膨胀性非牛顿流体多孔介质流动的模拟分析

在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程．采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间．对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效．分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度．当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑．压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显．     

Modelling of Power-law Fluid Flow Through Porous Media Using Smoothed Particle Hydrodynamics

The flow of non-Newtonian fluids through two-dimensional porous media is analyzed at the pore scale using the smoothed particle hydrodynamics (SPH) method. A fully explicit projection method is used to simulate incompressible flow. This study focuses on a shear-thinning power-law model (n < 1), though the method is sufficiently general to include other stress-shear rate relationships. The capabilities of the proposed method are demonstrated by analyzing a Poiseuille problem at low Reynolds numbers. Two test cases are also solved to evaluate validity of Darcy’s law for power-law fluids and to investigate the effect of anisotropy at the pore scale. Results show that the proposed algorithm can accurately simulate non-Newtonian fluid flows in porous media.     

On Optimal Designs Using Topology Optimization for Flow Through Porous Media Applications

Transport in Porous Media - Topology optimization (TopOpt) is a mathematical-driven design procedure to realize optimal material architectures. This procedure is often used to automate the design...     

Homogenization of Wall-Slip Gas Flow Through Porous Media

The permeability of reservoir rocks is most commonly measured with an atmospheric gas. Permeability is greater for a gas than for a liquid. The Klinkenberg equation gives a semi-empirical relation between the liquid and gas permeabilities. In this paper, the wall-slip gas flow problem is homogenized. This problem is described by the steady state, low velocity Navier–Stokes equations for a compressible gas with a small Knudsen number. Darcy's law with a permeability tensor equal to that of liquid flow is shown to be valid to the lowest order. The lowest order wall-slip correction is a local tensorial form of the Klinkenberg equation. The Klinkenberg permeability is a positive tensor. It is in general not symmetric, but may under some conditions, which we specify, be symmetric. Our result reduces to the Klinkenberg equation for constant viscosity gas flow in isotropic media.     

Network Modeling of Non-Darcy Flow Through Porous Media

Darcy's law is inadequate for describing high-velocity gas flow in porous media, which occurs in the near well-bore region of high capacity gas and condensate reservoirs. This study is directed at understanding the non-Darcy flow behavior. A pore-level network model has been developed to describe high velocity flow. The inputs to the model are pore size distributions and network coordination numbers. The outputs are permeability, non-Darcy coefficient, tortuousity and porosity. The additional pressure gradient term is found to be proportional to the square of the velocity in accordance with the Forchheimer's equation. The correlation between the non-Darcy coefficient and other flow properties (the permeability, the porosity and the tortuousity) is found to depend on the morphological parameters being changed. General correlations are derived between these flow properties.     

Inertia Effects in High-Rate Flow Through Heterogeneous Porous Media

The paper deals with the effects of large scale permeability–heterogeneity on flows at high velocities through porous media. The media is made of a large number of homogeneous blocks where the flow is assumed to be governed by the Forchheimer equation with a constant inertial coefficient. By assuming the validity of the Forchheimer equation at the large scale, an effective inertial coefficient is deduced from numerical simulations. Different media are investigated: serial-layers, parallel-layers and correlated media. The numerical results show that: (i) for the serial-layers, the effective inertial coefficient is independent of the Reynolds number and decreases when the variance and the mean permeability ratio increases; (ii) for the parallel-layers and the correlated media, the effective inertial coefficient is function of the Reynolds number and increases when the variance and the mean permeability ratio increases. Theoretical relationships are proposed for the inertial coefficient as function of the Reynolds number and the characteristics of the media.     

Pore-Scale Simulation of Fluid Flow Through Deformable Porous Media Using Immersed Boundary Coupled Lattice Boltzmann Method

Transport in Porous Media - The interaction between the fluid flow and the deformable porous media is crucial in the applications of adsorption/absorption. The immersed boundary coupled lattice...     

可变形多孔介质中的一维非定常耦合渗流

在Ｂｉｏｔ理论的基础上,考虑到可变形多孔介质的渗透系数依赖于孔隙变形的特点,建立了耦合渗流问题的基本方程;用初始层校正法求出了一维非定常耦合渗流问题的摄动解;实例计算表明,耦合分析与非耦合分析之间的判别较大,因此耦合效应不能忽略。     

