Todd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
|
|
- Charity Kellie Gibson
- 5 years ago
- Views:
Transcription
1 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 and National Center for Theoretical Sciences National Sun Yat-sen University (Taiwan)
2 Outline 1. The PDE s of Transport Problems and Local Conservation Principles 2. Characteristic Methods for Approximating the solution 3. Satisfying the Local Volume Conservation Principle Discretely 4. Numerical Results 5. Conclusions
3 Transport Problems and Local Conservation
4 Suppose Conservative Fluid Flow ξ is a conserved quantity ξ (mass/volume) v is the fluid velocity (length/time) ξv is the flux of ξ (mass/area/time) Q is an external source or sink of fluid (mass/volume/time) Within a region of space R, the total amount of ξ changes in time by d dt ξ dx = ξv ν da(x) + Q dx } R{{}} R {{}} R {{} Change in R Flow across R Sources/sinks = conservation locally on R R ξ t dx = (ξv) dx } R {{} Divergence Theorem This is true for each region R, so in fact ξ t + (ξv) = Q + R Q dx
5 A Transport Problem 1 One incompressible fluid (tracer) flowing miscibly in another incompressible fluid, within an incompressible medium. Velocity of the bulk fluid. Conservation of bulk fluid mass (ξ = φρ) gives ξ t + (ξv) = Q = u = q u is the (unknown) bulk fluid velocity (v = u/φ) φ is the porosity (constant in time) ρ is the (constant) density q is the source/sink (wells, Q = ρq) Simple Tracer Transport. Conservation of tracer mass (ξ = φc) gives φc t + (cu) = c I q + + cq q c (c) c is the (unknown) tracer concentration c I is the given concentration of injected fluid q + /q is q when positive/negative
6 A Transport Problem 2 However, transport is not the only process occurring! Mass flux. v = cu D c (Transported plus Diffusive Flux) where D is the diffusion/dispersion coefficient Chemical reactions. q = q c (c) + R(c) (Wells plus Reactions) where R is the reaction term Tracer Transport. Conservation of tracer mass gives φc t + (cu D c) = q c (c) + R(c)
7 Operator Splitting of Transport Equation 1 φc t + (cu D c) = q c (c) + R(c) Discretization in time: t > 0 and t n = n t. We want to solve the transport and reactive part of the equation explicitly and the diffusive part implicitly. Thus, we want φ cn+1 c n + (c n u) (D c n+1 ) = q c (c n ) + R(c n ) t This is equivalent to the three steps (Reaction) (Transport) cn φ c t = R(c n ) φc t = R(c) φĉ c t + (cn u) = q c (c n ) φc t + (cu) = q c (c) (Diffusion) φ cn+1 ĉ (D c n+1 ) = 0 φc t (D c) = 0 t with some intermediate c and ĉ.
8 Operator Splitting of Transport Equation 2 Nonlinear Ordinary Differential Equation part (Reaction) φc t = R(c) Linear Hyperbolic part (Transport) φc t + (cu) = q c (c) Linear Parabolic part (Diffusion/dispersion) φc t (D c) = 0 We discuss approximations of the transport step only.
9 Locally Conservative Methods A locally conservative method is one for which the approximate solution satisfies the conservation principle, but only over certain discrete regions. Normally, one would take the grid elements R and require ( φct + (cu) ) dx = q dx R R but we will need to be more general than this. Remark: The reactive and diffusive steps must also be solved by locally conservative methods, or local conservation will break down!
10 Characteristic Methods for Linear Transport
11 Characteristic Tracing of Points The characteristic trace-forward of the point x is denoted ˇx = ˇx(x; t). It satisfies the ordinary differential equation dˇx u(ˇx, t) = dt φ(ˇx), tn < t t n+1 ˇx(t n ) = x In the absence of sources/sinks and diffusion, fluid particles simply travel along the characteristics of the equation. Time t n+1 ˇx t n x Space The concentration is constant along this space-time path, since dc(ˇx, t) dt = c t dˇx + c dt = c t + c u φ = 1 ( φct + u c ) = 0 φ
12 Characteristic Trace-back of Points The characteristic trace-back of the point x is denoted ˆx = ˆx(x; t). It satisfies the (time backward) ordinary differential equation dˆx dt = u(ˆx, t) φ(ˆx), tn t < t n+1 ˆx(t n+1 ) = x In the absence of sources/sinks and diffusion, fluid particles simply travel along the characteristics of the equation. Time t n+1 x t n ˆx Space Again, the concentration is constant along this space-time path.
13 Modified Method of Characteristics (MMOC) (Douglas and Russell, 1982) Key idea: Use a finite difference approximation of the characteristic derivative dc dt c t(x, t) + u(x, t) This results in the approximation φ c(x, t) c(x, t + t) c(ˆx, t) t c(x, t + t) c(ˆx, t) φ t = (c I c)q + t n+1 c(x) at each grid point t n c(ˆx) Problems: Because the method is based on points, it violates local mass conservation constraints for both the bulk fluid and the tracer.
14 Characteristic Trace-back of Regions To obtain mass conservation,... Key idea: Trace regions rather than points! The particles in a grid element E trace back to a region Ê Ê = {ˆx Ω : ˆx = ˆx(x; t n ) for some x E}. In space and time, we actually trace a region E = E(E) given by E = {(ˆx, t) Ω [t n, t n+1 ] : ˆx = ˆx(x; t) for some x E}. t n+1 E t n Ê E
15 Local Mass Conservation of the Tracer φc t + (cu) = q c (c) Integrate in space-time over E and use the divergence theorem ( ) ( ) cu u x,t dx dt = E φc c ν E φ x,t dσ ( ) u = E φcn+1 dx Ê φcn dx + c ν S φ x,t dσ The( last ) term is integration on the space-time sides S of E, u but is orthogonal to ν φ x,t there! The local mass constraint: E φcn+1 dx = + Ê φcn dx E q c dx dt S B cu ν dσ t n+1 t n S B Ê E E ν x,t
16 Local Mass Conservation of the Bulk Fluid A similar local mass constraint holds for the bulk fluid (c 1) φ t + u = u = q Since we are dealing with incompressible fluids, we call this the local volume constraint. The local volume constraint: E φ dx = Ê φ dx + E q dx dt S B u ν dσ
17 Characteristics Mixed Method (CMM) (A., Chilakapati, and Wheeler, 1992; A. and Wheeler, 1995) Use lowest order Raviart-Thomas mixed finite elements. The scalar test function is a constant on each element in space. E φcn+1 dx = Ê φcn dx + q c dx dt cu ν dσ E S B Remark: Practical implementation requires that Ê be approximated by Ẽ Ê, a polygon. This is equivalent to modifying the velocity field, so tracer mass is still conserved locally by the above equation Ê Ẽ Vol(Ê) Vol(Ẽ)
18 Volume Error Problem: The local volume constraint may be violated for this perturbed velocity! This is because the volume constraint does not enter into the definition of the method. This leads to incorrect densities of the tracer, which leads, e.g., to bad approximation of reaction dynamics Relative volume errors around an injection well.
19 Local Volume Conservation in Characteristic Methods
20 Volume Conservation Key idea: To obtain volume conservation, define Ẽ Ê to be a simple shape that satisfies the volume constraint E φ dx Ẽ φ dx = E q dx dt S B u ν dσ Strategy: Suppose that E is a rectangle. Perturb the vertices and midpoints of Ê only a little so that we get a polygon Ẽ with 8 sides such that the above constraint is satisfied. We call this method the Volume Corrected Characteristics Mixed Method (VCCMM) Ê Ẽ Vol(Ê) = Vol(Ẽ) Problem: It is easy to introduce systematic bias into the flow field and thereby produce unphysical flows. We must do the adjustment very carefully!
21 An Example of Unphysical Flow Volume not conserved Volume conserved 10 years years Biased trace-back adjustment has introduced unphysical flow corresponding to a large, incorrect velocity channel.
22 Forward Trace of Injection Wells Most of the error is near injection wells. Characteristic tracing back in time traces into the well, which is difficult to approximate. Key idea 1: Trace the well forward (out of the well). Adjust region to satisfy the volume constraint (cf. Healy and Russell, 2000). The characteristic trace-forward of the point x is denoted ˇx = ˇx(x; t), and it satisfies the ordinary differential equation dˇx u(ˇx, t) = dt φ(ˇx), tn < t t n+1 ˇx(t n ) = x
23 Conservation Near Injection Wells Well volume constraint: Adjust W so W φ dx = W φ dx + E f q dx dt E Ẽ W W Element volume constraint: Adjust Ẽ so E φ dx = Ẽ φ dx+ Vol(E W) Vol( W) E q dx dt Transport: E φcn+1 dx = Ẽ φcn dx + Vol(E W) Vol( W) E q c dx dt
24 Inflow boundaries Like injection wells, inflow boundaries trace back out of the domain. Idea: Either trace inflow boundaries forward, or fold the time axis down to the xy-plane to create a ghost region. y y t t n t n+1 x t x u ν φ Ẽ Ω Ω Volume constraint: Replace φ by u ν in the ghost region: E Ẽ Ω φ dx φ dx + u ν dσ = S B E Ẽ φ dx φ dx = Mass constraint: Replace φc n by c n Iu ν in the ghost region: E φcn+1 dx Ẽ φcn dx = q c dx dt. E E q dx dt,
25 Trace-back Point Adjustment Key idea 2: Adjust points in the direction of the flow; that is, along the characteristics in time (cf. Douglas, Huang, and Pereira, 1999). To define x, for τ n t n, we solve (backwards) d x dt = u( x, t) φ( x), τn < t t n+1 x(t n+1 ) = x We convert space error into time error: t n+1 x t n τ n ˆx x
26 Trace-back (or Forward) Point Adjustment 1 Proceed away from injection wells and inflow boundaries by layers. For each layer, obtain volume conservation in two steps. 1. Volume conservation of the layer. Adjust the exterior contour of the entire layer along the characteristics until the volume of the layer is correct (within a small tolerance). That is, in the absence of other sources, inflow boundaries, and sinks, E in the layer E φ dx = E in the layer Ẽ φ dx Flow Adjusted point (fixed) Points adjusted simultaneously in the direction of the characteristic (we use a type of bisection algorithm) Points adjusted individually transverse to the flow in Step 2
27 Trace-back (or Forward) Point Adjustment 2 2. Element volume conservation. Within the layer, sequentially adjust the interior midpoint of each element transverse to the flow until the volume of the element is correct (within a small tolerance). Flow > x i,j < > xi,j+1/ x i,j 3 Adjusted point (fixed) Points individually adjusted transverse to the flow Remark: This is an extremely fast direct algorithm.
28 Numerical Results
29 Numerical Results Darcy s Law completes the equations: u = k µ p p is the fluid pressure k is the permeability µ is the fluid viscosity Measure the variability of k by the dimensionless coefficient of variation where the mean is C v = 1 M k ( 1 Vol(Ω) M k = 1 Vol(Ω) Ω (k(x) M k) 2 dx Ω k(x) dx ) 1/2
30 A Nuclear Contamination Problem 1 The permeability is log-normal and fractal. M k = cm 2 (about 20 md) C v = (varies over five orders of magnitude). 256 Inflow Y 64 1E-08 1E-09 1E-10 1E-11 1E-12 1E-13 Outflow Injection well
31 A Nuclear Contamination Problem 2 If we use a small time step of t = 1.5 years, we can trace back into the injection well CMM volume errors VCCMM (volume errors 10 9 )
32 A Nuclear Contamination Problem 3 Concentration at 30 years on a grid, with t = 1.5 yr 2.6CFL E E E E E E E E E E E E E E E E CMM VCCMM CMM overshoots the maximum concentration of 1E-5 by 34% up to 1.34E-5.
33 The Courant-Friedrichs-Lewy (CFL) Condition Explicit methods in 1-dimension have a time step restriction, known as the Courant-Friedrichs-Lewy (CFL) time-step, given by where h is the grid spacing. t t CFL,1-D = max x Ω hφ(x) u(x) In 2-dimensions, we should limit t to half this value, t t CFL,2-D = max x Ω hφ(x) 2 u(x) Godunov s method is a popular method, that is unstable if the CFL condition is violated. In principle, characteristic methods are not subject to this constraint, and large time steps can be used.
34 A Nuclear Contamination Problem 4 Concentration at 30 years on a grid E E E E E E E E E E E E E E Godunov t = yr = 1 CFL VCCMM-TF t = 3 yr = 5.1 CFL We use trace-forwarding near the well. No overshoot for either method. Less numerical diffusion for VCCMM-TF. 51 Godunov steps vs. 10 for VCCMM-TF.
35 A Nuclear Contamination Problem 5 Concentration at 30 years on a grid E E E E E E E E E E E E E E Godunov t = yr = 1 CFL VCCMM-TF t = 1 yr = 6.8 CFL We use trace-forwarding near the well. No overshoot for either method. Less numerical diffusion for VCCMM-TF. 205 Godunov steps vs. 30 for VCCMM-TF.
36 A Nuclear Contamination Problem 6 Concentration at 30 years on a grid E E E E E E E E E E E E E E Godunov t =.0366 yr = 1 CFL VCCMM-TF t = 0.5 yr = 13.7 CFL We use trace-forwarding near the well. No overshoot for either method. Less numerical diffusion for VCCMM-TF. 820 Godunov steps vs. 60 for VCCMM-TF.
37 A Quarter Five-Spot Problem 1 Geostatistically generated permeability. M k = 100 md C v = 2.58 (varies over four orders of magnitude). 5E-12 1E-12 5E-13 1E-13 5E-14 1E-14 5E-15 1E-15
38 A Quarter Five-Spot Problem 2 Concentration at 3.36 years using t = yr = CFL CMM VCCMM-TF CMM shows both overshoot and undershoot. Very large initial volume imbalances throughout the domain. If t = yr = CFL initially creates degenerate trace-back regions, which cannot be used.
39 A Linear Flood Problem 1 A test with an inflow boundary. Permeability field of the quarter five-spot problem Linear pressure drop across the domain in the x-direction. Concentration at 25 years using t = 0.18 year = 3 CFL CMM VCCMM CMM has severe overshoots up to c = Initial volume errors exceed 10% in the interior of the domain. If t = 9 CFL, initial relative volume errors are around 25%.
40 A Linear Flood Problem 2 Concentration at 25 years using trace-forwarding of the inflow boundary VCCMM-TF t = 0.18 yr = 3 CFL VCCMM-TF t = 0.36 yr = 6 CFL
41 A Fluvial Domain Problem 1 Domain feet 2 solved on a grid. Permeability of 3 values, M k = Darcy and C v = φ = 0.2. Wells in opposite corners, injecting 1 pore volume every 3 years. t = years = CFL The permeability, in Darcies
42 A Fluvial Domain Problem 2 VCCMM-TF concentration Time 1.05 years (step 70) Time 1.65 years (step 110) Remark: CMM alone produces negative concentrations on a grid with t = 0.01 year, indicating that the trace-back regions self intersect.
43 A Fluvial Domain Problem 3 CMM with trace-forwarding of wells only Time 1.05 years (step 70) Time 1.65 years (step 110)
44 Comparison with an Analytic Solution Comparison of VCCMM-TF concentration with an analytic solution for radial flow from a well in a horizontally infinite, uniform porous medium Analytic Analytic VCCMM-TF 0.6 VCCMM-TF Longitudinal dispersion 20 cm Longitudinal dispersion 100 cm The analytic solution is derived by a code from P.A. Hsieh (1986). VCCMM-TF uses t = 8 CFL.
45 Conclusions 1. A critical aspect of approximating hyperbolic transport problems is to conserve the mass of the tracer locally. 2. It is just as critical to the conserve locally the mass of the bulk fluid. 3. When adjusting trace-back points, care must be used to avoid introducing systematic bias into the transport computation. The key is to adjust along the characteristics when possible. 4. Trace-forwarding is needed around injection wells. 5. Inflow boundaries can be treated transparently through a space-time fold-down strategy, or with trace-forwarding. 6. In principle, the only restriction on the time step is that the traceback regions not degenerate or self-intersect (this problem can be alleviated by tracing more points of the boundary E). 7. VCCMM-TF allows us to take large time steps (such as 14 times the 2-D CFL limited time step) with no overshoots, which results in less numerical dispersion than Godunov s method.
Finite 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 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 informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationAn Eulerian-Lagrangian WENO Scheme for Nonlinear Conservation Laws
An Eulerian-Lagrangian WENO Scheme for Nonlinear Conservation Laws Chieh-Sen Huang a,, Todd Arbogast b, a Department of Applied Mathematics, National Sun Yat-sen University, Kaohsiung 84, Taiwan, R.O.C.
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 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 informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationProject Description: MODELING FLOW IN VUGGY MEDIA. Todd Arbogast, Mathematics
Project Description: MODELING FLOW IN VUGGY MEDIA Todd Arbogast, Mathematics Steve Bryant, Petroleum & Geosystems Engineering Jim Jennings, Bureau of Economic Geology Charlie Kerans, Bureau of Economic
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
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 informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More information10.34 Numerical Methods Applied to Chemical Engineering. Quiz 2
10.34 Numerical Methods Applied to Chemical Engineering Quiz 2 This quiz consists of three problems worth 35, 35, and 30 points respectively. There are 4 pages in this quiz (including this cover page).
More information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationAN INTRODUCTION TO: DIFFERENTIAL EQUATIONS AND THEIR NUMERICAL APPROXIMATION. Todd Arbogast
AN INTRODUCTION TO: DIFFERENTIAL EQUATIONS AND THEIR NUMERICAL APPROXIMATION Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
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 information1 Introduction to MATLAB
L3 - December 015 Solving PDEs numerically (Reports due Thursday Dec 3rd, carolinemuller13@gmail.com) In this project, we will see various methods for solving Partial Differential Equations (PDEs) using
More informationMonotonicity Conditions for Discretization of Parabolic Conservation Laws. Hilde Kristine Hvidevold
Monotonicity Conditions for Discretization of Parabolic Conservation Laws Master of Science Thesis in Applied Mathematics Hilde Kristine Hvidevold Department of Mathematics University of Bergen June 2,
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationRT3D Rate-Limited Sorption Reaction
GMS TUTORIALS RT3D Rate-Limited Sorption Reaction This tutorial illustrates the steps involved in using GMS and RT3D to model sorption reactions under mass-transfer limited conditions. The flow model used
More informationICES REPORT A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams
ICS RPORT 15-17 July 2015 A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams by Omar Al-Hinai, Mary F. Wheeler, Ivan Yotov The Institute for Computational ngineering
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
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 informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More informationAcoustic Wave Equation
Acoustic Wave Equation Sjoerd de Ridder (most of the slides) & Biondo Biondi January 16 th 2011 Table of Topics Basic Acoustic Equations Wave Equation Finite Differences Finite Difference Solution Pseudospectral
More informationA NUMERICAL APPROXIMATION OF NONFICKIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXX, 1(21, pp. 75 84 Proceedings of Algoritmy 2 75 A NUMERICAL APPROXIMATION OF NONFICIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA R. E. EWING, Y. LIN and J. WANG
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationNumerical fractional flow modelling of inhomogeneous
Numerical fractional flow modelling of inhomogeneous air sparging E.F. Kaasschieter\ G.J. Mulder^ & J.D. van der Werff ten Bosch^ Department of Mathematics and Computing Science, hoven University of Technology,
More informationNumerical Hydraulics
ETHZ, Fall 017 Numerical Hydraulics Assignment 4 Numerical solution of 1D solute transport using Matlab http://www.bafg.de/ http://warholian.com Numerical Hydraulics Assignment 4 ETH 017 1 Introduction
More informationProblem Set 4 Issued: Wednesday, March 18, 2015 Due: Wednesday, April 8, 2015
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0139.9 NUMERICAL FLUID MECHANICS SPRING 015 Problem Set 4 Issued: Wednesday, March 18, 015 Due: Wednesday,
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
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 informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2017 Randall J. LeVeque Applied Mathematics University of Washington January 4, 2017 http://faculty.washington.edu/rjl/classes/am574w2017
More informationW.A. van El. Modeling Dispersion and Mixing in EOR Processes
-06-015 W.A. van El Modeling Dispersion and Mixing in EOR Processes Title : Modeling Dispersion and Mixing in EOR Processes Author(s) : W.A. van El Date : June 015 Professor(s) : Prof. dr. W.R. Rossen
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationNew Mexico Tech Hyd 510
Aside: The Total Differential Motivation: Where does the concept for exact ODEs come from? If a multivariate function u(x,y) has continuous partial derivatives, its differential, called a Total Derivative,
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.306 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Monday March 12, 2018. Turn it in (by 3PM) at the Math.
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationChp 4: Non-linear Conservation Laws; the Scalar Case. By Prof. Dinshaw S. Balsara
Chp 4: Non-linear Conservation Laws; the Scalar Case By Prof. Dinshaw S. Balsara 1 4.1) Introduction We have seen that monotonicity preserving reconstruction and iemann solvers are essential building blocks
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
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 informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationICES REPORT Accuracy of WENO and Adaptive Order WENO Reconstructions for Solving Conservation Laws
ICES REPORT 17-29 October 2017 Accuracy of WENO and Adaptive Order WENO Reconstructions for Solving Conservation Laws by Todd Arbogast, Chieh-Sen Huang, Xikai Zhao The Institute for Computational Engineering
More informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationAdvection, Conservation, Conserved Physical Quantities, Wave Equations
EP711 Supplementary Material Thursday, September 4, 2014 Advection, Conservation, Conserved Physical Quantities, Wave Equations Jonathan B. Snively!Embry-Riddle Aeronautical University Contents EP711 Supplementary
More information2. Conservation of Mass
2 Conservation of Mass The equation of mass conservation expresses a budget for the addition and removal of mass from a defined region of fluid Consider a fixed, non-deforming volume of fluid, V, called
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
More informationSolving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions
International Workshop on MeshFree Methods 3 1 Solving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions Leopold Vrankar (1), Goran Turk () and Franc Runovc (3) Abstract:
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
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 informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationThe Orchestra of Partial Differential Equations. Adam Larios
The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline
More informationSimple Lagrangian Model in 1 D
Chapter Simple Lagrangian Model in D. Ice Dynamics as Fluid Dynamics At the scale of interest to us here, it is not possible to resolve ea individual ice floe. We, thus, have to find an approximation whi
More informationIntroduction to PDEs and Numerical Methods: Exam 1
Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently
More informationSplitting methods in the design of coupled flow and mechanics simulators
Splitting methods in the design of coupled flow and mechanics simulators Sílvia Barbeiro CMUC, Department of Mathematics, University of Coimbra PhD Program in Mathematics Coimbra, November 17, 2010 Sílvia
More informationFully implicit higher-order schemes applied to polymer flooding. K.-A. Lie(SINTEF), T.S. Mykkeltvedt(IRIS), X. Raynaud (SINTEF)
Fully implicit higher-order schemes applied to polymer flooding K.-A. Lie(SINTEF), T.S. Mykkeltvedt(IRIS), X. Raynaud (SINTEF) Motivation Challenging and complex to simulate water-based EOR methods, unresolved
More informationA semi-lagrangian finite difference WENO scheme for scalar nonlinear conservation laws
A semi-lagrangian finite difference WENO scheme for scalar nonlinear conservation laws Chieh-Sen Huang a,, Todd Arbogast b,, Chen-Hui Hung c a Department of Applied Mathematics, National Sun Yat-sen University,
More informationFlow and Transport. c(s, t)s ds,
Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationScalar Conservation Laws and First Order Equations Introduction. Consider equations of the form. (1) u t + q(u) x =0, x R, t > 0.
Scalar Conservation Laws and First Order Equations Introduction. Consider equations of the form (1) u t + q(u) x =, x R, t >. In general, u = u(x, t) represents the density or the concentration of a physical
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationSTUDY OF DISCRETE DUALITY FINITE VOLUME SCHEMES FOR THE PEACEMAN MODEL
SUDY OF DISCREE DUALIY FINIE VOLUME SCHEMES FOR HE PEACEMAN MODEL C. CHAINAIS-HILLAIRE, S. KRELL, AND A. MOUON Abstract. In this paper, we are interested in the finite volume approximation of a system
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationAnalysis of Two-Grid Methods for Nonlinear Parabolic Equations by Expanded Mixed Finite Element Methods
Advances in Applied athematics and echanics Adv. Appl. ath. ech., Vol. 1, No. 6, pp. 830-844 DOI: 10.408/aamm.09-m09S09 December 009 Analysis of Two-Grid ethods for Nonlinear Parabolic Equations by Expanded
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
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 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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
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 informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
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 informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
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 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationExamination paper for TMA4195 Mathematical Modeling
Department of Mathematical Sciences Examination paper for TMA4195 Mathematical Modeling Academic contact during examination: Espen R. Jakobsen Phone: 73 59 35 12 Examination date: December 16, 2017 Examination
More informationCurvature and the continuity of optimal transportation maps
Curvature and the continuity of optimal transportation maps Young-Heon Kim and Robert J. McCann Department of Mathematics, University of Toronto June 23, 2007 Monge-Kantorovitch Problem Mass Transportation
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS
ON THE MULTIPOINT MIXED FINITE VOLUME METHODS ON QUADRILATERAL GRIDS VINCENT FONTAINE 1,2 AND ANIS YOUNES 2 1 Laboratoire de Physique des Bâtiments et des Systèmes, Université de La Réunion, 15 avenue
More informationMA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20
MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20 MA 201 (2016), PDE 2 / 20 Vibrating string and the wave equation Consider a stretched string of length
More informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More information