AddDriftDiffusionAction

This Action automatically adds the necessary variables and objects to model the drift-diffusion equations within a domain block.

Overview

The DriftDiffusionAction object provides the necessary variables and objects to model the transport of plasma species using the drift-diffusion approximation.

The equation for the species density using the drift-diffusion approximation is:

$var element = document.getElementById("moose-equation-80ffb403-852c-4af1-aa2b-185c5d5c1678");katex.render("\\frac{d n_{j}}{d t} + \\nabla \\cdot \\left( \\text{sign}_{j} \\mu_{j} \\vec{E} n_{j} - D_{j} \\nabla n_{j} \\right) = S_{j}", element, {displayMode:true,throwOnError:false});$

Where:

the subscript $var element = document.getElementById("moose-equation-0e0b0348-d32b-4c16-b7a0-eed1e3ae24f3");katex.render("j", element, {displayMode:false,throwOnError:false});$ represents properties of the plasma species,
$var element = document.getElementById("moose-equation-84b4f967-f2d3-4aca-be8f-0aec851dedb8");katex.render("\\text{sign}", element, {displayMode:false,throwOnError:false});$ indicates the advection behavior ( $var element = document.getElementById("moose-equation-86e19a29-1740-4bd9-9062-f2c617d6103a");katex.render("\\text{+}1", element, {displayMode:false,throwOnError:false});$ for positively charged species, $var element = document.getElementById("moose-equation-7d9842a0-a0b8-4bc0-b695-50fc8e4e778c");katex.render("\\text{-}1", element, {displayMode:false,throwOnError:false});$ for negatively charged species, and $var element = document.getElementById("moose-equation-ec8bd9ac-ba1d-4334-989e-810eda0a6e48");katex.render("0", element, {displayMode:false,throwOnError:false});$ for neutral species),
$var element = document.getElementById("moose-equation-7ad9bf95-b505-4d0c-b119-69072e5d8c0f");katex.render("\\mu", element, {displayMode:false,throwOnError:false});$ is the mobility coefficient,
$var element = document.getElementById("moose-equation-c9a43127-a26f-44ea-81e6-e1f7ceda3bb9");katex.render("D", element, {displayMode:false,throwOnError:false});$ is the diffusion coefficient,
$var element = document.getElementById("moose-equation-3639e896-4745-438b-a805-542ebe003c8b");katex.render("n", element, {displayMode:false,throwOnError:false});$ is the density,
$var element = document.getElementById("moose-equation-ca167db5-bdcf-4d86-862e-d7c4f27c13e3");katex.render("\\vec{E}", element, {displayMode:false,throwOnError:false});$ is the electric field, and
$var element = document.getElementById("moose-equation-bc1ffc71-fd42-4fec-9312-d9801dcfd640");katex.render("S_{j}", element, {displayMode:false,throwOnError:false});$ is the source term.

In addition to the species density, Zapdos also tracks the mean energy density of the electrons using the equation:

$var element = document.getElementById("moose-equation-5171297a-54df-47f3-b9cc-8129e849e5d6");katex.render("\\frac{d n_{\\varepsilon}}{d t} + \\nabla \\cdot \\left( - \\mu_{\\varepsilon} \\vec{E} n_{\\varepsilon} - D_{\\varepsilon} \\nabla n_{\\varepsilon} \\right) - \\left( - \\mu_{e} \\vec{E} n_{e} - D_{e} \\nabla n_{e} \\right) \\cdot \\vec{E} = S_{\\varepsilon}", element, {displayMode:true,throwOnError:false});$

Where:

