jax_sbgeom.interfaces.plasma module

flux_surface_reaction_rates_simple(s_values: Array, n0: float, nalpha: float, Ti0: float, Tialpha: float)[source]

Function to create a simple neutron reaction rate profile as a function of the flux surface label s.

All profiles are defined as:

n(s) = n0 * (1 - s)**n_alpha T(s) = T0 * (1 - s)**Ti_alpha

Tritium and deuterium densities are assumed equal (1/2 of n(s)). Edge values are set to zero.

Parameters:

  • s_values (Array) – Normalized flux surface label [0, 1]

  • n0 (float) – Central ion density m^-3

  • nalpha (float) – Ion density profile exponent

  • Ti0 (float) – Central ion temperature in keV

  • Tialpha (float) – Ion temperature profile exponent

  • n_samples (int)

Returns:

reaction_rate – Reaction rate profile in m^-3 s^-1

Return type:

jnp.ndarray