File evolve_energy.hxx

Defines

struct EvolveEnergy : public Component 

#include <evolve_energy.hxx>

Evolves species internal energy in time

Mesh inputs

P<name>_src A source of pressure, in Pascals per second This can be over-ridden by the source option setting.

Public Functions

EvolveEnergy(std::string name, Options &options, Solver *solver)

Inputs

<name>
- bndry_flux Allow flows through radial boundaries? Default is true
- density_floor Minimum density floor. Default 1e-5 normalised units.
- diagnose Output additional diagnostic fields?
- evolve_log Evolve logarithm of pressure? Default is false
- hyper_z Hyper-diffusion in Z
- kappa_coefficient Heat conduction constant. Default is 3.16 for electrons, 3.9 otherwise
- kappa_limit_alpha Flux limiter, off by default.
- poloidal_flows Include poloidal ExB flows? Default is true
- precon Enable preconditioner? Note: solver may not use it even if enabled.
- thermal_conduction Include parallel heat conduction? Default is true
E<name> e.g. “Ee”, “Ed+”
- source Source of energy [W / s]. NOTE: This overrides mesh input P<name>_src
- source_only_in_core Zero the source outside the closed field-line region?
- neumann_boundary_average_z Apply Neumann boundaries with Z average?

virtual void transform(Options &state) override

Inputs

Sets

virtual void finally(const Options &state) override

Optional inputs

species
- <name>
  - velocity. Must have sound_speed or temperature
  - energy_source
  - collision_rate (needed if thermal_conduction on)
fields
- phi Electrostatic potential -> ExB drift

virtual void outputVars(Options &state) override: Add extra fields for output, or set attributes e.g docstrings.

virtual void precon(const Options &state, BoutReal gamma) override: Preconditioner

Private Members

bool bndry_flux: Heat capacity at constant volume (3/2 for ideal monatomic gas)

BoutReal kappa_coefficient: Leading numerical coefficient in parallel heat flux calculation