MOXOxygenToMetalRatio

MOX oxygen to metal ratio material.

Description

This material has the goal of computing an oxygen-to-metal ratio when the oxygen redistribution kernel is activated (MOXOxygenDiffusion). Moreover it can be useful for other modules that currently utilize just an initial oxygen-to-metal ratio. The equations are from Lassmann (1987). $var element = document.getElementById("moose-equation-1affd272-49ed-468c-b4e7-3a0224c787ed");katex.render("\\begin{aligned} & \\text{Hyper-stoichiometric case:}\\\\ \\text{ }\\\\ \\end{aligned}", element, {displayMode:true,throwOnError:false});$ $var element = document.getElementById("moose-equation-ab0018dd-68c7-40a4-b83a-d57e3a8bd6f5");katex.render("\\begin{aligned} O/M &= 2 + c_{i}\\\\ \\end{aligned}", element, {displayMode:true,throwOnError:false});$ where $var element = document.getElementById("moose-equation-c3422414-c1ad-4dc5-a832-e9c30c3727ff");katex.render("c_{i}", element, {displayMode:false,throwOnError:false});$ is the atomic fraction of interstitial atoms of oxygen. $var element = document.getElementById("moose-equation-a1cbfc48-3644-42bc-8d16-ab445a2f8481");katex.render("\\begin{aligned} & \\text{Hypo-stoichiometric case:}\\\\ \\text{ }\\\\ \\end{aligned}", element, {displayMode:true,throwOnError:false});$ $var element = document.getElementById("moose-equation-e876cbcb-18a4-4e06-90cf-ea615292f9cb");katex.render("\\begin{aligned} O/M &= 2 - 2c_{v}\\\\ \\end{aligned}", element, {displayMode:true,throwOnError:false});$ where $var element = document.getElementById("moose-equation-fac18207-7a8d-4629-adb9-d5880b910db2");katex.render("c_{v}", element, {displayMode:false,throwOnError:false});$ is the atomic fraction of oxygen vacancies.