the subscript $var element = document.getElementById("moose-equation-50e9656b-d963-4791-9525-f0c11a904fa0");katex.render("e", element, {displayMode:false,throwOnError:false});$ represents properties of the electrons,
$var element = document.getElementById("moose-equation-d2a4da21-36ab-44fa-9d82-dfae267786c0");katex.render("n_{\\varepsilon}", element, {displayMode:false,throwOnError:false});$ is the electron mean energy density (defined as $var element = document.getElementById("moose-equation-27df8e32-3672-4166-a306-6e988710934f");katex.render("n_{\\varepsilon} = \\varepsilon n_{e}", element, {displayMode:false,throwOnError:false});$ )
$var element = document.getElementById("moose-equation-3dbab072-40c6-4640-a83d-41eaa7b4116b");katex.render("\\varepsilon", element, {displayMode:false,throwOnError:false});$ is the electron mean energy in units of eV,
$var element = document.getElementById("moose-equation-7afb4cdb-dc79-4a76-864a-70dd3fccf929");katex.render("\\mu_{\\varepsilon}", element, {displayMode:false,throwOnError:false});$ is the electron mean energy density mobility coefficient (defined as $var element = document.getElementById("moose-equation-baea183a-3e4b-4019-a9a6-fc313de7d988");katex.render("\\mu_{\\varepsilon} = \\frac{5}{3} \\mu_{e}", element, {displayMode:false,throwOnError:false});$ by default),
$var element = document.getElementById("moose-equation-8e8e4b14-a6f4-4e5b-a7f5-a009c0815b83");katex.render("D_{\\varepsilon}", element, {displayMode:false,throwOnError:false});$ is the electron mean energy density diffusion coefficient (defined as $var element = document.getElementById("moose-equation-47c65652-61e9-4f55-bd0f-71e63b3cf6a0");katex.render("D_{\\varepsilon} = \\frac{5}{3} D_{e}", element, {displayMode:false,throwOnError:false});$ by default), and
$var element = document.getElementById("moose-equation-8e342416-7686-4dc0-af56-1f987df1e058");katex.render("S_{\\varepsilon}", element, {displayMode:false,throwOnError:false});$ is the electron mean energy source term.

DriftDiffusionAction provides the objects for the right hand side of the above equations, while CRANE's AddZapdosReactions action is often used to supply the chemistry network which defines $var element = document.getElementById("moose-equation-cbed25c4-c1e4-4919-9912-da1d29084746");katex.render("S_{j}", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-a992da8a-5eb8-423c-8916-42af062e9e04");katex.render("S_{\\varepsilon}", element, {displayMode:false,throwOnError:false});$ .

If using the electrostatic approximation for the electric field, DriftDiffusionAction supplies the Poisson's equation, in the form of:

$var element = document.getElementById("moose-equation-4b701a65-fde3-4c60-bc15-b930e23a850e");katex.render("\\nabla \\cdot \\left( - \\varepsilon_{0} \\nabla V \\right) = e \\left( \\left[\\sum{ \\text{sign}_{i} n_{i} } \\right] - n_{e} \\right)", element, {displayMode:true,throwOnError:false});$

Where:

$var element = document.getElementById("moose-equation-5750c57c-58ba-452a-9766-6d2f0a274a75");katex.render("V", element, {displayMode:false,throwOnError:false});$ is the electrostatic potential, and
$var element = document.getElementById("moose-equation-a810d503-b46a-4787-b86c-f9ec3cacefa6");katex.render("\\varepsilon_{0}", element, {displayMode:false,throwOnError:false});$ is the permittivity of free space.

DriftDiffusionAction can supply additional outputs for electron temperature, current, and electric field by supplying the auxkernels ElectronTemperature, Current, and Efield respectively.

Example Input File Syntax

[DriftDiffusionAction<<<{"href": "../../syntax/DriftDiffusionAction/index.html"}>>>]
  [Plasma]
    electrons<<<{"description": "User given variable name for energy dependent electrons"}>>> = em
    charged_particle<<<{"description": "User given variable name for energy independent charged particle"}>>> = Ar+
    Neutrals<<<{"description": "The names of the neutrals that should be added"}>>> = Ar*
    mean_energy<<<{"description": "The gives the mean energy a variable name"}>>> = mean_en
    field<<<{"description": "The gives the field a variable name"}>>> = potential
    Is_field_unique<<<{"description": "Is the field unique to this block?If not, then the field variable should be defined in the Variable Block."}>>> = true
    using_offset<<<{"description": "Is the LogStabilizationMoles Kernel being used"}>>> = false
    position_units<<<{"description": "Units of position"}>>> = ${dom0Scale}
    Additional_Outputs<<<{"description": "Current list of available ouputs options in this action: Current, ElectronTemperature, EField"}>>> = 'ElectronTemperature Current EField'
  []
