Modeling of Solute Transport in a 3D Rough-Walled Fracture–Matrix System

Fluid flow and solute transport in a 3D rough-walled fracture–matrix system were simulated by directly solving the Navier–Stokes equations for fracture flow and solving the transport equation for the whole domain of fracture and matrix with considering matrix diffusion. The rough-walled fracture–matrix model was built from laser-scanned surface tomography of a real rock sample, by considering realistic features of surfaces roughness and asperity contacts. The numerical modeling results were compared with both analytical solutions based on simplified fracture surface geometry and numerical results by particle tracking based on the Reynolds equation. The aim is to investigate impacts of surface roughness on solute transport in natural fracture–matrix systems and to quantify the uncertainties in application of simplified models. The results show that fracture surface roughness significantly increases heterogeneity of velocity field in the rough-walled fractures, which consequently cause complex transport behavior, especially the dispersive distributions of solute concentration in the fracture and complex concentration profiles in the matrix. Such complex transport behaviors caused by surface roughness are important sources of uncertainty that needs to be considered for modeling of solute transport processes in fractured rocks. The presented direct numerical simulations of fluid flow and solute transport serve as efficient numerical experiments that provide reliable results for the analysis of effective transmissivity as well as effective dispersion coefficient in rough-walled fracture–matrix systems. Such analysis is helpful in model verifications, uncertainty quantifications and design of laboratorial experiments.     

Two approaches to modeling unstable flow and mixing of variable density fluids in porous media

Density variation of aqueous phase fluids flowing in a porous medium, resulting from spatial and temporal variation of solute concentration, often gives rise to unstable flow, and therefore has a significant effect on solute transport. Studies on simulating unstable flow and mixing of variable density fluids in seemingly homogeneous porous media are rare. In this study, a three-dimensional (3-D) and a one-dimensional (1-D) model were developed to simulate unstable flow and mixing in a vertical, nominally 1-D system. In the 3-D model, the fluid flow and solute transport equations were solved numerically with a very fine spatial discretization. The 1-D numerical model was derived from a theoretical model to simulate the flow and mixing of fluids with variable density and viscosity at the field scale. To evaluate the models, simulated results were compared with experimental data from displacement experiments in a vertical sand column. The results show that the 1-D model provides fairly good prediction of breakthrough curves and that the 3-D model is able to qualitatively simulate breakthrough curves for highly unstable flow and mixing.Contribution from the Alabama Agric. Exp. Sta. as AAES Journal No. 3-955037.     

Propagation of a pre-existing turbulent fluid fracture

The propagation of rough and smooth wall pre-existing turbulent fluid fractures is investigated. The laminar fluid fracture is included as a special case for comparison. Lubrication theory is assumed to apply in the fracture and turbulence is introduced through the wall shear stress. The Perkins–Kern–Nordgren approximation is made in which the fluid pressure is proportional to the half-width of the fracture. The fracture half-width satisfies a non-linear diffusion equation. By using a linear combination of the Lie point symmetries of the non-linear diffusion equation a group invariant solution for the fracture length, volume and half-width is derived. The evolution of the length, half-width and mean flow velocity is analysed for a range of working conditions at the fracture entry. It is found that the mean flow velocity increases approximately linearly along the fracture.     

A thermodynamic approach to the rheology of highly interactive filler-polymer mixtures: Part I – Theory

The rheological behavior of highly interactive filler-polymer mixtures is simulated utilizing a double network created by the entangled polymer matrix and the adsorbed polymer. Both networks are represented by a nonlinear viscoelastic constitutive equation. The dependence of rheological properties on filler concentration is taken into account through the bridging density resulting from polymer-filler interactions and a hydrodynamic reinforcement. The relative contribution of both networks is computed through the energy balance consistent with the thermodynamics of the polymer-filler chemical interactions and fluid mechanics. This self-consistent approach allows one to calculate the strain dependence of dynamic properties under oscillatory flow and shear rate dependence of stresses under steady simple shear flow and upon start up and cessation of shear flow. Received: 11 May 2000 Accepted: 8 March 2001     

Reactive Solute Transport in a Nonstationary,Fractured Porous Medium: A Dual-Porosity Approach

