首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Liquid distributions in unsaturated porous media under different gravitational accelerations and corresponding macroscopic gaseous diffusion coefficients were investigated to enhance understanding of plant growth conditions in microgravity. We used a single-component, multiphase lattice Boltzmann code to simulate liquid configurations in two-dimensional porous media at varying water contents for different gravity conditions and measured gas diffusion through the media using a multicomponent lattice Boltzmann code. The relative diffusion coefficients (D rel) for simulations with and without gravity as functions of air-filled porosity were in good agreement with measured data and established models. We found significant differences in liquid configuration in porous media, leading to reductions in D rel of up to 25% under zero gravity. The study highlights potential applications of the lattice Boltzmann method for rapid and cost-effective evaluation of alternative plant growth media designs under variable gravity.  相似文献   

2.
The finite-element method based on a Galerkin technique was used to formulate the problem of simulating the two-dimensional (cross-sectional) transient movement of water and solute in saturated or partially saturated nonuniform porous media. The numerical model utilizes linear triangular elements. Nonreactive, as well as reactive solutes whose behaviour can be described by a distribution coefficient or first-order reaction term were considered. The flow portion of the model was tested by comparison of the model results with experimental and finite-difference results for transient flow in an unsaturated sand column and the solute transport portion of the model was tested by comparison with analytical solution results. The model was applied to a hypothetical case involving movement of water and solutes in tile-drained soils. The simulation results showed the development of distinct solute leaching patterns in the soil as drainage proceeded. Although applied to a tile drainage problem in this study, the model should be equally useful in the study of a wide range of two-dimensional water and solute migration problems.  相似文献   

3.
Diffusive mass exchange into immobile water regions within heterogeneous porous aquifers influences the fate of solutes. The percentage of immobile water is often unidentified in natural aquifers though. Hence, the mathematical prediction of solute transport in such heterogeneous aquifers remains challenging. The objective of this study was to find a simple analytical model approach that allows quantifying properties of mobile and immobile water regions and the portion of immobile water in a porous system. Therefore, the Single Fissure Dispersion Model (SFDM), which takes into account diffusive mass exchange between mobile and immobile water zones, was applied to model transport in well‐defined saturated dual‐porosity column experiments. Direct and indirect model validation was performed by running experiments at different flow velocities and using conservative tracer with different molecular diffusion coefficients. In another column setup, immobile water regions were randomly distributed to test the model applicability and to determine the portion of immobile water. In all setups, the tracer concentration curves showed differences in normalized maximum peak concentration, tailing and mass recovery according to their diffusion coefficients. These findings were more pronounced at lower flow rates (larger flow times) indicating the dependency of diffusive mass exchange into immobile water regions on tracers' molecular diffusion coefficients. The SFDM simulated all data with high model efficiency. Successful model validation supported the physical meaning of fitted model parameters. This study showed that the SFDM, developed for fissured aquifers, is applicable in porous media and can be used to determine porosity and volume of regions with immobile water. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

4.
A macroscopic form of Ohm's law is obtained for isotropic porous media saturated with an electrically conductive fluid by using volumetric averaging concepts. Closure of the macroscopic charge transport equation is aided by approximative modelling of the average geometric structures of three different types of isotropic porous media, namely foamlike materials, granular media and crossflow over prismatic bundles. Modelling of the microscopic charge transport necessitated the introduction of a representative interstitial flux of charge carriers and required quantification of the geometric tortuosity applicable to transport phenomena in general. Deterministic expressions for the formation factor are obtained and compare favourably with experimental results.  相似文献   

5.
In the dispersion theory, a linear relationship has been verified between the coefficient of hydrodynamic dispersion and water velocity, both in saturated and in unsaturated porous media. But for unsaturated soils the variability of flow directions and microscopic velocities can be larger than in saturated soils because of the lower degree of water saturation. This leads to an increased dispersion. Therefore, relationships between water content and relative water velocity fluctuations and water content together with the coefficient of dispersivity in unsaturated porous media respectively have been investigated systematically by displacement experiments in glass beads and coarse-textured sandy soil columns. The breakthrough curves (BTCs) of chloride showed that an increase of solute mixing with a decrease of water content was caused by an increase of flow velocity fluctuations for different pathways. In order to explain the observed tailing effect in unsaturated flow, two mathematical models were used to fit theoretically derived nonlinear functions of water content dependent dispersivities for both porous media. The close agreement between the observed and computed results suggests that the theoretical model of hydrodynamic dispersion can be extended to transport in unsaturated porous media, providing that BTCs of the effluent water are used to estimate representative dispersivity parameters of soils.  相似文献   