[]

Input Parameters

position_unitsUnits of position
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Units of position
potential_unitsUnits of potential (V or kV)
C++ Type:std::string
Controllable:No
Description:Units of potential (V or kV)
use_molesFalseWhether to convert from units of moles to \#.
Default:False
C++ Type:bool
Controllable:No
Description:Whether to convert from units of moles to \#.

Required Parameters

Additional_OutputsCurrent list of available ouputs options in this action: Current, ElectronTemperature, EField
C++ Type:std::vector<std::string>
Controllable:No
Description:Current list of available ouputs options in this action: Current, ElectronTemperature, EField
First_DriftDiffusionAction_in_blockTrueIs this the first DriftDiffusionAction for this block?If not, then the field solver kernels and Position auxkernel will NOT be supplied by this action.
Default:True
C++ Type:bool
Controllable:No
Description:Is this the first DriftDiffusionAction for this block?If not, then the field solver kernels and Position auxkernel will NOT be supplied by this action.
Is_field_uniqueFalseIs the field unique to this block?If not, then the field variable should be defined in the Variable Block.
Default:False
C++ Type:bool
Controllable:No
Description:Is the field unique to this block?If not, then the field variable should be defined in the Variable Block.
NeutralsThe names of the neutrals that should be added
C++ Type:std::vector<NonlinearVariableName>
Unit:(no unit assumed)
Controllable:No
Description:The names of the neutrals that should be added
active__all__ If specified only the blocks named will be visited and made active
Default:__all__
C++ Type:std::vector<std::string>
Controllable:No
Description:If specified only the blocks named will be visited and made active
blockThe subdomain that this action applies to.
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The subdomain that this action applies to.
charged_particleUser given variable name for energy independent charged particle
C++ Type:std::vector<NonlinearVariableName>
Unit:(no unit assumed)
Controllable:No
Description:User given variable name for energy independent charged particle
eff_fieldsThe effective fields that only affect their respective secondary charged particle.
C++ Type:std::vector<NonlinearVariableName>
Unit:(no unit assumed)
Controllable:No
Description:The effective fields that only affect their respective secondary charged particle.
eff_fields_property_namesName of the solver interface material property for the effective fields.
C++ Type:std::vector<std::string>
Controllable:No
Description:Name of the solver interface material property for the effective fields.
electronsUser given variable name for energy dependent electrons
C++ Type:NonlinearVariableName
Unit:(no unit assumed)
Controllable:No
Description:User given variable name for energy dependent electrons
familyLAGRANGESpecifies the family of FE shape functions to use for this variable
Default:LAGRANGE
C++ Type:MooseEnum
Options:LAGRANGE, MONOMIAL, HERMITE, SCALAR, HIERARCHIC, CLOUGH, XYZ, SZABAB, BERNSTEIN, L2_LAGRANGE, L2_HIERARCHIC, NEDELEC_ONE, LAGRANGE_VEC, MONOMIAL_VEC, RAVIART_THOMAS, RATIONAL_BERNSTEIN, SIDE_HIERARCHIC, L2_HIERARCHIC_VEC, L2_LAGRANGE_VEC, L2_RAVIART_THOMAS
Controllable:No
Description:Specifies the family of FE shape functions to use for this variable
fieldThe gives the field a variable name
C++ Type:std::string
Controllable:No
Description:The gives the field a variable name
field_property_namefield_solver_interface_propertyName of the solver interface material property.
Default:field_solver_interface_property
C++ Type:std::string
Controllable:No
Description:Name of the solver interface material property.
field_solverelectrostaticElectrostatic or electromagnetic field solver (default = electrostatic).
Default:electrostatic
C++ Type:MooseEnum
Options:electrostatic, electromagnetic
Controllable:No
Description:Electrostatic or electromagnetic field solver (default = electrostatic).
inactiveIf specified blocks matching these identifiers will be skipped.
C++ Type:std::vector<std::string>
Controllable:No
Description:If specified blocks matching these identifiers will be skipped.
initial_conditionSpecifies a constant initial condition for this variable
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Specifies a constant initial condition for this variable
initial_from_file_varGives the name of a variable for which to read an initial condition from a mesh file
C++ Type:std::string
Controllable:No
Description:Gives the name of a variable for which to read an initial condition from a mesh file
mean_energyThe gives the mean energy a variable name
C++ Type:NonlinearVariableName
Unit:(no unit assumed)
Controllable:No
Description:The gives the mean energy a variable name
offset20The offset parameter that goes into the exponential function
Default:20
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The offset parameter that goes into the exponential function
orderFIRSTSpecifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)
Default:FIRST
C++ Type:MooseEnum
Options:CONSTANT, FIRST, SECOND, THIRD, FOURTH
Controllable:No
Description:Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)
scalingSpecifies a scaling factor to apply to this variable
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Specifies a scaling factor to apply to this variable
secondary_charged_particlesThese are charged particles whose advection term in determined by an effective field.
C++ Type:std::vector<NonlinearVariableName>
Unit:(no unit assumed)
Controllable:No
Description:These are charged particles whose advection term in determined by an effective field.
using_offsetFalseIs the LogStabilizationMoles Kernel being used
Default:False
C++ Type:bool
Controllable:No
Description:Is the LogStabilizationMoles Kernel being used

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.

