A numerical simulation code for a double mach reflection problem.

dmr is a code of finite difference method for computational fluid dynamics, written in FORTRAN77. The code can simulate double mach reflection problem.

density contour (t=0.2) of double mach reflection problem

dmr was firstly (now) writtern in FORTRAN77.

The objective of the contributor is to evaluate the impact of each new feature in "Modern Fortran" on computation speed. "Modern Fortran" refers to Fortran 90/95 and later.

If you have any suggestion/knowledge/question, please post issues!

finite difference method

third order MUSCL interpolation with minmod limiter

SLAU

third order TVD Runge Kutta method

The double mach reflection problem is firstly proposed by Woodward and Collela (JCP, 1984)[1].

This problem is an important test case for the assessment of the resolution of Euler codes[2].

Only a fortran compiler is required to compile this code.

$ gfortran main.f

$./a.out

grid.xyz is a two dimensional grid file.

flowfield_?????.q are solution files.

File formats of these files are plot3d.

Any visualization softwares can be used that supports plot3d format. For example, Paraview, Fieldveiw, etc...

Automatically published documents generated by FORD are here.

This software is release under the MIT Licese, see LICENSE

Please report bugs and other issues through the issue tracker at:

https://github.com/toya42/dmr/issues

[1] Woodward and Collela, Journal of Computational Physics, 1984.

[2] Vevek et. al., Jornal of Scientific Computing, 2019.