Numerical fractional flow modelling of inhomogeneous
|
|
- Bryan Webster
- 5 years ago
- Views:
Transcription
1 Numerical fractional flow modelling of inhomogeneous air sparging E.F. Kaasschieter\ G.J. Mulder^ & J.D. van der Werff ten Bosch^ Department of Mathematics and Computing Science, hoven University of Technology, P.O. Box 513, 5600 MB hoven, The Netherlands, lands, Resource Analysis, Zuiderstraat 110, 2611 SJ Netherlands, jasper, v. nl Abstract Spillage of organic compounds into the subsurface environment can result in costly remediation. As a possibly effective remediation technique, the injection of air into the groundwater (air sparging) has gained attention. The injected air migrates towards the unsaturated zone, volatising contaminants from the groundwater and delivering oxygen to biota for biodegradation (see The two-phase flow of air and water turns out to be predominantly driven by convection. It is modelled using a fractional flow approach, which yields a hyperbolic governing equation for the air saturation. The Godunov method is evaluated for the numerical solution of the one-dimensional problem in an inhomogeneous porous medium. 1 Governing equations Air sparging is a three-dimensional two-phase flow process of air and water, starting in the saturated zone of a porous medium and passing to the unsaturated zone.
2 52 Computer Methods in Water Resources XII It is assumed that the flow of both air and water can be adequately described by Darcy's law. The porous medium is assumed to be isotropic. With respect to the flow, both phases are assumed to have constant and uniform viscosity, and to be incompressible. If the phases are immiscible, then the following equations are representative: OQ <^(x)-^- + V.vi = 0, l = a,w, (1) v; = -- -V(p, + m#4, f = o, w. (2) A Equation (1) expresses conservation of volume. Herein, </> [ ] denotes the porosity while Si [ ] and Vf [m s~*] respectively denote effective saturation and specific discharge. The specific discharge is given by Darcy's law, i.e. (2), in which p\ [Pa], p\ [kg m~^], /// [Pa s], k\ [ } are the fluid pressure, fluid density, viscosity and relative permeability. Clearly, the index / can either be a (air phase) or w (water phase). Also in (2), we have y [m s~^] denoting the acceleration due to gravity and K [m^] the absolute permeability. The two phases are linked by the identity & +,^ = l. (3) The following empirical relations are currently widely used to describe how capillary pressure PC [Pa] and relative permeability in a two-phase system depend on saturation (see [2]): a(x) (4) where in [-] and a [m *] are empirical soil parameters satisfying 0 < m < I and a > 0. Not included in (4) is the phenomenon of hysteresis. The system (l)-(4) is referred to as the mixed form of Richards equation (see [1]). 2 The fractional flow model A formulation that is often used for multi-phase flow problems is the so-called fractional flow model. It is equivalent to (l)-(4) and yields
3 Computer Methods in Water Resources XII 53 two separate equations for saturation and specific discharge. One of the major advantages of this formulation is that the model becomes more accessible to analysis. Adding the equations (2) for both phases and using (3) we obtain V-v = 0, in which v = v% + v%, is the total specific discharge. To obtain the equation for air saturation, (2) is rewritten as v/ / = &,%;, (5) where e~ is the unity vector in the positive ^-direction and the mobility of phase /, A/ [m^ Pa~* s~*], is defined as, Kki. A/ = -, / = a, w. Vl We subtract (5), / = w, from (5), / = a, to obtain v% = A[v + A^((/)^ - /)%)#ez - Vpc)], where the air fractional flow function fa [ ] is defined by A, fa = a + The systems (1) and (2) can now be rewritten in terms of air saturation and total specific discharge: <Xx) - (% + V [F(x, 5% v) - v A(x, / 2)Vpc(x, ^ \ yj ^)]? = 0, ^g^ V-v = 0, _, where The vector-valued function F [m s *] and the scalar function A [m^ Pa~^ s~^] are respectively referred to as the flux and mobility function. The flux function represents convective transport of the air phase and includes gravitational effects. Its non-linear behaviour is illustrated in Figure 1 in which a graph of F% is depicted
4 54 Computer Methods in Water Resources XII 4.5,x Figure 1: Typical flux function S for K = ^ nf and v, = mm s"*. The values used for the other relevant parameters can be found listed at the end of this article. Figure 1 shows that F? is not necessarily monotone in S. This is due to the presence of gravity in the system. FZ is an increasing function if, and only if, (see [5]). From numerical experiments it appears that our system is dominated by convection. Note that the saturation equation changes from a parabolic into a hyperbolic partial differential equation with the disappearance of the capillary pressure. This observation will enable us to identify appropriate numerical techniques. 3 One-dimensional modelling In this article the simple one-dimensional instationary case is studied in order to gain insight in the various physical and numerical problems
5 Computer Methods in Water Resources XII 55 encountered in modelling the air sparging process. We will consider vertical upward flow resulting from a constant vertical injection of air at a certain depth. Thus, consider the one-dimensional fractional flow model for air saturation resulting from (6): dv + F(Z, S) - X(z, S)-pc(z, S) = 0, (7) where Ffz,^) = A(^)t; + A(z,^)(^-/)J^. From (8) it is clear that v is uniform. Suppose that (7) holds for -A& < z < ht. The lower boundary z = hb is located at the position of the vertical air injection, i.e. the boundary condition at z = hf> is r\ f (z, f) - A(z, ^)g;pc(^ 61 = ^. (9) From this it follows that v = v^. We locate the upper boundary sufficiently high above the natural water table at z - 0. We assume the natural no-flow equilibrium as initial condition. Thus, initially the water table is located at z - 0, i.e. S(z) = 0 if z < 0. The noflow saturation distribution above the water table now follows as the solution of the one-dimensional analogon of (2). Thus, p^ = -p^gz and pa = -pagz for z > 0, and from thefirstequation of (4) it follows that 4 The Godunov scheme Since (7) is assumed to be dominated by convection, we will consider the hyperbolic conservation law ^W + A^) = o do) with the boundary condition F(z, S] - ^ at z = -h\>. To discretise (10) the domain [-/^,/^] is divided into M cells of equal width Az. This gives M nodes z, = -^ + (; - l/2)az, ; = 1,...,M, where Az - (hb + ht)/m. The nodes Zj are the centres of the grid cells. Using this block-centred approach, the air saturation S is approximated
6 56 Computer Methods in Water Resources XII at the cell centres, giving unique values per cell. After collecting those values, a piecewise constant approximation for the saturation profile is constructed. Time-sampling takes place at t^ =raatf,n =0,1,..., where A/ is the time step. We write S for a numerical approximation of the exact solution S(zj,tn)- The discretisation scheme considered is the Godunov scheme (see [6] and [3]). The scheme is explicit, which is a very attractive feature. However, the time step is limited by the Courant-Friedrichs-Lewy condition CFLj < 1, where max is the Courant number. The Godunov scheme is based on solving Riemann problems at the interfaces ~j+i/2 between grid cells. At each time level the numerical solution of the previous time level serves as initial data. The solutions of the Riemann problems are assembled to form a piecewise constant numerical solution at the new time level. Provided that CFLj < 1, j = 1,..., Af, the Godunov scheme is ; S7_i,S7))], (11) where <j>j = </>(zj) and Fj(S) = F(z^S). The function STZ(»; S,S+i) is denned by where W is the solution of the Riemann problem for the interface O TT7 ^ +1 Equation (11) results from averaging the juxtaposition 5^T*~* of the Riemann solutions ^( j 5^"*""^, S^^} over the corresponding grid cells at time level /nii i-e.
7 Computer Methods in Water Resources XII 57 / = 240 = Figure 2: Solution for the Godunov scheme It is possible to compute S'^ without solving Riemann problems explicitly. The construction of S +* has to be chosen differently in order to incorporate the inflow boundary condition (9):?n+l _ A/ 5 A numerical experiment Let A%z) = 5.3'10-^m^ifz> Oorz < - if 0 < z < 1, and the other parameters be chosen as listed at the end of this article. The CFL-condition results into the requirement Af < Az. We choose Az = m and Af = 7.5 s. The results are displayed in Figure 2.
8 58 Computer Methods in Water Resources XII Parameter default values symbol g ht m a Pa meaning acceleration due to gravity depth of computational domain height of computational domain Van Genuchten soil constant Darcy injection flux in the air phase Van Genuchten soil parameter viscosity of the air phase viscosity of the water phase density of the air phase density of the water phase porosity value 9.81 m s"2 6 m 4 m 2/ m! 2m-: s Pa s Pa s 1.24 kg m-3 1 # kg m References [1] Van Dijke, M.I.J., Van der Zee, S.E.A.T.M. & Van Duijn, C.J., Multi-phase flow modelling of air sparging, Advances in Water 18, pp , [2] Van Genuchten, M.Th., A closed form equation for predicting the hydraulic conductivity of unsaturated soils, Soil Science So- 4, pp , [3] Godlewski, E. & Raviart, P. -A., Hyperbolic Systems of Conservation Laws, Ellipses, Paris, [4] Hinchee, R.E., Miller, R.N. & Johnson, P.C., fn 5% Air Sparging, Bioventing, and Related Remediation Processes, Battelle Press, Columbus, [5] Kaasschieter, E.F., Van der Werff ten Bosch, J.D., & Mulder, G.J., A numerical fractional flow model for air sparging, Report RAN A of the Department of Mathematics and Computing Science, Eindhoven University of Technology, [6] LeVeque, R.J., Numerical Methods for Conservation Laws, Birkhauser Verlag, Basel, 1990.
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 informationTodd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
CONSERVATIVE CHARACTERISTIC METHODS FOR LINEAR TRANSPORT PROBLEMS Todd Arbogast Department of Mathematics and, (ICES) The University of Texas at Austin Chieh-Sen (Jason) Huang Department of Applied Mathematics
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 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 informationESTIMATION OF TRANSVERSE DISPERSTVITY IN THE MIXING ZONE OF FRESH-SALT GROUNDWATER
ModelCARE 90: Calibration and Reliability in Groundwater Modelling (Proceedings of the conference held in The Hague, September 1990). IAHS Publ. no. 195, 1990. ESTIMATION OF TRANSVERSE DISPERSTVITY IN
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
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 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 informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
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 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 informationCHARACTERIZATION OF FRACTURES IN GEOTHERMAL RESERVOIRS USING RESISTIVITY
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 01 SGP-TR-194 CHARACTERIZATION OF FRACTURES IN GEOTHERMAL RESERVOIRS
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 informationStreamline calculations. Lecture note 2
Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational
More informationApplications of Partial Differential Equations in Reservoir Simulation
P-32 Applications of Partial Differential Equations in Reservoir Simulation Deepak Singh Summary The solution to stochastic partial differential equations may be viewed in several manners. One can view
More informationFlood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum
Archives of Hydro-Engineering and Environmental Mechanics Vol. 56 (29), No. 3 4, pp. 121 137 IBW PAN, ISSN 1231 3726 Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum Dariusz
More information13.7 Power Applied by a Constant Force
13.7 Power Applied by a Constant Force Suppose that an applied force F a acts on a body during a time interval Δt, and the displacement of the point of application of the force is in the x -direction by
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 informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationB024 RESERVOIR STREAMLINE SIMULATION ACCOUNTING
1 B024 RESERVOIR STREAMLINE SIMULATION ACCOUNTING FOR EFFECTS OF CAPILLARITY AND WETTABILITY R.A. BERENBLYUM, A.A. SHAPIRO, E.H. STENBY IVC-SEP, Department of Chemical Engineering, Technical University
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 informationBuilding ground level
TMA4195 MATHEMATICAL MODELLING PROJECT 212: AQUIFER THERMAL ENERGY STORAGE 1. Introduction In the project we will study a so-called Aquifer Thermal Energy Storage (ATES) system with the aim of climitizing
More informationTwo Phase Transport in Porous Media
Two Phase Transport in Porous Media Lloyd Bridge Iain Moyles Brian Wetton Mathematics Department University of British Columbia www.math.ubc.ca/ wetton CRM CAMP Seminar, October 19, 2011 Overview Two phase
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 informationRupelton layer. 1 An introduction to the problem and some basic definitions and an introduction to the problem
Abstract In this work, we use the finite volume code pdelib/thermoha2 of Enchery and Fuhrmann [2] to understand the influence of layers of small thichness on the behaviour of the fluid. To this end consider
More informationCOMPUTER MODELLING OF HEAT AND MASS FLOW IN STEAMING GROUND AT KARAPITI THERMAL AREA, NEW ZEALAND
PROCEEDINGS, Twenty-Ninth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 26-28, 2004 SGP-TR-175 COMPUTER MODELLING OF HEAT AND MASS FLOW IN STEAMING GROUND
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 informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More information18.02 Multivariable Calculus Fall 2007
MIT OpenCourseWare http://ocw.mit.edu 18.02 Multivariable Calculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. V9. Surface Integrals Surface
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 informationMMsFEM and Streamlines
MMsFEM and Streamlines Combining a mixed multiscale FEM with streamline simulation for enhanced reservoir performance prediction on large grid models. Jørg E. Aarnes, Vegard Kippe and Knut-Andreas Lie
More informationEVALUATION OF CRITICAL FRACTURE SKIN POROSITY FOR CONTAMINANT MIGRATION IN FRACTURED FORMATIONS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationAbstract. This research is supported by US Army Corps of Engineers under Grant No. DACA K-0055 with Penn State.
Modeling Water Flow and Chemical and Sediment Transport in Watershed Systems Hwai-Ping Cheng and Gour-Tsyh Yeh Department of Civil and Environmental Engineering, The Pennsylvania State University, University
More 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 informationThermo hydro mechanical coupling for underground waste storage simulations
Thermo hydro mechanical coupling for underground waste storage simulations Clément Chavant - Sylvie Granet Roméo Frenandes EDF R&D 1 Outline Underground waste storage concepts Main phenomena and modelisation
More informationTwo-dimensional flow in a porous medium with general anisotropy
Two-dimensional flow in a porous medium with general anisotropy P.A. Tyvand & A.R.F. Storhaug Norwegian University of Life Sciences 143 Ås Norway peder.tyvand@umb.no 1 Darcy s law for flow in an isotropic
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 08-08 March2008 The complexe flux scheme for spherically symmetrie conservation laws by J.H.M. ten Thije Boonkkamp,
More informationJ. Environ. Res. Develop. Journal of Environmental Research And Development Vol. 8 No. 1, July-September 2013
SENSITIVITY ANALYSIS ON INPUT PARAMETERS OF 1-D GROUNDWATER FLOW GOVERNING EQUATION : SOLVED BY FINITE-DIFFERENCE AND THOMAS ALGORITHM IN MICROSOFT EXCEL Goh E.G.* 1 and Noborio K. 2 1. Department of Engineering
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 informationTravelling wave solutions for the Richards equation incorporating non-equilibrium effects in the capillarity pressure
UHasselt Computational Mathematics Preprint eries Travelling wave solutions for the Richards equation incorporating non-equilibrium effects in the capillarity pressure Hans van Duijn, Koondanibha Mitra
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 informationLecture 16 Groundwater:
Reading: Ch 6 Lecture 16 Groundwater: Today 1. Groundwater basics 2. inert tracers/dispersion 3. non-inert chemicals in the subsurface generic 4. non-inert chemicals in the subsurface inorganic ions Next
More informationUnsaturated Flow (brief lecture)
Physical Hydrogeology Unsaturated Flow (brief lecture) Why study the unsaturated zone? Evapotranspiration Infiltration Toxic Waste Leak Irrigation UNSATURATAED ZONE Aquifer Important to: Agriculture (most
More informationCranfield ^91. College of Aeronautics Report No.9007 March The Dry-Bed Problem in Shallow-Water Flows. E F Toro
Cranfield ^91 College of Aeronautics Report No.9007 March 1990 The Dry-Bed Problem in Shallow-Water Flows E F Toro College of Aeronautics Cranfield Institute of Technology Cranfield. Bedford MK43 OAL.
More informationI. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010
: an : an (Joint work with A. Fasano) Dipartimento di Matematica U. Dini, Università di Firenze (Italy) borsi@math.unifi.it http://web.math.unifi.it/users/borsi porous EMS SCHOOL ON INDUSTRIAL MATHEMATICS
More informationLocal Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf**
Controlled CO 2 Diversified fuels Fuel-efficient vehicles Clean refining Extended reserves Local Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf** *: Laboratoire J.L.Lions **:Université
More informationVisualization of Natural Convection in Enclosure. Filled with Porous Medium by Sinusoidally. Temperature on the One Side
Applied Mathematical Sciences, Vol., 2012, no. 97, 801-812 Visualization of Natural Convection in Enclosure Filled with Porous Medium by Sinusoidally Temperature on the One Side Paweena Khansila Department
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 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 informationDarcy s Law. Darcy s Law
Darcy s Law Last time Groundwater flow is in response to gradients of mechanical energy Three types Potential Kinetic Kinetic energy is usually not important in groundwater Elastic (compressional) Fluid
More informationThe RAMSES code and related techniques 4. Source terms
The RAMSES code and related techniques 4. Source terms Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term - Relaxation towards
More informationChapter 3 Permeability
3.2 Darcy s Law In 1856, Darcy investigated the flow of water through sand filters for water purification. His experimental apparatus is shown in Figure 3.11. By empirical observation Figure 3.11 Schematic
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
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 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 informationTutorial for the supercritical pressure pipe with STAR-CCM+
Tutorial for the supercritical pressure pipe with STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities with
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 informationPORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM
Excerpt from the Proceedings of the COMSOL Conference 2010 Paris PORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM Igor Bogdanov 1*, Sylvain Jardel 1, Anis Turki 1, Arjan Kamp 1 1 Open &
More informationNUMERICAL SOLUTION OF TWO-REGION ADVECTION-DISPERSION TRANSPORT AND COMPARISON WITH ANALYTICAL SOLUTION ON EXAMPLE PROBLEMS
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 130 137 NUMERICAL SOLUTION OF TWO-REGION ADVECTION-DISPERSION TRANSPORT AND COMPARISON WITH ANALYTICAL SOLUTION ON EXAMPLE PROBLEMS
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department ofmathematîcs and Computer Science. CAS A-Report April 2005
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department ofmathematîcs and Computer Science CAS A-Report 5-16 April 25 Analysis of wet pressing of paper: the three-phase model. Part I: constant air density by D.
More informationNUMERICAL TREATMENT OF TWO-PHASE FLOW IN CAPILLARY HETEROGENEOUS POROUS MEDIA BY FINITE-VOLUME APPROXIMATIONS
NUMERICAL TREATMENT OF TWO-PHASE FLOW IN CAPILLARY HETEROGENEOUS POROUS MEDIA BY FINITE-VOLUME APPROXIMATIONS HELMER ANDRÉ FRIISA AND STEINAR EVJE B,C Abstract. This paper examines two-phase flow in porous
More informationTemperature and Internal Heat Generation in a Porous Medium
TECHNISCHE MECHANIK, Band 21, Heft 4, (20011313318 Manuskripteingang: 20. September 200] Free Convection over 3 Vertical Flat Plate with a Variable Wall Temperature and Internal Heat Generation in a Porous
More informationConvergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates
Convergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates Xin Wen and Shi Jin Abstract We study the L 1 -error estimates
More informationDISCRETE ELEMENT SIMULATIONS OF WATER FLOW THROUGH GRANULAR SOILS
15th ASCE Engineering Mechanics Conference June 2-5, 2002, Columbia University, New York, NY EM 2002 DISCRETE ELEMENT SIMULATIONS OF WATER FLOW THROUGH GRANULAR SOILS Usama El Shamy 1, Student Member,
More informationEffect of Sorption/Curved Interface Thermodynamics on Pressure transient
PROCEEDINGS, Twentieth Workshop on Geothermal Rey~volr Englneerlng Stanford Unhrenlty, Stanfoni, Callfornla, January 2426 1995 SGP-m-150 Effect of Sorption/Curved Interface Thermodynamics on Pressure transient
More informationThe relative connectivity-tortuosity tensor for conduction of water in anisotropic unsaturated soils.
Final draft 20-01-2004 For submittal to Vadose Zone Journal The relative connectivity-tortuosity tensor for conduction of water in anisotropic unsaturated soils. Peter A.C. Raats 1, Z. Fred Zhang 2, Andy.
More informationAn Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term
An Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term Alexander Kurganov Department of Mathematics, Tulane University, 683 Saint Charles Avenue, New Orleans, LA 78,
More informationTransactions on Modelling and Simulation vol 13, 1996 WIT Press, ISSN X
The cooling of molten glass in a mould Ph. Simons, R. Mattheij Scientific Computing Group, Eindhoven University of Technology, 56 MB Eindhoven, The Netherlands Abstract The actual formation of glass jars
More informationExamination paper for TPG4150 Reservoir Recovery Techniques
1 Department of Petroleum Engineering and Applied Geophysics Examination paper for TPG4150 Reservoir Recovery Techniques Academic contact during examination: Jon Kleppe Phone: 91897300/73594925 Examination
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 informationNUMERICAL AND EXPERIMENTAL ANALYSIS OF SEEPAGE BENEATH A MODEL OF A GRAVITY DAM
Engineering Review Vol. 33, Issue 2, 75-84, 2013. 75 NUMERICAL AND EXPERIMENTAL ANALYSIS OF SEEPAGE BENEATH A MODEL OF A GRAVITY DAM T. Jelenkovi V. Travaš * Chair of Hydraulic Engineering, Faculty of
More informationAdvanced numerical methods for transport and reaction in porous media. Peter Frolkovič University of Heidelberg
Advanced numerical methods for transport and reaction in porous media Peter Frolkovič University of Heidelberg Content R 3 T a software package for numerical simulation of radioactive contaminant transport
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2011 Randall J. LeVeque Applied Mathematics University of Washington January 3, 2011 R.J. LeVeque, University of Washington AMath 574,
More information&
Inverse Problem for Transient Unsaturated Flow Identifiability and Non Uniqueness K. S. Hari Prasad & M. S. Ghidaoui Department of Civil and Structural Engineering The Hong Kong University of Science and
More informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More information1. INTRODUCTION TO CFD SPRING 2019
1. INTRODUCTION TO CFD SPRING 2019 1.1 What is computational fluid dynamics? 1.2 Basic principles of CFD 1.3 Stages in a CFD simulation 1.4 Fluid-flow equations 1.5 The main discretisation methods Appendices
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 informationA FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS
A FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS Ruben Juanes 1 and Knut-Andreas Lie 2 1 Stanford University, Dept. Petroleum Engineering, USA 2 SINTEF IKT, Dept., Norway ECMOR IX, August 30 September
More informationPaper E A pore network model for calculation of interfacial velocities
Paper E A pore network model for calculation of interfacial velocities Submitted to Advances in Water Resources, fall 21. A PORE NETWORK MODEL FOR CALCULATION OF INTERFACIAL VELOCITIES H.F. Nordhaug a,
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 informationMembrane processes selective hydromechanical diffusion-based porous nonporous
Membrane processes Separation of liquid or gaseous mixtures by mass transport through membrane (= permeation). Membrane is selective, i.e. it has different permeability for different components. Conditions
More informationEFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL Vahid SHOKRI 1*,Kazem ESMAEILI 2 1,2 Department of Mechanical Engineering, Sari Branch, Islamic Azad
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
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 informationSENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE. Juan E. Santos
SENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE Juan E. Santos Instituto del Gas y del Petróleo, Facultad de Ingeniería UBA and Department
More informationTutorial for the heated pipe with constant fluid properties in STAR-CCM+
Tutorial for the heated pipe with constant fluid properties in STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
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 informationCHAPTER 135 NUMERICAL SIMULATION OF WAVE MOTION ON AND IN COASTAL STRUCTURES
CHAPTER 135 NUMERICAL SIMULATION OF WAVE MOTION ON AND IN COASTAL STRUCTURES J.W. van der Meer 1 ', H.A.H. Petit 11, P. van den Bosch 1 ', G. Klopman 1 ' and R.D. Broekens 2 ' Abstract A 2-dimensional
More informationSemi-analytical solutions of contaminant transport equation with nonlinear sorption in 1D
Semi-analytical solutions of contaminant transport equation with nonlinear sorption in D Peter Frolkovič and Jozef Kačur November 9, 24 Abstract A new method to determine semi-analytical solutions of one-dimensional
More informationBasic equations of motion in fluid mechanics
1 Annex 1 Basic equations of motion in fluid mechanics 1.1 Introduction It is assumed that the reader of this book is familiar with the basic laws of fluid mechanics. Nevertheless some of these laws will
More informationUNIT II Real fluids. FMM / KRG / MECH / NPRCET Page 78. Laminar and turbulent flow
UNIT II Real fluids The flow of real fluids exhibits viscous effect that is they tend to "stick" to solid surfaces and have stresses within their body. You might remember from earlier in the course Newtons
More informationPore-scale modeling extension of constitutive relationships in the range of residual saturations
WATER RESOURCES RESEARCH, VOL. 37, NO. 1, PAGES 165 170, JANUARY 2001 Pore-scale modeling extension of constitutive relationships in the range of residual saturations Rudolf J. Held and Michael A. Celia
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationFluid Mechanics Introduction
Fluid Mechanics Introduction Fluid mechanics study the fluid under all conditions of rest and motion. Its approach is analytical, mathematical, and empirical (experimental and observation). Fluid can be
More informationDevelopment of a Thermo-Hydro-Geochemical Model for Low Temperature Geoexchange Applications
Development of a ThermoHydroGeochemical Model for Low Temperature Geoexchange Applications F. Eppner *1, P. Pasquier 1, and P. Baudron 1 1 Department of Civil, Geological and Mining Engineering, Polytechnique
More informationThe Deflation Accelerated Schwarz Method for CFD
The Deflation Accelerated Schwarz Method for CFD J. Verkaik 1, C. Vuik 2,, B.D. Paarhuis 1, and A. Twerda 1 1 TNO Science and Industry, Stieltjesweg 1, P.O. Box 155, 2600 AD Delft, The Netherlands 2 Delft
More informationGeneralised Separable Solution of Double Phase Flow through Homogeneous Porous Medium in Vertical Downward Direction Due to Difference in Viscosity
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 932-9466 Vol. 8, Issue (June 203), pp. 305-37 Applications and Applied Mathematics: An International Journal (AAM) Generalised Separable Solution
More informationLarge Time Step Scheme Behaviour with Different Entropy Fix
Proceedings of the Paistan Academy of Sciences: A. Physical and Computational Sciences 53 (): 13 4 (016) Copyright Paistan Academy of Sciences ISSN: 0377-969 (print), 306-1448 (online) Paistan Academy
More information