A numerical method of moment is developed for solute flux through a nonstationary, fractured porous medium. Solute flux is described as a space-time process where time refers to the solute flux breakthrough and space refers to the transverse displacement distribution at a control plane. A first-order mass diffusion model is applied to describe interregional mass diffusion between fracture (advection) and matrix (nonadvection) regions. The chemical is under linear equilibrium sorption in both fracture and matrix regions. Hydraulic conductivity in the fracture region is assumed to be a spatial random variable. In this study, the general framework of Zhang et al.(2000) is adopted for solute flux in a nonstationary flow field. A time retention function related to physical and chemical sorption in the dual-porosity medium is developed and coupled with solute advection along random trajectories. The mean and variance of total solute flux are expressed in terms of the probability density function of the parcel travel time and transverse displacement. The influences of various factors on solute transport are investigated. These factors include the interregional mass diffusion rate between fracture and matrix regions, chemical sorption coefficients in both regions, water contents in both regions, and location of the solute source. In comparison with solute transport in a one-region medium, breakthrough curves of the mean and variance of the total solute flux in a two-region medium have lower peaks and longer tails. As compared with the classical stochastic studies on solute transport in fractured media, the numerical method of moment provides an approach for applying the stochastic method to study solute transport in more complicated fractured media.     

Fluid Flow Through Single Fractures With Directional Shear Dislocations

This paper numerically investigates the fluid flow behavior through single fractures with directional shear dislocations. Synthetic fractures are generated with directional shear dislocations, and the lattice Boltzmann method is used to simulate the fracture flow. With an ignorance of tortuosity effect, a notable overestimation of hydraulic conductivity is observed when the simplified local cubic law is used. During the closure process, the decreasing rate of conductivity is found to be highly related to the roughness of fractures. The conductivity of smoother fractures decreases faster than that of rougher fractures. By conducting simulations on fractures with a constant shear displacement, the effective conductivity is found to vary with the shear directions. The results show that the conductivity of rougher fractures is less sensitive to the shear directions than that of smoother fractures. As fracture surfaces come into contact, a sharp decrease in effective conductivity is observed and the decreasing trend flattens as the contact ratio continues to increase. A new model is proposed based on the bottleneck model to predict the conductivity of sheared fractures. By integrating the tortuosity and channeling effects into the original model, the proposed new model shows a better performance in predicting the conductivity, especially for fractures with rougher surfaces.     

裂缝气液重力置换的流动实验及仿真

通过可视化模拟实验对真实裂缝下气液重力置换现象进行了模拟,得到了其气液界面现象与常规平板裂缝的差异;通过对压力、置换量等数据的处理证实了气液重力置换的机理;根据实验的液体漏失流量数据提出了钻井过程中漏失、置换、溢流区间;最后根据k-ε双方程模型进行了气液重力置换流动的仿真计算,得到三维裂缝流道中都是钻井液下部侵入,呈现一个近似的直角三角形形状。本文对现场钻井工程方案设计及安全钻井具有指导意义。     

On the Errors Associated with Two-Dimensional Stochastic Solute Transport Models

Many subsurface solute transport studies employ numerical modeling techniques to estimate solute arrival times. Simplifying assumptions must be made to define the modeling domain within a mathematical framework. One common assumption is that the vertical flow is negligible such that the flow field can be simulated with a two-dimensional model. Reducing the vertical dimension reduces the number of flow paths that a solute can take. In a heterogenous medium, artificially removing the 3rd dimension may lead to erroneous results. We investigate the error in the simulated solute breakthrough associated with a two-dimensional model. We also use a stochastic solution of solute arrival time to derive a transform of a two-dimensional ln (k) field so that solute transport more closely resembles three-dimensional transport behavior. The moment equations for two- and three-dimensional domains were solved simultaneously to calculate this transform. The results indicate that the removal of the vertical variability (3D 2D) introduces a 5–10% error in the predicted solute breakthrough. The error tends to increase with increased hydraulic conductivity variance. Numerical experiments confirm that the transform developed herein decreases the relative error of particle breakthrough curves.     

Mesoscopic Modeling of Capillarity-Induced Two-Phase Transport in a Microfluidic Porous Structure

Mesoscopic modeling at the pore scale offers great promise in exploring the underlying structure transport performance of flow through porous media. The present work studies the fluid flow subjected to capillarity-induced resonance in porous media characterized by different porous structure and wettability. The effects of porosity and wettability on the displacement behavior of the fluid flow through porous media are discussed. The results are presented in the form of temporal evolution of percentage saturation and displacement of the fluid front through porous media. The present study reveals that the vibration in the form of acoustic excitation could be significant in the mobilization of fluid through the porous media. The dependence of displacement of the fluid on physicochemical parameters like wettability of the surface, frequency along with the porosity is analyzed. It was observed that the mean displacement of the fluid is more in the case of invading fluid with wetting phase where the driving force strength is not so dominant.     

A Numerical Method of Moments for Solute Transport in Nonstationary Flow Fields

