Back to CFD

Study notes

LearnExamsInterviews

Qwestrum Engineering360 · Mechanical Engineering · CFD

Discretization Methods

Discretisation converts PDEs to algebraic equations at grid points; the central-difference approximation ∂φ/∂x ≈ (φ_{i+1} − φ_{i−1})/(2Δx) is second-order accurate. FDM, FVM, and FEM are the main schemes, per Versteeg & Malalasekera.
Exam tip: keep SI units consistent end-to-end, write the governing relation symbolically before substituting, and sanity-check magnitude and sign.

Key formulas & points

Skim these first — then read the full notes below.

  • FDM: differential form on structured grid
  • FVM: integral form, flux conservation guaranteed
  • FEM: weak form with shape functions (less common in CFD)

Topic details

Introduction

Discretisation is how CFD turns continuous equations into a computable system. Indian CFD courses cover finite-difference (FDM), finite-volume (FVM, the CFD standard), and finite-element (FEM) methods.

Scope in B.Tech and GATE syllabus

FDM approximates derivatives by Taylor-series differences on a structured grid; FVM integrates the conservation laws over control volumes, guaranteeing conservation and handling complex geometries; FEM uses weighted residuals with shape functions.

Why this topic matters in practice

Scheme choice affects accuracy (order of truncation error) and stability. Upwind schemes add numerical diffusion but are stable; central schemes are accurate but can oscillate. Deriving difference formulas and understanding accuracy/stability trade-offs are the exam skills.

Key relations & formulas

∂ϕ/∂x∣i≈(ϕi+1−ϕi−1)(2Δx)∂\phi/∂x|_i \approx \frac{(\phi_{i+1} - \phi_{i-1})}{(2\Delta x)}
(central difference, 2nd order)

Formulas (Indian textbook notation)

  • ∂2ϕ/∂x2∣i≈(ϕi+1−2ϕi+ϕi−1)Δx2∂^{2}\phi/∂x^{2}|_i \approx \frac{(\phi_{i+1} - 2\phi_{i} + \phi_{i-1})}{\Delta x^{2}}
FVM:∫CV∇⋅(ρVϕ)dV=∫CS(ρVϕ)⋅dSFVM: \int _CV ∇\cdot (\rho V\phi)dV = \int _CS (\rho V\phi)\cdot dS
(divergence theorem)
CourantnumberC=VΔtΔx<1Courant number C = \frac{V\Delta t}{\Delta x} < 1
(explicit stability)

Notation and sign conventions

Relation 1 —
∂ϕ/∂x∣i≈∂\phi/∂x|_i \approx
∂ϕ/∂x∣i≈(ϕi+1−ϕi−1)(2Δx)∂\phi/∂x|_i \approx \frac{(\phi_{i+1} - \phi_{i-1})}{(2\Delta x)}
(central difference, 2nd order)
Write this relation with symbols exactly as in Computational Fluid Dynamics — John Anderson before substituting numbers. Examiners award partial marks for a correct setup even when arithmetic slips.
Relation 2 —
∂2ϕ/∂x2∣i≈∂^{2}\phi/∂x^{2}|_i \approx

Formulas (Indian textbook notation)

  • ∂2ϕ/∂x2∣i≈(ϕi+1−2ϕi+ϕi−1)Δx2∂^{2}\phi/∂x^{2}|_i \approx \frac{(\phi_{i+1} - 2\phi_{i} + \phi_{i-1})}{\Delta x^{2}}
