共查询到20条相似文献，搜索用时 187 毫秒
1.
Conventional modelling of transport problems for porous media usually assumes that the Darcy flow velocities are steady. In certain practical situations, the flow velocity can exhibit time‐dependency, either due to the transient character of the flow process or time dependency in the boundary conditions associated with potential flow. In this paper, we consider certain one‐ and three‐dimensional problems of the advective transport of a chemical species in a fluid‐saturated porous region. In particular, the advective flow velocity is governed by the piezo‐conduction equation that takes into account the compressibilities of the pore fluid and the porous skeleton. Time‐ and/or mesh‐refining adaptive schemes used in the computational modelling are developed on the basis of a Fourier analysis, which can lead to accurate and optimal solutions for the advective transport problem with time‐ and space‐dependent advective flow velocity distributions. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
2.
This paper considers advective transport of a soluble contaminant through saturated soil with non-linear sorption of the contaminant onto a stationary porous media. The non-linear sorption isotherms considered in the transport analysis are the Langmuir and Freundlich sorption isotherms. A special case of the Freundlich sorption isotherm is the linear sorption isotherm, and it is shown that in this case transport through a homogeneous soil results in the initial concentration profile simply being translated in the direction of the groundwater flow. However, when the sorption isotherm is non-linear the initial concentration profile distorts as it is translated with the groundwater flow, leading to the development of concentration shock fronts and rarefactions. Analytic solutions to the non-linear first-order hyperbolic equations are developed for a number of contaminant transport problems of practical significance. It is shown that in the case of the Langmuir sorption isotherms, shock fronts develop at the leading edge of the concentration profile while for the Freundlich sorption isotherm shock fronts may develop at either the leading or trailing edge of the concentration profile. Copyright © 1999 John Wiley & Sons Ltd. 相似文献
3.
The present paper develops an analytical approach to the problem of heat-induced moisture movement in the vicinity of a spherical heat source embedded in an undeformable, moist porous solid of infinite extent. A transient-state distribution of temperature within the infinite medium is assumed to induce the moisture transport process. The numerical results, presented in the paper, illustrate the influence of the moisture transport characteristics on the time-dependent distribution of moisture within the porous medium. © 1998 John Wiley & Sons, Ltd. 相似文献
4.
This paper presents a model for the analysis of plane waves diffraction at a cavity in an infinite homogeneous poroelastic saturated medium, lined by a lining composed of four equal segments. An elastic boundary layer is placed between the cavity lining and the infinite porous medium. The boundary layer is simulated by ‘elastic boundary conditions’ in which the bulk matrix stress is proportional to the relative displacement between the lining and the surrounding medium matrix boundary. In addition, fluid impermeability through the intermediate layer is assumed. For the frequencies, that differ from the pseudoresonanse frequencies, the problem was reduced to the problem of an ideal elastic medium. A closed‐form analytical solution of the problem was obtained using Fourier–Bessel series, the convergence of which was proven. It was shown that the number of series terms required to obtain a desired level of accuracy can be determined in advance. The influence of the medium porosity on the medium dynamic stress concentration was studied. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
5.
The leakage effect in porous fissured media has been considered in a general sense by introducing a new expression of the leakage term in this paper. The double porosity concept is employed and the related expressions are formulated using the upwind finite element approach. Considering the infinite extension of the problem domain, a mapped transient infinite element has been presented to simulate the far field of the infinite medium. Since the mass transfer function of the present mapped transient infinite element is dependent on both space and time variables, the mechanism of transient contaminant migration problems in infinite porous fractured media can be rigorously simulated because the property matrices of the element are evaluated at any time of interest. By comparing the current numerical results with the analytical ones, the accuracy, correctness and effectiveness of the present method have been established. Three different time discretization schemes were examined and it was found that either the central difference or the backward difference approximation is suitable for the upwind finite element simulation of transient contaminant migration problems. 相似文献
6.
Using the numerical model presented in the first paper of this research,1 a parametric study has been carried out in this paper to investigate the effect of several important parameters on the transient contaminant transport in infinite porous fractured media. From the related numerical results, it has been demonstrated that: (1) transmissive coefficient between the porous block and the fissured network has a significant influence on the value of the concentration but has little effect on the speed of contaminant transport; (2) porosities in the porous block and fissured network have a significant influence on the maximum value of the concentration; (3) average linear velocity of flow has a significant influence on both the concentration distribution and speed of contaminant transport; (4) dispersion coefficient of the medium affects not only the shape of the concentration versus time curve but also the peak value of the concentration. 相似文献
7.
V. M. Murthy 《Journal of Earth System Science》1999,108(4):317-320
Surface processes play an important role in the simulation of desertification and climate change. The present study shows
that enhancement of evaporation at a given place is due to surface temperature and wetness inhomogeneities compared to its
surroundings. Between any adjacent dry and wet strips the induced advective heat transport changes the available net radiation
between sensible and latent heat fluxes. It also accommodates a redistribution of surface energy between two adjacent inhomogeneous
surface strips. A combined model of Tarpley (1994) and Ya Guo and Schuepp (1994) about the advective heat transport over dry
and wet strips has been tested over Andhra Pradesh. Taking the Ananthpur district as a relatively drier region, the induced
transport over the rest of the districts is estimated. This study has been done for two monsoon seasons of 1990 and 1991.
Attempts are made to identify an indexing criteria based on this study to estimate the magnitude of the heat transport. 相似文献
8.
This study presents an exact solution for drained spherical cavity expansion subjected to constant stress and zero displacement boundary conditions in finite medium. The solution is exact as no simplification is involved in the solution procedure in the plastic region and can be solved as an initial value problem. The effect of finite radial extent on the results of cavity expanding to a certain radius is accounted for through the initial values at the elastic-plastic boundary. A critical state constitutive model with a nonassociated flow rule is adopted. The model parameters from literature are then used in generating the results for cavity expansion in infinite and finite radial extent to highlight the difference. Also, the results are presented in a way that can be used to account for the boundary effect for the interpretation of cone penetration tests conducted in calibration chambers. 相似文献
9.
山东淄博市大武水源地裂隙岩溶水中污染物运移的数值研究 总被引：11，自引：0，他引：11
在分析研究淄博市大武水源地裂隙岩溶含水层的水力性质和污染物运移特征的基础上 ,对裂隙岩溶水的水头和污染物运移进行数值研究。目前国内外对裂隙岩溶水进行数值计算时 ,通常用等价多孔介质模型 ,但裂隙岩溶介质和多孔介质有很大不同。裂隙岩溶介质的储水和导水空间为裂隙网络 ,导水系数大 ,地下水的实际平均流速比孔隙水大得多 ,但给水度和贮水系数小。当用等价多孔介质模型进行模拟时应考虑这些特点。对于污染物运移的模拟 ,要同时求解水头方程和对流弥散方程 ,可采用MODFLOW和MT3D软件进行模拟。研究区裂隙岩溶水水头的数值计算表明 ,等效多孔介质模型水头的拟合误差能满足国标GB/T144 97- 93的要求。各时段地下水水量均衡计算的精度也满足要求。对流弥散方程的数值计算 ,由于Peclet数高达 95 .6 7,对流占绝对优势 ,可能存在数值弥散和数值振荡 ,因而采用多种方法进行了比较。对于同一问题 ,同时采用上游有限差分法 (UFDM) ,混合的欧拉拉格朗日方法 (特征线法MOC、改进特征线法MMOC和混合特征线法HMOC) ,总变异消减法(TVD)进行计算 ,并比较其结果。结果表明 ,混合特征线法 (HMOC)和总变异消减法 (TVD)比较适合于对流占优势的运移问题计算。由于渗透系数K和有效孔隙度θ对溶质运移结果的影响很大 ,? 相似文献
10.
D. W. Smith 《国际地质力学数值与分析法杂志》2000,24(8):693-722
Contaminant migration through soil is usually modelled mathematically using the dispersion–advection equation. This type of model finds application when planning the remediation of contaminated land, predicting the movement of polluted groundwater and designing engineered landfills. Usually the analysis assumes that the porous media through which the contaminant migrates is stationary. However, the construction of landfills on clay soils means that the soil beneath the landfill will undergo time‐dependent deformation as the soil consolidates. To date, there are no published data on the effect a deforming porous media may have on contaminant transport beneath a landfill; indeed, there appears to be no theory of contaminant migration through a deforming soil. In this paper, a one‐dimensional theory of contaminant migration through a saturated deforming porous media is developed based on a small and large strain analysis of a consolidating soil and conservation of contaminant mass. By selection of suitable parameters, the new transport equation reduces to the familiar one‐dimensional dispersion–advection equation for a saturated soil with linear, reversible, equilibrium controlled sorption of the contaminant onto the soil skeleton. Analytic solutions to a quasi‐steady‐state contaminant transport problem for a deforming media are presented, and a preliminary assessment made of the potential importance of soil deformation on the results of a contaminant migration analysis. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
11.
根据饱和多孔介质中考虑释放效应的一维渗透作用的可溶性污染物的迁移控制方程,通过Laplace变换和Fourier变换及其逆变换求得相应的通解形式。根据半无限体表面可溶性污染物点源注入情形下的基本解,通过积分方法得到半无限体表面圆形区域上作用循环浓度污染物后,在多孔介质内部污染物迁移过程的求解方法。作为一个典型算例,对渗透作用下循环浓度污染源引起的多孔介质中的迁移过程进行分析。算例表明,当污染源浓度为周期变化时,多孔介质内部污染物的浓度随时间增长由不稳定的周期变化过程逐渐过渡到稳定的循环过程,而某一深度处污染物浓度的相位则相应滞后。此时,稳定后的污染物浓度周期与污染源浓度周期相同。实际上,随多孔介质表面污染源浓度的周期变化,在靠近多孔介质表面一定深度范围内的污染物浓度在深度方向也呈增大或减小的交替变化过程。另一方面,在污染源循环变化过程中污染物不断向深处推进,而最终其影响范围限定在某一深度内。 相似文献
12.
The goal of this paper is to present an analytical solution to have a first insight of the impact of ice formation on the surrounding porous rock on underground cavities like reservoir, pipes, tunnels or wellbores. Among the other analytical solutions found in the literature on this topic, the originality of this work resides in the rigorous theoretical framework of poromechanics, which considers the coupling between liquid water and ice crystal under thermodynamic equilibrium. Liquid water transport, thermal conduction, and elastic properties of the phases are also considered. Two analytical solutions are presented, based on a linearization of the system of governing equations. The first one deals with a spherical cavity within an infinite porous medium leading to an exact analytical solution. It allows validating the Stehfest’s algorithm on the numerical inversion of Laplace Transform, used in the second analytical solution, which considers a cylindrical excavation. The validity of this solution is assessed by comparing its results to that issued from a numerical resolution of the nonlinear system of equations. The analytical solution is then ultimately used to identify the influence of key parameters like the thermal/hydraulic conductivities, the amount of ice formed and the thermal dilatation coefficients on the mechanical response of a cylindrical cavity submitted to an internal frost. 相似文献
13.
Fatiha Bouchelaghem 《Mathematical Geosciences》2011,43(1):55-73
This paper presents an analytical study of fluid flow in a porous medium presenting pores of two different length scales:
at the smallest or microscopic scale, the presence of connected voids confers a porous medium structure to the material investigated,
while at the upper or mesoscopic scale, macro-pores are present. This microstructure is employed to represent the progressive
opening of inter-aggregate pore spaces observed in natural compacted montmorillonites polluted by heavy metal ions. Three-dimensional
analytical expressions are rigorously derived for pore fluid velocity and excess pore fluid pressure within the porous matrix,
around an occluded ellipsoidal inter-aggregate void. The eccentricity ratio is employed to characterize the geometrical shape
of the ellipsoidal void, while its orientation with respect to the inflow in the far field is determined by the dip angle
θ. As an application, we investigate the flow focusing effect for varying eccentricity ratios and dip angles. 相似文献
14.
A previously developed two-dimensional numerical model is further developed for simulating the transport of dissolved contaminants originating from dissolution of a coal tar pool in a stratified, saturated porous medium. The model is used to simulate contaminant transport resulting from a rectangular-shaped coal-tar-pool dissolution experiment conducted in a large-scale experimental aquifer. The experimental porous medium consists of two sand strata, a high-hydraulic-conductivity upper stratum and a low-hydraulic-conductivity bottom stratum. The experiment was conducted to a time of 354 days and the groundwater velocity was changed several times within this duration. Model simulations show good agreement against observed contaminant concentrations, and simulations show that dissolved solute below the pool migrated deeper into the bottom stratum as compared to the upper stratum. Furthermore, simulations also suggest that contaminant concentrations in the lower stratum never reached quasi steady-state during the experimental time frame. 相似文献
15.
This paper presents a new numerical tool to model the two-dimensional contaminant transport through saturated porous media
using a meshfree method, called radial point interpolation method (RPIM) with polynomial reproduction. In RPIM, an approximate
solution is constructed entirely in terms of a set of nodes and no characterisation of the interrelationship of the nodes
is needed. The advection–dispersion equation with sorption is considered to illustrate the applicability of the RPIM. The
Galerkin weak form of the governing equation is formulated using 2D meshfree shape functions constructed using thin plate
spline radial basis functions. MATLAB code is developed to obtain the numerical solution. Three numerical examples are presented
and the results are compared with those obtained from the finite element method and analytical solutions. In order to test
the practical applicability and performance of the RPIM, two case studies of contaminant transport through landfill liners
are presented. A good agreement is obtained between the results of the RPIM and the field investigation data. 相似文献
16.
Mariusz Kaczmarek 《国际地质力学数值与分析法杂志》2001,25(8):757-770
In this paper a chemically induced deformation of porous material taking place during advective–dispersive transport of a chemical is considered. Linearized governing equations are derived and analytical solutions of 2 one‐dimensional problems for a homogeneous layer with drained boundaries are developed. Numerical results for a particular clayey material and a chemical migrating through the layer showing distributions of concentration of chemical, changes in porosity of the material and pore fluid pressure, and evolution of settlement of the layer as functions of time are discussed. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
17.
非平衡-非线性吸附情况下填埋场污染物运移分析 总被引：10，自引：0，他引：10
采用Langmuir等温吸附线方程描述非线性吸附性能, 基于改进的混合元方法，通过数值计算与分析，探讨了非平衡吸附条件下污染物运移过程及其机理。计算结果表明: 当考虑非平衡、非线性吸附性能时, 污染物穿透曲线即浓度的时程变化曲线尖锐而狭窄、峰值点前移, “拖长尾”现象不明显, 由此说明，土颗粒对污染物的非平衡、非线性吸附使得污染物的穿透能力增强, 滞留能力下降。进一步的变动参数比较分析表明：Langmuir等温线方程中的参数B、压实粘土衬里的渗透性及地下水渗流速度对污染物运移过程具有显著的影响。 相似文献
18.
A double-porosity model is used to describe the oscillatory gas motion and associated contaminant transport induced by cyclical variations in the barometric pressure at the surface of a fractured porous medium. Flow along the fractures and within the permeable matrix blocks is locally one-dimensional. The interaction between fractures and blocks includes seepage of fluid as well as diffusion of contaminant. To guard against artificial numerical diffusion, the FRAM filtering remedy and methodology of Chapman is used in calculating the advective fluxes along fractures and within blocks. The entire system of equations, including the fracture-matrix interaction terms, is solved by a largely implicit non-iterative algorithm which remains stable and conservative even when the computational time step is large compared to the cross-block transit time of pressure waves. The numerical accuracy is tested by comparison with exact solutions for oscillatory and unidirectional flows, some of which include diffusion interaction between the fracture and the matrix. The method is used to estimate the rate of vertical transport of radioactive gases through the rubblized chimney produced by an underground nuclear explosion. 相似文献
19.
Analytical solutions for advective–dispersive solute transport in double‐layered finite porous media
Peter John Cleall 《国际地质力学数值与分析法杂志》2011,35(4):438-460
Analytical solutions for advection and dispersion of a conservative solute in a one‐dimensional double‐layered finite porous media are presented. The solutions are applicable to five scenarios that have various combinations of fixed concentration, fixed flux and zero concentration gradient conditions at the inlet and outlet boundaries that provide a wide number of options. Arbitrary initial solute concentration distributions throughout the media can be considered via explicit formulations or numerical integration. The analytical solutions presented have been verified against numerical solutions from a finite‐element‐based approach and an existing closed‐form solution for double‐layered media with an excellent correlation being found in both cases. A practical application pertaining to advective transport induced by consolidation of underlying sediment layers on contaminant movement within a capped contaminated sediment system is presented. Comparison of the calculated concentrations and fluxes with alternative approaches clearly illustrates the need to consider advection processes. Consideration of the different features of contaminant transport due to varying pore‐water velocity fields in primary consolidation and secondary consolidation stages is achieved via the use of non‐uniform initial concentration distributions within the proposed analytical solutions. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
20.
Solute transport through a porous medium is typically modelled assuming the porous medium is rigid. However, many applications exist where the porous medium is deforming, including, municipal landfill liners, mine tailings dams, and land subsidence. In this paper, mass balance laws are used to derive the flow and transport equations for a deforming porous medium. The equations are derived in both spatial and material co‐ordinate systems. Solute transport through an engineered landfill liner is used as an illustrative example to show the differences between the theory for a rigid porous medium, and small and large deformation analysis of a deforming porous medium. It is found that the large deformation model produces shorter solute breakthrough times, followed by the small deformation model, and then the rigid porous medium model. It is also found that it is important to include spatial and temporal void ratio variations in the large deformation analysis. It is shown that a non‐linear large deformation model may greatly reduce the solute breakthrough time, compared to a standard transport analysis typically employed by environmental engineers. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献