Dynamic Rupture in a 3-D Particle-based Simulation of a Rough Planar Fault     

Steefen Abe  Shane Latham  Peter Mora 《Pure and Applied Geophysics》2006,163(9):1881-1892

An appreciation of the physical mechanisms which cause observed seismicity complexity is fundamental to the understanding of the temporal behaviour of faults and single slip events. Numerical simulation of fault slip can provide insights into fault processes by allowing exploration of parameter spaces which influence microscopic and macroscopic physics of processes which may lead towards an answer to those questions. Particle-based models such as the Lattice Solid Model have been used previously for the simulation of stick-slip dynamics of faults, although mainly in two dimensions. Recent increases in the power of computers and the ability to use the power of parallel computer systems have made it possible to extend particle-based fault simulations to three dimensions. In this paper a particle-based numerical model of a rough planar fault embedded between two elastic blocks in three dimensions is presented. A very simple friction law without any rate dependency and no spatial heterogeneity in the intrinsic coefficient of friction is used in the model. To simulate earthquake dynamics the model is sheared in a direction parallel to the fault plane with a constant velocity at the driving edges. Spontaneous slip occurs on the fault when the shear stress is large enough to overcome the frictional forces on the fault. Slip events with a wide range of event sizes are observed. Investigation of the temporal evolution and spatial distribution of slip during each event shows a high degree of variability between the events. In some of the larger events highly complex slip patterns are observed.  相似文献   

Lattice preferred orientation in CaIrO<Subscript>3</Subscript> perovskite and post-perovskite formed by plastic deformation under pressure     

Ken Niwa  Takehiko Yagi  Kenya Ohgushi  Sébastien Merkel  Nobuyoshi Miyajima  Takumi Kikegawa 《Physics and Chemistry of Minerals》2007,34(9):679-686

Lattice preferred orientations (LPO) developed in perovskite and post-perovskite structured CaIrO₃ were studied using the radial X-ray diffraction technique combined with a diamond anvil cell. Starting materials of each phase were deformed from 0.1 MPa to 6 GPa at room temperature. Only weak LPO was formed in the perovskite phase, whereas strong LPO was formed in the post-perovskite phase with an alignment of the (010) plane perpendicular to the compression axis. The present result suggests that the (010) is a dominant slip plane in the post-perovskite phase and it is in good agreement with the crystallographic prediction, dislocation observations via transmission electron microscopy, and a recent result of simple shear deformation experiment at 1 GPa–1,173 K. However, the present result contrasts markedly from the results on MgGeO₃ and (Mg,Fe)SiO₃, which suggested that the (100) or (110) is a dominant slip plane with respect to the post-perovskite structure. Therefore it is difficult to discuss the behavior of the post-perovskite phase in the Earth’s deep interior based on existing data of MgGeO₃, (Mg,Fe)SiO₃ and CaIrO₃. The possible sources of the differences between MgGeO₃, (Mg,Fe)SiO₃ and CaIrO₃ are discussed.  相似文献   

The influence of pressure on the structure and dynamics of hydrogen bonds in zoisite and clinozoisite     

Björn Winkler  Julian D. Gale  Keith Refson  Dan J. Wilson  Victor Milman 《Physics and Chemistry of Minerals》2008,35(1):25-35

Density functional theory calculations have been used to study the pressure-induced changes of the hydrogen bond of Fe-free orthozoisite and clinozoisite and the concomitant shifts of the OH-stretching frequencies. Two independent parameter-free lattice dynamical calculations have been employed. One was based on a plane-wave basis set in conjunction with norm-conserving pseudopotentials and a density functional perturbation theory approach, while the other used a localised basis set and a finite displacement algorithm for the lattice dynamical calculations. Both models confirm the unusually large pressure-induced red-shift found experimentally (−33.89 cm⁻¹/GPa) in orthozoisite, while the pressure-induced shifts in clinozoisite are much smaller (−5 to −9 cm⁻¹/GPa). The atomistic model calculations show that in orthozoisite the nearly linear O–H⋯O arrangement is compressed by about 8% on a pressure increase to 10 GPa, while concomitantly the O–H distance is significantly elongated (by 2.5% at 10 GPa). In clinozoisite, the O–H⋯O arrangement is kinked at ambient conditions and remains kinked at high pressures, while the O-H distance is elongated by only 0.5% at 10 GPa. The current calculations confirm that correlations between the distances and dynamics of hydrogen bonds, which have been established at ambient conditions, cannot be used to infer hydrogen positions at high pressures.  相似文献   

