A Multilevel Method for the Solution of Time Dependent Optimal Transport

In this paper we present a new computationally efficient numerical scheme for the minimizing flow for the computation of the optimal $L_2$mass transport mapping using the fluid approach. We review the method and discuss its numerical properties. We then derive a new scaleable, efficient discretization and a solution technique for the problem and show that the problem is equivalent to a mixed form formulation of a nonlinear fluid flow in porous media. We demonstrate the effectiveness of our approach using a number of numerical experiments.     

Discontinuous Galerkin relaxation algorithm for solute transport in porous media

We consider a degenerate parabolic convection dominated equation which models the transport of contaminant in porous media. The numerical scheme is fulfilled by combining the DG – Discontinuous Galerkin Method with and an efficient relaxation algorithm that recently developed. Numerical results show the efficiency of our scheme. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence of finite volume schemes for triangular systems of conservation laws

We consider non-strictly hyperbolic systems of conservation laws in triangular form, which arise in applications like three-phase flows in porous media. We device simple and efficient finite volume schemes of Godunov type for these systems that exploit the triangular structure. We prove that the finite volume schemes converge to weak solutions as the discretization parameters tend to zero. Some numerical examples are presented, one of which is related to flows in porous media. The research of K. H. Karlsen was supported by an Outstanding Young Investigators Award from the Research Council of Norway.     

Two‐phase flows in karstic geometry

Multiphase flow phenomena are ubiquitous. Common examples include coupled atmosphere and ocean system (air and water), oil reservoir (water, oil, and gas), and cloud and fog (water vapor, water, and air). Multiphase flows also play an important role in many engineering and environmental science applications. In some applications such as flows in unconfined karst aquifers, karst oil reservoir, proton membrane exchange fuel cell, multiphase flows in conduits, and in porous media must be considered together. Geometric configurations that contain both conduit (or vug) and porous media are termed karstic geometry. Despite the importance of the subject, little work has been performed on multiphase flows in karstic geometry. In this paper, we present a family of phase–field (diffusive interface) models for two‐phase flow in karstic geometry. These models together with the associated interface boundary conditions are derived utilizing Onsager's extremum principle. The models derived enjoy physically important energy laws. A uniquely solvable numerical scheme that preserves the associated energy law is presented as well. Copyright © 2013 John Wiley & Sons, Ltd.     

Laminar-turbulent reactive flows in porous media

We state some recent results concerning liquid-vapor phase transitions for a fluid flow through a porousmedium. The focus is on the friction exerted by the porous medium, which is modeled in such a way to include both laminar and turbulent flows. In this way we obtain a hyperbolic system of three balance laws with a forcing term that is discontinuous in the state variables. Existence, uniqueness and qualitative behavior of traveling waves is proved by a novel regularization technique.     

不可压饱和多孔弹性杆动力响应的多辛方法

研究了不可压饱和多孔弹性杆的一维动力响应问题.基于多孔介质理论,在流相和固相微观不可压、固相骨架小变形的假定下,建立了不可压流体饱和多孔弹性杆一维轴向动力响应的数学模型.利用Hamilton空间体系的多辛理论,构造了不可压饱和多孔弹性杆轴向振动方程的多辛形式及其多种局部守恒律.采用中点Box离散方法得到轴向振动方程的多辛离散格式和局部能量守恒律以及局部动量守恒律的离散格式;数值模拟了不可压饱和多孔弹性杆的轴向振动过程,记录了每一时间步的局部能量数值误差和局部动量数值误差.结果表明,已构造的多辛离散格式具有很高的精确性和较长时间的数值稳定性,这为解决饱和多孔介质的动力响应问题提供了新的途径.     

Une méthode particulaire déterministe pour des équations diffusives non linéaires

We study in this Note a deterministic particle method for heat (or Fokker–Planck) equations or for porous media equations. This method is based upon an approximation of these equations by nonlinear transport equations and we prove the convergence of that approximation. Finally, we present some numerical experiments for the heat equation and for an example of porous media equations.     

Non-isothermal,compressible gas flow for the simulation of an enhanced gas recovery application

In this work, we present a framework for numerical modeling of CO₂ injection into porous media for enhanced gas recovery (EGR) from depleted reservoirs. Physically, we have to deal with non-isothermal, compressible gas flows resulting in a system of coupled non-linear PDEs. We describe the mathematical framework for the underlying balance equations as well as the equations of state for mixing gases. We use an object-oriented finite element method implemented in C++. The numerical model has been tested against an analytical solution for a simplified problem and then applied to CO₂ injection into a real reservoir. Numerical modeling allows to investigate physical phenomena and to predict reservoir pressures as well as temperatures depending on injection scenarios and is therefore a useful tool for applied numerical analysis.     

