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

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

Generalized Darcy's Law and the Structure of the Phase and Relative Phase Permeabilities for Two-Phase Flows through Anisotropic Porous Media

For two-phase immiscible fluid flows a generalized Darcy's law is written in invariant tensor form for crystallographic point symmetry groups and anisotropic textures. The representation of the phase permeability coefficient tensors and the structure of the expressions for the relative phase permeabilities are analyzed for all symmetry groups. The relation between the phase and absolute permeability coefficient tensors is specified by a fourth-rank tensor with the external symmetry coinciding with external symmetry of the phase permeability tensors. It is shown that the external symmetry of the phase permeability coefficient tensors can differ from the external symmetry of the absolute permeability tensor. For triclinic and monoclinic symmetry groups it is shown that the phase permeability coefficient tensors may not be coaxial with each other and with the absolute permeability tensor; moreover, the directions of the principal axes of the phase permeability coefficient tensors can depend on the saturation.     

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.     

Coriolis Effects on Filtration Law in Rotating Porous Media

We investigate the filtration law of incompressible viscous Newtonian fluids in rigid non-inertial porous media, for example, rotating porous media. The filtration law is obtained by upscaling the flow at the pore scale. We use the method of multiple scale expansions which gives rigorously the macroscopic behaviour without any prerequisite on the form of the macroscopic equations. For finite Ekman numbers the filtration law is shown to resemble a Darcy's law, but with a non-symmetric permeability tensor which depends on the angular velocity of the porous matrix. We obtain the filtration analog of the Hall effect. For large Ekman numbers the filtration law is a small correction to the classical Darcy's law. The corrector is antisymmetric. In this case we recover a structure of law which is similar to phenomenological laws introduced in the literature, but with a dissimilar effective coefficient.     

Infiltration of Liquid Metals in Porous Compacts: Modeling of Permeabilities During Reactive Melt Infiltration

Reactive infiltration is a fast and cost-effective technique for manufacturing ceramic-matrix composites (CMCs). CMCs are used in elevated temperature applications like rocket engine casings, jet nozzles, gas turbine blades and nuclear cladding. There is an urgent need for minimizing experimental costs as well as optimizing process parameters during manufacture, so that we have minimized manufacturing costs and reduced infiltration times. Towards this end, the objective of this research was to develop an integrated micro-macro model of reactive flow of molten silicon in a porous preform consisting of carbon-coated silicon carbide fibers and then optimize process parameters computationally. The overall objective of the research was to arrive at a modified equation of Darcy's law for flow through a porous medium with the help of numerical/computational modeling. This paper deals with the flow of silicon through porous carbon at the macro level. The macro flow of silicon was integrated with an available micro model by determining the transient porosity from the micro model and using it in Darcy's law written for the macro flow of silicon. From the results of this study, we recommend suitable process parameters such as initial temperature of the solid reactant and the specific kind of reactants to be used for achieving complete infiltration. These conclusions are drawn after observation of the rate of decrease of permeability with more reaction.     

Flows of constant stretch history for polymeric materials with power-law distributions of relaxation times

It is herein shown that for separable integral constitutive equations with power-law distributions of relaxation times, the streamlines in creeping flow are independent of flow rate.For planar flows of constant stretch history, the stress tensor is the sum of three terms, one proportional to the rate-of-deformation tensor, one to the square of this tensor, and the other to the Jaumann derivative of the rate-of-deformation tensor. The three tensors are the same as occur in the Criminale-Ericksen-Filbey Equation, but the coefficients of these tensors depend not only on the second invariant of the strain rate, but also on another invariant which is a measure of flow strength. With the power-law distribution of relaxation times, each coefficient is equal to the second invariant of the strain rate tensor raised to a power, times a function that depends only on strength of the flow. Axisymmetric flows of constant stretch history are more complicated than the planar flows, because three instead of two nonzero normal components appear in the velocity gradient tensor. For homogeneous axisymmetric flows of constant stretch history, the stress tensor is given by the sum of the same three terms. The coefficients of these terms again depend on the flow strength parameter, but in general the dependences are not the same as in planar flow.     

Motion through spherical droplet with non-homogenous porous layer in spherical container

