Simulation of Single-Phase Multicomponent Flow Problems in Gas Reservoirs by Eulerian–Lagrangian Techniques

Over the past two decades most discussions of the simulation of miscible displacement in porous media were related to incompressible flow problems; recently, however, attention has shifted to compressible problems. The first goal of this paper is the derivation of the governing equations (mathematical models) for a hierarchy of miscible isothermal displacements in porous media, starting from a very general single-phase, multicomponent, compressible flow problem; these models are then compared with previously proposed models. Next, we formulate an extension of the modified method of characteristics with adjusted advection to treat the transport and dispersion of the components of the miscible fluid; the fluid displacement must be coupled in a two-stage operator-splitting procedure with a pressure equation to define the Darcy velocity field required for transport and dispersion, with the outer stage incorporating an implicit solution of the nonlinear parabolic pressure equation and an inner stage for transport and diffussion in which the mass fraction equations are solved sequentially by first applying a globally conservative Eulerian–Lagrangian scheme to solve for transport, followed by a standard implicit procedure for including the diffusive effects. The third objective is a careful investigation of the underlying physics in compressible displacements in porous media through several high resolution numerical experiments. We consider real binary gas mixtures, with realistic thermodynamic correlations, in homogeneous and heterogeneous formations.     

An Analytical Solution to the One-Dimensional Solute Advection-Dispersion Equation in Multi-Layer Porous Media

An analytical solution to the one-dimensional solute advection-dispersion equation in multi-layer porous media is derived using a generalized integral transform method. The solution was derived under conditions of steady-state flow and arbitrary initial and inlet boundary conditions. The results obtained by this solution agree well with the results obtained by numerically inverting Laplace transform-generated solutions previously published in the literature. The analytical solution presented in this paper provides more flexibility with regard to the inlet conditions. The numerical evaluation of eigenvalues and matrix exponentials required in this solution technique can be accurately and efficiently computed using the sign-count method and eigenvalue evaluation methods commonly available. The illustrative calculations presented herein have shown how an analytical solution can provide insight into contaminant distribution and breakthrough in transport through well defined layered column systems. We also note that the method described here is readily adaptable to two and three-dimensional transport problems.     

Variational Lagrangian-thermodynamics of nonisothermal finite strain mechanics of porous solids and thermomolecular diffusion

The nonisothermal finite strain dynamics of a porous solid containing a viscous fluid is developed on the basis of a new thermodynamics of open systems and irreversible processes. The same theory is applicable to the mechanics of a nonporous solid with thermomolecular diffusion of a substance in solution. New fundamental concepts of “thermobaric” and “convective” potentials are presented in the context of porous solids. Field equations and Lagrangian equations with generalized coordinates are derived directly from a variational principle of “virtual dissipation”. Inclusion of nonlinear viscoelasticity and plastic behavior is indicated. Partial saturation of pore fluid is discussed. The theory is applicable to the mechanics of a non porous solid with thermolecular diffusion of several molecular species in solution, and under certain conditions to the analogous case of a porous solid containg a fluid mixture. It is shown how the Lagrangian equations provide the foundation of finite element methods.     

Predicting transport by Lagrangian coherent structures with a high-order method

Recent developments in identifying Lagrangian coherent structures from finite-time velocity data have provided a theoretical basis for understanding chaotic transport in general flows with aperiodic dependence on time. As these theoretical developments are extended and applied to more complex flows, an accurate and general numerical method for computing these structures is needed to exploit these ideas for engineering applications. We present an unstructured high-order hp/spectral-element method for solving the two-dimensional compressible form of the Navier–Stokes equations. A corresponding high-order particle tracking method is also developed for extracting the Lagrangian coherent structures from the numerically computed velocity fields. Two different techniques are used; the first computes the direct Lyapunov exponent from an unstructured initial particle distribution, providing easier resolution of structures located close to physical boundaries, whereas the second advects a small material line initialized close to a Lagrangian saddle point to delineate these structures. We demonstrate our algorithm on simulations of a bluff-body flow at a Reynolds number of Re = 150 and a Mach number of M = 0.2 with and without flow forcing. We show that, in the unforced flow, periodic vortex shedding is predicted by our numerical simulations that is in stark contrast to the aperiodic flow field in the case with forcing. An analysis of the Lagrangian structures reveals a transport barrier that inhibits cross-wake transport in the unforced flow. The transport barrier is broken with forcing, producing enhanced transport properties by chaotic advection and consequently improved mixing of advected scalars within the wake.     

Lagrangian Gas Dynamics in Two Dimensions and Lagrangian systems

