toya42/dmr

A numerical simulation code for a double mach reflection problem.

dmr

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

Motivation of this repository

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!

Numerical method

spatial discretization

finite difference method

interpolation and limiter

third order MUSCL interpolation with minmod limiter

Riemann solver

SLAU

time integration

third order TVD Runge Kutta method

Double Mach Reflection problem

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].

Dependency

Only a fortran compiler is required to compile this code.

Compile & Execute

Compile

$ gfortran main.f

Execute

$./a.out

Visualize

Output file format

grid.xyz is a two dimensional grid file.

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

File formats of these files are plot3d.

Visualize

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

Document

Automatically published documents generated by FORD are here.

License

This software is release under the MIT Licese, see LICENSE

Contact

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

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

References

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

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