谈谈渗流

渗流是流体在多孔介质中的流动,渗流现象广泛地存在于自然界、工程材料、动物、植物中。多孔介质种类繁多,包括岩石(含各类矿藏)、土壤、生物材料和人工多孔介质材料等。渗流理论已经成为人类开发地下水、地热、石油、天然气、煤炭与煤层气等诸多地下资源的重要理论基础。本文从渗流的基本概念、渗流的分类、渗流的影响因素、渗流的特征以及渗流的研究意义等方面进行了阐述。     

ON FLOW IN POROUS MEDIA

渗流是流体在多孔介质中的流动,渗流现象广泛地存在于自然界、工程材料、动物、植物中。多孔介质种类繁多,包括岩石(含各类矿藏)、土壤、生物材料和人工多孔介质材料等。渗流理论已经成为人类开发地下水、地热、石油、天然气、煤炭与煤层气等诸多地下资源的重要理论基础。本文从渗流的基本概念、渗流的分类、渗流的影响因素、渗流的特征以及渗流的研究意义等方面进行了阐述。     

Nonlinear fluid flow laws for orthotropic porous media

Nonlinear fluid flow laws for orthotropic porous media are written in invariant tensor form. As usual in the theory of fluid flow through porous media [1, 2], the equations contain the flow velocity up to the second power. Expressions that determine the nonlinear resistances to fluid flow are presented and it is shown that, on going over from linear to nonlinear flow laws, the asymmetry effect may manifest itself, that is, the fluid flow characteristics may differ along the same straight line in the positive and negative directions. It is shown that, as compared with the linear fluid flow law for orthotropic media when for three symmetry groups a single flow law is sufficient, in nonlinear laws the anisotropy manifestations are much more variable and each symmetry group must be described by specific equations. A system of laboratory measurements for finding the nonlinear flow characteristics for orthotropic porous media is considered.     

Homotopy Solution for the Channel Flow of a Third Grade Fluid

The solution for the flow of a third grade fluid bounded by two parallel porous plates is given using homotopy analysis method (HAM). A comparison is made with the exact numerical solution for the various values of the physical parameters. It is found that a proper choice of the auxiliary parameter occurring in HAM solution gives very close results.     

Numerical and analytical study to estimate the effect of two length scales upon the permeability of a fibrous porous medium

A fibrous porous medium with two length scales is modeled as a bed of porous cylinders aligned perpendicular to the flow of viscous fluid. The flow behavior is described using Stokes and Darcy flow equations in the regions around (higher length scale) and within the cylinders (lower length scale) respectively. The typical ratio of higher and lower length-scale regions enable us to invoke lubrication approximation and simplify the equations to develop a closed form solution for the overall permeability of this dual-scale porous medium. A parametric analysis is performed to explore the dependence of permeability on factors such as the volumetric ratio of higher and lower length-scale regions, permeability and size of inclusions in the smaller length-scale region. The analytical model is compared with the numerical results and the trend is compared with the experiments.     

A Computational Method for Free-Surface Viscous Flows Through Porous Media

In this paper we have investigated free-surface unsteady viscous flow through porous media using a finite difference technique via primitive formulation. A transformation has been used to convert the irregular flow domain into a rectangular domain. Two illustrations have been presented in case of a standing wave problem and a cavity flow problem with a free surface. The method is found satisfactory for numerical solution for Reynolds number upto 200.     

Finite Element Analysis of Convection Flow Through a Porous Medium in a Horizontal Channel

