Contour Dynamics of Incompressible 3-D Fluids in a Porous Medium with Different Densities

We consider the problem of the evolution of the interface given by two incompressible fluids through a porous medium, which is known as the Muskat problem and in two dimensions it is mathematically analogous to the two-phase Hele–Shaw cell. We focus on a fluid interface given by a jump of densities, being the equation of the evolution obtained using Darcy’s law. We prove local well-posedness when the smaller density is above (stable case) and in the unstable case we show ill-posedness.     

Lattice Boltzmann方法模拟多孔介质Klinkenberg效应

格子Boltzmann数值模拟方法是研究复杂的多孔介质结构特别是Klinkenberg效应的有效方法之一，对处理复杂边值问题尤其有效，用格子Boltzmann方法研究了气流穿越多孔介质问题，并将数值计算结果与实验结果进行了比较，结果表明格子Boltzmann方法是数值模拟气流穿越多孔介质问题的有效方法之一。     

Heat and mass transfer on MFD flow through porous media over an accelerating surface in the presence of suction and blowing

The heat and mass transfer of electrically conducting fluid through porous media over an accelerating surface subject both to power law surface temperature and power law heat flux variations with a temperature-dependent heat source in the presence of a transverse uniform magnetic field is studied. A series solution to the energy and species concentration equation in terms of Kummer’s function is studied. The effect of Prandtl number and Schmidt number is studied with the help of graphs.     

Influence of the surrounding pressure on acoustic properties of slightly compressible porous media permeated with air-filled bubbles

The effective wave velocity, attenuation, and nonlinear properties of slightly compressible porous media permeated with air-filled bubbles are studied numerically by employing the nonlinear Hooke’s law for different surrounding pressures. Numerical simulations show that the acoustic properties of porous media are greatly affected by the surrounding pressure if the shear modulus of the elastic medium is very small due to the fact that the acoustic wave propagation in porous media are strongly influenced by the nonlinear oscillation of bubbles; moreover, the oscillation of a bubble depends on the equilibrium bubble radius, which is affected by the surrounding pressures. Published in Russian in Akusticheskiĭ Zhurnal, 2006, Vol. 52, No. 4, pp. 490–496. The text was submitted by the authors in English.     

变形双重介质分形油藏非达西流动分析

考虑与实际生产相符的介质的双重特性和分形特征,并考虑介质的变形,引入双重分形介质渗透率模数,建立应力敏感地层双重分形介质系统的流动方程.采用Douglas-Jones预估-校正法获得无限大地层定产量生产时变形双重分形介质模型的数值解和无限大地层定压生产时分形介质双孔模型的数值解,作出了典型的压力曲线图版.     

随机多孔介质的蒙特卡罗重构与渗流特性模拟

根据多孔介质微观结构的分形尺度标度特征,采用蒙特卡罗方法分别重构随机多孔介质的微观颗粒和孔隙结构,并基于分形毛管束模型研究多尺度多孔介质的气体渗流特性,建立多孔介质微观结构和宏观渗流特性的定量关系。结果表明：分形蒙特卡罗重构的多孔介质微细结构接近真实介质结构,气体渗流特性的计算结果与格子玻尔兹曼模拟数据较为吻合; 多孔介质气体渗透率随着克努森数的增加而增大,孔隙分形维数对于气体渗流的微尺度效应具有显著影响,而迂曲度分形维数对于表观渗透率和固有渗透率的比值影响可以忽略。提出的分形蒙特卡罗方法具有收敛速度快且计算误差与维数无关的优点,有利于深入理解多尺度多孔介质的渗流机理。     

Influence of gas law on ultrasonic behaviour of porous media under pressure

This paper deals with the influence of gas law on ultrasonic behaviour of porous media when the saturating fluid is high pressured. Previous works have demonstrated that ultrasonic transmission through a porous sample with variations of the static pressure (up to 18 bars) of the saturating fluid allows the characterization of high damping materials. In these studies, the perfect gas law was used to link static pressure and density, which is disputable for high pressures. This paper compares the effects of real and perfect gas laws on modeled transmission coefficient for porous foams at these pressures. Direct simulations and a mechanical parameters estimation from minimization show that results are very similar in both cases. The real gas law is thus not necessary to describe the acoustic behaviour of porous media at low ultrasonic frequencies (100 kHz) up to 20 bars.     

Underground Storage of Natural Gas in Hydrate State: Primary Injection Stage

The paper is devoted to simulation of the initial stage of natural gas hydrate underground storage: gas injection into aquifer just below permafrost rocks. It is based on the mathematical model of multiphase non-isothermal real gas and water flow in porous media. The model takes into account the transformation of gas and water into hydrate at certain temperature, which depends on gas flow pressure. The dynamics of hydrate and water saturation as well as the pressure and temperature fields in a reservoir with given porosity, permeability, and initial values of pressure, temperature and water saturation were studied. An implicit finite-difference scheme is used to approximate the original boundary-value problem. The finite-difference equations are solved using simple iteration and sweeping algorithms. Several examples of calculations corresponding to real cases are given. Calculations have revealed that the final result strongly depends on the combination of porosity and permeability of a reservoir.     