About the Porous Media Flow Through Circular and Elliptical Anisotropic Inhomogeneities

In this article, we describe horizontal groundwater flow due to a uniform flow at infinity around a cylindrical or elliptical inhomogeneity, where the permeability inside the inhomogeneity is anisotropic and different from the isotropic domain outside the inhomogeneity. The orientation of the uniform flow with respect to the orientation of the ellipse is arbitrary as well as the orientation of the anisotropy inside the ellipse. We derive an expression for the ratio of the flow through the ellipse with respect to the flow in the homogeneous case.     

Transient Heat Transport in Gas Flow Through Granular Porous Media

Models are developed to describe the time-evolution of gas and solid temperature profiles in a class of granular porous media; this time-dependence being the result of a heat source supplying energy uniformly to the gas as it flows into the medium. The solid-phase is treated as a fixed, axi-symmetric bed of randomly packed spheres of uniform size and material properties. An analytic expression for the locally averaged surface temperature of the solid phase is formulated in terms of the time history of the corresponding local gas temperature. This approach avoids the need to assume locally uniform temperatures within the solid phase, and provides a coupled model for the thermal energy transport in the system. A numerical scheme for treating the resulting transport equations is outlined, and results presented. A quasi-steady approximation is proposed, and this approximation is assessed by reference to numerical results obtained from the numerical scheme. One application of the work is to fixed-bed catalytic reactors and absorbers, and results are presented indicating how the regeneration times of such systems depend on operating parameters.     

Dynamics of Flow Through Porous Media with Unsteady Phase Permeabilities

The influence of nonequilibrium effects developing due to the formation of an emulsion of each phase (gas bubbles in the water and water dust in the gas) on the flow dynamics is investigated with reference to the displacement of water by a gas. The nonequilibrium effects manifest themselves in a change in the shape of the phase permeability curves (they become convex) and the threshold phase saturations in the course of flow through the porous medium. A kinetic equation in which the relaxation time is proportional to the seepage rate is used to describe such effects. The case in which the liquid displaced by the gas is itself gassed and the volume concentration of the gas bubbles is constant is considered.     

New General Maximum Entropy Model for Flow Through Porous Media

Transport in Porous Media - New experimental and numerical techniques constitute the major recent advancements in the study of flow through porous media. However, a model that duly links the micro-...     

Shock Propagation in a One-Dimensional Flow Through Deformable Porous Media

We consider a one-dimensional incompressible flow through a porous medium undergoing deformations such that the porosity and the hydraulic conductivity can be considered to be functions of the flux intensity. The medium is initially dry and we neglect capillarity, so that a sharp wetting front proceeds into the medium. We consider the open problem of the continuation of the solution possibly in the case of onset of singularities, which can be interpreted as a local collapse of the medium. In particular we analyze the case in which the boundary pressure has a piecewise constant derivative.     

Thermodynamic Automaton Simulations of Fluid Flow and Diffusion in Porous Media

A thermodynamic automaton model of fluid flow in porous media is presented. The model is a nonrelativistic version of a Lorentz invariant lattice gas model constructed by Udey et al. (1998). In the previous model it was shown that the energy momentum tensor and the relativistic Boltzman equation can be rigorously derived from the collision and propagation rules. In the present paper we demonstrate that this nonrelativistic model can be used to accurately simulate well known results involving single phase flow and diffusion in porous media. The simulation results show that (1) one-phase flow simulations in porous media are consistent with Darcy's law; (2) the apparent diffusion coefficient decreases with a decrease in permeability; (3) small scale heterogeneity does not affect diffusion significantly in the cases considered.     

Nonlinear Pressure Diffusion in Flow of Compressible Liquids Through Porous Media