Simulating flows in multi-layered and spatially-variable permeability media via a new Gray Lattice Boltzmann model     

《Computers and Geotechnics》2015

The simulation of flow in porous and fibrous permeable media is of high importance in many scientific and industrial applications. Although the finite element models at the representative elementary volume scales are used to solve a huge amount of scientific and engineering problems, they are hardly used to efficiently simulate pore-fluid flow problems at the particle scales. This encourages the development of numerical models to match the needs of such studies. In this paper, we propose a new Gray Lattice Boltzmann numerical model for simulating fluid flow in permeable media. Unlike most previous models, our proposed model has the ability to simulate multi-layers and space-variable permeability while preserving the continuity of the macroscopic velocity field. The model is verified with the available analytical solutions and a derived analytical expression for the case of variable porosity. In addition, we examine the importance of introducing a transition layer with a defined porosity function near the boundaries and interfaces. If this layer exists in practice, then the numerical results reveal that it cannot be neglected, and its impact is significant on the obtained velocity distribution. Finally, in the light of the obtained results, we can state that the proposed model has great potential to simulate complex and heterogeneous media with smoothness and accuracy, so that it may enrich the research content of the emerging computational geosciences.  相似文献   

A lattice Boltzmann model coupled with a Large Eddy Simulation model for flows around a groyne     

Xin YU  Xuelin TANG  Wuchang WANG  Fujun WANG  Zhicong CHEN  Xiaoyan SHI 《国际泥沙研究》2010,25(3):271-282

This present paper proposes a two-dimensional lattice Boltzmann model coupled with a Large Eddy Simulation (LES) model and applies it to flows around a non-submerged groyne in a channel. The LES of shallow water equations is efficiently performed using the Lattice Boltzmann Method (LBM) and the turbulence can be taken into account in conjunction with the Smagorinsky Sub-Grid Stress (SGS) model. The bounce-back scheme of the non-equilibrium part of the distribution function is used to determine the unknown distribution functions at inflow boundary, the zero gradient of the distribution function is set normal to outflow boundary to obtain the unknown distribution functions here and the bounce-back scheme, which states that an incoming particle towards the boundary is bounced back into fluid, is applied to the solid wall to ensure non-slip boundary conditions. The initial flow field is defined firstly and then is used to calculate the local equilibrium distributions as initial conditions of the distribution functions. These coupled models successfully predict the flow characteristics, such as circulating flow, velocity and water depth distributions. The comparisons between the simulated results and the experimental data show that the model scheme has the capacity to solve the complex flows in shallow water with reasonable accuracy and reliability.  相似文献   

Modeling granular soils liquefaction using coupled lattice Boltzmann method and discrete element method     

《Soil Dynamics and Earthquake Engineering》2014

In this paper, a novel coupled pore-scale model of pore-fluid interacting with discrete particles is presented for modeling liquefaction of saturated granular soil. A microscale idealization of the solid phase is achieved using the discrete element method (DEM) while the fluid phase is modeled at a pore-scale using the lattice Boltzmann method (LBM). The fluid forces applied on the particles are calculated based on the momentum exchange between the fluid and particles. The presented model is based on a first principles formulation in which pore-pressure develops due to actual changes in pore space as particles׳ rearrangement occurs during shaking. The proposed approach is used to model the response of a saturated soil deposit subjected to low and large amplitude seismic excitations. Results of conducted simulations show that at low amplitude shaking, the input motion propagates following the theory of wave propagation in elastic solids. The deposit response to the strong input motion indicates that liquefaction took place and it was due to reduction in void space during shaking that led to buildup in pore-fluid pressure. Soil liquefaction was associated with soil stiffness degradation and significant loss of interparticle contacts. Simulation results also indicate that the level of shaking-induced shear strains and associated volumetric strains play a major role in the onset of liquefaction and the rate of pore-pressure buildup.  相似文献   

A new lattice Boltzmann model for interface reactions between immiscible fluids     

《Advances in water resources》2015

In this paper, we describe a lattice Boltzmann model to simulate chemical reactions taking place at the interface between two immiscible fluids. The phase-field approach is used to identify the interface and its orientation, the concentration of reactant at the interface is then calculated iteratively to impose the correct reactive flux condition. The main advantages of the model is that interfaces are considered part of the bulk dynamics with the corrective reactive flux introduced as a source/sink term in the collision step, and, as a consequence, the model’s implementation and performance is independent of the interface geometry and orientation. Results obtained with the proposed model are compared to analytical solution for three different benchmark tests (stationary flat boundary, moving flat boundary and dissolving droplet). We find an excellent agreement between analytical and numerical solutions in all cases. Finally, we present a simulation coupling the Shan Chen multiphase model and the interface reactive model to simulate the dissolution of a collection of immiscible droplets with different sizes rising by buoyancy in a stagnant fluid.  相似文献   

