Macroscopic modelling of transport phenomena in porous media. 2: Applications to mass,momentum and energy transport

In this second paper, the averaging rules presented in Part 1 are employed in order to develop a general macroscopic balance equation and particular equations for mass, mass of a component, momentum and energy, all of a phase in a porous medium domain. These balance equations involve averaged fluxes. Then macroscopic equations are developed for advective, dispersive and diffusive fluxes, all in terms of averaged state variables of the system. These are combined with the macroscopic balance equations to yield field equations that serve as the core of the mathematical models that describe the transport of extensive quantities in a porous medium domain. It is shown that the methodology of averaging leads to a better understanding of the effective stress concept employed in dealing with transport phenomena in deformable porous media.     

Scale-dependent Macroscopic Balance Equations Governing Transport Through Porous Media: Theory and Observations

A theory is developed providing a rational framework for spatial scale- dependent fluid’s flow and heat transfer, and mass of a component migrating with it through porous media. Introducing the assumption of a non-Brownian type motion and referring to asymptotic expansion in powers of a small defined parameter, we develop a novel approach associated with macroscopic balance equations obtained by averaging over a Representative Elementary Volume (REV). We prove that these equations can be decomposed into a primary part that refers to the REV length scale and a secondary part valid at a length scale smaller than that of the corresponding REV length. Further to our previous development, we obtain two general forms of the primary and secondary macroscopic balance equations. One is based on the assumption that the advective flux of the extensive quantity is dominant over that of the dispersive flux, whereas the other disregards this assumption. Moreover we also introduce the primary and secondary macroscopic forms for the fluid heat- transfer equation. Considering a Newtonian fluid, the resulting primary Navier–Stokes equation can vary from a nonlinear wave equation to a drag-dominant equation at the fluid–solid interface (Darcy’s law). The secondary momentum balance equation describes a wave equation governing the concurrent propagation of the intensive momentum and the dispersive momentum flux, deviating from their corresponding average terms. The primary macroscopic fluid heat-transfer equation accounts for advective and dispersive heat fluxes and the secondary macroscopic heat-transfer equation involves the simultaneous advection of heat deviating from its corresponding intensive average quantity. The primary macroscopic solute mass balance equation accounts for advection and hydrodynamic dispersion. The secondary macroscopic component mass balance equation is in the form of pure advection governing migration of the deviation from the average component concentration. At this stage, we focus on establishing the viability of the developed theory. We do this by arguing that field observations of motion at small spatial scales are coherent with the hyperbolic characteristics of the secondary balance equations. Field observations under natural gradient flow conditions show excessive high concentration (average of 50 mg/L) of colloids under land irrigated by sewage effluents. We argue that this displacement of condensed colloidal parcels manifests the theoretical findings for the smaller spatial scale. Further evidence show the accumulation of particles moving behind the front of an emitted shockwave. We consider this as an experimental proof reinforcing the argument that colloidal migration is subject to the action of a shockwave in the fluid and pure advection transport, governed by the respective suggested hyperbolic macroscopic balance equations of fluid momentum and component mass at the smaller spatial scale.     

Advective fluxes in multiphase porous media under nonisothermal conditions

We consider the case in which more than one fluid phase occupies the void space of a porous medium. The advective flux law is formulated for a fluid phase, under nonisothermal conditions and with the presence of solutes in the fluid phases. The derivation of the flux laws is based on an approximated version of the averaged balance equation for linear momentum. Taking into account momentum transfer through the interface between the fluid phases, leads to coupling between the flow in adjacent phases. Fluxes are also shown to depend on the surface tension at the interface between the adjacent fluid phases. Since the latter depends on temperature and solute concentration in the two phases, the advective flux is shown to depend on both temperature and solute concentration gradients in the two phases. A preliminary order of magnitude analysis gives conditions under which the coupling phenomena are not negligible. The approach is applied to the unsaturated zone, as a typical example of a multiphase porous medium.The main conclusion is that the well known Darcy law for single phase flow, may have to be modified for a multi fluid phase system, especially when temperature and solute concentration are not uniform.     

Effects of Material Symmetry on the Coefficients of Transport in Anisotropic Porous Media