6.
This pore-scale modeling study in saturated porous media shows that compound-specific effects are important not only at steady-state and for the lateral displacement of solutes with different diffusivities but also for transient transport and solute breakthrough. We performed flow and transport simulations in two-dimensional pore-scale domains with different arrangement of the solid grains leading to distinct characteristics of flow variability and connectivity, representing mildly and highly heterogeneous porous media, respectively. The results obtained for a range of average velocities representative of groundwater flow (0.1–10 m/day), show significant effects of aqueous diffusion on solute breakthrough curves. However, the magnitude of such effects can be masked by the flux-averaging approach used to measure solute breakthrough and can hinder the correct interpretation of the true dilution of different solutes. We propose, as a metric of mixing, a transient flux-related dilution index that allows quantifying the evolution of solute dilution at a given position along the main flow direction. For the different solute transport scenarios we obtained dilution breakthrough curves that complement and add important information to traditional solute breakthrough curves. Such dilution breakthrough curves allow capturing the compound-specific mixing of the different solutes and provide useful insights on the interplay between advective and diffusive processes, mass transfer limitations, and incomplete mixing in the heterogeneous pore-scale domains. The quantification of dilution for conservative solutes is in good agreement with the outcomes of mixing-controlled reactive transport simulations, in which the mass and concentration breakthrough curves of the product of an instantaneous transformation of two initially segregated reactants were used as measures of reactive mixing.  相似文献   

7.
Wave-induced flow is observed as the dominated factor for P wave propagation at seismic frequencies. This mechanism has a mesoscopic scale nature. The inhomogeneous unsaturated patches are regarded larger than the pore size, but smaller than the wavelength. Surface wave, e.g., Rayleigh wave, which propagates along the free surface, generated by the interfering of body waves is also affected by the mesoscopic loss mechanisms. Recent studies have reported that the effect of the wave-induced flow in wave propagation shows a relaxation behavior. Viscoelastic equivalent relaxation function associated with the wave mode can describe the kinetic nature of the attenuation. In this paper, the equivalent viscoelastic relaxation functions are extended to take into account the free surface for the Rayleigh surface wave propagation in patchy saturated poroelastic media. Numerical results for the frequency-dependent velocity and attenuation and the time-dependent dynamical responses for the equivalent Rayleigh surface wave propagation along an interface between vacuum and patchy saturated porous media are reported in the low-frequency range (0.1–1,000 Hz). The results show that the dispersion and attenuation and kinetic characteristics of the mesoscopic loss effect for the surface wave can be effectively represented in the equivalent viscoelastic media. The simulation of surface wave propagation within mesoscopic patches requires solving Biot’s differential equations in very small grid spaces, involving the conversion of the fast P wave energy diffusion into the Biot slow wave. This procedure requires a very large amount of computer consumption. An efficient equivalent approach for this patchy saturated poroelastic media shows a more convenient way to solve the single phase viscoelastic differential equations.  相似文献   

8.
Water retention characteristic or water retention curve (WRC) is an important constitutive feature of porous media, and also meanwhile is an indispensable requirement in hydraulic transport modelling. Previous experiments have indicated that the specific surface area of porous media has effects on the WRC. It has also been observed that a linear relationship generally exists between the air–water interface area and the water saturation within unsaturated porous media. However it seems that no detailed study on their internal linkage with the WRC has been reported yet. This paper, at first gives a review of the development of WRC modelling, then it tries to explain the water retention characteristic according to the physical and chemical behaviours of the phases involved in unsaturated porous media. Using the traditional capillary theory, the volume averaging theorem and the advances in physical chemistry of interfacial surfaces, this paper then derives out a formula which represents the water retention characteristic of porous media. This formula demonstrates the internal linkage of the WRC to the specific surface area of porosities. It also shows agreements with the experimental observations mathematically. Based on this formula, a fitting model is proposed for the static WRC of porous media. Finally, this model is tested to fit the WRC measurements of a wide range of porous materials. Comparison with other main models is presented.  相似文献   

