This test case represents a low-enthalpy flow for which the assumption of a perfect gas has been shown to be valid. Experimental data for surface pressure and heat transfer rates exists. Numerical work comparing approximate Riemann solvers with flux-split methods has matched the experimental results. Numerical results are often characterized by carbuncle phenomena, which can be suppressed through tuning of an entropy correction factor for some schemes.
The computational mesh contains two 2D blocks. Both of the block sizes are 79 ´ 71. The mesh geometry is shown in Figure 2.
Both mesh files and project input files can be accessed below. Remember to place the grid files in a subfolder with the set up file /Hyper.
Change the directory to the subfolder with the selected grid and spatial scheme. Start the simulation by
mpirun –np 2 mpiaeroflo.exe < hyperbb.afl
The calculation starts from a small time step size Δt = 1.0×10-6 for 10,000 steps. Then calculate the simulation with a larger time step size Δt = 5.0×10-6 for 50,000 steps. Finally the simulation is calculated with a time step size Δt = 1.0×10-5 for 1,000,000 steps.
Figures 3, 4 and 5 are the nondimensionalized pressure contours for the MUSCL, WENO 33 and WENO 34 calculations, respectively.
Figure 6 shows the surface pressure distribution for different spatial schemes. The results are also compared with experimental results. Note that both of the WENO 33 and WENO 34 schemes blow out during the calculation.
Figure 7 compares the convergence performance of the three spatial schemes.
© 2025 All Rights Reserved.
