Saturation-Dependent Hydraulic Conductivity Anisotropy for Multifluid Systems in Porous Media

2007 ◽  
Vol 6 (4) ◽  
pp. 925-934 ◽  
Author(s):  
Z. F. Zhang ◽  
M. Oostrom ◽  
A.L. Ward
2011 ◽  
Vol 133 (6) ◽  
Author(s):  
O. Ivanchenko ◽  
V. Ivanchenko

Convection-enhanced delivery (CED) is a drug delivery technique used to target specific regions of the central nervous system (CNS) for the treatment of neurodegenerative diseases and cancer while bypassing the blood–brain barrier (BBB). The application of CED is limited by low volumetric flow rate infusions in order to prevent the possibility of backflow. Consequently, a small convective flow produces poor drug distribution inside the treatment region, which can render CED treatment ineffective. Novel catheter designs and CED protocols are needed in order to improve the drug distribution inside the treatment region and prevent backflow. In order to develop novel backflow-free catheter designs, the impact of the micro-fluid injection into deformable porous media was investigated experimentally as well as numerically. Fluid injection into the porous media has a considerable effect on local transport properties such as porosity and hydraulic conductivity because of the local media deformation. These phenomena not only alter the bulk flow velocity distribution of the micro-fluid flow due to the changing porosity, but significantly modify the flow direction, and even the volumetric flow distribution, due to induced local hydraulic conductivity anisotropy. These findings help us to design backflow-free catheters with safe volumetric flow rates up to 10 μl/min. A first catheter design reduces porous media deformation in order to improve catheter performance and control an agent volumetric distribution. A second design prevents the backflow by reducing the porosity and hydraulic conductivity along a catheter’s shaft. A third synergistic catheter design is a combination of two previous designs. Novel channel-inducing and dual-action catheters, as well as a synergistic catheter, were successfully tested without the occurrence of backflow and are recommended for future animal experiments.


2021 ◽  
Author(s):  
Dylan R. Harp ◽  
Vitaly Zlotnik ◽  
Charles J. Abolt ◽  
Brent D. Newman ◽  
Adam L. Atchley ◽  
...  

Abstract. The pathways and timing of drainage from inundated ice-wedge polygon centers in a warming climate have important implications for carbon flushing, advective heat transport, and transitions from carbon dioxide to methane dominated emissions. This research provides intuition on this process by presenting the first in-depth analysis of drainage from a single polygon based on fundamental hydrogeological principles. We use a recently developed analytical solution to provide a baseline for the effects of polygon aspect ratios (radius to thawed depth) and hydraulic conductivity anisotropy (horizontal to vertical hydraulic conductivity) on drainage pathways and temporal depletion of ponded water heights of inundated ice-wedge polygon centers. By varying the polygon aspect ratio, we evaluate the effect of polygon size (width), inter-annual increases in active layer thickness, and seasonal increases in thaw depth on drainage. One of the primary insights from the model is that most inundated ice-wedge polygon drainage occurs along an annular region of the polygon center near the rims. This implies that inundated polygons are most intensely flushed by drainage in an annular region along their horizontal periphery, with implications for transport of nutrients (such as dissolved organic carbon) and advection of heat towards ice wedge tops. The model indicates that polygons with large aspect ratios and high anisotropy will have the most distributed drainage. Polygons with large aspect ratio and low anisotropy will have their drainage most focused near the their periphery and will drain most slowly. Polygons with small aspect ratio and high anisotropy will drain most quickly. Our results, based on idealized scenarios, provide a baseline for further research considering geometric and hydraulic complexities of ice-wedge polygons.


2007 ◽  
Vol 7 (1) ◽  
pp. 59-66 ◽  
Author(s):  
D.E. Rolston