For just the resolution of the oxygen diffusion equation (MOXOxygenDiffusion) the indices $var element = document.getElementById("moose-equation-4bcbf0d2-532e-4b69-ace4-b901ef0e2b96");katex.render("i", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-93687f0a-7837-48be-9ae1-a3d9583075ad");katex.render("v", element, {displayMode:false,throwOnError:false});$ are omitted like in Lassmann (1987).

Input Parameters

oxygenCoupled oxygen
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:Coupled oxygen

Required Parameters

O_M_initial1.98Initial oxygen to metal ratio
Default:1.98
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Initial oxygen to metal ratio
blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
boundaryThe list of boundaries (ids or names) from the mesh where this object applies
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The list of boundaries (ids or names) from the mesh where this object applies
computeTrueWhen false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the MaterialBase via MaterialBasePropertyInterface::getMaterialBase(). Non-computed MaterialBases are not sorted for dependencies.
Default:True
C++ Type:bool
Controllable:No
Description:When false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the MaterialBase via MaterialBasePropertyInterface::getMaterialBase(). Non-computed MaterialBases are not sorted for dependencies.
constant_onNONEWhen ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped
Default:NONE
C++ Type:MooseEnum
Options:NONE, ELEMENT, SUBDOMAIN
Controllable:No
Description:When ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped
declare_suffixAn optional suffix parameter that can be appended to any declared 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 declared properties. The suffix will be prepended with a '_' character.

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.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
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
Controllable:No
Description:Determines whether this object is calculated using an implicit or explicit form
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
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
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

output_propertiesList of material properties, from this material, to output (outputs must also be defined to an output type)
C++ Type:std::vector<std::string>
Controllable:No
Description:List of material properties, from this material, to output (outputs must also be defined to an output type)
outputsnone Vector of output names where you would like to restrict the output of variables(s) associated with this object
Default:none
C++ Type:std::vector<OutputName>
Controllable:No
Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object

Outputs Parameters

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
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.

Material Property Retrieval Parameters

Input Files

(assessment/MOX/JOYO/B14/PTM001/analysis/b14_ptm001_1D_sample3_ox_noAm.i)
(assessment/MOX/JOYO/B14/PTM001/analysis/b14_ptm001_1D_sample3_ox.i)
(test/tests/mox_oxygen_transport/mox_oxygen_transport_hyper.i)
(test/tests/mox_oxygen_transport/mox_oxygen_transport_hypo_V_Pu_low3.3.i)
(test/tests/mox_oxygen_to_metal_ratio/mox_OM_ratio_hypo.i)
(test/tests/mox_oxygen_transport/mox_oxygen_transport_hypo_V_Pu_up3.3.i)
(test/tests/mox_oxygen_to_metal_ratio/mox_OM_ratio_hyper.i)

References

K. Lassmann. The OXIRED model for redistribution of oxygen in nonstoichiometric uranium-plutonium oxides. Journal of Nuclear Materials, 150:10–16, 1987.

@ARTICLE{lassmann_1987_oxired,
    author = "Lassmann, K.",
    title = "The {OXIRED} model for redistribution of oxygen in nonstoichiometric uranium-plutonium oxides",
    journal = "Journal of Nuclear Materials",
    year = "1987",
    volume = "150",
    pages = "10-16"
}

(assessment/MOX/JOYO/B14/PTM001/analysis/b14_ptm001_1D_sample3_ox_noAm.i)

# Sample at midplane

initial_fuel_density = 11026.4

[GlobalParams]
  density = ${initial_fuel_density}
  family = LAGRANGE
[]
[Problem]
  type = ReferenceResidualProblem
  extra_tag_vectors = 'ref'
  reference_vector = 'ref'
[]
[Mesh]
  coord_type = RZ
  [layered1D_mesh]
    type = Layered1DMeshGenerator
    fuel_height = 10e-3
    pellet_outer_radius = 0.0027
    pellet_bottom_coor = 0.0
    pellet_mesh_density = customize
    nx_p = 200
    elem_type = EDGE2
    slices_per_block = 1
    include_plenum = false
    include_clad = false
  []
[]
[UserObjects]
  [pin_geometry]
    type = Layered1DFuelPinGeometry
    include_clad = false
    mesh_generator = layered1D_mesh
  []
[]
[Variables]
  [oxygen]
    initial_condition = 0.009 #nominal value
    scaling = 1e-20
  []
[]
[AuxVariables]
  [temp]
  []
  [fission_rate]
  []
  [burnup]
  []
  [oxygen_to_metal_ratio]
    order = CONSTANT
    family = MONOMIAL
    initial_condition = 1.982
  []
  [saved_t]
  []
  [saved_o]
  []
[]
[Functions]
  [power_history]
    type = PiecewiseLinear
    x = '-200 0  72000  158040  160200  246600  248400  249000.012  251280'
    y = '0 0  39456.16 39456.16  43890.71  43890.71  53442.06  53442.06  0'
  []
[]
[Kernels]
  [oxygen_time_derivative]
    type = TimeDerivative
    variable = oxygen
  []
  [oxygen]
    type = MOXOxygenDiffusion
    variable = oxygen
    temperature = temp
    burnup = burnup
    oxygen_to_metal_ratio = oxygen_to_metal_ratio
    O_M_initial = 1.982
    q0 = 0.3
    extra_vector_tags = 'ref'
  []
[]
[AuxKernels]
  [fission_rate]
    type = FissionRateGeneral
    fission_rate_formulation = MOX
    variable = fission_rate
    block = fuel
    initial_porosity = 0.1372
    rod_ave_lin_pow = power_history
    pellet_diameter = 0.0054
    execute_on = timestep_begin
    energy_per_fission = 3.2e-11
  []
  [burnup]
    type = BurnupAux
    block = fuel
    fission_rate = fission_rate
    variable = burnup
    execute_on = timestep_begin
  []
  [OM]
    type = MaterialRealAux
    variable = oxygen_to_metal_ratio
    property = oxygen_to_metal_ratio
    execute_on = timestep_end
  []
[]
[Materials]
  [fuel_thermal]
    type = MAMOXThermal
    temperature = temp
    porosity = 0.1372
    block = fuel
    Am_content = 0.0
    oxy_to_metal_ratio = 1.982
  []
  [fuel_density]
    type = ParsedMaterial
    block = fuel
    property_name = density
    expression = ${initial_fuel_density}
  []
  [O_M_ratio]
    type = MOXOxygenToMetalRatio
    oxygen = oxygen
    output_properties = 'oxygen_to_metal_ratio'
    O_M_initial = 1.982
  []
[]
[Preconditioning]
  [SMP]
    type = SMP
    full = true
  []
[]
[Dampers]
  [bound]
    type = BoundingValueNodalDamper
    max_value = 1
    min_value = 0
    variable = oxygen
  []
[]
[Executioner]
  type = Transient
  solve_type = 'PJFNK'
  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-pc_type -sub_pc_type'
  petsc_options_value = 'asm lu'
  line_search = 'none'
  l_max_its = 50
  l_tol = 8e-3
  nl_max_its = 50
  nl_rel_tol = 1e-5
  nl_abs_tol = 1e-5
  start_time = -200
  n_startup_steps = 1
  end_time = 251280
  dtmax = 10000
  dtmin = 0.25
  [TimeStepper]
    type = IterationAdaptiveDT
    dt = 1e1
    optimal_iterations = 15
    iteration_window = 2
    linear_iteration_ratio = 100
    growth_factor = 2
    cutback_factor = .5
    force_step_every_function_point = true
    timestep_limiting_function = power_history
  []
[]

[PerformanceMetricOutputs]
[]

[Outputs]
  perf_graph = true
  exodus = false
  color = true
  csv = false
  [console]
    type = Console
    max_rows = 25
  []
  [line_plot]
    type = CSV
    execute_on = 'TIMESTEP_END'
    time_step_interval =  1
    file_base = 1d
    create_final_symlink = true
  []
[]

(assessment/MOX/JOYO/B14/PTM001/analysis/b14_ptm001_1D_sample3_ox.i)

# Sample at midplane

initial_fuel_density = 11057.75

[GlobalParams]
  density = ${initial_fuel_density}
  family = LAGRANGE
[]
[Problem]
  type = ReferenceResidualProblem
  extra_tag_vectors = 'ref'
  reference_vector = 'ref'
[]
[Mesh]
  coord_type = RZ
  [layered1D_mesh]
    type = Layered1DMeshGenerator
    fuel_height = 10e-3
    pellet_outer_radius = 0.0027
    pellet_bottom_coor = 0.0
    pellet_mesh_density = customize
    nx_p = 200
    elem_type = EDGE2
    slices_per_block = 1
    include_plenum = false
    include_clad = false
  []
[]
[UserObjects]
  [pin_geometry]
    type = Layered1DFuelPinGeometry
    include_clad = false
    mesh_generator = layered1D_mesh
  []
[]
[Variables]
  [oxygen]
    initial_condition = 0.009 #nominal value
    scaling = 1e-20
  []
[]
[AuxVariables]
  [temp]
  []
  [fission_rate]
  []
  [burnup]
  []
  [oxygen_to_metal_ratio]
    order = CONSTANT
    family = MONOMIAL
    initial_condition = 1.982
  []
  [saved_t]
  []
  [saved_o]
  []
[]
[Functions]
  [power_history]
    type = PiecewiseLinear
    x = '-200 0  72000  158040  160200  246600  248400  249000.012  251280'
    y = '0 0  39456.16 39456.16  43890.71  43890.71  53442.06  53442.06  0'
  []
[]
[Kernels]
  [oxygen_time_derivative]
    type = TimeDerivative
    variable = oxygen
  []
  [oxygen]
    type = MOXOxygenDiffusion
    variable = oxygen
    temperature = temp
    burnup = burnup
    oxygen_to_metal_ratio = oxygen_to_metal_ratio
    O_M_initial = 1.982
    q0 = 0.31
    extra_vector_tags = 'ref'
  []
[]
[AuxKernels]
  [fission_rate]
    type = FissionRateGeneral
    fission_rate_formulation = MOX
    variable = fission_rate
    block = fuel
    initial_porosity = 0.1372
    rod_ave_lin_pow = power_history
    pellet_diameter = 0.0054
    execute_on = timestep_begin
    energy_per_fission = 3.2e-11
  []
  [burnup]
    type = BurnupAux
    block = fuel
    fission_rate = fission_rate
    variable = burnup
    execute_on = timestep_begin
  []
  [OM]
    type = MaterialRealAux
    variable = oxygen_to_metal_ratio
    property = oxygen_to_metal_ratio
    execute_on = timestep_end
  []
[]
[Materials]
  [fuel_thermal]
    type = MAMOXThermal
    temperature = temp
    porosity = 0.1372
    block = fuel
    Am_content = 0.0237
    oxy_to_metal_ratio = 1.982
  []
  [fuel_density]
    type = ParsedMaterial
    block = fuel
    property_name = density
    expression = ${initial_fuel_density}
  []
  [O_M_ratio]
    type = MOXOxygenToMetalRatio
    oxygen = oxygen
    output_properties = 'oxygen_to_metal_ratio'
    O_M_initial = 1.982
  []
[]
[Preconditioning]
  [SMP]
    type = SMP
    full = true
  []
[]
[Dampers]
  [bound]
    type = BoundingValueNodalDamper
    max_value = 1
    min_value = 0
    variable = oxygen
  []
[]
[Executioner]
  type = Transient
  solve_type = 'PJFNK'
  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-pc_type -sub_pc_type'
  petsc_options_value = 'asm lu'
  line_search = 'none'
  l_max_its = 50
  l_tol = 8e-3
  nl_max_its = 50
  nl_rel_tol = 1e-5
  nl_abs_tol = 1e-5
  start_time = -200
  n_startup_steps = 1
  end_time = 251280
  dtmax = 10000
  dtmin = 0.25
  [TimeStepper]
    type = IterationAdaptiveDT
    dt = 1e1
    optimal_iterations = 15
    iteration_window = 2
    linear_iteration_ratio = 100
    growth_factor = 2
    cutback_factor = .5
    force_step_every_function_point = true
    timestep_limiting_function = power_history
  []
[]

[PerformanceMetricOutputs]
[]

[Outputs]
  perf_graph = true
  exodus = false
  color = true
  csv = false
  [console]
    type = Console
    max_rows = 25
  []
  [line_plot]
    type = CSV
    execute_on = 'TIMESTEP_END'
    time_step_interval =  1
    file_base = 1d
    create_final_symlink = true
  []
[]

(test/tests/mox_oxygen_transport/mox_oxygen_transport_hyper.i)

# This test is designed to exercise the MOXOxygen diffuion kernel. The kernel contains Fickian
# and Soret diffusion terms and a diffusion coefficient that is a function of temperature.
# There doesn't seem to be a simple analytical solution to this pde, so only general observations
# can be made as to whether or not the result is reasonable. The solution should be a form of
# exponential. The solution's shape agrees with Fig.11.20 from "Fundamental ascpects of nuclear reactor fuel elements", by Olander.
# HYPERstochiometric case
[GlobalParams]
  density = 10431.0
  energy_per_fission = 3.2e-11 # J/fission
[]

[Problem]
  type = FEProblem
[]

[Mesh]
  coord_type = RZ
  [mesh]
    type = GeneratedMeshGenerator
    dim = 1
    xmax = 2.794e-3 #pellet radius
    nx = 100
    elem_type = EDGE
  []
[]

[Variables]
  [temp]
    initial_condition = 1400.0
  []
  [oxygen]
    order = FIRST
    family = LAGRANGE
    initial_condition = 0.02
    scaling = 1e-20
  []
[]
[AuxVariables]
  [fission_rate]
  []
  [burnup]
  []
[]
[Functions]
  [power_history]
    type = PiecewiseLinear
    x = '0 10000'
    y = '0 50000'
  []
[]
[Kernels]
  [heat]
    type = HeatConduction
    variable = temp
  []
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temp
  []
  [heat_source]
    type = NeutronHeatSource
    variable = temp
    fission_rate = fission_rate
  []
  [oxygen_time_derivative]
    type = TimeDerivative
    variable = oxygen
  []
  [oxygen]
    type = MOXOxygenDiffusion
    variable = oxygen
    temperature = temp
    burnup = burnup
    O_M_initial = 2.02
  []
[]

[AuxKernels]
  [fission_rate]
    type = FissionRateGeneral
    fission_rate_formulation = MOX
    variable = fission_rate
    rod_ave_lin_pow = power_history
    execute_on = timestep_begin
    pellet_diameter = 0.005588
  []
  [burnup]
    type = BurnupAux
    fission_rate = fission_rate
    variable = burnup
    execute_on = timestep_begin
  []
[]

[BCs]
  [temp_outside]
    type = DirichletBC
    variable = temp
    boundary = right
    value = 1400
  []
[]

[Materials]
  [fuel_thermal]
    type = MAMOXThermal
    temperature = temp
    oxy_to_metal_ratio = 2.02
    Am_content = 0.0
    Np_content = 0.0
    output_properties = 'thermal_conductivity'
  []
  [fuel_density]
    type = ParsedMaterial
    property_name = density
    expression = 10431.0
  []
  [O_M_ratio]
    type = MOXOxygenToMetalRatio
    oxygen = oxygen
    output_properties = 'oxygen_to_metal_ratio'
    O_M_initial = 2.02
    outputs = all
  []
[]

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

[Executioner]
  type = Transient

  solve_type = 'PJFNK'

  # The -snes_converged_reason option is used here to diagnose some
  # intermittent solver failures we have been seeing with this test
  petsc_options = '-snes_ksp_ew -snes_converged_reason'
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu       superlu_dist'

  line_search = 'none'

  l_tol = 8e-3
  nl_abs_tol = 1e-9
  nl_rel_tol = 1e-9
  n_startup_steps = 1
  end_time = 10000
  dt = 200
  dtmin = 200
[]

[Postprocessors] #Oxygen to metal ratio is an easier variable to look at for the user
  [ave_ox]
    type = ElementAverageValue
    variable = oxygen
  []
  [max_ox]
    type = NodalExtremeValue
    value_type = max
    variable = oxygen
  []
  [min_ox]
    type = NodalExtremeValue
    value_type = min
    variable = oxygen
  []
  [ave_om_ratio]
    type = ElementAverageValue
    variable = oxygen_to_metal_ratio
  []
  [max_om_ratio]
    type = ElementExtremeValue
    value_type = max
    variable = oxygen_to_metal_ratio
  []
  [min_om_ratio]
    type = ElementExtremeValue
    value_type = min
    variable = oxygen_to_metal_ratio
  []
[]
[VectorPostprocessors]
  [radial_oxygen]
    type = LineValueSampler
    variable = oxygen
    start_point = '0.0 0.0 0.0'
    end_point = '2.794e-3 0.0 0.0'
    num_points = 200
    execute_on = final
    sort_by = x
    outputs = line_plot
  []
  [radial_oxygen-to-metal-ratio]
    type = LineValueSampler
    variable = oxygen_to_metal_ratio
    start_point = '0.0 0.0 0.0'
    end_point = '2.794e-3 0.0 0.0'
    num_points = 200
    execute_on = final
    sort_by = x
    outputs = line_plot
  []
[]
[Outputs]
  exodus = false
  csv = true
  [console]
    type = Console
  []
  [line_plot]
    type = CSV
    execute_on = 'FINAL'
    time_step_interval =  1
  []
[]

[Debug]
  show_var_residual_norms = true
[]

(test/tests/mox_oxygen_transport/mox_oxygen_transport_hypo_V_Pu_low3.3.i)

# This test is designed to exercise the MOXOxygen diffuion kernel. The kernel contains Fickian
# and Soret diffusion terms and a diffusion coefficient that is a function of temperature.
# There doesn't seem to be a simple analytical solution to this pde, so only general observations
# can be made as to whether or not the result is reasonable. The solution should be a form of
# exponential. The solution's shape agrees with Fig.11.21 from "Fundamental ascpects of nuclear reactor fuel elements", by Olander.
# HYPOstochiometric case with a Plutonium Valence lower than 3.3
[GlobalParams]
  density = 10431.0
  energy_per_fission = 3.2e-11 # J/fission
[]

[Mesh]
  coord_type = RZ
  [mesh]
    type = GeneratedMeshGenerator
    dim = 1
    xmax = 2.794e-3 #pellet radius
    nx = 100
    elem_type = EDGE
  []
[]

[Variables]
  [temp]
    initial_condition = 1400.0
  []
  [oxygen]
    order = FIRST
    family = LAGRANGE
    initial_condition = 0.044
    scaling = 1e-20
  []
[]
[AuxVariables]
  [fission_rate]
  []
  [burnup]
  []
[]
[Functions]
  [power_history]
    type = PiecewiseLinear
    x = '0 10000'
    y = '0 50000'
  []
[]
[Kernels]
  [heat]
    type = HeatConduction
    variable = temp
  []
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temp
  []
  [heat_source]
    type = NeutronHeatSource
    variable = temp
    fission_rate = fission_rate
  []
  [oxygen_time_derivative]
    type = TimeDerivative
    variable = oxygen
  []
  [oxygen]
    type = MOXOxygenDiffusion
    variable = oxygen
    temperature = temp
    burnup = burnup
    oxygen_to_metal_ratio = oxygen_to_metal_ratio
    O_M_initial = 1.912
  []
[]

[AuxKernels]
  [fission_rate]
    type = FissionRateGeneral
    fission_rate_formulation = MOX
    variable = fission_rate
    rod_ave_lin_pow = power_history
    execute_on = timestep_begin
    pellet_diameter = 0.005588
  []
  [burnup]
    type = BurnupAux
    fission_rate = fission_rate
    variable = burnup
    execute_on = timestep_begin
  []
[]

[BCs]
  [temp_outside]
    type = DirichletBC
    variable = temp
    boundary = right
    value = 1400
  []
[]

[Materials]
  [fuel_thermal]
    type = MAMOXThermal
    temperature = temp
    oxy_to_metal_ratio = 1.912
    Am_content = 0.0
    Np_content = 0.0
    output_properties = 'thermal_conductivity'
  []
  [fuel_density]
    type = ParsedMaterial
    property_name = density
    expression = 10431.0
  []
  [O_M_ratio]
    type = MOXOxygenToMetalRatio
    oxygen = oxygen
    output_properties = 'oxygen_to_metal_ratio'
    O_M_initial = 1.912
    outputs = all
  []
[]

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

[Executioner]
  type = Transient

  solve_type = 'PJFNK'

  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu       superlu_dist'

  line_search = 'none'

  l_tol = 8e-3
  nl_abs_tol = 1e-9
  nl_rel_tol = 1e-9
  n_startup_steps = 1
  end_time = 10000
  dt = 200
  dtmin = 200
[]

[Postprocessors] #Oxygen to metal ratio is an easier variable to look at for the user
  [ave_ox]
    type = ElementAverageValue
    variable = oxygen
  []
  [max_ox]
    type = NodalExtremeValue
    value_type = max
    variable = oxygen
  []
  [min_ox]
    type = NodalExtremeValue
    value_type = min
    variable = oxygen
  []
  [ave_om_ratio]
    type = ElementAverageValue
    variable = oxygen_to_metal_ratio
  []
  [max_om_ratio]
    type = ElementExtremeValue
    value_type = max
    variable = oxygen_to_metal_ratio
  []
  [min_om_ratio]
    type = ElementExtremeValue
    value_type = min
    variable = oxygen_to_metal_ratio
  []
[]
[VectorPostprocessors]
  [radial_oxygen]
    type = LineValueSampler
    variable = oxygen
    start_point = '0.0 0.0 0.0'
    end_point = '2.794e-3 0.0 0.0'
    num_points = 200
    execute_on = final
    sort_by = x
    outputs = line_plot
  []
  [radial_oxygen-to-metal-ratio]
    type = LineValueSampler
    variable = oxygen_to_metal_ratio
    start_point = '0.0 0.0 0.0'
    end_point = '2.794e-3 0.0 0.0'
    num_points = 200
    execute_on = final
    sort_by = x
    outputs = line_plot
  []
[]
[Outputs]
  exodus = false
  csv = true
  [console]
    type = Console
  []
  [line_plot]
    type = CSV
    execute_on = 'FINAL'
    time_step_interval =  1
  []
[]

[Debug]
  show_var_residual_norms = true
[]

(test/tests/mox_oxygen_to_metal_ratio/mox_OM_ratio_hypo.i)

# This test is designed to exercise the MOXOxygenToMetalRatio Material.
# Here the HYPOstoichiometric case is tested
[GlobalParams]
[]

[Mesh]
  coord_type = RZ
  [mesh]
    type = GeneratedMeshGenerator
    dim            = 1
    xmax           = 2.794e-3 #pellet radius
    nx             = 100
    elem_type      = EDGE
  []
[]

[Variables]
  [oxygen]
  []
[]

[Kernels]
  [diff]
    type = Diffusion
    variable = oxygen
  []
[]

[BCs]
[ox_right]
  type = DirichletBC
  variable = oxygen
  boundary = right
  value = 0.0
[]
[ox_left]
  type = DirichletBC
  variable = oxygen
  boundary = left
  value = 0.05
[]
[]

[Materials]
  [O_M_ratio]
    type = MOXOxygenToMetalRatio
    oxygen = oxygen
    output_properties = 'oxygen_to_metal_ratio'
    outputs = all
  []
[]

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

[Executioner]
  type = Steady
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Postprocessors]
[ave_om_ratio]
  type = ElementAverageValue
  variable = oxygen_to_metal_ratio
