Initial and Boundary Conditions in Computational Fluid Dynamics

1. Introduction

Initial and boundary conditions are fundamental in Computational Fluid Dynamics (CFD) as they are crucial for defining a unique solution to Partial Differential Equations (PDEs) that govern fluid flow and heat transfer processes. The manner in which these conditions are implemented numerically can significantly affect the stability and convergence of the computational solution. In numerical methods, it may sometimes be necessary to specify more conditions on a boundary than required by the analytical formulation of the problem.

2. Classification of PDEs and Required Conditions

The classification of PDEs—hyperbolic, parabolic, or elliptic—determines the initial and/or boundary conditions necessary for a well-posed problem.

2.1 Elliptic PDEs: Equilibrium Problems

• Nature: Associated with equilibrium (steady-state) problems.
• Domain of Dependence: Entire boundary influences every point inside the domain.
• Examples: Laplace's equation, steady-state heat conduction.
• Boundary Conditions:
  • Dirichlet: Prescribes the function value (e.g., temperature T).
  • Neumann: Prescribes the derivative (e.g., heat flux ∂T/∂n).
  • Robin: Combines both function and its derivative.
• Ill-posed Example: Solving Laplace’s equation on an open domain with boundary conditions on one line only.

2.2 Parabolic PDEs: Marching Problems

• Nature: Govern diffusion processes (e.g., unsteady heat conduction).
• Domain of Dependence: Entire domain up to a time t₁.
• Examples: Transient heat conduction in solids.
• Conditions: Requires both initial and boundary conditions.

2.3 Hyperbolic PDEs: Wave Problems

• Nature: Describe wave propagation and transport phenomena.
• Domain of Dependence: Limited to regions bounded by characteristic curves.
• Examples: Wave equation, Euler equations.
• Conditions: Need initial values and first derivatives.
• Constraint: Initial curve must not align with a characteristic curve.
• Boundary Conditions: Determined by the direction of characteristics.

3. Numerical Implementation

3.1 General Approaches

• Taylor Series: Approximates derivatives near boundaries.
• Polynomial Fitting: Used near walls for gradients.
• Finite-Volume Method: Applies conservation laws over control volumes.

3.2 Boundary Condition Scenarios

Solid Walls (No-Slip)

• Set u = v = 0 at walls.
• Implemented using:
  • Direct assignment
  • Ghost cells or reflection
  • Modified conservation equations

Symmetry Lines

• Conditions: ∂u/∂y = 0, v = 0 at the centerline.

Inflow/Outflow Boundaries

• Supersonic Inflow: All variables specified.
• Subsonic Inflow: Some variables specified (e.g., density).
• Supersonic Outflow: Use extrapolation from interior.
• Subsonic Outflow: One variable specified (e.g., pressure).

Shock Waves

• Treated as discontinuities with special handling.
• Rankine–Hugoniot relations used to relate variables across the shock.

Irregular Meshes and Curved Boundaries

• Use:
  • Fine mesh with zig-zag approximation
  • Interpolation techniques
  • Body-fitted coordinates
  • Finite-volume on irregular cells

Thermal Boundary Conditions

• Dirichlet: Fixed temperature.
• Neumann: Fixed heat flux.
• Robin: Mixed condition (e.g., convective boundary).

Wakes

• Velocity potential has a jump; circulation enforced.
• Normal derivatives may be extrapolated.

Far-Field Conditions

• Applied at distant boundaries.
• Kutta condition may be used for circulation (e.g., airfoils).