&
|
|
- Chad Skinner
- 5 years ago
- Views:
Transcription
1 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 Technology hari@wcmailust.hk & GHIDAOUI@usthk.ust.hk Abstract In the present study, the identifiability and non uniqueness of soil hydraulic properties in unsaturated flow systems are discussed. Necessary and sufficient conditions for identifiability are derived. A parameter estimation algorithm is developed by combining the Adjoint State Theory with Conjugate Gradient Method and Levenberg-Marquardt method. These two methods are compared through few numerical examples. It is shown that Levenberg-Marquardt algorithm is superior to Conjugate Gradient Method. However, Levenberg-Marquardt algorithm also fails to converge to the global minimum when the initial guess parameters are far from the true values. Introduction Moisture flow through the unsaturated zone is usually analysed by solving Richards Equation [4]. Solution of Richards equation requires knowledge of soil hydraulic conductivity and water content versus water pressure functions. In recent years, many researchers have investigated the feasibility of determining these functions simultaneously from transient flow by numerical inversion of the governing initial value problem. In such approaches, the unknown parameters are estimated by minimizing deviations between the observed and predicted flow attributes such as water contents, pressure heads or infiltration rates. The inverse problem is often ill-posed which results in non unique and unstable solutions. Application of parameter estimation in unsaturated
2 34 Computer Methods in Water Resources XII flow are few, consisting analysis of inflow and outflow experiments on laboratory soil cores and some field studies. Kool et al.[3] reviewed the parameter estimation methods in the unsaturated zone. The objective of the present study is to (i) develop theoretical analysis for identifiability of parameters used in different soil hydraulic relationships (ii) to develop a parameter estimation algorithm using Adjoint State Method and Conjugate Gradient Method and Levenberg - Marquardt methods (hi) performing numerical studies to test the robustness of these methods. Governing Equation The pressure head based formulation of Richards equation for one dimensional, single phase, non hysteretic transient vertical flow of water through a rigid and homogeneous unsaturated soil can be written as (Celiaetal.[l]) where ij./ is the pressure head, K is the hydraulic conductivity, C is the soil moisture capacity, 9 is the volumetric moisture content, z is the vertical co-ordinate taken positive upwards, t is the time co-ordinate and L is the length of the soil. In order to solve eqn(l), constitutive relationships between vj/ and nonlinear terms K and C have to be specified. The empirical relationships proposed by Gardner[2] and Van Genuchten[7] are considered in the present study. Gardner's Relationship : Van Genuchten's Relationship : (2) (3) for H/<0 ^ = for > > 0
3 Computer Methods in Water Resources XII 35 (5) In eqns (2)-(5), (3, a and n are unsaturated soil parameters with m = - l/«, IQat is the saturated hydraulic conductivity and S^ is the effective saturation defined as where 6^ is the saturated moisture content and 0, is the residual moisture content. Identifiability Identifiability refers to the forward problem. The parameters are said to be unidentifiable if more than one set of model parameters produce identical outputs. If two sets (K^CJ and (K^A) produce the same responses i then The parameters are not identifiable, if 5K * 0 and 5C * 0. A necessary, but not sufficient condition for identifiability is 5K = 0 and 5C = 0. Case : dj/dt = 0 and (dy/dz + ) = 0 In this case, i = -z, which corresponds to a no flow problem. In such a case (8K,5C) need not necessarily be equal to (0,0) making K and C unidentifiable. Since, there is no flow in and out of the system, no information can be obtained from such a data for identifying K and C. Case 2 : dj/dt = 0 and (dy/dz + ) * 0 In this case 8K = 0 and 6C need not be zero indicating that the steady state data is not capable of identifying the parameter C. However, the parameter K can be identified from such type of data.
4 36 Computer Methods in Water Resources XII Case 3 : dv //dt # 0 and (di/dz + ) = 0 In this case v / = -z, which usually does not arise in transient unsaturated flow. Hence this type of data is also not useful for identifying the parameters. Case 4 : dy/dt # 0 and (df/dz + ) * 0 In this case, the data is the transient flow data. In such a case, for identifiability, the necessary condition is 5K = 0 and 5C = 0. Thus, the transient data can be used to identify the parameters K and C. However, is 5K - 0 and 5C = 0 a sufficient condition for identifiability? To illustrate this consider the following two sub cases. Case 4a Let us assume that K and C depend on i and two parameters PI and P2- Then 5K = 8P, + 6P? = 0 5P, ' 8P; ' dc dc 8C = ^ W^k 8P, + <^TTk 5P,=0 / (8) / Generally, (9K/9Pi,9K/9P2) and (9C/9Pi,6C/9P2) are two linearly independent vectors. For eqn(8) to be true, we should have 5Pi = 0 and 8?2 = 0, i.e., PI and?2 are identifiable. In the case of Gardner's model K and C depend upon two parameters K^t and p. Let the ratios of dc/dk^ and 9K/(3Ksat & 9C/3(3 and 3K/3(3 be denoted by aj and a2 respectively. If (5C/5Ksat,3C/ap) and (SK/dK^SK/dp) are linearly independent, then a2 = Aai, where A is a constant. It is clear from eqn(2) that aj is zero. In figs.(l) and (2), the ratios a^ and &2 are plotted against pressure head for different values of K^ and P respectively. It is clear from these figures that at pressures very near to saturation, these two curves tend to become parallel to each other, indicating the linear dependence. Hence, if the data consists of pressure heads near to saturation, K^t and P may not be estimated accurately. Case 4b Let us assume that K and C depend upon j and three parameters PI,?2 and Pg. In such a case 0P, a?2 6P; 6K - 8P, + 5?2 + SPj = 0 (9) 5P ' a* ^ ap ^
5 Computer Methods in Water Resources XII 37 This system of equations will either have no solution or will have many possible solutions. Hence, in this case, all the three parameters are not identifiable simultaneously. However, fixing one of these parameters enables one to identify the other two parameters. Van Genuchten's relationship has three parameters K^t, a and n for describing the soil hydraulic properties. These three parameters can't be identified simultaneously from transient flow data. Inverse Problem The commonly used objective function G (Russo et al.[5])can be written as : LT 00 ^" (io) where L is the length of the soil profile, T is the total time of simulation, Qm and qc are the measured and computed infiltration rates. The parameter estimation involves the estimation of the parameters K^t, a and n for Van Genuchten relationship. The sensitivity coefficients at each iteration of the parameter estimation are obtained using Adjoint State Theory. Following Sun and Yeh[6], the adjoint equation for eqns () and (0) can be obtained as +AK(^LO,t.. = 0, t = T,0<z<L; X = 0, t < T, z = 0..^ T () The adjoint problem consists of solving for A, after solving the direct problem for j. Conjugate Gradient and Levenberg - Marquardt methods are used for the minimization. To test the robustness of these methods, several numerical tests are performed. A hypothetical data consisting of infiltration rates is generated by solving eqn(l), with the following set of parameters; K^t = 5cm/hr, a = 0.02 cm"% % = 2.3, 6^ = 0.38, 8, = The length L of the soil profile is considered as m and the initial moisture content(8j throughout the soil profile is taken as 0.2. The moisture content at the bottom is held at 0.2 while a moisture content equal to 9^ is applied at the top. Infiltration data are obtained by solving the direct problem (eqn(l)) for a period of about 8 minutes. This data
6 38 Computer Methods in Water Resources XII are used as measured infiltration rates in all the test runs. 8g and 9,. are considered as constants during the parameter estimation. With only one parameter (K^t, oc or n) as an unknown, both Conjugate Gradient method and Levenberg - Marquardt algorithm converge to the true solution. A comparison of Conjugate Gradient method and Levenberg-Marquardt algorithm is presented in Table.. Table. Comparison of Conjugate Gradient and Levenberg - Marquardt Methods Iteration Conjugate Gradient Method Levenberg - Marquardt Method Ksatrinn = 9cm/hr, K^truei = 5cm/hr tt(ini) ~ 0.07 Cm"ttftrue) ~ 0.02 Cm" "(inn = 3.5 "(true) = It is clear from Table. that the Levenberg - Marquardt algorithm converges to the true solution in fewer iterations as compared to Conjugate Gradient Method. For the case, when two parameters are unknown, Conjugate Gradient method fails to converge to global minimum where as Levenberg - Marquardt algorithm converges to the true values. However, Levenberg - Marquardt algorithm also fails to converge, when the initial parameters are far away from the true parameters. Fig.3 shows the parametric space through which the global
7 Computer Methods in Water Resources XII 39 minimum can be obtained for the case in which K^t and a are considered as unknown parameters. It is clear from Fig. 3 that the non uniqueness is quite pronounced in the estimation of parameters in unsaturated flow systems. Conclusions It has been shown that, with the unsteady flow data, the soil hydraulic model with two parameters can be identified uniquely, where as a model with three parameters are unidentfiable. Adjoint equations are derived for the Richards equation and are coupled with Conjugate Gradient and Levenberg - Marquardt algorithm for the estimation of soil parameters. Levenberg - Marquardt algorithm is found to be superior to Conjugate Gradient Method in estimating the parameters. References. Celia, M.A., Bouloutas, E.T. and Zarba, R.L., A General Mass Conservative Numerical Solution for the Unsaturated Flow Equation, z, 26(7), pp , Gardner, W. R., Some Steady State Solutions of Unsaturated Moisture Flow Equations with Application to Evaporation from a Watertable, Soil Science, 85, pp , Kool, J.B., Parker, P.C. and Van Genuchten, M.T., Parameter Estimation for Unsaturated Flow and Transport Models - A Review, Journal of Hydrology, 9, pp Richards, L.A., Capillary Conduction of Liquids Through Porous Medium, Physics,, pp , Russo, D., Bresler, E., Shani,U. and Parker, J.C., Analysis of Infiltration Events in Relation to Determining Soil Hydraulic Properties by Inverse Problem Methodology, Water Resources z, 27(6),pp , 99. Sun, N.Z. and Yeh, W.W.G., Coupled Inverse Problems in Groundwater Modeling,. Sensitivity Analysis and Parameter Identification, Water Resources Research, 26(0), pp , 990.
8 320 Computer Methods in Water Resources XII 7. Van Genuchten, M. T., A Closed Form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils, 6W &v'gmce 6bcz'e(y of 'ca JowrW, 44, pp , tt H ^ in ct re -02- nt m -... i I a, (=0) (3 0.0cm"'' P 0.02 cm"^ 0.03 cm'' 0.04 cm'' p = 0.05 cm'' A " > ' / W' / - /'?^ : // / / V / Pressure Head (cm) Fig. : Variation of a, and a% with pressure head - K^ = 3cm/hr Pressure Head (cm) Fig. 2 : Variation of a, and a% with pressure head - p = 0.03 cm"'' Global minimum Fig. 3 : Parametric Space of K^ and a - Global and Local Convergence
Comparison 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 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 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 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 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 informationIn situ estimation of soil hydraulic functions using a multistep soil-water extraction technique
WATER RESOURCES RESEARCH, VOL. 34, NO. 5, PAGES 1035 1050, MAY 1998 In situ estimation of soil hydraulic functions using a multistep soil-water extraction technique M. Inoue Arid Land Research Center,
More informationAnalytical solutions for water flow and solute transport in the unsaturated zone
Models for Assessing and Monitoring Groundwater Quality (Procsedines of a Boulder Symposium July 1995). IAHS Publ. no. 227, 1995. 125 Analytical solutions for water flow and solute transport in the unsaturated
More informationDirect Hydraulic Parameter and Function Estimation for Diverse Soil Types Under Infiltration and Evaporation
Transp Porous Med (2017) 116:797 823 DOI 10.1007/s11242-016-0801-0 Direct Hydraulic Parameter and Function Estimation for Diverse Soil Types Under Infiltration and Evaporation Jianying Jiao 1 Ye Zhang
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 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 informationIn 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 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 informationSolute dispersion in a variably saturated sand
WATER RESOURCES RESEARCH, VOL. 39, NO. 6, 1155, doi:10.1029/2002wr001649, 2003 Solute dispersion in a variably saturated sand Takeshi Sato Department of Civil Engineering, Gifu University, Yanagido, Gifu,
More informationNumerical Simulations of Two-Phase Flow in Rigid Porous Media
University of Colorado, Boulder CU Scholar Civil Engineering Graduate Theses & Dissertations Civil, Environmental, and Architectural Engineering Spring 1-1-2017 Numerical Simulations of Two-Phase Flow
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 informationCoupling TRIGRS and TOPMODEL in shallow landslide prediction. 1 Presenter: 王俊皓 Advisor: 李錫堤老師 Date: 2016/10/13
Coupling TRIGRS and TOPMODEL in shallow landslide prediction 1 Presenter: 王俊皓 Advisor: 李錫堤老師 Date: 016/10/13 Outline Introduction Literature review Methodology Pre-result Future work Introduction 3 Motivation
More informationPrinciples of soil water and heat transfer in JULES
Principles of soil water and heat transfer in JULES Anne Verhoef 1, Pier Luigi Vidale 2, Raquel Garcia- Gonzalez 1,2, and Marie-Estelle Demory 2 1. Soil Research Centre, Reading (UK); 2. NCAS-Climate,
More informationAbstract. This research is supported by US Army Corps of Engineers under Grant No. DACA K-0055 with Penn State.
Modeling Water Flow and Chemical and Sediment Transport in Watershed Systems Hwai-Ping Cheng and Gour-Tsyh Yeh Department of Civil and Environmental Engineering, The Pennsylvania State University, University
More information6.6 Solute Transport During Variably Saturated Flow Inverse Methods
6.6 Solute Transport During Variably Saturated Flow Inverse Methods JIÌÍ ŠIMæNEK AND MARTINUS TH. VAN GENUCHTEN, USDA-ARS, George E. Brown, Jr. Salinity Laboratory, Riverside, California DIEDERIK JACQUES,
More informationAn investigation of temporal adaptive solution of Richards equation for sharp front problems
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn:2319-765x. Volume 10, Issue 2 Ver. II (Mar-Apr. 2014), PP 106-117 An investigation of temporal adaptive solution of Richards equation for
More informationStability Analysis of Landslide Dam under Rainfall
Stability Analysis of Landslide Dam under Rainfall Pei-Hsun Tsai, Zheng-Yi Feng 2, Fan-Chieh Yu 3 and Jian-Han Lin 4 Associate Professor, Department of Construction Engineering, Chaoyang University of
More informationComparison of conductivity averaging methods for one-dimensional unsaturated flow in layered soils
A. Szymkiewicz a R. Helmig b Comparison of conductivity averaging methods for one-dimensional unsaturated flow in layered soils Stuttgart, October 1 a Faculty of Civil and Environmental Engineering, Gdansk
More informationSecond Order Linear Equations
Second Order Linear Equations Linear Equations The most general linear ordinary differential equation of order two has the form, a t y t b t y t c t y t f t. 1 We call this a linear equation because the
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 informationHydraulic tomography: Development of a new aquifer test method
WATER RESOURCES RESEARCH, VOL. 36, NO. 8, PAGES 2095 2105, AUGUST 2000 Hydraulic tomography: Development of a new aquifer test method T.-C. Jim Yeh and Shuyun Liu Department of Hydrology and Water Resources,
More informationDetermination of Parameters for Bimodal Hydraulic Functions by Inverse Modeling
Determination of Parameters for Bimodal Hydraulic Functions by Inverse Modeling Torsten Zurmühl and Wolfgang Durner* ABSTRACT content,, with respect to the logarithm of matric head, To determine the hydraulic
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 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 informationAn application of nested Newton-type algorithm for finite difference method solving Richards' equation
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn:239-765x. Volume, Issue Ver. II. (Feb. 24), PP 2-32 An application of nested Newton-type algorithm for finite dference method solving Richards'
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 informationIdentification of Parameters of the Basic Hydrophysical Characteristics of Soil
Identification of Parameters of the Basic Hydrophysical Characteristics of Soil Elena S. Zasukhina Dorodnicyn Computing Centre, FRC CSC RAS Vavilov st. 40, 119333 Moscow, Russia. elenazet2017@yandex.ru
More informationOne-Dimensional Steady State Infiltration in Heterogeneous Soils
WATER RESOURCES RESEARCH, VOL. 25, NO. 10, PAGES 2149-2158, OCTOBER 1989 One-Dimensional Steady State Infiltration in Heterogeneous Soils T.-C. JIM YEI Department of Hydrology and Water Resources, University
More information12 SWAT USER S MANUAL, VERSION 98.1
12 SWAT USER S MANUAL, VERSION 98.1 CANOPY STORAGE. Canopy storage is the water intercepted by vegetative surfaces (the canopy) where it is held and made available for evaporation. When using the curve
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 informationSecond-Order Linear ODEs (Textbook, Chap 2)
Second-Order Linear ODEs (Textbook, Chap ) Motivation Recall from notes, pp. 58-59, the second example of a DE that we introduced there. d φ 1 1 φ = φ 0 dx λ λ Q w ' (a1) This equation represents conservation
More informationA hybrid Marquardt-Simulated Annealing method for solving the groundwater inverse problem
Calibration and Reliability in Groundwater Modelling (Proceedings of the ModelCARE 99 Conference held at Zurich, Switzerland, September 1999). IAHS Publ. no. 265, 2000. 157 A hybrid Marquardt-Simulated
More informationRecent Advances in Profile Soil Moisture Retrieval
Recent Advances in Profile Soil Moisture Retrieval Jeffrey P. Walker, Garry R. Willgoose and Jetse D. Kalma Department of Civil, Surveying and Environmental Engineering The University of Newcastle, Callaghan,
More information2 integration methods estimate and control the temporal error through adapting both the temporal step size and the order of the temporal approximation
1 Higher order, locally conservative temporal integration methods for modeling multiphase flow in porous media C.T. Miller aλ, M.W. Farthing ay, C.E. Kees a z, and C.T. Kelley bx a Center for the Advanced
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 informationFlow in unsaturated random porous media, nonlinear numerical analysis and comparison to analytical stochastic models
Advances in Water Resources Vol. 22, No. 3, pp. 257 272, 1998 1998 Elsevier Science Ltd Printed in Great Britain. All rights reserved PII:S0309-1708(98)00010-4 0309-1708/98/$ - see front matter Flow in
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 informationHYDRAULIC CONDUCTIVITY ESTIMATION IN PARTIALLY SATURATED SOILS USING THE ADJOINT METHOD
HYDRAULIC CONDUCTIVITY ESTIMATION IN PARTIALLY SATURATED SOILS USING THE ADJOINT METHOD J. SANTOS, Y. EFENDIEV, AND L. GUARRACINO Abstract An iterative algorithm based on the adjoint method for the estimation
More informationNumerical investigations of hillslopes with variably saturated subsurface and overland flows
Numerical investigations of hillslopes with variably saturated subsurface and overland flows ARC DYNAS H. Beaugendre, T. Esclaffer, A. Ern and E. Gaume DYNAS Workshop 06-08/12/04 DYNAS Workshop 06-08/12/04
More informationAPPLICATION OF 1D HYDROMECHANICAL COUPLING IN TOUGH2 TO A DEEP GEOLOGICAL REPOSITORY GLACIATION SCENARIO
PROCEEDINGS, TOUGH Symposium 2015 Lawrence Berkeley National Laboratory, Berkeley, California, September 28-30, 2015 APPLICATION OF 1D HYDROMECHANICAL COUPLING IN TOUGH2 TO A DEEP GEOLOGICAL REPOSITORY
More information1 Robust Solution of Richards' Equation for Non-Uniform Porous Media Cass T. Miller and Glenn A. Williams Department of Environmental Sciences and Eng
1 Robust Solution of Richards' Equation for Non-Uniform Porous Media Cass T. Miller and Glenn A. Williams Department of Environmental Sciences and Engineering, University of North Carolina, Chapel Hill,
More informationApplication of Fuzzy Logic and Uncertainties Measurement in Environmental Information Systems
Fakultät Forst-, Geo- und Hydrowissenschaften, Fachrichtung Wasserwesen, Institut für Abfallwirtschaft und Altlasten, Professur Systemanalyse Application of Fuzzy Logic and Uncertainties Measurement in
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 informationApplied Mathematics 205. Unit I: Data Fitting. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit I: Data Fitting Lecturer: Dr. David Knezevic Unit I: Data Fitting Chapter I.4: Nonlinear Least Squares 2 / 25 Nonlinear Least Squares So far we have looked at finding a best
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #4: Continuity and Flow Nets
1.7, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #4: Continuity and Flow Nets Equation of Continuity Our equations of hydrogeology are a combination of o Conservation of mass o Some empirical
More informationA THIRD-ORDER NUMERICAL SCHEME WITH UPWIND WEIGHTING FOR SOLVING THE SOLUTE TRANSPORT EQUATION
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, VOL. 40, 1623 1637 (1997) A THIRD-ORDER NUMERICAL SCHEME WITH UPWIND WEIGHTING FOR SOLVING THE SOLUTE TRANSPORT EQUATION KANGLE HUANG, JIR[ I
More informationNumerical fractional flow modelling of inhomogeneous
Numerical fractional flow modelling of inhomogeneous air sparging E.F. Kaasschieter\ G.J. Mulder^ & J.D. van der Werff ten Bosch^ Department of Mathematics and Computing Science, hoven University of Technology,
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 informationChapter Seven. For ideal gases, the ideal gas law provides a precise relationship between density and pressure:
Chapter Seven Horizontal, steady-state flow of an ideal gas This case is presented for compressible gases, and their properties, especially density, vary appreciably with pressure. The conditions of the
More information1D Verification Examples
1 Introduction 1D Verification Examples Software verification involves comparing the numerical solution with an analytical solution. The objective of this example is to compare the results from CTRAN/W
More informationDOCUMENTATION FOR PREPARING THE INPUT FILE OF THE UPDATED HYDRUS PACKAGE FOR MODFLOW-2000
DOCUMENTATION FOR PREPARING THE INPUT FILE OF THE UPDATED HYDRUS PACKAGE FOR MODFLOW-2000 Updated May 2008 Navin Kumar C. Twarakavi 1, Hyeyoung Sophia Seo 2 and Jirka Simunek 1 1 University of California,
More informationSuction Potential and Water Absorption from Periodic Channels in Different Types of Homogeneous Soils
I. Solekhudin and K.C. Ang/ Electronic Journal of Boundary Elements, Vol. 1, No. 2, pp. 42-55 (212) Suction Potential and Water Absorption from Periodic Channels in Different Types of Homogeneous Soils
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 informationEffective unsaturated hydraulic conductivity for one-dimensional structured heterogeneity
WATER RESOURCES RESEARCH, VOL. 41, W09406, doi:10.1029/2005wr003988, 2005 Effective unsaturated hydraulic conductivity for one-dimensional structured heterogeneity A. W. Warrick Department of Soil, Water
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 informationJournal of Computational and Applied Mathematics. An adaptive time discretization of the classical and the dual porosity model of Richards equation
Journal of Computational and Applied Mathematics 233 (2010) 3167 3177 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
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 informationSOLUTE TRANSPORT. Renduo Zhang. Proceedings of Fourteenth Annual American Geophysical Union: Hydrology Days. Submitted by
THE TRANSFER FUNCTON FOR SOLUTE TRANSPORT Renduo Zhang Proceedings 1994 WWRC-94-1 1 n Proceedings of Fourteenth Annual American Geophysical Union: Hydrology Days Submitted by Renduo Zhang Department of
More informationNumerical Issues Arising in the Simulations of Transient Water Flow. in Layered Unsaturated Soils. Ruowen Liu
Numerical Issues Arising in the Simulations of Transient Water Flow in Layered Unsaturated Soils by Ruowen Liu A Dissertation Presented in Partial Fulfillment of the Requirements for the Degree Doctor
More informationPROBABILISTIC PREDICTIONS OF INFILTRATION INTO HETEROGENEOUS MEDIA WITH UNCERTAIN HYDRAULIC PARAMETERS
International Journal for Uncertainty Quantification, 1(1):35 47, 211 PROBABILISTIC PREDICTIONS OF INFILTRATION INTO HETEROGENEOUS MEDIA WITH UNCERTAIN HYDRAULIC PARAMETERS Peng Wang & Daniel M. Tartakovsky
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
More informationENGCOMP - Computational Engineering
Coordinating unit: Teaching unit: Academic year: Degree: ECTS credits: 2015 250 - ETSECCPB - Barcelona School of Civil Engineering 751 - DECA - Department of Civil and Environmental Engineering MASTER'S
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 informationWorkshop Opening Presentation of the DYNAS Project
Workshop Opening Presentation of the DYNAS Project Alexandre Ern, Eric Gaume, Cyril Kao, Jérôme Jaffré CERMICS École Nationale des Ponts et Chausées Marne-la-Vallée, France DYNAS Workshop 06-08/12/04 p.1
More informationA finite analytic algorithm for modeling variably saturated flows. Duan a,b, Li Chen a,b. Chang an University, Ministry of Education, P. R.
1 2 3 A finite analytic algorithm for modeling variably saturated flows Zaiyong Zhang a,b,c, Wenke Wang a,b,*, Chengcheng Gong a,b, Zhoufeng Wang a,b, Lei Duan a,b, Li Chen a,b 4 5 6 7 8 9 10 a Key Laboratory
More informationUpscaled flow and transport properties for heterogeneous unsaturated media
WATER RESOURCES RESEARCH, VOL. 38, NO. 5, 1053, 10.1029/2000WR000072, 2002 Upscaled flow and transport properties for heterogeneous unsaturated media Raziuddin Khaleel Fluor Federal Services, Richland,
More informationLEAST SQUARES PROBLEMS WITH INEQUALITY CONSTRAINTS AS QUADRATIC CONSTRAINTS
LEAST SQUARES PROBLEMS WITH INEQUALITY CONSTRAINTS AS QUADRATIC CONSTRAINTS JODI MEAD AND ROSEMARY A RENAUT. Introduction. The linear least squares problems discussed here are often used to incorporate
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 informationAn Equivalent Circuit Formulation of the Power Flow Problem with Current and Voltage State Variables
An Equivalent Circuit Formulation of the Power Flow Problem with Current and Voltage State Variables David M. Bromberg, Marko Jereminov, Xin Li, Gabriela Hug, Larry Pileggi Dept. of Electrical and Computer
More informationFluid Mechanics Introduction
Fluid Mechanics Introduction Fluid mechanics study the fluid under all conditions of rest and motion. Its approach is analytical, mathematical, and empirical (experimental and observation). Fluid can be
More informationNon-polynomial Least-squares fitting
Applied Math 205 Last time: piecewise polynomial interpolation, least-squares fitting Today: underdetermined least squares, nonlinear least squares Homework 1 (and subsequent homeworks) have several parts
More informationEffects of Soil Heterogeneity on One-Dimensional Infiltration
Utah State University DigitalCommons@USU Reports Utah Water Research Laboratory January 1977 Effects of Soil Heterogeneity on One-Dimensional Infiltration Charles Courette Roland W. Jeppson Follow this
More informationSimulation of Nondifferentiable Models for Groundwater Flow and Transport
1 Simulation of Nondifferentiable Models for Groundwater Flow and Transport C. T. Kelley a K. R. Fowler b C. E. Kees a a Department of Mathematics and Center for Research in Scientific Computation North
More informationUNIT 4 CHEMICAL KINETICS
Concentration Ankit Gupta Classes UNIT 4 CHEMICAL KINETICS MARK QUESTIONS Q.. In the reaction A B, if the concentration of A is plotted against time, the nature of the curve obtained will be as shown.
More informationOn Parameter Estimation for Neuron Models
On Parameter Estimation for Neuron Models Abhijit Biswas Department of Mathematics Temple University November 30th, 2017 Abhijit Biswas (Temple University) On Parameter Estimation for Neuron Models November
More informationA MATLAB Implementation of the. Minimum Relative Entropy Method. for Linear Inverse Problems
A MATLAB Implementation of the Minimum Relative Entropy Method for Linear Inverse Problems Roseanna M. Neupauer 1 and Brian Borchers 2 1 Department of Earth and Environmental Science 2 Department of Mathematics
More informationCapillary effect on watertable fluctuations in unconfined aquifers
Capillary effect on watertable fluctuations in unconfined aquifers Jun Kong 1, Cheng-Ji Shen, Pei Xin, Zhiyao Song 3, Ling Li,1,#, D.A. Barry 4, D.-S. Jeng 5, F. Stagnitti 6, D.A. Lockington and J.-Y.
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 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 informationPoints to Learn. τ C. τ R. conservative substance reactive substance. τ W V Q. out. out
Pots to Learn Steady State: defition Assimilation capacity - Concept - Mathematical defition: Transfer Coefficient, β =C /C Residence time: hydraulic conservative substance reactive substance τ R = a =
More informationpt)x+ x control volume
Appendix I DERIVATION OF EQUATIONS OF STEADY-STATE GROUNDWATER FLOW Consider a unit volume of saturated porous media (Figure AI. 1). In fluid mechanics, such a volume is called a control volume. The boundaries
More informationGEOL 474/674 Practice Exam #1 Fall This exam counts 20% of your grade for this course; your instructions are as follows:
Write your name here: This exam counts 20% of your grade for this course; your instructions are as follows: 1) You have 75 minutes to finish this exam. 2) No books, notes, or discussion is allowed. If
More information6. GRID DESIGN AND ACCURACY IN NUMERICAL SIMULATIONS OF VARIABLY SATURATED FLOW IN RANDOM MEDIA: REVIEW AND NUMERICAL ANALYSIS
Harter Dissertation - 1994-132 6. GRID DESIGN AND ACCURACY IN NUMERICAL SIMULATIONS OF VARIABLY SATURATED FLOW IN RANDOM MEDIA: REVIEW AND NUMERICAL ANALYSIS 6.1 Introduction Most numerical stochastic
More informationPollution. Elixir Pollution 97 (2016)
42253 Available online at www.elixirpublishers.com (Elixir International Journal) Pollution Elixir Pollution 97 (2016) 42253-42257 Analytical Solution of Temporally Dispersion of Solute through Semi- Infinite
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 informationMETHODOLOGY FOR DRENABILITY STUDIES OF STACKED MINING COARSE TAILINGS
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
More informationON THE CONSERVATION OF MASS AND ENERGY IN HYGROTHERMAL NUMERICAL SIMULATION WITH COMSOL MULTIPHYSICS
ON THE CONSERVATION OF MASS AND ENERGY IN HYGROTHERMAL NUMERICAL SIMULATION WITH COMSOL MULTIPHYSICS Michele Bianchi Janetti, Fabian Ochs, and Wolfgang Feist,2 Unit for Energy Efficient Buildings, University
More informationIncremental identification of transport phenomena in wavy films
17 th European Symposium on Computer Aided Process Engineering ESCAPE17 V. Plesu and P.S. Agachi (Editors) 2007 Elsevier B.V. All rights reserved. 1 Incremental identification of transport phenomena in
More informationThe matric flux potential: history and root water uptake
The matric flux potential: history and root water uptake Marius Heinen, Peter de Willigen (Alterra, Wageningen-UR) Jos van Dam, Klaas Metselaar (Soil Physics, Ecohydrology and Groundwater Management Group,
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 informationThe most common methods to identify velocity of flow are pathlines, streaklines and streamlines.
4 FLUID FLOW 4.1 Introduction Many civil engineering problems in fluid mechanics are concerned with fluids in motion. The distribution of potable water, the collection of domestic sewage and storm water,
More informationWMS 9.0 Tutorial GSSHA Modeling Basics Infiltration Learn how to add infiltration to your GSSHA model
v. 9.0 WMS 9.0 Tutorial GSSHA Modeling Basics Infiltration Learn how to add infiltration to your GSSHA model Objectives This workshop builds on the model developed in the previous workshop and shows you
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 informationOn variational data assimilation for 1D and 2D fluvial hydraulics
On variational data assimilation for D and D fluvial hydraulics I. Gejadze, M. Honnorat, FX Le Dimet, and J. Monnier LJK - MOISE project-team, Grenoble, France. Contact: Jerome.Monnier@imag.fr Civil Engineering
More informationResearch Scholar, PACIFIC University, Udaipur, Rajasthan, India. Department of Mathematics, UGC Visiting Professor, HNGU Patan, Gujarat, India
2018 IJSRSET Volume 4 Issue 4 Print ISSN: 2395-1990 Online ISSN : 2394-4099 Themed Section : Engineering and Technology Adomian Decomposition Study of Unidimensional Flow in Unsaturated Porous Media Khatri
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 information