Crossover from Nonlinearity Controlled to Heterogeneity Controlled Mixing in Two-Phase Porous Media Flows

We use high resolution Monte Carlo simulations to study the dispersive mixing in two-phase, immiscible, porous media flow that results from the interaction of the nonlinearities in the flow equations with geologic heterogeneity. Our numerical experiments show that distinct dispersive regimes occur depending on the relative strength of nonlinearity and heterogeneity. In particular, for a given degree of multiscale heterogeneity, controlled by the Hurst exponent which characterizes the underlying stochastic model for the heterogeneity, linear and nonlinear flows are essentially identical in their degree of dispersion, if the heterogeneity is strong enough. As the heterogeneity weakens, the dispersion rates cross over from those of linear heterogeneous flows to those typical of nonlinear homogeneous flows.     

A multiphase approach for evaluating the horizontal and rocking impedances of pile group foundations

This paper advocates the use of a multiphase model, already developed for static or quasi‐static geotechnical engineering problems, for simulating the behaviour of piled raft foundations subject to horizontal as well as rocking dynamic solicitations. It is shown that such a model, implemented in a FEM code, yields appropriate predictions for the foundation impedance characteristics, provided that shear and bending effects in the piles are taken into account, thus corroborating the findings of the asymptotic homogenization theory. Besides, it is notably pointed out that such a multiphase‐based computational tool makes it possible to assess the dynamic behaviour of pile groups in a much quicker way than when using direct numerical simulations, which may face oversized problems when large pile groups are concerned. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical simulation of the response of a reinforced wall to a high speed train passage

The dynamic response of a mechanically stabilized earth wall to the passing of a high‐speed train is modelled using the finite element method. A three‐dimensional analysis is carried out, using a specific framework that allows performing the analysis with a moderate computational effort. In the first place, a so‐called multiphase approach is used to take into account the reinforcing strips. The moving load is taken into account by performing the calculation in a mobile referential using the properties of symmetry of the train cars and a simplifying assumption of periodicity for the whole train. We also assume a steady state. A partial validation of the approach is obtained by means of a comparison with an analytical solution. The quick increase in displacements induced by the train passing when the speed comes close to the celerity of Rayleigh waves clearly appears in the results. The vertical displacements, vertical stresses in the backfill, tensile forces in the strips and the influence of the stiffness of the soil are discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Seismic stability analysis of piled embankments: A multiphase approach

The seismic stability analysis of an embankment lying over a soft foundation soil reinforced by a group of vertical piles is performed within the framework of the upper bound kinematic approach of yield design. The analysis is based on a previously developed ‘multiphase’ model of the reinforced ground, which explicitly accounts for the shear and bending resistances of the piles. Making use of appropriate failure mechanisms involving shear zones across which the reinforcements are continuously deforming, along with ‘plastic hinge’ surfaces, upper bound estimates to the critical seismic coefficient of the structure are derived. The results, which are confirmed by the simulations obtained from a finite element elastoplastic code, give clear evidence of the key role played by the bending strength capacities of the piles in ensuring the stability of the pile reinforced embankment. Copyright © 2009 John Wiley & Sons, Ltd.     

Element mobility during continental collision: insights from polymineralic metamorphic vein within UHP eclogite in the Dabie orogen

Metamorphic dehydration and partial melting are two important processes during continental collision. They have significant bearing on element transport at the slab interface under subduction‐zone P–T conditions. Petrological and geochemical insights into the two processes are provided by a comprehensive study of leucocratic veins in ultrahigh‐pressure (UHP) metamorphic rocks. This is exemplified by this study of a polymineralic vein within phengite‐bearing UHP eclogite in the Dabie orogen. The vein is primarily composed of quartz, kyanite, epidote and phengite, with minor accessory minerals such as garnet, rutile and zircon. Primary multiphase solid inclusions occur in garnet and epidote from the both vein and host eclogite. They are composed of quartz ± K‐feldspar ± plagioclase ± K‐bearing glass and exhibit irregular to negative crystal shapes that are surrounded by weak radial cracks. This suggests their precipitation from solute‐rich metamorphic fluid/melt that involved the reaction of phengite breakdown. Zircon U–Pb dating for the vein gave two groups of concordant ages at 217 ± 2 and 210 ± 2 Ma, indicating two episodes of zircon growth in the Late Triassic. The same minerals from the two rocks give consistent δ¹⁸O and δD values, suggesting that the vein‐forming fluid was directly derived from the host UHP eclogite. The vein is much richer in phengite and epidote than the host eclogite, suggesting that the fluid is associated with remarkable concentration of such water‐soluble elements as LILE and LREE migration. Garnet and rutile in the vein exhibit much higher contents of HREE (2.2–5.7 times) and Nb–Ta (1.8–2.0 times) than those in the eclogite, indicating that these normally water‐insoluble elements became mobile and then were sunken in the vein minerals. Thus, the vein‐forming agent would be primarily composed of the UHP aqueous fluid with minor amounts of the hydrous melt, which may even become a supercritical fluid to have a capacity to transport not only LILE and LREE but also HREE and HFSE at subduction‐zone metamorphic conditions. Taken together, significant amounts of trace elements were transported by the vein‐forming fluid due to the phengite breakdown inside the UHP eclogite during exhumation of the deeply subducted continental crust.     