The objective of this article is to highlight certain features of a number of coefficients that appear in models of phenomena of transport in anisotropic porous media, especially the coefficient of dispersion the second-rank tensor D _ij , and the dispersivity coefficient, the fourth-rank tensor a _ijkl , that appear in models of solute transport. Although we shall focus on the transport of mass of a dissolved chemical species in a fluid phase that occupies the void space, or part of it, the same discussion is also applicable to transport coefficients that appear in models that describe the advective mass flux of a fluid and the diffusive transport of other extensive quantities, like heat. The case of coupled processes, e.g. the simultaneous transport of heat and mass of a chemical species, are also considered. The entire discussion will be at the macroscopic level, at which a porous medium domain is visualized as a homogenized continuum.     

Approximate Analytical Solutions for Flow of a Third-Grade Fluid Through a Parallel-Plate Channel Filled with a Porous Medium

The flow of a non-Newtonian fluid through a porous media in between two parallel plates at different temperatures is considered. The governing momentum equation of third-grade fluid with modified Darcy’s law and energy equation have been derived. Approximate analytical solutions of momentum and energy equations are obtained by using perturbation techniques. Constant viscosity, Reynold’s model viscosity, and Vogel’s model viscosity cases are treated separately. The criteria for validity of approximate solutions are derived. A numerical residual error analysis is performed for the solutions. Within the validity range, analytical and numerical solutions are in good agreement.     

A numerical method for flows in porous and homogenous fluid domains coupled at the interface by stress jump

A numerical method was developed for flows involving an interface between a homogenous fluid and a porous medium. The numerical method is based on the finite volume method with body‐fitted and multi‐block grids. A generalized model, which includes Brinkman term, Forcheimmer term and non‐linear convective term, was used to govern the flow in the porous medium region. At its interface, a shear stress jump that includes the inertial effect was imposed, together with a continuity of normal stress. Furthermore, the effect of the jump condition on the diffusive flux was considered, additional to that on the convective part which has been usually considered. Numerical results of three flow configurations are presented. The method is suitable for coupled problems with regions of homogeneous fluid and porous medium, which have complex geometries. Copyright © 2006 John Wiley & Sons, Ltd.     

A theoretical analysis of forced convection in a porous-saturated circular tube: Brinkman–Forchheimer model

Fully developed forced convection inside a circular tube filled with saturated porous medium and with uniform heat flux at the wall is investigated on the basis of a Brinkman–Forchheimer model. The matched asymptotic expansion method is applied at small Darcy numbers. For large Darcy numbers, the solution for the Brinkman–Forchheimer momentum equation is found in terms of an asymptotic expansion. Once the velocity distribution is determined, the energy equation is solved using the same asymptotic technique. The results for the two limiting cases of clear fluid and Darcy flow conditions show good agreement with those available in the literature.     

Frictional passage of fluid through inhomogeneous porous system: a variational principle

A Fermat-like principle of minimum time is formulated for nonlinear steady paths of fluid flow in inhomogeneous isotropic porous media where fluid streamlines are curved by a location dependent hydraulic conductivity. The principle describes an optimal nature of nonlinear paths in steady Darcy’s flows of fluids. An expression for the total path resistance leads to a basic analytical formula for an optimal shape of a steady trajectory. In the physical space an optimal curved path ensures the maximum flux or shortest transition time of the fluid through the porous medium. A sort of “law of bending” holds for the frictional fluid flux in Lagrange coordinates. This law shows that—by minimizing the total resistance—a ray spanned between two given points takes the shape assuring that a relatively large part of it resides in the region of lower flow resistance (a ‘rarer’ region of the medium).     

Heat and mass transfer on a MHD third grade fluid with partial slip flow past an infinite vertical insulated porous plate in a porous medium

The influence of third grade, partial slip and other thermophysical parameters on the steady flow, heat and mass transfer of viscoelastic third grade fluid past an infinite vertical insulated plate subject to suction across the boundary layer has been investigated. The space occupying the fluid is porous. The momentum equation is characterized by a highly nonlinear boundary value problem in which the order of the differential equation exceeds the number of available boundary conditions. An efficient numerical scheme of midpoint technique with Richardson’s extrapolation is employed to solve the governing system of coupled nonlinear equations of momentum, energy and concentration. Numerical calculations were carried out for different values of various interesting non-dimensional quantities in the slip flow regime with heat and mass transfer and were shown with the aid of figures. The values of the wall shear stress, the local rate of heat and mass transfers were obtained and tabulated. The analysis shows that as the fluid becomes more shear thickening, the momentum boundary layer decreases but the thermal boundary layer increases; the magnetic field strength is found to decrease with an increasing temperature distribution when the porous plate is insulated. The consequences of increasing the permeability parameter and Schmidt number decrease both the momentum and concentration boundary layer thicknesses respectively whereas an increase in the thermal Grashof number gives rise to the thermal boundary layer thickness.     