Does pore water velocity affect the reaction rates of adsorptive solute transport in soils? Demonstration with pore-scale modelling     

Xiaoxian Zhang  John W. Crawford  Iain M. Young 《Advances in water resources》2008

Modelling adsorptive solute transport in soils needs a number of parameters to describe its reaction kinetics and the values of these parameters are usually determined from batch and displacement experiments. Some experimental results reveal that when describing the adsorption as first-order kinetics, its associated reaction rates are not constants but vary with pore water velocity. Explanation of this varies but an independent verification of each explanation is difficult because simultaneously measuring the spatiotemporal distributions of dissolved and adsorbed solutes in soils is formidable. Pore-scale modelling could play an important role to address this gap and has received increased attention over the past few years. This paper investigated the transport of adsorptive solute in a simple porous medium using pore-scale modelling. Fluid flow through the void space of the medium was assumed to be laminar and in saturated condition, and solute transport consisted of advection and molecular diffusion; the sorption and desorption occurring at the fluid–solid interface were modelled as linear first-order kinetics. Based on the simulated spatiotemporal distribution of dissolved and adsorbed solutes at pore scale, volumetric-average reaction kinetics at macroscopic scale and its associated reactive parameters were measured. Both homogeneous adsorption where the reaction rates at microscopic scale are constant, and heterogeneous adsorption where the reaction rates vary from site to site, were investigated. The results indicate that, in contrast to previously thought, the macroscopic reaction rates directly measured from the pore-scale simulations do not change with pore velocity under both homogeneous and heterogeneous adsorptions. In particular, we found that for the homogeneous adsorption, the macroscopic adsorption remains first-order kinetic and can be described by constant reaction rates, regardless of flow rate; whilst for the heterogeneous adsorption, the macroscopic adsorption kinetics continues not to be affected by flow rate but is no longer first-order kinetics that can be described by constant reaction rates. We discuss how these findings could help explain some contrary literature reports over the dependence of reaction rates on pore water velocity.  相似文献   

Dynamic Evolution of Precursory Field and Its Relation With the Group of Strong Earthquakes     

Li Li  Zhang Guomin 《中国地震研究》2000,14(2):122-135

Spatial scanning is done for two regions in Chinese Mainland,where displayed a denseprecursory network during 1994～1998.The two regions are the mid-southern segment of theNorth-south seismic belt(20°～35°N,95°～110°E)and North China(36°～42°N,110°～120°E).We took 0.5°×0.5°as a spatial window with a step of 0.25°and 4 months as atemporal window with a step of 1 month.For the two regions,the anomaly density is scannedfrom 1994 and 1995 respectively in the two regions.The precursory anomalies are all fromthe Division of Seismic Trend in China and the Division of Seismic Trend in the Capital Area,Center for Analysis and Prediction,China Seismological Bureau.A seismogenic tectonicmodel is introduced to explain the scanning results.In the model,the frictional strength ofthe focal sources is distributed randomly.After the boundary plate motion rate and all othergeological parameters are given,the stress of the sources in the system changesinhomogenously due to the variation of the frictional streng  相似文献   

The structure of grain boundaries in granite-origin ultramylonite studied by high-resolution electron microscopy     

T. Hiraga  T. Nagase  M. Akizuki 《Physics and Chemistry of Minerals》1999,26(8):617-623

The structure of grain boundaries in a granite-origin ultramylonite, composed mainly of fine-grained feldspar and quartz, was studied by high-resolution electron microscopy (HREM). At most of the boundaries, not only between the same minerals but also between different minerals, lattice fringes in adjacent grains meet at the interface with no other appreciable phases. In these boundaries, some of the straight segments correspond to a low-index plane of one of the connected grains. Boundaries containing voids, with a spheroidal shape elongated along the boundaries, were observed only between quartz grains. It is suggested that these boundaries were formed by healing of microcracks. The structural width of major boundaries, deduced from lattice-fringe imaging, is less than about 0.5 nm. Received: 15 September 1998 / Revised, accepted: 8 April 1999  相似文献