The problem of the creeping flow through a spherical droplet with a nonhomogenous porous layer in a spherical container has been studied analytically. Darcy's model for the flow inside the porous annular region and the Stokes equation for the flow inside the spherical cavity and container are used to analyze the flow. The drag force is exerted on the porous spherical particles enclosing a cavity, and the hydrodynamic permeability of the spherical droplet with a non-homogeneous porous layer is calculated. Emphasis is placed on the spatially varying permeability of a porous medium, which is not covered in all the previous works related to spherical containers. The variation of hydrodynamic permeability and the wall effect with respect to various flow parameters are presented and discussed graphically. The streamlines are presented to discuss the kinematics of the flow. Some previous results for hydrodynamic permeability and drag forces have been verified as special limiting cases.     

Averaged model with capillary nonequilibrium effects for two-phase flow through a highly heterogeneous porous medium

Two-phase flow through a medium with two porosities in which the absolute permeabilities and the capillary pressure functions of the components differ by an order of magnitude is investigated. A classification and diagram of the elementary flows are proposed at the single cell level. An averaged model is developed for a single class of systems in which source-type capillary-dispersion flow predominates in the blocks. This model contains a nonlinear kinetic relation between the average values of the capillary pressure functions. An expansion of the effective phase permeability tensor allowing it to be calculated efficiently is proposed. The capillary relaxation time is explicitly determined. Examples of calculations of the averaged phase permeability tensor and the capillary relaxation time are given. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 93–103, May–June, 1998. The work was carried out with the support from the Russian Foundation for Basic Research (project No. 95-01-01179a).     

On the Validity of Darcy's Law for Stable High-Concentration Displacements in Granular Porous Media

Recently, it has been suggested that Darcy's Law might not be applicable for modelling miscible, density-dependent flow in porous media. To investigate this, three sets of careful laboratory column experiments were performed on coarse and medium sands, consisting of upward displacement of water by sodium chloride solutions with concentrations ranging from 5 to 200g/l. Data on salt concentrations and water pressures were collected in horizontal transects along the flow direction. Salt concentration data were also collected in the influent and exit lines. The experimental data were analysed using a simplified approach based on Darcy's Law alone, applied with the assumption of a sharp interface. Darcy's Law was used to estimate porous medium permeability by fitting predictions to experimental data. Consistent estimates of permeability were obtained for each set of experiments. The results indicate that Darcy's Law adequately describes high concentration displacements through saturated coarse- and medium-grained porous media.     

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.     

Thermomechanics of Swelling Biopolymeric Systems

A two-scale theory for the swelling biopolymeric media is developed. At the microscale, the solid polymeric matrix interacts with the solvent through surface contact. The relaxation processes within the polymeric matrix are incorporated by modeling the solid phase as viscoelastic and the solvent phase as viscous at the mesoscale. We obtain novel equations for the total stress tensor, chemical potential of the solid phase, heat flux and the generalized Darcy's law all at the mesoscale. The constitutive relations are more general than those previously developed for the swelling colloids. The generalized Darcy's law could be used for modeling non-Fickian fluid transport over a wide range of liquid contents. The form of the generalized Fick's law is similar to that obtained in earlier works involving colloids. Using two-variable expansions, thermal gradients are coupled with the strain rate tensor for the solid phase and the deformation rate tensor for the liquid phase. This makes the experimental determination of the material coefficients easier and less ambiguous.     

Analytical Solutions to Gas Flow Problems in Low Permeability Porous Media

In most of conventional porous media the flow of gas is basically controlled by the permeability and the contribution of gas flow due to gas diffusion is ignored. The diffusion effect may have significant impact on gas flow behavior, especially in low permeability porous media. In this study, a dual mechanism based on Darcy flow as well as diffusion is presented for the gas flow in homogeneous porous media. Then, a novel form of pseudo pressure function was defined. This study presents a set of novel analytical solutions developed for analyzing steady-state and transient gas flow through porous media including effective diffusion. The analytical solutions are obtained using the real gas pseudo pressure function that incorporates the effective diffusion. Furthermore, the conventional assumption was used for linearizing the gas flow equation. As application examples, the new analytical solutions have been used to design new laboratory and field testing method to determine the porous media parameters. The proposed laboratory analysis method is also used to analyze data from steady-state flow tests of three core plugs. Then, permeability (k) and effective diffusion coefficient (D _e) was determined; however, the new method allows one to analyze data from both transient and steady-state tests in various flow geometries.     

The influence of Hall current on the rotating oscillating flows of an Oldroyd-B fluid in a porous medium

