cfd-institute / euler2d-kfvs

2D solver for Euler equations in quadrilateral grid, using kinetic flux vector splitting scheme, written in OOP F2003

euler2D-kfvs-Fortran2003

Euler2D-kfvs-Fortran2003 is a compressible Euler equations solver using kinetic flux vector splitting (KFVS) scheme. It is the successor of my Fortran 95 code euler-2D-kfvs

Some features:

• Cell centered finite volume method.

• Unstructured quadrilateral meshes.

• Kinetic flux vector splitting (KFVS) scheme.

• Written in Modern Fortran (Fortran 2003). Object-oriented programming has been implemented, mainly in the mesh reader. The solver is kept in Fortran 95 now, to retain its execution speed, which is the major advantage of the Fortran programming language. For further development, maybe I will upgrade the solver to Fortran 2003, in order to integrate new numerical schemes.

• My code makes use of GMSH as its mesh generator (quadridlateral grid). You need to read the mesh reader of the code (msh_read_gmsh.f90) to understand how to generate a mesh which is suitable for my code.

Euler2D-kfvs-Fortran2003 is inspired by the structured finite volume solver of Prof. Luc Mieussens

More details

Donations:

If you think that any information you obtained here is worth of some money and are willing to pay for it, feel free to send any amount through paypal.

euler2d-kfvs's Issues

Build errors

Truong,

I am trying to compile the code using the PGI 18.7 compilers and I'm running into an error. I use pgfortran as the compiler and get the following output:

make
pgfortran -c mod_vector_algebra.f90
pgfortran -c mod_read_gmsh.f90
pgfortran -c mod_point_util.f90
pgfortran -c mod_cell_2D.f90
pgfortran -c mod_detect_nearest_neighbor.f90
PGF90-S-0155-detect_neighbor is use associated and cannot be redeclared (mod_detect_nearest_neighbor.f90: 7)
0 inform, 0 warnings, 1 severes, 0 fatal for detect_neighbor
make: *** [mod_detect_nearest_neighbor.o] Error 2

It looks like mod_cell_2D.f90 and mod_detect_nearest_neighbor.f90 both have subroutines "detect_neighbor":

grep detect_neighbor *.f90
mod_cell_2D.f90: procedure :: find_neighbor => detect_neighbor
mod_cell_2D.f90: subroutine detect_neighbor(this)
mod_detect_nearest_neighbor.f90: subroutine detect_neighbor

Is this a bug in the code?

Thank you!

Jeff Layton