Upscaling of Nonwetting Phase Residual Transport in Porous Media: A Network Approach

Based on the volume averaging method, a macroscopic model is developed for the upscaling of NAPL transport in a porous medium idealised by a network model. Under the assumption of local mass non-equilibrium, a macroscopic equation involving a dispersion tensor, additional convective terms and a linear form for the interfacial mass flux is obtained. The resolution of the two local closure problems obtained allow the determination of the local properties without adjustable parmeters. These problems are solved in a semi-analytical, semi-numerical manner on the network. The originality of this work is the association of the upscaling by volume averaging method with the network approach. The local properties, including the dispersion tensor and the mass exchange coefficient, can therefore be calculated over a large number of pore-bodies and pore-throats in a computationaly tractable manner, thus leading to more significant results. Results are presented for 3D, spatially periodic models of porous media.     

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.     

Generalized Taylor-Aris moment analysis of the transport of sorbing solutes through porous media with spatially-periodic retardation factor

Taylor-Aris dispersion theory, as generalized by Brenner, is employed to investigate the macroscopic behavior of sorbing solute transport in a three-dimensional, hydraulically homogeneous porous medium under steady, unidirectional flow. The porous medium is considered to possess spatially periodic geochemical characteristics in all three directions, where the spatial periods define a rectangular parallelepiped or a unit-element. The spatially-variable geochemical parameters of the solid matrix are incorporated into the transport equation by a spatially-periodic distribution coefficient and consequently a spatially-periodic retardation factor. Expressions for the effective or large-time coefficients governing the macroscopic solute transport are derived for solute sorbing according to a linear equilibrium isotherm as well as for the case of a first-order kinetic sorption relationship. The results indicate that for the case of a chemical equilibrium sorption isotherm the longitudinal macrodispersion incorporates a second term that accounts for the eflect of averaging the distribution coefficient over the volume of a unit element. Furthermore, for the case of a kinetic sorption relation, the longitudinal macrodispersion expression includes a third term that accounts for the effect of the first-order sorption rate. Therefore, increased solute spreading is expected if the local chemical equilibrium assumption is not valid. The derived expressions of the apparent parameters governing the macroscopic solute transport under local equilibrium conditions agreed reasonably with the results of numerical computations using particle tracking techniques.     

On the spatial-temporal averaging method for modeling transport in porous media

The spatial-temporal averaging procedure is considered with a nonhomogeneous distribution of elementary domains in the spatial-temporal space and the probabilistic interpretation of the ST-averaging is also given. Several averaging theorems and corollaries about the averages of spatial and temporal derivatives are presented and rigorously proved which allow elementary domain to vary in space and time. The macroscopic transport equation in the most general condition and the simplified macroscopic equation under the special form of distributions are developed which may be reduced to the classical macroscopic transport equation as the spatial-temporal average degenerates into the volume average.     

Investigation of the scattering of antiplane shear waves by two griffith cracks using the non-local theory

In this paper, the scattering of harmonic antiplane shear waves by two finite cracks is studied using the non-local theory. The Fourier transform is applied and a mixed boundary value problem is formulated. Then a set of triple integral equations is solved using a new method, namely Schmidt's method. This method is more exact and more reasonable than Eringen's for solving this kind of problem. The result of the stress near the crack tip was obtained. Contrary to the classical elasticity solution, it is found that no stress singularity is present at the crack tip, which can explain the problem of macroscopic and microscopic mechanics.     

Method of homogenization applied to dispersion in porous media

The theory of homogenization which is a rigorous method of averaging by multiple scale expansions, is applied here to the transport of a solute in a porous medium. The main assumption is that the matrix has a periodic pore structure on the local scale. Starting from the pores with the Navier-Stokes equations for the fluid motion and the usual convective-diffusion equation for the solute, we give an alternative derivation of the three-dimensional macroscale dispersion tensor for solute concentration. The original result was first found by Brenner by extending Brownian motion theory. The method of homogenization is an expedient approach based on conventional continuum equations and the technique of multiple-scale expansions, and can be extended to more complex media involving three or more contrasting scales with periodicity in every but the largest scale.     

Maxwell's Equations in Two-Phase Systems II: Two-Equation Model

The method of volume averaging is used to derive the two-equation model for Maxwell's equations in a two-phase system. The analysis provides a set of macroscopic transport equations for the electric and magnetic fields that are completely coupled in terms of the closure problem. When the closure problem is quasi-steady, a formal solution is obtained and estimates are developed for the differences between the averaged fields in the individual phases. These estimates lead to constraints for the condition of local electrodynamic equilibrium.     

Coupled Upscaling Approaches For Conduction, Convection, and Radiation in Porous Media: Theoretical Developments