In this article, analysis is presented to study the effect of Hall current on the rotating flow of a non-Newtonian fluid in a porous medium taking into consideration the modified Darcy's law. The Oldroyd-B fluid model is used to characterize the non-Newtonian fluid behavior. The governing equations for unsteady rotating flow have been modeled in a porous medium. The analysis includes the flows induced by general periodic oscillations and elliptic harmonic oscillations of a plate. The effect of the various emerging parameters is discussed on the velocity distribution. The analytical results are confirmed mathematically by giving comparison with previous studies in the literature. It is observed that the velocity distribution increases with an increase of Hall parameter. The behavior of permeability is similar to that of the Hall parameter.     

Flow of Maxwell fluids in porous media

In this paper we analyze the flow of a Maxwell fluid in a rigid porous medium using the method of volume averaging. We first present the local volume averaged momentum equation which contains Darcy-scale elastic effects and undetermined integrals of the spatial deviations of the pressure and velocity. A closure problem is developed in order to determine the spatial deviations and thus obtain a closed form of the momentum equation that contains a time-dependent permeability tensor. To gain some insight into the effects of elasticity on the dynamics of flow in porous media, the entire problem is transformed to the frequency domain through a temporal Fourier transform. This leads to a dynamic generalization of Darcy's law. Analytical results are provided for the case in which the porous medium is modeled as a bundle of capillary tubes, and a scheme is presented to solve the transformed closure problem for a general microstructure.     

A multi-scale theory of swelling porous media: I. Application to one-dimensional consolidation

A theory is developed which describes flow in multi-scale, saturated swelling media. To upscale information, both the hybrid theory of mixtures and the homogenization technique are employed. In particular, a model is formulated in which vicinal water (water adsorbed to the solid phase) is treated as a separate phase from bulk (non-vicinal) water. A new form of Darcy's law governing the flow of both vicinal and bulk water is derived which involves an interaction potential to account for the swelling nature of the system. The theory is applied to the classical one-dimensional consolidation problem of Terzaghi and to verify Low's empirical, exponential, swelling result for clay at the macroscale.     

On the transient non-Fickian dispersion theory

The Fickian dispersion equation is the basic relationship used to describe the nonconvective mass flux of a solute in a porous medium. This equation prescribes a linear relationship between the dispersive mass flux and the concentration gradient. An important characteristic of the Fickian relationship is that it is independent of the history of dispersion (e.g. the time rate of change of the dispersion flux). Also, the dispersivities are supposed to be medium constants and invariant with temporal and spatial scales of observation. It is believed that in general these restrictions do not hold. A number of authors have proposed various alternative relationships. For example, differential equations have been employed that prescribe a relationship between the dispersion flux and its time and space derivatives. Also, stochastic theories result in integro-differential equations in which dispersion tensor grow asymptotically with time or distance. In this work, three different approaches, which lead to three different non-Fickian equations with a transient character, are discussed and their primary features and differences are highlighted. It is shown that an effective dispersion tensor defined in the framework of the transient non-Fickian theory, grows asymptotically with time and distance; a result which also follows from stochastic theories. Next, principles of continuum mechanics are employed to provide a solid theoretical basis for the non-Fickian transient dispersion theory. The equation of motion of a solute in a porous medium is used to provide a rigorous derivation of various dispersion relationships valid under different conditions. Under various simplifying assumptions, the generalized theory is found to agree with the conventional Fickian theory as well as several other non-Fickian relationships found in the literature. Moreover, it is shown that for nonconservative solutes, the traditional dispersion tensor is affected by the rate of mass exchange of the solute.Also with National Institute of Public Health and Environmental Protection (RIVM), PO Box 1; 3720BA Bilthoven, The Netherlands     

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

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

Large-scale averaging analysis of multiphase flow in fractured reservoirs

The method of large-scale averaging is applied to derive and analyze a dual-porosity model of multiphase flow in naturally fractured reservoirs. The dual-porosity model contains the usual equations based on Darcy's law, and the coupling terms representing the fluid transfer between the matrix and the fractures. Both quasisteady and transient closure schemes are considered to obtain and analyze the fracture and matrix permeability tensors and the fluid transfer terms. The techniques developed here are not restricted to regular geometric fractures. Computational work aimed at showing the implications of the theory behind the derivation of the present dual-porosity model is also described. In particular, comparisons among the dual-porosity model, the single porosity model, and other dual-porosity models are presented through numerical experiments.