Aerodynamics, Aerothermodynamics, and Acoustics Competency Logo

Hypersonic Airbreathing Propulsion Branch
Hypersonic Airbreathing Propulsion Branch

VULCAN-CFD: Roman God of the ForgeVULCAN-CFD, named after the Roman god of fire, is a Computational Fluid Dynamics software package (available for serial and parallel computational platforms) for turbulent reacting and non-reacting flows at conditions ranging from subsonic to hypersonic speeds. The computational cost of propulsion flow analysis is reduced through the use of special turbulent wall treatments, multi-grid methods for elliptic and space marching schemes, and conditioning of the governing equations to reduce numerical stiffness. VULCAN-CFD can simulate two-dimensional, axi-symmetric, or three-dimensional problems on structured multi-block grid systems. A variety of PDE-based turbulence models are available (including explicit algebraic Reynolds stress models) for Reynolds-averaged simulations, one-equation and algebraic sub-grid closures (including dynamic variants) are available for Large-eddy simulations, and hybrid Reynolds-Averaged / Large Eddy simulation capabilities are also present. A variety of robust upwind-biased algorithms are for Reynolds-averaged simulations, and a low-dissipation numerical framework is available for scale-resolving simulations. VULCAN-CFD offers significant geometric flexibility; boundary conditions can be imposed on any boundary or boundary subset, and the code has (C0)continuous as well as non-C(0) continuous block-to-block interface capabilities. VULCAN-CFD also utilizes a flexible thermodynamic and kinetic model specification framework based on user-specified databases (many of them are standard formats utilized by the propulsion community). The working fluid can be simulated as a calorically perfect single component gas, a mixture of thermally perfect gases (with or without chemical reactions), or a mixture of gases in thermodynamic non-equilibrium (with or without chemical reactions).

For a detailed description of the current capabilities of VULCAN-CFD see the VULCAN-CFD USER MANUAL


Multiple Unix and Linux platform support for either
     1) Block level parallelization using MPI or MPICH for multi-processor machines
     2) Serial execution for single CPU machines
2-D, axi-symmetric, or 3-D structured multi-block grid topologies
     - Automated 3-point, 5-point, and 6-point topology singularity detection
Block-to-block interface options:
     1) Arbitrary block-to-block C(0) continuous connectivity
     2) Arbitrary block-to-block non-C(0) continuous connectivity
Steady-state algorithms for spatially elliptic/hyperbolic and parabolic/hyperbolic equations
     1) Runge-Kutta with implicit residual smoothing
     2) Diagonalized approximate factorization
     3) Incomplete LU
Unsteady algorithms
     1) Runge-Kutta
     2) Diagonalized approximate with dual time-stepping
     3) Incomplete LU with dual time-stepping
Convergence acceleration options:
     1) Multi-grid
     2) Mesh sequencing
     3) Low Mach no. preconditioning
Gas models
     1) Single component calorically and thermally perfect gases
     2) Arbitrary multi-component mixtures of thermally perfect gases
     3) Arbitrary multi-component mixtures of gases in thermal non-equilibrium
Chemistry models
     1) Frozen flow
     2) Arbitrary, non-equilibrium, finite-rate chemical kinetics
Inviscid flux reconstruction algorithms
     1) MUSCL kappa scheme of van Leer
     2) PPM scheme of Colella and Woodward
     3) WENO scheme of Carpenter and Fisher
     4) 2nd, 4th, and 6th order symmetric (central) schemes
Inviscid flux schemes
     1) Roe flux difference split scheme with entropy fixes
     2) Edwards low dissipation flux split scheme
     3) Toro HLLC scheme
     4) Carpenter & Fisher α - split flux scheme
Robust low-dissipation inviscid flux schemes
     - Hybridization of dissipative and non-dissipative schemes via a discontinuity sensor
         a) Blend of dissipative (MUSCL, PPM, or WENO) with symmetric reconstruction
         b) Blend of dissipative (MUSCL, PPM, or WENO) with the α - split flux scheme
Viscous flux treatments for laminar or turbulent flow
     1) Full Navier-Stokes
     2) Thin-layer Navier-Stokes
     3) Parabolized Navier-Stokes
Turbulence model options for Reynolds-Averaged Simulations
     1) Spalart-Allmaras (conservative form)
     2) Menter k-omega (BSL)
     3) Menter k-omega (SST)
     4) Wilcox k-omega (1998) (low and high Reynolds number forms)
     5) Wilcox k-omega (2006) (low and high Reynolds number forms)
     6) Gatski-Speziale EASM k-omega
Turbulence model options for Large Eddy Simulations
     1) Smagorinsky with van Driest wall damping
     2) Vreman
     3) Dynamic Smagorinsky
     4) Dynamic Vreman
Surface turbulence boundary condition treatment options
     1) Wilcox compressible wall matching (k-omega based models)
     2) Solve to wall (all models)
Hybrid Reynolds-Averaged / Large Eddy Simulation options
     1) Detached Eddy Simulation (two-equation based model of Strelets)
     2) Hybrid RAS/LES model of Edwards and Baurle
Turbulence/Chemistry interaction closure models
     1) Eddy break-up model of Magnussen and Hjertager
     2) Temperature fluctuations: average reaction rate coefficient
         a) Assumed gaussian probability density function (PDF)
         b) Assumed beta probability density function (PDF)
     3) Species fluctuations: average product of concentration
         - Assumed multi-variate beta probability density function (PDF)


Hybrid structured/unstructured simulation capability
Compressible Flamelet Generated Manifold capability for efficient chemical kinetics treatment
Filtered Density Function capability for improved turbulence-chemistry interaction closure

Home Page
Graphical User Interface
Input Manual
Software Request


NASA Official Responsible For Content: Robert A. Baurle

This site's most recent update took place: 06/02/2014

Site Curator: Robert A. Baurle (

Privacy Policy and Important Notices Copyright Notice