Weak Zone Related Seismic Cycles in Progressive Failure Leading to Collapse in Brittle Crust

Until quite recently our understanding of the basic mechanical process responsible for earthquakes and faulting was not well known. It can be argued that this was partly a consequence of the complex nature of fracture in crust and in part because evidence of brittle phenomena in the natural laboratory of the earth is often obliterated or obscured by other geological processes. While it is well understood that the spatial and temporal complexity of earthquakes and the fault structures emerge from geometrical and material built-in heterogeneities, one important open question is how the shearing becomes localized into a band of intense fractures. Here the authors address these questions through a numerical approach of a tectonic plate by considering rockmass heterogeneity both in microscopic scale and in mesoscopic scale. Numerical simulations of the progressive failure leading to collapse under long-range slow driving forces in the far-field show earthquake-like rupture behavior. $En Echelon$ crack-arrays are reproduced in the numerical simulation. It is demonstrated that the underlying fracturing induced acoustic emissions (or seismic events) display self-organized criticality------from disorder to order. The seismic cycles and the geometric structures of the fracture faces, which are found greatly depending on the material heterogeneity (especially on the macroscopic scale), agree with that observed experimentally in real brittle materials. It is concluded that in order to predict a main shock, one must have extremely detailed knowledge on very minor features of the earth's crust far from the place where the earthquake originated. If correct, the model proposed here seemingly provides an explanation as to why earthquakes to date are not predicted so successfully. The reason is not that the authors do not understand earthquake mechanisms very well but that they still know little about our earth's crust.

Instability from steady and axisymmetric to steady and asymmetric floating half zone convection in a fat liquid bridge of larger Prandtl number

The linear instability analysis of the present paper shows that the thermocapillary convection in a half floating zone of larger Prandtl number has a steady instability mode w(i) = 0 and m = 1 for a fat liquid bridge V = 1.2 with small geometrical aspect ratio A = 0.6. This conclusion is different from the usual idea of hydrothermal instability, and implies that the instability of the system may excite a steady and axial asymmetric state before the onset of oscillation in the ease of large Prandtl number.

Two Bifurcation Transition Processes in Floating Half Zone Convection of Larger Prandtl Number Fluid

Processes of the onset oscillation in the thermocapillaxy convection under the Earth's gravity are investigated by the numerical simulation and experiments in a floating half zone of large Prandtl number with different volume ratio. Both computational and experimental results show that the steady and axisymmetric convection turns to the oscillatory convection of m=1 for the slender liquid bridge, and to the oscillatory convection before a steady and 3D asymmetric state for the case of a fat liquid bridge. It implies that, there are two critical Marangoni numbers related, respectively, to these two bifurcation transitions for the fat liquid bridge. The computational results agree with the results of ground-based experiments.

Influence of liquid bridge volume on instability of floating half zone convection

Thermocapillary instabilities on floating half zone convection in microgravity environment were investigated by linear instability analysis method. The critical Marangoni numbers were obtained and compared with the experimental ones. The influences of the liquid bridge volume and the aspect ratio on the critical Marangoni number were analyzed. It is found that the liquid bridge volume and the aspect ratio have great influence on the critical Marangoni number. There was a gap region where the oscillatory convection will not be observed in present analyses and in experiments in the curve of the critical Marangoni number vs the liquid bridge volume for the case of large Prandtl number and small aspect ratio.

Effects on Dense-phase Zone Temperature of Fluidized Bed with Uneven Air Distribution

由于采用非均匀布风,内旋流流化床的移动区空气量不足,导致燃烧不充分,温度较低。当移动区未流化时,密相区内存在较明显的温度不均匀性。随着移动区流速的提高,温度差迅速减小。当移动区流速超过2.0#mu#m后,密相区温度基本均匀一致。流动区流速对密相区温度均匀有一定的影响,流速越高,温度越均匀。

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection III. Three-dimensional model

An unsteady and three-dimensional model of the floating-half-zone convection on the ground is studied by the direct numerical simulation for the medium of 10 cSt silicon oil, and the influence of the liquid bridge volume on the critical applied temperature difference is especially discussed. The marginal curves for the onset of oscillation are separated into two branches related, respectively, to the slender liquid bridge and the fat liquid bridge. The oscillatory features of the floating-half-zone convection are also discussed.

Influence of Liquid Bridge Volume on the Floating Zone Convection

The liquid bridge volume is a critical geometrical parameter in addition to the aspect ratio for onset of oscillation in the floating zone convection. The oscillatory features are generally divided into two characteristic regions: slender liquid bridge region and fat liquid bridge region. The oscillatory modes in two regions are discussed in the present paper.

Statistical analysis of synthetic earthquake catalogs generated by models with various levels of fault zone disorder

The stress release model, a stochastic version of the elastic rebound theory, is applied to the large events from four synthetic earthquake catalogs generated by models with various levels of disorder in distribution of fault zone strength (Ben-Zion, 1996) They include models with uniform properties (U), a Parkfield-type asperity (A), fractal brittle properties (F), and multi-size-scale heterogeneities (M). The results show that the degree of regularity or predictability in the assumed fault properties, based on both the Akaike information criterion and simulations, follows the order U, F, A, and M, which is in good agreement with that obtained by pattern recognition techniques applied to the full set of synthetic data. Data simulated from the best fitting stress release models reproduce, both visually and in distributional terms, the main features of the original catalogs. The differences in character and the quality of prediction between the four cases are shown to be dependent on two main aspects: the parameter controlling the sensitivity to departures from the mean stress level and the frequency-magnitude distribution, which differs substantially between the four cases. In particular, it is shown that the predictability of the data is strongly affected by the form of frequency-magnitude distribution, being greatly reduced if a pure Gutenburg-Richter form is assumed to hold out to high magnitudes.