We propose a new and canonical way of writing the equations of gas dynamics in Lagrangian coordinates in two dimensions as a weakly hyperbolic system of conservation laws. One part of the system is called the physical part and contains physical variables; the other part is the geometrical part. We show that the physical part is symmetrizable. We show that the weak hyperbolicity is due to shear contact discontinuities. Free divergence constraints play an important role in the system. We prove the L² stability of the physical part of the system. Based on this formulation, we derive a new conservative and entropy-consistent finite-volume numerical scheme. We prove the stability of the numerical scheme. Numerical results show the potential interest of this approach. Various examples (Born-Infeld, MHD, 3D lagrangian gas dynamics) can be written using the same abstract formalism.     

饱和-非饱和土壤中污染物运移过程的数值模拟

本文提出了一个模拟饱和 非饱和土壤中溶和污染物运移过程的数值模型．模拟的控制污染物运移的物理 化学现象包括：对流，机械逸散，分子弥散，吸附，蜕变，不动水效应．发展了一个修正的特征线Ｇａｌｅｒｋｉｎ方法以离散污染物运移过程的控制方程并导出了一个用于有限元方程求解的显式算法．数值例题结果表明所提出模型和算法的功能     

Torsion of cylindrical bodies with varying strain properties

We analyze the results of experimental studies of effective strain properties of damaged, porous, and other inhomogeneous materials and study the main laws of their behavior under strain. We consider the possible versions of constitutive relations taking account of the dependence of the properties of the media under study on the loading conditions or the strain conditions and the relations between the shear and bulk strains. Since the traditional statement of the torsion problems for bodies with such properties cannot be used, we analyze the strain consistency equations and the relations between the strains and displacements in cylindrical coordinates and obtain expressions for the displacements in an appropriate generalized form, which can be used not only for the torsion problems. We study how the distributions of displacements, strains, and stresses under torsion depend on the parameter characterizing the susceptibility of the material strain properties to variations in the stress state type. We show that, in the case of torsion of a cylinder of circular cross-section, there is no deplanation of the cross-section, just as in the classical solution, but the distributions of displacements, strains, and stresses significantly differ from the well-known solutions.     

<Emphasis Type="Italic">A Priori</Emphasis> Parameter Estimation for the Thermodynamically Constrained Averaging Theory: Species Transport in a Saturated Porous Medium

The thermodynamically constrained averaging theory (TCAT) has been used to develop a simplified entropy inequality (SEI) for several major classes of macroscale porous medium models in previous works. These expressions can be used to formulate hierarchies of models of varying sophistication and fidelity. A limitation of the TCAT approach is that the determination of model parameters has not been addressed other than the guidance that an inverse problem must be solved. In this work we show how a previously derived SEI for single-fluid-phase flow and transport in a porous medium system can be reduced for the specific instance of diffusion in a dilute system to guide model closure. We further show how the parameter in this closure relation can be reliably predicted, adapting a Green’s function approach used in the method of volume averaging. Parameters are estimated for a variety of both isotropic and anisotropic media based upon a specified microscale structure. The direct parameter evaluation method is verified by comparing to direct numerical simulation over a unit cell at the microscale. This extension of TCAT constitutes a useful advancement for certain classes of problems amenable to this estimation approach.     

Variational irreversible thermodynamics of physical-chemical solids with finite deformation

A new thermodynamics of open thermochemical systems and a variational principle of virtual dissipation are applied to the finite deformation of a solid coupled to thermomolecular diffusion and chemical reactions. A variational derivation is obtained of the field differential equations as well as Lagrangian equations with generalized coordinates. New formulas for the affinity and a new definition of the chemical potential are presented. An outline is given of an unusually large field of applications, such as active transport in biological systems, finite element methods, plastic properties as analogous to chemical reactions, phase changes and recrystalization, porous solids, heredity and initially stressed solids. A new and unified insight is thus provided in highly diversified problems.     

Stochastic Flow Simulation and Particle Transport in a 2D Layer of Random Porous Medium

A stochastic numerical method is developed for simulation of flows and particle transport in a 2D layer of porous medium. The hydraulic conductivity is assumed to be a random field of a given statistical structure, the flow is modeled in the layer with prescribed boundary conditions. Numerical experiments are carried out by solving the Darcy equation for each sample of the hydraulic conductivity by a direct solver for sparse matrices, and tracking Lagrangian trajectories in the simulated flow. We present and analyze different Eulerian and Lagrangian statistical characteristics of the flow such as transverse and longitudinal velocity correlation functions, longitudinal dispersion coefficient, and the mean displacement of Lagrangian trajectories. We discuss the effect of long-range correlations of the longitudinal velocities which we have found in our numerical simulations. The related anomalous diffusion is also analyzed.     

Factorization of Transport Coefficients in Macroporous Media

