hy2Foam is an open-source two-temperature computational fluid dynamics (CFD) solver that has been developed to tackle the highly complex flow physics of the hypersonic planetary atmospheric entry. Implemented within the OpenFOAM framework, the code has the capability to model physical phenomena relative to the high-speed chemically-reacting environment surrounding a spacecraft. The core of the solver initially relied on OpenFOAM solvers rhoCentralFoam and reactingFoam and it has been complemented with the many new features listed below

THERMAL NON-EQUILIBRIUM
  • addition of the electronic and electron energy modes
  • implementation of Park’s two-temperature model for a single vibro-electronic energy pool
  • modelling of energy transfers between the different energy pools
    • Vibrational-Translational (V—T): Landau-Teller
    • Heavy-particle — electron (H—e): Appleton-Bray
  • modifications to the Smoluchowski temperature jump and Maxwell velocity slip boundary conditions
  • computation of the mean free path and breakdown parameter  
CHEMICAL NON-EQUILIBRIUM
  • finite-rate volume chemistry
  • customizable chemistry databases (Park 1993, Park 1994, Dunn & Kang 1973, Scanlon 2015) including dissociation, electron impact dissociation, electron impact ionization, associative ionization, exchange and charge exchange reactions
  • chemistry-vibration coupling: Park TTv model, coupled vibration-dissociation-vibration (CVDV) model
  • polyatomics species (e.g., CO2 for the Mars atmospheric entry)  
TRANSPORT
  • species viscosity: Blottner, power law, Sutherland, constant
  • species thermal conductivity: Eucken
  • mixing rules: Wilke, Armaly & Sutton
  • species diffusion models: Lewis number, generalised Fick’s law, SCEBD model, inclusion of pressure gradient effects
  • computation of the convective & diffusive wall heat flux  
TURBULENCE
  • no changes made to OF-v1706: laminar, k-omega SST, k-Epsilon, Spalart Allmaras, etc  
MAGNETOHYDRODYNAMICS (MHD)
  • low magnetic Reynolds number model
  • MHD source terms added to the Navier-Stokes-Fourier equations: Lorentz force & Joule heating
  • various electrical conductivity models: Bush, Chapman-Cowling, Raizer, Spitzer-Harm
  • Hall effect  
MISCELLANEOUS
  • inlet boundary condition to gradually increase the free-stream flow velocity
  • dictionaries that can be re-read on-the-fly (e.g., for high-performance computing)
  • adaptive mesh refinement (DyM) based on the Mach number gradient or normalised gradients of density or pressure
  • layer functions for the strongly-coupled hybrid CFD-DSMC code hyperFoam

 

This work was initiated at the University of Strathclyde (Glasgow, UK) and was funded by the Engineering and Physical Sciences Research Council (EPSRC) from February 2014 until early 2017.
The first version of the MHD module was developed at the Ioffe Institute & Peter the Great St. Petersburg Polytechnic University (St Petersburg, Russia) in 2017-2018.
Volunteering work from mid-2017 onwards.