An analytical solution for describing the transient temperature distribution in an aquifer thermal energy storage system     

Kuang‐Yi Li  Shaw‐Yang Yang  Hund‐Der Yeh 《水文研究》2010,24(25):3676-3688

A mathematical model is developed for predicting the temperature distribution in an aquifer thermal energy storage (ATES) system, which consists of a confined aquifer bounded from above and below by the rocks of different geological properties. The main transfer processes of heat include the conduction and advection in the aquifer and the conduction in the rocks. The semi‐analytical solution in dimensionless form for the model is developed by Laplace transforms and its corresponding time‐domain solution is evaluated by the modified Crump method. Field geothermal property data are used to simulate the temperature distribution in an ATES system. The results show that the heat transfer in the aquifer is fast and has a vast effect on the vicinity of the wellbore. However, the aquifer temperature decreases with increasing radial and vertical distances. The temperature in the aquifer may be overestimated when ignoring the effect of thermal conductivity. The temperature distribution in an ATES system depends on the vertical thermal conduction in the rocks and the horizontal advection and thermal conduction in the aquifer. The present solution is useful in designing and simulating the heat injection facility in the ATES systems. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Chronostratigraphic constraints on Middle Pleistocene faunal assemblages and Acheulian industries from the Cretone lacustrine basin,central Italy          下载免费PDF全文

Fabrizio Marra  Piero Ceruleo  Brian Jicha  Luca Pandolfi  Carmelo Petronio  Leonardo Salari  Biagio Giaccio  Gianluca Sottili 《第四纪科学杂志》2016,31(7):641-658

  相似文献   

Ring foundation on elastic subgrade: an analytical solution for computer modelling using the Lagrangian multiplier method          下载免费PDF全文

Shahin Nayyeri  Masood Hajali  Caesar Abi Shdid 《国际地质力学数值与分析法杂志》2016,40(14):2017-2030

In the practice of geotechnical engineering, the case of a ring footing carrying a set of concentrated point loads is a common problem. At times, the induced vertical and angular displacements for the ring footing need to be evaluated at a relatively precise level. By making use of the governing set of equations derived for the case of a general curved beam, expressions that can be easily implemented in modern computing software are derived for the vertical and angular displacements of a ring footing of rectangular cross section as functions of the radial position. The loading case considered is a vertical point load, and the soil is modelled as elastic. Estimates of the displacements have been shown for a common range of practical applications. The behaviour for a set of concentrated loads may be evaluated using the derived equations through direct superposition. Nonlinear finite element analysis is used to evaluate the vertical deflection and angular twist of the ring foundation. Numerical analysis performed for three ring foundations with different radii and cross sections is reported to validate the accuracy of the derived analytical solution. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

A constitutive model for mechanical and chemo‐mechanical compaction in sedimentary basins and finite element analysis          下载免费PDF全文

A. Brüch  S. Maghous  F.L.B. Ribeiro  L. Dormieux 《国际地质力学数值与分析法杂志》2016,40(16):2238-2270

Compaction and associated fluid flow are fundamental processes in sedimentary basin deformation. Purely mechanical compaction originates mainly from pore fluid expulsion and rearrangement of solid particles during burial, while chemo‐mechanical compaction results from Intergranular Pressure‐Solution (IPS) and represents a major mechanism of deformation in sedimentary basins during diagenesis. The aim of the present contribution is to provide a comprehensive 3D framework for constitutive and numerical modeling of purely mechanical and chemo‐mechanical compaction in sedimentary basins. Extending the concepts that have been previously proposed for the modeling of purely mechanical compaction in finite poroplasticity, deformation by IPS is addressed herein by means of additional viscoplastic terms in the state equations of the porous material. The finite element model integrates the poroplastic and poroviscoplastic components of deformation at large strains. The corresponding implementation allows for numerical simulation of sediments accretion/erosion periods by progressive activation/deactivation of the gravity forces within a fictitious closed material system. Validation of the numerical approach is assessed by means of comparison with closed‐form solutions derived in the context of a simplified compaction model. The last part of the paper presents the results of numerical basin simulation performed in one dimensional setting, demonstrating the ability of the modeling to capture the main features in elastoplastic and viscoplastic compaction. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

A coupled discrete element model for the simulation of soil and water flow through an orifice          下载免费PDF全文

Yao Tang  Dave H. Chan  David Z. Zhu 《国际地质力学数值与分析法杂志》2017,41(14):1477-1493

Soil erosion around defective underground pipes can cause ground collapses and sinkholes in urban areas. Most of these soil erosion events are caused by fluidization of the surrounding soil with subsequent washing into defective sewer pipes. In this study, this soil erosion process is simplified as the gradual washout of sand particles mixed with water through an orifice. The discrete element method is used to simulate the large deformation behavior of the sand particles, and the Darcy fluid model is coupled with this approach to simulate fluid flow through porous sand media. A coupled 3D discrete element model is developed and implemented based on this scheme. To simulate previous experiments using this coupled model considering the current computing capacity, we incorporated a ‘supply layer’ to study the continuous erosion process. The coupled model can predict the erosion flow rates of sand and water and the shape of erosion void. Thus, the model can be used as an effective and efficient tool to investigate the soil erosion process around defective pipes. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

