File thermal_force.hxx

Defines

THERMAL_FORCE_H

struct ThermalForce : public Component

#include <thermal_force.hxx>

Simple calculation of the thermal force

Important: This implements a quite crude approximation, which is intended for initial development and testing. The expressions used are only valid for trace heavy ions and light main ion species, and would not be valid for Helium impurities in a D-T plasma, for example. For this reason only collisions where one ion has an atomic mass < 4, and the other an atomic mass > 10 are considered. Warning messages will be logged for species combinations which are not calculated.

Options used:

• <name>

  • electron_ion : bool Include electron-ion collisions?

  • ion_ion : bool Include ion-ion elastic collisions?

Public Functions

inline ThermalForce(std::string name, Options &alloptions, Solver*)

virtual void transform(Options &state) override

Inputs

• species

  • e [ if electron_ion true ]

    • charge

    • density

    • temperature

  • <species>

    • charge [ Checks, skips species if not set ]

    • AA

    • temperature [ If AA < 4 i.e. “light” species ]

Outputs

• species

  • e

    • momentum_source [ if electron_ion true ]

  • <species> [ if AA < 4 (“light”) or AA > 10 (“heavy”) ]

    • momentum_source

Private Members

bool electron_ion

Include electron-ion collisions?

bool ion_ion

Include ion-ion elastic collisions?

bool first_time = {true}

True the first time transform() is called.