Homogenization and numerical Upscaling. Unsaturated flow and two-phase flow
|
|
- Bertram Taylor
- 6 years ago
- Views:
Transcription
1 Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart
2 Outline Block 1: Introduction and Repetition Homogenization of the Richards equation Multiscale methods for the Richards equation Block 2: Homogenization of the one-dimensional two-phase flow equation Multiscale methods for the two-phase flow problem
3 Literature Homogenization of the Richards equation: J. Lewandowska and J.-P. Laurent. Moisture transfer in an unsaturated heterogeneous porous medium. Transport in Porous Media, 45: , J. Lewandowska, A. Szymkiewicz, K. Burzynski, and M. Vauclin. Modeling of unsaturated water flow in double-porosity soils by the homogenization approach. Advances in Water Resources, 27: , I. Neuweiler and O. Cirpka. Homogenization of richards equation in permeability fields with different connectivities. Water Resources Research, 41(2): doi: /2004wr003329, I. Neuweiler and H. Eichel. Effective parameter functions for richards equation in layered porous media. Vadose Zone Journal, accepted for publication, 2006.
4 Literature Numerical Homogenization of the Richards equation: Y. Efendiev, T. Hou and V. Ginting: Multiscale finite element methods for nonlinear problems and their application, Comm. Math. Sci. 2(4), , Homogenization of viscous dominated flow: R. Mauri. Heat and mass transport in nonhomogeneous random velocity fields. Physical Review E, 68:doi: /PhysRevE , C. J. VanDuijn, A. Mikelic, and I. S. Pop. Effective equations for twophase flow with trapping on the micro scale. SIAM Journal of Applied Mathematics, 62(5): , 2002.
5 Literature Multiscale methods for viscous dominated two-phase flow: P. Jenny, S.H. Lee and H. A. Tchelepi: Adaptive multiscale finite/volume method for multiphase flow and transport in porous media, Multiscale Model. Simul. 3(1), 50-64, P. Jenny, S.H. Lee and H. A. Tchelepi: Multi-scale finite-volume method for elliptic problems in subsurface flow simulation, Journal of Computational Physics 187, 47-67, I. Lunati and P. Jenny, Multiscale finite-volume method for compressible multiphase flow in porous media, Journal of Computational Physics 216, , L. J. Durlofsky, Numerical calculation of equivalent grid block permeability tensors for heterogeneous porous media, Water Resources Research 27(5), , Y. Efendiev and L. J. Durlofsky: Numerical modeling of subgrid heterogeneity in two-phase flow simulations, Water Resources Research 38(8), doi: /2000WR000190, 2002.
6 Block 1 Introduction
7 Introduction and Repetition Homogenization Theory: Starting point: Model with space-dependent parameters Derivation of a model for the spatially averaged parameters Requirements: Space scales in the medium have to be separated Assumptions: Solution has two separated scales, solution is stationary on the small scale Homogenized model is valid for: Fixed typical time scale (usually fast time scales are neglected), fixed dimensionless numbers (problematic for nonlinear problems, where they change), fixed parameter contrasts (also problematic for nonlinear problems) Treated for the elliptic problem
8 Introduction and Repetition Fields with scale separation: Counter example: Length scale L REV Length scale l Unit cell l????
9 Introduction and Repetition Derivation of an upscaled, homogeneous model, where the heterogeneity is no longer seen Consider the limit...
10 Introduction and Repetition Hypothesis: local periodicity l
11 Procedure: Introduction and Repetition Adaption of a dimensionless form of the equation Scaling of the dimensionless numbers and parameter contrasts with ε Choice of a reference system Separating the space variables into small scale and large scale Expansion of the variables in orders of ε Solving for different orders of ε separately Limit ε >0, lowest order equation (space averaged over the unit cell) is the upscaled equation Effective parameters of the model are defined from cell problems Proof of convergence, existence, (-> mathematicians)
12 Introduction and Repetition Credo: Homogenization is a very useful method and improves our understanding of processes in heterogeneous media. It can help to develop models for processes in heterogeneous media. However, we make a lot of inherent assumptions and approximations. As long as we are aware of them: FINE!
13 Block 1 Homogenization of Richards equation
14 Contaminant River Infiltration Root uptake Unsaturated zone: air + water Scale: 10 0 m Aquifer
15 Richards equation Mass balance of the water phase assuming air is always at atmospheric pressure Flux: Darcys law Saturation of water Flux Water pressure head Saturated conductivity Relative permeability
16 Air pressure is constant: Richards equation Capillary pressure head Capillary pressure saturation relation Relative permeability
17 Capillary pressure saturation relation in porous media Unique relation between water saturation and capillary pressure: (Pressure on the macroscale, averaged over the pore space) Entry pressure (capillary entry pressure of the largest pore)
18 Relative permeability
19 Richards equation But also
20 z above groundwater [m] Homogenization: Richards equation Typical staedy state profiles: Groundwaterlevel: Constant head Water-Content Profile in Silt as Function of Recharge [mm/year] Gravity dominated 0 0,1 0,2 0,3 0,4 0,5 θ [-] Capillary-Head Profile in Silt as Function of Recharge [mm/year] Gravity dominated h c [m] Constant flux hydrostatic
21 1. Step: Dimensionless equation Variables with stars are dimensionless Typical values for: Pressure head (needed only for the parameter functions) Conductivity (e.g. geometric mean in 2d) Length scale (there are two, here: large length scale) Time scale Dimensionless number:
22 Choice of the time scale Time scale due to large scale pressure gradients: Time scale due gravity, large scale: Time scale due to small scale pressure gradients: Time scale due gravity, small scale: In principle a scale separation of time would be necessary. This is not done, usually a large time scale is chosen.
23 Problem: estimation of pressure head fluctuations Could be due to large scale boundary conditions Could be due to fluctuations of the solution in a heterogeneous medium Could be due a solution in a saturation regime, where the P c -S curve is steep Difficult to estimate a priori
24 Example: Capillary-Head Profile in Silt as Function of Recharge [mm/year] Due to small scale 14 fluctuations Due to small large scale solution h c [m] hydrostatic
25 Choice of the parameter contrast, e.g in: In a two-material composite, what determines the flow is With typical values: Problem: during a transient flow process this ratio can change a lot, as it depends strongly on the saturation. For the estimation of the dimensionless numbers, time scale and parameter contrast we have in principle to choose the saturation dependent values!! Means: We have to have an idea about the solution already.
26 2. Step: Two-scale Expansion Reference system: Large scale L Separation of scales -> space variable can be made dimensionless with two different length scales -> Two different variables: Transformation of the gradient:
27 Expansion of the variables in terms of ε
28 First case: Homogenization: Richards equation -> moderate heterogeneities E.g. lognormally distributed field -> small variances -> Typical pressure head fluctuations are in the same order of magnitude as the large length scale
29 Example for such a scenario: Slow changes in a scenario such as: 6 Silt with recharge rate: 10 mm / year Material 1 Material 2 z [m] 0 0 0,5 1 1,5 2 2,5 3 3,5 h C [m] Or horizontal flow
30 Expansion of everything in terms of ε The whole equation with these parameters:
31 Solving for different orders of ε separately: Order Boundary conditions: Is -periodic is in average constant on Problem corresponds to gravitational free unsaturated flow with periodic boundary conditions without source term -> only solution is that is constant on the small scale (would otherwise be a perpetuum mobile).
32 This condition is equivalent to capillary equilibrium The typical time scale is too large to appear in this condition -> quasi-instantaneous The pressure head (equal to negative capillary pressure head) is constant over the unit cell -> not influenced by the heterogeneities Implications for the paramter functions: Not constant on the small scale, but instantaneously in equilibrium -> not dynamics
33 Unit cell Heterogeneous capillary entry pressure head: Water Saturation: Unsaturated conductivity:
34 Order Coupling of large and small scale -> decoupling via
35 Cell problem: Boundaries: is periodic. As in single phase flow, as depends only explicitely on,but not indirectly, as is constant in -> Problem boils down to the single phase flow problem with a fixed value of
36 Upscaled problem: Order Same problem as the Richards equation with effective retention function and effective unsaturated conductivity function
37 For an isotropic field:
38 3. Step: Derivation of the effective retention function Choose a value for the pressure head h (0) Calculate the saturation field at this pressure head according to the local retention curves Average the saturation over the unit cell This yields a point on the effective retention curve Repeat the procedure for more values h (0).
39 Exercise 1: Example: 2d Checkerboardfield Porosity: Unit cell Material 1: Material 2: Calculate the values of the effective retention curve at
40 Solution: Θ
41 4. Step: Derivation of the effective unsaturated conductivity Cell problem In each component:
42 That means: Boundaries: is periodic. This is equivalent to a singlephase pressure problem on the unit cell with a unit pressure gradient in both directions. If the conductivity is isotropic, it is sufficient to solve one problem, the effective unsaturated conductivity is the resulting averaged flux (average pressure gradient is unit). As this is equivalent to single phase flow: Can be estimated with the methods we know for this problem.
43 Exercise 2: Example: 2d Checkerboardfield Porosity: Unit cell Material 1: Material 2: Calculate the values of the effective conductivity curve at Use the geometric mean as approximation for the effective conductivity at one value
44 Solution:
45
46 Analytic results for log-normally distributed K field and Brooks Corey model assuming Gaussian fields Effective retention curve: (no approximation) From: Neuweiler and Cirpka, WRR, 2005
47 Analytic results for log-normally distributed K field and Brooks Corey model assuming Gaussian fields Effective permeability curve: (second order in variance of total permeability) From: Neuweiler and Cirpka, WRR, 2005
48 Second case: Strong parameter contrast in a two-material composite (Lewandowska et al., AWR, 2004) -> Typical pressure head fluctuations are in the same order of magnitude as the large length scale
49 In material 1: In material 2:
50 Variable of the upscaled system: Was: Is now: Consider only: Small capacity Large capacity Coupled to the inclusion via the flux boundaries
51 Order in material 1: -> Capillary equilibrium in material 1 Order in material 2: -> dynamics on the small scale
52 Order in material 1:. As in the Richards problem before, only this time we consider only material 1. The inclusion sphere of material two is treated as impermeable. -> Effective unsaturated conductivity from the cell problem:
53 Order in material 1: Before: Averaging over the unit cell cancelled the derivative with respect to the small scale. Now: Leads to coupling to the inclusion
54 Leads to the averaged equation of order in material 1: Effective parameters as in the usual Richards problem with the cell problem Dynamic effect: Requires a parametrization.
55 Physical interpretation: Typical time scale of gravity on the large scale is equal to the typical time scale of capillary forces on the small scale. -> Appears as an additional effect
56 Third case: Gravity dominated flow (Neuweiler and Eichel, Vadose Zone Journal, in press)
57 Order : No longer capillary equilibrium, equation has to be solved numerically.
58 Example vertical flow in a layered medium: Local head distribution
59 Upscaled problem in leading order: -> Effective curves can be approximated reasonably well with the capillary equilibrium approach
60 Block 1 Numerical homogenization of Richards equation
61 Numerical Homogenization of Richards equation Aim: Improvement of fine scale FE methods Idea: making use of two scale expansions Richards Equ. -in variational formulation - on the coarse mesh
62 Numerical Homogenization of Richards equation Macro/micro approach: use a coarse mesh for defining the nodal values of φ( x j ) and a fine mesh for computing the basis functions ϕ j ( x). The problem dimension is that of the coarse mesh Multiscale FEM are defined for non-periodic problems but their quantitative convergence analysis is made in the periodic case. Main references for this example: Hou, Efendiev, Wu, Babuska, Matache, Schwab, Brizzi...
63 Numerical Homogenization of Richards equation ε [ K u ( φ )( φ + e )] ε = 0 i i z The construction of the basis functions neglects gravity effects!
64 Comparison of fine scale solution with MsFEM: Numerical Homogenization of Richards equation We follow infiltration fronts moving from top to bottom in an a)isotropic medium b)anisotroic medium
Unsaturated 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 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 informationMultiscale Computation for Incompressible Flow and Transport Problems
Multiscale Computation for Incompressible Flow and Transport Problems Thomas Y. Hou Applied Mathematics, Caltech Collaborators: Y. Efenidev (TAMU), V. Ginting (Colorado), T. Strinopolous (Caltech), Danping
More informationHomogenization Theory
Homogenization Theory Sabine Attinger Lecture: Homogenization Tuesday Wednesday Thursday August 15 August 16 August 17 Lecture Block 1 Motivation Basic Ideas Elliptic Equations Calculation of Effective
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 informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part IV Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
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 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 informationMultiscale Finite Element Methods. Theory and
Yalchin Efendiev and Thomas Y. Hou Multiscale Finite Element Methods. Theory and Applications. Multiscale Finite Element Methods. Theory and Applications. October 14, 2008 Springer Dedicated to my parents,
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 informationA GENERALIZED CONVECTION-DIFFUSION MODEL FOR SUBGRID TRANSPORT IN POROUS MEDIA
MULTISCALE MODEL. SIMUL. Vol. 1, No. 3, pp. 504 526 c 2003 Society for Industrial and Applied Mathematics A GENERALIZED CONVECTION-DIFFUSION MODEL FOR SUBGRID TRANSPORT IN POROUS MEDIA Y. EFENDIEV AND
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 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 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 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 informationGG655/CEE623 Groundwater Modeling. Aly I. El-Kadi
GG655/CEE63 Groundwater Modeling Model Theory Water Flow Aly I. El-Kadi Hydrogeology 1 Saline water in oceans = 97.% Ice caps and glaciers =.14% Groundwater = 0.61% Surface water = 0.009% Soil moisture
More informationSeminar zu aktuellen Themen der Numerik im Wintersemester 2010/2011
Seminar zu aktuellen Themen der Numerik im Wintersemester 2010/2011 Modeling Two-Phase-Two-Component Processes in Porous Media Using an Upscaling-Multi-Scale Method Elin Solberg 27 January 2011 Contents
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 informationAnalysis of oil displacement by water in oil reservoirs with horizontal wells
Analysis of oil displacement by water in oil reservoirs with horizontal wells Paulo Dore Fernandes, Thiago Judson L. de Oliveira and Rodrigo A. C. Dias Problem Description This work involves near-well
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 informationSimulation of Imbibition Phenomena in Fluid Flow through Fractured Heterogeneous Porous Media with Different Porous Materials
Journal of Applied Fluid Mechanics, Vol. 10, No. 5, pp. 1451-1460, 2017. Available online at.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. DOI: 10.169/acadpub.jafm.73.242.2721 Simulation of Imbibition
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 informationTreatment of tensorial relative permeabilities with multipoint flux approximation
M. Wolff a B. Flemisch a R. Helmig a I. Aavatsmark b Treatment of tensorial relative permeabilities with multipoint flux approximation Stuttgart, March 2011 a Department of Hydromechanics and Modeling
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 informationNumerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs
Transport in Porous Media 49: 209 240, 2002. 2002 Kluwer Academic Publishers. Printed in the Netherlands. 209 Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured
More informationModeling of 1D Anomalous Diffusion In Fractured Nanoporous Media
LowPerm2015 Colorado School of Mines Low Permeability Media and Nanoporous Materials from Characterisation to Modelling: Can We Do It Better? IFPEN / Rueil-Malmaison - 9-11 June 2015 CSM Modeling of 1D
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 informationModeling moisture transport by periodic homogenization in unsaturated porous media
Vol. 2, 4. 377-384 (2011) Revue de Mécanique Appliquée et Théorique Modeling moisture transport by periodic homogenization in unsaturated porous media W. Mchirgui LEPTIAB, Université de La Rochelle, France,
More informationA MULTISCALE METHOD FOR MODELING TRANSPORT IN POROUS MEDIA ON UNSTRUCTURED CORNER-POINT GRIDS
A MULTISCALE METHOD FOR MODELING TRANSPORT IN POROUS MEDIA ON UNSTRUCTURED CORNER-POINT GRIDS JØRG E. AARNES AND YALCHIN EFENDIEV Abstract. methods are currently under active investigation for the simulation
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 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 informationAPPENDIX Tidally induced groundwater circulation in an unconfined coastal aquifer modeled with a Hele-Shaw cell
APPENDIX Tidally induced groundwater circulation in an unconfined coastal aquifer modeled with a Hele-Shaw cell AaronJ.Mango* Mark W. Schmeeckle* David Jon Furbish* Department of Geological Sciences, Florida
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 informationTREATMENT OF TENSORIAL RELATIVE PERMEABILITIES WITH MULTIPOINT FLUX APPROXIMATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number 3, Pages 725 744 c 2012 Institute for Scientific Computing and Information TREATMENT OF TENSORIAL RELATIVE PERMEABILITIES WITH
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 informationHydraulic properties of porous media
PART 5 Hydraulic properties of porous media Porosity Definition: Void space: n V void /V total total porosity e V void /V solid Primary porosity - between grains Secondary porosity - fracture or solution
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 informationEvaporation-driven soil salinization
Evaporation-driven soil salinization Vishal Jambhekar 1 Karen Schmid 1, Rainer Helmig 1,Sorin Pop 2 and Nima Shokri 3 1 Department of Hydromechanics and Hydrosystemmodeling, University of Stuttgart 2 Department
More informationAnalysis of Multiphase Flow under the Ground Water
Analysis of Multiphase Flow under the Ground Water Pramod Kumar Pant Department of Mathematics, Bhagwant University, Ajmer, Rajasthan, India Abstract The single-phase fluid flow through a porous medium
More informationMULTISCALE FINITE ELEMENT METHODS FOR NONLINEAR PROBLEMS AND THEIR APPLICATIONS
MULTISCALE FINITE ELEMENT METHODS FOR NONLINEAR PROBLEMS AND THEIR APPLICATIONS Y. EFENDIEV, T. HOU, AND V. GINTING Abstract. In this paper we propose a generalization of multiscale finite element methods
More informationMULTISCALE FINITE ELEMENT METHODS FOR NONLINEAR PROBLEMS AND THEIR APPLICATIONS
COMM. MATH. SCI. Vol. 2, No. 4, pp. 553 589 c 2004 International Press MULTISCALE FINITE ELEMENT METHODS FOR NONLINEAR PROBLEMS AND THEIR APPLICATIONS Y. EFENDIEV, T. HOU, AND V. GINTING Abstract. In this
More informationQuantifying shallow subsurface flow and salt transport in the Canadian Prairies
Quantifying shallow subsurface flow and salt transport in the Canadian Prairies Andrew Ireson GIWS, University of Saskatchewan www.usask.ca/water Uri Nachshon Garth van der Kamp GIWS, University of Saskatchewan
More informationSoils, Hydrogeology, and Aquifer Properties. Philip B. Bedient 2006 Rice University
Soils, Hydrogeology, and Aquifer Properties Philip B. Bedient 2006 Rice University Charbeneau, 2000. Basin Hydrologic Cycle Global Water Supply Distribution 3% of earth s water is fresh - 97% oceans 1%
More informationII. Dealing with Physical Laws: A Stochastic Flowpath Analysis Example
II. Dealing with Physical Laws: A Stochastic Flowpath Analysis Example In this chapter we will examine a physical problem where the significance of observing the underlying principles will be displayed.
More informationDNAPL migration through interbedded clay-sand sequences
Groundwater Quality: Natural and Enhanced Restoration of Groundwater Pollution (Proceedings ofthe Groundwater Quality 2001 Conference held al Sheffield. UK. June 2001). IAHS Publ. no. 275. 2002. 455 DNAPL
More informationTENSOR RELATIVE PERMEABILITIES: ORIGINS, MODELING AND NUMERICAL DISCRETIZATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number 3, Pages 701 724 c 2012 Institute for Scientific Computing and Information TENSOR RELATIVE PERMEABILITIES: ORIGINS, MODELING AND
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 information1 FLUID-MECHANICAL INTERACTION SINGLE FLUID PHASE
FLUID-MECHANICAL INTERACTION SINGLE FLUID PHASE 1-1 1 FLUID-MECHANICAL INTERACTION SINGLE FLUID PHASE 1.1 Introduction FLAC models the flow of fluid (e.g., groundwater) through a permeable solid, such
More informationXingye Yue. Soochow University, Suzhou, China.
Relations between the multiscale methods p. 1/30 Relationships beteween multiscale methods for elliptic homogenization problems Xingye Yue Soochow University, Suzhou, China xyyue@suda.edu.cn Relations
More informationMULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION
MULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION P. DOSTERT, Y. EFENDIEV, AND T.Y. HOU Abstract. In this paper, we study multiscale
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 informationMacro-Scale Dynamic Effects in Homogeneous and Heterogeneous Porous Media
Transp Porous Med (2005) 58:121 145 Springer 2005 DOI 10.1007/s11242-004-5472-6 Macro-Scale Dynamic Effects in Homogeneous and Heterogeneous Porous Media SABINE MANTHEY 1,, S. MAJID HASSANIZADEH 2 and
More informationNumerical Simulation of Flows in Highly Heterogeneous Porous Media
Numerical Simulation of Flows in Highly Heterogeneous Porous Media R. Lazarov, Y. Efendiev, J. Galvis, K. Shi, J. Willems The Second International Conference on Engineering and Computational Mathematics
More informationField Scale Modeling of Local Capillary Trapping during CO 2 Injection into the Saline Aquifer. Bo Ren, Larry Lake, Steven Bryant
Field Scale Modeling of Local Capillary Trapping during CO 2 Injection into the Saline Aquifer Bo Ren, Larry Lake, Steven Bryant 2 nd Biennial CO 2 for EOR as CCUS Conference Houston, TX October 4-6, 2015
More informationPRECONDITIONING MARKOV CHAIN MONTE CARLO SIMULATIONS USING COARSE-SCALE MODELS
PRECONDITIONING MARKOV CHAIN MONTE CARLO SIMULATIONS USING COARSE-SCALE MODELS Y. EFENDIEV, T. HOU, AND W. LUO Abstract. We study the preconditioning of Markov Chain Monte Carlo (MCMC) methods using coarse-scale
More informationComparison of Heat and Mass Transport at the Micro-Scale
Comparison of Heat and Mass Transport at the Micro-Scale E. Holzbecher, S. Oehlmann Georg-August Univ. Göttingen *Goldschmidtstr. 3, 37077 Göttingen, GERMANY, eholzbe@gwdg.de Abstract: Phenomena of heat
More informationSunday September 28th. Time. 06:00 pm 09:00 pm Registration (CAAS) IP: Invited Presentation (IS or SP) 55 mn. CP: Contributed Presentation 25 mn
Sunday September 28th 06:00 pm 09:00 pm Registration (CAAS) IP: Invited Presentation (IS or SP) 55 mn CP: Contributed Presentation 25 mn Monday September 29th 07:30 am 08:30 am Registration 08:30 am 09:00
More informationMulti-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media
Multi-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media Jan Tecklenburg a,, Insa Neuweiler a, Jesus Carrera b, Marco Dentz b a Institute of Fluid Mechanics
More informationOn the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2006; 22:741 751 Published online 13 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cnm.846
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 informationModeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids
Modeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids Jørg Espen Aarnes and Vera Louise Hauge SINTEF ICT, Deptartment of Applied Mathematics Applied Mathematics
More informationMixed-hybrid finite element method for modelling two-phase flow in porous media
Journal of Math-for-Industry, Vol. 3 (211C-2), pp. 9 19 Mixed-hybrid finite element method for modelling two-phase flow in porous media Radek Fučík and Jiří Mikyška Revised on September 23, 211 Abstract.
More informationModeling Reactive Flows in Porous Media
Modeling Reactive Flows in Porous Media Peter Lichtner (lead PI), Los Alamos National Laboratory Glenn Hammond, Pacific Northwest National Laboratory Richard Tran Mills, Oak Ridge National Laboratory NCCS
More informationAssessment of Hydraulic Conductivity Upscaling Techniques and. Associated Uncertainty
CMWRXVI Assessment of Hydraulic Conductivity Upscaling Techniques and Associated Uncertainty FARAG BOTROS,, 4, AHMED HASSAN 3, 4, AND GREG POHLL Division of Hydrologic Sciences, University of Nevada, Reno
More informationSurface Processes Focus on Mass Wasting (Chapter 10)
Surface Processes Focus on Mass Wasting (Chapter 10) 1. What is the distinction between weathering, mass wasting, and erosion? 2. What is the controlling force in mass wasting? What force provides resistance?
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part II Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
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 informationTwo-phase gravity currents in porous media
J. Fluid Mech. (211), vol. 678, pp. 248 27. c Cambridge University Press 211 doi:1.117/jfm.211.11 Two-phase gravity currents in porous media MADELEINE J. GOLDING, JEROME A. NEUFELD, MARC A. HESSE AND HERBERT
More informationOperator Upscaling for the Wave Equation
Operator Upscaling for the Wave Equation Tetyana Vdovina Susan E. Minkoff UMBC), Oksana Korostyshevskaya Department of Computational and Applied Mathematics Rice University, Houston TX vdovina@caam.rice.edu
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 information5. Which surface soil type has the slowest permeability rate and is most likely to produce flooding? A) pebbles B) sand C) silt D) clay A) B) C) D)
1. During a heavy rainstorm, soil samples A and B both became saturated with water. However, 10 minutes after the storm ended, the soils appeared as shown below. Which statement best explains the observed
More informationNUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 NUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES Steven
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 informationC) D) 3. Which graph best represents the relationship between soil particle size and the rate at which water infiltrates permeable soil?
1. Which earth material covering the surface of a landfill would permit the least amount of rainwater to infiltrate the surface? A) silt B) clay C) sand D) pebbles 2. Which graph best represents the relationship
More informationEffects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
Transport in Porous Media (2006) 64: 1 14 Springer 2006 DOI 10.1007/s11242-005-1126-6 Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
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 information1 Modeling Immiscible Fluid Flow in Porous Media
Excerpts from the Habilitation Thesis of Peter Bastian. For references as well as the full text, see http://cox.iwr.uni-heidelberg.de/people/peter/pdf/bastian_habilitationthesis.pdf. Used with permission.
More informationINTERPOLATION AND UPDATE IN DYNAMIC DATA-DRIVEN APPLICATION SIMULATIONS
INTERPOLATION AND UPDATE IN DYNAMIC DATA-DRIVEN APPLICATION SIMULATIONS CRAIG C. DOUGLAS University of Kentucky, Department of Computer Science, 325 McVey Hall, Lexington, KY 4561 and Yale University,
More informationpifreeze A Freeze / Thaw Plug-in for FEFLOW User Guide
pifreeze A Freeze / Thaw Plug-in for FEFLOW User Guide MIKE 2016 DHI headquarters Agern Allé 5 DK-2970 Hørsholm Denmark +45 4516 9200 Telephone +45 4516 9333 Support +45 4516 9292 Telefax mike@dhigroup.com
More informationFracture-Matrix Flow Partitioning and Cross Flow: Numerical Modeling of Laboratory Fractured Core Flood
Fracture-Matrix Flow Partitioning and Cross Flow: Numerical Modeling of Laboratory Fractured Core Flood R. Sanaee *, G. F. Oluyemi, M. Hossain, and M. B. Oyeneyin Robert Gordon University *Corresponding
More informationpressure-saturations relationships for two-phase flow in 3D porous media: implications micro-heterogeneities
Loughborough University Institutional Repository Dynamic effects in capillary pressure-saturations relationships for two-phase flow in 3D porous media: implications of micro-heterogeneities This item was
More informationMultiscale Methods for Multiphase Flow in Porous Media
Multiscale Methods for Multiphase Flow in Porous Media Jan M. Nordbotten Department of Mathematics, University of Bergenjan.nordbotten@math.uib.no Department of Civil and Environmental Engineering, Princeton
More informationInvestigations in Geologic Carbon Sequestration: Multiphase Flow of CO2 and Water in Reservoir Rocks. Annual Report 2015
Investigations in Geologic Carbon Sequestration: Multiphase Flow of CO2 and Water in Reservoir Rocks Annual Report 2015 Sally M. Benson, David Cameron, Ferdinand Hingerl, Andrew Gyenis, Boxiao Li, Christin
More informationCivil Engineering Department College of Engineering
Civil Engineering Department College of Engineering Course: Soil Mechanics (CE 359) Lecturer: Dr. Frederick Owusu-Nimo FREQUENCY CE 260 Results (2013) 30 25 23 25 26 27 21 20 18 15 14 15 Civil Geological
More informationA Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach
SPE-172990-MS A Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach Benjamin Ganis, Kundan Kumar, and Gergina Pencheva, SPE; Mary F. Wheeler, The
More informationTwo-Scale Wave Equation Modeling for Seismic Inversion
Two-Scale Wave Equation Modeling for Seismic Inversion Susan E. Minkoff Department of Mathematics and Statistics University of Maryland Baltimore County Baltimore, MD 21250, USA RICAM Workshop 3: Wave
More informationR.M. Manasipov and J.E. Mindel. Montanuniversitaet Leoben. October 3, 2014
Comparative Analysis of Compositional Two-Phase Flow Modeling in Heterogeneous Media Between the Discrete Event Simulation Method Coupled to a Split Node Formulation and Classical Timestep-Driven Approaches
More informationModeling two-phase flow in strongly heterogeneous porous media
Presented at the COMSOL Conference 2010 China COMSOL 2010 I«^rc Modeling two-phase flow in strongly heterogeneous porous media Zhaoqin Huang Research Center for Oil & Gas Flow in Reservoir, School of Petroleum
More informationA Numerical Simulation Framework for the Design, Management and Optimization of CO 2 Sequestration in Subsurface Formations
A Numerical Simulation Framework for the Design, Management and Optimization of CO 2 Sequestration in Subsurface Formations Investigators Hamdi Tchelepi, Associate Professor of Petroleum Engineering; Lou
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 informationAN UNFITTED METHOD FOR TWO-PHASE FLOW IN FRACTURED POROUS MEDIA
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 AN UNFIED MEHOD FOR WO-PHASE FLOW IN FRACURED POROUS MEDIA A. Fumagalli and A. Scotti
More informationAvailable online at ScienceDirect
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 114 (17) 3343 3352 www.elsevier.com/locate/procedia 13th International Conference on Greenhouse Gas Control Technologies, GHGT-13,
More informationCO 2 storage capacity and injectivity analysis through the integrated reservoir modelling
CO 2 storage capacity and injectivity analysis through the integrated reservoir modelling Dr. Liuqi Wang Geoscience Australia CO 2 Geological Storage and Technology Training School of CAGS Beijing, P.
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 informationNatural Convection in Vertical Channels with Porous Media and Adiabatic Extensions
Natural Convection in Vertical Channels with Porous Media and Adiabatic Extensions Assunta Andreozzi 1,a, Bernardo Buonomo 2,b, Oronzio Manca 2,c and Sergio Nardini 2,d 1 DETEC, Università degli Studi
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 informationOnset of convection of a reacting fluid layer in a porous medium with temperature-dependent heat source
AMERICAN JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH 2011, Science Huβ, http://www.scihub.org/ajsir ISSN: 2153-649X doi:10.5251/ajsir.2011.2.6.860.864 Onset of convection of a reacting fluid layer in
More informationEvaporation-driven transport and precipitation of salt in porous media: A multi-domain approach
Evaporation-driven transport and precipitation of salt in porous media: A multi-domain approach Vishal Jambhekar Karen Schmid, Rainer Helmig Department of Hydromechanics and Hydrosystemmodeling EGU General
More information2. Governing Equations. 1. Introduction
Multiphysics Between Deep Geothermal Water Cycle, Surface Heat Exchanger Cycle and Geothermal Power Plant Cycle Li Wah Wong *,1, Guido Blöcher 1, Oliver Kastner 1, Günter Zimmermann 1 1 International Centre
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 information