We prove the fundamental theorem about factorization of the phenomenological coefficients for transport in macroporous media. By factorization we mean the representation of the transport coefficients as products of geometric parameters of the porous medium and the parameters characteristic of the multicomponent fluid saturating the porous space. The two permeabilities of the porous medium, the convective and the diffusional ones, are separated. A similarity between the diffusional permeability and the porosity–tortuosity factor of the Kozeny–Carman theory is demonstrated. We do not make any specific assumption about stochastic or deterministic structure of the porous medium. The fluxes in fluid on the pore level are described by general relations of the non-equilibrium thermodynamics.     

二维柱坐标系中守恒的保球对称的中心型拉氏方法

提出了一种Godunov型中心型拉氏方法,用于求解二维柱坐标系中的可压缩多介质Euler方程组,该方法完全在体积控制体上离散,不仅保证质量、动量和总能量守恒,且该方法在二维柱坐标系中保一维球对称;并且对一维球对称问题在球对称网格划分下,精度测试表明该方法具有一阶精度,算例显示方法非常有效。     

各向异性展状介质中弹性波的传播矩阵解法

基于广义胡克定律及混和变量弹性波方程，解析求得各层介质内位移、应力传递矩阵，给出了直角坐标系下各向异性层状介质中弹性波的传播矩阵解法．该方法适用于非轴对称各向异性和点源作用，较好地解决了数值计算中有效数字精度损失问题．数值结果表明，计算效率、准确性及稳定性均较好．     

ONSET CONDITION OF STRAIN LOCALIZATION IN MATRIX OF SATURATED POROUS MEDIA

Introduction Strainlocalizationofgeomaterialsisoneofmostpopularfailuretypesinnature,which canbeshowedaslandslidesandmudflowsinmountainousareasunderincessantorheavy raining,especiallythevegetationisseverelydamagedbywoodsharvest;pipingeffect,a typeoflocalfa…     

A symplectic integrator for dynamic coupling between nonlinear vessel motion with variable cross-section and bottom topography and interior shallow-water sloshing

The coupled motion between shallow-water sloshing in a moving vessel with variable cross-section and bottom topography, and the vessel dynamics is considered, with the vessel dynamics restricted to horizontal motion governed by a nonlinear spring. The coupled fluid and vessel equations in Eulerian coordinates are transformed to the Lagrangian particle path setting which leads to a formulation with nice properties for numerical simulation. In the Lagrangian representation, a simple and fast numerical algorithm based on the Störmer–Verlet method, is implemented. The numerical scheme conserves the total energy in the system, as well as giving the partition of energy between the fluid and vessel. Numerical simulations of the coupled nonlinear dynamics are presented.     

非均匀介质中流动和传热问题的有限分析法

数值求解非均匀介质中的输运问题广泛应用于科学计算和工程领域．介质的强非均匀性给相关问题的准确求解带来极大的困难．近年来，本课题组将有限分析法拓展到该领域，建立了非均匀介质中输运问题的有限分析法．该算法基于网格奇点邻域内类拉普拉斯方程局部解析解构建，算法具有很高的精度，且不依赖于介质的非均匀性强度．不管相邻网格传导率差异如何，仅需对原始网格进行很少地细分就可以获得非常准确的计算结果，因此与其他传统数值算法相比，可以大幅提高计算精度和效率．该算法可广泛应用于求解非均匀多孔介质中的渗流、复合材料中的热传导及电场分布等问题．     

Conditional Stochastic Analysis for Multiphase Trsnsport in Randomly Heterogeneous,variably Saturated Media

Groundwater contamination of organics has recently become a problem of growing concern over the resulting health and environmental problems. In general, the multiphase system of nonaqueous phase liquid (NAPL), water and air has to be studied in order to realistically describe the movement of such materials in the subsurface. Numerous models have been developed to study multiphase flow and/or multispecies transport in porous media. However, using models to study the influence of medium heterogeneity on such flow and transport is only a recent event. It has been demonstrated for single-phase flow and transport in saturated and unsaturated media that the study of medium heterogeneity is amenable to stochastic analysis. In this paper, we extend our Eulerian–Lagrangian stochastic theory for single-phase transport to the problem of multiphase–multispecies transport in randomly heterogeneous media under the conditions that the flow is steady-state and the phases are in local chemical equilibrium. We present theoretical expressions to describe the first two conditional moments of the random concentration of any species in any phase. Though they reveal some of the fundamental properties and help gaining insight into the nature of the problem, these expressions cannot be evaluated without either high resolution Monte Carlo simulation or approximation (closure). Therefore, we propose two sets of workable approximations, one being a weak approximation and the other being a linearized pseudo-Fickian approximation. The former yields a nonlinear integro-differential equation for the first conditional moment and the latter yields a linear differential equation. Then the second moments can be computed from explicit expressions from either the weak or pseudo-Fickian approximation.