File adas_reaction.hxx

Defines

ADAS_REACTION_H
struct OpenADASRateCoefficient
#include <adas_reaction.hxx>

Represent a 2D rate coefficient table (T,n) Reads data from a file, then interpolates at required values.

Public Functions

OpenADASRateCoefficient(const std::string &filename, int level)

Read the file, extracting data for the given ionisation level

Parameters:

  • filename – The file to read. Path relative to run working directory

  • level – The first index in the log coefficient array (ionisation level)

BoutReal evaluate(BoutReal T, BoutReal n)

Inputs:

Parameters:

  • n – Electron density in m^-3

  • T – Electron temperature in eV

Returns:

rate in units of m^3/s or eV m^3/s

Public Members

std::vector<std::vector<BoutReal>> log_coeff
std::vector<BoutReal> log_temperature
std::vector<BoutReal> log_density
BoutReal Tmin
BoutReal Tmax

Range of T [eV].

BoutReal nmin
BoutReal nmax

Range of density [m^-3].

struct OpenADAS : public Component
#include <adas_reaction.hxx>

Read in and perform calculations with OpenADAS data https://open.adas.ac.uk/

Uses the JSON files produced by: https://github.com/TBody/OpenADAS_to_JSON

Subclassed by ADASCarbonIonisation< level >, ADASCarbonRecombination< level >, ADASLithiumIonisation< level >, ADASLithiumRecombination< level >, ADASNeonIonisation< level >, ADASNeonRecombination< level >

Public Functions

inline OpenADAS(const Options &units, const std::string &rate_file, const std::string &radiation_file, int level, BoutReal electron_heating)

Inputs

Notes

  • The rate and radiation file names have “json_database/” prepended

Parameters:

  • units – Options tree containing normalisation constants

  • rate_file – A JSON file containing reaction rate <σv> rates (e.g. SCD, ACD)

  • radiation_file – A JSON file containing radiation loss rates (e.g. PLT, PRB)

  • level – The lower ionisation state in the transition e.g. 0 for neutral -> 1st ionisation and 1st -> neutral recombination

  • electron_heating – The heating of the electrons per reaction [eV] This is the ionisation energy, positive for recombination and negative for ionisation

void calculate_rates(Options &electron, Options &from_ion, Options &to_ion)

Perform the calculation of rates, and transfer of particles/momentum/energy

Parameters:

  • electron – The electron species e.g. state[“species”][“e”]

  • from_ion – The ion on the left of the reaction

  • to_ion – The ion on the right of the reaction

Private Members

OpenADASRateCoefficient rate_coef

Reaction rate coefficient.

OpenADASRateCoefficient radiation_coef

Energy loss (radiation) coefficient.

BoutReal electron_heating

Heating per reaction [eV].

BoutReal Tnorm
BoutReal Nnorm
BoutReal FreqNorm

Normalisations.

struct OpenADASChargeExchange : public Component
#include <adas_reaction.hxx>

Subclassed by ADASCarbonCX< level, Hisotope >, ADASLithiumCX< level, Hisotope >, ADASNeonCX< level, Hisotope >

Public Functions

inline OpenADASChargeExchange(const Options &units, const std::string &rate_file, int level)
void calculate_rates(Options &electron, Options &from_A, Options &from_B, Options &to_A, Options &to_B)

Perform charge exchange

from_A + from_B -> to_A + to_B

from_A and to_A must have the same atomic mass from_B and to_B must have the same atomic mass The charge of from_A + from_B must equal the charge of to_A + to_B

Private Members

OpenADASRateCoefficient rate_coef

Reaction rate coefficient.

BoutReal Tnorm
BoutReal Nnorm
BoutReal FreqNorm

Normalisations.