HagelaarEnergyAdvectionBC

Kinetic advective electron energy boundary condition (Based on Hagelaar et al. (2000))

Overview

HagelaarEnergyAdvectionBC is an electric field driven outflow boundary condition with the addition of ion induced secondary electron energy. HagelaarEnergyAdvectionBC assumes the electrostatic approximation for the electric field.

The outflow is defined as

$var element = document.getElementById("moose-equation-76852ce6-8aa7-4a0f-bfc8-982118fb86a8");katex.render("a = \\begin{cases} 1, & \\mu_{e} \\ \\nabla (V) \\cdot \\textbf{n} > 0\\\\ 0, & \\mu_{e} \\ \\nabla (V) \\cdot \\textbf{n} \\leq 0\\\\ \\end{cases} \\\\[10pt] \\Gamma_{\\epsilon} \\cdot \\textbf{n} = \\frac{1-r_{e}}{1+r_{e}} \\left[ -(2 a_{e}-1) \\frac{5}{3} \\mu_{e} \\text{-} \\nabla (V) \\cdot \\textbf{n} \\ n_{e} \\ \\epsilon \\right] - \\frac{5}{3} \\epsilon_{\\gamma} (1 - a_{e}) \\gamma \\left[ \\mu_{+} \\ \\left( \\text{-} \\nabla (V) \\right) \\cdot \\textbf{n} \\ n_{+} \\right]", element, {displayMode:true,throwOnError:false});$

Where $var element = document.getElementById("moose-equation-9f493a3d-2758-4ca5-bce1-ece16d5e63e3");katex.render("\\Gamma_{\\epsilon} \\cdot \\textbf{n}", element, {displayMode:false,throwOnError:false});$ is the outflow of the electron energy density normal to the boundary, $var element = document.getElementById("moose-equation-699c916c-6016-4733-ad86-c9500e24eb3c");katex.render("\\textbf{n}", element, {displayMode:false,throwOnError:false});$ is the normal vector of the boundary, $var element = document.getElementById("moose-equation-dd89bab0-f0ce-4e18-a4b3-5eb7e1f3a569");katex.render("n_{e}", element, {displayMode:false,throwOnError:false});$ is the electron density, $var element = document.getElementById("moose-equation-ab46ca52-421a-4d79-80b7-efd5c6baf8e6");katex.render("n_{+}", element, {displayMode:false,throwOnError:false});$ is the mean energy density, $var element = document.getElementById("moose-equation-4bc38486-2fb4-47e7-a8ab-593c37749c25");katex.render("\\mu_{j}", element, {displayMode:false,throwOnError:false});$ is the mobility coefficient, $var element = document.getElementById("moose-equation-4bd85c11-df0c-43b3-ae40-76572f7f5e41");katex.render("V", element, {displayMode:false,throwOnError:false});$ is the electric potential, $var element = document.getElementById("moose-equation-424cddf3-b1dd-4b38-87c9-9577cdaeb510");katex.render("\\epsilon_{\\gamma}", element, {displayMode:false,throwOnError:false});$ is the energy of the secondary electron, and $var element = document.getElementById("moose-equation-4d3684e6-31c7-4a97-81d7-c56dfedac337");katex.render("\\gamma", element, {displayMode:false,throwOnError:false});$ is the ion induced secondary electron coefficient. $var element = document.getElementById("moose-equation-bc45d01b-c517-4e22-90d6-c444d41f129d");katex.render("a", element, {displayMode:false,throwOnError:false});$ is defined such that the outflow is only defined when the drift velocity is directed towards the wall and zero otherwise. $var element = document.getElementById("moose-equation-0866d44f-4288-4c81-8cee-604953b9dfcc");katex.render("r", element, {displayMode:false,throwOnError:false});$ is defined as the fraction of particles reflected by the surface. When converting the density to log form and applying a scaling factor of the mesh, the strong form for HagelaarEnergyAdvectionBC is defined as

$var element = document.getElementById("moose-equation-a71c71f0-96d8-4a78-b5f6-cfa658b0dc5e");katex.render("a = \\begin{cases} 1, & \\mu_{e} \\ \\nabla (V) \\cdot \\textbf{n} > 0\\\\ 0, & \\mu_{e} \\ \\nabla (V) \\cdot \\textbf{n} \\leq 0\\\\ \\end{cases} \\\\[10pt] \\Gamma_{\\epsilon} \\cdot \\textbf{n} = \\frac{1-r_{e}}{1+r_{e}} \\left[ -(2 a_{e}-1) \\frac{5}{3} \\mu_{e} \\left( \\text{-} \\nabla (V / l_{c}) \\right) \\cdot \\textbf{n} \\ \\exp(N_{e}) \\ \\epsilon \\right] - \\frac{5}{3} \\epsilon_{\\gamma} (1 - a_{e}) \\gamma \\left[ \\mu_{+} \\ \\left( \\text{-} \\nabla (V / l_{c}) \\right) \\cdot \\textbf{n} \\ \\exp(N_{+}) \\right]", element, {displayMode:true,throwOnError:false});$