In the flow of liquids through porous media, nonlinear effects arise from the dependence of the fluid density, porosity, and permeability on pore pressure, which are commonly approximated by simple exponential functions. The resulting flow equation contains a squared gradient term and an exponential dependence of the hydraulic diffusivity on pressure. In the limiting case where the porosity and permeability moduli are comparable, the diffusivity is constant, and the squared gradient term can be removed by introducing a new variable y, depending exponentially on pressure. The published transformations that have been used for this purpose are shown to be special cases of the Cole–Hopf transformation, differing in the choice of integration constants. Application of Laplace transformation to the linear diffusion equation satisfied by y is considered, with particular reference to the effects of the transformation on the boundary conditions. The minimum fluid compressibilities at which nonlinear effects become significant are determined for steady flow between parallel planes and cylinders at constant pressure. Calculations show that the liquid densities obtained from the simple compressibility equation of state agree to within 1% with those obtained from the highly accurate Wagner-Pru?  equation of state at pressures to 20 MPa and temperatures approaching 600 K, suggesting possible applications to some geothermal systems.     

Casson Fluid Flow Due to Non-Coaxial Rotation of a Porous Disk and the Fluid at Infinity Through a Porous Medium

The flow of the Casson fluid due to non-coaxial rotation of a disk and the fluid at infinity is investigated. Partial differential equations are made dimensionless and coupled. The exact solution of the resultant nonlinear initial-boundary-value problem is solved by applying the Laplace transform. The shear stresses at the disk surface and the steady state stresses are computed. The effects of dimensionless parameters on the dimensionless primary and secondary velocities are analyzed.     

Modeling Fluid Flow in Fractured Porous Media with the Interfacial Conditions Between Porous Medium and Fracture

One of the most popular models that has been applied to predict the fluid velocity inside the fracture with impermeable walls is the cubic law. It highlights that the mean flux along the fracture is proportional to the cubic of fracture aperture. However, for a fractured porous medium, the normal and tangential interface conditions between the fracture and porous matrix can change the velocity profile inside the fracture. In this paper, a correction factor is introduced for flow equation along the fracture by imposing the continuity of normal and tangential components of velocity at the interface between the fracture and porous matrix. As a result, the mean velocity inside the fracture depends not only on the fracture aperture, but also on a set of non-dimensional numbers, including the matrix porosity, the ratio of intrinsic permeability of fracture to that of matrix, the wall Reynolds number, and the ratio of normal velocity on one wall to the other. Finally, the introduced correction factor is employed within the extended finite element method, which is widely used for numerical simulation of fluid flow within the fractured porous media. Several numerical results are presented for the fluid flow through a specimen containing single fracture, in order to investigate the deviation from the cubic law in different case studies.

     

Free Fluid Flow Through a Visco-Deformable Porous Medium with a Moving Boundary

One-dimensional Darcy-law flow through a porous matrix representing a high-viscosity liquid is investigated. The flow develops in a region which depends on time due to sedimentation. The problem considered simulates the geological process of sedimentation in a basin. In accordance with geological data, the permeability and viscosity coefficients of the matrix are assumed to depend nonlinearly on the porosity. The asymptotic properties of the flow are described for large times. The agreement between the results of asymptotic and numerical solutions is satisfactory at intermediate times and good at large times under the realistic sedimentary basin conditions. The simplicity of the asymptotic solution obtained makes it possible to vary the problem parameters and determine the porosity, pressure, and velocities for particular geological conditions by means of simple calculations.     

Dynamic Network Model of Mobility Control in Emulsion Flow Through Porous Media

Modeling the flow of emulsion in porous media is extremely challenging due to the complex nature of the associated flows and multiscale phenomena. At the pore scale, the dispersed phase size can be of the same order of magnitude of the pore length scale and therefore effective viscosity models do not apply. A physically meaningful macroscopic flow model must incorporate the transport of the dispersed phase through the porous material and the changes on flow resistance due to drop deformation as it flows through pore throats. In this work, we present a dynamic capillary network model that uses experimentally determined pore-level constitutive relationships between flow rate and pressure drop in constricted capillaries to obtain representative transient macroscopic flow behavior emerging from microscopic emulsion flow at the pore level. A parametric analysis is conducted to study the effect of dispersed phase droplet size and capillary number on the flow response to both emulsion and alternating water/emulsion flooding in porous media. The results clearly show that emulsion flooding changes the continuous-phase mobility and consequently flow paths through the porous media, and how the intensity of mobility control can be tuned by the emulsion characteristics.