[]
[max_om_ratio]
  type = ElementExtremeValue
  value_type = max
  variable = oxygen_to_metal_ratio
[]
[min_om_ratio]
  type = ElementExtremeValue
  value_type = min
  variable = oxygen_to_metal_ratio
[]
[]
[VectorPostprocessors]
  [radial_oxygen-to-metal-ratio]
    type = LineValueSampler
    variable = oxygen_to_metal_ratio
    start_point = '0.0 0.0 0.0'
    end_point = '2.794e-3 0.0 0.0'
    num_points = 200
    execute_on = final
    sort_by = x
    outputs = line_plot
  []
[]
[Outputs]
  exodus = false
  csv = true
  [console]
    type = Console
  []
  [line_plot]
    type = CSV
    execute_on = 'FINAL'
    time_step_interval =  1
  []
[]

[Debug]
 show_var_residual_norms = true
[]

(test/tests/mox_oxygen_transport/mox_oxygen_transport_hypo_V_Pu_up3.3.i)

# This test is designed to exercise the MOXOxygen diffuion kernel. The kernel contains Fickian
# and Soret diffusion terms and a diffusion coefficient that is a function of temperature.
# There doesn't seem to be a simple analytical solution to this pde, so only general observations
# can be made as to whether or not the result is reasonable. The solution should be a form of
# exponential. The solution's shape agrees with Fig.11.21 from "Fundamental ascpects of nuclear reactor fuel elements", by Olander.
# HYPOstochiometric case with a Plutonium Valence between 3.3 and 4
[GlobalParams]
  density = 10431.0
  energy_per_fission = 3.2e-11 # J/fission