Where $var element = document.getElementById("moose-equation-634e0417-0810-4efb-bf49-4afcdf0166a4");katex.render("N_{j}", element, {displayMode:false,throwOnError:false});$ is the molar density of the species in log form and $var element = document.getElementById("moose-equation-5bc014e2-2bb8-43ca-8ed7-ae3536c02f42");katex.render("l_{c}", element, {displayMode:false,throwOnError:false});$ is the scaling factor of the mesh.

Example Input File Syntax

[BCs]
  [mean_en_physical_left]
    type = HagelaarEnergyAdvectionBC
    variable = mean_en
    boundary = 'left'
    potential = potential
    ions = Arp
    r = 0
    position_units = ${dom0Scale}
    secondary_electron_energy = 3
    emission_coeffs = 0.02
  []
[]

Input Parameters

boundaryThe list of boundary IDs from the mesh where this object applies
C++ Type:std::vector<BoundaryName>
Unit:(no unit assumed)
Controllable:No
Description:The list of boundary IDs from the mesh where this object applies
emission_coeffsThe species-dependent secondary electron emission coefficients for this boundary
C++ Type:std::vector<std::string>
Unit:(no unit assumed)
Controllable:No
Description:The species-dependent secondary electron emission coefficients for this boundary
ionsA list of ion densities in log form
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:A list of ion densities in log form
position_unitsUnits of position.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Units of position.
potentialThe electric potential
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:The electric potential
rThe reflection coefficient
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The reflection coefficient
secondary_electron_energyThe secondary electron energy in eV
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The secondary electron energy in eV
variableThe name of the variable that this residual object operates on
C++ Type:NonlinearVariableName
Unit:(no unit assumed)
Controllable:No
Description:The name of the variable that this residual object operates on

Required Parameters

displacementsThe displacements
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:The displacements
prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Optional Parameters

absolute_value_vector_tagsThe tags for the vectors this residual object should fill with the absolute value of the residual contribution
C++ Type:std::vector<TagName>
Unit:(no unit assumed)
Controllable:No
Description:The tags for the vectors this residual object should fill with the absolute value of the residual contribution
extra_matrix_tagsThe extra tags for the matrices this Kernel should fill
C++ Type:std::vector<TagName>
Unit:(no unit assumed)
Controllable:No
Description:The extra tags for the matrices this Kernel should fill
extra_vector_tagsThe extra tags for the vectors this Kernel should fill
C++ Type:std::vector<TagName>
Unit:(no unit assumed)
Controllable:No
Description:The extra tags for the vectors this Kernel should fill
matrix_tagssystemThe tag for the matrices this Kernel should fill
Default:system
C++ Type:MultiMooseEnum
Unit:(no unit assumed)
Options:nontime, system
Controllable:No
Description:The tag for the matrices this Kernel should fill
vector_tagsnontimeThe tag for the vectors this Kernel should fill
Default:nontime
C++ Type:MultiMooseEnum
Unit:(no unit assumed)
Options:nontime, time
Controllable:No
Description:The tag for the vectors this Kernel should fill

Tagging Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Unit:(no unit assumed)
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
diag_save_inThe name of auxiliary variables to save this BC's diagonal jacobian contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
C++ Type:std::vector<AuxVariableName>
Unit:(no unit assumed)
Controllable:No
Description:The name of auxiliary variables to save this BC's diagonal jacobian contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Unit:(no unit assumed)
Controllable:Yes
Description:Set the enabled status of the MooseObject.
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:Determines whether this object is calculated using an implicit or explicit form
save_inThe name of auxiliary variables to save this BC's residual contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
C++ Type:std::vector<AuxVariableName>
Unit:(no unit assumed)
Controllable:No
Description:The name of auxiliary variables to save this BC's residual contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Unit:(no unit assumed)
Controllable:No
Description:The seed for the master random number generator
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

Input Files

(test/tests/reflections/Schottky_300_V_5_um/Input.i)

References

GJM Hagelaar, FJ De Hoog, and GMW Kroesen. Boundary conditions in fluid models of gas discharges. Physical Review E, 62(1):1452, 2000. doi:10.1103/PhysRevE.62.1452.

