共查询到20条相似文献，搜索用时 31 毫秒
1.
-- In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by Dieterich and Kilgore (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance Dc on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments. 相似文献
2.
Despite the insight gained from 2-D particle models, and given that the dynamics of crustal faults occur in 3-D space, the question remains, how do the 3-D fault gouge dynamics differ from those in 2-D? Traditionally, 2-D modeling has been preferred over 3-D simulations because of the computational cost of solving 3-D problems. However, modern high performance computing architectures, combined with a parallel implementation of the Lattice Solid Model (LSM), provide the opportunity to explore 3-D fault micro-mechanics and to advance understanding of effective constitutive relations of fault gouge layers. In this paper, macroscopic friction values from 2-D and 3-D LSM simulations, performed on an SGI Altix 3700 super-cluster, are compared. Two rectangular elastic blocks of bonded particles, with a rough fault plane and separated by a region of randomly sized non-bonded gouge particles, are sheared in opposite directions by normally-loaded driving plates. The results demonstrate that the gouge particles in the 3-D models undergo significant out-of-plane motion during shear. The 3-D models also exhibit a higher mean macroscopic friction than the 2-D models for varying values of interparticle friction. 2-D LSM gouge models have previously been shown to exhibit accelerating energy release in simulated earthquake cycles, supporting the Critical Point hypothesis. The 3-D models are shown to also display accelerating energy release, and good fits of power law time-to-failure functions to the cumulative energy release are obtained. 相似文献
3.
In this study, 3-D Lattice Solid Model (LSMearth or LSM) was extended by introducing particle-scale rotation. In the new model,
for each 3-D particle, we introduce six degrees of freedom: Three for translational motion, and three for orientation. Six
kinds of relative motions are permitted between two neighboring particles, and six interactions are transferred, i.e., radial,
two shearing forces, twisting and two bending torques. By using quaternion algebra, relative rotation between two particles
is decomposed into two sequence-independent rotations such that all interactions due to the relative motions between interactive
rigid bodies can be uniquely decided. After incorporating this mechanism and introducing bond breaking under torsion and bending
into the LSM, several tests on 2-D and 3-D rock failure under uni-axial compression are carried out. Compared with the simulations
without the single particle rotational mechanism, the new simulation results match more closely experimental results of rock
fracture and hence, are encouraging. Since more parameters are introduced, an approach for choosing the new parameters is
presented. 相似文献
4.
Ryszard Staroszczyk 《Acta Geophysica》2011,59(5):833-857
A multi-grain model for a migration recrystallization process in polar ice is presented. The model is based on the Sachs-Reuss
approximation of the stress homogeneity in a polycrystalline aggregate. An individual crystal of ice is treated as a transversely
isotropic and incompressible medium which deforms by viscous creep. The highly anisotropic viscous behaviour of the ice crystal
is described by a constitutive law expressing microscopic strain-rate in terms of the deviatoric stress and three fluidity
parameters that define different viscous resistances of the crystal in different glide directions. It is assumed that the
recrystallization occurs in those crystals in the aggregate which are most slowly deforming, and new crystals are nucleated
at orientations which favour the crystal deformation by basal glide. The model predictions are illustrated by results of numerical
simulations of simple flows, showing the evolution of the microscopic structure of ice and the variation of macroscopic viscosities
with increasing deformations. 相似文献
5.
--The earthquake generation cycle consists of tectonic loading, quasi-static rupture nucleation, dynamic rupture propagation and stop, and subsequent stress redistribution and fault restrengthening. From a macroscopic point of view, the entire process of earthquake generation cycles should be consistently described by a coupled nonlinear system of a slip-response function, a fault constitutive law and a driving force. On the basis of such a general idea, we constructed a realistic 3-D simulation model for earthquake generation cycles at a transcurrent plate boundary by combining the viscoelastic slip-response function derived for a two-layered elastic-viscoelastic structure model, the slip- and time-dependent fault constitutive law that has an inherent mechanism of fault restrengthening, and the steady relative plate motion as a driving force into a single closed system. With this model we numerically simulated the earthquake generation cycles repeated in a seismogenic region on a plate interface, and examined space-time changes in shear stress, slip deficits and fault constitutive properties during one complete cycle in detail. The occurrence of unstable dynamic slip brings about decrease both in fault strength and shear stress to a constant residual level. After the arrest of dynamic slip, the breakdown strength drop j†p of fault is restored rapidly and the process of stress accumulation resumes in the seismogenic region. On the other hand, the restoration of the critical weakening displacement Dc proceeds gradually with time through the interseismic period. The restoration of Dc can be regarded as the macroscopic manifestation of the microscopic recovery process of fractal fault surface structure. Through numerical simulation with a multi-segmented fault model, we examined the effects of viscoelastic fault-to-fault interaction. The effect of transient viscoelastic stress transfer through the asthenosphere is significant as well as the direct effect of elastic stress transfer, and it possibly explains the time lag of the sequential occurrence of large events along a plate boundary. 相似文献
6.
实际地震勘探中,储层物性参数的差异是导致地震波响应特征发生变化的根本原因,而建立储层物性参数与地震响应特征之间的联系,需要跨越微观孔隙尺度、介观测井尺度以及宏观地震尺度等三个不同尺度空间.本文基于已知井的岩石物理实验数据和测井数据,利用复杂多孔隙介质理论将微观尺度孔隙岩石粗化到介观测井尺度,利用Backus平均理论将介观测井尺度的模型进一步粗化到宏观地震尺度,最终,得到地震尺度裂缝-多孔隙介质模型.其数值计算结果与测井数据和地震数据的对比表明：基于两级尺度粗化算法的裂缝多孔隙介质模型在给定参数下是有效的,且基于该模型的地震响应特征分析方法能够对储层的地震响应特征随物性参数的变化进行分析. 相似文献
7.
Experimental and numerical modelling of sedimentation in a rectangular shallow basin 总被引：9，自引：6，他引：3
Sameh A. KANTOUSH Erik BOLLAERT and Anton J. SCHLEISS Dr. Ecole Polytechnique Fédérale de Lausanne 《国际泥沙研究》2008,23(3):212-232
Numerical simulation of flows in shallow reservoirs has to be checked for its consistency in predicting real flow conditions and sedimentation patterns. Typical flow patterns may exhibit flow separation at the inlet, accompanied by several recirculation and stagnation areas all over the reservoir surface. The aim of the present research project is to study the influence of the geometry of a reservoir on sediment transport and deposition numerically and experimentally, focusing on a prototype reservoir depth between 5 and 15 m as well as suspended sediment transport.
A series of numerical simulations is presented and compared with scaled laboratory experiments, with the objective of testing the sensitivity to different flow and sediment parameters and different turbulence closure schemes. Different scenarios are analyzed and a detailed comparison of preliminary laboratory tests and some selected simulations are presented.
The laboratory experiments show that suspended sediment transport and deposition are determined by the initial flow pattern and by the upstream and downstream boundary conditions. In the experiments, deposition in the rectangular basin systematically developed along the left bank, although inflow and outflow were positioned symmetrically along the centre of the basin. Three major horizontal eddies developed influencing the sediment deposition pattern. Although asymmetric flow patterns are privileged, a symmetric pattern can appear from time to time. This particular behaviour could also be reproduced by a two-dimensional depth-averaged flow and sediment transport model （CCHE2D）. The paper presents numerical simulations using different turbulence closure schemes （k-ε and eddy viscosity models）. In spite of the symmetric setup, these generally produced an asymmetric flow pattern that can easily switch sides depending on the assumptions made for the initial and boundary conditions. When using the laboratory experiment as a reference, the most reliable numerical results have been obtai 相似文献
A series of numerical simulations is presented and compared with scaled laboratory experiments, with the objective of testing the sensitivity to different flow and sediment parameters and different turbulence closure schemes. Different scenarios are analyzed and a detailed comparison of preliminary laboratory tests and some selected simulations are presented.
The laboratory experiments show that suspended sediment transport and deposition are determined by the initial flow pattern and by the upstream and downstream boundary conditions. In the experiments, deposition in the rectangular basin systematically developed along the left bank, although inflow and outflow were positioned symmetrically along the centre of the basin. Three major horizontal eddies developed influencing the sediment deposition pattern. Although asymmetric flow patterns are privileged, a symmetric pattern can appear from time to time. This particular behaviour could also be reproduced by a two-dimensional depth-averaged flow and sediment transport model （CCHE2D）. The paper presents numerical simulations using different turbulence closure schemes （k-ε and eddy viscosity models）. In spite of the symmetric setup, these generally produced an asymmetric flow pattern that can easily switch sides depending on the assumptions made for the initial and boundary conditions. When using the laboratory experiment as a reference, the most reliable numerical results have been obtai 相似文献
8.
-- The seismic records from significant earthquakes are profoundly affected by 3-D variations in crustal structure both in the source zone itself and in propagation to some distance. Even in structurally complex zones such as Japan and Mexico relatively coherent arrivals are found associated with different classes of propagation paths. The presence of strong lateral variations can disrupt the arrivals, and impose significant variations in propagation characteristics for different directions from the source as illustrated by observations for the 1995 Kobe and 2000 Tottori-ken Seibu earthquakes in western Japan. Such effects can be modelled in 3 dimensions using a hybrid scheme with a pseudospectral representation for horizontal coordinates and finite differences in depth. This arrangement improves parallel implementation by minimising communication costs. For a realistic 3-D model for the structure in western Japan the 3-D simulations to frequencies close to 1 Hz provide a good representation of the observations from subduction zones events such as the 1946 Nankai earthquake and the 2000 Tottori-ken Seibu earthquake. The model can therefore be used to investigate the pattern of ground motion expected for future events e.g., in current seismic gaps. 相似文献
9.
Analytical and numerical analysis of pressure drawdown in a poroelastic reservoir with complete overburden effect considered 总被引：1，自引：0，他引：1
Starting from an analytical reservoir model that incorporates full interaction with an elastic overburden, a new hybrid mathematical approach is developed by combining two numerical discretization methods. A tabular reservoir (petroleum reservoir or an aquifer) in an infinite or semi-infinite domain is viewed as a macroscopic displacement discontinuity, allowing use of the efficient displacement discontinuity mathematical method to calculate stresses and displacements that arise because of pressure changes. A 3-D finite element method using a poroelastic formulation is used to discretize the reservoir itself. By coupling the displacement discontinuity and finite element methods, a 3-D large-scale poroelastic reservoir can be simulated within an infinite or semi-infinite domain. The numerical model has been verified through comparison to known solutions, and some time-dependent pressure drawdown problems are analyzed. Results indicate that including the complete overburden (reservoir surroundings) response has a significant effect on pressure drawdown in a poroelastic reservoir during pumping, and should be incorporated in appropriate applications such as well test equations and subsidence analyses. 相似文献
10.
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. 相似文献
11.
Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model 总被引：4，自引：0，他引：4
— The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes. 相似文献
12.
-- Fault zones are thought to consist of regions with reduced seismic velocity. When sources are located in or close to these low-velocity zones, guided seismic head and trapped waves are generated which may be indicative of the structure of fault zones at depth. Observations above several fault zones suggest that they are common features of near fault radiation, yet their interpretation may be highly ambiguous. Analytical methods have been developed to calculate synthetic seismograms for sources in fault zones as well as at the material discontinuities. These solutions can be used for accurate modeling of wave propagation in plane-parallel layered fault zone structures. However, at present it is not clear how modest deviations from such simplified geometries affect the generation efficiency and observations of trapped wave motion. As more complicated models cannot be solved by analytical means, numerical methods must be employed. In this paper we discuss 3-D finite-difference calculations of waves in modestly irregular fault zone structures. We investigate the accuracy of the numerical solutions for sources at material interfaces and discuss some dominant effects of 3-D structures. We also show that simple mathematical operations on 2-D solutions generated with line sources allow accurate modeling of 3-D wave propagation produced by point sources. The discussed simulations indicate that structural discontinuities of the fault zone (e.g., fault offsets) larger than the fault zone width affect significantly the trapping efficiency, while vertical properly gradients, fault zone narrowing with depth, small-scale structures, and moderate geometrical variations do not. The results also show that sources located with appropriate orientations outside and below a shallow fault zone layer can produce considerable guided wave energy in the overlying fault zone layer. 相似文献
13.
-- A new technique for the parallel computing of 3-D seismic wave propagation simulation is developed by hybridizing the Fourier pseudospectral method (PSM) and the finite-difference method (FDM). This PSM/FDM hybrid offers a good speed-up rate using a large number of processors. To show the feasibility of the hybrid, a numerical 3-D simulation of strong ground motion was conducted for the 1999 Chi-Chi, Taiwan earthquake (Mw 7.6). Comparisons between the simulation results and observed waveforms from a dense strong ground motion network in Taiwan clearly demonstrate that the variation of the subsurface structure and the complex fault slip distribution greatly affect the damage during the Chi-Chi earthquake. The directivity effect of the fault rupture produced large S-wave pulses along the direction of the rupture propagation. Slips in the shallow part of the fault generate significant surface waves in Coastal Plain along the western coast. A large velocity gradient in the upper crust can propagate seismic waves to longer distances with minimum attenuation. The S waves and surface waves were finally amplified further by the site effect of low-velocity sediments in basins, and caused the significant disasters. 相似文献
14.
H. Wiechen 《Annales Geophysicae》1999,17(5):595-603
Magnetic reconnection is a process of fundamental importance for the dynamics of the Earth’s plasma sheet. In this context, the development of thin current sheets in the near-Earth plasma sheet is a topic of special interest because they could be a possible cause of microscopic fluctuations acting as collective non-idealness from a macroscopic point of view. Simulations of the near-Earth plasma sheet including boundary perturbations due to localized inflow through the northern (or southern) plasma sheet boundary show developing thin current sheets in the near-Earth plasma sheet about 8–10 RE tailwards of the Earth. This location is largely independent from the localization of the perturbation. The second part of the paper deals with the problem of the macroscopic non-ideal consequences of microscopic fluctuations. A new model is presented that allows the quantitative calculation of macroscopic non-idealness without considering details of microscopic instabilities or turbulence. This model is only based on the assumption of a strongly fluctuating, mixing dynamics on microscopic scales in phase space. The result of this approach is an expression for anomalous non-idealness formally similar to the Krook resistivity but now describing the macroscopic consequences of collective microscopic fluctuations, not of collisions. 相似文献
15.
Displacement capacity of masonry piers: parametric numerical analyses versus international building codes 总被引：1，自引：1，他引：0
Maurizio Orlando Luca Salvatori Paolo Spinelli Mario De Stefano 《Bulletin of Earthquake Engineering》2016,14(8):2259-2271
The nonlinear behaviour of masonry piers loaded in their plane is investigated by parametric numerical simulations. Each pier has a cantilever scheme, is loaded by a constant axial load and is subjected to an increasing horizontal displacement at the top. The macro-modelling approach is used to perform numerical analyses, adopting two different constitutive laws: a total strain crack model and a plastic model. The numerical model is calibrated on a block-masonry type for which experimental tests are available in literature. Parametric numerical simulations are performed by varying the aspect-ratio and the compression level, in order to assess the influence of such parameters on both shear strength and displacement capacity. By comparing numerical results with formulas of international codes, a good agreement for the shear strength is obtained, while significant differences are observed for the displacement capacity, which is influenced by both parameters. The authors propose a simple empirical formula for the displacement capacity, obtained by fitting the numerical results. The expression can be useful in the practical design for considering the influence of aspect-ratio and compression level, currently neglected by building codes. 相似文献
16.
用数值方法和岩石力学实验模拟地震断层的破裂过程时, 模型或试件的尺寸选取不合适将影响断层附近的应力场. 本文利用二维弹性有限元方法, 计算了一个含椭圆孔的岩石试件在不同尺寸下的应力场, 并与相应的无限域椭圆孔问题的解析解进行了比较. 数值分析结果表明, 含椭圆孔的岩石试件在单轴压缩情况下, 其边界效应对内部应力场的影响不容忽视; 试件的临界尺寸长宽比约为3:2, 孔的尖端到边界的距离与孔的长轴的比值约为7.3. 这种数值分析试件尺寸的方法, 不仅可以为岩石力学实验提供理论依据, 同时也为确定数值模拟断层错动等问题中几何模型的大小提供了参考价值. 相似文献
17.
--We present imaging results in geophysics based on using multidimensional Gaussian wavelets as a filter in a 2-D Cartesian domain. Besides decomposing the field into various distinct lengthscales, we have also constructed the 2-D maps describing the spatial distributions of the maximum of the wavelet-transformed L2-norm Emax (x,y) and its corresponding local wavenumber kmax (x,y), where x and y are the Cartesian coordinates. For geoid anomalies, using a wavelet filter extending to 90 degrees, we have discerned the distinct outlines of convergent and divergent tectonic zones and have conducted a quantitative comparison of the short-wavelength gravitational anomalies at those wavelengths between two different geographical locations. We have also compared the wavelet results with a nonlinear bandpass filter in the spectral domain where a Gaussian filter with the logarithm of the degree l acting as the argument has been employed. A wavelet solution, with a length-scale corresponding to 256 degrees, would need a filter with over 400 spherical harmonics centering around l=157 for an optimal spatial fit. The computational effort with the bandpass filter technique greatly exceeds those associated with wavelets. We have also shown the ability of the wavelets to analyze the vastly different scales present in high Rayleigh number convection and the mixing of passive heterogeneities driven by thermal convection. Wavelets will be a useful tool for rapid analyzing of the large multidimensional fields to be captured in many other geophysical endeavors, such as the upcoming gravity satellite missions and satellite radar interferometry images. 相似文献
18.
《水文科学杂志》2013,58(3)
Abstract Two multi-tracer tests were performed in fissured rocks accessible in underground laboratories to examine a new fluorescent dye: pyrene-1,3,6,8-tetra sulphonic acid (PTS). The first test was carried out at the Lindau Rock Laboratory (LRL), Germany, in a highly permeable ore dike, and the second, at the Grimsel Test Site (GTS), Switzerland, in a heterogeneous granite fault zone (AU 126). At the LRL new tracer was injected together with uranine in a convergent flow field (monopole test), and slightly different tracer breakthrough curves were observed according to different diffusion coefficients of both tracers. The matrix porosity calculated with the aid of the one-dimensional (1-D) single-fissure dispersion model (SFDM) agrees well with that found in earlier tracer tests and with measurements performed on core samples. At the GTS, the PTS tracer was applied together with pyranine in two-well injection–withdrawal (dipole) tests. Both tracers yielded identical tracer concentration curves, which confirm their conservative behaviour. Mathematical simulations performed with the aid of a 3-D numerical model (FRAC3DVS) yielded equally good fits for different sets of parameters, independent of whether matrix porosity was included or neglected. That lack of unique solution and the difficulty in observing the influence of matrix diffusion result from a wide distribution of the transit times of particular streamlines, which is characteristic for injection–withdrawal tests. However, both tracer tests clearly indicated that the new tracer (PTS) behaves conservatively at high pH values and can be successfully used for groundwater labelling. 相似文献
19.
第四纪松散沉积物中活断层滑动面的显微构造研究方法 总被引：3，自引：0，他引：3
第四纪松散深积物中的活断层滑动面是重要的古地震遗迹，也是鉴定古地震事件和确定古地震时代的重要标志。以往，以第四纪松散沉积物中活断层滑动面主要借助宏观观察进行研究，但遇到了一些难以解决的问题，如隐性断层问题。作者利用显微构造分析，找到了研究这些问题的途径，其方法是对第四纪松散沉积物进行定向原状样品采集，通过室内样品固结，显微薄片制作和显微构造现象观察等，从显微域里研究各种变形现象。用此方法，可以在没有宏观资料支持的情况下确定断面的运动性质；发现并描述第四纪松散沉积物中的隐性断层，研究断层的尖来机制，还可以对工程建设场地发现的不同松散沉积物相接触的现象进行活断层鉴别。 相似文献
20.
Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based on the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods. 相似文献