A Lagrangian perturbation approach has been applied to develop the method of moments for predicting mean and variance of solute flux through a three-dimensional nonstationary flow field. The flow nonstationarity may stem from medium nonstationarity, finite domain boundaries, and/or fluid pumping and injecting. The solute flux is described as a space–time process where time refers to the solute flux breakthrough and space refers to the transverse displacement distribution at the control plane. The analytically derived moment equations for solute transport in a nonstationary flow field are too complicated to solve analytically, a numerical finite difference method is implemented to obtain the solutions. This approach combines the stochastic model with the flexibility of the numerical method to boundary and initial conditions. This method is also compared with the numerical Monte Carlo method. The calculation results indicate the two methods match very well when the variance of log-conductivity is small, but the method of moment is more efficient in computation.     

Vortices and wall shear stresses in asymmetric and symmetric channels with sinusoidal wavy walls for pulsatile flow at low Reynolds numbers

Numerical simulation and flow visualization were performed to study the dynamical behavior of vortices generated in channels with two different geometries, i.e., a periodically converging–diverging channel and serpentine channel, both having sinusoidal wavy walls. This system for pulsatile flow is used to enhance heat and mass transfer in very viscous liquids. The numerical results predict well the dynamical behavior of vortices and agree with the flow visualizations. For both channels, the vortex expands in each furrow of the channel walls during the deceleration phase and shrinks during the acceleration phase, which leads to fluid exchange between the vortex and the mainstream. The time-averaged vortex strength and wall shear stresses increase, as the frequency of fluid oscillation increases under a fixed oscillatory fraction of the flow rate. However, above a certain value of the frequency, they reversely decrease due to viscous effects. This frequency for the serpentine channel is smaller than that for the converging–diverging channel. The channel geometries are found to have an important effect of the flow characteristics.     

明渠层流近床面颗粒拖曳力研究

泥沙颗粒受到的拖曳力是泥沙运动的主要驱动力，而当前应用于计算流体力学-离散颗粒法(CFD-DPM)耦合模型进行水沙运动模拟的泥沙颗粒拖曳力公式均没有考虑明渠流底床边壁作用的影响。求解不可压缩Navier-Stokes方程，对明渠层流不同雷诺数条件下床面附近不同高度处颗粒所受拖曳力进行了模拟，根据模拟结果变化规律，提出了综合考虑床面和水流惯性对标准拖曳力影响的修正拖曳力计算公式。与常用的单颗粒标准拖曳力公式和考虑遮蔽效应的多颗粒拖曳力公式相比，采用本文修正公式得到的水沙作用力更接近高精度数值解，应用于CFD-DPM输沙模拟获得的输沙结果与输沙率公式结果一致，应用分析表明输沙模拟应当采用粗糙底床边界。     

Pore-scale network model for drainage-dominated three-phase flow in porous media

Drainage displacements in three-phase flow under strongly wetting conditions are completely described by a simple generalisation of well understood two-phase drainage mechanisms. As in two-phase flow, the sequence of throat invasions in three-phase flow is determined by fluid connectivity and threshold capillary pressure for the invading interface. Flow through wetting and intermediate spreading films is important in determining fluid recoveries and the progress of the displacement in three-phase flow. Viscous pressure drops associated with flow through films give rise to multiple filling and emptying of pores. A three-phase, two-dimensional network model based on the pore-scale fluid distributions and displacement mechanisms reported by Øren et al. and which accounts for flow through both wetting and intermediate fluid films is shown to correctly predict all the important characteristics of three-phase flow observed in glass micromodel experiments.     

The Dual-Reciprocity Boundary Element Analysis for Hydraulically Fractured Shale Gas Reservoirs Considering Diffusion and Sorption Kinetics

This work presents a new application of boundary element method (BEM) to model fluid transport in unconventional shale gas reservoirs with discrete hydraulic fractures considering diffusion, sorption kinetics and sorbed-phase surface diffusion. The fluid transport model consists of two governing partial differential equations (PDEs) written in terms of effective diffusivities for free and sorbed gases, respectively. Boundary integral formulations are analytically derived using the fundamental solution of the Laplace equation for the governing PDEs and Green’s second identity. The domain integrals arising due to the time-dependent function and nonlinear terms are transformed into boundary integrals employing the dual-reciprocity method. This transformation retains the domain-integral-free, boundary-integral-only character of standard BEM approaches. In the proposed solution, the free- and sorbed-gas flow in the shale matrix is solved simultaneously after coupling the fracture flow equation of free gas. Well production performance under the effect of relaxation phenomenon due to delayed responses of sorbed gas under nonequilibrium sorption condition is rigorously captured by imposing the zero-flux condition at fracture–matrix interface for the sorbed-gas transport equation. The validity of proposed solution is verified using several case studies through comparison against a commercial finite-element numerical simulator.

     

A model for solute transport following an abrupt pressure impact in saturated porous media

