Weighted average flux method and flux limiters for the numerical simulation of shock waves in rigid porous media

The one‐dimensional flow field generated by the passage of a shock wave in a rigid, thermoelastic porous foam has been simulated using a two‐phase mathematical model. The work presented here makes use of the weighted average flux method to solve the system of six equations that govern the problem. Spurious oscillations are eliminated through the application of total variation diminishing limiting methods. Four different limiters were tested: van Leer, SuperA, MinA and van Albada. Numerical tests were carried out to verify the performance of each flux limiter in terms of accuracy. The results were compared to analytical and previously obtained data to assess the performance of the mathematical model. Excellent agreement was obtained. Copyright © 2002 John Wiley & Sons, Ltd.     

Horizontal boundary-layer natural convection in a porous medium saturated with a gas

Boundary-layer analysis is performed for free convection flow over a hot horizontal surface embedded in a porous medium saturated with a gas of variable properties. The variable gas properties are accounted for via the assumption that thermal conductivity and dynamic viscosity are proportional to temperature. A similarity solution is shown to exist for the case of constant surface temperature. Numerical results for the stream function, horizontal velocity, and temperature profiles within the boundary layer as well as for the mass of entrained gas, surface slip velocity, and heat transfer rate at different values of the wall-temperature parameter are presented. Asymptotic solutions for large heating are also available to support the numerical work.     

Interaction of evaporation fronts with a formation interface in a porous medium

The evolutionarity and structure of water-vapor phase discontinuities formed in a geothermal reservoir on the interface between permeable formations with different properties are considered. In the short-wave approximation a graphic method is proposed for solving the problem of breakdown of an arbitrary discontinuity in a geothermal reservoir consisting of two formations with different properties.     

Transient non-Darcy free convection between parallel vertical plates in a fluid saturated porous medium

Transient non-Darcy free convection between two parallel vertical plates in a fluid saturated porous medium is investigated using the generalized momentum equation proposed by Vafai and Tien. The effects of porous inertia and solid boundary are considered in addition to the Darcy flow resistance. Exact solutions are found for the asymptotic states at small and large times. The large time solutions reveal that the velocity profiles are rather sensitive to the Darcy number Da when Da<1. It has also been found that boundary friction alters the velocity distribution near the wall, considerably. Finite difference calculations have also been carried out to investigate the transient behaviour at the intermediate times in which no similarity solutions are possible. This analytical and numerical study reveals that the transient free convection between the parallel plates may well be described by matching the two distinct asymptotic solutions obtained at small and large times.Nomenclature C empirical constant for the Forchheimer term - f velocity function for the small time solution - F velocity function for the large time solution - g acceleration due to gravity - Gr* micro-scale Grashof number - H a half distance between two infinite plates - K permeability - Nu Nusselt number - Pr Prandtl number - t time - T temperature - u, v Darcian velocity components - x, y Cartesian coordinates - effective thermal diffusivity - coefficient of thermal expansion - porosity - dimensionless time - similarity variable - dimensionless temperature - viscosity - kinematic viscosity - density - the ratio of heat capacities     

Studying the Darcy-Stefan problem on phase transition in a saturated porous soil

The Cauchy problem for the Darcy-Stefan model, which describes the process of freezing (thawing) of a saturated porous soil with allowance for liquid-phase filtration, is considered. The model includes the Darcy law, the equation of liquid-phase incompressibility, the equation of absence of solid-phase motion, the equation of energy balance in the porous soil-saturating continuous medium system, and also the Stefan condition and the condition of continuity of the normal components of the velocity field at the interface boundary. The existence of generalized solutions of the problem satisfying an additional condition of entropy nondecreasing in a thermomechanical system (i.e., the second law of thermodynamics) is proved by the method of the kinetic equation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 4, pp. 81–93, July–August, 2008.     

Compressibility and anisotropy effects in the processes of two-phase flow through a porous medium