This study deals with macroscopic modeling of heat transfer in porous media subjected to high temperature. The derivation of the macroscopic model, based on thermal non-equilibrium, includes coupling of radiation with the other heat transfer modes. In order to account for non-Beerian homogenized phases, the radiation model is based on the generalized radiation transfer equation and, under some conditions, on the radiative Fourier law. The originality of the present upscaling procedure lies in the application of the volume averaging method to local energy conservation equations in which radiation transfer is included. This coupled homogenization mainly raises three challenges. First, the physical natures of the coupled heat transfer modes are different. We have to deal with the coexistence of both the material system (where heat conduction and/or convection take place) and the non-material radiation field composed of photons. This radiation field is homogenized using a statistical approach leading to the definition of radiation properties characterized by statistical functions continuously defined in the whole volume of the porous medium. The second difficulty concerns the different scales involved in the upscaling procedure. Scale separation, required by the volume averaging method, must be compatible with the characteristic length scale of the statistical approach. The third challenge lies in radiation emission modeling, which depends on the temperature of the material system. For a semi-transparent phase, this temperature is obtained by averaging the local-scale temperature using a radiation intrinsic average while a radiation interface average is used for an opaque phase. This coupled upscaling procedure is applied to different combinations of opaque, transparent, or semi-transparent phases. The resulting macroscopic models involve several effective transport properties which are obtained by solving closure problems derived from the local-scale physics.     

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.     

Wave propagation in anisotropic media with non-local elasticity

In an effort to understand and quantify the effect of non-local elasticity on the wave propagation response of laminated composite layered media, a frequency-wavenumber domain based finite element method is employed. The developed elements are based on the exact solution in the transformed domain and thus exactly represent the dynamics of a layer. This feature enables to model a layer of any thickness by a single element and drastically reduces the cost of computation. The effect of non-locality on the dispersion relation and in turn on the wave response is compared with local (classical) elasticity solutions. A procedure and sample example is outlined to estimate the magnitude of the non-locality parameter by comparing the dispersion relation with lattice dynamics. The effect of non-locality, in terms of the mode-shift and appearance of dispersion on the modes of Lamb waves is further demonstrated.     

A stochastic complex model with random imaginary noise

In this paper, we study a stochastic complex beam–beam interaction model subjected to random imaginary noise. The general procedure is presented to obtain the Fokker–Planck–Kolmogorov equation (FPK) using stochastic averaging method in the case of a special example. The exact stationary probability of FPK is examined theoretically under certain conditions and then the first and second moments for the amplitude are expressed analytically. Finally, a numerical simulation is performed to verify the theoretical results of moments and excellent agreement can be observed between these two results.     

A Solute Transport in Stratified Media

Two-dimensional and steady solute transport in a stratified porous formation is analysed under assumption that the effect of pore-scale dispersion is negligible. The longitudinal dispersion produced as a result of the vertical variation of hydraulic conductivity is analysed by averaging the variability of a solute flux concentration and conductivity. The evolution of the solute flux concentration is expressed with respect to the correlated variable, that is the travel (arrival) time at a fixed location and the averaging procedure is constructed to satisfy the boundary condition where the inlet concentration is a known function of time. In such a statement, a velocity-averaged solute flux concentration is described by a conventional dispersion model (CDM) with a dispersion coefficient which is a function of the arrival time. It is demonstrated that such CDM satisfies the assumption that hydraulic conductivity of the layers is gamma distributed with the parameter of distribution which is chosen to represent a reasonable value of the field scale solute dispersion. The overall behaviour of the model is illustrated by several examples of two-dimensional mass transport.     

Direct numerical simulation of turbulence over resolved and modeled rough walls with irregularly distributed roughness

To simulate turbulent flow over a rough wall without resolving complicated rough geometries, a macroscopic rough wall model is developed based on spatial (plane) averaging theory. The plane-averaged drag force term, which arises through averaging the Navier–Stokes equations in a plane parallel to a rough wall, can be modeled using a plane porosity and a plane hydraulic diameter. To evaluate the developed model, direct and macroscopic model simulations for turbulence over irregularly distributed semi-spheres at Reynolds number of 300 are carried out using the D3Q27 multiple-relaxation time lattice Boltzmann method. The results show that the developed model can be used to predict rough wall skin friction. The results agree quantitatively with standard turbulence statistics such as mean velocity and Reynolds stress profiles with the fully resolved DNS data. Since velocity dispersion occurs inside the rough wall and is found to contribute to turbulence energy dissipation, which the developed model cannot account for, the developed model fails to reproduce dispersion-related turbulence energy dissipation. However, it is found that the plane-averaged drag force term can successfully recover the deficiency of dispersion-related turbulence energy dissipation.     

Kinetic models for dilute solutions of dumbbells in non-homogeneous flows revisited