Study of a Numerical Approach for a Transient Flow Problem in Porous Media

In this work, we deal with the numerical study of the new approximation method proposed in [7] for a transient flow problem in porous media. The stationary problem, obtained from a time discretization of this transient problem, is considered as an optimal shape design formulation. We prove the existence of the solution of the discrete optimal shape problem obtained from finite element discretization. We study the convergence and give numerical results showing the efficiency of the proposed approach.     

A multiphase single relaxation time lattice Boltzmann model for heterogeneous porous media

The lattice Boltzmann (LB) method has been shown to be a highly efficient numerical method for solving fluid flow in confined domains such as pipes, irregularly shaped channels or porous media. Traditionally the LB method has been applied to flow in void regions (pores) and no flow in solid regions. However, in a number of scenarios, this may not suffice. That is partial flow may occur in semi-porous regions. Recently gray-scale LB methods have been applied to model single phase flow in such semi-porous materials. Voxels are no longer completely void or completely solid but somewhere in between. We extend the single relaxation time LB method to model multiphase, immiscible flow (e.g., gas and liquid or water and oil) in a semi-porous medium. We compare the solution to test cases and find good agreement of the model as compared to analytical solutions. We then apply the model to real porous media and recover both capillary and viscous flow regimes. However, some deficiencies in the single relaxation time LB method applied to multiphase flow are uncovered and we describe methods to overcome these limitations.     

A Hierarchy of Models for Two-Phase Flows

Summary. We derive a hierarchy of models for gas-liquid two-phase flows in the limit of infinite density ratio, when the liquid is assumed to be incompressible. The starting model is a system of nonconservative conservation laws with relaxation. At first order in the density ratio, we get a simplified system with viscosity, while at the limit we obtain a system of two conservation laws, the system of pressureless gases with constraint and undetermined pressure. Formal properties of this constraint model are provided, and sticky blocks solutions are introduced. We propose numerical methods for this last model, and the results are compared with the two previous models. Received April 20, 2000; accepted September 12, 2000 %%Online publication November 15, 2000 Communicated by Gérard Iooss     

Minimization of grid orientation effects in simulation of oil recovery processes through use of an unsplit higher order scheme

We present an unsplit second-order finite difference algorithm for hyperbolic conservation laws in several variables. Although the method can be directly implemented for general hyperbolic systems, we focus in this article on reducing grid orientation effects in porous media flow. In particular, we consider miscible and immiscible displacement processes. Our main concern is to develop a scheme that can easily be implemented into existing standard finite-difference-based reservoir simulators as an option to be used if grid orientation effects occur. The principle of the scheme is to build a higher order scheme to reduce numerical dispersion and that does not split the spatial operator to reduce the effect of the grid orientation. Numerical results are presented, which show that the method presented here reduces the effect of the numerical dispersion to a level that minimizes the grid orientation effects in a computationally efficient manner. © 1996 John Wiley & Sons, Inc.     

Variational Phase Field Modeling of Fracture Caused by Fluid Transport in Poroelastic Solids

The numerical modeling of failure mechanisms due to fracture based on sharp crack discontinuities is extremely demanding and suffers in situations with complex crack topologies. This drawback can be overcome by recently developed diffusive crack modeling concepts, which are based on the introduction of a crack phase field. Such an approach is conceptually in line with gradient-extended continuum damage models which include internal length scales. In this paper, we extend our recently outlined mechanical framework [1–3] towards the phase field modeling of fracture in the coupled problem of fluid transport in deforming porous media. Here, extremely complex crack patterns may occur due to drying or hydraulic induced fracture, the so called fracking. We develop new variational potentials for Biot-type fluid transport in porous media at finite deformations coupled with phase field fracture. It is shown, that this complex coupled multi-field problem is related to an intrinsic mixed variational principle for the evolution problem. This principle determines the rates of deformation, fracture phase field and fluid content along with the fluid potential. We develop a robust computational implementation of the coupled problem based on the potentials mentioned above and demonstrate its performance by the numerical simulation of complex fracture patterns. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Propagation of nonlinear thermoelastic waves in porous media within the theory of heat conduction with memory: physical derivation and exact solutions

In this paper, we study a mathematical model of nonlinear thermoelastic wave propagation in fluid‐saturated porous media, considering memory effect in the heat propagation. In particular, we derive the governing equations in one dimension by using the Gurtin–Pipkin theory of heat flux history model and specializing the relaxation function in such a way to obtain a fractional Erdélyi–Kober integral. In this way, we obtain a nonlinear model in the framework of time‐fractional thermoelasticity, and we find an explicit analytical solution by means of the invariant subspace method. A second memory effect that can play a significant role in this class of models is parametrized by a generalized time‐fractional Darcy law. We study the equations obtained also in this case and find an explicit traveling wave type solution. Copyright © 2016 John Wiley & Sons, Ltd.     