Functions of the relative phase permeabilities which take into account the tensor nature of the relations between the absolute and phase permeabilities are constructed on the basis of experimental data. Problems of the displacement of gas by water and of water by gas in an orthotropic pore space are studied. The relative phase permeabilities obtained are compared with functions that depend only on the local saturation. The effect of the phase compressibility and the anisotropy of the absolute permeability tensor is considered.     

A model for solute transport following an abrupt pressure impact in saturated porous media

A mathematical model is developed of an abrupt pressure impact applied to a compressible fluid with solute, flowing through saturated porous media. Nondimensional forms of the macroscopic balance equations of the solute mass and of the fluid mass and momentum lead to dominant forms of these equations. Following the onset of the pressure change, we focus on a sequence of the first two time intervals at which we obtain reduced forms of the balance equations. At the very first time period, pressure is proven to be distributed uniformly within the affected domain, while solute remains unaffected. During the second time period, the momentum balance equation for the fluid conforms to a wave form, while the solute mass balance equation conforms to an equation of advective transport. Fluid's nonlinear wave equation together with its mass balance equation, are separately solved for pressure and velocity. These are then used for the solution of solute's advective transport equation. The 1-D case, conforms to a pressure wave equation, for the solution of fluid's pressure and velocity. A 1-D analytical solution of the transport problem, associates these pressure and velocity with an exponential power which governs solute's motion along its path line.     

A numerical study of fluids with pressure‐dependent viscosity flowing through a rigid porous medium

Much of the work on flow through porous media, especially with regard to studies on the flow of oil, are based on ‘Darcy's law’ or modifications to it, such as Darcy–Forchheimer or Brinkman models. While many theoretical and numerical studies concerning flow through porous media have taken into account the inhomogeneity and anisotropy of the porous solid, they have not taken into account the fact that the viscosity of the fluid and drag coefficient could depend on the pressure in applications, such as enhanced oil recovery (EOR). Experiments clearly indicate that the viscosity varies exponentially with respect to the pressure and the viscosity can change, in some applications, by several orders of magnitude. The fact that the viscosity depends on pressure immediately implies that the ‘drag coefficient’ would also depend on the pressure. In this paper we consider modifications to Darcy's equation wherein the drag coefficient is a function of pressure, which is a realistic model for technological applications, such as EOR and geological carbon sequestration. We first outline the approximations behind Darcy's equation and the modifications that we propose to Darcy's equation, and derive the governing equations through a systematic approach using mixture theory. We then propose a stabilized mixed finite element formulation for the modified Darcy's equation. To solve the resulting nonlinear equations we present a solution procedure based on the consistent Newton–Raphson method. We solve representative test problems to illustrate the performance of the proposed stabilized formulation. One of the objectives of this paper is to show that the dependence of viscosity on the pressure can have a significant effect both on the qualitative and quantitative nature of the solution. Copyright © 2010 John Wiley & Sons, Ltd.     

A boundary layer analysis of combined heat and mass transfer by natural convection from a concentrated source in a saturated porous medium

A boundary layer analysis was carried out to investigate the coupled phenomena of heat and mass transfer by natural convection from concentrated heat and mass sources embedded in saturated porous media. Both line and point source problems were treated. The boundary layer equations based on Darcy's law and Boussinesq approximation were solved by means of similarity transformation to obtain the details of velocity, temperature and concentration distributions above a concentrated heat source. Two important parameters, namely the Lewis number Le and the buoyancy ratioN were identified to conduct a series of numerical integrations. For the case of small Le, a substance diffuses further away from the plume centerline, such that the mass transfer influences both velocity and temperature profiles over a wide range. For large Le, on the other hand, the substance diffuses within a narrow range along the centerline. Naturally, the influence of mass transfer is limited to the level of the centerline velocity, so that a peaky velocity profile appears for positiveN whereas a velocity defect emerges along the centerline for negativeN. For such cases of large Le, the temperature profiles are found to be fairly insensitive to Le.     

