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.  相似文献   

Engineering characteristics and uses of duricrusts in Australia     

G. H. McNally 《Australian Journal of Earth Sciences》2013,60(6):535-547

The uses and shortcomings of duricrusts (ferricrete, calcrete and silcrete) in engineering construction (as used for road‐base, aggregate, foundation materials and aquifers) are reviewed. Australian production of these materials represents about one‐third of all unprocessed road‐base and they are especially important as pavement courses for lightly trafficked, low‐cost rural and outback roads. However, duricrusts are regarded as marginal materials at best because of their typically poor grading, particle unsoundness, high fines plasticity, and absorption of water and bitumen. These materials are used because they are available locally in areas that are otherwise lacking in hard rock materials, such as deeply weathered and sedimentary rock terrains. The weathering profiles of which they form part are characterized by high permeability (despite being clay‐rich), variable cementation, low compressibility and a tendency to become weaker with depth. Although the more indurated layers are unrippable, they are also difficult to blast. Pedogenic (nodular) duricrusts were formerly sought for natural road‐base, because they occur widely and require only rudimentary processing. However, well‐cemented groundwater (vadose) duricrusts are now the preferred deposits, even though they require crushing and screening. Ferricrete is the most widely exploited of the duricrusts for engineering purposes, especially in northern and southwestern Australia, although calcrete is important in South Australia and in the Murray Basin. Silcrete is only a minor source of aggregate and road‐base, mainly in western Queensland.  相似文献   

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.  相似文献   

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.  相似文献   

Coupled finite‐element simulation of injection well testing in unconsolidated oil sands reservoir     

Bin Xu  Ron C.K. Wong 《国际地质力学数值与分析法杂志》2013,37(18):3131-3149

This paper presents a finite‐element (FE) model for simulating injection well testing in unconsolidated oil sands reservoir. In injection well testing, the bottom‐hole pressure (BHP) is monitored during the injection and shut‐in period. The flow characteristics of a reservoir can be determined from transient BHP data using conventional reservoir or well‐testing analysis. However, conventional reservoir or well‐testing analysis does not consider geomechanics coupling effects. This simplified assumption has limitations when applied to unconsolidated (uncemented) oil sands reservoirs because oil sands deform and dilate subjected to pressure variation. In addition, hydraulic fracturing may occur in unconsolidated oil sands when high water injection rate is used. This research is motivated in numerical modeling of injection well testing in unconsolidated oil sands reservoir considering the geomechanics coupling effects including hydraulic fracturing. To simulate the strong anisotropy in mechanical and hydraulic behaviour of unconsolidated oil sands induced by fluid injection in injection well testing, a nonlinear stress‐dependent poro‐elasto‐plastic constitutive model together with a strain‐induced anisotropic permeability model are formulated and implemented into a 3D FE simulator. The 3D FE model is used to history match the BHP response measured from an injection well in an oil sands reservoir. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Palaeoproterozoic mid‐basin inversion in the northern Mt Isa terrane,Queensland     

P. G. Betts 《Australian Journal of Earth Sciences》2013,60(5):735-748

The Mellish Park Syncline is located in the northern part of the Mt Isa terrane. It has an axial trace that transects the remnants of the unconformity‐bounded Palaeoproterozoic Leichhardt and Isa Superbasins. The syncline is separated into a lower and upper component based upon variation in fold geometry across the basin‐bounding unconformity. The lower syncline, in the Leichhardt Superbasin, is tight and has an inclined west‐dipping axial plane. The upper syncline, in the Isa Superbasin, is open and upright. The geometry of the lower syncline is a consequence of a period of shortening and basin inversion which post‐dated the Leichhardt Rift Event (ca 1780–1740 Ma) and pre‐dated the Mt Isa Rift Event (ca 1710–1655 Ma), forming an open and upright north‐oriented syncline. Subsequent southeast tilting and half‐graben development during the Mt Isa Rift Event resulted in the lower syncline being tilted into its inclined geometry. Sequences of the Isa Superbasin were then deposited onto the eroded syncline. The geometry of the upper syncline reflects regional east‐west shortening during the Isan Orogeny (ca 1590–1500 Ma). The position of the upper syncline was largely controlled by the pre‐existing lower syncline. At this time the lower syncline was reactivated and tightened by flexural slip folding.  相似文献   

Analytical and computational procedure for solving the consolidation problem of layered soils     

S. Sadiku 《国际地质力学数值与分析法杂志》2013,37(16):2789-2800

One‐dimensional consolidation analysis of layered soils conventionally entails solving a system of differential equations subject to the flow conditions at the bounding upper and lower surfaces, as well as the continuity conditions at the interface of every pair of contiguous layers. Formidable computational efforts are required to solve the ensuing transcendental equations expressing the matching conditions at the interfaces, using this method. In this paper, the jump discontinuities in the flow parameters upon crossing from one layer to the other have been systematically built into a single partial differential equation governing the space–time variation of the excess pore pressure in the entire composite medium, by the use of the Heaviside distribution. Despite the presence of the discontinuities in the coefficients of the differential equation, a closed‐form solution in the sense of an infinite generalized Fourier series is obtained, in addition to which is the development of a Green's function for the differential problem. The eigenfunctions of the composite medium are the coordinate functions of the series, obtained computationally through the application of the extended equations of Galerkin. The analysis has been illustrated by solving the consolidation problem of a four‐layer composite, and the results obtained agree very well with the results obtained by previous researchers. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Curve‐Number/Green–Ampt mixed procedure for streamflow predictions in ungauged basins: Parameter sensitivity analysis     

S. Grimaldi  A. Petroselli  N. Romano 《水文研究》2013,27(8):1265-1275

A procedure combining the Soil Conservation Service‐Curve Number (SCS‐CN) method and the Green–Ampt (GA) infiltration equation was recently developed to overcome some of the drawbacks of the classic SCS‐CN approach when estimating the volume of surface runoff at a sub‐daily time resolution. The rationale of this mixed procedure, named Curve Number for Green–Ampt (CN4GA), is to use the GA infiltration model to distribute the total volume of the net hyetograph (rainfall excess) provided by the SCS‐CN method over time. The initial abstraction and the total volume of rainfall given by the SCS‐CN method are used to identify the ponding time and to quantify the hydraulic conductivity parameter of the GA equation. In this paper, a sensitivity analysis of the mixed CN4GA parameters is presented with the aim to identify conditions where the mixed procedure can be effectively used within the Prediction in Ungauged Basin perspective. The effects exerted by changes in selected input parameters on the outputs are evaluated using rectangular and triangular synthetic hyetographs as well as 100 maximum annual storms selected from synthetic rainfall time series. When applied to extreme precipitation events, which are characterized by predominant peaks of rainfall, the CN4GA appears to be rather insensitive to the input hydraulic parameters of the soil, which is an interesting feature of the CN4GA approach and makes it an ideal candidate for the rainfall excess estimation at sub‐daily temporal resolution at ungauged sites. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献