Advanced Parameters

(test/tests/DriftDiffusionAction/RF_Plasma_actions.i)

dom0Scale = 25.4e-3

[GlobalParams]
  potential_units = kV
  use_moles = true
[]

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

[Problem]
  type = FEProblem
[]

#Action the supplies the drift-diffusion equations
#This action also adds JouleHeating and the ChargeSourceMoles_KV Kernels
[DriftDiffusionAction]
  [Plasma]
    electrons = em
    charged_particle = Ar+
    Neutrals = Ar*
    mean_energy = mean_en
    field = potential
    Is_field_unique = true
    using_offset = false
    position_units = ${dom0Scale}
    Additional_Outputs = 'ElectronTemperature Current EField'
  []
[]

#The Kernels supply the sources terms
[Kernels]
  #Net electron production from ionization
  [em_ionization]
    type = EEDFReactionLog
    variable = em
    electrons = em
    target = Ar
    mean_energy = mean_en
    reaction = 'em + Ar -> em + em + Ar+'
    coefficient = 1
  []
  #Net electron production from step-wise ionization
  [em_stepwise_ionization]
    type = EEDFReactionLog
    variable = em
    electrons = em
    target = Ar*
    mean_energy = mean_en
    reaction = 'em + Ar* -> em + em + Ar+'
    coefficient = 1
  []
  #Net electron production from metastable pooling
  [em_pooling]
    type = ReactionSecondOrderLog
    variable = em
    v = Ar*
    w = Ar*
    reaction = 'Ar* + Ar* -> Ar+ + Ar + em'
    coefficient = 1
  []

  #Net ion production from ionization
  [Ar+_ionization]
    type = EEDFReactionLog
    variable = Ar+
    electrons = em
    target = Ar
    mean_energy = mean_en
    reaction = 'em + Ar -> em + em + Ar+'
    coefficient = 1
  []
  #Net ion production from step-wise ionization
  [Ar+_stepwise_ionization]
    type = EEDFReactionLog
    variable = Ar+
    electrons = em
    target = Ar*
    mean_energy = mean_en
    reaction = 'em + Ar* -> em + em + Ar+'
    coefficient = 1
  []
  #Net ion production from metastable pooling
  [Ar+_pooling]
    type = ReactionSecondOrderLog
    variable = Ar+
    v = Ar*
    w = Ar*
    reaction = 'Ar* + Ar* -> Ar+ + Ar + em'
    coefficient = 1
  []

  #Net excited Argon production from excitation
  [Ar*_excitation]
    type = EEDFReactionLog
    variable = Ar*
    electrons = em
    target = Ar
    mean_energy = mean_en
    reaction = 'em + Ar -> em + Ar*'
    coefficient = 1
  []
  #Net excited Argon loss from step-wise ionization
  [Ar*_stepwise_ionization]
    type = EEDFReactionLog
    variable = Ar*
    electrons = em
    target = Ar*
    mean_energy = mean_en
    reaction = 'em + Ar* -> em + em + Ar+'
    coefficient = -1
  []
  #Net excited Argon loss from superelastic collisions
  [Ar*_collisions]
    type = EEDFReactionLog
    variable = Ar*
    electrons = em
    target = Ar*
    mean_energy = mean_en
    reaction = 'em + Ar* -> em + Ar'
    coefficient = -1
  []
  #Net excited Argon loss from quenching to resonant
  [Ar*_quenching]
    type = EEDFReactionLog
    variable = Ar*
    electrons = em
    target = Ar*
    mean_energy = mean_en
    reaction = 'em + Ar* -> em + Ar_r'
    coefficient = -1
  []
  #Net excited Argon loss from  metastable pooling
  [Ar*_pooling]
    type = ReactionSecondOrderLog
    variable = Ar*
    v = Ar*
    w = Ar*
    reaction = 'Ar* + Ar* -> Ar+ + Ar + em'
    coefficient = -2
    _v_eq_u = true
    _w_eq_u = true
  []
  #Net excited Argon loss from two-body quenching
  [Ar*_2B_quenching]
    type = ReactionSecondOrderLog
    variable = Ar*
    v = Ar*
    w = Ar
    reaction = 'Ar* + Ar -> Ar + Ar'
    coefficient = -1
    _v_eq_u = true
  []
  #Net excited Argon loss from three-body quenching
  [Ar*_3B_quenching]
    type = ReactionThirdOrderLog
    variable = Ar*
    v = Ar*
    w = Ar
    x = Ar
    reaction = 'Ar* + Ar + Ar -> Ar_2 + Ar'
    coefficient = -1
    _v_eq_u = true
  []

  #Energy loss from ionization
  [Ionization_Loss]
    type = EEDFEnergyLog
    variable = mean_en
    electrons = em
    target = Ar
    reaction = 'em + Ar -> em + em + Ar+'
    threshold_energy = -15.7
  []
  #Energy loss from excitation
  [Excitation_Loss]
    type = EEDFEnergyLog
    variable = mean_en
    electrons = em
    target = Ar
    reaction = 'em + Ar -> em + Ar*'
    threshold_energy = -11.56
  []
  #Energy loss from step-wise ionization
  [Stepwise_Ionization_Loss]
    type = EEDFEnergyLog
    variable = mean_en
    electrons = em
    target = Ar*
    reaction = 'em + Ar* -> em + em + Ar+'
    threshold_energy = -4.14
  []
  #Energy gain from superelastic collisions
  [Collisions_Loss]
    type = EEDFEnergyLog
    variable = mean_en
    electrons = em
    target = Ar*
    reaction = 'em + Ar* -> em + Ar'
    threshold_energy = 11.56
  []