Plane steady-state incompressible fluid flow through a porous medium with a nonlinear resistance laws in the presence of a sink

Plane nonlinear fluid flows through a porous medium which simulate a sink located at the same distance from the roof and floor of the stratum for two nonlinear flow laws are constructed. The following flow laws are taken: a power law and a law of special form reducing to analytic functions in the hodograph plane.     

A new numerical scheme for convective dominated heat transfer in a porous medium with strong temperature gradients

A new numerical scheme, theimplicit correction scheme, has been developed for heat transfer in a porous medium with strong temperature gradients. The scheme includes diffusion, convection and transverse heat transfer processes. By using correction coefficients which are based on transverse heat transfer, the effects of convection enthalpy flow and diffusion are modified. Under suitable limiting conditions, the implicit correction scheme can be reduced to the central-difference, upwind, or power-law scheme. The correction scheme is shown to be especially useful in calculations of the thermal effectiveness of the regenerator in Stirling cycle refrigeration.     

Numerical investigation of the first transition in the three-dimensional problem of convective flow in a porous medium

The branching off of steady-state regimes from mechanical equilibrium is studied for the problem of filtration convection in a parallelepiped. The conditions for the geometric parameters under which stable continuous families of steady-state regimes develop are found. The stability of equilibria of the family with respect to three-dimensional perturbations is analyzed in a numerical experiment using a finite-difference method.     

剧变截面圆管内渗流的数值计算方法

对于剧变截面圆管的渗流问题写出不可压缩渗流的基本方程组,对直接求解原始变量(速度和压力)的数值计算方法作出改进。先由非主流方向的运动方程计算压力,后由主流方向的运动方程计算主流方向的速度分量,再由连续性方程计算非主流方向的速度分量。这样可以避免在一般的求解原始变量方法中由连续性方程计算压力时出现的困难和麻烦。根据本方法和剧变截面圆管的特点,采用半交错不等距非正交贴体混合网格系。本文详细写出差分方程和迭代计算公式,对剧变截面圆管内的渗流算例进行数值计算。本方法的优点是简单和实用,在工程上具有较大的应用价值。     

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability简

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.     

覆水饱和双相介质圆弧场地P波激励下散射问题的理论解

采用波函数展开法对平面P波入射复杂水域地形的空间变异性地震动场进行研究,该水域地形具有覆水层、饱和双相介质、场地非平坦以及第二类分层(场地跨越分层界面)等属性.首先,依据地震波反射和透射特性推导直角坐标系下的自由波场分布;然后,根据场地属性并引入大圆弧法分析极坐标系下的含有待定系数的散射波场;进而,结合土-水分界面和饱和土层分界面边界条件,求解散射波场中的待定系数;最后,通过自由波场和散射波场得到覆水饱和双相介质圆弧场地波函数理论解.基于理论解,通过算例验证了理论推导的合理性及可靠性,分析了地表位移在不同入射条件下的差异性.结果 表明,相对于均匀介质,饱和双相介质会显著影响地表位移分布.此外,入射波频率和角度对地震地面运动特性也有较大的影响.     

流体饱和多孔介质动力问题的显式时域解法

为了简化分析,Zienkiewicz等基于Biot理论,在忽略流体相对于土骨架运动的加速度条件下,建立了以土骨架位移u和孔隙流体压力p为基本变量的u-p格式饱和两相介质动力方程。针对该u-p方程,在空间上,采用伽辽金法有限元离散,并结合对角化形式的质量矩阵和流体压缩矩阵,忽略相邻结点间的惯性和流体压缩量间的耦合作用。在时域内,基于杜修力等提出的显式算法和Euler预估-校正法,建立了一种具有二阶精度的全显式时域积分法。采用一维饱和土模型,对比提出算法的数值解与Simon方法的解析解,发现两者吻合良好,验证了本文方法的正确性。并分析了饱和土二维动力问题,以及渗透系数和排水条件对饱和土动力响应的影响。