Write this relation with symbols exactly as in Computational Fluid Dynamics — John Anderson before substituting numbers. Examiners award partial marks for a correct setup even when arithmetic slips.
Relation 3 —
FVM:∫CV∇⋅FVM: \int _CV ∇\cdot
FVM:∫CV∇⋅(ρVϕ)dV=∫CS(ρVϕ)⋅dSFVM: \int _CV ∇\cdot (\rho V\phi)dV = \int _CS (\rho V\phi)\cdot dS
(divergence theorem)
Write this relation with symbols exactly as in Computational Fluid Dynamics — John Anderson before substituting numbers. Examiners award partial marks for a correct setup even when arithmetic slips.
Relation 4 —
CourantnumberC=VΔtΔx<1Courant number C = \frac{V\Delta t}{\Delta x} < 1
CourantnumberC=VΔtΔx<1Courant number C = \frac{V\Delta t}{\Delta x} < 1
(explicit stability)
Write this relation with symbols exactly as in Computational Fluid Dynamics — John Anderson before substituting numbers. Examiners award partial marks for a correct setup even when arithmetic slips.

Fundamentals and definitions

A Taylor expansion gives derivative approximations: forward/backward differences are first-order (O(Δx)), while the central difference (φ_{i+1} − φ_{i−1})/(2Δx) is second-order (O(Δx²)) — more accurate for the same grid.

Governing relations in practice

The finite-volume method integrates the governing equations over each control volume, so fluxes across faces are conserved exactly — mass, momentum, and energy are conserved discretely, which is why FVM dominates CFD.

Design and analysis considerations

Convective terms need care: central differencing is accurate but can produce unphysical oscillations at high cell Peclet number, while upwind differencing is stable but introduces numerical (false) diffusion that smears gradients. Higher-order schemes (QUICK, second-order upwind) balance the two.

Advanced theory and extensions

Time discretisation (explicit vs implicit) affects stability and cost: explicit schemes are simple but limited by the CFL condition, implicit schemes allow larger time steps at the cost of solving a system each step. Selecting scheme, order, and stability treatment is the essence of discretisation.

Assumptions and validity limits

State assumptions explicitly before using any relation for discretization methods — steady state, uniform properties, linear elastic material, ideal gas, incompressible flow, etc., as applicable.
Wrong assumptions invalidate the entire solution even when the formula is correct. In CFD viva and GATE descriptive questions, listing valid assumptions often earns separate marks.

Step-by-step problem approach

1. Read the question and list given data with SI units (common in CFD papers).
2. Draw a neat labelled diagram where applicable — examiners in Indian universities award diagram marks even when arithmetic slips.
3. Identify which relation from this topic applies to discretization methods.
4. Use equation 1:
∂ϕ/∂x∣i≈∂\phi/∂x|_i \approx
.
5. Use equation 2:
∂2ϕ/∂x2∣i≈∂^{2}\phi/∂x^{2}|_i \approx
.
6. Substitute values, compute, and verify units and sign (direction).
7. State conclusion in one line — e.g. safe/unsafe, stable/unstable, feasible/infeasible.

Applications & exam relevance

Discretization Methods appears in aero, HVAC ducts, and turbomachinery. In Indian mechanical curricula this topic is tested because it connects theory to computational fluid flow simulation.
GATE and semester exams often combine discretization methods with earlier units — revise prerequisites before attempting mixed problems.
Industry interview panels sometimes ask: "Where did you use discretization methods?" — answer with a lab, mini-project, or plant visit example if possible.

Common mistakes in exams

• Assuming a central-difference (2nd-order) result is stable at high Peclet number
• Confusing the order of accuracy of forward (1st) and central (2nd) differences
• Forgetting numerical diffusion introduced by upwind schemes
• Violating the CFL stability limit with an explicit time scheme

Quick revision checklist

Before attempting discretization methods problems, confirm you can:
1. FDM: differential form on structured grid
2. FVM: integral form, flux conservation guaranteed
3. FEM: weak form with shape functions (less common in CFD)
Revise the solved examples in Computational Fluid Dynamics — John Anderson and one previous-year GATE or university paper for this unit.

Worked examples

Try the problem first — open the solution when you are ready to check.

Central-difference derivative

Problem

Given φ values 10, 14, 22 at equally spaced points i−1, i, i+1 with Δx = 0.5, estimate ∂φ/∂x at i by central difference.

Solution

∂φ/∂x ≈ (φ_{i+1} − φ_{i−1})/(2Δx) = (22 − 10)/(2 × 0.5) = 12/1 = 12.