Extended Hydrodynamics from Enskog’s Equation for a Two-Dimensional System General Formalism

Balance equations are derived from Enskog’s kinetic equation for a two-dimensional system of hard disks using Grad’s moment expansion method. This set of equations constitute an extended hydrodynamics for moderately dense bi-dimensional fluids. The set of independent hydrodynamic fields in the present formulations are: density, velocity, temperature and also—following Grad’s original idea—the symmetric and traceless pressure tensor p _ij and the heat flux vector q ^k . An approximation scheme similar in spirit to one made by Grad in his original work is made. Once the hydrodynamics is derived it is used to discuss the nature of a simple one-dimensional heat conduction problem. It is shown that, not too far from equilibrium, the nonequilibrium pressure in this case only depends on the density, temperature and heat flux vector. PACS: 51.10.+y, 05.20.Jj, 44.10.+i, 05.70.Ln     

Validity of the permeability Carman-Kozeny equation: A volume averaging approach

A volume averaging approach is used to estimate the porous media permeability. Contrary to traditional methods that rely on solving the Navier-Stokes equations for laminar flow, this approach has the advantage that it does not require the specification of some physical conditions and parameters (pressure drop and viscosity). Numerical results on synthetic models of porous media showed that (i) the local porous medium configuration has an important effect on the permeability value, and (ii) the Carman-Kozeny equation cannot describe the permeability behavior as a function of porosity and characteristic lengths. In turn, our results indicate that simple empirical equations, commonly used in practice, are unable to describe the permeability functionalities over a broad range of porous media configurations.     

Equation of State for Parallel Rigid Spherocylinders

The pair distribution function of monodisperse rigid spherocylinders is calculated by Shinomoto’s method, which was originally proposed for hard spheres. The equation of state is derived by two different routes: Shinomoto’s original route, in which a hard wall is introduced to estimate the pressure exerted on it, and the virial route. The pressure from Shinomoto’s original route is valid only when the length-to-width ratio is less than or equal to 0.25 (i.e., when the spherocylinders are nearly spherical). The virial equation of state is shown to agree very well with the results of numerical simulations of spherocylinders with length-to-width ratio greater than or equal to 2.     

A multiscale poromicromechanical approach to wave propagation and attenuation in bone

Ultrasonics is an important diagnostic tool for bone diseases, as it allows for non-invasive assessment of bone tissue quality through mass density–elasticity relationships. The latter are, however, quite complex for fluid-filled porous media, which motivates us to develop a rigorous multiscale poromicrodynamics approach valid across the great variety of different bone tissues. Multiscale momentum and mass balance, as well as kinematics of a hierarchical double porous medium, together with Darcy’s law for fluid flow and micro–poro-elasticity for the solid phase of bone, give access to the so-called dispersion relation, linking the complex wave numbers to corresponding wave frequencies. Experimentally validated results show that 2.25 MHz acoustical signals transmit healthy cortical bone (exhibiting a low vascular porosity) only in the form of fast waves, agreeing very well with experimental data, while both fast and slow waves transmit highly osteoporotic as well as trabecular bone (exhibiting a large vascular porosity). While velocities and wavelengths of both fast and slow waves, as well as attenuation lengths of slow waves, are always monotonously increasing with the permeability of the bone sample, the attenuation length of fast waves shows a minimum when considered as function of the permeability.     

Computation of collective modes and acoustic investigations at different temperatures of vitrous silica

This paper is mainly concerned with elastic and acoustic properties of vitrous silica besides the computation of phonon frequencies. Thus the phonon frequencies of vitrous silica have been calculated assuming the electronic bulk modulus,K _e, as equal to zero. New equations have been derived to relate the pressure derivatives of second order elastic constants to the acoustic Gruneisen’s parameters using both Bhatia-Singh’s parameters and Schofield’s equations. The calculated longitudinal and transverse Gruneisen’s parameters and the predicted absorption band spectra from Nagendranath’s equation and Bhatia Singh’s parameters are in good agreement with experiment. The calculated mean acoustic mode Gruneisen’s parameter evaluated from the pressure derivative of Nagendranath’s equation is also in good agreement with experiment. An erratum to this article is available at .     

On the theory of acoustothermometry of waterlike media: Effects of the quasi-static field,strong absorption,and radiation pattern

An integral equation relating the thermal noise’s pressure measured by an acoustothermometer to the one-dimensional temperature profile of a waterlike medium is derived by applying the fluctuation-dissipative theorem to hydrodynamic equations. On the basis of this equation, effects of the quasi-stationary field of thermal radiation, a strong absorption of acoustic waves in the medium, and a finite beam width of the receiving antenna are investigated. Conditions under which the solution to the equation coincides with the result of the transfer theory, which ignores the aforementioned effects, are determined. The role of the effects under study in acoustothermometry of biological media is investigated. A method is proposed for controlling the depth from which the radiation is received. This method allows the retrieval of the temperature profile of the medium from the data of acoustic sounding.     

