File neutral_boundary.hxx#

Defines

NEUTRAL_BOUNDARY_H#
struct NeutralBoundary : public Component#
#include <neutral_boundary.hxx>

Per-species boundary condition for neutral particles at sheath (Y) boundaries.

Sets boundary conditions:

  • Free boundary conditions on logarithm of density, temperature and pressure

  • No-flow boundary conditions on velocity and momentum.

Adds an energy sink corresponding to a flux of heat to the walls.

Heat flux into the wall is q = gamma_heat * n * T * v_th

where v_th = sqrt(eT/m) is the thermal speed

Public Functions

NeutralBoundary(std::string name, Options &options, Solver*)#
virtual void transform(Options &state) override#

state

  • species

    • <name>

      • density Free boundary

      • temperature Free boundary

      • pressure Free boundary

      • velocity [if set] Zero boundary

      • momentum [if set] Zero boundary

      • energy_source Adds wall losses

virtual void outputVars(Options &state) override#

Add extra fields for output, or set attributes e.g docstrings.

Private Members

std::string name#

Short name of species e.g “d”.

BoutReal Tnorm#
BoutReal target_energy_refl_factor#
BoutReal sol_energy_refl_factor#
BoutReal pfr_energy_refl_factor#

Fraction of energy retained after reflection.

BoutReal target_fast_refl_fraction#
BoutReal sol_fast_refl_fraction#
BoutReal pfr_fast_refl_fraction#

Fraction of neutrals undergoing fast reflection.

Field3D target_energy_source#
Field3D wall_energy_source#

Diagnostic for power loss.

bool diagnose#
bool lower_y#

Save diagnostic variables?

Boundary condition at lower y?

bool upper_y#

Boundary condition at upper y?

bool sol#

Boundary condition at sol?

bool pfr#

Boundary condition at pfr?