@article{hagelaar2000boundary,
    author = "Hagelaar, GJM and De Hoog, FJ and Kroesen, GMW",
    title = "Boundary conditions in fluid models of gas discharges",
    journal = "Physical Review E",
    volume = "62",
    number = "1",
    pages = "1452",
    year = "2000",
    publisher = "APS",
    doi = "10.1103/PhysRevE.62.1452"
}

(test/tests/reflections/Schottky_300_V_5_um/Input.i)

dom0Scale = 1
dom0Size = 5E-6 #m
vhigh = -400E-3 #kV

[GlobalParams]
  potential_units = kV
  #        potential_units = V
  use_moles = true
[]

[Mesh]
  [file]
    type = FileMeshGenerator
    file = 'Geometry.msh'
  []
  [add_left]
    type = SideSetsFromNormalsGenerator
    normals = '-1 0 0'
    new_boundary = 'left'
    input = file
  []
  [add_right]
    type = SideSetsFromNormalsGenerator
    normals = '1 0 0'
    new_boundary = 'right'
    input = add_left
  []
[]

[Problem]
  type = FEProblem
  #        kernel_coverage_check = false
[]

[Preconditioning]
  [smp]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient
  #       line_search = none
  end_time = 10E6

  steady_state_detection = 1
  steady_state_tolerance = 1E-15
  steady_state_start_time = 1E-6

  petsc_options = '-snes_converged_reason -snes_linesearch_monitor'
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_shift_type -pc_factor_shift_amount -ksp_type -snes_linesearch_minlambda'
  petsc_options_value = 'lu NONZERO 1.e-10 preonly 1e-3'
  nl_rel_tol = 1e-9
  nl_abs_tol = 1e-9

  dtmin = 1e-16
  dtmax = 0.1E-7
  nl_max_its = 100
  [TimeSteppers]
    [Adaptive]
      type = IterationAdaptiveDT
      cutback_factor = 0.4
      dt = 1e-13
      growth_factor = 1.2
      optimal_iterations = 100
    []
  []
[]

[Outputs]
  perf_graph = true
  print_linear_residuals = false
  [out]
    type = Exodus
    #               execute_on = 'final'
  []
[]

[Debug]
  show_var_residual_norms = true
[]

[UserObjects]
  [data_provider]
    type = ProvideMobility
    electrode_area = 5.02e-7 # Formerly 3.14e-6
    ballast_resist = 1e0
    e = 1.6e-19
  []
[]

[Kernels]
  ## Stabilization
  [Arp_log_stabilization]
    type = LogStabilizationMoles
    variable = Arp
    offset = 45
    block = 0
  []
  [em_log_stabilization]
    type = LogStabilizationMoles
    variable = em
    offset = 45
    block = 0
  []
  [mean_en_log_stabilization]
    type = LogStabilizationMoles
    variable = mean_en
    block = 0
    offset = 45
  []
  #       [mean_en_advection_stabilization]
  #               type = EFieldArtDiff
  #               variable = mean_en
  #               potential = potential
  #               block = 0
  #       []

  [em_time_deriv]
    type = ElectronTimeDerivative
    variable = em
    block = 0
  []
  [em_advection]
    type = EFieldAdvection
    variable = em
    potential = potential
    block = 0
    position_units = ${dom0Scale}
  []
  [em_diffusion]
    type = CoeffDiffusion
    variable = em
    block = 0
    position_units = ${dom0Scale}
  []
  [em_ionization]
    type = ElectronsFromIonization
    em = em
    variable = em
    potential = potential
    mean_en = mean_en
    block = 0
    position_units = ${dom0Scale}
  []

  [potential_diffusion_dom1]
    type = CoeffDiffusionLin
    variable = potential
    block = 0
    position_units = ${dom0Scale}
  []

  [Arp_charge_source]
    type = ChargeSourceMoles_KV
    variable = potential
    charged = Arp
    block = 0
  []
  [em_charge_source]
    type = ChargeSourceMoles_KV
    variable = potential
    charged = em
    block = 0
  []

  [Arp_time_deriv]
    type = ElectronTimeDerivative
    variable = Arp
    block = 0
  []
  [Arp_advection]
    type = EFieldAdvection
    variable = Arp
    potential = potential
    position_units = ${dom0Scale}
    block = 0
  []
  [Arp_diffusion]
    type = CoeffDiffusion
    variable = Arp
    block = 0
    position_units = ${dom0Scale}
  []
  [Arp_ionization]
    type = IonsFromIonization
    variable = Arp
    potential = potential
    em = em
    mean_en = mean_en
    block = 0
    position_units = ${dom0Scale}
  []

  [mean_en_time_deriv]
    type = ElectronTimeDerivative
    variable = mean_en
    block = 0
  []
  [mean_en_advection]
    type = EFieldAdvection
    variable = mean_en
    potential = potential
    block = 0
    position_units = ${dom0Scale}
  []
  [mean_en_diffusion]
    type = CoeffDiffusion
    variable = mean_en
    block = 0
    position_units = ${dom0Scale}
  []
  [mean_en_joule_heating]
    type = JouleHeating
    variable = mean_en
    potential = potential
    em = em
    block = 0
    position_units = ${dom0Scale}
  []
  [mean_en_ionization]
    type = ElectronEnergyLossFromIonization
    variable = mean_en
    potential = potential
    em = em
    block = 0
    position_units = ${dom0Scale}
  []
  [mean_en_elastic]
    type = ElectronEnergyLossFromElastic
    variable = mean_en
    potential = potential
    em = em
    block = 0
    position_units = ${dom0Scale}
  []
  [mean_en_excitation]
    type = ElectronEnergyLossFromExcitation
    variable = mean_en
    potential = potential
    em = em
    block = 0
    position_units = ${dom0Scale}
  []
