Exponential stability of stochastic generalized porous media equations with jump

Stochastic generalized porous media equation with jump is considered. The aim is to show the moment exponential stability and the almost certain exponential stability of the stochastic equation.     

On the cross-effects of coupled macroscopic transport equations in porous media

In this paper, the derivation of macroscopic transport equations for this cases of simultaneous heat and water, chemical and water or electrical and water fluxes in porous media is presented. Based on themicro-macro passage using the method of homogenization of periodic structures, it is shown that the resulting macroscopic equations reveal zero-valued cross-coupling effects for the case of heat and water transport as well as chemical and water transport. In the case of electrical and water transport, a nonsymmetrical coupling was found.Notations b mobility - c concentration of a chemical - D rate of deformation tensor - D molecular diffusion coefficient - D _ij ^eff macroscopic (or effective) diffusion tensor - electric field - E ₀ initial electric field - k _ij molecular tensor - j, j ^*, current densities - K _ij macroscopic permeability tensor - l characteristic length of the ERV or the periodic cell - L characteristic macroscopic length - L _ijkl coupled flows coefficients - n _i unit outward vector normal to - p pressure - q _t ,q _t ⁺ , heat fluxes - q _c ,q _c ⁺ , chemical fluxes - s specific entropy or the entropy density - S entropy per unit volume - t time variable - t _ij local tensor - T absolute temperature - v _i velocity - V ₀ initial electric potential - V electric potential - x macroscopic (or slow) space variable - y microscopic (or fast) space variable - _i local vectorial field - _i local vectorial field - electric charge density on the solid surface - , bulk and shear viscosities of the fluid - _ij local tensor - _ij local tensor - _i local vector - _ij molecular conductivity tensor - _ij ^eff effective conductivity tensor - homogenization parameter - fluid density - ₀ ion-conductivity of fluid - _ij dielectric tensor - _i ¹ , _i ² , _i ³ local vectors - ⁴ local scalar - _S solid volume in the periodic cell - _L volume of pores in the periodic cell - boundary between _S and _L - ^s rate of entropy production per unit volume - total volume of the periodic cell - _l volume of pores in the cell On leave from the Politechnika Gdanska; ul. Majakowskiego 11/12, 80-952, Gdask, Poland.     

Averaging methods for numerical simulations of flows through heterogeneous porous media

Many natural rock systems contain small patches of different permeability which affect the flow of fluids through them. As these heterogeneities become smaller and more numerous, they become harder to model numerically. We consider how to reduce the computational effort required in simulations by incorporating their effects in the boundary conditions at the edges of each grid block. This is in contrast with current methods which involve often arbitrary changes in the fluid properties. The method is restricted to the case of widely-spaced patches, which simplifies interaction effects. The system then reduces to an array of dipoles, and two averaging methods are proposed for finite grid blocks. Several infinite systems, including vertical and horizontal bands, are also considered as further approximations. There is a great wealth of existing results from different fields which lead to identical mathematical problems and which can be used in these cases. Finally, we consider how to use these techniques when the precise configuration of the grid block is not known, but only its statistical properties. This can lead to results which are very different from the deterministic case.     

A new stochastic method of reconstructing porous media

We present a new stochastic method of reconstructing porous medium from limited morphological information obtained from two-dimensional micro- images of real porous medium. The method is similar to simulated annealing method in the capability of reconstructing both isotropic and anisotropic structures of multi-phase but differs from the latter in that voxels for exchange are not selected completely randomly as their neighborhood will also be checked and this new method is much simpler to implement and program. We applied it to reconstruct real sandstone utilizing morphological information contained in porosity, two-point probability function and linear-path function. Good agreement of those references verifies our developed method’s powerful capability. The existing isolated regions of both pore phase and matrix phase do quite minor harm to their good connectivity. The lattice Boltzmann method (LBM) is used to compute the permeability of the reconstructed system and the results show its good isotropy and conductivity. However, due to the disadvantage of this method that the connectivity of the reconstructed system’s pore space will decrease when porosity becomes small, we suggest the porosity of the system to be reconstructed be no less than 0.2 to ensure its connectivity and conductivity.     

Analytical solution of partial differential equations for radial transport of a solute in double porous media