[]

[Mesh]
  coord_type = RZ
  [mesh]
    type = GeneratedMeshGenerator
    dim = 1
    xmax = 2.794e-3 #pellet radius
    nx = 100
    elem_type = EDGE
  []
[]

[Variables]
  [temp]
    initial_condition = 1400.0
  []
  [oxygen]
    order = FIRST
    family = LAGRANGE
    initial_condition = 0.01
    scaling = 1e-20
  []
[]
[AuxVariables]
  [fission_rate]
  []
  [burnup]
  []
[]
[Functions]
  [power_history]
    type = PiecewiseLinear
    x = '0 10000'
    y = '0 50000'
  []
[]
[Kernels]
  [heat]
    type = HeatConduction
    variable = temp
  []
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temp
  []
  [heat_source]
    type = NeutronHeatSource
    variable = temp
    fission_rate = fission_rate
  []
  [oxygen_time_derivative]
    type = TimeDerivative
    variable = oxygen
  []
  [oxygen]
    type = MOXOxygenDiffusion
    variable = oxygen
    temperature = temp
    burnup = burnup
    oxygen_to_metal_ratio = oxygen_to_metal_ratio
  []
[]

[AuxKernels]
  [fission_rate]
    type = FissionRateGeneral
    fission_rate_formulation = MOX
    variable = fission_rate
    rod_ave_lin_pow = power_history
    execute_on = timestep_begin
    pellet_diameter = 0.005588
  []
  [burnup]
    type = BurnupAux
    fission_rate = fission_rate
    variable = burnup
    execute_on = timestep_begin
  []
