Finite Volume Method
|
|
- Vivien Goodman
- 5 years ago
- Views:
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 Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
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 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 informationAdjoint 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 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 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 informationNumerical 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 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 informationBurgers 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 informationAE/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 informationParabolic 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 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 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 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 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 informationZonal 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 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 informationINTRODUCTION 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 informationPoisson 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 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 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 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 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 informationIntroduction 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 information0.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 informationConstrained 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,
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 informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationComparison 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 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 informationFirst, 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 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 informationAE/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 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 informationWeighted 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 informationRelevant 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 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 informationFinite 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 informationAA214B: 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 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 information0.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 informationConservation 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 informationMAC2313 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 informationModelling 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 informationA 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 informationScientific 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 informationChapter 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 informationSistemas 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 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 informationModule 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 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 informationLocal 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 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 informationSpace-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 informationEntropy 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 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 informationFundamentals 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 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 informationIntroduction 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 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 informationSolving 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 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 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 informationNumerical 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 informationBasic 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 informationMAE 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 informationIMPLEMENTATION 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 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 informationExtension 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 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 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 informationUNIVERSITY 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 informationNumerical 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 informationAProofoftheStabilityoftheSpectral 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 informationA 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 informationTHE 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 informationWeighted 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 informationNon-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 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 informationTwo-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 informationVISCOUS 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 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 information1 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 informationV (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 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 informationA 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 informationConvergence 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 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 informationInternational 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 informationActive 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 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 information2.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 informationLecture 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 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 informationAntony 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 informationRhie-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 informationLecture: 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 informationLecture 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 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 informationCFD 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