Modules

Flux splitting scheme: AUSM+ module ausmP !< Reference: Liou, M. S., “A sequel to AUSM: AUSM+,” Journal of Computational Physics, vol. 129, pp. 364–382, 1996

Setup boundary condition for the domain

Reconstruct the boundary face in case of 4th and 5th order higher order face state reconstruction method. Since the limited information is available at the boundaries, the boundary face is limiter to 3rd order accurate and is reconstructed using MUSCL Scheme even when rest of the domain is using WENO or PPM

Check for the input from the output_control.md file

A module contains generalized subroutine to copy variable in ghost cells

Interpolation of primitive variable from cell centers to faces

The geometry module calculates various geometrical quantities like face-normals, face-areas and cell-volumes to be used in computations.

Declare all the constants used by SA turbulence model

Allocate memory to laminar gradients if flow is viscous and allocate memory to tubulence gradients base upon the model being used

This module handles the MPI Communication calls for interface boundary conditions

Reference: Edwards, J.R., A low-diffusion flux-splitting scheme for Navier-Stokes calculations, Computers & Fluids, vol. 26, no. 6, pp.635-659, 1997

Setup the indicies map at interface between two blocks

Reference: Kitamura, K., Shima, E., Fujimoto, K. and Wang, Z.J., Performance of low-dissipation Euler fluxes and preconditioned LU-SGS at low speeds, Communications in Computational Physics, vol. 10 no. 1, pp.90-119, 2011

This module read input control files which include: 1. control.md 2. fvscheme.md 3. flow.md 4. res_control.md 5. state_read_write_control.md

Read output files from the restart folder

Read the restart file in the vtk format

Inviscid flux calculation through faces

Setup, run, and destroy the solver allocate/deallcoate memory, initialize, iterate

Allocate memory to the state variables and initialize them The state of the system is defined using the density, velocity and pressure (primitive variables qp), and trubulent and transition variables at the cell-center points.

This module march the solution is time.

Reference: Van Leer, B., Flux-vector splitting for the Euler equation. In Upwind and High-Resolution Schemes, Springer, Berlin, Heidelberg, pp. 80-89, 1997

Setup, destroy, calculate molecular and turbulence viscosity

Detect all the grid points on the wall boundary condition and store them in a single file

Reference: 1 Shu, C.-W., “High-order Finite Difference and Finite Volume WENO Schemes and Discontinuous Galerkin Methods for CFD,” International Journal of Computational Fluid Dynamics, vol. 17, 2003, pp. 107–118. Reference: 2 Huang, W. F., Ren, Y. X., and Jiang, X., “A simple algorithm to improve the performance of the WENO scheme on non-uniform grids,” Acta Mechanica Sinica/Lixue Xuebao, 2017, pp. 1–11.

Open/close and call other modules for writing solution based on the input: type of file, either vtk or tecplot modules are called

Writing solution in the output file in tecplot format with node data instead of cell-center data.