Homogenization and numerical Upscaling. Unsaturated flow and two-phase flow

Size: px
Start display at page:

Download "Homogenization and numerical Upscaling. Unsaturated flow and two-phase flow"

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)

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 information

1. Water in Soils: Infiltration and Redistribution

1. 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 information

Multiscale Computation for Incompressible Flow and Transport Problems

Multiscale 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 information

Homogenization Theory

Homogenization 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 information

In all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.

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 information

Introduction to Aspects of Multiscale Modeling as Applied to Porous Media

Introduction 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 information

Simulation of Unsaturated Flow Using Richards Equation

Simulation 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 information

RATE OF FLUID FLOW THROUGH POROUS MEDIA

RATE 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 information

Multiscale Finite Element Methods. Theory and

Multiscale 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 information

Temperature dependent multiphase flow and transport

Temperature 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 information

A GENERALIZED CONVECTION-DIFFUSION MODEL FOR SUBGRID TRANSPORT IN POROUS MEDIA

A 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 information

Modelling 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 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.

ψ 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 information

Darcy s Law, Richards Equation, and Green-Ampt Equation

Darcy 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

Water in Soil Sections in Craig

Water 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 information

GG655/CEE623 Groundwater Modeling. Aly I. El-Kadi

GG655/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 information

Seminar zu aktuellen Themen der Numerik im Wintersemester 2010/2011

Seminar 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 information

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 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 information

Analysis of oil displacement by water in oil reservoirs with horizontal wells

Analysis 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 information

Comparison of conductivity averaging methods for one-dimensional unsaturated flow in layered soils

Comparison 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 information

Simulation of Imbibition Phenomena in Fluid Flow through Fractured Heterogeneous Porous Media with Different Porous Materials

Simulation 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 information

dynamics of f luids in porous media

dynamics 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 information

Treatment of tensorial relative permeabilities with multipoint flux approximation

Treatment 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 information

Upscaling of Richards equation for soil moisture dynamics to be utilized in mesoscale atmospheric models

Upscaling 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 information

Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs

Numerical 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 information

Modeling of 1D Anomalous Diffusion In Fractured Nanoporous Media

Modeling 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 information

16 Rainfall on a Slope

16 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 information

Modeling moisture transport by periodic homogenization in unsaturated porous media

Modeling 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 information

A 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 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 information

Stochastic geometry and porous media

Stochastic 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

Flow in unsaturated random porous media, nonlinear numerical analysis and comparison to analytical stochastic models

Flow 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 information

APPENDIX 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 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 information

Darcy's Law. Laboratory 2 HWR 531/431

Darcy'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 information

TREATMENT OF TENSORIAL RELATIVE PERMEABILITIES WITH MULTIPOINT FLUX APPROXIMATION

TREATMENT 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 information

An 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 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 information

Hydraulic properties of porous media

Hydraulic 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 information

Analytical approach predicting water bidirectional transfers: application to micro and furrow irrigation

Analytical 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 information

Evaporation-driven soil salinization

Evaporation-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 information

Analysis of Multiphase Flow under the Ground Water

Analysis 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 information

MULTISCALE FINITE ELEMENT METHODS FOR NONLINEAR PROBLEMS AND THEIR APPLICATIONS

MULTISCALE 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 information

MULTISCALE FINITE ELEMENT METHODS FOR NONLINEAR PROBLEMS AND THEIR APPLICATIONS

MULTISCALE 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 information

Quantifying shallow subsurface flow and salt transport in the Canadian Prairies

Quantifying 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 information

Soils, Hydrogeology, and Aquifer Properties. Philip B. Bedient 2006 Rice University

Soils, 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 information

II. Dealing with Physical Laws: A Stochastic Flowpath Analysis Example

II. 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 information

DNAPL migration through interbedded clay-sand sequences

DNAPL 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 information

TENSOR RELATIVE PERMEABILITIES: ORIGINS, MODELING AND NUMERICAL DISCRETIZATION

TENSOR 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 information

Chapter Seven. For ideal gases, the ideal gas law provides a precise relationship between density and pressure:

Chapter 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 information