[]

[BCs]
  [temp_outside]
    type = DirichletBC
    variable = temp
    boundary = right
    value = 1400
  []
[]

[Materials]
  [fuel_thermal]
    type = MAMOXThermal
    temperature = temp
    oxy_to_metal_ratio = 1.98
    Am_content = 0.0
    Np_content = 0.0
    output_properties = 'thermal_conductivity'
  []
  [fuel_density]
    type = ParsedMaterial
    property_name = density
    expression = 10431.0
  []
  [O_M_ratio]
    type = MOXOxygenToMetalRatio
    oxygen = oxygen
    output_properties = 'oxygen_to_metal_ratio'
    outputs = all
  []
[]

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

[Executioner]
  type = Transient

  solve_type = 'PJFNK'

  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu       superlu_dist'

  line_search = 'none'

  l_tol = 8e-3
  nl_abs_tol = 1e-9
  nl_rel_tol = 1e-9
  n_startup_steps = 1
  end_time = 10000
  dt = 200
  dtmin = 200
[]

[Postprocessors] #Oxygen to metal ratio is an easier variable to look at for the user
  [ave_ox]
    type = ElementAverageValue
    variable = oxygen
  []
  [max_ox]
    type = NodalExtremeValue
    value_type = max
    variable = oxygen
  []
  [min_ox]
    type = NodalExtremeValue
    value_type = min
    variable = oxygen
  []
  [ave_om_ratio]
    type = ElementAverageValue
    variable = oxygen_to_metal_ratio
  []
  [max_om_ratio]
    type = ElementExtremeValue
    value_type = max
    variable = oxygen_to_metal_ratio
  []
  [min_om_ratio]
    type = ElementExtremeValue
    value_type = min
    variable = oxygen_to_metal_ratio
  []