The science of soil-water physics and contaminant transport in porous media began a little more than a century ago. The first equation to quantify the flow of water is attributed to Darcy. The next major development for unsaturated media was made by Buckingham in 1907. Buckingham quantified the energy state of soil water based on the thermodynamic potential energy. Buckingham then introduced the concept of unsaturated hydraulic conductivity, a function of water content. The water flux as the product of the unsaturated hydraulic conductivity and the total potential gradient has become the accepted Buckingham-Darcy law. Two decades later, Richards applied the continuity equation to Buckingham's equation and obtained a general partial differential equation describing water flow in unsaturated soils. For combined water and solute transport, it had been recognized since the latter half of the 19th century that salts and water do not move uniformly. It wasn't until the middle of the 20th century that scientists began to understand the complex processes of diffusion, dispersion, and convection and to develop mathematical formulations for solute transport. Knowledge on water flow and solute transport processes has expanded greatly since the early part of the 20th century to the present.


Author(s):  
William G. Gray ◽  
Michael A. Celia

The mathematical study of flow in porous media is typically based on the 1856 empirical result of Henri Darcy. This result, known as Darcy’s law, states that the velocity of a single-phase flow through a porous medium is proportional to the hydraulic gradient. The publication of Darcy’s work has been referred to as “the birth of groundwater hydrology as a quantitative science” (Freeze and Cherry, 1979). Although Darcy’s original equation was found to be valid for slow, steady, one-dimensional, single-phase flow through a homogeneous and isotropic sand, it has been applied in the succeeding 140 years to complex transient flows that involve multiple phases in heterogeneous media. To attain this generality, a modification has been made to the original formula, such that the constant of proportionality between flow and hydraulic gradient is allowed to be a spatially varying function of the system properties. The extended version of Darcy’s law is expressed in the following form: qα=-Kα . Jα (2.1) where qα is the volumetric flow rate per unit area vector of the α-phase fluid, Kα is the hydraulic conductivity tensor of the α-phase and is a function of the viscosity and saturation of the α-phase and of the solid matrix, and Jα is the vector hydraulic gradient that drives the flow. The quantities Jα and Kα account for pressure and gravitational effects as well as the interactions that occur between adjacent phases. Although this generalization is occasionally criticized for its shortcomings, equation (2.1) is considered today to be a fundamental principle in analysis of porous media flows (e.g., McWhorter and Sunada, 1977). If, indeed, Darcy’s experimental result is the birth of quantitative hydrology, a need still remains to build quantitative analysis of porous media flow on a strong theoretical foundation. The problem of unsaturated flow of water has been attacked using experimental and theoretical tools since the early part of this century. Sposito (1986) attributes the beginnings of the study of soil water flow as a subdiscipline of physics to the fundamental work of Buckingham (1907), which uses a saturation-dependent hydraulic conductivity and a capillary potential for the hydraulic gradient.


Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-13 ◽  
Author(s):  
Chi Yao ◽  
Chen He ◽  
Jianhua Yang ◽  
Qinghui Jiang ◽  
Jinsong Huang ◽  
...  

An original 3D numerical approach for fluid flow in fractured porous media is proposed. The whole research domain is discretized by the Delaunay tetrahedron based on the concept of node saturation. Tetrahedral blocks are impermeable, and fluid only flows through the interconnected interfaces between blocks. Fractures and the porous matrix are replaced by the triangular interface network, which is the so-called equivalent matrix-fracture network (EMFN). In this way, the three-dimensional seepage problem becomes a two-dimensional problem. The finite element method is used to solve the steady-state flow problem. The big finding is that the ratio of the macroconductivity of the whole interface network to the local conductivity of an interface is linearly related to the cubic root of the number of nodes used for mesh generation. A formula is presented to describe this relationship. With this formula, we can make sure that the EMFN produces the same macroscopic hydraulic conductivity as the intact rock. The approach is applied in a series of numerical tests to demonstrate its efficiency. Effects of the hydraulic aperture of fracture and connectivity of the fracture network on the effective hydraulic conductivity of fractured rock masses are systematically investigated.


Sign in / Sign up

Export Citation Format

Share Document