[]

[Variables]
  [potential]
  []
  [em]
    block = 0
  []
  [Arp]
    block = 0
  []
  [mean_en]
    block = 0
  []

[]

[AuxVariables]
  [e_temp]
    block = 0
    order = CONSTANT
    family = MONOMIAL
  []
  [x]
    order = CONSTANT
    family = MONOMIAL
  []
  [x_node]
  []
  [rho]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [em_lin]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [Arp_lin]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [Efield]
    order = CONSTANT
    family = MONOMIAL
  []
  [Current_em]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [Current_Arp]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [tot_gas_current]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [EFieldAdvAux_em]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [DiffusiveFlux_em]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [PowerDep_em]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [PowerDep_Arp]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [ProcRate_el]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [ProcRate_ex]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [ProcRate_iz]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
[]

[AuxKernels]
  [PowerDep_em]
    type = ADPowerDep
    density_log = em
    potential = potential
    art_diff = false
    potential_units = kV
    variable = PowerDep_em
    position_units = ${dom0Scale}
    block = 0
  []
  [PowerDep_Arp]
    type = ADPowerDep
    density_log = Arp
    potential = potential
    art_diff = false
    potential_units = kV
    variable = PowerDep_Arp
    position_units = ${dom0Scale}
    block = 0
  []
  [ProcRate_el]
    type = ADProcRate
    em = em
    potential = potential
    proc = el
    variable = ProcRate_el
    position_units = ${dom0Scale}
    block = 0
  []
  [ProcRate_ex]
    type = ADProcRate
    em = em
    potential = potential
    proc = ex
    variable = ProcRate_ex
    position_units = ${dom0Scale}
    block = 0
  []
  [ProcRate_iz]
    type = ADProcRate
    em = em
    potential = potential
    proc = iz
    variable = ProcRate_iz
    position_units = ${dom0Scale}
    block = 0
  []
  [e_temp]
    type = ElectronTemperature
    variable = e_temp
    electron_density = em
    mean_en = mean_en
    block = 0
  []
  [x_g]
    type = Position
    variable = x
    position_units = ${dom0Scale}
    block = 0
  []
  [x_ng]
    type = Position
    variable = x_node
    position_units = ${dom0Scale}
    block = 0
  []
  [rho]
    type = ParsedAux
    variable = rho
    coupled_variables = 'em_lin Arp_lin'
    expression = 'Arp_lin - em_lin'
    execute_on = 'timestep_end'
    block = 0
  []
  [tot_gas_current]
    type = ParsedAux
    variable = tot_gas_current
    coupled_variables = 'Current_em Current_Arp'
    expression = 'Current_em + Current_Arp'
    execute_on = 'timestep_end'
    block = 0
  []
  [em_lin]
    type = Density
    #               convert_moles = true
    variable = em_lin
    density_log = em
    block = 0
  []
  [Arp_lin]
    type = Density
    #               convert_moles = true
    variable = Arp_lin
    density_log = Arp
    block = 0
  []
  [Efield_g]
    type = Efield
    component = 0
    potential = potential
    variable = Efield
    position_units = ${dom0Scale}
    block = 0
  []
  [Current_em]
    type = ADCurrent
    potential = potential
    density_log = em
    variable = Current_em
    art_diff = false
    block = 0
    position_units = ${dom0Scale}
  []
  [Current_Arp]
    type = ADCurrent
    potential = potential
    density_log = Arp
    variable = Current_Arp
    art_diff = false
    block = 0
    position_units = ${dom0Scale}
  []
  [EFieldAdvAux_em]
    type = ADEFieldAdvAux
    potential = potential
    density_log = em
    variable = EFieldAdvAux_em
    block = 0
    position_units = ${dom0Scale}
  []
  [DiffusiveFlux_em]
    type = ADDiffusiveFlux
    density_log = em
    variable = DiffusiveFlux_em
    block = 0
    position_units = ${dom0Scale}
  []
