Simulation of the micro-physics of rocks using LSMearth

The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (MORA and PLACE, 1994; PLACE and MORA, 1999). A new modular and flexible LSM approach has been developed that allows different microphysics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.

Predicting cave initiation and propagation

The most widely used method for predicting the onset of continuous caving is Laubscher's caving chart. A detailed examination of this method was undertaken which concluded that it had limitations which may impact on results, particularly when dealing with stronger rock masses that are outside current experience. These limitations relate to inadequate guidelines for adjustment factors to rock mass rating (RMR), concerns about the position on the chart of critical case history data, undocumented changes to the method and an inadequate number of data points to be confident of stability boundaries. A review was undertaken on the application and reliability of a numerical method of assessing cavability. The review highlighted a number of issues, which at this stage, make numerical continuum methods problematic for predicting cavability. This is in particular reference to sensitivity to input parameters that are difficult to determine accurately and mesh dependency. An extended version of the Mathews method for open stope design was developed as an alternative method of predicting the onset of continuous caving. A number of caving case histories were collected and analyzed and a caving boundary delineated statistically on the Mathews stability graph. The definition of the caving boundary was aided by the existence of a large and wide-ranging stability database from non-caving mines. A caving rate model was extrapolated from the extended Mathews stability graph but could only be partially validated due to a lack of reliable data.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

A review of drainage and spontaneous rupture in free standing thin films with tangentially immobile interfaces

A review of spontaneous rupture in thin films with tangentially immobile interfaces is presented that emphasizes the theoretical developments of film drainage and corrugation growth through the linearization of lubrication theory in a cylindrical geometry. Spontaneous rupture occurs when corrugations from adjacent interfaces become unstable and grow to a critical thickness. A corrugated interface is composed of a number of waveforms and each waveform becomes unstable at a unique transition thickness. The onset of instability occurs at the maximum transition thickness, and it is shown that only upper and lower bounds of this thickness can be predicted from linear stability analysis. The upper bound is equivalent to the Freakel criterion and is obtained from the zeroth order approximation of the H-3 term in the evolution equation. This criterion is determined solely by the film radius, interfacial tension and Hamaker constant. The lower bound is obtained from the first order approximation of the H-3 term in the evolution equation and is dependent on the film thinning velocity A semi-empirical equation, referred to as the MTR equation, is obtained by combining the drainage theory of Manev et al. [J. Dispersion Sci. Technol., 18 (1997) 769] and the experimental measurements of Radoev et al. [J. Colloid Interface Sci. 95 (1983) 254] and is shown to provide accurate predictions of film thinning velocity near the critical thickness of rupture. The MTR equation permits the prediction of the lower bound of the maximum transition thickness based entirely on film radius, Plateau border radius, interfacial tension, temperature and Hamaker constant. The MTR equation extrapolates to Reynolds equation under conditions when the Plateau border pressure is small, which provides a lower bound for the maximum transition thickness that is equivalent to the criterion of Gumerman and Homsy [Chem. Eng. Commun. 2 (1975) 27]. The relative accuracy of either bound is thought to be dependent on the amplitude of the hydrodynamic corrugations, and a semiempirical correlation is also obtained that permits the amplitude to be calculated as a function of the upper and lower bound of the maximum transition thickness. The relationship between the evolving theoretical developments is demonstrated by three film thickness master curves, which reduce to simple analytical expressions under limiting conditions when the drainage pressure drop is controlled by either the Plateau border capillary pressure or the van der Waals disjoining pressure. The master curves simplify solution of the various theoretical predictions enormously over the entire range of the linear approximation. Finally, it is shown that when the Frenkel criterion is used to assess film stability, recent studies reach conclusions that are contrary to the relevance of spontaneous rupture as a cell-opening mechanism in foams. (C) 2003 Elsevier Science B.V. All rights reserved.

Simultaneous FT-NIR and ESR analyses to yield propagation rate coefficients for polymerization of methyl methacrylate based monomers

A hyphenated instrumental approach has been used to obtain reliable values for the propagation rate coefficients as a function of conversion for polymerizations of methyl methacrylate (MMA) and a mixture of MMA and ethyleneglycol dimethacrylate (EGDMA) with a 1:1 concentration of double bonds, from near the onset of the Trommsdorf region into the glass region. ESR spectroscopy was used to measure the radical concentration while FT-NIR fibre-optic spectroscopy was employed to measure instantaneously the double-bond concentration within the temperature-controlled cavity of the ESR instrument during polymerization. The advantage of this approach to the measurement of the rate coefficient is that it is equally applicable to branching and linear polymerizations. For the polymerization of methyl methacrylate, the values of the rate coefficient at the lowest conversions at which reliable values could be obtained were in agreement with recently reported values obtained by the PLP-SEC method. For the lowest conversions, the values obtained were 403 1 mol(-1) s(-1) at 306 K for MMA and 5201 mol(-1) s(-1) at 310 K for a 1:1 mixture of MMA and EGDMA. (C) 2003 Society of Chemical Industry.