[]

[AuxVariables]
  [x_node]
  []

  [rho]
    order = CONSTANT
    family = MONOMIAL
  []

  [Ar]
  []

  [emRate]
    order = CONSTANT
    family = MONOMIAL
  []
  [exRate]
    order = CONSTANT
    family = MONOMIAL
  []
  [swRate]
    order = CONSTANT
    family = MONOMIAL
  []
  [deexRate]
    order = CONSTANT
    family = MONOMIAL
  []
  [quRate]
    order = CONSTANT
    family = MONOMIAL
  []
  [poolRate]
    order = CONSTANT
    family = MONOMIAL
  []
  [TwoBRate]
    order = CONSTANT
    family = MONOMIAL
  []
  [ThreeBRate]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[AuxKernels]
  [emRate]
    type = ProcRateForRateCoeff
    variable = emRate
    v = em
    w = Ar
    reaction = 'em + Ar -> em + em + Ar+'
  []
  [exRate]
    type = ProcRateForRateCoeff
    variable = exRate
    v = em
    w = Ar*
    reaction = 'em + Ar -> em + Ar*'
  []
  [swRate]
    type = ProcRateForRateCoeff
    variable = swRate
    v = em
    w = Ar*
    reaction = 'em + Ar* -> em + em + Ar+'
  []
  [deexRate]
    type = ProcRateForRateCoeff
    variable = deexRate
    v = em
    w = Ar*
    reaction = 'em + Ar* -> em + Ar'
  []
  [quRate]
    type = ProcRateForRateCoeff
    variable = quRate
    v = em
    w = Ar*
    reaction = 'em + Ar* -> em + Ar_r'
  []
  [poolRate]
    type = ProcRateForRateCoeff
    variable = poolRate
    v = Ar*
    w = Ar*
    reaction = 'Ar* + Ar* -> Ar+ + Ar + em'
  []
  [TwoBRate]
    type = ProcRateForRateCoeff
    variable = TwoBRate
    v = Ar*
    w = Ar
    reaction = 'Ar* + Ar -> Ar + Ar'
  []
  [ThreeBRate]
    type = ProcRateForRateCoeffThreeBody
    variable = ThreeBRate
    v = Ar*
    w = Ar
    x = Ar
    reaction = 'Ar* + Ar + Ar -> Ar_2 + Ar'
  []

  [x_ng]
    type = Position
    variable = x_node
    position_units = ${dom0Scale}
  []

  [Ar_val]
    type = ConstantAux
    variable = Ar
    # value = 3.22e22
    value = -2.928623
    execute_on = INITIAL
  []