[]

[BCs]
  ## Potential boundary conditions ##
  [potential_left]
    type = NeumannCircuitVoltageMoles_KV
    variable = potential
    boundary = left
    function = potential_bc_func
    ions = Arp
    data_provider = data_provider
    electrons = em
    electron_energy = mean_en
    r = 0
    position_units = ${dom0Scale}
    emission_coeffs = 0.02
  []

  [potential_dirichlet_right]
    type = DirichletBC
    variable = potential
    boundary = right
    value = 0
  []

  ## Electron boundary conditions ##
  [Emission_left]
    type = SchottkyEmissionBC
    #               type = SecondaryElectronBC
    variable = em
    boundary = 'left'
    potential = potential
    ions = Arp
    electron_energy = mean_en
    r = 1
    position_units = ${dom0Scale}
    tau = 100E-9
    relax = true
    emission_coeffs = 0.02
  []

  # [em_physical_left]
  #       type = HagelaarElectronBC
  #       variable = em
  #       boundary = 'left'
  #       potential = potential
  #       electron_energy = mean_en
  #       r = 0
  #       position_units = ${dom0Scale}
  # []

  [em_physical_right]
    type = HagelaarElectronAdvectionBC
    variable = em
    boundary = right
    potential = potential
    r = 0
    position_units = ${dom0Scale}
  []

  ## Argon boundary conditions ##
  #       [Arp_physical_left_diffusion]
  #               type = HagelaarIonDiffusionBC
  #               variable = Arp
  #               boundary = 'left'
  #               r = 0
  #               position_units = ${dom0Scale}
  #       []
  [Arp_physical_left_advection]
    type = HagelaarIonAdvectionBC
    variable = Arp
    boundary = 'left'
    potential = potential
    r = 0
    position_units = ${dom0Scale}
  []

  #       [Arp_physical_right_diffusion]
  #               type = HagelaarIonDiffusionBC
  #               variable = Arp
  #               boundary = right
  #               r = 0
  #               position_units = ${dom0Scale}
  #       []
  [Arp_physical_right_advection]
    type = HagelaarIonAdvectionBC
    variable = Arp
    boundary = right
    potential = potential
    r = 0
    position_units = ${dom0Scale}
  []

  ## Mean energy boundary conditions ##
  [mean_en_physical_left]
    type = HagelaarEnergyAdvectionBC
    variable = mean_en
    boundary = 'left'
    potential = potential
    ions = Arp
    r = 0
    position_units = ${dom0Scale}
    secondary_electron_energy = 3
    emission_coeffs = 0.02
  []

  [mean_en_physical_right]
    type = HagelaarEnergyBC
    variable = mean_en
    boundary = right
    potential = potential
    electrons = em
    r = 0
    position_units = ${dom0Scale}
  []
[]

[ICs]
  [potential_ic]
    type = FunctionIC
    variable = potential
    function = potential_ic_func
  []

  [em_ic]
    type = ConstantIC
    variable = em
    value = -42
    block = 0
  []

  [Arp_ic]
    type = ConstantIC
    variable = Arp
    value = -45
    block = 0
  []

  [mean_en_ic]
    type = ConstantIC
    variable = mean_en
    value = -36
    block = 0
  []
[]

[Functions]
  [potential_bc_func]
    type = ParsedFunction
    symbol_names = 'VHigh'
    symbol_values = '${vhigh}'
    expression = 'VHigh'
  []
  [potential_ic_func]
    type = ParsedFunction
    expression = '-${vhigh} * (${dom0Size} - x) / ${dom0Size}'
  []
[]

[Materials]
  [gas_block]
    type = Gas
    interp_trans_coeffs = true
    interp_elastic_coeff = true
    ramp_trans_coeffs = false
    em = em
    potential = potential
    ip = Arp
    mean_en = mean_en
    user_se_coeff = 0.02
    user_work_function = 4.55 # eV
    user_field_enhancement = 55
    user_Richardson_coefficient = 80E4
    user_cathode_temperature = 1273
    property_tables_file = td_argon_mean_en.txt
    block = 0
  []