Influence of leaf retention on cutting propagation of Lavandula dentata L.

Vegetative propagation of lavender offers several advantages over sexual propagation, among them crop homogeneity and yield of higher quality essential oil. However, Lavandula species have been propagated mostly by seeds and are said to be recalcitrant to rooting when propagated by cuttings. During cutting propagation, one of the important variables that influence the rooting capacity of cuttings is the leaf retention. The objective of this work was to evaluate the influence of leaf retention on rooting of L. dentata cuttings. Apical cuttings of L. dentata of 10 cm in length, keeping approximately 1/3, 1/2 or 2/3 of their leaves were planted in commercial substrate Plantmax HT® under intermittent mist. After two months, averages of root number, length of the longest root, root fresh and dry weights, and the survival percentage were evaluated. Root length and fresh weight were statistically greater with 2/3 of leaf retention and when fewer leaves were kept on the cuttings, lower means of root dry weight was observed. Under the conditions applied in this study, greater leaf retention was better for rooting of L. dentata cuttings.

Chemical sterilization in in vitro propagation of Arundina bambusifolia Lindl. and Epidendrum ibaguense Kunth

There is a great demand for simpler and less costly laboratory techniques and for more accessible procedures for orchid breeders who do not have the necessary theoretical basis to use the traditional seed and clone production methods of orchids in vitro. The aim of this study was to assess the use of sodium hypochlorite (NaClO) as a decontaminant in the process of inoculating adult orchid explants of Arundina bambusifolia and Epidendrum ibaguenses. Solutions of NaClO (1.200, 2.400, 3.600, 4.800 and 6.000 mg L-1 - equivalent to 50, 100, 150, 200 and 250 mL L-1 of commercial bleach - CB) were sprayed on the explants (1.0 mL) and the culture medium (GB5), in the presence or absence of activated charcoal (2 g L-1). The explants used were nodal segments of field-grown adult plants. The procedures for inoculating the explants were conducted outside the laminar flow chamber (LFC), except for the control treatment (autoclaved medium and explant inoculation inside the LFC). The best results for fresh weight yield, height and number of shoots were obtained using NaClO in solution at 1.200 mg L-1 (equivalent to 50 mL L-1 commercial bleach) with activated charcoal in the culture medium. Fresh weight figures were 1.10 g/jar for Arundina bambusifolia and 0.16 g/jar for Epidendrum ibaguenses. Spraying the NaClO solutions controls the contamination of the culture medium already inoculated with the explants.

Rooting of hardwood cuttings of Roxo de Valinhos fig (Ficus carica L.) with different propagation strategies

The objective of this study was to evaluate the substrate, cuttings collection time, the position and the cutting depth, and the propagation environment on rooting of 'Purple Valinhos' fig tree cuttings in Southwestern Paraná, Brazil. Two experiments were carried out at UTFPR, Câmpus Dois Vizinhos, with hardwoods cuttings from Roxo de Valinhos fig tree. The first experiment used a randomized block design, in 3 x 3 x 2 factorial (substrate x environment x collection time), with four replications of 10 cuttings per plot. The cuttings were collected in the first fifteen days of July and August. The substrates were sand, soil and the mixture of these [1:1 (v / v)]. The environments used were open sky, tunnel with plastic cover and tunnel with half-shade black net cover. The second experiment used a randomized block design, 2 x 2 x 3 factorial (shoot cutting position x soil cover x shoot cutting depth), with four replications of 12 cuttings per plot. In the factor position, the vertically (0 º inclination) and inclined (45 º inclination) shoot cuttings were evaluated. Soil cover was tested with mulching plastic cover or not. The tested depths were 1/3, 1/2 and 2/3 in relation to the total length of the shoot cutting. In both experiments, the following were analyzed: rooting and mortality indices, number of leaves and primary shoots, length of the three largest roots per cutting. It was conclude that, the protected environment with plastic cover on sand as substrate must recommended for the rooting of fig estaca, collecting them in the first half of July. The inclination position and cutting depth of the estaca and the substrate coverage with plastic mulching did not influence the results.

Evaluating Tsunami Impact on the Gulf of Cadiz Coast (Northeast Atlantic).

