U3Si2TricubicInterpolationUserObject

Performs tricubic interpolation in temperature, temperature gradient, and burnup (fission density) to determine the degradation to the thermal conductivity and gaseous swelling of U3Si2 fuel.

Description

The U3Si2TricubicInterpolationUserObject performs tricubic interpolation using the data contained in the input comma separated value (.csv) files for calculations utilized by SilicideFuelThermal and U3Si2VolumetricSwellingEigenstrain.

Required CSV Files

The user must specify six input csv files for the algorithm.

The first three represent the grid points of the tricubic calculation for which values are known. These grid points are for temperature, temperature gradient, and fission density. These files must be in ROW format consisting of a single row of data with each column delimited. The units of the values in the files should be K for temperature, K/mm for temperature gradient and fissions/cm $var element = document.getElementById("moose-equation-7074187a-cdb3-4614-88e1-ad1bdd60afe9");katex.render("^3", element, {displayMode:false,throwOnError:false});$ with the power of 10 $var element = document.getElementById("moose-equation-8f215278-d18f-409f-9fcc-219d17cec1d3");katex.render("^{21}", element, {displayMode:false,throwOnError:false});$ removed (i.e., a fission density of 2.0 $var element = document.getElementById("moose-equation-2a977509-b3d3-4478-ba9d-941d216adf94");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-9c22fc3a-5fc8-4cb3-87b6-9fddd3416eca");katex.render("^{21}", element, {displayMode:false,throwOnError:false});$ would be input as 2.0).

The last three files represent the values of grain boundary coverage (FCOV), gaseous swelling strain due to intragranular bubbles (GSWb), and total gaseous swelling (GSW) for a given combination of temperature, temperature gradient and fission density. These files must have as many rows as the product of the number of columns in the temperature_grid_points file and number of columns in the temperature_gradient_grid_points file. The number of columns in these three files should equal the number of columns in the fission_density_grid_points file.

Example Input Syntax

[UserObjects<<<{"href": "../../syntax/UserObjects/index.html"}>>>]
  [thermal_conductivity_degradation]
    type = U3Si2TricubicInterpolationUserObject<<<{"description": "Performs tricubic interpolation in temperature, temperature gradient, and burnup (fission density) to determine the degradation to the thermal conductivity and gaseous swelling of U3Si2 fuel.", "href": "U3Si2TricubicInterpolationUserObject.html"}>>>
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = 'initial nonlinear'
    temperature_grid_points_file<<<{"description": "File holding data for temperature grid points to be used with the tricubic interpolation."}>>> = temperature.csv
    temperature_gradient_grid_points_file<<<{"description": "File holding data for the temperature gradient grid points to be used with the tricubic interpolation."}>>> = temperature_gradient.csv
    fission_density_grid_points_file<<<{"description": "File holding data for the fission density (burnup) grid points to be used with the tricubic interpolation. Units are fissions/cm^3 with the trailing 10^21 removed (i.e., 2.5e20^21 would be 2.5)."}>>> = fission_density.csv
    grain_boundary_coverage_file<<<{"description": "File holding data for the grain boundary coverage (FCOV) to be used with the tricubic interpolation calculation."}>>> = fcov.csv
    intragranular_gaseous_swelling_file<<<{"description": "File holding data for the gaseous swelling strain due to intragranular bubbles (GSWb) to be used with the tricubic interpolation calculation."}>>> = gswb.csv
    total_gaseous_swelling_file<<<{"description": "File holding data for the total gaseous swelling strain (GSW) to be used with the tricubic interpolation calculation."}>>> = gsw.csv
  []
[]

Input Parameters

fission_density_grid_points_fileFile holding data for the fission density (burnup) grid points to be used with the tricubic interpolation. Units are fissions/cm^3 with the trailing 10^21 removed (i.e., 2.5e20^21 would be 2.5).
C++ Type:FileName
Controllable:No
Description:File holding data for the fission density (burnup) grid points to be used with the tricubic interpolation. Units are fissions/cm^3 with the trailing 10^21 removed (i.e., 2.5e20^21 would be 2.5).
grain_boundary_coverage_fileFile holding data for the grain boundary coverage (FCOV) to be used with the tricubic interpolation calculation.
C++ Type:FileName
Controllable:No
Description:File holding data for the grain boundary coverage (FCOV) to be used with the tricubic interpolation calculation.
intragranular_gaseous_swelling_fileFile holding data for the gaseous swelling strain due to intragranular bubbles (GSWb) to be used with the tricubic interpolation calculation.
C++ Type:FileName
Controllable:No
Description:File holding data for the gaseous swelling strain due to intragranular bubbles (GSWb) to be used with the tricubic interpolation calculation.
temperature_gradient_grid_points_fileFile holding data for the temperature gradient grid points to be used with the tricubic interpolation.
C++ Type:FileName
Controllable:No
Description:File holding data for the temperature gradient grid points to be used with the tricubic interpolation.
temperature_grid_points_fileFile holding data for temperature grid points to be used with the tricubic interpolation.
C++ Type:FileName
Controllable:No
Description:File holding data for temperature grid points to be used with the tricubic interpolation.
total_gaseous_swelling_fileFile holding data for the total gaseous swelling strain (GSW) to be used with the tricubic interpolation calculation.
C++ Type:FileName
Controllable:No
Description:File holding data for the total gaseous swelling strain (GSW) to be used with the tricubic interpolation calculation.

Required Parameters

allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Controllable:No
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:XFEM_MARK, NONE, INITIAL, LINEAR, NONLINEAR_CONVERGENCE, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM
Controllable:No
Description:The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
Default:0
C++ Type:int
Controllable:No
Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in POSTAUX
force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREAUX
force_preicFalseForces the UserObject to be executed in PREIC during initial setup
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREIC during initial setup

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

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

(test/tests/solid_mechanics/u3si2_eigenstrains/u3si2_vswelling/argonne_rz.i)
(test/tests/thermalSilicideFuel/thermalU3Si2_degradation.i)

(test/tests/thermalSilicideFuel/thermalU3Si2_degradation.i)

#
# This input file is used to test the thermal conductivity degradation
# model for U3Si2 fuel in an axisymmetric RZ geometry.
#
# The thermal conductivity model was developed using rate theory calculations at
# the Argonne National Laboratory:
#
# Y. Miao, K. A. Gamble, D. Andersson, B. Ye, Z. Mei, G. Hofman and A. M. Yacout,
# "Gaseous swelling of U3Si2 during steady-state LWR operation: A rate theory
# investigation," Nuclear Engineering and Design, 322, 2017, p. 336-344.
#
# Note that for this test the input CSV files are made up for testing the
# algorithm.  The full CSV data set developed using the rate theory calculations
# from Argonne National Laboratory can be obtained by contacted a BISON
# developer.
#
# For this problem, a 1x1 mm square is simulated. The temperature is fixed at
# 450 K and the burnup is fixed at 9.086425e-3 FIMA which corresponds to a
# fission density of 0.25e21 fissions/cm^3.
#
# The thermal conductivity degradation calculation is complex and outlined
# in the documentation.  Given the supplied input CSVs the parameters FCOV and GSWb
# calculated from the tricubic interpolation userobject should be 7.5 and 0.00375
# respectively.  Then, using these values and the given temperature and burnup
# above the analytically calculated thermal conductivity is given as: 9.79462 W/m-K
#
# This value is confirmed when looking at the avg_th_cond postprocessor.
#

[GlobalParams]
  displacements = 'disp_x disp_y'
  volumetric_locking_correction = true
[]

[Mesh]
  coord_type = RZ
  [mesh]
    type = GeneratedMeshGenerator
    dim = 2
    xmin = 0
    xmax = 1e-3
    ymin = 0
    ymax = 1e-3
    nx = 1
    ny = 1
  []
[]

[UserObjects]
  [thermal_conductivity_degradation]
    type = U3Si2TricubicInterpolationUserObject
    execute_on = 'initial nonlinear'
    temperature_grid_points_file = temperature.csv
    temperature_gradient_grid_points_file = temperature_gradient.csv
    fission_density_grid_points_file = fission_density.csv
    grain_boundary_coverage_file = fcov.csv
    intragranular_gaseous_swelling_file = gswb.csv
    total_gaseous_swelling_file = gsw.csv
  []
[]

[Variables]
  [temp]
    order = FIRST
    family = LAGRANGE
    initial_condition = 450.0
  []
[]

[AuxVariables]
  [th_cond]
    order = CONSTANT
    family = MONOMIAL
  []
  [burnup]
    order = FIRST
    family = LAGRANGE
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = FINITE
        block = 0
        incremental = true
        add_variables = true
      []
    []
  []
[]

[Kernels]
  [heat]
    type = HeatConduction
    variable = temp
  []
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temp
  []
[]

[Functions]
  [burnupFunc]
    type = ConstantFunction
    value = 9.086425e-3  # Corresponds to a fission density of 0.25e21 fissions/cm^3
  []
[]

[AuxKernels]
  [th_cond]
    type = MaterialRealAux
    variable = th_cond
    property = thermal_conductivity
    execute_on = 'initial timestep_end'
  []
  [burnup]
    type = FunctionAux
    variable = burnup
    function = burnupFunc
    execute_on = timestep_begin
  []
[]

[BCs]
  [left_x]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0.0
  []
  [bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0.0
  []
  [left_temp]
    type = DirichletBC
    variable = temp
    boundary = left
    value = 450.0
  []
  [right_temp]
    type = DirichletBC
    variable = temp
    boundary = right
    value = 450.0
  []
[]

[Materials]
  [stress]
    type = ComputeFiniteStrainElasticStress
  []
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = 0
    youngs_modulus = 1.0e11
    poissons_ratio = 0.17
  []
  [thermal]
    type = SilicideFuelThermal
    block = 0
    thermal_conductivity_degradation = 'thermal_conductivity_degradation'
    thermal_conductivity_model = WHITE
    burnup = burnup
    temperature = temp
  []
  [density]
    type = ParsedMaterial
    block = 0
    property_name = density
    expression = 11590.0
  []
[]

[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'
  l_max_its = 60
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-9
  l_tol = 1e-5

  start_time = 0.0
  end_time = 1.0
  num_steps = 1
[]

[Postprocessors]
  [ave_burnup]
    type = ElementAverageValue
    variable = burnup
    block = 0
  []
  [avg_th_cond]
    type = ElementAverageValue
    variable = th_cond
    execute_on = 'initial timestep_end'
  []
  [avg_temp]
    type = AverageNodalVariableValue
    variable = temp
    block = 0
  []
[]

[Outputs]
  exodus = true
[]

(test/tests/solid_mechanics/u3si2_eigenstrains/u3si2_vswelling/argonne_rz.i)

#
# This input file is used to test the Argonne gaseous swelling model for U3Si2 fuel
# in an axisymmetric RZ geometry.
#
# The gaseous swelling model was developed using rate theory calculations at
# the Argonne National Laboratory:
#
# Y. Miao, K. A. Gamble, D. Andersson, B. Ye, Z. Mei, G. Hofman and A. M. Yacout,
# "Gaseous swelling of U3Si2 during steady-state LWR operation: A rate theory
# investigation," Nuclear Engineering and Design, 322, 2017, p. 336-344.
#
# Note that for this test the input CSV files are made up for testing the
# algorithm.  The full CSV data set developed using the rate theory calculations
# from Argonne National Laboratory can be obtained by contacted a BISON
# developer.
#
# For this problem, a 1x1 mm square is simulated. The temperature is fixed at
# 450 K and the burnup is fixed at 9.086425e-3 FIMA which corresponds to a
# fission density of 0.25e21 fissions/cm^3.
#
# Given the input csv used in the U3Si2TricubicInterpolationUserObject the
# interpolated value for gaseous swelling should be 0.0075. Comparing to the
# value reported by the gas_swell postprocessor confirms the correct implementation
# of the tricubic interpolation scheme.
#

initial_fuel_density = 11590.0

[GlobalParams]
  density = ${initial_fuel_density}
  displacements = 'disp_x disp_y'
  volumetric_locking_correction = true
[]

[Mesh]
  coord_type = RZ
  [mesh]
    type = GeneratedMeshGenerator
    dim = 2
    xmin = 0
    xmax = 1e-3
    ymin = 0
    ymax = 1e-3
    nx = 1
    ny = 1
  []
[]

[Variables]
  [temp]
    order = FIRST
    family = LAGRANGE
    initial_condition = 450.0
  []
[]

[UserObjects]
  [gaseous_swelling]
    type = U3Si2TricubicInterpolationUserObject
    execute_on = 'initial nonlinear'
    temperature_grid_points_file = temperature.csv
    temperature_gradient_grid_points_file = temperature_gradient.csv
    fission_density_grid_points_file = fission_density.csv
    grain_boundary_coverage_file = fcov.csv
    intragranular_gaseous_swelling_file = gswb.csv
    total_gaseous_swelling_file = gsw.csv
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = FINITE
        block = 0
        incremental = true
        add_variables = true
        eigenstrain_names = 'swell'
      []
    []
  []
[]

[Kernels]
  [heat]         # gradient term in heat conduction equation
    type = HeatConduction
    variable = temp
  []

  [heat_ie]       # time term in heat conduction equation
    type = HeatConductionTimeDerivative
    variable = temp
  []
[]

[AuxVariables]
 [burnup]
    order = FIRST
    family = LAGRANGE
 []
 [gaseous_swelling]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [burnupFunc]
    type = ConstantFunction
    value = 9.086425e-3  # Corresponds to a fission density of 0.25e21 fissions/cm^3
  []
[]

[AuxKernels]
  [burnup]
    type = FunctionAux
    variable = burnup
    function = burnupFunc
    execute_on = timestep_begin
  []
  [gaseous_swelling]
    type = MaterialRealAux
    property = gaseous_swelling
    variable = gaseous_swelling
    block = 0
    execute_on = timestep_end
  []
[]

[BCs]
  [left_x]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0.0
  []
  [bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0.0
  []
  [left_temp]
    type = DirichletBC
    variable = temp
    boundary = left
    value = 450.0
  []
  [right_temp]
    type = DirichletBC
    variable = temp
    boundary = right
    value = 450.0
  []
[]

[Materials]
  [fuel_elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = 0
    youngs_modulus = 1e11
    poissons_ratio = 0.17
  []
  [fuel_swelling]
    type = U3Si2VolumetricSwellingEigenstrain
    block = 0
    temperature = temp
    burnup = burnup
    complete_burnup = 10
    gaseous_swelling_type = ARGONNE
    gaseous_swelling_object = gaseous_swelling
    eigenstrain_name = swell
    include_densification = false
    include_solid_swelling = false
  []
  [fuel_stress]
    type = ComputeFiniteStrainElasticStress
    block = 0
  []
  [thermal]
    type = SilicideFuelThermal
    block = 0
    temperature = temp
  []
  [density]
    type = StrainAdjustedDensity
    block = 0
    strain_free_density = ${initial_fuel_density}
  []
[]

[Executioner]
  type = Transient

  solve_type = 'PJFNK'

  petsc_options = '-snes_ksp_ew '
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
  petsc_options_value = 'lu       superlu_dist                  51'
  l_max_its = 60
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-10
  l_tol = 1e-5

  start_time = 0.0
  end_time = 1.0
  num_steps = 1
[]

[Postprocessors]
  [ave_burnup]
    type = ElementAverageValue
    variable = burnup
    block = 0
  []
  [gas_swell]
    type = ElementAverageValue
    variable = gaseous_swelling
    block = 0
  []
  [avg_temp]
    type = AverageNodalVariableValue
    variable = temp
    block = 0
  []
[]

[Outputs]
  exodus = true
[]

(test/tests/thermalSilicideFuel/thermalU3Si2_degradation.i)

#
# This input file is used to test the thermal conductivity degradation
# model for U3Si2 fuel in an axisymmetric RZ geometry.
#
# The thermal conductivity model was developed using rate theory calculations at
# the Argonne National Laboratory:
#
# Y. Miao, K. A. Gamble, D. Andersson, B. Ye, Z. Mei, G. Hofman and A. M. Yacout,
# "Gaseous swelling of U3Si2 during steady-state LWR operation: A rate theory
# investigation," Nuclear Engineering and Design, 322, 2017, p. 336-344.
#
# Note that for this test the input CSV files are made up for testing the
# algorithm.  The full CSV data set developed using the rate theory calculations
# from Argonne National Laboratory can be obtained by contacted a BISON
# developer.
#
# For this problem, a 1x1 mm square is simulated. The temperature is fixed at
# 450 K and the burnup is fixed at 9.086425e-3 FIMA which corresponds to a
# fission density of 0.25e21 fissions/cm^3.
#
# The thermal conductivity degradation calculation is complex and outlined
# in the documentation.  Given the supplied input CSVs the parameters FCOV and GSWb
# calculated from the tricubic interpolation userobject should be 7.5 and 0.00375
# respectively.  Then, using these values and the given temperature and burnup
# above the analytically calculated thermal conductivity is given as: 9.79462 W/m-K
#
# This value is confirmed when looking at the avg_th_cond postprocessor.
#

[GlobalParams]
  displacements = 'disp_x disp_y'
  volumetric_locking_correction = true
[]

[Mesh]
  coord_type = RZ
  [mesh]
    type = GeneratedMeshGenerator
    dim = 2
    xmin = 0
    xmax = 1e-3
    ymin = 0
    ymax = 1e-3
    nx = 1
    ny = 1
  []
[]

[UserObjects]
  [thermal_conductivity_degradation]
    type = U3Si2TricubicInterpolationUserObject
    execute_on = 'initial nonlinear'
    temperature_grid_points_file = temperature.csv
    temperature_gradient_grid_points_file = temperature_gradient.csv
    fission_density_grid_points_file = fission_density.csv
    grain_boundary_coverage_file = fcov.csv
    intragranular_gaseous_swelling_file = gswb.csv
    total_gaseous_swelling_file = gsw.csv
  []
[]

[Variables]
  [temp]
    order = FIRST
    family = LAGRANGE
    initial_condition = 450.0
  []
[]

[AuxVariables]
  [th_cond]
    order = CONSTANT
    family = MONOMIAL
  []
  [burnup]
    order = FIRST
    family = LAGRANGE
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = FINITE
        block = 0
        incremental = true
        add_variables = true
      []
    []
  []
[]

[Kernels]
  [heat]
    type = HeatConduction
    variable = temp
  []
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temp
  []
[]

[Functions]
  [burnupFunc]
    type = ConstantFunction
    value = 9.086425e-3  # Corresponds to a fission density of 0.25e21 fissions/cm^3
  []
[]

[AuxKernels]
  [th_cond]
    type = MaterialRealAux
    variable = th_cond
    property = thermal_conductivity
    execute_on = 'initial timestep_end'
  []
  [burnup]
    type = FunctionAux
    variable = burnup
    function = burnupFunc
    execute_on = timestep_begin
  []
[]

[BCs]
  [left_x]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0.0
  []
  [bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0.0
  []
  [left_temp]
    type = DirichletBC
    variable = temp
    boundary = left
    value = 450.0
  []
  [right_temp]
    type = DirichletBC
    variable = temp
    boundary = right
    value = 450.0
  []
[]

[Materials]
  [stress]
    type = ComputeFiniteStrainElasticStress
  []
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = 0
    youngs_modulus = 1.0e11
    poissons_ratio = 0.17
  []
  [thermal]
    type = SilicideFuelThermal
    block = 0
    thermal_conductivity_degradation = 'thermal_conductivity_degradation'
    thermal_conductivity_model = WHITE
    burnup = burnup
    temperature = temp
  []
  [density]
    type = ParsedMaterial
    block = 0
    property_name = density
    expression = 11590.0
  []
[]

[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'
  l_max_its = 60
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-9
  l_tol = 1e-5

  start_time = 0.0
  end_time = 1.0
  num_steps = 1
[]

[Postprocessors]
  [ave_burnup]
    type = ElementAverageValue
    variable = burnup
    block = 0
  []
  [avg_th_cond]
    type = ElementAverageValue
    variable = th_cond
    execute_on = 'initial timestep_end'
  []
  [avg_temp]
    type = AverageNodalVariableValue
    variable = temp
    block = 0
  []
[]

[Outputs]
  exodus = true
[]

Description
Example Input Syntax
Input Parameters
Input Files