[]

(test/tests/reflections/Schottky_300_V_5_um/Input.i)

dom0Scale = 1
dom0Size = 5E-6 #m
vhigh = -400E-3 #kV

[GlobalParams]
  potential_units = kV
  #        potential_units = V
  use_moles = true
[]

[Mesh]
  [file]
    type = FileMeshGenerator
    file = 'Geometry.msh'
  []
  [add_left]
    type = SideSetsFromNormalsGenerator
    normals = '-1 0 0'
    new_boundary = 'left'
    input = file
  []
  [add_right]
    type = SideSetsFromNormalsGenerator
    normals = '1 0 0'
    new_boundary = 'right'
    input = add_left
  []
[]

[Problem]
  type = FEProblem
  #        kernel_coverage_check = false
[]

[Preconditioning]
  [smp]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient
  #       line_search = none
  end_time = 10E6

  steady_state_detection = 1
  steady_state_tolerance = 1E-15
  steady_state_start_time = 1E-6

  petsc_options = '-snes_converged_reason -snes_linesearch_monitor'
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_shift_type -pc_factor_shift_amount -ksp_type -snes_linesearch_minlambda'
  petsc_options_value = 'lu NONZERO 1.e-10 preonly 1e-3'
  nl_rel_tol = 1e-9
  nl_abs_tol = 1e-9

  dtmin = 1e-16
  dtmax = 0.1E-7
  nl_max_its = 100
  [TimeSteppers]
    [Adaptive]
      type = IterationAdaptiveDT
      cutback_factor = 0.4
      dt = 1e-13
      growth_factor = 1.2
      optimal_iterations = 100
    []
  []
[]

[Outputs]
  perf_graph = true
  print_linear_residuals = false
  [out]
    type = Exodus
    #               execute_on = 'final'
  []
[]

[Debug]
  show_var_residual_norms = true
[]

[UserObjects]
  [data_provider]
    type = ProvideMobility
    electrode_area = 5.02e-7 # Formerly 3.14e-6
    ballast_resist = 1e0
    e = 1.6e-19
  []
[]

[Kernels]
  ## Stabilization
  [Arp_log_stabilization]
    type = LogStabilizationMoles
    variable = Arp
    offset = 45
    block = 0
  []
  [em_log_stabilization]
    type = LogStabilizationMoles
    variable = em
    offset = 45
    block = 0
  []
  [mean_en_log_stabilization]
    type = LogStabilizationMoles
    variable = mean_en
    block = 0
    offset = 45
  []
  #       [mean_en_advection_stabilization]
  #               type = EFieldArtDiff
  #               variable = mean_en
  #               potential = potential
  #               block = 0
  #       []

  [em_time_deriv]
    type = ElectronTimeDerivative
    variable = em
    block = 0
  []
  [em_advection]
    type = EFieldAdvection
    variable = em
    potential = potential
    block = 0
    position_units = ${dom0Scale}
  []
  [em_diffusion]
    type = CoeffDiffusion
    variable = em
    block = 0
    position_units = ${dom0Scale}
  []
  [em_ionization]
    type = ElectronsFromIonization
    em = em
    variable = em
    potential = potential
    mean_en = mean_en
    block = 0
    position_units = ${dom0Scale}
  []

  [potential_diffusion_dom1]
    type = CoeffDiffusionLin
    variable = potential
    block = 0
    position_units = ${dom0Scale}
  []

  [Arp_charge_source]
    type = ChargeSourceMoles_KV
    variable = potential
    charged = Arp
    block = 0
  []
  [em_charge_source]
    type = ChargeSourceMoles_KV
    variable = potential
    charged = em
    block = 0
  []

  [Arp_time_deriv]
    type = ElectronTimeDerivative
    variable = Arp
    block = 0
  []
  [Arp_advection]
    type = EFieldAdvection
    variable = Arp
    potential = potential
    position_units = ${dom0Scale}
    block = 0
  []
  [Arp_diffusion]
    type = CoeffDiffusion
    variable = Arp
    block = 0
    position_units = ${dom0Scale}
  []
  [Arp_ionization]
    type = IonsFromIonization
    variable = Arp
    potential = potential
    em = em
    mean_en = mean_en
    block = 0
    position_units = ${dom0Scale}
  []

  [mean_en_time_deriv]
    type = ElectronTimeDerivative
    variable = mean_en
    block = 0
  []
  [mean_en_advection]
    type = EFieldAdvection
    variable = mean_en
    potential = potential
    block = 0
    position_units = ${dom0Scale}
  []
  [mean_en_diffusion]
    type = CoeffDiffusion
    variable = mean_en
    block = 0
    position_units = ${dom0Scale}
  []
  [mean_en_joule_heating]
    type = JouleHeating
    variable = mean_en
    potential = potential
    em = em
    block = 0
    position_units = ${dom0Scale}
  []
  [mean_en_ionization]
    type = ElectronEnergyLossFromIonization
    variable = mean_en
    potential = potential
    em = em
    block = 0
    position_units = ${dom0Scale}
  []
  [mean_en_elastic]
    type = ElectronEnergyLossFromElastic
    variable = mean_en
    potential = potential
    em = em
    block = 0
    position_units = ${dom0Scale}
  []
  [mean_en_excitation]
    type = ElectronEnergyLossFromExcitation
    variable = mean_en
    potential = potential
    em = em
    block = 0
    position_units = ${dom0Scale}
  []