Proper orthogonal decomposition of oscillatory Marangoni flow in half-zone liquid bridges of low-Pr fluids

Proper orthogonal decomposition (POD) using method of snapshots was performed on three different types of oscillatory Marangoni flows in half-zone liquid bridges of low-Pr fluid (Pr = 0.01). For each oscillation type, a series of characteristic modes (eigenfunctions) have been extracted from the velocity and temperature disturbances, and the POD provided spatial structures of the eigenfunctions, their oscillation frequencies, amplitudes, and phase shifts between them. The present analyses revealed the common features of the characteristic modes for different oscillation modes: four major velocity eigenfunctions captured more than 99% of the velocity fluctuation energy form two pairs, one of which is the most energetic. Different from the velocity disturbance, one of the major temperature eigenfunctions makes the dominant contribution to the temperature fluctuation energy. On the other hand, within the most energetic velocity eigenfuction pair, the two eigenfunctions have similar spatial structures and were tightly coupled to oscillate with the same frequency, and it was determined that the spatial structures and phase shifts of the eigenfunctions produced the different oscillatory disturbances. The interaction of other major modes only enriches the secondary spatio-temporal structures of the oscillatory disturbances. Moreover, the present analyses imply that the oscillatory disturbance, which is hydrodynamic in nature, primarily originates from the interior of the liquid bridge. (C) 2007 Elsevier B.V. All rights reserved.

On the influence zone and the prediction of tensile strength of particulate polymeric composites

An intended numerical investigation is carried out. The results indicate that, even if a perfect adhesive bond is preserved between the particles and matrix materials, the two-phase element cell model is unable to predict the strength increment of the particulate polymeric composites (PPC). To explore the main reinforcing mechanism, additional microscopic experiment is performed. An ''influence zone'' was observed around each particle which is measured about 2 to 10 micrometers in thickness for a glass-polyethylene mixture. Then, an improved computational model is presented to include the ''influence zone'' effect and several mechanical behaviors of PPC are well simulated through this new model.

Transition of oscillatory floating half zone convection from Earth's gravity to microgravity

Oscillatory features of floating half zone convection were experimentally studied by using the drop shaft facility of Japan Microgravity Center which supported microgravity period of 10 s. Coordinated measurements including free surface deformation and oscillation, temperature and flow pattern in both 1-g and micro-g environment were obtained. The oscillatory frequency and amplitude in micro-g condition were lower and larger than the ones in l-g condition, respectively. The results gave, at first time, the oscillatory features such as free surface wave in micro-g, coordinated measurements of more than two physical quantities in the micro-g, and transition of thermocapillary oscillatory convection from I-g to micro-g.

Investigation on a simulation model of floating half zone convection .2. Experiment

A simulation model of floating half zone with non-uniform temperature distribution at the upper rod and uniform temperature distribution at lower rod was discussed by numerical investigation in a previous paper. In the present paper, the experimental investigation of the simulation model is given generally. The results of the present model show that the temperature profile is quite different and the critical applied temperature difference is lower than the one of usual model with same geometrical parameters in most cases. The features of critical Marangoni number depending on the liquid bridge volume are also different from the ones of usual model.

Experimental study of thermocapillary convection on a liquid bridge consisting of fluid with low Prandtl number in a floating half-zone

A device of mercury liquid bridge of floating half-zone is designed to experimentally explore thermocapillary convection and its instability of a low Prandtl number liquid. Noncontacted diagnostic techniques were developed to monitor surface flow and surface deformation. The surface flow and the influence of a growing surface film (or skin) on the flow were observed. It is shown that the film is a key factor in changing the behavior associated with the thermocapillary convection. The experiment indicates that the critical Marangoni number should be much higher than that expected by the numerical simulation. The condition and process of surface film growth are discussed. The surface oscillation of the mercury bridge wrapped with ''dirt-film'' was observed, and the characteristics and the frequency associated with this oscillation are given.

Numerical investigation on a simulation model of floating zone convection

A simulation model with adiabatic condition at the upper rod and constant temperature at the lower rod is studied numerically in this paper. The temperature distribution in a simulation model is closer to the one in the half part of a floating full zone in comparison with the one in a usual floating half zone model with constant temperature at both rods, because the temperature distribution of a floating full zone is symmetric for the middle plane in a microgravity environment. The results of the simulation model show that the temperature profiles and the how patterns are different from those of the usual floating half zone model. Another type of half zone model, with a special non-uniform temperature distribution at the upper rod and constant temperature at the lower rod, has been suggested by recent experiments. The temperature boundary condition of the upper rod has a maximum value in the center and a lower value near the free surface. This modified simulation model is also simulated numerically in the present paper. Copyright (C)1996 Elsevier Science Ltd.

A unified model for dislocation nucleation, dislocation emission and dislocation free zone

In this paper, a unified model for dislocation nucleation, emission and dislocation free zone is proposed based on the Peierls framework. Three regions are identified ahead of the crack tip. The emitted dislocations, located away from the crack tip in the form of an inverse pileup, define the plastic zone. Between that zone and the cohesive zone immediately ahead of the crack tip, there is a dislocation free zone. With the stress field and the dislocation density field in the cohesive zone and plastic zone being, respectively, expressed in the first and second Chebyshev polynomial series, and the opening and slip displacements in trigonometric series, a set of nonlinear algebraic equations can be obtained and solved with the Newton-Raphson Method. The results of calculations for pure shearing and combined tension and shear loading after dislocation emission are given in detail. An approximate treatment of the dynamic effects of the dislocation emission is also developed in this paper, and the calculation results are in good agreement with those of molecular dynamics simulations.