/Compressible-Flow

Riemann solver for compressible flows.

Primary LanguageFortran

Roe Approximate Riemann Solver

Program Files

File Purpose
main.f90 Main program containing functions and subroutines for 2D Euler flow solver.
read.py Post-processing program to create Mach animation and $C_p$ plot.
airfoil.py NACA airfoil generation program to create .geo file for Gmsh.
airfoil.geo File containing geometric information to generate mesh in Gmsh.
airfoil.su2 File containing mesh information to input to main program.
check.py Mesh visualisation program with Matplotlib to verify mesh.

Program Verification

NACA Airfoil

animation.mp4

image

Diamond Airfoil

animation.mp4

image

Solver Theory

image

image

HLL Approximate Riemann Solver

movie.mp4

Domain Parameters

Parameter Value
$L_x$ $x \in$ [-1.0, 1.0]
$L_y$ $y \in$ [-1.0, 1.0]
$L_t$ $t \in$ [0.0, 10.0]
$D_x$ $d_x \in$ [-0.60, -0.40]
$D_y$ $d_y \in$ [-0.01, 0.01]
$n_x$ 200
$n_y$ 200

Initial Conditions

Parameter Value
$\tilde{\rho}_\infty$ 1.0
$\tilde{u}_\infty$ $M_\infty^x$ = 1.4
$\tilde{v}_\infty$ $M_\infty^y$ = 0.0
$\tilde{p}_\infty$ 1.0 / $\gamma$

image

HLLC Approximate Riemann Solver

movie.mp4

Domain Parameters

Parameter Value
$L_x$ $x \in$ [-1.0, 1.0]
$L_y$ $y \in$ [-1.0, 1.0]
$L_t$ $t \in$ [0.0, 10.0]
$D_x$ $d_x \in$ [-0.60, -0.40]
$D_y$ $d_y \in$ [-0.01, 0.01]
$n_x$ 200
$n_y$ 200

Initial Conditions

Parameter Value
$\tilde{\rho}_\infty$ 1.0
$\tilde{u}_\infty$ $M_\infty^x$ = 1.4
$\tilde{v}_\infty$ $M_\infty^y$ = 0.0
$\tilde{p}_\infty$ 1.0 / $\gamma$

image