[]

[Variables]
  [potential]
  []
  [em]
    block = 0
  []
  [Arp]
    block = 0
  []
  [mean_en]
    block = 0
  []

[]

[AuxVariables]
  [e_temp]
    block = 0
    order = CONSTANT
    family = MONOMIAL
  []
  [x]
    order = CONSTANT
    family = MONOMIAL
  []
  [x_node]
  []
  [rho]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [em_lin]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [Arp_lin]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [Efield]
    order = CONSTANT
    family = MONOMIAL
  []
  [Current_em]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [Current_Arp]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [tot_gas_current]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [EFieldAdvAux_em]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [DiffusiveFlux_em]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [PowerDep_em]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [PowerDep_Arp]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [ProcRate_el]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [ProcRate_ex]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [ProcRate_iz]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
[]

[AuxKernels]
  [PowerDep_em]
    type = ADPowerDep
    density_log = em
    potential = potential
    art_diff = false
    potential_units = kV
    variable = PowerDep_em
    position_units = ${dom0Scale}
    block = 0
  []
  [PowerDep_Arp]
    type = ADPowerDep
    density_log = Arp
    potential = potential
    art_diff = false
    potential_units = kV
    variable = PowerDep_Arp
    position_units = ${dom0Scale}
    block = 0
  []
  [ProcRate_el]
    type = ADProcRate
    em = em
    potential = potential
    proc = el
    variable = ProcRate_el
    position_units = ${dom0Scale}
    block = 0
  []
  [ProcRate_ex]
    type = ADProcRate
    em = em
    potential = potential
    proc = ex
    variable = ProcRate_ex
    position_units = ${dom0Scale}
    block = 0
  []
  [ProcRate_iz]
    type = ADProcRate
    em = em
    potential = potential
    proc = iz
    variable = ProcRate_iz
    position_units = ${dom0Scale}
    block = 0
  []
  [e_temp]
    type = ElectronTemperature
    variable = e_temp
    electron_density = em
    mean_en = mean_en
    block = 0
  []
  [x_g]
    type = Position
    variable = x
    position_units = ${dom0Scale}
    block = 0
  []
  [x_ng]
    type = Position
    variable = x_node
    position_units = ${dom0Scale}
    block = 0
  []
  [rho]
    type = ParsedAux
    variable = rho
    coupled_variables = 'em_lin Arp_lin'
    expression = 'Arp_lin - em_lin'
    execute_on = 'timestep_end'
    block = 0
  []
  [tot_gas_current]
    type = ParsedAux
    variable = tot_gas_current
    coupled_variables = 'Current_em Current_Arp'
    expression = 'Current_em + Current_Arp'
    execute_on = 'timestep_end'
    block = 0
  []
  [em_lin]
    type = Density
    #               convert_moles = true
    variable = em_lin
    density_log = em
    block = 0
  []
  [Arp_lin]
    type = Density
    #               convert_moles = true
    variable = Arp_lin
    density_log = Arp
    block = 0
  []
  [Efield_g]
    type = Efield
    component = 0
    potential = potential
    variable = Efield
    position_units = ${dom0Scale}
    block = 0
  []
  [Current_em]
    type = ADCurrent
    potential = potential
    density_log = em
    variable = Current_em
    art_diff = false
    block = 0
    position_units = ${dom0Scale}
  []
  [Current_Arp]
    type = ADCurrent
    potential = potential
    density_log = Arp
    variable = Current_Arp
    art_diff = false
    block = 0
    position_units = ${dom0Scale}
  []
  [EFieldAdvAux_em]
    type = ADEFieldAdvAux
    potential = potential
    density_log = em
    variable = EFieldAdvAux_em
    block = 0
    position_units = ${dom0Scale}
  []
  [DiffusiveFlux_em]
    type = ADDiffusiveFlux
    density_log = em
    variable = DiffusiveFlux_em
    block = 0
    position_units = ${dom0Scale}
  []
