共查询到20条相似文献，搜索用时 31 毫秒
1.
H. R. Thomas 《国际地质力学数值与分析法杂志》1985,9(6):573-588
A theoretical model is presented to predict simultaneous transient coupled heat and moisture transfer in partly saturated soils. The formulation is in terms of volumetric moisture content, is two dimensional, includes gravitational flow and takes into account latent heat of vaporization effects. The numerical solution of the problem is accomplished by means of a finite element solution algorithm. Predictions from the numerical model are used to investigate the importance of gravitational flow, for the case of a soil stratum subjected to evaporation losses at the surface. The results achieved show good qualitative agreement with expected behaviour. 相似文献
2.
H. R. Thomas 《国际地质力学数值与分析法杂志》1988,12(1):31-44
The subject of this paper is the nonlinear analysis of heat and moisture transfer in partly saturated soil. In particular, an extension of the previous work of the author to two-dimensional applications is presented. Problems of greater complexity and engineering significance can hence be solved. The theoretical basis of the problem is first outlined before the two-dimensional formulation of the numerical solution is given. The finite element method is employed to give the spatial distributions of the variables, with a finite difference scheme being used to predict the temporal variations. The method allows the nonlinear nature of the soil parameters to be modelled, predictions being given of the coupled transient simultaneous transfer of heat and volumetric moisture content. The two-dimensional nature of the work is illustrated for the case of surface evaporative soil moisture losses from sloping ground. Subsoil conditions of a uniform deposit of loam are considered with values of representative material parameters assumed after a literature search. The coupled nature of heat and mass transfer is illustrated by means of an analysis of the moisture migration patterns that occur. Two-dimensional moisture flow, out of and into the soil, is shown to take place. 相似文献
3.
A parallel numerical finite difference model, employing the self-implicit method, for coupled heat and moisture transfer in unsaturated soil is presented. The model is programmed in Occam and executed on a parallel computing network of transputers. An assessment of the model was achieved via the simulation of a laboratory experiment. A very good correlation between experimental and numerical results was obtained. Comparison of results with those obtained from a parallel explicit method is also illustrated showing no significant difference. The computational time employing the new method was, however, found to be half of that obtained using the explicit method. The computational efficiency of the approach was also found to be very high. © 1997 John Wiley & Sons, Ltd. 相似文献
4.
非饱和土热湿耦合传输问题的数值解——理论及一维问题的解 总被引：1，自引：0，他引：1
本文介绍了求解非饱和土壤中热量和水分耦合传输问题的一种数值方法——积分有限差分方法(IFDM)。基于菲利普-迪弗瑞斯(Philip-De Vries)多孔介质热湿耦合流动模型,采用积分有限差分方法,编制了求解热湿传输问题的计算程序 HM1,该程序可用于求解各类工程中遇到的多孔介质一维传热传湿问题。文中还给出了计算结果与实测结果相比较的实例。 相似文献
5.
The results of a numerical simulation of the observed field-drying behaviour of an unsaturated clay are described. In particular, the simulation of soil drying during the summer of 1983 is addressed. Richards' pressure-based unsaturated flow formulation is adopted. Transient numerical solutions are achieved by making use of the finite element method coupled with a finite difference timestepping scheme. The application of the algorithm to the simulation of the problem in hand is described. The initial moisture content distribution represents saturated conditions as indicated by field-measured data. Boundary conditions are determined from an interpretation of meteorological data. A simple method of representing plant root extraction of soil moisture is proposed. Comparisons of numerical results with field-measured data collected by British Gas are presented. Good correlation is obtained. The numerical simulation is shown to be capable of representing the observed field-drying behaviour. The complexity of modelling soil drying, as opposed to infiltration, is clearly illustrated by the work presented. 相似文献
6.
A parallel numerical model, employing a finite difference explicit scheme for the analysis of coupled heat and moisture transfer in unsaturated soil, is employed to simulate a laboratory experiment of heating of medium sand. The model, written in a two-dimensional polar co-ordinate formulation, is programmed in the concurrent language Occam and executed on a parallel computing network of transputers. Parallelization is adopted as a means of overcoming computing difficulties, which limited numerical solutions to those at steady state, to enable transient behaviour to be simulated. The parallel algorithm was found to be very efficient, enabling a full solution of transient behaviour to be obtained. An investigation of the ability of the model to accurately simulate the complex, interrelated coupled nature of both two-dimensional transient and steady-state behaviour yielded very good correlation between experimental and numerical results. It can therefore be concluded that overall the results obtained provide confidence in the validity of the approach proposed. 相似文献
7.
8.
On the development of a model of the thermo-mechanical-hydraulic behaviour of unsaturated soils 总被引：3，自引：0，他引：3
The development of a model of the thermo-mechanical-hydraulic behaviour of unsaturated soil is described. A step by step approach is adopted, taking as a starting point a potential-based model of coupled heat and moisture transfer. Extensions of this work to include the effect of air transfer and the deformation characteristics of the soil are presented. Constitutive relationships which accommodate non-linear elasticity, thermoelasticity and elasto-plasticity are considered. Numerical solutions of the various versions of the model are presented. Validation of the versions of the model, at their differing stages of development, is considered in turn and results presented in support of each case. The need for further, on-going research is also suggested, with a view towards resolving remaining problems of concern in safety assessment. 相似文献
9.
An analysis of coupled heat and moisture movement in unsaturated soil in terms of the fundamental potentials for flow is examined. The approach adopted is based on the assumption that the total potential for liquid flow consists of two components, the elevation and the capillary potential. The fundamental potentials employed in the work are, therefore, temperature and capillary potential. The full theoretical formulation of the problem is presented, together with full details of the solution algorithm employed. Spatial discretization is achieved via the use of the finite element method, with the time-varying behaviour described by a finite difference technique. Soil parameter variations as functions of both temperature and moisture content are included in a one-dimensional approach. The work is applied to a practical engineering problem, namely heat and mass transfer in the upper layers of a soil stratum. This problem is of importance to the utilities, since many services are buried in this zone. Material parameters obtained from an associated programme of experimental work are employed. Moisture content and temperature profiles indicating the extent and rate of penetration of drying and heating fronts are produced. 相似文献
10.
Nadia Laredj Hanifi Missoum Karim Bendani Mustapha Maliki 《Arabian Journal of Geosciences》2011,5(5):935-942
Knowledge of transport processes of heat and moisture in soils of arid zones is vital to understanding the environmental and economic impacts of many activities: agriculture, waste disposal, geoenvironmental practices and earth sciences. Through extensive review and study on the different aspects of coupled transfer processes in swelling porous media, a general mathematical model for coupled heat, moisture, air flow and deformation problems in clayey soils is proposed in a consistent and unified manner. The model is characterized by the presence of a deformable solid matrix filled with two fluid phases (liquid water and air). In the proposed model, both pore water and air transfers are assumed to be governed by the generalized Darcy’s law. Fully coupled, non-linear partial differential equations are established and then solved by using a Galerkin weighted residual approach in space domain and an implicit integrating scheme in time domain. The obtained model has been finally validated by means of some case tests for the prediction of the thermo-hydro-mechanical behaviour of unsaturated swelling soils. The calculated relative errors between experimental and numerical results are 3% for temperature and 7% for stresses. Consequently, the developed numerical model predicts satisfactory results, compared to experimental test measures. The model is applicable to two-dimensional problems with various initial and boundary conditions; non-linear soil parameters can be easily included in this model. 相似文献
11.
In this paper, a fully coupled thermo-hydro-mechanical model is presented for two-phase fluid flow and heat transfer in fractured/fracturing porous media using the extended finite element method. In the fractured porous medium, the traction, heat, and mass transfer between the fracture space and the surrounding media are coupled. The wetting and nonwetting fluid phases are water and gas, which are assumed to be immiscible, and no phase-change is considered. The system of coupled equations consists of the linear momentum balance of solid phase, wetting and nonwetting fluid continuities, and thermal energy conservation. The main variables used to solve the system of equations are solid phase displacement, wetting fluid pressure, capillary pressure, and temperature. The fracture is assumed to impose the strong discontinuity in the displacement field and weak discontinuities in the fluid pressure, capillary pressure, and temperature fields. The mode I fracture propagation is employed using a cohesive fracture model. Finally, several numerical examples are solved to illustrate the capability of the proposed computational algorithm. It is shown that the effect of thermal expansion on the effective stress can influence the rate of fracture propagation and the injection pressure in hydraulic fracturing process. Moreover, the effect of thermal loading is investigated properly on fracture opening and fluids flow in unsaturated porous media, and the convective heat transfer within the fracture is captured successfully. It is shown how the proposed computational model is capable of modeling the fully coupled thermal fracture propagation in unsaturated porous media. 相似文献
12.
A poroelastic numerical model is presented to evaluate three-dimensional consolidation due to groundwater withdrawal from desaturating anisotropic porous media. This numerical model is developed based on the fully coupled governing equations for groundwater flow in deforming variably saturated porous media and the Galerkin finite element method. Two different cases of unsaturated aquifers are simulated for the purpose of comparison: a cross-anisotropic soil aquifer, and a corresponding isotropic soil aquifer composed of a geometrically averaged equivalent material. The numerical simulation results show that the anisotropy has a significant effect on the shapes of three-dimensional hydraulic head distribution and displacement vector fields. Such an effect of anisotropy is caused by the uneven partitioning of the hydraulic pumping stress between the vertical and horizontal directions in both groundwater flow field and solid skeleton deformation field. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
13.
This paper focuses attention on the development of a numerical model of the hydro/thermo/mechanical behaviour of unsaturated clay and its consequent verification and validation. The work presented describes on-going collaboration between the Cardiff School of Engineering and Atomic Energy of Canada. The model development, which was carried out at Cardiff, can be described as being based on a mechanistic approach to coupled heat, moisture and air flow. This is then linked to a deformation analysis of the material within a ‘consolidation’ type of model. The whole is solved via the finite element method to yield a computer software code named COMPASS (COde for Modelling PArtly Saturated Soil). Some aspects of verification and validation of the model have been addressed in-house. However, the purpose of current AECL work is to provide an independent, rigorous, structured programme of validation and the paper will also explore the further validation of COMPASS within this context. © 1998 by John Wiley & Sons, Ltd. 相似文献
14.
冻土水热力耦合研究现状及进展 总被引：5，自引：4，他引：1
开展冻土水热力三场耦合研究对解决寒区工程问题具有重要的理论指导意义。归纳了冻土水热力耦合的理论基础,认为目前的水分迁移驱动力假说仍然不能很好地解释水分迁移现象,分凝冰的形成机制及判据仍需进行深入的研究。分类和评价了常见的正冻土水热力耦合模型,发现流体动力学模型虽然能够很好地描述水热迁移现象,但未考虑非连续冰透镜体;而较复杂的刚冰模型虽然考虑了冻结缘内水热迁移耦合现象,但是参数众多;热力学模型从微观角度描述了冻土水热力并考虑孔隙吸力,但仍存在参数众多的问题。同时,对预融膜理论在冻土水热力耦合问题中的应用进行了分析和展望,认为可以借助预融膜理论对冻土水热力耦合中的能量、水分迁移驱动力以及迁移速率等进行描述。最后,基于冻土水热力三场耦合研究现状及存在的问题,提出了冻土水热力耦合研究的总体构想：研究与实际情况相符同时适用于稳态及非稳态的通用数学表达式,开展冻土物理学各个参数的动态变化研究,纳入非饱和土体在冻融过程中的水热力相互作用研究,实现水热力在真正意义上的耦合,同时,加强预融膜理论在大尺度、陆面过程以及水热边界等方面的应用研究。 相似文献
15.
Vegetation contributes to weak soil stabilisation through reinforcement of the soil, dissipation of excess pore pressure and increasing the shear strength by induced matric suction. This paper describes the way vegetation influences soil matric suction, shrinkage and ground settlement in the vadose zone through transpiration. A mathematical model for the rate of root water uptake, including the root growth rate considering ground conditions, type of vegetation and climatic parameters, has been developed. A finite element approach is employed to solve the transient coupled flow-deformation equations. The finite element mesh is built using partially saturated soil elements capable of representing the salient aspects of unsaturated permeability and the soil water characteristic curve. The model formulation is based on the effective stress theory of unsaturated soils. Based on this proposed model, the distribution of the ground matric suction profile adjacent to the tree is numerically analysed. Current field measurements of soil matric suction and moisture content collected from Miram site located in Victoria State, Australia by the authors are compared with the numerical predictions. The results indicate that the proposed root water uptake model incorporated in the numerical analysis can be used for prediction of ground properties influenced by tree roots. 相似文献
16.
Prediction of unsaturated soil behavior during earthquake loading has received increasing attention in geotechnical engineering research and practice in recent years. Development of a fully coupled analysis procedure incorporating a coupled hydromechanical elastoplastic constitutive model for dynamic analysis of unsaturated soils has, however, been limited. This paper presents the implementation of a coupled hydromechanical elastoplastic constitutive model into a fully coupled dynamic analysis procedure and its validation using a centrifuge test. First, the fully coupled finite element equations governing the dynamic behavior of unsaturated soils with the solid skeleton displacement, pore water pressure, and pore air pressure as nodal unknowns are briefly presented. The closest point projection method is then utilized to implement the coupled hydromechanical elastoplastic constitutive model into the finite element equations. The constitutive model includes hysteresis in soil–water characteristic curves, cyclic elastoplasticity of the solid skeleton, and the coupling mechanisms between the SWCCs and the solid skeleton. Finally, the analysis procedure is validated using the results from a dynamic centrifuge test on an embankment constructed of compacted unsaturated silt subjected to base shaking. Reasonable comparisons between the predicted and measured accelerations, settlements, and deformed shapes are obtained.
相似文献17.
In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic properties of the porous medium, i.e. saturation, relative permeability and capillary pressure. It is worth mentioning that the imposition of various boundary conditions is feasible notwithstanding the choice of the primary variables. The modified Pastor–Zienkiewicz generalized constitutive model is introduced into the mathematical formulation to simulate the mechanical behavior of the unsaturated soil. The accuracy of the proposed mathematical model for analyzing coupled fluid flows in porous media is verified by the resolution of several numerical examples for which previous solutions are known. Finally, the performance of the computational algorithm in modeling of large-scale porous media problems including the large elasto-plastic deformations is demonstrated through the fully coupled analysis of the failure of two earth and rockfill dams. Furthermore, the three-phase model is compared to its simplified one which simulates the unsaturated porous medium as a two-phase one with static air phase. The paper illustrates the shortcomings of the commonly used simplified approach in the context of seismic analysis of two earth and rockfill dams. It is shown that accounting the pore air as an independent phase significantly influences the unsaturated soil behavior. 相似文献
18.
在考虑相变的热能平衡方程和非饱和水分迁移质量控制方程的基础上,建立温度场-水分场的耦合模型,并采用一种无网格粒子算法（SPH）进行数值求解。其中,耦合方程中考虑了水流传热以及温度势对水流的直接驱动,在不考虑相变的情况下,该耦合模型可退化为常温下的水-热耦合模型,故可用于模拟冻融循环的相关问题。从求解热能平衡方程中的含冰量出发,实现解耦并对半无限单向冻结条件下介质内非稳态温度场和体积含水率分布场进行模拟,将耦合作用下的温度场与不耦合的解析解进行对比,反映出水分迁移对温度场存在较大影响。最后,求解了路基边坡在季节性周期温度边界下,温度场、水分场分布的演变规律,并评估了边坡阴阳面受热不均对水热两场分布的影响。计算结果基本能反映土冻结相变的实际物理过程,光滑粒子算法可以用于尝试解决冻土领域的其他相关问题。 相似文献
19.
Energy piles are bi-functional foundation elements used as structural support as well as ground heat exchangers for shallow geothermal energy systems. Because they are relatively short, energy piles may be partially embedded in unsaturated soils. Saturation conditions influence the thermal properties of the ground and therefore the heat exchange rate, which in turn affects the efficiency of energy piles. This paper combines analytical, experimental and numerical investigations to evaluate the heat exchange rate of energy piles partially or fully embedded in unsaturated soils. The proposed analytical solution is based on the cylindrical heat source theory that treats the soil as a semi-infinite, homogeneous, and isotropic medium. The solution from this theory is multiplied by a function developed analytically in this paper and the outcome is the heat exchange rate for energy piles in unsaturated soils. The proposed function depends on soil saturation, soil and pile thermal properties, and pile geometry. The analytical solution was compared against a finite element solution; which was in turn validated against results from laboratory experiments. Very satisfactory agreements between the analytical, numerical and experimental outputs were observed. The proposed method can be used for a quick and simple evaluation of the efficiency of energy piles in unsaturated conditions. The proposed analytical solution can also be a useful tool for the verification of numerical codes developed for the design of energy piles in unsaturated soils. 相似文献
20.
This paper presents the results of the application of a numerical model of flow of water in unsaturated soil to simulate real field behaviour. Numerical predictions are compared with field-monitored results as part of an assessment of the approach adopted. The numerical approach is based on a finite element solution of Richards' theoretical formulation, adopting a finite difference recurrence relationship to model the transient nature of the problem. The field results have been collected by British Gas in Kimmeridge clay at a site in Swindon. The simulation addressed the autumn wetting of the soil during 1983. Comparisons of volumetric moisture content readings from the field and numerical predictions give good correlation. It is concluded that the numerical model has adequately represented field behaviour over the depth of interest. Taken in conjunction with previous results, it is claimed that the model should now be capable of use in a predictive mode. 相似文献