Finite Volume Method

Size: px
Start display at page:

Download "Finite Volume Method"

Transcription

1 Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65

2 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 2 / 65

3 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 3 / 65

4 Divergence theorem In Cartesian coordinates, q = (u, v, w) q = u x + v y + w z Divergence theorem dv n V q dv = V V q ˆn ds Praveen. C (CTFD, NAL) FVM CMMACS 4 / 65

5 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 5 / 65

6 Conservation laws Basic laws of physics are conservation laws: mass, momentum, energy, charge Control volume n Streamlines Rate of change of φ in V = -(net flux across V ) φ dv = F ˆn ds t V V Praveen. C (CTFD, NAL) FVM CMMACS 6 / 65

7 Conservation laws If F is differentiable, use divergence theorem φ dv = t F dv or V [ ] φ t + F V dv = 0, V φ t + F = 0 In Cartesian coordinates, F = (f, g, h) φ t + f x + g y + h z = 0 for every control volume V Praveen. C (CTFD, NAL) FVM CMMACS 7 / 65

8 Differential and integral form Differential form φ t + F = 0 (conserved quantity) + divergence(flux) = 0 t Integral form φ dv = F ˆn ds t V V Praveen. C (CTFD, NAL) FVM CMMACS 8 / 65

9 Mass conservation equation φ = ρ, F = ρ q F ˆn ds = amount of mass flowing across ds per unit time ρ dv + ρ q ˆn ds = 0 t V V Using divergence theorem ρ dv + t or V V V [ ] ρ t + (ρ q) (ρ q) ds = 0 dv = 0 Praveen. C (CTFD, NAL) FVM CMMACS 9 / 65

10 Mass conservation equation Differential form ρ t + (ρ q) = 0 In Cartesian coordinates, q = (u, v, w) ρ t + x (ρu) + y (ρv) + z (ρw) = 0 Non-conservative form ρ t + u ρ x + v ρ y + w ρ ( u z + ρ x + v y + w ) = 0 z or in vector notation ρ t + q ρ + ρ q = 0 Praveen. C (CTFD, NAL) FVM CMMACS 10 / 65

11 Navier-Stokes equations Mass: Momentum: Energy: ρ t + (ρ q) = 0 t (ρ q) + (ρ q q) + p E t + (E + p) q = τ = ( q τ) + Q Praveen. C (CTFD, NAL) FVM CMMACS 11 / 65

12 Convection and diffusion Mass: Momentum: Energy: ρ t + (ρ q) = 0 t (ρ q) + (ρ q q) + p E t + (E + p) q = τ = ( q τ) + Q Praveen. C (CTFD, NAL) FVM CMMACS 12 / 65

13 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 13 / 65

14 Convection-diffusion equation Convection-diffusion equation u t + a u x = µ 2 u x 2 a,µ are constants, µ > 0 Conservation law form u t + f x f f c f d = 0 = f c + f d = au = µ u x Praveen. C (CTFD, NAL) FVM CMMACS 14 / 65

15 Convection equation Scalar convection equation u t + a u x = 0 Exact solution u(x, t + t) = u(x a t, t) Wave-like solution: convected at speed a Wave moves in a particular direction Solution can be discontinuous Hyperbolic equation Praveen. C (CTFD, NAL) FVM CMMACS 15 / 65

16 Scalar convection equation a > 0 t Praveen. C (CTFD, NAL) FVM CMMACS 16 / 65

17 Scalar convection equation a > 0 a t t t+ t Praveen. C (CTFD, NAL) FVM CMMACS 17 / 65

18 Convection equation Scalar convection equation u t + a u x = 0 Exact solution u(x, t + t) = u(x a t, t) Wave-like solution: convected at speed a Wave moves in a particular direction Solution can be discontinuous Hyperbolic equation Praveen. C (CTFD, NAL) FVM CMMACS 18 / 65

