DownloadFosite is a generic framework for the numerical solution of hyperbolic conservation laws in generalized orthogonal coordinates. Its main purpose is the simulation of compressible flows in accretion disks. The underlying numerical solution method belongs to the family of unsplit conservative finite volume TVD schemes. The method is 2nd order accurate in space and uses high order Runge-Kutta and multistep schemes for time evolution. In addition to the pure advection code several source terms have been implemented including viscous diffusion and gravitational acceleration.
Fosite is written with object-oriented patterns in Fortran 2003 and follows the Structure of Arrays (SoA) layout, operating on generic field datatypes. This allows for high performance on modern architectures (SIMD). It is parallelized and vectorized. The software is thereby optimized for the NEC SX-Aurora TSUBASA Vector Engine .
Features
- conservative 2nd order finite volume scheme
- solution of transport problems for compressible viscous flows
- conservation of linear and angular momentum, energy and rothalpy
- gravitational source terms including poisson solver to account for selfgravity
- various curvilinear orthogonal grids
- fast orbital advection algorithm (FARGO) for 2D polar geometries
- high-order time integration with adaptive step size control
- supported binary output formats: VTK, XDMF