共查询到20条相似文献，搜索用时 46 毫秒
1.
The goal of the research was to demonstrate the impact of thin porous interfacial transition zones (ITZs) between aggregates and cement matrix on fluid flow in unsaturated concrete caused by hydraulic/capillary pressure. To demonstrate this impact, a novel coupled approach to simulate the two-phase (water and moist air) flow of hydraulically and capillary-driven fluid in unsaturated concrete was developed. By merging the discrete element method (DEM) with computational fluid dynamics (CFD) under isothermal settings, the process was numerically studied at the meso-scale in two-dimensional conditions. A flow network was used to describe fluid behaviour in a continuous domain between particles. Small concrete specimens of a simplified particle mesostructure were subjected to fully coupled hydro-mechanical simulation tests. A simple uniaxial compression test was used to calibrate the pure DEM represented by bonded spheres, while a permeability and sorptivity test for an assembly of spheres was used to calibrate the pure CFD. For simplified specimens of the pure cement matrix, cement matrix with aggregate, and cement matrix with aggregate and ITZ of a given thickness, DEM/CFD simulations were performed sequentially. The numerical results of permeability and sorptivity were directly compared to the data found in the literature. A satisfactory agreement was achieved. Porous ITZs in concrete were found to reduce sorption by slowing the capillary-driven fluid flow, and to speed the full saturation of pores when sufficiently high hydraulic water pressures were dominant. 相似文献
2.
R. Jeffrey Dunn 《Engineering Geology》1985,21(3-4):215-223
Advective fluid flow, which depends directly on hydraulic conductivity, controls contaminant migration through soils in many cases. The state of practice in laboratory measurements of hydraulic conductivity of soils was reviewed and found to be extremely variable, with the result that measurements on a given soil material often give variations of several orders of magnitude. An experimental program was completed using compacted clay soils to study several facets of laboratory testing practice, and a suitable test procedure for consistent results was developed. Recommendations on several important aspects of procedure are presented. 相似文献
3.
Over the last decade, waste disposal has become a particularly sensitive issue in Algeria. New legislation concerning landfill
liner design has been adopted. Traditional methods of landfill liner characterization involve soil sampling and chemical analysis,
which are costly, destructive and time-consuming. New techniques are currently being investigated that aim to provide nondestructive
liner characterisation. This paper details technical aspects associated with electrical conductivity measurements within landfill
liners and presents experimental work to show the direct application of electrical techniques to track ionic movement through
a sand bentonite liner under chemically induced flow. Samples of sand bentonite were mixed and compacted with NaCl electrolytes
at different concentrations. The electrical conductivities of compacted specimens were measured with a two-electrode cell.
The effects of frequency and electrolyte concentration on the conductivity measurement were explored. The relationship between
the soil electrical conductivity and the NaCl electrolyte concentration in interstitial pore fluid was determined. The conductivity
measurements were used to quantify the pore fluid concentration and effective diffusion coefficient of sand bentonite liners.
It is concluded here that the electrical conductivity of compacted specimens depends mainly on the salt concentration in the
pore fluid, and that this approach could therefore be used to track ionic movement through liners during diffusion. 相似文献
4.
Guillermo A. Narsilio Olivier Buzzi Stephen Fityus Tae Sup Yun David W. Smith 《Computers and Geotechnics》2009,36(7):1200-1206
The accurate estimation of hydraulic conductivity is important for many geotechnical engineering applications, as the presence of fluids affects all aspects of soil behaviour, including its strength. Darcy’s law is the key experimental (or phenomenological) equation employed to model ground water flow. Yet, this phenomenological equation can be linked to a more fundamental microscale model of flow through the pore spaces of the porous material. This paper provides an experimental verification of the relationships between Darcy’s law (macroscale) and the Navier–Stokes equations (microscale) for actual complex pore geometries of a granular material. The pore geometries are experimentally obtained through state-of-the-art X-ray computer assisted micro-tomography. From the numerical modelling of the microscale flow based on actual pore geometries, it is possible to quantify and visualize the development of pore-scale fluid preferential flow-paths through the porous material, and to assess the importance of pore connectivity in soil transport properties. 相似文献
5.
Sellafield in West Cumbria was a potential site for the location of the UK's first underground repository for radioactive, intermediate level waste (ILW). The repository was to lie around 650 m beneath the ground surface within rocks of the Borrowdale volcanic group (BVG), a thick suite of SW dipping, fractured, folded and metamorphosed Ordovician meta-andesites and ignimbrites. These are overlain by an onlapping sequence of Carboniferous and Permo-Triassic sediments. In situ borehole measurements showed that upward trending fluid pressure gradients exist in the area of the potential repository site, and that there are three distinct fluid types in the subsurface; fresh, saline and brine (at depth, to the west of the site). Simulations of fluid flow in the Sellafield region were undertaken with a 2D, steady-state, coupled fluid and heat flow simulation code (OILGEN). In both simplified and geologically complex models, topographically driven flow dominated the regional hydrogeology. Fluids trended persistently upwards through the potential repository site. The dense brine to the west of the site promoted upward deflection of topographically driven groundwaters. The inclusion in hydrogeological models of faults and variably saline sub-surface fluids was essential to the accurate reproduction of regional hydraulic head variations. Sensitivity analyses of geological variables showed that the rate of groundwater flow through the potential repository site was dependent upon the hydraulic conductivity of the BVG, and was unaffected by the hydraulic conductivity of other hydrostratigraphic units. Calibration of the model was achieved by matching simulated subsurface pressures to those measured in situ. Simulations performed with BVG hydraulic conductivity 100 times the base case median value provided the “best-fit” comparison between the calculated equivalent freshwater head and that measured in situ, regardless of the hydraulic conductivity of other hydrostratigraphic units. Transient mass transport simulations utilising the hydraulic conductivities of this “best fit” simulation showed that fluids passing through the potential repository site could reach the surface in 15 000 years. Simple safety case implications drawn from the results of the study showed that the measured BVG hydraulic conductivity must be less than 0.03 m year−1 to be simply declared safe. Recent BVG hydraulic conductivity measurements showed that the maximum BVG hydraulic conductivity is around 1000 times this safety limit. 相似文献
6.
An unsaturated hydraulic conductivity model for compacted GMZ01 bentonite with consideration of temperature 总被引：1，自引：0，他引：1
W. M. Ye M. Wan B. Chen Y. G. Chen Y. J. Cui J. Wang 《Environmental Earth Sciences》2014,71(4):1937-1944
As one of the most important properties of compacted bentonite used as buffer/backfill materials, hydraulic conductivity is influenced by various factors including temperature, microstructure and suction (or degree of saturation), etc. Based on the readily available results of both temperature-controlled water-retention tests and unsaturated infiltration tests under confined (constant volume) conditions, influences of temperature and microstructure variations on unsaturated hydraulic conductivity of the compacted Gaomiaozi (GMZ01) bentonite were analyzed. Then, a revised unsaturated hydraulic conductivity model considering temperature effects and microstructure changes was developed. With this proposed model, prediction and comparison were made on the unsaturated hydraulic conductivity of the compacted GMZ01 bentonite at 20 °C. Results show that water-retention capacity of compacted GMZ01 bentonite decreases as temperature increases and the degree of the temperature influence depends on suction. Under confined conditions, influence of hydration on microstructure of compacted GMZ01 bentonite depends on pore size. The proposed model can well describe the variations of unsaturated hydraulic conductivity with suction at different temperatures. However, further improvement of the proposed model is needed to account for the phenomenon of inter-aggregate pores clogging that occurred at the initial stage of hydration of compacted GMZ01 bentonite under confined conditions. 相似文献
7.
Gap-graded granular soils are used as construction materials worldwide, and their hydraulic conductivity depends on their relative content of coarse and fine grains, initial conditions, and particle shape. In this study, a series of constant head hydraulic conductivity tests were performed on gap-graded granular soils with different initial relative densities, fine contents, and particle shapes. The test results show that the hydraulic conductivity decreases with an increase in fine fraction and then remains approximately constant beyond the “transitional fine content.” The role of the structural effect on the hydraulic conductivity is different from that on the mechanical properties (such as stiffness and shear strength). This can be attributed to the degree of filling within inter-aggregate voids, disturbance of soil structure, and densified fine bridges between coarse aggregates. The equivalent void ratio concept was introduced into the Kozeny–Carman formula to capture the effect of fines (aggregates) on the “coarse-dominated” (“fine-dominated”) structure, and a simple model is proposed to capture the change of hydraulic conductivity of gap-granular soils. The model incorporates a structural variable to capture the effect of fines on “coarse-dominated” structure and coarse aggregates on “fine-dominated” structure. The performance of the model was verified with experimental data from this study and previously reported data compiled from the literature. The results reveal that the proposed model is simple yet effective at capturing the hydraulic conductivity of gap-graded granular soils with a wide range of fine contents, initial conditions, and particle shapes.
相似文献8.
Internal erosion by suffusion can change dramatically the constitutive behavior of granular materials by modifying the fabric of granular materials. In this study, the effect of an internal fluid flow on granular materials is investigated at the material point scale using the numerical coupling between a discrete element method (DEM) and a pore-scale finite volume (PFV) coupling scheme. The influence of the stress state and the hydraulic loading (direction and intensity) on the occurrence of grain transport in dense widely graded granular samples is thus investigated and interpreted in terms of micromechanics. In particular, it is shown that grain transport is increased when the macroscopic flow direction is aligned with the privileged force chain orientation. The stress-induced microstructure modifications are shown to influence the transport distances by controlling the number of rattlers. 相似文献
9.
Numerical simulation of the quicksand phenomenon by a 3D coupled Discrete Element ‐ Lattice Boltzmann hydromechanical model 
下载免费PDF全文
下载免费PDF全文 This paper deals with the numerical simulation of the quicksand phenomenon using a coupled Discrete Elements – Lattice Boltzmann hydromechanical model. After the presentation of the developed numerical model, simulations of ascending fluid flow through granular deposits are performed. The simulations show that the quicksand actually triggers for a hydraulic gradient very close to the critical hydraulic gradient calculated from the global analysis of classical soil mechanics, that is, when the resultant of the applied external pressure balances submerged weight of the deposit. Moreover, they point out that the quicksand phenomenon does not occur only for hydraulic gradients above the critical hydraulic gradient, but also in some cases with slightly lower gradients. In such cases, a more permeable zone is first gradually built at the bottom of the deposit through a grain rearrangement, which increases the hydraulic gradient in the upper zones and triggers the phenomenon. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
10.
三轴试验的结果表明,由于应用于水利工程的沥青混凝土自身的孔隙率小,侧向作用效应十分明显,体积变形较小,且主要表现为剪胀特性,目前已有的非线性和弹塑性模型难以正确地反映复杂应力状态下的这种力学行为。考虑沥青混凝土不遵守关联流动法则和Drucker公设,其弹性系数伴随塑性变形过程而改变,建议了一个在弹塑性耦合基础上的双屈服面弹塑性模型,与室内三轴试验结果吻合良好。 相似文献
11.
考虑颗粒破碎影响的粗粒土本构模型 总被引：2，自引：1，他引：1
颗粒破碎直接改变了粗粒土本身结构,对粗粒土的剪胀、内摩擦角、峰值强度、渗透系数都会产生影响。为了能够准确地描述粗粒土的应力-应变关系,特别是高应力条件下出现显著颗粒破碎时的应力-应变关系,亟待建立考虑颗粒破碎的粗粒土本构模型。根据三轴试验数据,建立考虑颗粒破碎耗能的应力-应变关系,采用相关联流动法则导出考虑颗粒破碎的粗粒土本构模型。所建立的本构模型考虑了颗粒破碎对粗粒土剪胀、内摩擦角的影响。通过变异粒子群优化算法确定模型参数,拟合试验曲线。模型计算结果与试验曲线拟合较好,能够很好地描述粗粒土在不同围压下的体积剪胀、剪缩和应力硬化、软化现象。 相似文献
12.
考虑颗粒破碎的粗粒土剪胀性统一本构模型 总被引：2，自引：0，他引：2
粗粒土作为无黏性散粒状材料具有状态依赖特性,土体的剪切特性受密度和应力水平影响。易破碎是粗粒土的另一个特点,颗粒破碎影响粗粒土的剪胀、内摩擦角、峰值强度和渗透系数。为了能够准确地描述粗粒土的应力-应变关系,采用初始状态参量描述粗粒土的内部状态,根据三轴试验数据建立考虑颗粒破碎耗能的应力-应变关系,采用相关联流动法则推导考虑颗粒破碎的粗粒土剪胀性“统一本构模型”,并建立初始状态参量与模型参数之间的关系。所建立的统一本构模型既考虑了颗粒破碎对剪胀、内摩擦角的影响,又考虑了剪切特性对土体初始状态的依赖。采用变异粒子群算法拟合试验曲线,确定模型参数。模型计算结果能够很好地拟合试验曲线。采用同一组参数对假定的初始状态进行模拟计算,计算结果表明,模型能够模拟不同初始密度和应力水平下粗粒土变形的一般规律。 相似文献
13.
离散裂隙渗流方法与裂隙化渗透介质建模 总被引：4，自引：1，他引：4
流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。 相似文献
14.
In this article, we present a numerical method to deal with fluid–solid interactions and simulate particle–fluid systems as encountered in soils. This method is based on a coupling between two methods, now widely used in mechanics of granular media and fluid dynamics respectively: the discrete element (DE) method and the lattice Boltzmann (LB) method. The DE method is employed to model interactions between particles, whereas the LB method is used to describe an interstitial Newtonian fluid flow. The coupling presented here is a full one in the sense that particle motions act on fluid flow and reciprocally. This article presents in details each of the two methods and the principle of the coupling scheme. Determination of hydrodynamic forces and torques is also detailed, and the treatment of boundaries is explained. The coupled method is finally illustrated on a simple example of piping erosion, which puts in evidence that the combined LB–DE scheme constitutes a promising tool to study coupled problems in geomechanics. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
15.
Granular materials can behave as harmless sand dunes or as devastating landslides. A granular avalanche marks the transition between these distinct solid-like and fluid-like states. The solid-like state is typically described using plasticity models from critical state theory. In the fluid regime, granular flow is commonly captured using a visco-plastic model. However, due to our limited understanding of the mechanism governing the solid–fluid-like transition, characterizing the material behavior throughout the life cycle of an avalanche remains an open challenge. Here, we employ laboratory experiments of transient avalanches spontaneously generated by a rotating drum. We report measurements of dilatancy and grain kinematics before, during, and after each avalanche. Those measurements are directly incorporated into a rate-dependent plasticity model that quantitatively predicts the granular flow measured in experiments. Furthermore, we find that dilatancy in the solid-like state controls the triggering of granular avalanches and therefore plays a key role in the solid–fluid-like transition. With the proposed approach, we demonstrate that the life cycle of a laboratory avalanche, from triggering to run out, can be fully explained. Our results represent an important step toward a unified understanding of the physical phenomena associated with transitional behavior in granular media. 相似文献
16.
Thermochronometry based on radiogenic noble gases is critically dependent upon accurate knowledge of the kinetics of diffusion. With few exceptions, complex natural crystals are represented by ideal geometries such as infinite sheets, infinite cylinders, or spheres, and diffusivity is assumed to be isotropic. However, the physical boundaries of crystals generally do not conform to ideal geometries and diffusion within some crystals is known to be anisotropic. Our failure to incorporate such complexities into diffusive models leads to inaccuracies in both thermal histories and diffusion parameters calculated from fractional release data. To address these shortcomings we developed a code based on the lattice Boltzmann (LB) method to model diffusion from complex 3D geometries having isotropic, temperature-independent anisotropic, and temperature-dependent anisotropic diffusivity. In this paper we outline the theoretical basis for the LB code and highlight several advantages of this model relative to more traditional finite difference approaches. The LB code, along with existing analytical solutions for diffusion from simple geometries, is used to investigate the affect of intrinsic crystallographic features (e.g., crystal topology and diffusion anisotropy) on calculated diffusion parameters and a novel method for approximating thermal histories from crystals with complex topologies and diffusive anisotropy is presented. 相似文献
17.
基于人工交叉裂隙模型,通过室内透水试验,利用电荷耦合元件（CCD）照相机可视化技术,对流体在裂隙交叉点内的非线性流动特性进行研究。建立两种离散裂隙网络（DFN）模型,考虑两种边界条件,改变模型进口和出口之间的压力,直接求解Navier-Stokes（简称N-S）方程,对DFN的非线性渗流特性进行研究。结果表明,室内试验可以观测到与出口3相连的裂隙单元内发生了明显的非线性流动,且通过模型的流量Q和模型两端的压力P具有非线性关系。数值计算结果也表明,在水力梯度J较大时（比如J > 0.1）,通过DFN的Q和P具有非线性关系,而当J较小时（比如J < 10-4）,Q与P线性相关;根据文中的算例,建议利用局部立方定律求解DFN内每条裂隙的渗流特性的临界条件为J ≤10-4;裂隙表面粗糙会造成通过DFN渗流量的降低,但对相对流量误差的影响可忽略不计。 相似文献
18.
Liu Zhang-Rong Ye Wei-Min Cui Yu-Jun Zhu He-Hua Wang Qiong Chen Yong-Gui 《Acta Geotechnica》2021,16(10):3145-3160
Bentonite pellets are recognized as good buffer/backfill materials for sealing technological voids in high-level radioactive waste (HLW) repository. Compared to that of a traditional compacted bentonite block, one of the most important particularities of this material is the initially discrete pellets and the inevitable heterogeneous porosity formed, leading to a distinctive water retention behaviour. In this paper, water retention and mercury intrusion porosimetry (MIP) tests were conducted on pellet mixture (constant volume), single pellet (free swelling) and compacted block (constant volume) of GMZ bentonite, water retention properties and pore structure evolutions of the specimens were comparatively investigated. Results show that the water retention properties of the three specimens are almost similar to each other in the high suction range (>?10 MPa), while the water retention capacity of pellet mixture is lower than those of the compacted block and single pellet in the low suction range (<?10 MPa). Based on the capillary water retention theory (the Young–Laplace equation), a new concept ‘saturated void ratio’ that was positively related to water content and dependent on pore size distribution of the specimen was defined. Then, according to the product of saturated void ratio and water density in saturated void, differences of water retention properties for the three specimens at low suctions were explained. Meanwhile, MIP tests indicate that as suction decreases, the micro- and macrovoid ratios of pellet mixture and compacted block decrease as the mesovoid ratio increases, while all the void ratios of single pellets increase. This could be explained that upon wetting, water is successively adsorbed into the inter-layer, inter-particle and inter-pellet voids, leading to the subdivision of particles and swelling of aggregates and pellets. Under constant volume condition, aggregates and pellets tend to swell and fill into the inter-aggregates or inter-pellets voids. While under free swelling condition, the particles and aggregates in a single pellet tend to swell outward rather than squeezing into the inter-aggregate voids, leading to the expansion of the pores and even formation of cracks. Results including the effects of initial conditions (initial dry density and fabric) and constraint conditions (constant volume or free swelling) on the water retention capacity and pore structure evolution reached in this work are of great importance in designing of engineering barrier systems for the HLW repository.
相似文献19.
A program for the simulation of two‐dimensional (2‐D) fluid flow at the microstructural level of a saturated anisotropic granular medium is presented. The program provides a numerical solution to the complete set of Navier–Stokes equations without a priori assumptions on the viscous or convection components. This is especially suited for the simulation of the flow of fluids with different density and viscosity values and for a wide range of granular material porosity. The analytical solution for fluid flow in a simple microstructure of porous medium is used to verify the computer program. Subsequently, the flow field is computed within microscopic images of granular material that differ in porosity, particle size and particle shape. The computed flow fields are shown to follow certain paths depending on air void size and connectivity. The permeability tensor coefficients are derived from the flow fields, and their values are shown to compare well with laboratory experimental data on glass beads, Ottawa sand and silica sands. The directional distribution of permeability is expressed in a functional form and its anisotropy is quantified. Permeability anisotropy is found to be more pronounced in the silica sand medium that consists of elongated particles. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
20.
This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献