[]

[BCs]
  #Voltage Boundary Condition, same as in paper
  [potential_left]
    type = FunctionDirichletBC
    variable = potential
    boundary = 'left'
    function = potential_bc_func
    preset = false
  []
  [potential_dirichlet_right]
    type = DirichletBC
    variable = potential
    boundary = 'right'
    value = 0
    preset = false
  []

  #New Boundary conditions for electons, same as in paper
  [em_physical_right]
    type = LymberopoulosElectronBC
    variable = em
    boundary = 'right'
    emission_coeffs = 0.01
    #emission_coeffs = 1
    ks = 1.19e5
    #ks = 0.0
    ions = Ar+
    position_units = ${dom0Scale}
  []
  [em_physical_left]
    type = LymberopoulosElectronBC
    variable = em
    boundary = 'left'
    emission_coeffs = 0.01
    #emission_coeffs = 1
    ks = 1.19e5
    #ks = 0.0
    ions = Ar+
    position_units = ${dom0Scale}
  []

  #New Boundary conditions for ions, should be the same as in paper
  [Ar+_physical_right_advection]
    type = LymberopoulosIonBC
    variable = Ar+
    boundary = 'right'
    position_units = ${dom0Scale}
  []
  [Ar+_physical_left_advection]
    type = LymberopoulosIonBC
    variable = Ar+
    boundary = 'left'
    position_units = ${dom0Scale}
  []

  #New Boundary conditions for ions, should be the same as in paper
  #(except the metastables are not set to zero, since Zapdos uses log form)
  [Ar*_physical_right_diffusion]
    type = LogDensityDirichletBC
    variable = Ar*
    boundary = 'right'
    value = 100
  []
  [Ar*_physical_left_diffusion]
    type = LogDensityDirichletBC
    variable = Ar*
    boundary = 'left'
    value = 100
  []

  #New Boundary conditions for mean energy, should be the same as in paper
  [mean_en_physical_right]
    type = ElectronTemperatureDirichletBC
    variable = mean_en
    electrons = em
    value = 0.5
    boundary = 'right'
  []
  [mean_en_physical_left]
    type = ElectronTemperatureDirichletBC
    variable = mean_en
    electrons = em
    value = 0.5
    boundary = 'left'
  []

[]

[ICs]
  [em_ic]
    type = FunctionIC
    variable = em
    function = density_ic_func
  []
  [Ar+_ic]
    type = FunctionIC
    variable = Ar+
    function = density_ic_func
  []
  [Ar*_ic]
    type = FunctionIC
    variable = Ar*
    function = density_ic_func
  []
  [mean_en_ic]
    type = FunctionIC
    variable = mean_en
    function = energy_density_ic_func
  []

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

[Functions]
  [potential_bc_func]
    type = ParsedFunction
    expression = '0.100*sin(2*3.1415926*13.56e6*t)'
  []
  [potential_ic_func]
    type = ParsedFunction
    expression = '0.100 * (25.4e-3 - x)'
  []
  [density_ic_func]
    type = ParsedFunction
    expression = 'log((1e13 + 1e15 * (1-x/1)^2 * (x/1)^2)/6.022e23)'
  []
  [energy_density_ic_func]
    type = ParsedFunction
    expression = 'log(3./2.) + log((1e13 + 1e15 * (1-x/1)^2 * (x/1)^2)/6.022e23)'
  []
