Compress2d

Purpose:

Two dimensional compressible Navier-Stokes solver.

Main routines:

PreProcess
Conv_Step
Save_Conv
Diff_Step_NS
Restore_Conv
cIntegrate
Set_Dir_BCs
Analyser

Descriptions:

PreProcess

Process command line arguements.
Read parameters from .rea file.
Read mesh from .rea file.
Read boundary conditions.
Create numbering for vertices and edges.
Generate matrices for elliptic solves.
Create multistep storage.

Conv_Step

Calculate Euler upwind fluxes at boundaries of elements. Also calculate divergence of Euler flux terms on volume quadrature of elements.

Save_Conv

Store results from Conv_Step in temporary storage so that the viscous terms can be calculated using the element storage.

Diff_Step_NS

Calculate viscous terms, using DGM derivatives to calculate auxiliary variables and their derivatives.

Restore_Conv

Add Euler and Viscous terms together.

cIntegrate

Use Adams-Bashforth time integrator to update state vector.

Set_Dir_BCs

Correct fields at boundary to take into account viscous boundary conditions.

Analyser

Output fields, history data and force data.