/horses3d

HORSES3D: A high-order discontinuous Galerkin solver for flow simulations and multi-physics applications

Primary LanguageFortranMIT LicenseMIT

HORSES3D High-Order (DG) Spectral Element Solver

MIT License

Copyright (c) 2021 NUMATH https://numath.dmae.upm.es

Synopsis

HORSES3D is a multiphysics environment where the compressible Navier-Stokes equations, the incompressible Navier–Stokes equations, the Cahn–Hilliard equation and entropy–stable variants are solved. Arbitrary high–order, p–anisotropic discretisations are used, including static and dynamic p–adaptation methods (feature-based and truncation error-based). Explicit and implicit time-steppers for steady and time-marching solutions are available, including efficient multigrid and preconditioners. Numerical and analytical Jacobian computations with a coloring algorithm have been implemented. Multiphase flows are solved using a diffuse interface model: Navier–Stokes/Cahn–Hilliard. Turbulent models implemented include RANS: Spalart-Allmaras and LES: Smagorinsky, Wale, Vreman; including wall models. Immersed boundary methods can be used, to avoid creating body fitted meshes. Acoustic propagation can be computed using Ffowcs-Williams and Hawkings models.

HORSES3D supports curvilinear, hexahedral, conforming meshes in GMSH, HDF5 and SpecMesh/HOHQMesh format. A hybrid CPU-based parallelisation strategy (shared and distributed memory) with OpenMP and MPI is followed.

External libraries

The following external routines/libraries can be used with HORSES3D, but are not necessary: METIS, MPI, HDF5, MKL, PETSc.

Compilers and third-party software

HORSES3D is an object-oriented Fortran 2008 solver, that can be compiled using gcc and the Intel compiler, in Unix-based operating systems.

  • We recommend using recent versions of such compilers (2019 or newer).

  • Make is necessary (e.g., Gnu's version, which is included in most linux distributions).

  • Supported meshes are in GMSH, HDF5 (HOPR) and SpecMesh/HOHQMesh format.

  • Post processing can be performed in tecplot or paraview.

Compiling & Running

  1. Go to the Solver folder and configure the project

    cd Solver
./configure

  2. Build the solvers using make (see manual in the /doc folder)

    make clean
make all [options]

    with the desired options (defaults are bold):

    • PLATFORM=MACOSX/LINUX

    • MODE=DEBUG/RELEASE

    • COMPILER=ifort/gfortran

    • COMM=PARALLEL/SEQUENTIAL

    • ENABLE_THREADS=NO/YES

    • WITH_PETSC=YES/NO

    • WITH_METIS=YES/NO

    • WITH_HDF5=YES/NO

    • WITH_MKL=YES/NO

    For example:

    make all COMPILER=ifort COMM=PARALLEL

  3. Run the solver for the parameter file file.control (see manual in the /doc folder and examples in /test)

    ./horses3d.ns file.control

  4. Test cases for various physics are provided in the folder /test

Additional libraries

  • PETSc:

    export PETSC_DIR=path_to_petsc

  • METIS:

    export METIS_DIR=path_to_metis

  • HDF5:

    export HDF5_DIR=path_to_hdf5