Advanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
|
|
- Egbert Bailey
- 5 years ago
- Views:
Transcription
1 Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg
2 Content Motivation and background R 3 T Numerical modelling advection advection + retardation + reaction advection + nonlinear retardation advective level set equation 2
3 Motivation and Background UG software toolbox - Unstructured Grids... to simplify the implementation of parallel adaptive multigrid method on unstructured grid for complex engineering applications. P. Bastian et. al
4 Motivation and Background Locally adapted multilevel grid conforming multilevel grid structure coarsening possible 4
5 Motivation and Background D 3 F application based on UG ( ) Distributed Density Driven Flow numerical modelling of gravity induced flows near saltdomes Frolkovic, De Schepper: Numerical modeling of convection dominated transport coupled with density driven flow in porous media; Advances in Water Resources, 2 5
6 Motivation and Background R 3 T application based on UG(999-24) Reaction Retardation Radionuclides Transport numerical modelling of radioactive contaminant transport F., Lampe, Wittum: r3t - software package for numerical simulations of radioactive contaminant transport in groundwater; WiR 25 6
7 R 3 T RRRT - Radionuclides Reactions (Decay) Np U Pu U 233 U 234 decay chains of up to 4 nuclides t C i = k λki C k λ ij C i 7
8 R 3 T RRRT - Transport Nuclides in flowing groundwater Np U Pu U U convection-dispersion-diffusion PDEs (up to 4) t C i + V C i D i ( V ) C i =... 8
9 R 3 T RRRT - Retardation of transport sorption Nuclides in flowing groundwater immobilization Np U Pu U U up to 2 additional ordinary differential equations t ( R i C i) + k i ( K i C i C i ad)
10 R 3 T Illustrative example (see video on my homepage) t ( R i C i) + k i ( K i C i C i ad) +...
11 R 3 T Linear case sorption Nuclides in flowing groundwater immobilization Np U Pu U U t ( R i C i) + k i ( K i C i C i ad) +...
12 R 3 T Nonlinear case sorption Nuclides in flowing groundwater immobilization Np U Pu 238 U U t ( R i (C)C i) + k i ( K i (C)C i C i ad)
13 R 3 T Sparsity of differential equations U Pu tu 238 t P 238 t U U 234 T T T U 238 P 238 U T := u D U 238 P 238 U 234 = 3
14 R 3 T Sparsity of discrete equations T ii T ij T ik T ji T jj T jk T ki T kj T kk T ii T ij T ik T ji T jj T jk T ki T kj T kk T ii T ij T ik T ji T jj T jk T ki T kj T kk U 238 i U 238 j U 238 k P 238 i P 238 j P 238 k U 234 i U 234 j U 234 k local stiff matrix for a triangle finite element 4
15 R 3 T Sparse matrix storage method (Neuss, 999) Jsparse:Daa=" * * ***"; Jsparse:Taa=" a a a"; U U Pu 5
16 R 3 T - numerical modelling Finite volume methods x j " ij! i T e x i " e ik x k Grid - unstructured numerical solution given pointwise gradient easily obtained from FE interpolation vertex-centred finite volume method (FVM) finite volume mesh dual to finite elements 6
17 R 3 T - numerical modelling Numerical solution piecewise linear, continuous c(t n, x) = c n i + T e c n (x x i ), x T e piecewise constant, discontinuous c(t n, x) = c n i, x Ω i piecewise linear reconstruction, discontinuous c(t n, x) = c n i + Ω i c n (x x i ), x Ω i 7
18 Motivation and Background Numerical modelling numerical algorithms fit analytical model preserving physical properties,... stable, consistent,... available, simple and good in general: unstructured grids robust for rough data,... (st order schemes) precise for smooth parts,... (2nd order schemes) 8
19 Advection-Diffusion-Dispersion Model equation t c + J =, J = V c D c 9
20 Advection-Diffusion-Dispersion Model equation t c + J =, J = V c D c FVM Ω i c n+ i = Ω i c n i t n J n+/2 ij exact integral formulation: t n+ Ω i c(t n+ ) = Ω i c(t n ) t n Ω i Ω j n J 2
21 Advection-Diffusion-Dispersion Model equation t c + J =, J = V c D c FVM Ω i c n+ i = Ω i c n i t n J n+/2 ij physical property - mass Ω i c n i : Ω i c(t n, x) dx 2
22 Advection-Diffusion-Dispersion Model equation t c + J =, J = V c D c FVM Ω i c n+ i = Ω i c n i t n J n+/2 ij physical property - conservation law J n+/2 ij = J n+/2 ji 22
23 Advection-Diffusion-Dispersion Model equation t c + J =, J = V c D c FVM Ω i c n+ i = Ω i c n i t n J n+/2 ij fully coupled implicit discretization J n+/2 ij = J ij (c n+ i, c n+ j, ) multigrid linear solver,... 23
24 Advection Model equation ( V c ) t c + = 24
25 Advection Motivation - exact simulation (see video on my homepage) 25
26 Advection Model equation ( V c ) t c + = 26
27 Advection Model equation ( V c ) t c + = FVM Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij 27
28 Advection Model equation ( V c ) t c + = FVM Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij physical property - mass Ω i c n i : Ω i c(t n, x) dx 28
29 Advection Model equation ( V c ) t c + = FVM Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij physical property - conservation law V ij = V ji, c n+/2 ij = c n+/2 ji 29
30 Advection Model equation ( V c ) t c + = FVM Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij physical property - characteristic curves c n+/2 ij :=? 3
31 Advection - st order scheme Model equation ( V c ) t c + = FVM Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij Piecewise constant numerical solution { c n i V ij > c n+/2 ij = c n j V ij < 3
32 Advection - st order scheme Model equation ( V c ) t c + = FVM Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij physical property - residence time = Ω i τ i max{, Vij } 32
33 Advection - st order scheme Model equation ( V c ) t c + = FVM Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij physical property - CFL condition t n τ i 33
34 Advection - st order scheme Courant number = Courant number > Courant number < 34
35 Advection - st order scheme Flux-based method of characteristics } DistributeMass(j, t, ˆτ, q) { t = t + τ j ; if (t t n+ ) then { } b j = b j + ˆτ q ; return; if (t + ˆτ > t n+ ) then { } j m = j ; b j = b j + (ˆτ (t n+ t )) q ; ˆτ = t n+ t ; for (j m 2 Λ out j m ) DistributeMass(j m 2, t, ˆτ, v j m j m 2 v jm q) ; return ; F.: Flux-based method of characteristics for transport in porous media; CVS, 22 35
36 Advection and Reaction and Retardation Courant number 5 Computation time 2.5 hours 36
37 Advection and Reaction and Retardation Courant number 5 Computation time.7 hours 37
38 Advection and Reaction and Retardation Example of 3 radionuclides R=, R2=3, R3=9, small physical dispersion V = (,), small dispersion, linear decay chain initially only st component non-zero 38
39 Advection and Reaction and Retardation Example of 3 radionuclides R=, R2=3, R3=9, small physical dispersion 2nd order Godunov method with many time steps 39
40 Advection and Reaction and Retardation Example of 3 radionuclides R=, R2=3, R3=9, small physical dispersion standard operator splitting method, 2 time steps 2nd order Godunov method with many time steps 4
41 Advection and Reaction and Retardation Example of 3 radionuclides R=, R2=3, R3=9, small physical dispersion flux-based method of characteristics F.: Flux-based method of characteristics for coupled system of transport equations in in porous media; CVS, 22 4
42 Advection Godunov method use exact solution of related simpler problem D Riemann s problem justified by numerical hyperbolic equations e.g., D advection => st order upwind m. High-resolution FVM piecewise linear numerical solution structured grid - Leveque 22 unstructured grid? (e.g., Sonar 993) 42
43 Advection and retardation Model equation ( V c ) t (Rφc) + = Fast sorption (equilibrium) R := + φ φ ρk linear case R = R(x) 43
44 Advection and retardation Example - Henry isotherm (see video on my homepage) R = 2 44
45 Advection and retardation Model equation ( V c ) t (Rφc) + = Fast sorption R := + φ φ ρk nonlinear case R = R(x, c) 45
46 Advection and retardation Example - Freundlich isotherm (see video on my homepage) R = + u p 46
47 Advection and retardation Nonlinear hyperbolic equation ( ) V c t θ + = θ = θ(c), c = θ (c) shocks.8 correct speed sharp also with diffusion rarefaction waves F., Kačur: Semi-analytical solutions of contaminant transport equation with nonlinear sorption in D; Comp. Geosciences, 26, to appear x 47
48 Advection and retardation Implementation (see video on my homepage) linear sorption 48 nonlinear sorption
49 Advection - st order method Trivial example t c + V c =, V c(, x) const c n+ i = c n i t n const 49
50 Advection - st order method Consistent for structured grid? t c + V c =, V c(, x) const Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij 5
51 Advection - st order method Consistent for structured grid! t c + V c =, V c(, x) const Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij 5
52 Advection - st order method Nonconsistent for unstructured grid! t c + V c =, V c(, x) const Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij 52
53 Advection - st order method (2x2) t = t =.9635 t = t =.5895 t =.7854 t = t =.78.5 t = t =
54 Advection Level set equation t c + V c =, V = 54
55 Advection Level set equation FVM t c + V c =, V = Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij 55
56 Advection Level set equation FVM t c + V c =, V = Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij physical property - value c n i : c(t n, x i ) 56
57 Advection Level set equation FVM t c + V c =, V = Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij physical property - characteristic curves c n+/2 ij := c(t n, X ij (t n )) 57
58 Advection - 2nd order scheme Level set equation FVM t c + V c =, V = Ω i c n+ i = Ω i c n i t n V ij c n+/2 ij physical property - characteristic curves c n+/2 ij := c(t n, X ij (t n )) c n+/2 ij := c n ij tn 2 V i c n i 58
59 Advection - st versus 2nd order method.9 t =.8.7,8.6,6.5 Y.4,4.3.2, ,2,4,6,8 X 59
60 Advection - st order method (2x2) t = t =.9635 t = t =.5895 t =.7854 t = t =.78.5 t = t =
61 Advection - 2nd order method (2x2) t = t =.9635 t = t =.5895 t =.7854 t = t =.78.5 t = t =
62 Advection - 2nd order method (2x2) t = t =.3927 t = t =.78 t =.578 t = t = t = t =
63 Advection Flux-based level set method (see video on my homepage) F., Mikula: High resolution flux-based level set method; 25 63
64 Nonlinear advective level set equation Example with topological changes (see video on my homepage) F., Mikula: Flux-based level set method: finite volume emthod for evolving interfaces; 22 64
65 Conclusions Numerical modelling numerical algorithms fit analytical model preserving physical properties,... stable, consistent,... available, simple and good in general unstructured grids robust for rough data,... (st order schemes) precise for smooth parts,... (2nd order schemes) 65
Advanced 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 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 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 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 informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
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 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 informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
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 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 informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
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 informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
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 informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
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 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 informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
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 informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationFinite volume method for two-phase flows using level set formulation
Finite volume method for two-phase flows using level set formulation Peter Frolkovič,a,1, Dmitry Logashenko b, Christian Wehner c, Gabriel Wittum c a Department of Mathematics and Descriptive Geometry,
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationStable Semi-Discrete Schemes for the 2D Incompressible Euler Equations
Stable Semi-Discrete Schemes for the D Incompressible Euler Equations Doron Levy Department of Mathematics Stanford University dlevy@math.stanford.edu http://math.stanford.edu/ dlevy Maryland, May 004
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 informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationMixed and Mixed Hybrid Finite Element Methods: Theory, Implementation and Applications
Mixed and Mixed Hybrid Finite Element Methods: Theory, Implementation and Applications Max-Planck Institute for Mathematics in the Sciences Leipzig, Germany mailto:fradu@mis.mpg.de Introduction Objective
More informationPalindromic Discontinuous Galerkin Method
Palindromic Discontinuous Galerkin Method David Coulette, Emmanuel Franck, Philippe Helluy, Michel Mehrenberger, Laurent Navoret To cite this version: David Coulette, Emmanuel Franck, Philippe Helluy,
More informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationThe Center for Astrophysical Thermonuclear Flashes. FLASH Hydrodynamics
The Center for Astrophysical Thermonuclear Flashes FLASH Hydrodynamics Jonathan Dursi (CITA), Alan Calder (FLASH) B. Fryxell, T. Linde, A. Mignone, G. Wiers Many others! Mar 23, 2005 An Advanced Simulation
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationMultigrid solvers for equations arising in implicit MHD simulations
Multigrid solvers for equations arising in implicit MHD simulations smoothing Finest Grid Mark F. Adams Department of Applied Physics & Applied Mathematics Columbia University Ravi Samtaney PPPL Achi Brandt
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
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 informationRESEARCH HIGHLIGHTS. WAF: Weighted Average Flux Method
RESEARCH HIGHLIGHTS (Last update: 3 rd April 2013) Here I briefly describe my contributions to research on numerical methods for hyperbolic balance laws that, in my view, have made an impact in the scientific
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Pseudo-Time Integration 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 10 Outline 1
More informationReactive transport in porous media: formulations, non-linear solvers and preconditioners
Reactive transport in porous media: formulations, non-linear solvers and preconditioners Michel Kern with L. Amir, B. Gueslin, A. Taakili Institut National de Recherche en Informatique et Automatique High
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
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 informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationarxiv: v2 [math.na] 23 Apr 2016
Improved ZZ A Posteriori Error Estimators for Diffusion Problems: Conforming Linear Elements arxiv:508.009v2 [math.na] 23 Apr 206 Zhiqiang Cai Cuiyu He Shun Zhang May 2, 208 Abstract. In [8], we introduced
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationModeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids
Modeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids Jørg Espen Aarnes and Vera Louise Hauge SINTEF ICT, Deptartment of Applied Mathematics Applied Mathematics
More informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
More informationOn the design of higher-order FEM satisfying the discrete maximum principle
On the design of higher-order FEM satisfying the discrete maximum principle Dmitri Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund, Germany
More informationSpace-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries
Space-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries Laurence Halpern, Caroline Japhet, and Jérémie Szeftel 3 LAGA, Université Paris XIII, Villetaneuse
More informationNumerical modelling of phase change processes in clouds. Challenges and Approaches. Martin Reitzle Bernard Weigand
Institute of Aerospace Thermodynamics Numerical modelling of phase change processes in clouds Challenges and Approaches Martin Reitzle Bernard Weigand Introduction Institute of Aerospace Thermodynamics
More informationOrder of Convergence of Second Order Schemes Based on the Minmod Limiter
Order of Convergence of Second Order Schemes Based on the Minmod Limiter Boan Popov and Ognian Trifonov July 5, 005 Abstract Many second order accurate non-oscillatory schemes are based on the Minmod limiter,
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 informationKrylov single-step implicit integration factor WENO methods for advection-diffusion-reaction equations
Accepted Manuscript Krylov single-step implicit integration factor WENO methods for advection diffusion reaction equations Tian Jiang, Yong-Tao Zhang PII: S0021-9991(16)00029-2 DOI: http://dx.doi.org/10.1016/j.jcp.2016.01.021
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 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 informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
More informationEQUADIFF 9. Mária Lukáčová-Medviďová Numerical solution of compressible flow. Terms of use:
EQUADIFF 9 Mária Lukáčová-Medviďová Numerical solution of compressible flow In: Zuzana Došlá and Jaromír Kuben and Jaromír Vosmanský (eds.): Proceedings of Equadiff 9, Conference on Differential Equations
More informationSunday September 28th. Time. 06:00 pm 09:00 pm Registration (CAAS) IP: Invited Presentation (IS or SP) 55 mn. CP: Contributed Presentation 25 mn
Sunday September 28th 06:00 pm 09:00 pm Registration (CAAS) IP: Invited Presentation (IS or SP) 55 mn CP: Contributed Presentation 25 mn Monday September 29th 07:30 am 08:30 am Registration 08:30 am 09:00
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
More informationFrancis X. Giraldo,
1 Time-Integrators Francis X. Giraldo, giraldo@nrlmry.navy.mil, www.nrlmry.navy.mil/~giraldo/projects/nseam.html 1.1 Introduction Roughly speaking, there are 2 classes of TIs: 1. EulerianMethods(fixed-frame-e.g.,arockatthebottomofaflowing
More informationNUMERICAL METHODS IN ASTROPHYSICS An Introduction
-1 Series in Astronomy and Astrophysics NUMERICAL METHODS IN ASTROPHYSICS An Introduction Peter Bodenheimer University of California Santa Cruz, USA Gregory P. Laughlin University of California Santa Cruz,
More informationA general well-balanced finite volume scheme for Euler equations with gravity
A general well-balanced finite volume scheme for Euler equations with gravity Jonas P. Berberich, Praveen Chandrashekar, Christian Klingenberg Abstract We present a second order well-balanced Godunov-type
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More informationAdaptive algorithm for saddle point problem for Phase Field model
Adaptive algorithm for saddle point problem for Phase Field model Jian Zhang Supercomputing Center, CNIC,CAS Collaborators: Qiang Du(PSU), Jingyan Zhang(PSU), Xiaoqiang Wang(FSU), Jiangwei Zhao(SCCAS),
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 informationComparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows
Comparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows Dr Argiris I. Delis Dr Ioannis K. Nikolos (TUC) Maria
More informationPast, present and space-time
Past, present and space-time Arnold Reusken Chair for Numerical Mathematics RWTH Aachen Utrecht, 12.11.2015 Reusken (RWTH Aachen) Past, present and space-time Utrecht, 12.11.2015 1 / 20 Outline Past. Past
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationCalculation of Sound Fields in Flowing Media Using CAPA and Diffpack
Calculation of Sound Fields in Flowing Media Using CAPA and Diffpack H. Landes 1, M. Kaltenbacher 2, W. Rathmann 3, F. Vogel 3 1 WisSoft, 2 Univ. Erlangen 3 inutech GmbH Outline Introduction Sound in Flowing
More informationOn efficient implicit upwind schemes
On efficient implicit upwind schemes Jostein R. Natvig Knut Andreas Lie June 30, 2008 Abstract Although many advanced methods have been devised for the hyperbolic transport problems that arise in reservoir
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationWell-Balanced Schemes for the Euler Equations with Gravity
Well-Balanced Schemes for the Euler Equations with Gravity Alina Chertock North Carolina State University chertock@math.ncsu.edu joint work with S. Cui, A. Kurganov, S.N. Özcan and E. Tadmor supported
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationSpatial discretization scheme for incompressible viscous flows
Spatial discretization scheme for incompressible viscous flows N. Kumar Supervisors: J.H.M. ten Thije Boonkkamp and B. Koren CASA-day 2015 1/29 Challenges in CFD Accuracy a primary concern with all CFD
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 informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More informationNumerical resolution of a two-component compressible fluid model with interfaces
Numerical resolution of a two-component compressible fluid model with interfaces Bruno Després and Frédéric Lagoutière February, 25 Abstract We study a totally conservative algorithm for moving interfaces
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
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 informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationNumerical schemes for short wave long wave interaction equations
Numerical schemes for short wave long wave interaction equations Paulo Amorim Mário Figueira CMAF - Université de Lisbonne LJLL - Séminaire Fluides Compréssibles, 29 novembre 21 Paulo Amorim (CMAF - U.
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 informationHIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS
June 6, 7 :7 WSPC/Lecture Notes Series: 9in x 6in chapter HIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS Chi-Wang Shu Division of Applied Mathematics, Brown University Providence,
More information