[]

[Materials]
  [GasBasics]
    type = GasElectronMoments
    interp_trans_coeffs = false
    interp_elastic_coeff = false
    ramp_trans_coeffs = false
    user_p_gas = 133.322
    em = em
    mean_en = mean_en
    user_electron_mobility = 30.0
    user_electron_diffusion_coeff = 119.8757763975
    property_tables_file = Argon_reactions_paper_RateCoefficients/electron_moments.txt
  []
  [gas_species_0]
    type = ADHeavySpecies
    heavy_species_name = Ar+
    heavy_species_mass = 6.64e-26
    heavy_species_charge = 1.0
    mobility = 0.144409938
    diffusivity = 6.428571e-3
  []
  [gas_species_1]
    type = ADHeavySpecies
    heavy_species_name = Ar*
    heavy_species_mass = 6.64e-26
    heavy_species_charge = 0.0
    diffusivity = 7.515528e-3
  []
  [gas_species_2]
    type = ADHeavySpecies
    heavy_species_name = Ar
    heavy_species_mass = 6.64e-26
    heavy_species_charge = 0.0
  []
  [reaction_0]
    type = ZapdosEEDFRateConstant
    mean_energy = mean_en
    property_file = 'Argon_reactions_paper_RateCoefficients/ar_excitation.txt'
    reaction = 'em + Ar -> em + Ar*'
    electrons = em
  []
  [reaction_1]
    type = ZapdosEEDFRateConstant
    mean_energy = mean_en
    property_file = 'Argon_reactions_paper_RateCoefficients/ar_ionization.txt'
    reaction = 'em + Ar -> em + em + Ar+'
    electrons = em
  []
  [reaction_2]
    type = ZapdosEEDFRateConstant
    mean_energy = mean_en
    property_file = 'Argon_reactions_paper_RateCoefficients/ar_deexcitation.txt'
    reaction = 'em + Ar* -> em + Ar'
    electrons = em
  []
  [reaction_3]
    type = ZapdosEEDFRateConstant
    mean_energy = mean_en
    property_file = 'Argon_reactions_paper_RateCoefficients/ar_excited_ionization.txt'
    reaction = 'em + Ar* -> em + em + Ar+'
    electrons = em
  []
  [reaction_4]
    type = GenericRateConstant
    reaction = 'em + Ar* -> em + Ar_r'
    #reaction_rate_value = 2e-13
    reaction_rate_value = 1.2044e11
  []
  [reaction_5]
    type = GenericRateConstant
    reaction = 'Ar* + Ar* -> Ar+ + Ar + em'
    #reaction_rate_value = 6.2e-16
    reaction_rate_value = 373364000
  []
  [reaction_6]
    type = GenericRateConstant
    reaction = 'Ar* + Ar -> Ar + Ar'
    #reaction_rate_value = 3e-21
    reaction_rate_value = 1806.6
  []
  [reaction_7]
    type = GenericRateConstant
    reaction = 'Ar* + Ar + Ar -> Ar_2 + Ar'
    #reaction_rate_value = 1.1e-42
    reaction_rate_value = 398909.324
  []
[]

#New postprocessor that calculates the inverse of the plasma frequency
[Postprocessors]
  [InversePlasmaFreq]
    type = PlasmaFrequencyInverse
    variable = em
    use_moles = true
    execute_on = 'INITIAL TIMESTEP_BEGIN'
  []
[]

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

  [fdp]
    type = FDP
    full = true
  []
[]

[Executioner]
  type = Transient
  end_time = 3e-7
  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 fgmres 1e-3'

  dtmin = 1e-14
  l_max_its = 20

  scheme = bdf2
  dt = 1e-9
[]

[Outputs]
  file_base = 'RF_out'
  perf_graph = true
  [out]
    type = Exodus
  []
[]

Overview
Example Input File Syntax
Input Parameters