1 FLUID-MECHANICAL INTERACTION SINGLE FLUID PHASE

1 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 information

Xingye Yue. Soochow University, Suzhou, China.

Xingye 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 information

MULTISCALE 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 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 information

Advanced 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 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 information

Macro-Scale Dynamic Effects in Homogeneous and Heterogeneous Porous Media

Macro-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 information

Numerical Simulation of Flows in Highly Heterogeneous Porous Media

Numerical 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 information

Field 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 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 information

PRECONDITIONING MARKOV CHAIN MONTE CARLO SIMULATIONS USING COARSE-SCALE MODELS

PRECONDITIONING 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 information

Comparison of Heat and Mass Transport at the Micro-Scale

Comparison 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 information

Sunday 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. 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 information

Multi-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 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 information

On the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method

On 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 information

Earth dam steady state seepage analysis

Earth 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 information

Modeling 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 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 information

Mixed-hybrid finite element method for modelling two-phase flow in porous media

Mixed-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 information

Modeling Reactive Flows in Porous Media

Modeling 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 information

Assessment of Hydraulic Conductivity Upscaling Techniques and. Associated Uncertainty

Assessment 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 information

Surface Processes Focus on Mass Wasting (Chapter 10)

Surface 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 information

Introduction to Aspects of Multiscale Modeling as Applied to Porous Media

Introduction 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 information

Numerical Simulations of Two-Phase Flow in Rigid Porous Media

Numerical 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 information

Two-phase gravity currents in porous media

Two-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 information

Operator Upscaling for the Wave Equation

Operator 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 information

Instructor : Dr. Jehad Hamad. Chapter (7)

Instructor : 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 information

5. 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)

5. 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 information

NUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES

NUMERICAL 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 information

Effective unsaturated hydraulic conductivity for one-dimensional structured heterogeneity

Effective 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 information

C) D) 3. Which graph best represents the relationship between soil particle size and the rate at which water infiltrates permeable soil?

C) 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 information

Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium

Effects 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 information

EXAMPLE PROBLEMS. 1. Example 1 - Column Infiltration

EXAMPLE 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 information

1 Modeling Immiscible Fluid Flow in Porous Media

1 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 information

INTERPOLATION AND UPDATE IN DYNAMIC DATA-DRIVEN APPLICATION SIMULATIONS

INTERPOLATION 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 information

pifreeze A Freeze / Thaw Plug-in for FEFLOW User Guide

pifreeze 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 information

Fracture-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 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 information

pressure-saturations relationships for two-phase flow in 3D porous media: implications micro-heterogeneities

pressure-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 information

Multiscale Methods for Multiphase Flow in Porous Media

Multiscale 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 information

Investigations 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 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 information

Civil Engineering Department College of Engineering

Civil 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 information

A Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach

A 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 information

Two-Scale Wave Equation Modeling for Seismic Inversion

Two-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 information

R.M. Manasipov and J.E. Mindel. Montanuniversitaet Leoben. October 3, 2014

R.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 information

Modeling two-phase flow in strongly heterogeneous porous media

Modeling 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 information

A 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 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 information

Numerical investigations of hillslopes with variably saturated subsurface and overland flows

Numerical 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 information

AN UNFITTED METHOD FOR TWO-PHASE FLOW IN FRACTURED POROUS MEDIA

AN 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 information

Available online at ScienceDirect

Available 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 information

CO 2 storage capacity and injectivity analysis through the integrated reservoir modelling

CO 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 information

Analytical solutions for water flow and solute transport in the unsaturated zone

Analytical 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 information

Natural Convection in Vertical Channels with Porous Media and Adiabatic Extensions

Natural 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 information

Seepage. c ZACE Services Ltd. August 2011

Seepage. 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 information

Onset of convection of a reacting fluid layer in a porous medium with temperature-dependent heat source

Onset 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 information

Evaporation-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 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 information

2. Governing Equations. 1. Introduction

2. 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 information

Numerical Solution of the Two-Dimensional Time-Dependent Transport Equation. Khaled Ismail Hamza 1 EXTENDED ABSTRACT

Numerical 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