9.
While the tortuosity coefficient is commonly estimated using an expression based on total porosity, this relationship is demonstrated to not be applicable (and thus is often misapplied) over a broad range of soil textures. The fundamental basis for a correlation between the apparent diffusion tortuosity coefficient and hydraulic conductivity is demonstrated, although such a relationship is not typically considered. An empirical regression for estimating the tortuosity coefficient based on hydraulic conductivity for saturated, unconsolidated soil is derived based on results from 14 previously reported diffusion experiments performed with a broad range of soil textures. Analyses of these experimental results confirm that total porosity is a poor predictor for the tortuosity coefficient over a large range of soil textures. The apparent diffusion tortuosity coefficient is more reliably estimated based on hydraulic conductivity.  相似文献   

10.
Three-dimensional analytical solutions for solute transport in saturated, homogeneous porous media are developed. The models account for three-dimensional dispersion in a uniform flow field, first-order decay of aqueous phase and sorbed solutes with different decay rates, and nonequilibrium solute sorption onto the solid matrix of the porous formation. The governing solute transport equations are solved analytically by employing Laplace, Fourier and finite Fourier cosine transform techniques. Porous media with either semi-infinite or finite thickness are considered. Furthermore, continuous as well as periodic source loadings from either a point or an elliptic source geometry are examined. The effect of aquifer boundary conditions as well as the source geometry on solute transport in subsurface porous formations is investigated.  相似文献   

11.
The knowledge of hydraulic and electric properties of porous media and the relations between them is essential for the quantitative evaluation of electric well logs and the solution of other reservoir engineering problems. The state of the art in this field is not yet satisfying. Theories still show considerable discrepancies with practice due to oversimplified model approaches. Empirical relations are either too coarse, not fully determined, or valid under specific geological and geographical conditions only. This article deals with the development of a general theory of the electric and hydraulic resistance behavior of porous media on the basis of a very general statistical network model. A general solution of the relations between formation factor, permeability, and porosity is presented by means of a rigorous mathematical treatment of two limiting cases of such a network. The solution shows that the product of the formation factor and the permeability can be expressed in the expectation values and the variation coefficients of pore channel cross section and shape factor and by a network factor, that depends on the mesh texture of the network. This network factor is in the range zero to one. It is further shown that the path length increase enters both the electric and the hydraulic tortuosity by its square.  相似文献   

12.
We present a new field measurement and numerical interpretation method (combined termed “test”) to parameterize the diffusion of trichloroethene (TCE) and its biodegradation products (DPs) from the matrix of sedimentary rock. The method uses a dual-packer system to interrogate a low-permeability section of the rock matrix adjacent to a previously contaminated borehole and uses the borehole monitoring history to establish the pretest condition. TCE and its DPs are removed from the groundwater between the packers at the onset of the testing. The parameters estimated by fitting a radial diffusion model to the concentration history and borehole concentration data, also termed back diffusion, are the tortuosity factor and sorption coefficients of TCE and DPs in the rock matrix and the TCE and DP biodegradation rate coefficients in the borehole. We demonstrate the equipment design and the interpretive method using a borehole accessing the gray mudstone at a TCE contaminated site in the Newark Basin. In this test, both nonreactive (bromide) and reactive (trichlorofluoroethene) tracers are used to constrain the estimated parameters; however, the bromide tracer was not needed to estimate the parameters in this test. The parameters estimated from the field test are consistent with values measured independently in laboratory experiments using field samples of similar lithology. From the interpretation, we compute the TCE and DP concentration distributions in the rock matrix prior to the test to illustrate how the results can be used to enhance understanding of contaminant distribution in the rock matrix.  相似文献   

