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?