← All projects

Computational Fluids

Airfoil Flow, Three Methods

The flow around NACA airfoils in four phases: ANSYS Fluent (RANS) with a grid study, two physics-informed neural networks, and a source–sink potential-flow code in MATLAB.

Course
Fluid Mechanics II
Term
Fall 2025
Team
Individual project
Tools
  • Fluent
  • PyTorch
  • MATLAB

Overview

How do a RANS solver, physics-informed neural networks and a source–sink potential-flow method describe the flow around cambered and symmetric NACA airfoils, and when is each solution independent of its numerical settings?

Approach

  1. Phase 1: NACA 4421 at Re = 2×10⁶ and 3×10⁶ in ANSYS Fluent on structured C-meshes, with a grid study.
  2. Phase 2: Navier–Stokes PINN for NACA 4421 at Re = 100, α = 0°, with Fourier features and Sobol sampling.
  3. Phase 4: PINN for the potential flow around the symmetric NACA 16-006 at α = 5°.
  4. Phase 5: source–sink distribution for symmetric airfoils of five thicknesses in MATLAB.

Results

Fluent grid independence, five C-meshes (25 480 → 105 690 cells)C_L < 0.1 % and C_D ≈ 0.8 % between the two finest grids
PINN, NACA 4421 at Re = 100, surface C_p1.5 % to 25 % from the experimental data at four chord stations
PINN set-up independenceunchanged for n_int ≥ 16 000, width ≥ 128, ≥ 6 hidden layers
Source–sink method, symmetric airfoil at α = 0°C_L, circulation and C_D converge towards zero as the number of sources grows
Retained Fluent mesh78 950 cells: C_L = 0.497, C_D = 0.0216
PINN network6 × 128 tanh MLP, 16 Fourier features
PINN trainingAdam 30 000 iterations + L-BFGS
Source–sink method160 sources for the field plots

From the report

PINN velocity magnitude around NACA 4421 at Re = 100.
PINN velocity magnitude around NACA 4421 at Re = 100.
Structured C-mesh used in Fluent.
Structured C-mesh used in Fluent.
Source–sink method: streamlines with 160 sources.
Source–sink method: streamlines with 160 sources.