The mathematical model for radial transport of a solute is summed up in this paper. The action of non-equilibrium linear adsorption, the double property of porous media and the decay of solute are considered. With the first kind of boundary condition, one finds the analytical solution of these equations by Laplace transform and calculates the dimensionless solution by FORTRAN program with DJS-040. The distribution and change of solute are evaluated and the solution under various limit cases is given. By numerical analysis, one obtains some valuable conclusions.     

Finite element modeling of transport in porous media

Modern computational techniques enable, in principle, the modeling of transport in porous media, involving convection, adsorption and dispersion. Implementation of the techniques for practical problems leads to various difficulties, however. One of these is the difference in horizontal and vertical scales in natural situations; other difficulties encountered are numerical dispersion and the flow near singularities. In order to overcome these difficulties a two-dimensional flow model has been adapted to incorporate three-dimensional velocity components. This procedure takes into account that in regional flow fields the horizontal flow components in aquifers are much larger than the vertical components, and yet it enables to observe transport in vertical direction. Numerical dispersion is suppressed by particle tracking.     

Experimental determination of the flow transport coefficients in the coupled equations of two-phase flow in porous media

The transport coefficients in the coupled equations of two-phase flow are defined if the pressure gradient in one of the two flowing fluids is equal to zero. This definition has been used in experiments with oil and water in a sandpack and the four transport coefficients have been measured over wide water saturation ranges. The values of the cross coefficients were found to be significant as they ranged from 10 to 35% of the value of the effective permeability to water and from 5 to 15% of the effective permeability to oil, respectively.     

Fundamental solution of equations of flow in porous media

A fundamental solution is found to the system of equations which describes the slow flow of an incompressible viscous liquid through a porous material under the influence of a field of body forces.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 169–171, November–December, 1984.     

3D simulation of radionuclide transport in porous media

We present an efficient and easily implementable finite volume method simulating radionuclide transport through highly heterogeneous grounds in three space dimensions. The numerical concentration of the transported chemicals are proved to remain nonnegative and stable. Then, we run a realistic test case in which some radioactive iodine I¹²⁹ particles are released from a leak in an underground nuclear waste disposal site. The question of whether the radionuclide invades the underground and reach the ground surface is investigated. Because of the 3D nature of the problem, a particular emphasis is made on the control of CPU time. Copyright © 2009 John Wiley & Sons, Ltd.     

A scale-dependent equation for solute transport in porous media

The governing equation describing solute transport in porous media is reformulated using standard volume averaging techniques. The alternative formulation is based on a modified definition of the deviation, which allows for variation of macroscopic velocity across the REV. The new equation contains additional scale-dependent terms which are functions of the size of the averaging volume (REV). This result indicates that the scale-dependent nature of the dispersion phenomenon is inherent even at the scale of the REV.     

Deletion of nondominant effects in modeling transport in porous media

A methodology for eliminating nondominant effects in models that describe transport phenomena in porous media is presented. The methodology is based on the introduction of dimensionless numbers and on a proper evaluation of the order of magnitude of terms. These dimensionless numbers are redefined as characteristics of transport and transformation phenomena in porous media. It is shown that different time scales and different length scales may have to be employed for different variables. A method for evaluating the order of magnitude of the error of prediction when terms are deleted, is presented.     

Percolation theory approach to transport phenomena in porous media

The percolation theory approach to static and dynamic properties of the single- and two-phase fluid flow in porous media is described. Using percolation cluster scaling laws, one can obtain functional relations between the saturation fraction of a given phase and the capillary pressure, the relative permeability, and the dispersion coefficient, in drainage and imbibition processes. In addition, the scale dependency of the transport coefficient is shown to be an outcome of the fractal nature of pore space and of the random flow pattern of the fluids or contaminant.     

Determining equations of two-phase flows through anisotropic porous media

A new interpretation of the concept of relative phase permeability is given. Relative phase permeabilities are represented in the form of fourth-rank tensors. It is shown that in the case of anisotropic porous media functions depending not only on the saturation but also on the anisotropy parameters represented in the form of ratios of the principal values of the absolute permeability coefficient tensor correspond to the classical representation of the relative phase permeabilities. For a two-phase flow in anisotropic porous media with orthotropic and transversely-isotropic symmetry a generalized two-term Darcy’s law is analyzed. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 87–94, March–April, 1998. The work was carried out with support from the Russian Foundation for Fundamental Research (project No. 96-01-00623).     

