We developed a didactic wind tunnel based on a simplified Computational Fluid Dynamics model for incompressible flows, designed for educational and exploratory purposes.
The solver enforces mass conservation and momentum balance:
- ∇·u = 0
- ∂u/∂t + (u·∇)u = −∇p/ρ + ∇·[(ν+νₜ)∇u
The following modeling and numerical assumptions are adopted:
- Uniform staggered grid (velocities stored on cell faces; pressures at cell centers);
- Explicit time integration;
- Pressure projection method with iterative Poisson solvers: over-relaxed Gauss–Seidel (CPU mode) and Jacobi or Black–Red Gauss–Seidel (GPU mode);
- Semi-Lagrangian advection;
- Laminar viscosity with simplified local Smagorinsky LES (turbulent viscosity computed locally from the strain-rate tensor).
Below an example