Conceptual check — Discretization Methods

Problem

In a CFD semester or GATE paper you are asked: "State the main assumption, the governing relation, and one practical consequence of discretization methods." What should a complete answer include?

Practice questions

Most-asked interview and GATE questions for this topic — expand any item for a model answer.

  1. 1
    What is Discretization Methods, and why does it appear in B.Tech / GATE syllabi?

    Model answer

    Discretisation converts PDEs to algebraic equations at grid points; the central-difference approximation ∂φ/∂x ≈ (φ_{i+1} − φ_{i−1})/(2Δx) is second-order accurate. FDM, FVM, and FEM are the main schemes, per Versteeg & Malalasekera.
  2. 2
    State the relation ∂φ/∂x|_i ≈ and name each symbol.

    Model answer

    The governing relation is ∂ϕ/∂x∣i≈∂\phi/∂x|_i \approx. Write every symbol with SI units before substituting numbers.
  3. 3
    State the relation ∂²φ/∂x²|_i ≈ and name each symbol.

    Model answer

    The governing relation is ∂2ϕ/∂x2∣i≈∂^{2}\phi/∂x^{2}|_i \approx. Write every symbol with SI units before substituting numbers.
  4. 4
    State the relation FVM: ∫_CV ∇· and name each symbol.

    Model answer

    The governing relation is FVM:∫CV∇⋅FVM: \int _CV ∇\cdot. Write every symbol with SI units before substituting numbers.
  5. 5
    State the relation Courant number C = VΔt/Δx < 1 and name each symbol.

    Model answer

    The governing relation is CourantnumberC=VΔtΔx<1Courant number C = \frac{V\Delta t}{\Delta x} < 1. Write every symbol with SI units before substituting numbers.
  6. 6
    Explain: FDM: differential form on structured grid

    Model answer

    FDM: differential form on structured grid — state the assumption range and one exam trap linked to this point.
  7. 7
    Explain: FVM: integral form, flux conservation guaranteed

    Model answer

    FVM: integral form, flux conservation guaranteed — state the assumption range and one exam trap linked to this point.
  8. 8
    Explain: FEM: weak form with shape functions (less common in CFD)

    Model answer

    FEM: weak form with shape functions (less common in CFD) — state the assumption range and one exam trap linked to this point.
  9. 9
    How would you correct this error in a viva: Assuming a central-difference (2nd-order) result is stable at high Peclet number?

    Model answer

    Identify the wrong assumption or unit mix-up, rewrite the correct relation, and recompute with a one-line sanity check.
  10. 10
    How would you correct this error in a viva: Confusing the order of accuracy of forward (1st) and central (2nd) differences?

    Model answer

    Identify the wrong assumption or unit mix-up, rewrite the correct relation, and recompute with a one-line sanity check.
  11. 11
    How would you correct this error in a viva: Forgetting numerical diffusion introduced by upwind schemes?

    Model answer

    Identify the wrong assumption or unit mix-up, rewrite the correct relation, and recompute with a one-line sanity check.
  12. 12
    How would you correct this error in a viva: Violating the CFL stability limit with an explicit time scheme?

    Model answer

    Identify the wrong assumption or unit mix-up, rewrite the correct relation, and recompute with a one-line sanity check.

Exams & GATE

  • 1
    Anderson Ch. 4 — upwind differencing for stability at high Peclet number.
  • 2
    Avoid: Assuming a central-difference (2nd-order) result is stable at high Peclet number
  • 3
    Avoid: Confusing the order of accuracy of forward (1st) and central (2nd) differences
  • 4
    Avoid: Forgetting numerical diffusion introduced by upwind schemes

📖 Standard books (India)

  • Computational Fluid Dynamics — John Anderson

    Read: Syllabus unit

    CFD fundamentals for aerospace and ME

Explore related topics

See real mechanical engineering careers

After exams and interviews, see how engineers actually built careers — milestones and decisions from people in the field.

Browse Mechanical Engineering journeysPlan your path