Evolution of governing mass and momentum balances following an abrupt pressure impact in a porous medium

A mathematical model is developed of an abrupt pressure impact applied to a compressible fluid flowing through a porous medium domain. Nondimensional forms of the macroscopic fluid mass and momentum balance equations yield two new scalar numbers relating storage change to pressure rise. A sequence of four reduced forms of mass and momentum balance equations are shown to be associated with a sequence of four time periods following the onset of a pressure change. At the very first time period, pressure is proven to be distributed uniformly within the affected domain. During the second time interval, the momentum balance equation conforms to a wave form. The behavior during the third time period is governed by the averaged Navier-Stokes equation. After a long time, the fourth period is dominated by a momentum balance similar to Brinkman's equation which may convert to Darcy's equation when friction at the solid-fluid interface dominates.     

Linear and Non-linear Double Diffusive Convection in a Fluid-Saturated Anisotropic Porous Layer with Cross-Diffusion Effects

The double diffusive convection in a horizontal anisotropic porous layer saturated with a Boussinesq binary fluid, which is heated and salted from below in the presence of Soret and DuFour effects is studied analytically using both linear and non-linear stability analyses. The linear analysis is based on the usual normal mode technique, while the non-linear analysis is based on a minimal representation of double Fourier series. The generalized Darcy model including the time derivative term is employed for the momentum equation. The critical Rayleigh number, wavenumbers for stationary and oscillatory modes, and frequency of oscillations are obtained analytically using linear theory. The effects of anisotropy parameter, solute Rayleigh number, and Soret and DuFour parameters on the stationary, oscillatory convection, and heat and mass transfer are shown graphically. Some known results are recovered as special cases of the present problem.     

The Effect of Rotation on the Onset of Double Diffusive Convection in a Sparsely Packed Anisotropic Porous Layer

The effect of rotation on the onset of double diffusive convection in a sparsely packed anisotropic porous layer, which is heated and salted from below, is investigated analytically using the linear and nonlinear theories. The Brinkman model that includes the Coriolis term is employed for the momentum equation. The critical Rayleigh number, wavenumber for stationary and oscillatory modes and a dispersion relation are obtained analytically using linear theory. The effect of anisotropy parameters, Taylor number, Darcy number, solute Rayleigh number, Lewis number, Darcy–Prandtl number, and normalized porosity on the stationary, oscillatory and finite amplitude convection is shown graphically. It is found that contrary to its usual influence on the onset of convection in the absence of rotation, the mechanical anisotropy parameter show contrasting effect on the onset criterion at moderate and high rotation rates. The nonlinear theory based on the truncated representation of Fourier series method is used to find the heat and mass transfers. The effect of various parameters on heat and mass transfer is shown graphically. Some of the convection systems previously reported in the literature is shown to be special cases of the system presented in this study.     

Linear and nonlinear double diffusive convection in a rotating sparsely packed porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and nonlinear stability analyses. The Brinkman model that includes the Coriolis term is employed as the momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for the energy equation. The onset criterion for stationary, oscillatory, and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal diffusion, solute diffusion, and rotation that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number, and Taylor number on the stability of the system is investigated. The nonlinear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out. Some of the convection systems previously reported in the literature is shown to be special cases of the system presented in this study.     

The effect of a transition layer between a fluid and a porous medium: shear flow in a channel

Flow in a three-layer channel is modeled analytically. The channel consists of a transition layer sandwiched between a porous medium and a fluid clear of solid material. Within the transition layer, the reciprocal of the permeability varies linearly across the channel. The Brinkman model is used for the momentum equations for the porous medium layer and the transition layer. The velocity profile is obtained in closed form in terms of Airy, exponential, and polynomial functions. The overall volume flux and boundary friction factors are calculated and compared with values obtained with a two-layer model employing the Beavers–Joseph condition at the interface between a Darcy porous medium and a clear fluid.     