19 Diffusion equation Scalar diffusion equation u t = µ 2 u x 2 Solution u(x, t) e µk 2t f(x; k) Solution is smooth and decays with time Elliptic equation Praveen. C (CTFD, NAL) FVM CMMACS 19 / 65

20 Diffusion equation: Initial condition t Praveen. C (CTFD, NAL) FVM CMMACS 20 / 65

21 Diffusion equation t t+ t Praveen. C (CTFD, NAL) FVM CMMACS 21 / 65

22 Diffusion equation Scalar diffusion equation u t = µ 2 u x 2 Solution u(x, t) e µk 2t f(x; k) Solution is smooth and decays with time Elliptic equation Praveen. C (CTFD, NAL) FVM CMMACS 22 / 65

23 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 23 / 65

24 FVM in 1-D: Conservation law 1-D conservation law Initial condition u t + f(u) = 0, a < x < b, t > 0 x u(x, 0) = g(x) Boundary conditions at x = a and/or x = b Praveen. C (CTFD, NAL) FVM CMMACS 24 / 65

25 FVM in 1-D: Grid Divide computational domain into N cells a = x 1/2, x 1+1/2, x 2+1/2,..., x N 1/2, x N+1/2 = b h i i=1 i i=n x=a i 1/2 i+1/2 x=b Cell number i Cell size (width) C i = [x i 1/2, x i+1/2 ] h i = x i+1/2 x i 1/2 Praveen. C (CTFD, NAL) FVM CMMACS 25 / 65

26 FVM in 1-D Conservation law applied to cell C i xi+1/2 ( u t + f ) x x i 1/2 dx = 0 C i Cell-average value Integral form i 3/2 i 1/2 h i+1/2 i+3/2 i u i (t) = 1 xi+1/2 u(x, t) dx h i x i 1/2 h i du i dt + f(x i+1/2, t) f(x i 1/2, t) = 0 Praveen. C (CTFD, NAL) FVM CMMACS 26 / 65

27 FVM in 1-D i=1 i 1 i i+1 i=n x=a i 1/2 i+1/2 x=b Cell-average value is dicontinuous at the interface What is the flux? Numerical flux function f(x i+1/2, t) F i+1/2 = F(u i, u i+1 ) Praveen. C (CTFD, NAL) FVM CMMACS 27 / 65

28 FVM in 1-D Finite volume update equation (ODE) h i du i dt + F i+1/2 F i 1/2 = 0, i = 1,..., N Telescopic collapse of fluxes N ( ) du i h i + F dt i+1/2 F i 1/2 = 0 i=1 d dt N h i u i + F N+1/2 F 1/2 = 0 i=1 Rate of change of total u = -(Net flux across the domain) h i i=1 i i=n x=a i 1/2 i+1/2 x=b Praveen. C (CTFD, NAL) FVM CMMACS 28 / 65

29 Time integration At time t = 0 set the initial condition u o i := u i (0) = 1 xi+1/2 g(x) dx h i x i 1/2 Choose a time step t (stability and accuracy) or M Break time (0, T) into M intervals t = T M, tn = n t, n = 0, 1, 2,...,M For n = 0, 1, 2,...,M, compute u n i = u i (t n ) u(x i, t n ) u n+1 i = u n i t h i (F n i+1/2 F n i 1/2 ) Praveen. C (CTFD, NAL) FVM CMMACS 29 / 65

30 Numerical Flux Function Take average of the left and right cells U i+1 U i Central difference scheme du i dt i+1/2 F i+1/2 = 1 2 [f(u i) + f(u i+1 )] + f i+1 f i 1 h i = 0 Suitable for elliptic equations (diffusion phenomena) Unstable for hyperbolic equations (convection phenomena) Praveen. C (CTFD, NAL) FVM CMMACS 30 / 65

31 Convection equation: Centered flux Convection equation, a = constant u t + a u x = 0, Centered flux, F i+1/2 = (au i + au i+1 )/2 f(u) = au du i dt + a u i+1 u i 1 h i = 0 U i+1 U i i+1/2 Praveen. C (CTFD, NAL) FVM CMMACS 31 / 65