An approach to the modeling of viscoelastic damage. Application to the long‐term creep of gypsum rock materials     

B. Nedjar  R. Le Roy 《国际地质力学数值与分析法杂志》2013,37(9):1066-1078

A three‐dimensional phenomenological model is developed to describe the long‐term creep of gypsum rock materials. The approach is based on the framework of continuum damage mechanics where coupling with viscoelasticity is adopted. Specifically, a local damage model based on the concept of yield surface is proposed and deeply investigated. Among the many possibilities, we choose in this work its coupling with a generalized Kelvin–Voigt rheological model to formulate the whole behavior. Long‐term as well as short‐term relaxation processes can be integrated in the model by means of as many as necessary viscoelastic processes. The numerical discretization is described for an easy integration within a finite element procedure. Finally, a set of numerical simulations is given to show the possibilities of the presented model. It shows good agreement with some experimental results found in the literature. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Water waves generated by instantaneous disturbances at the bed of a sloping beach     

Rumpa Chakraborty 《地球物理与天体物理流体动力学》2013,107(5):481-496

The two-dimensional problem of the generation of water waves due to instantaneous disturbances prescribed at the bed of a beach sloping at an arbitrary angle is studied here. It is formulated in terms of an initial-boundary-value problem for the velocity potential describing the motion in the fluid region assuming the linear theory. Using the Laplace transform in time and the Mellin transform in distance, the problem is reduced to solving a difference equation whose method of solution is of considerable importance in the literature. The form of the free surface is obtained in terms of a multiple infinite integral that is evaluated by the method of steepest-descent. For some prescribed forms of the disturbance at the bed of the beach, the free surface is depicted in a number of figures for different beach angles. It is observed that as the beach angle decreases, the maximum wave height increases, which is plausible.  相似文献   

Calculating actual crop evapotranspiration under soil water stress conditions with appropriate numerical methods and time step     

Songhao Shang 《水文研究》2012,26(22):3338-3343

Calculation of actual crop evapotranspiration under soil water stress conditions is crucial for hydrological modeling and irrigation water management. Results of actual evapotranspiration depend on the estimation of water stress coefficient from soil water storage in the root zone, which varies with numerical methods and time step used. During soil water depletion periods without irrigation or precipitation, the actual crop evapotranspiration can be calculated by an analytical method and various numerical methods. We compared the results from several commonly used numerical methods, including the explicit, implicit and modified Euler methods, the midpoint method, and the Heun's third‐order method, with results of the analytical method as the bench mark. Results indicate that relative errors of actual crop evapotranspiration calculated with numerical methods in one time step are independent of the initial soil water storage in the range of soil water stress. Absolute values of relative error decrease with the order of numerical methods. They also decrease with the number of time step, which can ensure the numerical stability of successive simulation of soil water balance. Considering the calculation complexity and calculation errors caused by numerical approximation for different time step and maximum crop evapotranspiration, the explicit Euler method is recommended for the time step of 1 day (d) or 2 d for maximum crop evapotranspiration less than 5 mm/d, the midpoint method or the modified Euler method for the time step of up to one week or 10 d for maximum crop evapotranspiration less than 5 mm/d, and the Heun's third‐order method for the time step of up to 15 d. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Drained cavity expansions in soils of finite radial extent subjected to two boundary conditions     

Mohammad Pournaghiazar  Adrian R. Russell  Nasser Khalili 《国际地质力学数值与分析法杂志》2013,37(4):331-352

The problem of drained cavity expansion in soils of finite radial extent is investigated. Cylindrical and spherical cavities expanded from zero radius subjected to either constant stress or zero displacement at the finite boundary are considered. The generalised analytical solution procedure presented enables more advanced constitutive models to be implemented than have been possible in previous studies. Results generated for Sydney quartz sand highlight substantial differences between cavity limit pressures for the two boundary conditions and for boundaries of finite and infinite radial extent. This is of significance in accounting for boundary effects when interpreting cone penetration tests conducted in calibration chambers. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Particle manifold method (PMM): A new continuum‐discontinuum numerical model for geomechanics     

L. Sun  G.F. Zhao  J. Zhao 《国际地质力学数值与分析法杂志》2013,37(12):1711-1736

Particle manifold method (PMM) is a new extension of the numerical manifold method (NMM). PMM uses a mathematical cover system to describe the motion and deformation of a particle‐based physical domain. By introducing the concept of particle into NMM, PMM takes the advantages of easy topological and contact operations with particles. In this article, the methodology, formulations and implementation of the method are presented, together with modelling examples for validation. It is found that good solutions for both continuous and discontinuous problems are obtained by the new developed PMM. Due to the underlying coupled continuum‐discontinuum property of PMM, it has great potential for modelling of geomechanical problems. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献