The Gulf of Cadiz coasts are exposed to tsunamis. Emergency planning tools are now taking into account this fact, especially because a series of historical occurrences were strikingly significant, having left strong evidence behind, in the mareographic records, the geological evidence or simply the memory of the populations. The study area is a strip along the Algarve coast, south Portugal, an area known to have been heavily impacted by the 1 November 1755 event. In this study we use two different tsunami scenarios generated by the rupture of two thrust faults identified in the area, corresponding to 8.1-8.3 magnitude earthquakes. Tsunami propagation and inundation computation is performed using a non-linear shallow water code with bottom friction. Numerical modeling results are presented in terms of flow depth and current velocity with maximum values of 7 m and 8 m/s for inundation depth and flow speed, respectively. These results constitute a valuable tool for local authorities, emergency and decision planners to define the priority zones where tsunami mitigation measures must be implemented and to develop tsunami-resilient communities.

From local to global importance measures of uncertainty propagation in composite structures

The influence of uncertainties of input parameters on output response of composite structures is investigated in this paper. In particular, the effects of deviations in mechanical properties, ply angles, ply thickness and on applied loads are studied. The uncertainty propagation and the importance measure of input parameters are analysed using three different approaches: a first-order local method, a Global Sensitivity Analysis (GSA) supported by a variance-based method and an extension of local variance to estimate the global variance over the domain of inputs. Sample results are shown for a shell composite laminated structure built with different composite systems including multi-materials. The importance measures of input parameters on structural response based on numerical results are established and discussed as a function of the anisotropy of composite materials. Needs for global variance methods are discussed by comparing the results obtained from different proposed methodologies. The objective of this paper is to contribute for the use of GSA techniques together with low expensive local importance measures.

Uncertainty propagation in inverse reliability-based design of composite structures

An approach for the analysis of uncertainty propagation in reliability-based design optimization of composite laminate structures is presented. Using the Uniform Design Method (UDM), a set of design points is generated over a domain centered on the mean reference values of the random variables. A methodology based on inverse optimal design of composite structures to achieve a specified reliability level is proposed, and the corresponding maximum load is outlined as a function of ply angle. Using the generated UDM design points as input/output patterns, an Artificial Neural Network (ANN) is developed based on an evolutionary learning process. Then, a Monte Carlo simulation using ANN development is performed to simulate the behavior of the critical Tsai number, structural reliability index, and their relative sensitivities as a function of the ply angle of laminates. The results are generated for uniformly distributed random variables on a domain centered on mean values. The statistical analysis of the results enables the study of the variability of the reliability index and its sensitivity relative to the ply angle. Numerical examples showing the utility of the approach for robust design of angle-ply laminates are presented.

A VLF to UHF propagation model for evaluation of high-speed railway electromagnetic interference in RF monitoring base-stations

Considering that recent european high-speed railway system has a traction power system of kV 50 Hz, which causes electromagnetic emission for the outside world, it is important to dimension the railway system emissions, using a frequency/distance dependent propagation model. This paper presents an enhanced theoretical model for VLF to UHF propagation, railway system oriented. It introduces the near field approach (crucial in low frequency propagation) and also considers the source characteristics and type of measuring antenna. Simulations are presented, and comparisons are set with earlier far field models. Using the developed model, a real case study was performed in partnership with Refer Telecom (portuguese telecom operator for railways). The new propagation model was used in order to predict the future high-speed railway electromagnetic emissions in the Lisbon north track. The results show the model's prediction capabilities and also its applicability to realistic scenarios.

Halting the fuse discharge propagation using optical fiber microwires

We report and analyze the halting of the fuse effect propagation in optical fiber microwires. The increase of the mode field diameter in the tapered region decreases the optical intensity resulting in the extinction of the fuse effect. This fiber element presents a low insertion loss and can be introduced in the optical network in order to protect the active equipment from the damage caused by the fuse effect.

Halting the fuse discharge propagation Using optical fiber microwires

We report and analyze the halting of the fuse effect propagation in optical fiber microwires. The increase of the mode field diameter in the tapered region decreases the optical intensity resulting in the extinction of the fuse effect. This fiber element presents a low insertion loss and can be introduced in the optical network in order to protect the active equipment from the damage caused by the fuse effect. (C) 2012 Optical Society of America

Fractional Dynamics of Computer Virus Propagation

We propose a fractional model for computer virus propagation. The model includes the interaction between computers and removable devices. We simulate numerically the model for distinct values of the order of the fractional derivative and for two sets of initial conditions adopted in the literature. We conclude that fractional order systems reveal richer dynamics than the classical integer order counterpart. Therefore, fractional dynamics leads to time responses with super-fast transients and super-slow evolutions towards the steady-state, effects not easily captured by the integer order models.