A mathematical model is developed of an abrupt pressure impact applied to a compressible fluid with solute, flowing through saturated porous media. Nondimensional forms of the macroscopic balance equations of the solute mass and of the fluid mass and momentum lead to dominant forms of these equations. Following the onset of the pressure change, we focus on a sequence of the first two time intervals at which we obtain reduced forms of the balance equations. At the very first time period, pressure is proven to be distributed uniformly within the affected domain, while solute remains unaffected. During the second time period, the momentum balance equation for the fluid conforms to a wave form, while the solute mass balance equation conforms to an equation of advective transport. Fluid's nonlinear wave equation together with its mass balance equation, are separately solved for pressure and velocity. These are then used for the solution of solute's advective transport equation. The 1-D case, conforms to a pressure wave equation, for the solution of fluid's pressure and velocity. A 1-D analytical solution of the transport problem, associates these pressure and velocity with an exponential power which governs solute's motion along its path line.     

On study of horizontal thin film flow of Sisko fluid due to surface tension gradient

The present paper is concerned with the steady thin film flow of the Sisko fluid on a horizontal moving plate, where the surface tension gradient is a driving mechanism. The analytic solution for the resulting nonlinear ordinary differential equation is obtained by the Adomian decomposition method (ADM). The physical quantities are derived including the pressure profile, the velocity profile, the maximum residue time, the stationary points, the volume flow rate, the average film velocity, the uniform film thickness, the shear stress, the surface tension profile, and the vorticity vector. It is found that the velocity of the Sisko fluid film decreases when the fluid behavior index and the Sisko fluid parameter increase, whereas it increases with an increase in the inverse capillary number. An increase in the inverse capillary number results in an increase in the surface tension which in turn results in an increase in the surface tension gradient on the Sisko fluid film. The locations of the stationary points are shifted towards the moving plate with the increase in the inverse capillary number, and vice versa locations for the stationary points are found with the increasing Sisko fluid parameter. Furthermore, shear thinning and shear thickening characteristics of the Sisko fluid are discussed. A comparison is made between the Sisko fluid film and the Newtonian fluid film.     

New Correlative Models to Improve Prediction of Fracture Permeability and Inertial Resistance Coefficient

Presence of fracture roughness and occurrence of nonlinear flow complicate fluid flow through rock fractures. This paper presents a qualitative and quantitative study on the effects of fracture wall surface roughness on flow behavior using direct flow simulation on artificial fractures. Previous studies have highlighted the importance of roughness on linear and nonlinear flow through rock fractures. Therefore, considering fracture roughness to propose models for the linear and nonlinear flow parameters seems to be necessary. In the current report, lattice Boltzmann method is used to numerically simulate fluid flow through different fracture realizations. Flow simulations are conducted over a wide range of pressure gradients through each fracture. It is observed that creeping flow at lower pressure gradients can be described using Darcy’s law, while transition to inertial flow occurs at higher pressure gradients. By detecting the onset of inertial flow and regression analysis on the simulation results with Forchheimer equation, inertial resistance coefficients are determined for each fracture. Fracture permeability values are also determined from Darcy flow as well. According to simulation results through different fractures, two parametric expressions are proposed for permeability and inertial resistance coefficient. The proposed models are validated using 3D numerical simulations and experimental results. The results obtained from these two proposed models are further compared with those obtained from the conventional models. The calculated average absolute relative errors and correlation coefficients indicate that the proposed models, despite their simplicity, present acceptable outcomes; the models are also more accurate compared to the available methods in the literature.     

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.     

Modeling interface shear behavior of granular materials using micro-polar continuum approach

Recently, the authors have focused on the shear behavior of interface between granular soil body and very rough surface of moving bounding structure. For this purpose, they have used finite element method and a micro-polar elasto-plastic continuum model. They have shown that the boundary conditions assumed along the interface have strong influences on the soil behavior. While in the previous studies, only very rough bounding interfaces have been taken into account, the present investigation focuses on the rough, medium rough and relatively smooth interfaces. In this regard, plane monotonic shearing of an infinite extended narrow granular soil layer is simulated under constant vertical pressure and free dilatancy. The soil layer is located between two parallel rigid boundaries of different surface roughness values. Particular attention is paid to the effect of surface roughness of top and bottom boundaries on the shear behavior of granular soil layer. It is shown that the interaction between roughness of bounding structure surface and the rotation resistance of bounding grains can be modeled in a reasonable manner through considered Cosserat boundary conditions. The influence of surface roughness is investigated on the soil shear strength mobilized along the interface as well as on the location and evolution of shear localization formed within the layer. The obtained numerical results have been qualitatively compared with experimental observations as well as DEM simulations, and acceptable agreement is shown.