We analyse the convection flow of a viscous fluid through a horizontal channel enclosing a fully saturated porous medium. The Galerkin finite element analysis is used to discuss the flow and heat transfer through the porous medium using serendipity elements. The velocity, the temperature distributions and the rate of heat transfer are analysed for variations in the governing parameters. The profiles at different vertical levels are asymmetric curves, exhibiting reversal flow everywhere except on the midplane. In a given porous medium, for fixed G or N, the temperature in the fluid region at any position in fluids with a higher Prandtl number, is much higher than in fluids with a lower Prandtl number. Likewise, other parameters being fixed, lesser the permeability of the medium, lower the temperature in the flow field. Nu reduces across the flow at all axial positions, while it enhances along the axial direction of the channel. Nu reduces with decrease in the Darcy parameter D, and thus lesser the permeability of the medium, lesser the rate of heat transfer across the boundary at any axial position of the channel.     

Mixed Convection in a Vertical Porous Channel

A numerical study of mixed convection in a vertical channel filled with a porous medium including the effect of inertial forces is studied by taking into account the effect of viscous and Darcy dissipations. The flow is modeled using the Brinkman–Forchheimer-extended Darcy equations. The two boundaries are considered as isothermal–isothermal, isoflux–isothermal and isothermal–isoflux for the left and right walls of the channel and kept either at equal or at different temperatures. The governing equations are solved numerically by finite difference method with Southwell–Over–Relaxation technique for extended Darcy model and analytically using perturbation series method for Darcian model. The velocity and temperature fields are obtained for various porous parameter, inertia effect, product of Brinkman number and Grashof number and the ratio of Grashof number and Reynolds number for equal and different wall temperatures. Nusselt number at the walls is also determined for three types of thermal boundary conditions. The viscous dissipation enhances the flow reversal in the case of downward flow while it counters the flow in the case of upward flow. The Darcy and inertial drag terms suppress the flow. It is found that analytical and numerical solutions agree very well for the Darcian model. An erratum to this article is available at .     

Magnetohydrodynamics convective flow in a vertical slot through a porous medium

An analysis is presented with magnetohydrodynamics natural convective flow of a viscous Newtonian fluid saturated porous medium in a vertical slot. The flow in the porous media has been modeled using the Brinkman model. The fully-developed two-dimensional flow from capped to open ends is considered for which a continuum of solutions is obtained. The influence of pertinent parameters on the flow is delineated and appropriate conclusions are drawn. The asymptotic behaviour and the volume flux are analyzed and incorporated graphically for the three-parameter family of solution.     

The Flow Problem of Fluids Flow in a Fractal Reservoir with Double Porosity

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆｌｏｗｔｈｅｏｒｙａｎｄｉｔｓａｐｐｌｉｃａｔｉｏｎｏｆｆｌｕｉｄｓｆｌｏｗｉｎａｆｒａｃｔａｌｒｅｓｅｒｖｏｉｒｈａｖｅｃｏｎｔｉｎｕａｌｌｙｇｏｎｅｄｅｅｐｉｎｔｏｓｔｕｄｙｓｉｎｃｅＣｈａｎｇａｎｄＹｏｒｔｓｏｓ[1]ｂｕｉｌｔｔｈｅｆｌｏｗｍｏｄｅｌｏｆｆｌｕｉｄｔｈｒｏｕｇｈａｆｒａｃｔａｌｒｅｓｅｒｖｏｉｒ.ＴＯＮＧＤｅｎｇ＿ｋｅ[2 ]ｐｒｅｓｅｎｔｅｄｔｈｅｅｘａｃｔｓｏｌｕｔｉｏｎａｎｄｉｔｓｐｒｅｓｓｕｒｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｆｏｒｔｈｅｖａ…     

An overview on nonlinear porous flow in low permeability porous media

This paper gives an overview on nonlinear porous flow in low permeability porous media, reveals the microscopic mechanisms of flows, and clarifies properties of porous flow fluids. It shows that, deviating from Darcy's linear law, the porous flow characteristics obey a nonlinear law in a low-permeability porous medium, and the viscosity of the porous flow fluid and the permeability values of water and oil are not constants. Based on these characters, a new porous flow model, which can better describe low permeability reservoir, is established. This model can describe various patterns of porous flow, as Darcy's linear law does. All the parameters involved in the model, having definite physical meanings, can be obtained directly from the experiments.     

