Two-species relaxation towards equilibrium, hotspot burn: "conduction_reaction"

Reaction model

As for `equilibrating', 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 `equilibrating'.

Input

As for `equilibrating', with in addition:

length_scale
hotspot_surface_area_per_material_volume(function of reacted volume fraction)

Notes

The length scale is the spacing between hotspots.
The geometry function is of WXC++ type Realfunction1d - see the WXC++ manual for possible types.
In order to model ~spherical burn out then burn in, with N hotspot sites per volume, a reasonable choice for the function is
If p and T are within the specified tolerances, the equilibration step is omitted. (This is for speed.)