13.
This work presents a rigorous numerical validation of analytical stochastic models of steady state unsaturated flow in heterogeneous porous media. It also provides a crucial link between stochastic theory based on simplifying assumptions and empirical field and simulation evidence of variably saturated flow in actual or realistic hypothetical heterogeneous porous media. Statistical properties of unsaturated hydraulic conductivity, soil water tension, and soil water flux in heterogeneous soils are investigated through high resolution Monte Carlo simulations of a wide range of steady state flow problems in a quasi-unbounded domain. In agreement with assumptions in analytical stochastic models of unsaturated flow, hydraulic conductivity and soil water tension are found to be lognormally and normally distributed, respectively. In contrast, simulations indicate that in moderate to strong variable conductivity fields, longitudinal flux is highly skewed. Transverse flux distributions are leptokurtic. the moments of the probability distributions obtained from Monte Carlo simulations are compared to modified first-order analytical models. Under moderate to strong heterogeneous soil flux conditions (σ2y≥1), analytical solutions overestimate variability in soil water tension by up to 40% as soil heterogeneity increases, and underestimate variability of both flux components by up to a factor 5. Theoretically predicted model (cross-)covariance agree well with the numerical sample (cross-)covarianaces. Statistical moments are shown to be consistent with observed physical characteristics of unsaturated flow in heterogeneous soils.©1998 Elsevier Science Limited. All rights reserved  相似文献   

14.
An important quantity in groundwater protection is the residence time of water in an aquifer. It relates to both the travel time of a pollutant to arrive at a well and the time span required for self-purification of a polluted aquifer after removal of pollutant inputs. Time scales for aquifers can be gained from artificial tracer experiments or from environmental tracer data, the latter offering the only realistic alternative if time scales of years or decades have to be taken into account.

Different tracers show different time scales due to their different transport mechanisms especially in the unsaturated zone. While solute tracers are moved advectively with the seepage water, gas tracers pass the unsaturated zone diffusively through the air phase. Depending on the properties of the unsaturated zone (hydraulic properties, thickness) this difference in behavior can be used to separate the subsurface transport process into the unsaturated and the saturated parts.

In a field study in Germany, SF6 and 3H were used as environmental tracers. Both have a relatively well-known input function. Interpretation of data from observation wells by a box model approach led to spatially and temporally varying residence times. This was an indication that the influence of the unsaturated zone could not be neglected. While the gas tracer SF6 shows only residence times in the saturated zone, the tracer 3H reflects the whole travel time of water including both the unsaturated and saturated zones. Using a one-dimensional plug-flow model for the unsaturated zone combined with a detailed two-dimensional flow and transport model for the saturated zone leads to a holistic and consistent interpretation of the measured tracer concentrations. The observed pattern of old water under thick loess cover and younger water under areas where the fractured basalt aquifer crops out is reproduced after adjusting only two parameters: the effective porosity of the saturated aquifer and the product of field capacity and thickness of the unsaturated zone. While the effective porosity of the saturated zone is adjusted by means of the SF6 data, the field capacity of the loess layer is adjusted by means of the 3H observations. The thickness of the unsaturated zone is deduced from geological and pedological maps. All flow data are obtained from a calibrated flow model, which is based on geological data, observed heads and pumping tests only.

The transport model for the saturated zone was calibrated by fitting the porosity by means of gaseous tracer concentrations (SF6). The combined saturated–unsaturated zone model was then calibrated by fitting the field capacity of the unsaturated zone by means of 3H concentrations. With this model it was possible to verify the observed NO3 concentrations at the drinking water wells and to develop predictions for their future development under various scenarios of fertilizer input reduction in specific areas.  相似文献   


15.
The partial differential equation which governs the seepage of water in unsaturated and saturated porous media is solved numerically by a generalized Newton iteration technique for two models, one ditch drainage model and one earth dam model. For each model, which is two-dimensional, a few hypothetical soils with different moisture retention curves are considered. In both models only drainage from an initially saturated soil occurs; thus, the problem of hysteresis is avoided. The results of the computations are compared with those of corresponding saturated (pure groundwater) models; solutions obtained earlier by this author and others.

Computational instability phenomena appear when the slope of the retention curves is made steep, i.e., for poorly-graded soils.  相似文献   