[]
[VectorPostprocessors]
  [radial_oxygen]
    type = LineValueSampler
    variable = oxygen
    start_point = '0.0 0.0 0.0'
    end_point = '2.794e-3 0.0 0.0'
    num_points = 200
    execute_on = final
    sort_by = x
    outputs = line_plot
  []
  [radial_oxygen-to-metal-ratio]
    type = LineValueSampler
    variable = oxygen_to_metal_ratio
    start_point = '0.0 0.0 0.0'
    end_point = '2.794e-3 0.0 0.0'
    num_points = 200
    execute_on = final
    sort_by = x
    outputs = line_plot
  []
[]
[Outputs]
  exodus = false
  csv = true
  [console]
    type = Console
  []
  [line_plot]
    type = CSV
    execute_on = 'FINAL'
    time_step_interval =  1
  []
[]

[Debug]
  show_var_residual_norms = true
[]

(test/tests/mox_oxygen_to_metal_ratio/mox_OM_ratio_hyper.i)

# This test is designed to exercise the MOXOxygenToMetalRatio Material.
# Here the HYPERstoichiometric case is tested
[GlobalParams]
[]

[Mesh]
  coord_type = RZ
  [mesh]
    type = GeneratedMeshGenerator
    dim            = 1
    xmax           = 2.794e-3 #pellet radius
    nx             = 100
    elem_type      = EDGE
  []
