METHODOLOGY FOR DRENABILITY STUDIES OF STACKED MINING COARSE TAILINGS
|
|
- Roy Bishop
- 6 years ago
- Views:
Transcription
1 MEHODOLOGY FOR DRENABILIY SUDIES OF SACKED MINING COARSE AILINGS J. C. Machado Júnior Departamento de Engenharia Civil, Universidade Federal de Ouro Preto W. L. Oliveira Filho Departamento de Engenharia Civil, Universidade Federal de Ouro Preto Abstract: his paper presents a methodology for drenability studies of coarse mining tailings. hese studies have a particular impact in the wall stability and construction methods of stacked mining coarse tailings structures (Gomes et al. 1999a). In that context, flow processes on unsaturated soil are the major complicated factor. For that purpose, a short review of this flow problem is presented, including theory and flow material functions needed in the analysis and their determination. Some common scenarios of flow problems related to tailings piling up are discussed and analyzed using a commercial code, SEEP/W, as a numerical tool. his paper also aims to make a contribution to the diffusion of unsaturated soil mechanics out of the academic circles. 1. INRODUCION Nowadays the mining tailings disposal management is often an activity of primary importance for the success of a mining operation. he economics and also the environmental issues associated with that activity have made the mining industries to optimize the disposal areas and to develop new processes to discard tailings. o pile up coarse tailings is one of these new techniques for the coarse fraction of the waste generated by the iron ore processing plants (Gomes et al. 1999a). In this regard, one of the key issues of the stability analysis of these structures is the question of drenability of the tailings. he project will be more economical and safe, depending on how precise a model can be made of the flow processes taking place in the tailings. However, a factor that complicates the flow analysis is that the process of drainage needs to deal with unsaturated porous media, and a complete new field of study for geotechnical engineers has to be faced. In this paper, an effort to popularize unsaturated flow analysis out of the academic circles is made by addressing a methodology to study a real problem on mining tailings disposal. 2. FUNDAMENALS OF UNSAURAED FLOW brief review of unsaturated flow theory is worth to be presented here since there is a lack of general understanding in that subject within our community of geotechnical engineers HE GOVERNING EQUAION he governing equation of the flow processes in soils is well known from textbooks of basic soil mechanics (Craig 1996, Lambe and Whitman 1979, Ortigão 1996). his equation can be expressed in a vectorial form like, dm ( ρ ) dx dydz w wv = (1), dt where, ρ w is the water density, v is the darcian velocity vector, is a differential operator so that, = [ / x / y / z ], M W is the water mass, and t is the time. 233
2 Equation (1) relates the mass balance in a three dimensional set, x, y and z, with the storage term on the right hand side. he latter can be expressed in terms of soil porosity, n, and degree of saturation, S, so that, M w = ρw S n dxdydz (2) and (1) becomes ( ρ Sn W ) ( ρ W v ) = (3). Depending on the way that the properties ρ w, S, and n are considered, different soil flow processes are described by (3). he usual simplification for most of soil engineers is to consider ρ w, S and n as constants, and therefore their derivatives are equal to zero. he right hand side then is equal to zero and this is the case of saturated, steady state flow analysis. ( v ) = 0 (4). In another case, when hypotheses of isothermal conditions and saturated porous media are introduced, the derivatives on ρ w and S are zero, and (3) becomes n ( v ) = (5), which is the general equation for (large strain) consolidation analysis. he less known form of (3), at least for geotechnical engineers, appears when the derivatives of ρ w and n are zero, and (3) becomes ( v ) = n S (6). Considering the Darcy s law, (6) becomes ( K h) = n S (7), where K is the hydraulic conductivity tensor and h is the total head. he governing equation in this form describes the transient analysis of unsaturated flow of incompressible porous media and isothermal conditions. 2.2 Richard's equation Further developments of the governing equation for unsaturated flow are possible. he concept of volumetric water content, θ, is very often needed in such cases. he relation between water volume in the voids, V w, and total soil volume, V, defines this quantity. hus, Vw θ = = V n S (8). Also, the total hydraulic head, h, can be divided in two other terms, the elevation head, z, and pressure head, ψ. hus, h = z + ψ (9). Equation (7) then becomes, θ ( Ke + K ψ) = (10), where e is a vector with the component z unitary, so that e =[ ]. he last development of the flow equation is to take into account that K and θ are functions of pressure head in unsaturated flow. It means that, K = K (ψ) (11) and θ = θ(ψ) (12), which implies, θ = c( ψ) ψ (13), where c(ψ)= θ/ ψ is called specific capacity retention. By considering these relationships, the primary dependent variable in the flow equation becomes ψ and the governing equation assumes the form 234
3 ( K( ψ) e + K( ψ) ψ) = C( ψ) ψ (14). his ψ-based equation of flow for transient flow through an unsaturated media is called Richard's equation (Freeze and Cherry 1979). It s a nonlinear second order partial differential equation. he nonlinearity on (14) is because the derivative coefficients are not constants but dependent on the primary variable, ψ. he solution of (14) must satisfy the problem boundary conditions. hese can be a prescribed pressure head boundary condition (Dirichlet condition) so that, ψ (x,t) = ψ on Γ 1 (15), or a prescribed normal flux boundary condition (Neuman condition) so that, [ K( ψ) e + K( ψ) ψ] v n = on Γ 2 (16), where Γ=Γ 1 +Γ 2, is the boundary of the volume, V; ψ and v are prescribed values, n is the unit normal vector on Γ 2 and x is a position vector. he solution of (14) must also satisfy the initial condition, ψ( x,0) = ψ 0 (17), where ψ 0 is a known value. o solve (14), flow material functions (11) and (12) are required. he determination of these functions are discussed in section 3. Closed form solutions of (14) exist only in very few cases what makes the use a numerical scheme a necessity. 2.3 Applications o illustrate the use of the flow equation, the analyses of two very common field situations involving flow in unsaturated media are presented Infiltration in a soil column o study this situation, an unsaturated soil layer is thought to receive on its top boundary a constant rain precipitation, R, with dimensions [L/]. Because the problem can be considered one-dimensional, the governing equation, (14), can be written as ψ ψ ψ ψ z k ( )( 1 + z c z ) = ( ) (18). From Darcy's law, v=ki, where i = h/ z and h/ z = 1+ ψ/ z, hence at top boundary R = + 1 ψ z k ( ψ) (19) or ψ R = 1 (20). z k( ψ) At the bottom boundary, a Dirichlet prescription is usually assumed, like ψ = 0 (21). (18), (20) and (21) constitute the boundary value problem for infiltration in a soil column Drainage in a soil column his situation can also be considered onedimensional and therefore the governing equation looks the same as in the previous case. ψ ψ ψ ψ z k ( )( 1 + z c z ) = ( ) (22). he top boundary is now impervious which implies R = 0 or ψ/ z = -1 (23). At the bottom, it is usually assumed again, ψ = 0 (24). 235
4 (22), (23) and (24) defines the boundary value problem for drainage in a soil column. 3. DEERMINAION OF HE FLOW MAERIAL FUNCIONS As it was mentioned, the determination of k= k(ψ) and θ = θ(ψ) is required for the solution of the governing equation. hose functions are called characteristic curves. he first is the hydraulic conductivity function and the second is the soil-water characteristic curve (SWCC). hese curves often have sigmoid shapes and present hysteresis for infiltration and drainage cycles. he determination of characteristic curves can be achieved by specific laboratory or field tests or even by correlation (especially the hydraulic conductivity function). Benson and Gribb (1997) present a comprehensive review of the procedures available. In the laboratory the water retention curve can be obtained through the pressure plate test according to the ASM D he setup of the test is schematically shown in Fig. 1. A chamber houses the soil specimen that is placed on a porous disc of high air entry value. he disc is in contact with a water reservoir at atmospheric pressure. he difference between the pore-air, u a, and pore-water, u w, pressures is called matric suction and it s created by pressurizing the chamber with compressed air. he amount of water that leaves the specimen with the increasing pressure is monitored, allowing the determination of a relationship between the degree of saturation and suction. Both drainage and infiltration cycles are run in the test. It is possible to measure the permeability function directly but intrinsic test difficulties and lack of accuracy very often prevent its large use and lead the practitioner to use indirect methods. he usual procedure is to estimate the hydraulic conductivity function from the SWCC. he greater reliability and short period of time required for the SWCC determination favor this approach. he correlation also has its strength in the intimacy of the two relationships. he SWCC describes the amount of water present in a soil under various matric suctions. his curve essentially indicates the space available for the water to flow through the soil at various matric suctions since water can only flow through the waterfilled pores. herefore the shape of the curve can be used to estimate the permeability function. he function decreases monotonically from the saturated permeability, k sat, at zero matric suction and assumes lower values as the matric suction increases. Figure 1 Schematic diagram of a pressure plate extractor (Rahardjo and Leong, 1997) here are many empirical equations that relate the hydraulic conductivity to the volumetric water content or degree of saturation of the soil (Gardner, 1958, Brooks and Corey, 1964; van Genuchten, 1980). he van Genuchten equation is given below: β 1/ 2 β 1 k r = θ e 1 1 θ e β where, k r relative hydraulic conductivity, β - constant, and hydraulic conductivity is given by (21) k = k sat k r (22). he β coefficient of (21) is a fitting parameter obtained from the experimental data of SWCC and the theoretical model given below also due to van Genuchten, θ e = 1 [ 1+ ( h) ] α β β 1 1 (23) 236
5 and θ e θ θr = θ θ s r (24) where, θ - volumetric water content, θ e equivalent volumetric water content, θ r residual volumetric water content, θ s saturated volumetric water content, h = ( u a - u w ) matric suction, α, β - constants. Figure 2 Soil water characteristic curve (SEEP/W user s guide) 4. NUMERIC ANALYSIS In the following subitems a complete example of a numeric simulation is given. he analyses were performed using a commercial computer program called SEEP/W. his is a finite element software used to model seepage problems, developed by GEO-SLOPE International Ltd. 4.1 Flow material functions he characteristic curves used in the analyses were obtained from the piece of literature (SEEP/W user s guide) and correspond to a fine sand, k sat = 4.3 x 10-6, n=0.22, d 10 = mm. Figure 2 and 3 illustrate the flow material functions used in the analyses. 4.2 Scenarios for the numeric analyses Several field situations can be analyzed within a tailings management strategy. In this paper three scenarios were envisioned for the study of tailings drenability. hese case studies are very often found in practice and the flow processes happen to be gravitational and onedimensional (Collins and Znidarcic, 1997). he boundary conditions for these situations are presented in terms of total head, H, total nodal flow, Q, and flow per unit length, q. In the first scenario, a one-meter thick tailings column, initially saturated, is underlain by a drainage layer, and is left to drain freely. his situation is shown in Figure 4(a). Figure 3 Hydraulic condutivity function (SEEP/W user s guide) he second scenario is a follow-up to the previous case with the initial conditions corresponding to the end of that analysis. In the sequence, a constant flow rate is imposed in the upper boundary, simulating rain precipitation over the drained layer. Figure 5 (a) is a schematic drawing of that case study. he last scenario allows the evaluation of the flow pattern considering multiple layers in different drainage stages. In the actual case, the behavior of a newly disposed layer (saturated) and its iteration with the underlain layer at certain drainage stage are analyzed. his scenario is sketched in Figure 6 (a). 4.3 Analysis Results he results of the numeric simulation of a draining soil column are shown in Figure 4 (b) in terms of pressure head x depth x time. he instant t = 0 reflects the initial conditions of the problem (a hydrostatic water column). It can be seen that for t > 0 the imposed boundary conditions are Q = 0 at top boundary and H = 0 237
6 at the bottom boundary. Hydrostatic conditions firstly develop at top portion of the column and progress downwards with time. At the end steady state conditions are reached with a hydrostatic suction profile throughout the entire column. he results of the infiltration case are shown in Figure 5 (b). At t = 0, the initial conditions corresponding to a hydrostatic suction profile. For t > 0, the boundary conditions are q = 4.3 x 10-6 m/s at the top boundary, and H = 0 at bottom boundary. At the beginning of the simulation, at the top portion, it can be seen that high gradients are necessary to overcome the low permeability values associated with high matric suctions and to allow seepage to take place. As the degree of saturation increases, the entrance gradients are smaller, tending towards a gravitational value of 1. Because q = k sat, the steady state solution gives a profile with constant suction from top to bottom. he corresponding degree of saturation profile would be similar with a constant S = 1.0 from top to bottom. he last case results are shown in Figure 6 (b). he initial conditions reflect two different situations. From 0 to 1m a suction profile corresponding to a layer at certain stage in a drainage process is shown. From 1 to 2 m, a hydrostatic positive porewater pressure profile is seen. his one corresponds to a newly disposed saturated layer. In the initial lapses of time, t > 0, a drainage process of the top layer is taking place while the bottom layer is experiencing a saturation phase. At a certain point in those flow processes, a typical steady state infiltration profile can be seen but only instantaneously, because for later lapses of time, a soil column drainage pattern governs. In the last stage, the two layers behave like one, in a very similar way to the results of the first case study. he boundary conditions in this case are Q = 0 at the top boundary and H=0 at the bottom. Figure 4 - (a) Finite element mesh and boundary conditions to scenario #1; (b) SEEP/W simulation results to scenario #1. 238
7 Figure 5 - (a) Finite element mesh and boundary conditions to scenario #2; (b) SEEP/W simulation results to scenario #2. Figure 6 - (a) Finite element mesh and boundary conditions to scenario #3; (b) SEEP/W simulation results to scenario #3. 239
8 4.4 Discussion Questions about flow rate pressure head and time required for water to flow can now be easily and rationally addressed. Some variations of the case studies described before can be envisioned. Even more sophisticated situations can also be considered such as twodimensional flow and different materials, without introducing any major change in the methodology here described. his kind of analysis can be used in the design of stacked mining coarse tailings structures in areas such as wall stability studies, monitoring programs and construction methods (Znidarcic, 1999). hese aspects will be object of future work. 5. CONCLUSION From the simple cases presented in this paper it is apparent that the methodology has great capability to expand the horizons to study flow processes in unsaturated soils. he requirements for this development include: to perform new but not sophisticated tests, to use a competent and reliable computer code such as SEEP/W and others, and most of all to get acquainted with unsaturated flow processes. An effort to bring about this issues has been made in this work. 6. REFERENCES Benson, C. H. and Gribb, M. M. (1997).Measuring Unsaturated Hydraulic Conductivity in Laboratory and Field, Unsaturated Soil Engineering Practice. Geothecnical Special Publication (68) ASCE, Brooks, R. H. and Corey, A.. (1964). Hydraulic Properties of Porous Medium. Hydrology Paper No.3, Civil Engineering dept., Colorado State University, Fort Collings, Colorado. Craig, R. F. (1997) Soil Mechanics, Chapman & Hall, sixth edition. Collins, B. D., and Znidarcic, D. (1997), riggering Mechanics of Rainfall Induced Debris Flows, Proceedings of 2 nd Pan- American Symposium on Landslides and 2 nd Brazilian Conference on Slope Stability, II PSL/COBRE, Vol.1, Freeze, R.A. & Cherry,J.A. (1979). Groundwater, Prentice-Hall, Inc., Englewood Cliffs, New Jersey. Gardner, W. R. (1958). Some Steady State Solutions of the Unsaturated moisture Flow Equation with Aplications to Eva poration from Water able. Soil Science, 85, Geoslope International Ltd., User Guide SEEP/W for finite element seepage analysis. Version 3 Gomes, R. C., Araújo,L. G., Oliveira, W. L.F., Ribeiro, S. G. S., Nogueira, C. L., (1999a). Concepção e Projeto Básico da Disposição em Pilhas de Rejeito de Minério de Ferro em Cava Exaurida de Mineração. VI Congresso Brasileiro de Geotecnia Ambiental, REGEO 99. Lambe,. W., and Whitman, R. V.,(1979), Soil Mechanics, SI Version, John Wiley & Sons Ortigão, J. A. R., (1995), Introdução àmecânica dos Solos dos Estados Críticos, Livros écnicos e Científicos S. A.. Rahardjo, H. and Leong, E. C. (1997).Soilwater Characteristic Curves and Flux Boundary Probems, Unsaturated Soil Engineering Practice. Geothecnical Special Publication (68) ASCE, Van Genuchten, M.. (1980). A Closed Form Equation for Predicting the Hydraulic conductivity of Unsaturated Soils. Soil Science Society of America Journal, 44, Znidarcic, D. (1998). Report on the Review of Germano Exhausted Open Pit Project 240
In all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.
Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate
More informationInfluence of Soil Characteristic and Rainfall Intensity on Matric Suction of Unsaturated Weathered Soil Slope
Journal of the Korean Society of Civil Engineers Kim, Yong Min*, Lee, Kwang Woo**, Kim, Jung Hwan*** Influence of Soil Characteristic and Rainfall Intensity on Matric Suction of Unsaturated Weathered Soil
More information1. Water in Soils: Infiltration and Redistribution
Contents 1 Water in Soils: Infiltration and Redistribution 1 1a Material Properties of Soil..................... 2 1b Soil Water Flow........................... 4 i Incorporating K - θ and ψ - θ Relations
More informationTHEORY. Water flow. Air flow
THEORY Water flow Air flow Does Suction Gradient Cause Flow? Coarse stone Fine ceramic Suction gradient to the right No suction gradient but still flow Does water content gradient cause Flow? Suction gradient
More informationUnsaturated Flow (brief lecture)
Physical Hydrogeology Unsaturated Flow (brief lecture) Why study the unsaturated zone? Evapotranspiration Infiltration Toxic Waste Leak Irrigation UNSATURATAED ZONE Aquifer Important to: Agriculture (most
More informationComparison of Averaging Methods for Interface Conductivities in One-dimensional Unsaturated Flow in Layered Soils
Comparison of Averaging Methods for Interface Conductivities in One-dimensional Unsaturated Flow in Layered Soils Ruowen Liu, Bruno Welfert and Sandra Houston School of Mathematical & Statistical Sciences,
More informationChapter 7 Permeability and Seepage
Permeability and Seepage - N. Sivakugan (2005) 1 7.1 INTRODUCTION Chapter 7 Permeability and Seepage Permeability, as the name implies (ability to permeate), is a measure of how easily a fluid can flow
More informationGEO-SLOPE International Ltd, Calgary, Alberta, Canada Wick Drain
1 Introduction Wick Drain This example is about modeling the behavior of a wick drain. The primary purpose here is to illustrate how interface elements can conveniently be used to include the effects of
More informationPredicting the soil-water characteristics of mine soils
Predicting the soil-water characteristics of mine soils D.A. Swanson, G. Savci & G. Danziger Savci Environmental Technologies, Golden, Colorado, USA R.N. Mohr & T. Weiskopf Phelps Dodge Mining Company,
More informationSTABILITY OF RESIDUAL SOIL SLOPES BASED ON SPATIAL DISTRIBUTION OF SOIL PROPERTIES. Harianto Rahardjo*, Alfrendo Satyanaga
STABILITY OF RESIDUAL SOIL SLOPES BASED ON SPATIAL DISTRIBUTION OF SOIL PROPERTIES Harianto Rahardjo*, Alfrendo Satyanaga * Professor, School of Civil and Environmental Engineering, Nanyang Technological
More informationSoil Mechanics for Unsaturated Soils
Soil Mechanics for Unsaturated Soils Delwyn G. Fredlund University of Saskatchewan Saskatoon, Sask., Canada and Harianto Rahardjo Nanyang Technological University Singapore John Wiley & Sons Notes by:
More informationDarcy s Law, Richards Equation, and Green-Ampt Equation
Darcy s Law, Richards Equation, and Green-Ampt Equation 1. Darcy s Law Fluid potential: in classic hydraulics, the fluid potential M is stated in terms of Bernoulli Equation (1.1) P, pressure, [F L!2 ]
More informationψ ae is equal to the height of the capillary rise in the soil. Ranges from about 10mm for gravel to 1.5m for silt to several meters for clay.
Contents 1 Infiltration 1 1a Hydrologic soil horizons...................... 1 1b Infiltration Process......................... 2 1c Measurement............................ 2 1d Richard s Equation.........................
More informationTeaching Unsaturated Soil Mechanics as Part of the Undergraduate Civil Engineering Curriculum
Teaching Unsaturated Soil Mechanics as Part of the Undergraduate Civil Engineering Curriculum Delwyn G. Fredlund, Visiting Professor Kobe University, Kobe, Japan Sapporo, Hokkaido, Japan February 15, 2005
More informationEXAMPLE PROBLEMS. 1. Example 1 - Column Infiltration
EXAMPLE PROBLEMS The module UNSATCHEM is developed from the variably saturated solute transport model HYDRUS-1D [Šimůnek et al., 1997], and thus the water flow and solute transport parts of the model have
More informationTransient Analysis on Infiltration and Stability for Unsaturated Soils in Busan Landslide Area
Transient Analysis on Infiltration and Stability for Unsaturated Soils in Busan Landslide Area FAUZI ACHMAD ZAKY, SEBOONG OH Department of Civil Engineering Yeungnam University 8 Daehak-Ro, Gyeongsan,
More informationSimulation of Unsaturated Flow Using Richards Equation
Simulation of Unsaturated Flow Using Richards Equation Rowan Cockett Department of Earth and Ocean Science University of British Columbia rcockett@eos.ubc.ca Abstract Groundwater flow in the unsaturated
More informationSoil Water Atmosphere Plant (SWAP) Model: I. INTRODUCTION AND THEORETICAL BACKGROUND
Soil Water Atmosphere Plant (SWAP) Model: I. INTRODUCTION AND THEORETICAL BACKGROUND Reinder A.Feddes Jos van Dam Joop Kroes Angel Utset, Main processes Rain fall / irrigation Transpiration Soil evaporation
More informationHYDROLOGICAL CONDITION LEADING TO LANDSLIDE INITIATION
HYDROLOGICAL CONDITION LEADING TO LANDSLIDE INITIATION Adrin Tohari Indonesian Institute of Sciences, Research Center for Geotechnology, Jln Sangkuriang, Bandung 413. Telp: 364. Email: adrin@cbn.net.id
More informationDistribution of pore water pressure in an earthen dam considering unsaturated-saturated seepage analysis
E3S Web of Conferences 9, 194 (16) DOI: 1.11/ e3sconf/169194 E-UNSAT 16 Distribution of pore water in an earthen dam considering unsaturated-saturated seepage analysis 1a Kumar Venkatesh, Siva Ram Karumanchi
More informationThe Use of COMSOL for Integrated Hydrological Modeling
The Use of for Integrated Hydrological Modeling Ting Fong May Chui *, David L. Freyerg Department of Civil and Environmental Engineering, Stanford University *Corresponding author: 380 Panama Mall, Terman
More information16 Rainfall on a Slope
Rainfall on a Slope 16-1 16 Rainfall on a Slope 16.1 Problem Statement In this example, the stability of a generic slope is analyzed for two successive rainfall events of increasing intensity and decreasing
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationTikrit University. College of Engineering Civil engineering Department CONSOILDATION. Soil Mechanics. 3 rd Class Lecture notes Up Copyrights 2016
Tikrit University CONSOILDATION College of Engineering Civil engineering Department Soil Mechanics 3 rd Class Lecture notes Up Copyrights 2016 Stresses at a point in a soil mass are divided into two main
More informationProcedia Earth and Planetary Science 9 ( 2014 ) The Third Italian Workshop on Landslides
Available online at www.sciencedirect.com ScienceDirect Procedia Earth and Planetary Science 9 ( 2014 ) 214 221 The Third Italian Workshop on Landslides Prediction of suction evolution of silty pyroclastic
More informationFilling Pond Head vs Volume Functions
Introduction Filling Pond Head vs Volume Functions The objective of this illustration is show how to model the filling of a pond where the water is seeping into the pond from the soil. The head in the
More informationCyclic Triaxial Behavior of an Unsaturated Silty Soil Subjected to Suction Changes
6 th International Conference on Earthquake Geotechnical Engineering 1-4 November 215 Christchurch, New Zealand Cyclic Triaxial Behavior of an Unsaturated Silty Soil Subjected to Suction Changes T. Nishimura
More informationAdvanced Hydrology Prof. Dr. Ashu Jain Department of Civil Engineering Indian Institute of Technology, Kanpur. Lecture 6
Advanced Hydrology Prof. Dr. Ashu Jain Department of Civil Engineering Indian Institute of Technology, Kanpur Lecture 6 Good morning and welcome to the next lecture of this video course on Advanced Hydrology.
More informationModelling of pumping from heterogeneous unsaturated-saturated porous media M. Mavroulidou & R.I. Woods
Modelling of pumping from heterogeneous unsaturated-saturated porous media M. Mavroulidou & R.I. Woods Email: M.Mavroulidou@surrey.ac.uk; R. Woods@surrey.ac.uk Abstract Practising civil engineers often
More informationPractical methodology for inclusion of uplift and pore pressures in analysis of concrete dams
Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams Michael McKay 1 and Francisco Lopez 2 1 Dams Engineer, GHD Pty 2 Principal Dams/Structural Engineer, GHD Pty
More informationModule 2 Lecture 9 Permeability and Seepage -5 Topics
Module 2 Lecture 9 Permeability and Seepage -5 Topics 1.2.7 Numerical Analysis of Seepage 1.2.8 Seepage Force per Unit Volume of Soil Mass 1.2.9 Safety of Hydraulic Structures against Piping 1.2.10 Calculation
More informationAir Flow Modeling. An Engineering Methodology. February 2012 Edition. GEO-SLOPE International Ltd.
Air Flow Modeling with AIR/W An Engineering Methodology February 2012 Edition GEO-SLOPE International Ltd. Copyright 2007-2012 by GEO-SLOPE International, Ltd. All rights reserved. No part of this work
More informationEvaluation of Flow Transmissibility of Rockfill Structures
Evaluation of Flow Transmissibility of Rockfill Structures Toshihiro MORII 1 and Takahiko TATEISHI 2 Abstract To predict the hydraulic conditions during and after the construction of such structures as
More informationEarth dam steady state seepage analysis
Engineering manual No. 32 Updated 3/2018 Earth dam steady state seepage analysis Program: FEM Water Flow File: Demo_manual_32.gmk Introduction This example illustrates an application of the GEO5 FEM module
More informationInstructor : Dr. Jehad Hamad. Chapter (7)
Instructor : Dr. Jehad Hamad Chapter (7) 2017-2016 Soil Properties Physical Properties Mechanical Properties Gradation and Structure Compressibility Soil-Water Relationships Shear Strength Bearing Capacity
More informationHydrological process simulation in the earth dam and dike by the Program PCSiWaPro
Fakultät Umweltwissenschaften, Fachrichtung Hydrowissenschaften. Hydrological process simulation in the earth dam and dike by the Program PCSiWaPro Jinxing Guo, Peter-Wolfgang Graeber Table of contents
More information*** ***! " " ) * % )!( & ' % # $. 0 1 %./ +, - 7 : %8% 9 ) 7 / ( * 7 : %8% 9 < ;14. " > /' ;-,=. / ١
١ ******!" #$ % & '!( ) % * ") +,-./ % 01. 3 ( 4 56 7/4 ) 8%9 % : 7 ;14 < 8%9 % : *7./ = ;-, >/'." Soil Permeability & Seepage ٢ Soil Permeability- Definition ٣ What is Permeability? Permeability is the
More informationSeepage. c ZACE Services Ltd. August 2011
Seepage c ZACE Services Ltd August 2011 1 / 50 2 / 50 Seepage analysis for fully/partially saturated media 1 Steady state v F k,k = 0 2 Transient v F k,k c p = 0 Darcy velocity v F i = k ij k r (S) ( p
More informationINTERPRETATION OF UNDRAINED SHEAR STRENGTH OF UNSATURATED SOILS IN TERMS OF STRESS STATE VARIABLES
INTERPRETATION OF UNDRAINED SHEAR STRENGTH OF UNSATURATED SOILS IN TERMS OF STRESS STATE VARIABLES S. K. Vanapalli and D.G. Fredlund Department of Civil Engineering University of Saskatchewan, Saskatoon
More informationNUMERICAL ANALYSIS OF DESICCATION, SHRINKAGE AND CRACKING IN LOW PLASTICITY CLAYEY SOILS
IS Numerical - Modelling analysis of Shrinkage of desiccation, Cracking shrinkage, Porous and Media cracking in low plasticity clayey soils XIV International Conference on Computational Plasticity. Fundamentals
More informationTable of Contents Chapter 1 Introduction to Geotechnical Engineering 1.1 Geotechnical Engineering 1.2 The Unique Nature of Soil and Rock Materials
Table of Contents Chapter 1 Introduction to Geotechnical Engineering 1.1 Geotechnical Engineering 1.2 The Unique Nature of Soil and Rock Materials 1.3 Scope of This Book 1.4 Historical Development of Geotechnical
More informationResearch Article One-Dimensional Vacuum Steady Seepage Model of Unsaturated Soil and Finite Difference Solution
Hindawi Mathematical Problems in Engineering Volume 2017, Article ID 9589638, 7 pages https://doi.org/10.1155/2017/9589638 Research Article One-Dimensional Vacuum Steady Seepage Model of Unsaturated Soil
More informationA void ratio dependent water retention curve model including hydraulic hysteresis
DOI: 10.1051/ e3sconf/0911010 A void ratio dependent water retention curve model including hydraulic hysteresis Amin Y. Pasha 1, Arman Khoshghalb 1 and Nasser Khalili 1,a 1 School of Civil and Environmental
More informationFreezing Around a Pipe with Flowing Water
1 Introduction Freezing Around a Pipe with Flowing Water Groundwater flow can have a significant effect on ground freezing because heat flow via convection is often more effective at moving heat than conduction
More informationC. Lanni(1), E. Cordano(1), R. Rigon(1), A. Tarantino(2)
Landslide Processes: from geomorphologic mapping to dynamic modelling (6-7 February, 2009 - Strasbourg, France) A tribute to Prof. Dr. Theo van Asch C. Lanni(1), E. Cordano(1), R. Rigon(1), A. Tarantino(2)
More informationExperiment and Modeling of Soil-Water Characteristic Curve of Unsaturated Soil in Collapsing Erosion Area
Pol. J. Environ. Stud. Vol. 25, No. 6 (2016), 2509-2517 DOI: 10.15244/pjoes/64307 Original Research Experiment and Modeling of Soil-Water Characteristic Curve of Unsaturated Soil in Collapsing Erosion
More informationTemperature dependent multiphase flow and transport
Temperature dependent multiphase flow and transport J.F. Sykes, A.G. Merry and J. Zhu Department of Civil Engineering, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1 E-mail: sykesj@uwaterloo.ca
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Title Pore-water pressure characteristics of two instrumented residual soil slopes( Main article ) Author(s)
More informationField instrumentation for performance assessment of Geobarrier System
Field instrumentation for performance assessment of Geobarrier System,a Harianto Rahardjo, Zhai Qian, Alfrendo Satyanaga, Eng Choon Leong, Chien Looi Wang, Johnny Liang Heng Wong School of Civil and Environmental
More informationSTABILITY ANALYSIS OF EARTH DAM SLOPES SUBJECTED TO EARTHQUAKE USING ERT RESULTS INTERPRETATION
STABILITY ANALYSIS OF EARTH DAM SLOPES SUBJECTED TO EARTHQUAKE USING ERT RESULTS INTERPRETATION Eko Andi Suryo Lecturer / Department of Civil Engineering, Faculty of Engineering / University of Brawijaya
More informationUpdating the Coupling Algorithm in HYDRUS Package for MODFLOW
Updating the Coupling Algorithm in HYDRUS Package for MODFLOW SAHILA BEEGUM Guided by Dr. K P Sudheer, Dr. Indumathi M Nambi & Dr. Jirka Šimunek Department of Civil Engineering, Indian Institute of Technology
More informationSeepage through a dam embankment
Elevation GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com Introduction Seepage through a dam embankment The objective of this example is look at a rather simple case of flow through
More informationAnalytical approach predicting water bidirectional transfers: application to micro and furrow irrigation
Advances in Fluid Mechanics VI 633 Analytical approach predicting water bidirectional transfers: application to micro and furrow irrigation D. Crevoisier Irrigation Research Unit, Cemagref Montpellier,
More informationGEO-SLOPE International Ltd, Calgary, Alberta, Canada Convective air flow in heated column
1 Introduction Convective air flow in heated column The objective of this illustration is show the how coupled thermal / air flow module in AIR/W and TEMP/W can result in convective cells of air movement
More informationResearch Article Performance Evaluation of Four-Parameter Models of the Soil-Water Characteristic Curve
e Scientific World Journal, Article ID 56985, 2 pages http://dx.doi.org/55/204/56985 Research Article Performance Evaluation of Four-Parameter Models of the Soil-Water Characteristic Curve Siti Jahara
More informationRem: Revista Escola de Minas ISSN: Universidade Federal de Ouro Preto Brasil
Rem: Revista Escola de Minas ISSN: 0370-4467 editor@rem.com.br Universidade Federal de Ouro Preto Brasil Portugal Menezes, Luciana; Lopes Oliveira Filho, Waldyr; Carvalho Silva, Cláudio Henrique Determination
More informationWater in Soil Sections in Craig
Water in Soil Sections 2.1-2.6 in Craig Outlines Introduction Darcy s Law Volume of water flowing per unit time Measuring K in laboratory Seepage Theory Flow Net Introduction All soils are permeable materials,
More informationFrozen Ground Containment Barrier
Frozen Ground Containment Barrier GEO-SLOPE International Ltd. www.geo-slope.com 1200, 700-6th Ave SW, Calgary, AB, Canada T2P 0T8 Main: +1 403 269 2002 Fax: +1 888 463 2239 Introduction Frozen soil barriers
More informationUpscaling of Richards equation for soil moisture dynamics to be utilized in mesoscale atmospheric models
Exchange Processes at the imnd Surface for a Range of Space and Time Scales (Proceedings of the Yokohama Symposium, July 1993). [AHS Publ. no. 212, 1993. 125 Upscaling of Richards equation for soil moisture
More information7. STRESS ANALYSIS AND STRESS PATHS
7-1 7. STRESS ANALYSIS AND STRESS PATHS 7.1 THE MOHR CIRCLE The discussions in Chapters and 5 were largely concerned with vertical stresses. A more detailed examination of soil behaviour requires a knowledge
More information(Refer Slide Time: 02:10)
Soil Mechanics Prof. B.V.S. Viswanathan Department of Civil Engineering Indian Institute of Technology, Bombay Lecture 24 Flow of water through soils-v Welcome to lecture five of flow of water through
More informationInvestigation of transient effects on the soil-water characteristic curve of different granular soils
Investigation of transient effects on the soil-water characteristic curve of different granular soils M. Milatz, T. Törzs & J. Grabe Institute of Geotechnical Engineering and Construction Management, Hamburg
More informationWUFI Workshop at NTNU /SINTEF Fundamentals
WUFI Workshop at NTNU /SINTEF 2008 Fundamentals Contents: From steady-state to transient Heat storage and -transport Moisture storage and -transport Calculation of coupled transport Model limitations 2
More informationCan we distinguish Richards and Boussinesq s equations for hillslopes?: The Coweeta experiment revisited
WATER RESOURCES RESEARCH, VOL. 35, NO. 2, PAGES 589 593, FEBRUARY 1999 Can we distinguish Richards and Boussinesq s equations for hillslopes?: The Coweeta experiment revisited T. S. Steenhuis, 1 J.-Y.
More informationSeepage Analysis for Shurijeh Reservoir Dam Using Finite Element Method. S. Soleymani 1, A. Akhtarpur 2
Seepage Analysis for Shurijeh Reservoir Dam Using Finite Element Method S. Soleymani 1, A. Akhtarpur 2 1 Group of Dam Construction, Toossab Company, P.O. Box 917751569, Mashhad City, Iran, PH (+98) 511-7684091;
More informationAnalysis of Transient Seepage Through Levees
Analysis of Transient Seepage Through Levees Matthew David Sleep Dissertation submitted to the faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements
More informationAn Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation. Ahsan Munir, PhD Tom Spirka, PhD
An Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation Ahsan Munir, PhD Tom Spirka, PhD Agenda Provide an overview of COMSOL 4.3b Our products, solutions and applications Subsurface
More informationThe process of consolidation and settlement
Consolidation Based on part of the GeotechniCAL reference package by Prof. John Atkinson, City University, London The process of consolidation and settlement One-dimensional consolidation theory The oedometer
More informationTIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL
TIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL Qassun S. Mohammed Shafiqu and Maarib M. Ahmed Al-Sammaraey Department of Civil Engineering, Nahrain University, Iraq
More informationCHAPTER 2. SOIL-WATER POTENTIAL: CONCEPTS AND MEASUREMENT
SSC107 Fall 2000 Chapter 2, Page - 1 - CHAPTER 2. SOIL-WATER POTENTIAL: CONCEPTS AND MEASUREMENT Contents: Transport mechanisms Water properties Definition of soil-water potential Measurement of soil-water
More informationAPPLICATION OF TWO-PHASE REGRESSION TO GEOTECHNICAL DATA. E. Stoimenova, M. Datcheva, T. Schanz 1
Pliska Stud. Math. Bulgar. 16 (2004), 245-257 STUDIA MATHEMATICA BULGARICA APPLICATION OF TWO-PHASE REGRESSION TO GEOTECHNICAL DATA E. Stoimenova, M. Datcheva, T. Schanz 1 A method for estimating a transition
More informationEffect of long duration rainstorm on stability of Red-clay slopes
Xue et al. Geoenvironmental Disasters (2016) 3:12 DOI 10.1186/s40677-016-0046-9 Geoenvironmental Disasters RESEARCH Open Access Effect of long duration rainstorm on stability of Red-clay slopes Kaixi Xue
More informationCHARACTERISTICS OF SEDIMENT TRANSPORT IN SWASH ZONE DUE TO SATURATED-UNSATURATED SLOPED BEACH
CHARACTERISTICS OF SEDIMENT TRANSPORT IN SWASH ZONE DUE TO SATURATED-UNSATURATED SLOPED BEACH Masashi Ochi 1, Makoto Miyatake 2 and Katsutoshi Kimura 3 The influence of saturated-unsaturated sloped beach
More informationSTUDY OF THE BARCELONA BASIC MODEL. INFLUENCE OF SUCTION ON SHEAR STRENGTH
STUDY OF TH BARCLONA BASIC MODL. INFLUNC OF SUCTION ON SHAR STRNGTH Carlos Pereira ABSTRACT The Barcelona Basic Model, BBM, is one of the most used elasto-plastic models for unsaturated soils. This summary
More informationStudies of rainfall-induced slope failures
Invited Lecture. National Seminar Slope 2000, 27 April 2002, Bandung, Indonesia. In Paulus P Rahardjo (edt.) Slope 2002. Proceedings of the National Seminar, Slope 2002. 27-April 2002. Bandung, Indonesia.
More informationTime Rate of Consolidation Settlement
Time Rate of Consolidation Settlement We know how to evaluate total settlement of primary consolidation S c which will take place in a certain clay layer. However this settlement usually takes place over
More informationExperimental Analysis on Soil-Water Characteristic Curve of CH3COO - Contaminated Clay
Experimental Analysis on Soil-Water Characteristic Curve of CH3COO - Contaminated Clay Liwen Cao, Tianyu Xu, Yong Wang and Pan Huo School of Resource and Earth Science, China University of Mining and Technology
More informationSubsurface Flow Modelling At The Hillslope Scale: Numerical And Physical Analysis. Oscar Cainelli
Subsurface Flow Modelling At The Hillslope Scale: Numerical And Physical Analysis Oscar Cainelli 2007 Based on the doctoral thesis in Environmental Engineering (XIX cycle) defended in February 2007 at
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationUnsaturated seepage behavior study using soil column test
Unsaturated seepage behavior study using soil column test Author Kim, Dong Hyun, Kim, S., Gratchev, Ivan Published 2013 Conference Title Proceedings of 18th Southeast Asian Geotechnical cuminaugural AGSSEA
More informationDarcy's Law. Laboratory 2 HWR 531/431
Darcy's Law Laboratory HWR 531/431-1 Introduction In 1856, Henry Darcy, a French hydraulic engineer, published a report in which he described a series of experiments he had performed in an attempt to quantify
More informationHydraulic conductivity of granular materials
3 r d International Conference on New Developments in Soil Mechanics and Geotechnical Engineering, Hydraulic conductivity of granular materials Namir K.S.Al-Saoudi Building and construction Eng. Dept.
More informationThermal and hydraulic modelling of road tunnel joints
Thermal and hydraulic modelling of road tunnel joints Cédric Hounyevou Klotoé 1, François Duhaime 1, Lotfi Guizani 1 1 Département de génie de la construction, École de technologie supérieure, Montréal,
More informationClimate effects on landslides
GEORAMP ONE DAY SYMPOSIUM Climate effects on landslides E. E. Alonso, M. Sondón, N. M. Pinyol Universitat Politècnica de Catalunya October 14th, 2016. UPC, Barcelona Infiltration (evaporation) and slope
More informationConsolidation. Verification Manual GT. Written by: The SoilVision Systems Ltd. Team. Last Updated: Saturday, October 13, 2018
Consolidation Verification Manual GT Written by: The SoilVision Systems Ltd. Team Last Updated: Saturday, October 13, 2018 SoilVision Systems Ltd. Saskatoon, Saskatchewan, Canada Software License The software
More informationStochastic geometry and porous media
Stochastic geometry and transport in porous media Hans R. Künsch Seminar für Statistik, ETH Zürich February 15, 2007, Reisensburg Coauthors Thanks to the coauthors of this paper: P. Lehmann, A. Kaestner,
More information&
Inverse Problem for Transient Unsaturated Flow Identifiability and Non Uniqueness K. S. Hari Prasad & M. S. Ghidaoui Department of Civil and Structural Engineering The Hong Kong University of Science and
More information1 BASIC CONCEPTS AND MODELS
1 BASIC CONCEPTS AND ODELS 1.1 INTRODUCTION This Volume III in the series of textbooks is focused on applications of environmental isotopes in surface water hydrology. The term environmental means that
More informationA Critical Review of the Methodologies Employed for Suction Measurement for Developing the SWCC
The 12 th International Conference of International Association for Computer Methods and Advances in Geomechanics (IACMAG) 1-6 October, 08 Goa, India A Critical Review of the Methodologies Employed for
More informationElectrical Resistivity of Compacted Kaolin and its Relation with Suction
Electrical Resistivity of Compacted Kaolin and its Relation with Suction Dias, Ana Sofia ana.sofia.dias@tecnico.ulisboa.pt Summary The electrical characteristics of compacted kaolin were studied and related
More informationRapid Drawdown Stability Analysis of San Luis Dam
Rapid Drawdown Stability Analysis of San Luis Dam Stark, T.D., tstark@illinois.edu Department of Civil and Environmental Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801 Jafari,
More informationInfluence of Surrounding Soil Conditions and Joint Sealing on Seepage Resistance of a Sheet Pile Wall, Three Dimensional Numerical Analyses
IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Influence of Surrounding Soil Conditions and Joint Sealing on Seepage Resistance of a Sheet Pile Wall, Three Dimensional Numerical
More informationA novel modeling approach for the simulation of soil water interaction in a highly plastic clay
Geomech. Geophys. Geo-energ. Geo-resour. (216) 2:77 95 DOI 1.17/s4948-16-23-5 ORIGINAL ARTICLE A novel modeling approach for the simulation of soil water interaction in a highly plastic clay Ramy Saadeldin.
More informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More informationNumerical Modeling of Rainfall-induced Slope Failure
Numerical Modeling of Rainfall-induced Slope Failure An Thi Phuong TRAN 1), Youngjong SIM 2), Ah-Ram KIM 3), Gye-Chun CHO 4), and Dae-Hong KIM 5) 1), 3), 4) Department of Civil and Environmental Engineering,
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Title Horizontal drains in residual soil slopes( Main article ) Author(s) Citation Santoso, Vera Amalia; Rahardjo,
More informationNumerical Solution of the Two-Dimensional Time-Dependent Transport Equation. Khaled Ismail Hamza 1 EXTENDED ABSTRACT
Second International Conference on Saltwater Intrusion and Coastal Aquifers Monitoring, Modeling, and Management. Mérida, México, March 3-April 2 Numerical Solution of the Two-Dimensional Time-Dependent
More informationExperimental measurement of parameters governing flow rates and partial saturation in paper-based microfluidic devices
Experimental measurement of parameters governing flow rates and partial saturation in paper-based microfluidic devices Dharitri Rath 1, Sathishkumar N 1, Bhushan J. Toley 1* 1 Department of Chemical Engineering
More informationRATE OF FLUID FLOW THROUGH POROUS MEDIA
RATE OF FLUID FLOW THROUGH POROUS MEDIA Submitted by Xu Ming Xin Kiong Min Yi Kimberly Yip Juen Chen Nicole A project presented to the Singapore Mathematical Society Essay Competition 2013 1 Abstract Fluid
More informationSlope Stability Model of the Questa Rock Pile Phase 2
2 Proceedings Tailings and Mine Waste 2011 Slope Stability Model of the Questa Rock Pile Phase 2 Murray Fredlund SoilVision Systems Ltd., Saskatoon, Canada Haihua Lu SoilVision Systems Ltd., Saskatoon,
More information