[]

[BCs]
  ## Potential boundary conditions ##
  [potential_left]
    type = NeumannCircuitVoltageMoles_KV
    variable = potential
    boundary = left
    function = potential_bc_func
    ions = Arp
    data_provider = data_provider
    electrons = em
    electron_energy = mean_en
    r = 0
    position_units = ${dom0Scale}
    emission_coeffs = 0.02
  []

  [potential_dirichlet_right]
    type = DirichletBC
    variable = potential
    boundary = right
    value = 0
  []

  ## Electron boundary conditions ##
  [Emission_left]
    type = SchottkyEmissionBC
    #               type = SecondaryElectronBC
    variable = em
    boundary = 'left'
    potential = potential
    ions = Arp
    electron_energy = mean_en
    r = 1
    position_units = ${dom0Scale}
    tau = 100E-9
    relax = true
    emission_coeffs = 0.02
  []

  # [em_physical_left]
  #       type = HagelaarElectronBC
  #       variable = em
  #       boundary = 'left'
  #       potential = potential
  #       electron_energy = mean_en
  #       r = 0
  #       position_units = ${dom0Scale}
  # []

  [em_physical_right]
    type = HagelaarElectronAdvectionBC
    variable = em
    boundary = right
    potential = potential
    r = 0
    position_units = ${dom0Scale}
  []

  ## Argon boundary conditions ##
  #       [Arp_physical_left_diffusion]
  #               type = HagelaarIonDiffusionBC
  #               variable = Arp
  #               boundary = 'left'
  #               r = 0
  #               position_units = ${dom0Scale}
  #       []
  [Arp_physical_left_advection]
    type = HagelaarIonAdvectionBC
    variable = Arp
    boundary = 'left'
    potential = potential
    r = 0
    position_units = ${dom0Scale}
  []

  #       [Arp_physical_right_diffusion]
  #               type = HagelaarIonDiffusionBC
  #               variable = Arp
  #               boundary = right
  #               r = 0
  #               position_units = ${dom0Scale}
  #       []
  [Arp_physical_right_advection]
    type = HagelaarIonAdvectionBC
    variable = Arp
    boundary = right
    potential = potential
    r = 0
    position_units = ${dom0Scale}
  []

  ## Mean energy boundary conditions ##
  [mean_en_physical_left]
    type = HagelaarEnergyAdvectionBC
    variable = mean_en
    boundary = 'left'
    potential = potential
    ions = Arp
    r = 0
    position_units = ${dom0Scale}
    secondary_electron_energy = 3
    emission_coeffs = 0.02
  []

  [mean_en_physical_right]
    type = HagelaarEnergyBC
    variable = mean_en
    boundary = right
    potential = potential
    electrons = em
    r = 0
    position_units = ${dom0Scale}
  []
[]

[ICs]
  [potential_ic]
    type = FunctionIC
    variable = potential
    function = potential_ic_func
  []

  [em_ic]
    type = ConstantIC
    variable = em
    value = -42
    block = 0
  []

  [Arp_ic]
    type = ConstantIC
    variable = Arp
    value = -45
    block = 0
  []

  [mean_en_ic]
    type = ConstantIC
    variable = mean_en
    value = -36
    block = 0
  []
[]

[Functions]
  [potential_bc_func]
    type = ParsedFunction
    symbol_names = 'VHigh'
    symbol_values = '${vhigh}'
    expression = 'VHigh'
  []
  [potential_ic_func]
    type = ParsedFunction
    expression = '-${vhigh} * (${dom0Size} - x) / ${dom0Size}'
  []
[]

[Materials]
  [gas_block]
    type = Gas
    interp_trans_coeffs = true
    interp_elastic_coeff = true
    ramp_trans_coeffs = false
    em = em
    potential = potential
    ip = Arp
    mean_en = mean_en
    user_se_coeff = 0.02
    user_work_function = 4.55 # eV
    user_field_enhancement = 55
    user_Richardson_coefficient = 80E4
    user_cathode_temperature = 1273
    property_tables_file = td_argon_mean_en.txt
    block = 0
  []
[]

Overview
Example Input File Syntax
Input Parameters
Input Files
References