FUNDAMENTALS OF FINITE DIFFERENCE METHODS
|
|
- Carol Gallagher
- 6 years ago
- Views:
Transcription
1 FUNDAMENTALS OF FINITE DIFFERENCE METHODS By Deep Gupta 3 rd Year undergraduate, Mechanical Engg. Deptt., IIT Bombay Supervised by: Prof. Gautam Biswas, IIT Kanpur
2 Acknowledgements It has been a pleasure preparing this lecture. I take this opportunity to thank my supervisor Prof. Gautam Biswas (Professor,IIT Kanpur) for his enthusiastic support. A major chunk of this lecture is inspired by his lecture notes. Excerpts from Prof. Atul Sharma's (Asstt. Professor, IIT Bombay) lectures have been included with his kind permission and profuse thanks to him for the same. Last but not the least invaluable reference from J.W. Anderson's "Computational Fluid Dynamics" was indispensable in the completion of this lecture.
3 Contents The presentation discusses: Partial Differential eqns: classification and equations involved in fluid mechanics with boundary/initial conditions. Concept of discretization. Finite Difference Quotients, starting from Taylor s series. Finite Difference eqns.: concept, basic aspects like consistency and convergence. Explicit and Implicit methods of discretization: a comparison. Errors and stability analysis, including analysis of parabolic and hyperbolic equations. Fundamentals of Fluid Flow Modeling: Transportive, Conservative properties and Artificial Viscosity. Upwind Scheme of Discretization (First and Second). 3
4 Partial Differential Equations (Classification) Consider a second order differential eqn., φ φ φ φ φ φ =, Linear & non-linear eqns. Quasi-linear eqns. Homogeneous or non-homogeneous. Can be further classified as: B AC = parabolic 4 0, B 4 AC < 0, elliptic B 4 AC > 0, hyperbolic ( ) A B C D E F G x y x x y y x y 4
5 Partial Differential Equations (Examples) Unsteady N-S eqns are elliptic in space and parabolic in time. Laplace and Poisson eqns are elliptical in nature. φ φ 0, forlaplace' s + = x y Constt., forpoisson' s Heat Conduction eqn.(-d) is parabolic, φ t = B φ x 5
6 Partial Differential Equations (examples) Fluid flow problems have non-linear terms advection (for momentum eqn.) or convection (for energy eqn.). φ φ φ φ φ + u + v = B + S + t x y x y Parabolic in time and elliptical in space, turns hyperbolic for high speed flows. 6
7 Boundary/Initial Conditions No. of B.C.s=order of highest derivatives in differential equations. Unsteady problems require initial conditions. Spatial B.C.s in flow/ heat transfer problems (types): Dirchlet:- Φ=Φ (r) A Neumann:- Φ n=φ (r) A Mixed:- a(r)φ+b(r) Φ n=φ 3 (r) A 3 7
8 Discretization A closed form mathematical expression, such as a partial differential equation, is approximated by analogous expressions which prescribe values at only a finite no. of discrete points or volumes in domain. Analytical solution: continuous. Numerical solution: Discrete. Discretization Techniques Finite Difference Finite Volume Finite Element 8
9 Finite Difference Method (FDM) :Grid Generation x and y need not be uniform, only for simplifying the programming. Uniform spacing: Transformed Computational plane Non-Uniform spacing: Physical plane. 9
10 FDM: Finite Difference Quotient Taylor Series Expansion Example: f(x)=sin x x=0.; x=0.0; f(x+ x)=0.983; f(x)=0.95; π f =f =0.95(3.76%error f 3 = (0.775% error) f 4 =0.984(0.0% error) 0
11 FDM:Finite Difference Quotient (Cont.) Taylor Series expansion I Order Forward Difference
12 FDM:Finite Difference Quotient (Cont.)
13 FDM:Finite Difference Quotient (Cont.) 3
14 FDM:Finite Difference Quotient (Cont.) II Order Central Difference for mixed Derivative xy φ + i+, j+ i, j i+, j i, j+ ij, φ φ φ φ ( ) ( ) = +Ο x, y 4 xy Δ Δ ΔΔ Higher-order-accurate difference equations Disadvantages: Require more grid points, thus require more computer time for each time wise or spatial step. Advantages: Require a smaller no. of total grid points in a flow solution to obtain overall accuracy. 4
15 FDM: Polynomial approach One Sided Finite Differences, other side points not known. Not restricted to be used at boundaries only. u y For example, assume, u = a+ by+ cy After subsequent solving, get u 3u + 4u u = +Ο Δ y Δy Similarly, for third order accuracy, ij, ( y) 3 u + 8u 9u + u ij, ij, + ij, + ij, + 3 = +Ο Δ 6Δy ( y) 3 5
16 Finite Difference Equations Concept: All partial derivatives finite-difference quotients. Consider unsteady,-d heat conduction: T = α t T x Time t is a marching variable. Using FTCS (Forward Time Central Space ) method of discretization: T T α = 0= t x T PartialDifferentialEq. n+ i DifferenceEq. ( n n n T ) i+ Ti + Ti n T α i TE.. + Δt Δx ( ) 6
17 Finite Difference Equations: Consistency and Convergence Consistency: As the mesh is refined, lim PDE FDE = lim ( TE) = 0 ( ) mesh 0 mesh 0 Depends on Differencing Scheme, for eg., DuFort-Frankel scheme is not consistent for D unsteady state heat conduction eq. Convergence: Approximate sol. approaches the exact one for each value of the independent variable as grid spacing tends to zero, i.e., n u = u x, t asδx, Δt 0 ( ) i i n 7
18 Explicit and Implicit Methods Explicit method: Consider -D unsteady state heat conduction equation, T = α t F.D.E. is (FTCS): T n+ i n Ti = Δt α T x ( n n n T ) i+ Ti + Ti Δx ( ) T n i The eq. can be written as: n n n α Δt T n i+ Ti + T + i = Ti + ( Δx) Being Parabolic,this eq. lends itself to a marching solution. ( )( ) 8
19 Explicit and Implicit Methods (Cont.) Implicit Method: Considering the same -D heat conduction eq., but with different discretization, This method is known as Crank-Nicolson implicit scheme T n+ i ( n n n n n n T ) i+ + Ti + Ti Ti + Ti + Ti n T α i = Δt Δx ( ) Obtained by expressing the spatial differences in terms of averages between n & n+ time levels. 9
20 Explicit and Implicit methods (Cont.) Implicit Method (cont.): The eq. can be written as, n n r ( α Δt + n+ n n+ n n+ n T ) where r= i Ti = Ti+ + Ti+ Ti Ti + Ti + Ti Δx or, ( ) ( ) n+ + r n+ n+ n r n n Ti + Ti Ti+ = Ti + Ti + rti+ r r applied at all grid pts.,i.e., from i= to i=k+ to get the tridiagonal matrix,which can be solved. n+ n B() 0 0 K 0 T ( C() + A) n n B() K T3 C() n+ n 0 B(3) K 0 T4 = C(3) M M M M M M M n+ n K Bk ( ) T k ( Ck ( ) + D) where, A & D= Boundary conditions C(k)=R.H.S. of the eq. at i=k+ 0
21 Explicit and Implicit Methods (Cont.) Implicit Method for -D conduction: T T T P.D.E.= = α + t x y Applying Crank-NicolsonScheme, F.D.E.= n+ n Ti, j Ti, j α ( δ )( n n ) x δ + = + y Ti, j + Ti, j Δt where, The system is not n n n T, n i+ j T i, j + T i, j tridiagonal (5 unknowns) δ x T i, j = Δ x δ ( ) T T + T n n n n i, j + i, j i, j y T i, j = ( Δ y ) T, T, T, T & T n+ n+ n+ n+ n+ ij, i+, j i, j ij, + ij,
22 Explicit and Implicit Methods (Cont.)
23 Explicit and Implicit Methods (Cont.) 3
24 Explicit and Implicit Methods: Comparison Implicit Method Advantage: Fewer time steps for calculations over a time interval. Disadvantage: - Complicated to set up the program. -Computer time required each step is much high. -For larger Δt, truncation error is high. Explicit method Advantage: Simple solution algorithm. Disadvantage: Requires many time steps to carry out the calculations over a given time interval, due to restrictions on Δt imposed by stability constraints. 4
25 Explicit and Implicit methods: The Choice Time marching broadly has purposes: To obtain a steady state solution by calculating the flow in steps of time, until a final steady-state flow is approached. No need for timewise accuracy approach the final steady state IMPLICIT Preferred. To obtain an accurate timewise solution of an inherently unsteady flow, Timewise accuracy is absolutely necessary less Truncation Error EXPLICIT Preferred. 5
26 Errors and Stability Analysis Errors: Given, A=analytical sol of P.D.E. D=exact sol. of F.D.E. N=numerical sol. from a real computer with finite accuracy. Then, Discretization Error=A-D =Truncation Error+error introduced due to treatment of boundary condition. Round-off Error=ε=N-D Numerical error introduced for a repetitive no. of calculations. or, N=D+ε 6
27 Error and Stability Analysis (Cont.) Considering the F.D.E. of -D heat conduction, Numerical sol. N must satisfy the above eq., thus, But D is the exact sol., therefore, Subtracting the above two,ε also satisfies the F.D.E. 7
28 Errors and Stability Analysis (Cont.) As solution progresses from step n to n+,it is: Stable: If the ε i s shrink, or remain the same. Unstable: If the ε i s grow larger. Condition for stability: 8
29 Error and Stability Analysis (Contd.) εi s are of the form: This random variation of εi s can be expressed as a Fourier series: Sequential Addition of sine and cosine functions with sequentially increasing wave no. where k m is wave no 9
30 Error and Stability Analysis (Cont.) Wave Number, k m : Consider the sine function, m is related with the no. of wave fitted in a given interval. Assuming the time-wise distribution is exponential in t, Since the F.D.E. is linear, can deal with one term only. 30
31 Errors and Stability Analysis (Cont.) Substituting the value of ε m (x,t) in the F.D.E., 3
32 Errors and Stability Analysis (Cont.) Substituting value of ε m (x,t) in the condition for stability, G is called Amplification Factor. can be possibilities: Called Von Neumann stability analysis. 3
33 Stability of Hyperbolic equations Considering the first order wave eq. u u + c = 0 t x Applying Von Neumann stability analysis to its FTCS discretization shows unconditional unstability. So Lax method of discretization is used, where u(t) is represented by an average value between grid pts. i+ and i-,in time derivative, i.e., ut () = u + u ( n n ) i+ i 33
34 Stability of Hyperbolic Equations (Contd.) Then, u = t u u + u Δt ( + ) n+ n n i i i n n n n u i u n+ ui+ + ui Δt + i ui = c Δx Assuming error of the same form, ( ) at ε m x t = e e, m Amplification factor becomes, ( ) ik x ( ) sin ( ) e at Δ = cos kδ x ic kδx m From stability requirement, finally get, C Δ t = c Δ x C is Courant Number. The condition is CFL (Courant-Freidrichs- Lewy) condition. m 34
35 Stability of hyperbolic Equations (Contd.) Consider second order eq., u u = c t x If v=δu/δt and w=c(δu/δx),then the system of eqns can be written as, u t which is a first order eqn. v 0 c where & A u= w u + = x [ A] 0 = c 0 Eigenvalues of matrix A : det[a-λi]=0, or λ = c 0 Roots are λ =+c & λ =-c, representing traveling waves with speeds c & -c. The characteristic lines are ct (right-running) x = ct (left-running) 35
36 Stability of Hyperbolic Equations (Contd.) Courant No. (C)< Numerical Domain includes all the Analytical Domain, hence STABLE Courant no. (C)> Numerical Domain does not include all the Analytical Domain, hence UNSTABLE 36
37 Stability Analysis-Real Perspective From the definition of Round Off error ε, an incorrect impression is formed that no instabilities would be there, if perfect computer is used with no round-off error. General concept of numerical stability is, in reality,based on the timewise behavior of the solution itself. Stability doesn t inherently depend on behavior of roundoff error per se. General Von Neumann analysis is, where solution itself is written as a Fourier series, i. e., Ux = ik x Ve m where Vm m ( ) m is the amplitude of the mth harmonic of the solution. 37
38 Fluid Flow Modeling Fluid Flow problems complex. G.E. s form a non-linear system. Model eqn. must have a convective, diffusive and time-dependent term. Burger s eqn satisfies all requirements, + u = υ t x x ζ ζ ζ For inviscid flow, it is replaced by Euler s eqn., ζ ζ + u = 0 t x 38
39 Fluid Flow Modeling :Conservative Property Consider the vorticity transport eqn., ω t = V. ω + υ ω ( ) Integrating over fixed region R, finally get, ω d R= V ω. nda + υ ω. nda ( ) ( ) t R A0 A0 R This implies that the time rate of accumulation of ω in equals net advective flux rate of ω across A 0 into Rand net diffusive flux rate of ω across A 0 into R 39
40 Conservative Property (Contd.) Concept of conservative property is to maintain this integral relation in finite difference representation. Considering the conservative form (which follows the said property) of Burger s eqn., i. e., ω = t x ( uω ) Using FTCS method, the F.D.E. is, n+ n n n n n ωi ωi ui+ ωi+ ui ωi = Δt Δx 40
41 Conservative Property (Contd.) Consider a region R from i=i to i=i.evaluating Δ t i = I the integral ω Δ x as, i = I I I I n+ n ( ) n ( ) n ωi Δ x ωiδ x = uω uω i i+ Δt i= I i= I i= I R 4
42 Conservative Property (Contd.) Summation of R.H.S. finally gives, ( uω ) ( uω ) = I I + Rate of accumulation of ω i in R is equal to net advective flux rate across the boundary of running from i=i to I. F.D.E. has preserved integral conservation relations, hence possesses conservative property. 4
43 Conservative Property (Contd.) Conservative Property depends on the form of continuum eqn. used. For e.g., consider non-conservative form : ω ω = u t x After using FTCS technique and then applying integration as earlier, finally yields, I n n n n The F.D.E. fails to preserve the integral Gauss-divergence property, i.e., the conservative property. i= I u ω u ω i i i i+ 43
44 The Upwind Scheme and Transportive Property FTCS to Burger s eqn, found to be unconditionally unstable. Thus, Upwind Scheme of discretization is used, n+ n n n ζi ζi uζi uζi = + viscous terms u>0 Δt Δx n+ n n n ζi ζi uζi + uζi = + viscous terms u<0 Δt Δx Transportive Property: If the effect of perturbation is convected only in the direction of velocity. F.D.E. formed by FTCS method violates the said property. 44
45 Transportive Property (Contd.) Consider a perturbation ε m =δ in ζ. For u>0, ε m =δ (perturbation at m th space location), all other ε=0. Now, checking the upwind scheme, at the down stream location (m+) n+ n ζm + ζm+ 0 uδ uδ = =+ Δt Δx Δx which is acceptable. 45
46 Transportive Property (Contd.) At point m of disturbance, i.e., perturbation is being transported out. At upstream (m-), ζ ζ n + n m ζ m u δ 0 u = = δ Δt Δx Δx n+ m n ζ m 0 0 = = Δt Δx no perturbation effect is carried upstream. Upwind Method maintains unidirectionality. 0 46
47 Upwind Scheme and Artificial Viscosity From Taylor Series, Substituting in F.D.E. obtained by upwind scheme, or, ζ ζ n + i n i ± ( Δt) ( Δ t ) n ζ ζ = ζ i + Δ t + + t t ( Δ x ) n n ζ ζ ζ i x x i x n i = ± Δ + + ( Δx) ζ ζ 3 u ζ ζ 3 Δ t + +Ο ( Δ t) = Δx +Ο ( Δx) Δt t t Δx x x +[diffusive terms] ζ ζ uδt ζ ζ = u + uδx + ν +Ο Δx t x x Δ x x n i n i ( ) L L 3 47
48 Artificial Viscosity (Contd.) The eqn. may be written as, ζ = u ζ + ν ζ e + ν ζ + higher-order terms t x x x where, ν e = ( ) and C (Courant number)=u t/ x uδx C Artificial Viscosity depends on discretization procedure. For steady state, i.e., ζ/ t=0, ν e =u x/. For a -D convective-diffusive eqn., following the similar steps, the eqn. obtained, ζ u ζ v ζ ( ) ζ = + νex + ν + ( ν ) ey + ν ζ Where, t x y x y ν ex = ( ), uδx C x with C x =u t/ x, C y =v t/ x ν ey = uδy( C y ) ; 48
49 Second Upwind Differencing (Hybrid Scheme) Usage of higher order upwind method of differencing, or method which improves the accuracy, is favored. According to the method of second upwind differencing, u ζ u ζ u ζ = x Δx ( ) i, j R R L L where, u R =(u i,j +u i+,j )/ ; u L =(u i,j +u i-,j )/ ; and, ζ R =ζ i,j for u R >0; ζ R =ζ i+,j for u R <0; ζ L =ζ i-,j for u L >0; ζ L =ζ i,j for u L <0; 49
50 Hybrid Scheme (Contd.) For u R >0 and u L >0,get, For unsteady x-direction x momentum eqn., i.e., ζ=u, and introducing a factor η, which can express eqn. as a weighted average of central and upwind differencing, u x i, j where 0<η<. ( uζ ) u, u+, u, u, ij+ i j ij+ i j = ζij, ζi, j x Δx ( ui, j ui+, j) η ( ui, j ui+, j) = + + 4Δx ( ) ( u ) i, j ui, j η ui, j ui, j + The accuracy of the scheme can always be increased by a suitable adjustment of η value. 50
51 Summary The Presentation showed: Finite Difference methods is one of the important discretization technique in CFD, and is based on replacing partial derivatives with difference quotients. Consistency and convergence of a F.D.E. depend on the Truncation Error & differencing scheme used. Explicit and Implicit methods are different CFD techniques, and their use depend on the type of problem being solved. The exact form of stability criterion depends on the form of differential eqn, e.g.,von Neumann stability criteria for parabolic and CFL for hyperbolic eqns Fluid flow can be modeled using a simpler eqn, called Burger s eqn. Conservative property depends on form of continuum eqn. used. Space centred differences are more accurate than upwind differences, as indicated by Taylor s series, the whole system is not more accurate if Transportive Property is looked upon. So, a combination is maintained using second upwind differencing. 5
52 Thank You for your kind Attention And Welcome for Any Questions, Comments or Suggestions. 5
Basics of Discretization Methods
Basics of Discretization Methods In the finite difference approach, the continuous problem domain is discretized, so that the dependent variables are considered to exist only at discrete points. Derivatives
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 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 informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationModule 3: BASICS OF CFD. Part A: Finite Difference Methods
Module 3: BASICS OF CFD Part A: Finite Difference Methods THE CFD APPROACH Assembling the governing equations Identifying flow domain and boundary conditions Geometrical discretization of flow domain Discretization
More informationLecture 17: Initial value problems
Lecture 17: Initial value problems Let s start with initial value problems, and consider numerical solution to the simplest PDE we can think of u/ t + c u/ x = 0 (with u a scalar) for which the solution
More informationFinite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More 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 informationq t = F q x. (1) is a flux of q due to diffusion. Although very complex parameterizations for F q
! Revised Tuesday, December 8, 015! 1 Chapter 7: Diffusion Copyright 015, David A. Randall 7.1! Introduction Diffusion is a macroscopic statistical description of microscopic advection. Here microscopic
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 1: Finite Difference Method
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture1/1_1.htm 1 of 1 6/19/2012 4:29 PM The Lecture deals with: Classification of Partial Differential Equations Boundary and Initial Conditions Finite Differences
More informationLecture 4.5 Schemes for Parabolic Type Equations
Lecture 4.5 Schemes for Parabolic Type Equations 1 Difference Schemes for Parabolic Equations One-dimensional problems: Consider the unsteady diffusion problem (parabolic in nature) in a thin wire governed
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationNumerical methods Revised March 2001
Revised March 00 By R. W. Riddaway (revised by M. Hortal) Table of contents. Some introductory ideas. Introduction. Classification of PDE's.3 Existence and uniqueness.4 Discretization.5 Convergence, consistency
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
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 informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
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 informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationChapter 4. Nonlinear Hyperbolic Problems
Chapter 4. Nonlinear Hyperbolic Problems 4.1. Introduction Reading: Durran sections 3.5-3.6. Mesinger and Arakawa (1976) Chapter 3 sections 6-7. Supplementary reading: Tannehill et al sections 4.4 and
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 informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 7:
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture7/7_1.htm 1 of 1 6/20/2012 12:26 PM The Lecture deals with: Errors and Stability Analysis file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture7/7_2.htm
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 informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationDiscretization of Convection Diffusion type equation
Discretization of Convection Diffusion type equation 10 th Indo German Winter Academy 2011 By, Rajesh Sridhar, Indian Institute of Technology Madras Guides: Prof. Vivek V. Buwa Prof. Suman Chakraborty
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More information2.2. Methods for Obtaining FD Expressions. There are several methods, and we will look at a few:
.. Methods for Obtaining FD Expressions There are several methods, and we will look at a few: ) Taylor series expansion the most common, but purely mathematical. ) Polynomial fitting or interpolation the
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationFinite volume method for CFD
Finite volume method for CFD Indo-German Winter Academy-2007 Ankit Khandelwal B-tech III year, Civil Engineering IIT Roorkee Course #2 (Numerical methods and simulation of engineering Problems) Mentor:
More informationThe Finite Difference Method
Chapter 5. The Finite Difference Method This chapter derives the finite difference equations that are used in the conduction analyses in the next chapter and the techniques that are used to overcome computational
More informationCharacteristic finite-difference solution Stability of C C (CDS in time/space, explicit): Example: Effective numerical wave numbers and dispersion
Spring 015 Lecture 14 REVIEW Lecture 13: Stability: Von Neumann Ex.: 1st order linear convection/wave eqn., F-B scheme Hyperbolic PDEs and Stability nd order wave equation and waves on a string Characteristic
More informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
More informationChapter 5 Types of Governing Equations. Chapter 5: Governing Equations
Chapter 5 Types of Governing Equations Types of Governing Equations (1) Physical Classification-1 Equilibrium problems: (1) They are problems in which a solution of a given PDE is desired in a closed domain
More informationNumerical Analysis and Computer Science DN2255 Spring 2009 p. 1(8)
Numerical Analysis and Computer Science DN2255 Spring 2009 p. 1(8) Contents Important concepts, definitions, etc...2 Exact solutions of some differential equations...3 Estimates of solutions to differential
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationNumerical Solution of partial differential equations
G. D. SMITH Brunei University Numerical Solution of partial differential equations FINITE DIFFERENCE METHODS THIRD EDITION CLARENDON PRESS OXFORD Contents NOTATION 1. INTRODUCTION AND FINITE-DIFFERENCE
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More information2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
3 Quantitative Properties of Finite Difference Schemes 31 Consistency, Convergence and Stability of FD schemes Reading: Tannehill et al Sections 333 and 334 Three important properties of FD schemes: Consistency
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More 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 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 informationChapter 2 Finite-Difference Discretization of the Advection-Diffusion Equation
Chapter Finite-Difference Discretization of the Advection-Diffusion Equation. Introduction Finite-difference methods are numerical methods that find solutions to differential equations using approximate
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationIntroduction to PDEs and Numerical Methods Tutorial 4. Finite difference methods stability, concsistency, convergence
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 4. Finite difference methods stability, concsistency, convergence Dr. Noemi Friedman,
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
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 informationClassification of partial differential equations and their solution characteristics
9 TH INDO GERMAN WINTER ACADEMY 2010 Classification of partial differential equations and their solution characteristics By Ankita Bhutani IIT Roorkee Tutors: Prof. V. Buwa Prof. S. V. R. Rao Prof. U.
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationECE539 - Advanced Theory of Semiconductors and Semiconductor Devices. Numerical Methods and Simulation / Umberto Ravaioli
ECE539 - Advanced Theory of Semiconductors and Semiconductor Devices 1 General concepts Numerical Methods and Simulation / Umberto Ravaioli Introduction to the Numerical Solution of Partial Differential
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 informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
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 informationFINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES
FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES SHRUTI JAIN B.Tech III Year, Electronics and Communication IIT Roorkee Tutors: Professor G. Biswas Professor S. Chakraborty ACKNOWLEDGMENTS I would like
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationFinite Differences: Consistency, Stability and Convergence
Finite Differences: Consistency, Stability and Convergence Varun Shankar March, 06 Introduction Now that we have tackled our first space-time PDE, we will take a quick detour from presenting new FD methods,
More informationAnswers to Exercises Computational Fluid Dynamics
Answers to Exercises Computational Fluid Dynamics Exercise - Artificial diffusion upwind computations.9 k=. exact.8 k=..7 k=..6 k=..5.4.3.2...2.3.4.5.6.7.8.9 x For k =.,. and., and for N = 2, the discrete
More informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationApproximations of diffusions. Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 3b Approximations of diffusions Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.1 04/10/2018 1 Learning objectives Become aware of the existence of stability conditions for the
More informationSOE3213/4: CFD Lecture 3
CFD { SOE323/4: CFD Lecture 3 @u x @t @u y @t @u z @t r:u = 0 () + r:(uu x ) = + r:(uu y ) = + r:(uu z ) = @x @y @z + r 2 u x (2) + r 2 u y (3) + r 2 u z (4) Transport equation form, with source @x Two
More informationarxiv: v1 [physics.comp-ph] 22 Feb 2013
Numerical Methods and Causality in Physics Muhammad Adeel Ajaib 1 University of Delaware, Newark, DE 19716, USA arxiv:1302.5601v1 [physics.comp-ph] 22 Feb 2013 Abstract We discuss physical implications
More information8. Introduction to Computational Fluid Dynamics
8. Introduction to Computational Fluid Dynamics We have been using the idea of distributions of singularities on surfaces to study the aerodynamics of airfoils and wings. This approach was very powerful,
More informationNumerical Solutions of the Burgers System in Two Dimensions under Varied Initial and Boundary Conditions
Applied Mathematical Sciences, Vol. 6, 22, no. 3, 563-565 Numerical Solutions of the Burgers System in Two Dimensions under Varied Initial and Boundary Conditions M. C. Kweyu, W. A. Manyonge 2, A. Koross
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
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 informationDissipation and Dispersion
Consider the problem with periodic boundary conditions Dissipation and Dispersion u t = au x 0 < x < 1, t > 0 u 0 = sin 40 πx u(0, t) = u(1, t) t > 0 If a > 0 then the wave is moving to the left and if
More informationFundamentals Physics
Fundamentals Physics And Differential Equations 1 Dynamics Dynamics of a material point Ideal case, but often sufficient Dynamics of a solid Including rotation, torques 2 Position, Velocity, Acceleration
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More informationProblem Set 4 Issued: Wednesday, March 18, 2015 Due: Wednesday, April 8, 2015
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0139.9 NUMERICAL FLUID MECHANICS SPRING 015 Problem Set 4 Issued: Wednesday, March 18, 015 Due: Wednesday,
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
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 informationNumerical Hydraulics
ETHZ, Fall 017 Numerical Hydraulics Assignment 4 Numerical solution of 1D solute transport using Matlab http://www.bafg.de/ http://warholian.com Numerical Hydraulics Assignment 4 ETH 017 1 Introduction
More informationNumerical methods for complex systems. S. V. Gurevich
Numerical methods for complex systems S. V. Gurevich December 16, 2008 2 Part I Partial Differential Equations 3 Chapter 1 Intorduction 1.1 Definition, Notation and Classification A differential equation
More informationNumerical Methods for Engineers and Scientists
Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL
More informationIntroduction to PDEs and Numerical Methods: Exam 1
Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationLecture 16: Relaxation methods
Lecture 16: Relaxation methods Clever technique which begins with a first guess of the trajectory across the entire interval Break the interval into M small steps: x 1 =0, x 2,..x M =L Form a grid of points,
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 informationComputational Fluid Dynamics
Computational Fluid Dynamics Dr.Eng. Reima Iwatsu Phone: 0355 69 4875 e-mail: iwatsu@las.tu-cottbus.de NACO Building Room 53-107 Time Summer Term Lecture: Tuesday 7:30-9:00 (every two weeks) LG4/310 Exercise:
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 Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More information3.3. Phase and Amplitude Errors of 1-D Advection Equation
3.3. Phase and Amplitude Errors of 1-D Advection Equation Reading: Duran section 2.4.2. Tannehill et al section 4.1.2. The following example F.D. solutions of a 1D advection equation show errors in both
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem in a 2D Box Region Let us consider a model problem of parabolic
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 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 informationNumerical Methods for PDEs
Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationJMBC Computational Fluid Dynamics I Exercises by A.E.P. Veldman
JMBC Computational Fluid Dynamics I Exercises by A.E.P. Veldman The exercises will be carried out on PC s in the practicum rooms. Several (Matlab and Fortran) files are required. How these can be obtained
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationIntro to Research Computing with Python: Partial Differential Equations
Intro to Research Computing with Python: Partial Differential Equations Erik Spence SciNet HPC Consortium 28 November 2013 Erik Spence (SciNet HPC Consortium) PDEs 28 November 2013 1 / 23 Today s class
More information