Application of generalised effective-medium theory to transport in porous media

The use of effective-medium treatments to estimate bulk properties pertaining to transport (of, for example, fluids, heat, particles or electricity) through random composite media (such as reservoir rocks), is widespread. This is because they are relatively simple, often reasonably accurate (on occasion, remarkably so) and in many cases yield closed-form expressions for the properties concerned. However, the single-bond effective-medium treatment (EMT) of random resistor networks that has been used to determine transport coefficients for various transport problems in pore networks is limited to some special isotropic networks with nearest-neighbour connections. We demonstrate here that transport through two different fracture system models, with stress-induced anisotropy, can be treated using an EMT originally applied to anisotropic resistor networks. The main purpose of the present contribution, however, is to present a new, more general effective medium formalism applicable to networks of arbitrary topology. This new generalised EMT is used to obtain a new criterion for percolation of an arbitrary conducting network under random dilution. A specific application to unsaturated flow through a pore network with nearest- and next-nearest-neighbour connections is also given.     

On equilibrium and primary variables in transport in porous media

Thermodynamic equilibrium, which involves mechanical, thermal, and chemical equilibria, in a multiphase porous medium, is defined and discussed, both at the microscopic level, and at the macroscopic one. Conditions are given for equilibrium in the presence of forces between the surface of the solid matrix and the fluid phases. The concept ofapproximate thermodynamic equilibrium is introduced and discussed, employing the definition of athermodynamic potential. This discussion serves as the basis for the methodology of determining the number of degrees of freedom in models of phenomena of transport (of mass, energy, and momentum) in porous media. Equilibrium and nonequilibrium cases are considered. The proposed expressions for the number of degrees of freedom in macroscopic transport models, represent the equivalent ofGibbs phase rule in thermodynamics. Based on balance considerations and thermodynamic relationships, it is shown that the number of degrees of freedom, NF, in a problem of transport in a deformable porous medium, involving NP fluid phases and NC components, under nonisothermal conditions, with equilibrium among all phases and components, is $${\text{NF = NC + NP + 4}}{\text{.}}$$ Under nonequilibrium conditions among the phases, the rule takes the form $${\text{NF = NC }} \times {\text{ NP + 2NP + NC + 4}}{\text{.}}$$ In both cases, when fluid phase velocities are determined by Darcy's law, NF is reduced by NP. When the solid matrix is nondeformable, NF is reduced by 3. The number of degrees of freedom is also determined for conditions of approximate chemical and thermal equilibria, and for conditions of equilibrium that prevail only among some of the phases present in the system. Examples of particular cases are presented to illustrate the proposed methodology.     

Qualitative properties of plasticity equations for porous media in plane deformation

The instantaneous stress-strain state of a porous rigid-plastic material obeying the cylindrical yield condition and the associated flow rule is considered in the case of plane deformation. It is shown that the type of the system of equations depends on the stress state. In the hyperbolic case, the equations of characteristics and relations along them are derived. An exact solution to the model boundary-value problem with the maximum friction law taken into account is obtained. An asymptotic analysis near the maximum friction surface is performed.     

Probabilistic approach for transport of contaminants through porous media

The channels formed between individual particles in porous media have variable dimensions and orientations. The porosity, permeability and its anisotropy exhibit random spatial distributions. The probabilistic approach can effectively describe the transport of contaminants through porous media and is analysed in this paper. Numerical results are obtained by considering (I) random dispersion coefficients without and with spatial structure, (II) random time distribution of concentration at the inlet boundary, (III) random velocity distribution in the flow field without and (IV) with variable dispersion coefficient, (V) non-linearity of the governing equation and (VI) anisotropy of the dispersion coefficient. Two methods are used for probabilistic predictions: (1) Gaussian field approach in conjunction with Monte Carlo method and (2) random walk method. The input random parameters are assumed to have normal and log-normal distributions according to available experimental data. The probability distribution functions of the contaminant concentration at different locations within the flow domain are calculated and compared with the input distributions as a function of the mean and fluctuation Peclet numbers. The one-dimensional case is analysed in detail and the illustrative numerical predictions are compared with analytical and experimental results. The extension to a two-dimensional domain is discussed in the last part of this paper.