[]

[Variables]
  [oxygen]
  []
[]

[Kernels]
  [diff]
    type = Diffusion
    variable = oxygen
  []
[]

[BCs]
[ox_right]
  type = DirichletBC
  variable = oxygen
  boundary = right
  value = 0.0
[]
[ox_left]
  type = DirichletBC
  variable = oxygen
  boundary = left
  value = 0.03
[]
[]

[Materials]
  [O_M_ratio]
    type = MOXOxygenToMetalRatio
    oxygen = oxygen
    output_properties = 'oxygen_to_metal_ratio'
    O_M_initial = 2.02
    outputs = all
  []
[]

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

[Executioner]
  type = Steady
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Postprocessors]
[ave_om_ratio]
  type = ElementAverageValue
  variable = oxygen_to_metal_ratio
[]
[max_om_ratio]
  type = ElementExtremeValue
  value_type = max
  variable = oxygen_to_metal_ratio
[]
[min_om_ratio]
  type = ElementExtremeValue
  value_type = min
  variable = oxygen_to_metal_ratio
[]
[]
[VectorPostprocessors]
  [radial_oxygen-to-metal-ratio]
    type = LineValueSampler
    variable = oxygen_to_metal_ratio
    start_point = '0.0 0.0 0.0'
    end_point = '2.794e-3 0.0 0.0'
    num_points = 200
    execute_on = final
    sort_by = x
    outputs = line_plot
  []
[]
[Outputs]
  exodus = false
  csv = true
  [console]
    type = Console
  []
  [line_plot]
    type = CSV
    execute_on = 'FINAL'
    time_step_interval =  1
  []
[]

[Debug]
 show_var_residual_norms = true
[]

Description
Input Parameters
Input Files
References