多孔介质中两相不可压缩可混溶驱动问题的动态混合元方法及其特征修正格式

许多依赖时间的问题涉及到局部化现象,如突出的前沿位置、激波、边界层等, 其位置随时间而变动.多孔介质中两相不可压缩可混溶驱动问题是一典型的、有代表性 的"局部化现象"问题,其数学模型为耦合非线性偏微分方程组的初边值问题.为减轻数 值解在局部前沿位置的数值振荡,提高解的精确性,本文给出了该问题的动态混合元格 式和沿特征线修正的动态混合元格式,证明了其收敛性,并给出了误差估计.     

Extended finite element modeling of deformable porous media with arbitrary interfaces

In this paper, an enriched finite element method is presented for numerical simulation of saturated porous media. The arbitrary discontinuities, such as material interfaces, are encountered via the extended finite element method (X-FEM) by enhancing the standard FEM displacements. The X-FEM technique is applied to the governing equations of porous media for the spatial discretization, followed by a generalized Newmark scheme used for the time domain discretization. In X-FEM, the material interfaces are represented independently of element boundaries and the process is accomplished by partitioning the domain with some triangular sub-elements whose Gauss points are used for integration of the domain of elements. Finally, several numerical examples are analyzed, including the dynamic analysis of the failure of lower San Fernando dam, to demonstrate the efficiency of the X-FEM technique in saturated porous soils.     

Convergence of relaxation schemes for conservation laws

We study the stability and the convergence for a class of relaxing numerical schemes for conservation laws. Following the approach recently proposed by S. Jin and Z. Xin we use a semilinear local relaxation approximation, with a stiff lower order term, and we construct some numerical first and second order accurate algorithms, which are uniformly bounded in the L^∞ and BV norms with respect to the relaxation parameter. The relaxation limit is also investigated.     

Generalization of a Relaxation Scheme for Systems of Forced Nonlinear Hyperbolic Conservation Laws with Spatially Dependent Flux Functions

A generalization of a finite difference method for calculating numerical solutions to systems of nonlinear hyperbolic conservation laws in one spatial variable is investigated. A previously developed numerical technique called the relaxation method is modified from its initial application to solve initial value problems for systems of nonlinear hyperbolic conservation laws. The relaxation method is generalized in three ways herein to include problems involving any combination of the following factors: systems of nonlinear hyperbolic conservation laws with spatially dependent flux functions, nonzero forcing terms, and correctly posed boundary values. An initial value problem for the forced inviscid Burgers' equation is used as an example to show excellent agreement between theoretical solutions and numerical calculations. An initial boundary value problem consisting of a system of four partial differential equations based on the two-layer shallow-water equations is solved numerically to display a more general applicability of the method than was previously known.     

A multiscale Darcy–Brinkman model for fluid flow in fractured porous media

The aim of this work is to present a reduced mathematical model for describing fluid flow in porous media featuring open channels or fractures. The Darcy’s law is assumed in the porous domain while the Stokes–Brinkman equations are considered in the fractures. We address the case of fractures whose thickness is very small compared to the characteristic diameter of the computational domain, and describe the fracture as if it were an interface between porous regions. We derive the corresponding interface model governing the fluid flow in the fracture and in the porous media, and establish the well-posedness of the coupled problem. Further, we introduce a finite element scheme for the approximation of the coupled problem, and discuss solution strategies. We conclude by showing the numerical results related to several test cases and compare the accuracy of the reduced model compared with the non-reduced one.     

Entropy solution theory for fractional degenerate convection–diffusion equations

We study a class of degenerate convection-diffusion equations with a fractional non-linear diffusion term. This class is a new, but natural, generalization of local degenerate convection-diffusion equations, and include anomalous diffusion equations, fractional conservation laws, fractional porous medium equations, and new fractional degenerate equations as special cases. We define weak entropy solutions and prove well-posedness under weak regularity assumptions on the solutions, e.g. uniqueness is obtained in the class of bounded integrable solutions. Then we introduce a new monotone conservative numerical scheme and prove convergence toward the entropy solution in the class of bounded integrable BV functions. The well-posedness results are then extended to non-local terms based on general Lévy operators, connections to some fully non-linear HJB equations are established, and finally, some numerical experiments are included to give the reader an idea about the qualitative behavior of solutions of these new equations.