16.
The main processes affecting the migration of a solute in a fissured aquifer will be advection and dispersion in the fissures, diffusion into the porous matrix; and adsorption. This paper considers solute transport in an idealized fissured aquifer consisting of slabs of saturated rock-matrix separated by equally spaced, planar fissures. The solution of the transport equations is developed as far as Laplace transforms of the solute concentrations in the fissure and matrix water. Numerical inversion of the transforms is used to investigate characteristic behaviour of the model for a number of special cases.  相似文献   

17.
二维格子气自动机模拟孔隙介质的电传输特性   总被引：14，自引：1，他引：13       下载免费PDF全文
本文应用2-D格子气自动机模拟饱和油水两相的多孔介质的导电特性。在油水两相界面处,引入反射与透射系数来决定粒子的运动状态,通过调节反射与透射系数就可以改变油水两相的导电性差异。用模拟结果考察Archie公式的地层因素F=R0/Rw和电阻率增大系数I=Rt/R0,其中,R0为百分之百饱和水时的岩石电阻率,Rw为水的电阻率,Rt为不同流体饮和度时的岩石电阻率。结果表明F与孔隙度φ间,I与含水饱和度间都存在幂关系,并可以表示为F=aφ-m,I=bS-nw模拟结果同时证实：公式中的参数a、m的变化反映了孔隙微观结构的变化,参数b、n主要受孔隙度大小和油相分布状态的影响。  相似文献   

18.
The macroscopic spreading and mixing of solute plumes in saturated porous media is ultimately controlled by processes operating at the pore scale. Whilst the conventional picture of pore-scale mechanical dispersion and molecular diffusion leading to persistent hydrodynamic dispersion is well accepted, this paradigm is inherently two-dimensional (2D) in nature and neglects important three-dimensional (3D) phenomena. We discuss how the kinematics of steady 3D flow at the pore scale generate chaotic advection—involving exponential stretching and folding of fluid elements—the mechanisms by which it arises and implications of microscopic chaos for macroscopic dispersion and mixing. Prohibited in steady 2D flow due to topological constraints, these phenomena are ubiquitous due to the topological complexity inherent to all 3D porous media. Consequently 3D porous media flows generate profoundly different fluid deformation and mixing processes to those of 2D flow. The interplay of chaotic advection and broad transit time distributions can be incorporated into a continuous-time random walk (CTRW) framework to predict macroscopic solute mixing and spreading. We show how these results may be generalised to real porous architectures via a CTRW model of fluid deformation, leading to stochastic models of macroscopic dispersion and mixing which both honour the pore-scale kinematics and are directly conditioned on the pore-scale architecture.  相似文献   

19.
Double-porosity is an important characteristic of microstructure in a large range of geomaterials. It designs porous media with connected fissures/fractures or aggregated soils. The origin of double-porosity can be natural or/and it can result from mechanical, chemical or biological damage. The presence of double-porosity can significantly affect the behaviour of geomaterials. In this paper we provide an experimental evidence of the double-porosity effects by performing laboratory experiments. Series of tracer dispersion experiments (in saturated and unsaturated steady-state water flow conditions) in a physical model of double-porosity geomaterial were carried out. For the comparative purposes, experiments of the same type were also performed in a singleporosity model medium. The results clearly showed that the double-porosity microstructure leads to the non-Fickian behaviour of the tracer (early breakthrough and long tail) in both saturated and unsaturated cases.  相似文献   

20.
As illustrated variously by wetting and drying scanning curves, flow in unsaturated porous media is inherently nonlocal. This nonlocality is also manifest in hysteresis in the classical Darcy conductivity. It is the authors' belief that most current theories of unsaturated/saturated flow are often inadequate, as they do not account for spatial nonlocality and memory. Here we provide a fundamental theory in which nonlocality of the flow constitutive theory is a natural consequence of force balances. The results are derived from general principles in statistical physics and under appropriate limiting conditions, the classical Darcy's Law is recovered for saturated flow. A notable departure in this theory from other nonlocal flow theories is that a classical Darcy type equation on a small scale need not exist.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号