Stress and temperature relaxation towards equilibrium, hotspot burn: "surface_reactions"

Reaction model

As for `reactive', except that an additional contribution is added to the reaction rate to reflect the different temperature and pressure near a hotspot.

Hotspots are modelled by giving the initial material a non-zero volume fraction of products. The initial state in the products should be chosen carefully to be at the same pressure and temperature as the reactant.

Equilibration method

As for `reactive'.

Input

As for `reactive', with in addition:

For each reaction:
   flame_velocity (WXC++ Realfunction2d function of density, temperature)
For each material i in mixture:
   For each material j in mixture:
      common_surface_area_per_volume (WXC++ Realfunction2d function of volume fraction i, volume fraction j)