We propose a two-fluid theory to model a dilute polymer solution assuming that it consists of two phases, polymer and solvent, with two distinct macroscopic velocities. The solvent phase velocity is governed by the macroscopic Navier–Stokes equations with the addition of a force term describing the interaction between the two phases. The polymer phase is described on the mesoscopic level using a dumbbell model and its macroscopic velocity is obtained through averaging. We start by writing down the full phase-space distribution function for the dumbbells and then obtain the inertialess limits for the Fokker–Planck equation and for the averaged friction force acting between the phases from a rigorous asymptotic analysis. The resulting equations are relevant to the modelling of strongly non-homogeneous flows, while the standard kinetic model is recovered in the locally homogeneous case.     

Method of moments for laminar dispersion in an oscillatory flow through curved channels with absorbing walls

The unsteady dispersion of a solute, when the fluid is driven through a curved channel with absorbing walls by an imposed pulsatile pressure gradient, is studied using the method of moments. The study examines the effect of oscillatory Reynolds number, amplitude/frequency of the pressure pulsation and boundary absorption on the longitudinal dispersion. The methodology involves a set of unsteady integral moment equations obtained by applying the Aris-Barton method of moments on the convective-diffusion equation for a curved channel. Central moments are obtained from the moment equations which are solved by a finite-difference implicit scheme. The effect of curvature and boundary absorption on the effective dispersion coefficient from the initial to the stationary stage of the oscillatory flow is studied. Amplitude of the effective dispersion coefficient is found to increase with curvature and decrease with frequency of the pressure pulsation. For large Peclet number and Schmidt number, the amplitude of the dispersion coefficient can be 1.6 times that in a straight channel at large times. Also, for large times, the amplitude of the dispersion coefficient is twice the amplitude of the dispersion coefficient as α, the frequency parameter changes from 0.5 to 1.0. The axial distributions of mean concentration are determined from the first four central moments by using the Hermite polynomial representation. The effect of curvature is to delay the stationary state and also the approach to normality of the concentration distribution. The study has importance in understanding the spreading of pollutants in tidal basins and natural current fields.     

Approximations for solute transport through porous media with flow transverse to layering

Analytical approaches for the prediction of solute transport in layered porous media are investigated for the case of flow perpendicular to the direction of layering. One approach involves the use of averaging techniques to treat the profile as an equivalent homogeneous medium. The method is demonstrated on hypothetical and laboratory-measured data sets and a criterion for validity of the method is given. The second approach involves the use of time convolution to predict breakthrough curves for layered systems on the assumption that layer interactions have no significant effect on transport. Accuracy criteria are derived by comparing moments of the exact and approximate solutions and it is found that the convolution method has broader applicability than the equivalent single-layer analysis. An extension of the convolution method to include consideration of nonequilibrium transport due to the presence of mobile-immobile regions is presented and demonstrated by analysis of laboratory breakthrough data from a two-layer system exhibiting mobile-immobile regions.     

Investigation of a Griffith crack subject to anti-plane shear by using the non-local theory

Field equations of the non-local elasticity are solved to determine the state of stress in a plate with a Griffith crack subject to the anti-plane shear. Then a set of dual-integral equations is solved using Schmidts method. Contrary to the classical elasticity solution, it is found that no stress singularity is present at the crack tip. The significance of this result is that the fracture criteria are unified at both the macroscopic and the microscopic scales.     

Analysis of One-Dimensional Advection–Diffusion Model with Variable Coefficients Describing Solute Transport in a Porous medium

This work presents the analytical solution and temporal moments of one-dimensional advection–diffusion model with variable coefficients. Two case studies along with the two different sets of boundary conditions are considered at the inlet and outlet of the domain. In the first case, a time-dependent solute dispersion in the homogeneous domain along uniform flow is taken into account, whereas in the second case, due to inhomogeneity of domain, velocity is taken spatially dependent and the dispersion is assumed proportional to the square of the velocity. The Laplace transform is used to obtain the analytical solutions. The analytical temporal moments are derived from the Laplace domain solutions. To verify the correctness of the analytical solutions, a high-resolution second-order finite volume scheme is applied. Different case studies are considered and discussed. Both analytical and numerical results are in good agreement with each other.     

An improved approach for random parametric response of dynamic systems with non-linear inertia

A non-Gaussian closure scheme is developed for determining the stationary response of dynamic systems including non-linear inertia and stochastic coefficients. Numerical solutions are obtained and examined for their validity based on the preservation of moments properties. The method predicts the jump phenomenon, for all response statistics at an excitation level very close to the threshold level of the condition of almost sure stability. In view of the increased degree of non-linearity, resulting from the non-Gaussian closure scheme, the mean square of the response displacement is found to be less than those values predicted by other methods such as the Gaussian closure or the first order stochastic averaging.