Skip Navigation | ANU Home | Search ANU | Directories
The Australian National University
Research School of Astronomy and Astrophysics
Fyris Alpha Simulation Code

Fyris Alpha

by Dr Ralph Sutherland, RSAA, ANU

Fyris Alpha is a high resolution, shock capturing, multi-phase, up-wind Godunov method hydrodynamics code, which includes a variable equation of state, and optional microphysics such as cooling, gravity and mutlitple tracer variables. Fyris Alpha is a successor to earlier generations of hydrodynamics codes which arose from the original work by Collella and Woodward (1984) and earlier pioneering work by van Leer (1979) on high resolution shock capturing hydrodynamics and the PPM algorithm. The Riemann solver used is a new solver, derived from the exact solver of Gottlieb and Groth (1988), which is significantly faster than previous solvers.

The new code (and it's framework ZLS) has been designed and developed at RSAA for use primarily in astrophysical applications, such as galactic and interstellar bubbles, hypersonic shocks, and a range of jet phenomena. Fyris Alpha boasts both higher performance and more detailed microphysics than it's predecessors, with the aim of producing output that is closer to the observational domain, such as emission line fluxes, and eventually, detailed spectral synthesis.

Fyris Alpha is approximately 85,000 lines of C code (not including third party libraries), and was begun during a study visit to UAO, supported by the Centre for Complex Dynamics, at the Ångström Labotoriet, Uppsala University, Sweden and the code is named after the river that flows through that town.

The code is conceptually made up of an overarching framework, called ZLS, which provides all the standard services needed for any simulation, plus provision for modules of specific algorithm code. Fyris Alpha is the first of these modules, and it encapsulates the split sweep semi-lagrangian remap PPM method used by ppmlr (in turn developed from VH1, Blondin et al. 1998) but with an improved Riemann solver. It also has a number of optimisations that have improved the speed so that additional calculations neeed for multi-phase simulations become practical. Future modules, including magnetohydrodynamics (Fyris Beta) are in development.

The The code paper documenting the hydrodynamic module Fyris Alpha has been accepted, and posted to astro-ph. A release version of the code is now available.

Much of the code was developed with funding from the ARC Discovery Project DP 0664434.


Fyris Alpha: 3D Kelvin Helmholtz Instability.