东川式铜矿伴生组分地质特征

东川式铜矿的成因，自否定岩浆期后热液成矿，转为以沉积—变质成矿为主流派观点以后，沉积—成岩期后成矿否定了客观存在的后期多元、多期叠加改造成矿的现象。本文由伴生组分特征研究入手，肯定了后期多元、多期叠加改造成矿的存在，希引起研究者重视，对找富铜矿及回收伴生的金、钴，均有现实的指导意义。     

Present situation of the research and application of engineering multiwave seismic prospecting in China

ＰｒｅｓｅｎｔｓｉｔｕａｔｉｏｎｏｆｔｈｅｒｅｓｅａｒｃｈａｎｄａｐｐｌｉｃａｔｉｏｎｏｆｅｎｇｉｎｅｅｒｉｎｇｍｕｌｔｉｗａｖｅｓｅｉｓｍｉｃｐｒｏｓｐｅｃｔｉｎｇｉｎＣｈｉｎａＸＩＡＯ－ＰＩＮＧＷＵ（吴小平），ＨＯＮＧ－ＲＵＺＨＡＯ（赵鸿儒）（Ｄ...     

Comparison of Galerkin-type discretization techniques for two-phase flow in heterogeneous porous media

A comparison of Standard Galerkin, Petrov-Galerkin, and Fully-Upwind Galerkin methods for the simulation of two-phase flow in heterogeneous porous media is presented. On the basis of the coupled pressure-saturation equations, a generalized formulation for all three finite element methods is derived and analysed. For flow in homogeneous media, the Petrov-Galerkin method gives excellent results. But this method fails miserably for problems with heterogeneous media. This is because it is not able to capture correctly processes that take place at interfaces when, for instance, the capillary pressure-saturation relationship after Brooks and Corey is assumed. The Fully-Upwind Galerkin method is superior to the Petrov-Galerkin approach because it is able to give correct results for flow in homogeneous and heterogeneous media for the two models of van Genuchten and Brooks-Corey. The widely used formulation which is correct for the homogeneous case cannot be used for heterogeneous media. Instead the straightforward approach of gradp_c in combination with a chord-slope technique must be utilized.     

Numerical simulation of NAPL flow in the subsurface

A three-dimensional, three-phase numerical model is presented for simulating the movement of immiscible fluids, including nonaqueous-phase liquids (NAPLs), through porous media. The model is designed to simulate soil flume experiments and for practical application to a wide variety of contamination scenarios involving light or dense NAPLs in heterogeneous subsurface systems. The model is derived for the three-phase flow of water, NAPL, and air in porous media. The basic governing equations are based upon the mass conservation of the constitutents within the phases. The descretization chosen to transform the governing equations into the approximating equations, although logically regular, is very general. The approximating equations are a set of simultaneous coupled nonlinear equations which are solved by the Newton-Raphson method. The linear system solutions needed for the Newton-Raphson method are obtained using a matrix of preconditioner/accelerator iterative methods. Because of the special way the governing equations are implemented, the model is capable of simulating many of the phenomena considered necessary for the sucessful simulation of field problems including entry pressure phenomena, entrapment, and preferential flow paths. The model is verified by comparing it with several exact analytic test solutions and three soil flume experiments involving the introduction and movement of light nonaqueous-phase liquid (LNAPL) or dense nonaqueous-phase liquid (DNAPL) in heterogeneous sand containing a watertable. This revised version was published online in August 2006 with corrections to the Cover Date.