Jump Condition for Diffusive and Convective Mass Transfer Between a Porous Medium and a Fluid Involving Adsorption and Chemical Reaction

In this paper, mass transfer at the fluid–porous medium boundaries is studied. The problem considers both diffusive and convective transport, along with adsorption and reaction effects in the porous medium. The result is a mass flux jump condition that is expressed in terms of effective transport coefficients. Such coefficients (a total dispersion tensor and effective reaction and adsorption coefficients) may be computed from the solution of the corresponding closure problem here stated and solved as a function of the Péclet number (Pe), the porosity and a local Thiele modulus. For the case of negligible convective transport (i.e., ), the closure problem reduces to the one recently solved by the authors for diffusion and reaction between a fluid and a microporous medium.     

Thermal radiation effect on non-Darcy natural convection with lateral mass transfer

 A boundary layer analysis has been presented to study the influence of thermal radiation and lateral mass flux on non-Darcy natural convection over a vertical flat plate in a fluid saturated porous medium. Forchheimer extension is considered in the flow equations, and the Rosseland approximation is used to describe the radiative heat flux in the energy equation. Similarity solution for the transformed governing equations is obtained and the combined effect of thermal radiation and fluid suction/injection on the heat transfer rate is discussed. Numerical results for the details of the velocity and temperature profiles as well as Nusselt number have been presented. Received on 7 July 1999     

Fully Developed Forced Convection Heat Transfer in a Porous Channel with Asymmetric Heat Flux Boundary Conditions

An analytical study is performed on steady, laminar, and fully developed forced convection heat transfer in a parallel plate channel with asymmetric uniform heat flux boundary conditions. The channel is filled with a saturated porous medium, and the lower and upper walls are subjected to different uniform heat fluxes. The dimensionless form of the Darcy–Brinkman momentum equation is solved to determine the dimensionless velocity profile, while the dimensionless energy equation is solved to obtain temperature profile for a hydrodynamically and thermally fully developed flow in the channel. Nusselt numbers for the lower and upper walls and an overall Nusselt number are defined. Analytical expressions for determination of the Nusselt numbers and critical heat flux ratio, at which singularities are observed for individual Nusselt numbers, are obtained. Based on the values of critical heat flux ratio and Darcy number, a diagram is provided to determine the direction of heat transfer between the lower or upper walls while the fluid is flowing in the channel.     

Theoretical derivation of the conservation equations for single phase flow in porous media: a continuum approach

This paper deals with the theoretical derivation of the conservation equations for single phase flow in a porous medium. The derivation is obtained within the framework of the continuum mechanics and classical thermodynamics. The adopted procedure provides the conservation equations of mass, momentum, mechanical energy, total energy, internal energy, entropy, temperature, enthalpy, Gibbs free energy and Helmholtz free energy. The obtained results highlight the connection between the basic equations of fluid mechanics and of fluid flow in porous media, as well as the restrictions and the limitations of Darcy’s law and Richards’ equation.     

Filtration in a Porous Granular Medium: 1. Simulation of Pore-Scale Particle Deposition and Clogging

This paper presents a numerical model for simulating the pore-scale transport and infiltration of dilute suspensions of particles in a granular porous medium under the action of hydrodynamic and gravitational forces. The formulation solves the Stokes’ flow equations for an incompressible fluid using a fixed grid, multigrid finite difference method and an embedded boundary technique for modeling particle–fluid coupling. The analyses simulate a constant flux of the fluid suspension through a cylindrical model pore. Randomly generated particles are collected within the model pore, initially through contact and attachment at the grain surface (pore wall) and later through mounding close to the pore inlet. Simple correlations have been derived from extensive numerical simulations in order to estimate the volume of filtered particles that accumulate in the pore and the differential pressure needed to maintain a constant flux through the pore. The results show that particle collection efficiency is correlated with the Stokes’ settling velocity and indirectly through the attachment probability with the particle–grain surface roughness. The differential pressure is correlated directly with the maximum mound height and indirectly with particle size and settling velocity that affect mound packing density. Simple modification factors are introduced to account for pore length and dip angle. These parameters are used to characterize pore-scale infiltration processes within larger scale network models of particle transport in granular porous media in a companion paper. Articlenote: Currently at GZA GeoEnvironmental Inc., 1 Edgewater Drive, Norwood, MA 02062, U.S.A.