Flow of a biomagnetic viscoelastic fluid: application to estimation of blood flow in arteries during electromagnetic hyperthermia,a therapeutic procedure for cancer treatment

The paper deals with the theoretical investigation of a fundamental problem of biomaguetic fluid flow through a porous medium subject to a magnetic field by using the principles of biomagnetic fluid dynamics （BFD）. The study pertains to a situation where magnetization of the fluid varies with temperature. The fluid is considered to be non-Newtonian, whose flow is governed by the equation of a second-grade viscoelastic fluid. The walls of the channel are assumed to be stretchable, where the surface velocity is proportional to the longitudinal distance from the origin of coordinates. The problem is first reduced to solving a system of coupled nonlinear differential equations involving seven parameters. Considering blood as a biomagnetic fluid and using the present analysis, an attempt is made to compute some parameters of the blood flow by developing a suitable numerical method and by devising an appropriate finite difference scheme. The computational results are presented in graphical form, and thereby some theoretical predictions are made with respect to the hemodynamical flow of the blood in a hyperthermal state under the action of a magnetic field. The results clearly indicate that the presence of a magnetic dipole bears the potential so as to affect the characteristics of the blood flow in arteries to a significant extent during the therapeutic procedure of electromagnetic hyperthermia. The study will attract the attention of clinicians, to whom the results would be useful in the treatment of cancer patients by the method of electromagnetic hyperthermia.     

MIXED TRANSFORM FINITE ELEMENT METHOD FOR SOLVING THE NON-LINEAR EQUATION FOR FLOW IN VARIABLY SATURATED POROUS MEDIA

A new computational method is developed for numerical solution of the Richards equation for flow in variably saturated porous media. The new method, referred to as the mixed transform finite element method, employs the mixed formulation of the Richards equation but expressed in terms of a partitioned transform. An iterative finite element algorithm is derived using a Newton–Galerkin weak statement. Specific advantages of the new method are demonstrated with applications to a set of one— dimensional test problems. Comparisons with the modified Picard method show that the new method produces more robust solutions for a broad range of soil– moisture regimes, including flow in desiccated soils, in heterogeneous media and in layered soils with formation of perched water zones. In addition, the mixed transform finite element method is shown to converge faster than the modified Picard method in a number of cases and to accurately represent pressure head and moisture content profiles with very steep fronts. © 1997 by John Wiley & Sons, Ltd.     

Effect of yield stress on the flow of a Casson fluid in a homogeneous porous medium bounded by a circular tube

The effect of yield stress on the flow characteristics of a Casson fluid in a homogeneous porous medium bounded by a circular tube is investigated by employing the Brinkman model to account for the Darcy resistance offered by the porous medium. The non-linear coupled implicit system of differential equations governing the flow is first transformed into suitable integral equations and are solved numerically. Analytical solution is obtained for a Newtonian fluid in the case of constant permeability, and the numerical solution is verified with that of the analytic solution. The effect of yield stress of the fluid and permeability of the porous medium on shear stress and velocity distributions, plug flow radius and flow rate are examined. The minimum pressure gradient required to start the flow is found to be independent of the permeability of the porous medium and is equal to the yield stress of the fluid.     

低渗透多孔介质中的非线性渗流理论

文中论述了低渗透性多孔介质中非线性渗流理论的几个问题,阐明了渗流流体的性质,指出了多孔介质对流体通过的选择性,提出了新的非线性渗流方程,用实验资料对其进行了验证,分析了该方程演变功能,表明它可以描述各种渗流规律.该方程的各项参数都可从实验中直接得到,应用方便,并且参数的物理意义明确.     

Large-eddy simulation of fluid mixing in tee with sintered porous medium

