Scribd
Upload
3 views11 pages

Numerical Analysis Techniques in CFD

The document outlines a series of tasks related to basic numerical analysis in computational fluid dynamics (CFD), including derivations, classifications, and simulations of various equations and methods. It emphasizes the importance of justifications and proper documentation of solutions, which should be submitted in LaTeX or a readable format. Key topics include ideal gas equations, PDE characteristics, stability conditions, numerical integration methods, and implicit time-stepping for advection and diffusion equations.

Uploaded by

ledecarte
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
3 views11 pages

Numerical Analysis Techniques in CFD

The document outlines a series of tasks related to basic numerical analysis in computational fluid dynamics (CFD), including derivations, classifications, and simulations of various equations and methods. It emphasizes the importance of justifications and proper documentation of solutions, which should be submitted in LaTeX or a readable format. Key topics include ideal gas equations, PDE characteristics, stability conditions, numerical integration methods, and implicit time-stepping for advection and diffusion equations.

Uploaded by

ledecarte
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

BASIC NUMERICAL ANALYSIS IN CFD.

Note-The solutions for the numerical part must be


submitted in Latex(preferably) or can also be submitted
as a pdf or in clearly readable written form on a blank A4
size paper and proper justifications and reasonings
should be provided as and where required.

[Link] the equations of state for e and T for an ideal gas (see Eq. 2-5 and Eq. 2-6)

[Link] A-Classify the following two PDE's

Part B-Find the characteristics of each of the following PDE's


Part C-Find the solution of the wave equation

with initial data 𝑢(𝑥, 0) = 1 and 𝑢𝑦 (𝑥, 0) = 0

Part [Link] the second derivative (Eq. 4-10) as a difference of first derivatives.

Part E
[Link] A-Derivation of the stability condition for the heat equation

Part B-Excel simulation of 1D heat equation

B.1 Excel simulation of 1D wave equation using central, forward,backward discretization


,Lax-Wendroff scheme ,MacCormack scheme and using 4 step Runge-Kutta scheme (central)

B.2 Derivation of CFL condition starting from the amplification function Eq. 6-22
B.3 Calculate the error of the upwind scheme for different angular frequencies (, 4 , 6 , 8 , 12 ,
24 ) and CFL numbers (0.75, 0.5, 0.25) (see 6.4.2 Example)
[Link] Integration of y‘= 2x using an explicit, implicit and trapezoidal method (Boundary
condition: y(0) = 0)

4.1-Convert the MacCormack algorithm Eq. 6-65 into the Lax-Wendroff algorithm Eq. 6-60
5.3

5.4 Derivation of the Newton algorithm and the derivation of the delta form from the implicit
discretization of Burgers equation (Eq. 6-84)

6. Part A-Calculate the 2nd row of the Jacobian matrix (𝐴) given below

Also Calculate: Eigenvalues, Eigenvectors, diagonalize 𝐴 (Attention: 𝑅 must be inverse of 𝐿, or


eigenvector matrices have to be normalized with the determinant) [Excel-Commands: minv,
mmult with Strg+Shift+Enter]

Part B-Derivation of the momentum equation for 1D- Euler flow (Eq. 7-24)

Part C-Derivation of the Lax-Wendroff scheme (Eq. 8-4)


Part D-Reproduce the ADI-factorization of the 2D Beam& Warming scheme (without numerical
dissipation) and find the additional 4th order term on the implicit side

[Link] A-Derivation of the relations for 𝑝 and 𝑢 at the boundaries (Eq. 10-5, Eq. 10-6).

Part-B Excel programming of the shock tube problem (eventually, comparison with the
MCC-Code): conditions left: 𝑝 = 1 𝑏𝑎𝑟, 𝜌 = 1 𝑘𝑔 𝑚3 ; conditions right: 𝑝 = 0.1 𝑏𝑎𝑟, 𝜌 = 0.125 𝑘𝑔
𝑚3

Part-C Show the identity of Eq. 11-42 and Eq. 11-43

Part-D Show the identity of Eq. 11-68, Eq. 11-69 and Eq. 11-70.
[Link]-A Determine the 𝑅 matrix for one dimensional flow from the general matrix (Eq. 11-87)

Part-B Calculation of the state vector at the cell interfaces (Eq. 12-6, Eq. 12-7) with the help of
the MUSCL approach (Eq. 12-1, Eq. 12-2).
Part-C Derive the NVF of the second order (Eq. 12-10) and third order (Eq. 12-8) MUSCL
schemes. Show that they are at least second order accurate!

Part -D Derive the differential 𝜕𝑇/ 𝜕𝑥 at (0,0) for a quadratic auxiliary cell with the vertex
coordinates (1,0), (0,1), (-1,0) and (0,-1) using the finite volume discretization based on the
Gauss theorem

Part-E Prove the identities Eq. 1 and Eq. 2 of the Favre averaging below-
Note-Please make a proper and detailed code and document each and every part
of the code with proper justification and to the point and exact way [Link]
graphs and outcomes should be properly documented and well [Link]
values assumed or taken shall be marked and specified.

[Link] Time stepping for the Advection Equation in 1D

on a uniform grid with m = 1000 points, using a method-of-lines (MOL) approach.


Spatially discretize this equation using the third-order upwind biased scheme , and integrate it
forward in time to a time T = 1 using a time step size corresponding to an advective CFL
number of ν = a∆t/∆x = 1, 10, 100 and the following implicit temporal integrators:
1. Backward Euler method.
2. Implicit trapezodial method (Crank-Nicolson).
Compare the three numerical solutions w(t = 1) to the exact solution u(x, 1) and comment on
your observations. Explain how you solved the linear system arising in the implicit temporal
discretization.
9.1 Implicit Timestepping for the Diffusion Equation

One Dimension

9.2 Two dimensions

You might also like