A Maximum Principle for the Muskat Problem for Fluids with Different Densities

We study the fluid interface problem given by two incompressible fluids with different densities evolving by Darcy’s law. This scenario is known as the Muskat problem for fluids with the same viscosities, being in two dimensions mathematically analogous to the two-phase Hele-Shaw cell. We prove in the stable case (the denser fluid is below) a maximum principle for the L ^∞ norm of the free boundary.     

Effective medium method of slightly compressible elastic media permeated with air-filled bubbles

An effective medium method is developed for the slightly compressible elastic media permeated with air-filled bubbles, according to the nonlinear oscillation of the bubble, which happens when compressional wave travels through the porous media. The effective Lame coefficients of the porous medium and the nonlinear elastic wave equation are deduced, based on the fact that the micro-unit of the effective medium should have the same stress and strain as the micro-unit of the porous media. The linearized properties obtained by this method are in good agreement with the results of Gaunaurd’s classic theory [Gaunaurd G.C. and überall H., J. Acoust. Soc. Am., 1978, 63: 1699–1711]. Furthermore, the nonlinear coefficient, which is an important property of the porous media, can also be acquired by this method. __________ Translated from Acta Acustica, 2006 (in Chinese) (in press)     

A fractional-order Darcy's law

By using spatial averaging methods, in this work we derive a Darcy's-type law from a fractional Newton's law of viscosity, which is intended to describe shear stress phenomena in non-homogeneous porous media. As a prerequisite towards this end, we derive an extension of the spatial averaging theorem for fractional-order gradients. The usage of this tool for averaging continuity and momentum equations yields a Darcy's law with three contributions: (i) similar to the classical Darcy's law, a term depending on macroscopic pressure gradients and gravitational forces; (ii) a fractional convective term induced by spatial porosity gradients; and (iii) a fractional Brinkman-type correction. In the three cases, the corresponding permeability tensors should be computed from a fractional boundary-value problem within a representative cell. Consistency of the resulting Darcy's-type law is demonstrated by showing that it is reduced to the classical one in the case of integer-order velocity gradients and homogeneous porous media.     

地层渗透参数反演计算

利用渗流源网格压力(及井底压力)对网格渗透率导数关系(称为渗透率压力场敏感系数)实现由三维渗流方程压强数值解计算井底压力对网格渗透参数导数的计算.基于渗透率压力场敏感系数的计算结果构造地层渗透参数反演方程,由多井系统井压监测数据实现反演地层渗透参数的计算.算例显示:反演方法对均匀(或分块均匀)油藏模型的计算取得了较好的结果,反演迭代10次后计算结果与模型参数的最大相对误差收敛于1%之内.该方法对油藏动态模拟研究及实现油藏参数动态跟踪具有实际意义.     

气孔分布特性对含水多孔介质气体扩散的影响

本文基于多孔介质的气孔分布特性,计算了多孔介质在含水状态下的扩散性能,并且比较了采用两种方式计算相对渗透率时的相对扩散性能。其结果表明,基于气孔分布的计算结果低于与气孔分布无关的计算结果。另外,疏水性含水多孔介质的扩散性能低于亲水性含水多孔介质的扩散性能,基于气孔分布计算含水多孔介质的气体扩散性能时,Wyllie公式并不适用。     

Characteristic regimes of premixed gas combustion in high-porosity micro-fibrous porous media

Dynamical behaviour of the premixed flame propagating in the inert high-porosity micro-fibrous porous media has been studied numerically. Effects of mixture filtration velocity, equivalence ratio and burner transverse size on the flame structure have been investigated and the regions of existence of different combustion regimes have been determined. It was found that the influence of the hydrodynamic instability on the flame dynamics is significant in the case of the moderate and high filtration velocities and this effect is negligible at the low velocities. At the moderate filtration velocities the effect of hydrodynamic instability manifests in the flame front deformation and in particular in the flame inclination. It was found that the flame can be stabilized within the whole interval of the filtration gas velocity, whereas in the ordinary porous media the standing wave is settled only at fixed value of gas filtration velocity. This finding is in line with recent experimental results on combustion in micro-fibrous porous media (Yang et al., Combust. Sci. Tech. 181 (2009), 1–16). Possible physical interpretation of the flame anchoring effect may be given on the base of present numerical analysis. At the high filtration velocities the hydrodynamic instability manifests itself in periodical appearance of the moving wrinkles on the flame front surface which forms non stationary high temperature trailing spots behind the leading part of the flame front. Such dynamics may be associated with splitting wave structures which were revealed in previous experiments (Yang et al., Combust. Sci. Tech. 181 (2009), 1–16).