32 Convection equation: Upwind flux Upwind flux U i+1 U i or F i+1/2 = F i+1/2 = a + a 2 i+1/2 { au i, if a > 0 au i+1, if a < 0 u i + a a u i+1 2 Praveen. C (CTFD, NAL) FVM CMMACS 32 / 65

33 Convection equation: Upwind flux Upwind scheme or du i dt du i dt du i dt + a u i u i 1 h i = 0, if a > 0 + a u i+1 u i h i = 0, if a < 0 + a + a u i u i 1 2 h i + a a u i+1 u i = 0 2 h i Praveen. C (CTFD, NAL) FVM CMMACS 33 / 65

34 Diffusion equation Diffusion equation u t + f x = 0, f = µ u x No directional dependence; use centered formula On uniform grid F i+1/2 = µ u i+1 u i x i+1 x i du i dt = µ u i+1 2u i + u i 1 h 2 Praveen. C (CTFD, NAL) FVM CMMACS 34 / 65

35 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 35 / 65

36 Choice of numerical parameters What scheme to use? What should be the grid size h? What should be the time step t? Praveen. C (CTFD, NAL) FVM CMMACS 36 / 65

37 Stability: Convection equation Take a > 0 CFL number Centered flux Upwind flux Solution at t = 1 λ = a t h u n+1 i = u n i λ 2 un i+1 + λ 2 un i 1 u n+1 i = (1 λ)u n i + λu n i 1 Praveen. C (CTFD, NAL) FVM CMMACS 37 / 65

38 Convection equation: Initial condition 1 Initial condition Praveen. C (CTFD, NAL) FVM CMMACS 38 / 65

39 Convection equation: Centered scheme, λ = 0.8 8e+10 6e+10 Exact solution At t=1 4e+10 2e e+10-4e+10-6e+10-8e Praveen. C (CTFD, NAL) FVM CMMACS 39 / 65

40 Convection equation: Upwind scheme, λ = Exact solution At t= Praveen. C (CTFD, NAL) FVM CMMACS 40 / 65

41 Convection equation: Upwind scheme, λ = Exact solution At t= Praveen. C (CTFD, NAL) FVM CMMACS 41 / 65

42 Stability: Convection equation Centered flux: Unstable Upwind flux: conditionally stable λ 1 or t h a Praveen. C (CTFD, NAL) FVM CMMACS 42 / 65

43 Stability: Diffusion equation Central scheme Stability parameter Update equation u n+1 i u n i t = µ un i+1 2un i + u n i 1 h 2 ν = µ t h 2 u n+1 i = νu n i+1 + (1 2ν)un i + νu n i 1 Stability condition ν 1 2 or t h2 2µ Praveen. C (CTFD, NAL) FVM CMMACS 43 / 65

44 Stability: Convection-diffusion equation Convection-diffusion equation u t + f x = 0, f = au µ u x Upwind scheme for convective term; central scheme for diffusive term F i+1/2 = a + a 2 u i + a a u i+1 µ u i+1 u i 2 x i+1 x i Update equation u n+1 i = + + ( a + t + µ t ) h h ( 2 1 a t µ t h h ( 2 a t h + µ t h 2 u n i 1 ) u n i ) u n i+1 Praveen. C (CTFD, NAL) FVM CMMACS 44 / 65

45 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 45 / 65

46 Non-linear conservation law Inviscid Burgers equation Finite difference scheme u t + u u x = 0 du i dt + u i + u i u i u i 1 2 h + u i u i u i+1 u i 2 h = 0 Gives wrong results Not conservative; cannot be put in FV form h du i dt + F i+1/2 F i 1/2 = 0 Praveen. C (CTFD, NAL) FVM CMMACS 46 / 65

47 Non-linear conservation law Conservation law u t + ( ) u 2 = 0 or x 2 Upwind finite volume scheme u t + f x = 0, u2 f(u) = 2 with h du i dt + F i+1/2 F i 1/2 = 0 { 1 F i+1/2 = 2 u2 i if u i+1/ u2 i+1 if u i+1/2 < 0 Approximate interface value, for example u i+1/2 = 1 2 (u i + u i+1 ) Praveen. C (CTFD, NAL) FVM CMMACS 47 / 65

48 Transonic flow over NACA Praveen. C (CTFD, NAL) FVM CMMACS 48 / 65

49 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 49 / 65

50 Computational domain (Josy) Praveen. C (CTFD, NAL) FVM CMMACS 50 / 65

51 Grids and finite volumes Praveen. C (CTFD, NAL) FVM CMMACS 51 / 65

52 Structured grid Praveen. C (CTFD, NAL) FVM CMMACS 52 / 65

53 Unstructured grid Praveen. C (CTFD, NAL) FVM CMMACS 53 / 65

54 Unstructured hybrid grid Praveen. C (CTFD, NAL) FVM CMMACS 54 / 65

55 Unstructured hybrid grid Praveen. C (CTFD, NAL) FVM CMMACS 55 / 65

56 Cartesian grid Saras (Josy) Praveen. C (CTFD, NAL) FVM CMMACS 56 / 65

57 Block-structured grid Fighter aircraft (Nair) Praveen. C (CTFD, NAL) FVM CMMACS 57 / 65

58 Cell-centered and vertex-centered Praveen. C (CTFD, NAL) FVM CMMACS 58 / 65

59 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 59 / 65

60 FVM in 2-D Conservation law Integral form A i du i dt u t + f x + g y = 0 + (fn x + gn y ) ds = 0 C i C i C i Praveen. C (CTFD, NAL) FVM CMMACS 60 / 65

61 FVM in 2-D Approximate flux across the cell face nij C j C i F ij S ij (fn x + gn y ) ds Finite volume approximation A i du i dt + j N(i) F ij S ij = 0 Praveen. C (CTFD, NAL) FVM CMMACS 61 / 65

62 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 62 / 65

63 Summary Integral form of conservation law Satisfies conservation of mass, momentum, energy, etc. Can capture discontinuities like shocks Can be applied on any type of grid Useful for complex geometry Praveen. C (CTFD, NAL) FVM CMMACS 63 / 65

64 Outline 1 Divergence theorem 2 Conservation laws 3 Convection-diffusion equation 4 FVM in 1-D 5 Stability of numerical scheme 6 Non-linear conservation law 7 Grids and finite volumes 8 FVM in 2-D 9 Summary 10 Further reading Praveen. C (CTFD, NAL) FVM CMMACS 64 / 65

65 For further reading I Blazek J., Computational Fluid Dynamics: Principles and Applications, Elsevier, Hirsch Ch., Numerical Computation of Internal and External Flows, Vol. 1 and 2, Wiley. LeVeque R. J., Finite Volume Methods for Hyperbolic Equations, Cambridge University Press, Wesseling P., Principles of Computational Fluid Dynamics, Springer, Ferziger J. H. and Peric M., Computational methods for fluid dynamics, 3 rd edition, Springer, Praveen. C (CTFD, NAL) FVM CMMACS 65 / 65

Numerical Solution of Partial Differential Equations governing compressible flows

Numerical Solution of Partial Differential Equations governing compressible flows Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore

More information

The one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:

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

Numerical Solution Techniques in Mechanical and Aerospace Engineering

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

Adjoint code development and optimization using automatic differentiation (AD)

Adjoint code development and optimization using automatic differentiation (AD) Adjoint code development and optimization using automatic differentiation (AD) Praveen. C Computational and Theoretical Fluid Dynamics Division National Aerospace Laboratories Bangalore - 560 037 CTFD

More information

Finite volume method on unstructured grids

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

Conservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion

Conservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear

More information

Math background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids

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

Burgers equation - a first look at fluid mechanics and non-linear partial differential equations

Burgers equation - a first look at fluid mechanics and non-linear partial differential equations Burgers equation - a first look at fluid mechanics and non-linear partial differential equations In this assignment you will solve Burgers equation, which is useo model for example gas dynamics anraffic

More information

AE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

AE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept. AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says

More information

Parabolic Flow in Parallel Plate Channel ME 412 Project 4

Parabolic Flow in Parallel Plate Channel ME 412 Project 4 Parabolic Flow in Parallel Plate Channel ME 412 Project 4 Jingwei Zhu April 12, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of this project is to develop and apply a computer

More information

An Introduction to the Discontinuous Galerkin Method

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

Numerical Solutions to Partial Differential Equations

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

Finite Volume Schemes: an introduction

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 information

2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13

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

Zonal modelling approach in aerodynamic simulation

Zonal modelling approach in aerodynamic simulation Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability

More information

Getting started: CFD notation

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

INTRODUCTION TO PDEs

INTRODUCTION TO PDEs INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial

More information

Poisson Equation in 2D

Poisson Equation in 2D A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method

More information

Block-Structured Adaptive Mesh Refinement

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

The Orchestra of Partial Differential Equations. Adam Larios

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

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

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

Introduction to the Finite Volumes Method. Application to the Shallow Water Equations. Jaime Miguel Fe Marqués

Introduction to the Finite Volumes Method. Application to the Shallow Water Equations. Jaime Miguel Fe Marqués Introduction to the Finite Volumes Method. Application to the Shallow Water Equations. Jaime Miguel Fe Marqués Contents Preliminary considerations 3. Study of the movement of a fluid................ 3.

More information

0.2. CONSERVATION LAW FOR FLUID 9

0.2. CONSERVATION LAW FOR FLUID 9 0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is

More information

Constrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics

Constrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics Project work at the Department of Mathematics, TUHH Constrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics Katja Baumbach April 4, 005 Supervisor:

More information

[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,

[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of, 269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,

More information

Computational Fluid Dynamics-1(CFDI)

Computational Fluid Dynamics-1(CFDI) بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to

More information

Comparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models

Comparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models Comparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models Julian Koellermeier, Manuel Torrilhon October 12th, 2015 Chinese Academy of Sciences, Beijing Julian Koellermeier,

More information

Math 660-Lecture 23: Gudonov s method and some theories for FVM schemes

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

First, Second, and Third Order Finite-Volume Schemes for Diffusion

First, Second, and Third Order Finite-Volume Schemes for Diffusion First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.

More information

LECTURE # 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) 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 information

AE/ME 339. K. M. Isaac. 9/22/2005 Topic 6 FluidFlowEquations_Introduction. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

AE/ME 339. K. M. Isaac. 9/22/2005 Topic 6 FluidFlowEquations_Introduction. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept. AE/ME 339 Computational Fluid Dynamics (CFD) 1...in the phrase computational fluid dynamics the word computational is simply an adjective to fluid dynamics.... -John D. Anderson 2 1 Equations of Fluid

More information

Advanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg

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

Weighted ENO Schemes

Weighted ENO Schemes Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic

More information

Relevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]

Relevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE] Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential

More information

Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2

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

Finite volume approximation of the relativistic Burgers equation on a Schwarzschild (anti-)de Sitter spacetime

Finite volume approximation of the relativistic Burgers equation on a Schwarzschild (anti-)de Sitter spacetime Turkish Journal of Mathematics http:// journals. tubitak. gov. tr/ math/ Research Article Turk J Math 2017 41: 1027 1041 c TÜBİTAK doi:10.906/mat-1602-8 Finite volume approximation of the relativistic

More information

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29

More information

Wave propagation methods for hyperbolic problems on mapped grids

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

0.3.4 Burgers Equation and Nonlinear Wave

0.3.4 Burgers Equation and Nonlinear Wave 16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave

More information

Conservation Laws of Surfactant Transport Equations

Conservation Laws of Surfactant Transport Equations Conservation Laws of Surfactant Transport Equations Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada Winter 2011 CMS Meeting Dec. 10, 2011 A. Cheviakov

More information

MAC2313 Final A. (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative.

MAC2313 Final A. (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative. MAC2313 Final A (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative. ii. The vector field F = 5(x 2 + y 2 ) 3/2 x, y is radial. iii. All constant

More information

Modelling of turbulent flows: RANS and LES

Modelling of turbulent flows: RANS and LES Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012

More information

A finite-volume algorithm for all speed flows

A finite-volume algorithm for all speed flows A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,

More information

Scientific Computing I

Scientific Computing I Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline

More information

Chapter 2: Fluid Dynamics Review

Chapter 2: Fluid Dynamics Review 7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading

More information

Sistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas

Sistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas Sistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas Pedro González-Casanova Henríquez Unidad de Investigación en Cómputo Aplicado DGSCA, UNAM Elvira Torondel y Ricard Garrido

More information

Conservation Laws and Finite Volume Methods

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

Module 2: Governing Equations and Hypersonic Relations

Module 2: Governing Equations and Hypersonic Relations Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume

More information

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

Local discontinuous Galerkin methods for elliptic problems

Local discontinuous Galerkin methods for elliptic problems COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn

More information

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

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

Space-time XFEM for two-phase mass transport

Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,

More information

Entropy stable schemes for compressible flows on unstructured meshes

Entropy stable schemes for compressible flows on unstructured meshes Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/

More information

A Multi-Dimensional Limiter for Hybrid Grid

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

Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics

Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/

More information

PDE Solvers for Fluid Flow

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

Introduction to Physical Acoustics

Introduction to Physical Acoustics Introduction to Physical Acoustics Class webpage CMSC 828D: Algorithms and systems for capture and playback of spatial audio. www.umiacs.umd.edu/~ramani/cmsc828d_audio Send me a test email message with

More information

Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract

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

Solving the Euler Equations!

Solving the Euler Equations! http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p

More information

A Very Brief Introduction to Conservation Laws

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

Chapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation

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

Numerical Methods for Hyperbolic Conservation Laws Lecture 4

Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/

More information

Basic Aspects of Discretization

Basic Aspects of Discretization Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)

More information

MAE 101A. Homework 7 - Solutions 3/12/2018

MAE 101A. Homework 7 - Solutions 3/12/2018 MAE 101A Homework 7 - Solutions 3/12/2018 Munson 6.31: The stream function for a two-dimensional, nonviscous, incompressible flow field is given by the expression ψ = 2(x y) where the stream function has

More information

IMPLEMENTATION OF A PARALLEL AMG SOLVER

IMPLEMENTATION OF A PARALLEL AMG SOLVER IMPLEMENTATION OF A PARALLEL AMG SOLVER Tony Saad May 2005 http://tsaad.utsi.edu - tsaad@utsi.edu PLAN INTRODUCTION 2 min. MULTIGRID METHODS.. 3 min. PARALLEL IMPLEMENTATION PARTITIONING. 1 min. RENUMBERING...

More information

5. FVM discretization and Solution Procedure

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

Extension to moving grids

Extension to moving grids Extension to moving grids P. Lafon 1, F. Crouzet 2 & F. Daude 1 1 LaMSID - UMR EDF/CNRS 2832 2 EDF R&D, AMA April 3, 2008 1 Governing equations Physical coordinates Generalized coordinates Geometrical

More information

Équation de Burgers avec particule ponctuelle

Équation de Burgers avec particule ponctuelle Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin

More information

A Study on Numerical Solution to the Incompressible Navier-Stokes Equation

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

UNIVERSITY OF OSLO. Faculty of Mathematics and Natural Sciences. MEK4300/9300 Viscous flow og turbulence

UNIVERSITY OF OSLO. Faculty of Mathematics and Natural Sciences. MEK4300/9300 Viscous flow og turbulence UNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences Examination in: Day of examination: Friday 15. June 212 Examination hours: 9. 13. This problem set consists of 5 pages. Appendices: Permitted

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis

More information

AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy

AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University

More information

A Stable Spectral Difference Method for Triangles

A Stable Spectral Difference Method for Triangles A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna

More information

THE STOKES SYSTEM R.E. SHOWALTER

THE STOKES SYSTEM R.E. SHOWALTER THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes

More information

Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods

Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu

More information

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso

More information

A numerical study of SSP time integration methods for hyperbolic conservation laws

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

Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5

Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Jingwei Zhu May 14, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of

More information

VISCOUS FLUX LIMITERS

VISCOUS FLUX LIMITERS VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,

More information

Notes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow

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

1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.

1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω. 1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

Computational Analysis of an Imploding Gas:

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

A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations

A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm

More information

Convergence and Error Bound Analysis for the Space-Time CESE Method

Convergence and Error Bound Analysis for the Space-Time CESE Method Convergence and Error Bound Analysis for the Space-Time CESE Method Daoqi Yang, 1 Shengtao Yu, Jennifer Zhao 3 1 Department of Mathematics Wayne State University Detroit, MI 480 Department of Mechanics

More information

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

International Engineering Research Journal

International Engineering Research Journal Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule

More information

Active Flux for Advection Diffusion

Active Flux for Advection Diffusion Active Flux for Advection Diffusion A Miracle in CFD Hiroaki Nishikawa National Institute of Aerospace! NIA CFD Seminar! August 25, 2015 In collaboration with the University of Michigan Supported by NASA

More information

A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws

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

2.20 Fall 2018 Math Review

2.20 Fall 2018 Math Review 2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more

More information

Lecture 5: Kinetic theory of fluids

Lecture 5: Kinetic theory of fluids Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities

More information

ICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws

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

Antony Jameson. AIAA 18 th Computational Fluid Dynamics Conference

Antony Jameson. AIAA 18 th Computational Fluid Dynamics Conference Energy Estimates for Nonlinear Conservation Laws with Applications to Solutions of the Burgers Equation and One-Dimensional Viscous Flow in a Shock Tube by Central Difference Schemes Antony Jameson AIAA

More information

Rhie-Chow interpolation in OpenFOAM 1

Rhie-Chow interpolation in OpenFOAM 1 Rhie-Chow interpolation in OpenFOAM 1 Fabian Peng Kärrholm Department of Applied Mechanics Chalmers Univesity of Technology Göteborg, Sweden, 2006 1 Appendix from Numerical Modelling of Diesel Spray Injection

More information

Lecture: Wave-induced Momentum Fluxes: Radiation Stresses

Lecture: Wave-induced Momentum Fluxes: Radiation Stresses Chapter 4 Lecture: Wave-induced Momentum Fluxes: Radiation Stresses Here we derive the wave-induced depth-integrated momentum fluxes, otherwise known as the radiation stress tensor S. These are the 2nd-order

More information

Lecture Notes on Computational Fluid Dynamics

Lecture Notes on Computational Fluid Dynamics AMS 260, Computational Fluid Dynamics Lecture Notes on Computational Fluid Dynamics Instructor: Prof. Dongwook Lee (dlee79@ucsc.edu) MWF, 2:00 pm 3:10 pm at Kresge classroom 319 Winter, 2015 https://courses.soe.ucsc.edu/courses/ams260/winter15/01

More information

Introduction LECTURE 1

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

CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing. Lakhdar Remaki

CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing. Lakhdar Remaki CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing Lakhdar Remaki Outline What we doing in CFD? CFD in Industry Shock-filter model for mesh adaptation

More information