Scribd
Upload
133 views

CFD Cheat Sheet: December 15, 2014

This document provides a cheat sheet on concepts in computational fluid dynamics (CFD). It defines key terms like the amplification factor and von Neumann stability analysis. It explains how to analyze the diffusion error and classify PDEs as hyperbolic, parabolic, or elliptic. It also summarizes common numerical methods like Jacobi and Gauss-Seidel iteration. The document stresses the relationship between stability, consistency, and convergence of numerical schemes. Finally, it defines common boundary conditions and describes how advection and diffusion terms affect solutions.

Uploaded by

Charles Roberson
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
133 views

CFD Cheat Sheet: December 15, 2014

This document provides a cheat sheet on concepts in computational fluid dynamics (CFD). It defines key terms like the amplification factor and von Neumann stability analysis. It explains how to analyze the diffusion error and classify PDEs as hyperbolic, parabolic, or elliptic. It also summarizes common numerical methods like Jacobi and Gauss-Seidel iteration. The document stresses the relationship between stability, consistency, and convergence of numerical schemes. Finally, it defines common boundary conditions and describes how advection and diffusion terms affect solutions.

Uploaded by

Charles Roberson
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

CFD Cheat Sheet

December 15, 2014

Part

1.1

Von Neumann Stability Analysis

n+1

Amplification factor |G| = | U n 1 (stable)


U

UIn = Uin eIi


eI(i1) = eIi eI
eI = cos() + I sin()
eI = cos() I sin()
=

at
x

Diffusion Error D =

|G|
,
|G|

where |G|=1

b2 4ac
Hyperbolic = postive (If first order PDE, always hyperbolic)
Parabolic = 0
Elliptic = negative
U n+1 U n

Ut = i t i
Jacobi = Dxk+1 = b Lxx U xk
Gauss-Seidel = (D + L)xk+1 = b U xk
Convergence Stability Consistency
If 2 characteristics are present, then all 3 are said to be true of the scheme
The solution will converge to an exact solution as t 0 and x 0
A numerical scheme is said to be stable if any component of the initial solution is amplified
without bound
As t and x 0 we should approach the exact differential equation.
Finite difference approximations can be checked quickly via:
1

U (x) =
=

Ui1 +Ui+1
2

Ui Uxi x+Ui +Uxi x


2

2Ui
2

= Ui and therefore consistent.


1D
2D
Dirichlet
u() = c
u() = f ()
Neumann
Robin

u
x ()

u() +

=d

u
x ()

u
x ()

=e

= g()

u() + u() +

u
x ()

= e()

For and equation Ut + aUx = 0, we can see that this will equal Ut = aUx , advection our
solution in a negative direction.
For and equation with a diffusive term Uxx this will diffusive our equation in a direction.

First Order Schemes: Diffusive (Smooth, Stable, Non-Accurate)


Second Order Schemes: Dispersive (Accurate, Oscillatory)

You might also like