Mixing processes of hot and cold fluids in a tee with and without sintered copper spheres are simulated by FLUENT using the large-eddy simulation (LES) turbulent flow model and the sub-grid scale (SGS) Smagorinsky-Lilly (SL) model with buoyancy. Comparisons of numerical results of the two cases with and without sintered copper spheres show that the porous medium significantly reduces velocity and temperature fluctuations because the porous medium can effectively restrict the fluid flow and enhance heat transfer. The porous medium obviously increases the pressure drop in the main duct. The porous medium reduces the power spectrum density (PSD) of temperature fluctuations in the frequency range from 1 Hz to 10 Hz.     

低渗透多孔介质中的非线性渗流理论简

文中论述了低渗透性多孔介质中非线性渗流理论的几个问题,阐明了渗流流体的性质,指出了多孔介质对流体通过的选择性,提出了新的非线性渗流方程,用实验资料对其进行了验证,分析了该方程演变功能,表明它可以描述各种渗流规律.该方程的各项参数都可从实验中直接得到,应用方便,并且参数的物理意义明确.     

Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs

A numerical method as well as a theoretical study of non-Darcy fluid flow through porous and fractured reservoirs is described. The non-Darcy behavior is handled in a three-dimensional, multiphase flow reservoir simulator, while the model formulation incorporates the Forchheimer equation for describing single-phase or multiphase non-Darcy flow and displacement. The non-Darcy flow through a fractured reservoir is handled using a general dual-continuum approach. The numerical scheme has been verified by comparing its results against those of analytical methods. Numerical solutions are used to obtain some insight into the physics of non-Darcy flow and displacement in reservoirs. In addition, several type curves are provided for well-test analyses of non-Darcy flow to demonstrate a methodology for modeling this type of flow in porous and fractured rocks, including flow in petroleum and geothermal reservoirs.     

A numerical method for forced convection in porous and homogeneous fluid domains coupled at interface by stress jump

A numerical method was developed for flows involving an interface between a homogeneous fluid and a porous medium. It is based on the finite volume method with body‐fitted and multi‐block grids. The Brinkman–Forcheimmer extended model was used to govern the flow in the porous medium region. At its interface, the flow boundary condition imposed is a shear stress jump, which includes the inertial effect, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The forced convection through a porous insert over a backward‐facing step is investigated. The results are presented with flow configurations for different Darcy numbers, 10^?2 to 10^?5, porosity from 0.2 to 0.8, Reynolds number from 10 to 800, and the ratio of insert length to channel height from 0.1 to 0.3. The heat transfer is improved by using porous insert. To enhance the heat transfer with minimal frictional losses, it is preferable to have a medium length of insert with medium Darcy number, and larger Reynolds number. The interfacial stress jump coefficients β and β₁ were varied from ?1 to 1, and within this range the average and local lower‐wall Nusselt numbers are not sensitive to the parameters. Copyright © 2007 John Wiley & Sons, Ltd.     

Forced convection from an isolated heat source in a channel with porous medium

A numerical study is made of flow and heat transfer characteristics of forced convection in a channel that is partially filled with a porous medium. The flow geometry models convective cooling process in a printed circuit board system with a porous insert.The channel walls are assumed to be adiabatic. Comprehensive numerical solutions are acquired to the governing Navier-Stokes equations, using the Brinkman-Forchheimer-extended Darcy model for the regions of porous media. Details of flow and thermal fields are examined over ranges of the principal parameters; i.e., the Reynolds number Re, the Darcy number Da (≡K/H²), the thickness of the porous substrate S, and the ratio of thermal conductivities R_k (≡k_eff/k). Two types of the location of the porous block are considered. The maximum temperature at the heat source and the associated pressure drop are presented for varying Re, Da, S, and R_k. The results illustrate that as S increases or Da decreases, the fluid flow rate increases. Also, as R_k increases for fixed Da, heat transfer rates are augmented. Explicit influences of Re on the flow and heat transport characteristics are also scrutinized. Assessment is made of the utility of using a porous insert by cross comparing the gain in heat transport against the increase in pressure drop.