Standard Metallic Fuel Rod Outputs Action System

Description

The StandardMetallicFuelRodOutputs action simplifies the input file of a metallic fuel simulation by creating up to 11 postprocessors, vector postprocessors, and user objects which are commonly used to analyze the results of Bison simulations. The standard names given to these postprocessors and vector postprocessors by this action aid in the generation of plots with a common plotting tool.

warning:For Metallic Fuel Use Only

This action is designed for use with only Metallic Fuel simulations run with a mesh generated by either the mesh script or the FuelPinMeshGenerator; this action relies on the use of the boundary condition naming conventions in these files.

The user may elect to use the StandardMetallicFuelRodOutputs action to generate output for only the cladding, only the fuel pellets, or both as appropriate to the simulation. All of the possible outputs, generated with the default parameter setting rod_component = Both, are shown in Table 1.

Table 1: Correspondence Among Action Functionality and MooseObjects for the StandardMetallicFuelRodOutputs action

Functionality	Generated Postprocessor Name	Created Classes
Average temperature of the cladding interior	`average_interior_clad_temperature`	SideAverageValue
Average centerline fuel pellet temperature	`average_centerline_fuel_temperature`	AxisymmetricCenterlineAverageValue
Average plenum temperature	`plenum_temperature`	SideAverageValue
Volume of the plenum gases	`plenum_volume`	InternalVolume
Volume of the fuel pellets	`pellet_volume`	InternalVolume
Maximum axial displacement of the clad interior	`maximum_clad_elongation`	NodalExtremeValue
Maximum axial displacement of the fuel pellets	`maximum_fuel_elongation`	NodalExtremeValue
Radial displacement of the outer cladding at the simulation end	input file name + `_clad_radial_displacement_FINAL`	NodalValueSampler
Radial displacement of the outer fuel pellet at the simulation end	input file name + `_fuel_outer_radial_displacement_FINAL`	NodalValueSampler
Max Total Hoop Strain	`max_total_hoop_strain`	ElementExtremeValue
Cladding Total Hoop Strain	`clad_total_hoop_strain`	LineValueSampler
Generated fission gas	`fission_gas_generated`	ElementIntegralMaterialProperty
Released fission gas	`fission_gas_released`	ElementIntegralMaterialProperty
Percentage of released fission gas	`fission_gas_released_percentage`	FGRPercent
Rod average burnup	`average_burnup`	ElementAverageValue

note:Plenum Pressure

Additionally, this action can create the PlenumPressure by including the initial_pressure. For simplicity, the gas constant for this situation is assumed to be R = 8.3143 which is in units of J $var element = document.getElementById("moose-equation-a4f30d8f-22bb-4dbd-8f17-cad040fccb11");katex.render("\\cdot", element, {displayMode:false,throwOnError:false});$ K $var element = document.getElementById("moose-equation-d2036ac9-980d-4372-a41f-ba960a2aa543");katex.render("^{-1}\\cdot", element, {displayMode:false,throwOnError:false});$ mol $var element = document.getElementById("moose-equation-a1caf35b-d733-40a5-af67-38d141fb1be8");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ . Additional unmeshed volumes and their corresponding temperatures can be included when calculating plenum pressure, e.g. the volume between dish and chamfer surfaces of adjacent pellets by using the parameters additional_volumes and temperature_of_additional_volumes.

note:Volume of the Plenum Gas

In calculation the volume of the plenum gas, additional volumes can be included via the parameter internal_volume_addition. This can apply to situations such as the addition of sodium within the internal volume.

note

The action accommodates different naming conventions. The amount of fission gas produced can be called fis_gas_prod or fgm_produced, while the amount of fission gas released can be named either fis_gas_rel or fgm_released.

Example Input File Syntax

The default behavior of the StandardMetallicFuelRodOutputs action will create a standard set of outputs for both the fuel mesh block component of the rod and the cladding:

[StandardMetallicFuelRodOutputs<<<{"href": "index.html"}>>>]
  initial_pressure<<<{"description": "The initial pressure of the plenum gas in Pa."}>>> = 0.084e6 # Pa
  fuel_pellet_blocks<<<{"description": "The list of names of the blocks (subdomains) for the fuel pellets."}>>> = 'pellet'
  plenum_boundary_name<<<{"description": "The name of the boundary to use for the plenum temperature and plenum volume postprocessor calculations. Necessary only for non-standard clad only meshes."}>>> = 'inside_surfaces'
[]

It is also possible to use the StandardMetallicFuelRodOutputs action to generate outputs for only the fuel or only the cladding, depending on the rod_component selected. Here an a similar example of the StandardLWRFuelRodOutputs action for a fuel pellet only simulation is shown:

[StandardLWRFuelRodOutputs<<<{"href": "../StandardLWRFuelRodOutputs/index.html"}>>>]
  temperature<<<{"description": "The temperature variable name."}>>> = temp
  rod_component<<<{"description": "The mesh blocks for which the standard outputs should be generated. Options are FUEL CLAD BOTH."}>>> = fuel
  fuel_pellet_blocks<<<{"description": "The list of names of the blocks (subdomains) for the fuel pellets."}>>> = fuel_pellet
[]

note:Set Boundary Names for Cladding Burst Cases

For cladding burst tests which are run only on a hollow tube mesh, the advanced parameters plenum_boundary_name and external_clad_boundary_name should be set to ensure the correct boundary or boundaries are used in the calculation of the plenum volume, plenum temperature, the average clad interior temperature, and the radial displacement of the cladding exterior.

Note that when plenum_boundary_name is set, the boundaries specified by arguments for this parameter are used for both the plenum quantities and the average clad interior temperature.

[StandardLWRFuelRodOutputs<<<{"href": "../StandardLWRFuelRodOutputs/index.html"}>>>]
  rod_component<<<{"description": "The mesh blocks for which the standard outputs should be generated. Options are FUEL CLAD BOTH."}>>> = clad
  cladding_blocks<<<{"description": "The list of names of the block (subdomain) for the cladding. Necessary only for non-standard clad only meshes."}>>> = 0
  plenum_boundary_name<<<{"description": "The name of the boundary to use for the plenum temperature and plenum volume postprocessor calculations. Necessary only for non-standard clad only meshes."}>>> = 'left'
  external_clad_boundary_name<<<{"description": "The name of the boundary to use for the cladding radial displacement vector postprocessor calculations. Necessary only for non-standard clad only meshes."}>>> = 'right'
[]

Input Parameters

displacementsThe nonlinear displacement variables
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:The nonlinear displacement variables

Required Parameters

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
additional_volumesThe name of the postprocessor(s) that hold additional volumes that are connected to the cavity but not meshed.
C++ Type:std::vector<PostprocessorName>
Unit:(no unit assumed)
Controllable:No
Description:The name of the postprocessor(s) that hold additional volumes that are connected to the cavity but not meshed.
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
fuel_pellet_blocksThe list of names of the blocks (subdomains) for the fuel pellets.
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of names of the blocks (subdomains) for the fuel pellets.
fuel_pin_geometryName of fuel pin geometry UserObject.
C++ Type:UserObjectName
Controllable:No
Description:Name of fuel pin geometry UserObject.
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_pressureThe initial pressure of the plenum gas in Pa.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The initial pressure of the plenum gas in Pa.
internal_volume_additionAn additional volume to be included in the internal volume calculation.
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:An additional volume to be included in the internal volume calculation.
layeredFalseWhether this analysis uses a layered 1D or layered 2D mesh.
Default:False
C++ Type:bool
Controllable:No
Description:Whether this analysis uses a layered 1D or layered 2D mesh.
out_of_plane_strainstrain_yyThe out-of-plane strain variable.
Default:strain_yy
C++ Type:VariableName
Unit:(no unit assumed)
Controllable:No
Description:The out-of-plane strain variable.
out_of_plane_strain_claddingcladding_strain_yyName of the UserObject that contains the layered average out_of_plane_strain in the cladding.
Default:cladding_strain_yy
C++ Type:UserObjectName
Controllable:No
Description:Name of the UserObject that contains the layered average out_of_plane_strain in the cladding.
out_of_plane_strain_fuelfuel_strain_yyName of the UserObject that contains the layered average out_of_plane_strain in the fuel.
Default:fuel_strain_yy
C++ Type:UserObjectName
Controllable:No
Description:Name of the UserObject that contains the layered average out_of_plane_strain in the fuel.
rod_componentBOTHThe mesh blocks for which the standard outputs should be generated. Options are FUEL CLAD BOTH.
Default:BOTH
C++ Type:MooseEnum
Options:FUEL, CLAD, BOTH
Controllable:No
Description:The mesh blocks for which the standard outputs should be generated. Options are FUEL CLAD BOTH.
temperaturetemperatureThe temperature variable name.
Default:temperature
C++ Type:VariableName
Unit:(no unit assumed)
Controllable:No
Description:The temperature variable name.
temperature_of_additional_volumesThe name of the postprocessor(s) that hold the temperatures of the additional volumes.
C++ Type:std::vector<PostprocessorName>
Unit:(no unit assumed)
Controllable:No
Description:The name of the postprocessor(s) that hold the temperatures of the additional volumes.
use_automatic_differentiationFalseFlag to use automatic differentiation (AD) objects when possible
Default:False
C++ Type:bool
Controllable:No
Description:Flag to use automatic differentiation (AD) objects when possible

Optional Parameters

cladding_blocksThe list of names of the block (subdomain) for the cladding. Necessary only for non-standard clad only meshes.
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of names of the block (subdomain) for the cladding. Necessary only for non-standard clad only meshes.
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.
external_clad_boundary_nameThe name of the boundary to use for the cladding radial displacement vector postprocessor calculations. Necessary only for non-standard clad only meshes.
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The name of the boundary to use for the cladding radial displacement vector postprocessor calculations. Necessary only for non-standard clad only meshes.
plenum_boundary_nameThe name of the boundary to use for the plenum temperature and plenum volume postprocessor calculations. Necessary only for non-standard clad only meshes.
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The name of the boundary to use for the plenum temperature and plenum volume postprocessor calculations. Necessary only for non-standard clad only meshes.

Advanced Parameters

Associated Actions

Available Actions

Bison App
StandardMetallicFuelRodOutputsActionSets up the Auxvariables, Auxkernels, and Postprocessors required to generate the standard set of fuel rod simulation outputs for a 2D-RZ metallic fuel simulation. Not currently suitable for use with layered 1D or 3D meshes.

(test/tests/standard_metallic_outputs_action/x441_mini_fuel_rod.i)

[GlobalParams]
  density = 15800.0
  order = SECOND
  family = LAGRANGE
  energy_per_fission = 3.2e-11 # J/fission
  volumetric_locking_correction = false
  displacements = 'disp_x disp_y'
  temperature = temp
[]

[Problem]
  type = ReferenceResidualProblem
  reference_vector = 'ref'
  extra_tag_vectors = 'ref'
[]

[Mesh]
  coord_type = RZ
  [smeared_pellet_mesh]
    # Nominal Design Geometric Parameters (X441)
    type = FuelPinMeshGenerator
    clad_thickness = 0.38e-03
    pellet_outer_radius = 2.195e-03
    pellet_height = 3.4e-2
    clad_top_gap_height = 2.7e-2
    clad_gap_width = 0.345e-3
    bottom_clad_height = 2.24e-3
    top_clad_height = 2.24e-3
    clad_bot_gap_height = 0.31e-3 # arbitrary
    # meshing parameters
    clad_mesh_density = customize
    pellet_mesh_density = customize
    nx_p = 1
    ny_p = 5
    nx_c = 1
    ny_c = 5
    ny_cu = 1
    ny_cl = 1
    pellet_quantity = 1
    elem_type = QUAD8
  []
  # mesh options
  patch_size = 2
  patch_update_strategy = auto
  partitioner = centroid
  centroid_partitioner_direction = y
[]

[Variables]
  [disp_x]
  []
  [disp_y]
  []
  [temp]
    initial_condition = 298
  []
[]

[AuxVariables]
  # Aux variables for output
  [effective_creep_strain]
    order = CONSTANT
    family = MONOMIAL
  []
  [gap_cond]
    order = CONSTANT
    family = MONOMIAL
  []
  [coolant_htc]
    order = CONSTANT
    family = MONOMIAL
  []
  [cumulative_damage_index]
    order = CONSTANT
    family = MONOMIAL
  []
  [element_failed]
    order = CONSTANT
    family = MONOMIAL
  []
  [solid_swell]
    block = pellet
    order = CONSTANT
    family = MONOMIAL
  []
  [gas_swell]
    block = pellet
    order = CONSTANT
    family = MONOMIAL
  []
  [volumetric_strain]
    block = pellet
    order = CONSTANT
    family = MONOMIAL
  []
  [hoop_stress]
    order = CONSTANT
    family = MONOMIAL
  []
  [hoop_creep_strain]
    order = CONSTANT
    family = MONOMIAL
  []
  [hoop_elastic_strain]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [power_history]
    type = PiecewiseLinear
    x = '0 5e3'
    y = '0 44722'
  []
  [coolant_press_ramp]
    type = PiecewiseLinear
    x = '0 5e3'
    y = '0.151e6 0.151e6'
  []
  [coolant_temp_ramp]
    type = PiecewiseLinear
    x = '0 5e3'
    y = '298.0  648.0'
  []
  [axial_peaking_factors]
    type = PowerPeakingFunction
    fit = EBRII_ROW_4
    pellet_length = 343.0e-3
    pellet_y_start = 8.1e-3
  []
  [engr_radial_strain_fuel]
    type = ParsedFunction
    expression = 'fuel_disp_rad / 2.195e-03'
    symbol_values = 'max_fuel_radial_disp'
    symbol_names = 'fuel_disp_rad'
  []
  [engr_axial_strain_fuel]
    type = ParsedFunction
    expression = 'fuel_disp_axial / 343.0e-3'
    symbol_values = 'max_fuel_elongation'
    symbol_names = 'fuel_disp_axial'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [fuel]
    strain = FINITE
    generate_output = 'stress_xx stress_yy stress_zz vonmises_stress
      hydrostatic_stress creep_strain_xx creep_strain_yy creep_strain_zz
      elastic_strain_xx elastic_strain_yy elastic_strain_zz strain_xx strain_yy
      strain_zz'
    block = pellet
    eigenstrain_names = 'fuel_thermal_strain gas_swelling_eigenstrain
     solid_swelling_eigenstrain'
    extra_vector_tags = 'ref'
  []
  [clad]
    strain = FINITE
    generate_output = 'stress_xx stress_yy stress_zz vonmises_stress
      hydrostatic_stress creep_strain_xx creep_strain_yy creep_strain_zz
      elastic_strain_xx elastic_strain_yy elastic_strain_zz strain_xx strain_yy
      strain_zz'
    block = clad
    eigenstrain_names = 'clad_thermal_eigenstrain'
    extra_vector_tags = 'ref'
  []
[]

[Kernels]
  # Define kernels for the various terms in the PDE system
  [gravity]
    type = Gravity
    variable = disp_y
    value = -9.81
    extra_vector_tags = 'ref'
  []
  [heat]
    type = HeatConduction
    variable = temp
    extra_vector_tags = 'ref'
  []
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temp
    extra_vector_tags = 'ref'
  []
  [heat_source]
    type = NeutronHeatSource
    variable = temp
    block = pellet
    fission_rate = fission_rate
    extra_vector_tags = 'ref'
  []
[]

[AuxKernels]
  [effective_creep_strain]
    type = MaterialRealAux
    property = effective_creep_strain
    variable = effective_creep_strain
  []
  [conductance]
    type = MaterialRealAux
    property = gap_conductance
    variable = gap_cond
    boundary = 10
  []
  [failed_element]
    boundary = 2
    type = MaterialRealAux
    property = failed
    variable = element_failed
  []
  [gas_swell]
    type = MaterialRealAux
    variable = gas_swell
    property = gas_swelling
    execute_on = timestep_end
  []
  [solid_swell]
    type = MaterialRealAux
    variable = solid_swell
    property = solid_swelling
    execute_on = timestep_end
  []
  [volumetric_strain]
    type = RankTwoScalarAux
    rank_two_tensor = total_strain
    variable = volumetric_strain
    scalar_type = VolumetricStrain
    execute_on = timestep_end
    block = pellet
  []
  [hoop_stress]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = hoop_stress
    index_j = 2
    index_i = 2
    execute_on = timestep_end
  []
  [hoop_creep_strain]
    type = RankTwoAux
    rank_two_tensor = creep_strain
    variable = hoop_creep_strain
    index_j = 2
    index_i = 2
    execute_on = timestep_end
    block = clad
  []
  [hoop_elastic_strain]
    type = RankTwoAux
    rank_two_tensor = elastic_strain
    variable = hoop_elastic_strain
    index_j = 2
    index_i = 2
    execute_on = timestep_end
    block = clad
  []
[]

[Contact]
  [pellet_clad_mechanical]
    primary = 5
    secondary = 10
    penalty = 1e12
    model = frictionless
    formulation = kinematic
    normalize_penalty = true
    tangential_tolerance = 1e-3
    normal_smoothing_distance = 0.1
  []
[]

[ThermalContact]
  [thermal_contact]
    type = GapHeatTransfer
    variable = temp
    primary = 5
    secondary = 10
    quadrature = true
    gap_conductivity = 61.0
    min_gap = 0.345e-3
  []
[]

[BCs]
  [no_x_all]
    type = DirichletBC
    variable = disp_x
    boundary = 12
    value = 0.0
  []
  [no_y_fuel]
    type = DirichletBC
    variable = disp_y
    boundary = 20
    value = 0.0
  []
  [no_y_clad]
    type = DirichletBC
    variable = disp_y
    boundary = 1
    value = 0.0
  []
  [Pressure]
    [coolantPressure]
      boundary = '1 2 3'
      function = coolant_press_ramp
    []
  []
[]

[CoolantChannel]
  [convective_clad_surface]
    boundary = '1 2 3'
    variable = temp
    inlet_temperature = coolant_temp_ramp
    inlet_pressure = coolant_press_ramp
    inlet_massflux = 5261.5 # kg/m^2-sec
    coolant_material = sodium
    rod_diameter = 5.84e-3 # m
    rod_pitch = 7.48e-3 # m (Pitch-to-diameter Ratio = 1.28)
    linear_heat_rate = power_history
    axial_power_profile = axial_peaking_factors
    subchannel_geometry = triangular
  []
[]

[Materials]
  [fission_rate]
    type = UPuZrFissionRate
    rod_linear_power = power_history
    axial_power_profile = axial_peaking_factors
    pellet_radius = 2.195e-03
    X_Zr = 0.225
    block = pellet
    outputs = all
  []
  [burnup]
    type = UPuZrBurnup
    initial_X_Zr = 0.225
    density = 15800
    block = pellet
    outputs = all
  []
  [fuel_elasticity_tensor]
    type = UPuZrElasticityTensor
    X_Zr = 0.225
    X_Pu = 0.0
    block = pellet
  []
  [fuel_elastic_stress]
    type = ComputeMultipleInelasticStress
    tangent_operator = nonlinear
    inelastic_models = 'fuel_upuzrcreep'
    block = pellet
  []
  [fuel_upuzrcreep]
    type = UPuZrCreepUpdate
    block = pellet
    porosity = porosity
    max_inelastic_increment = 2e-3
  []
  [fuel_thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = pellet
    thermal_expansion_coeff = 1.18e-5
    stress_free_temperature = 295.0
    eigenstrain_name = fuel_thermal_strain
  []
  [gas_swelling]
    type = UPuZrGaseousEigenstrain
    eigenstrain_name = gas_swelling_eigenstrain
    initial_porosity = 0.0
    bubble_number_density = 1e20
    interconnection_initiating_porosity = 0.23
    interconnection_terminating_porosity = 0.25
    anisotropic_factor = 0.4
    outputs = all
    output_properties = 'porosity gaseous_porosity'
    block = pellet
  []
  [solid_swelling]
    type = BurnupDependentEigenstrain
    eigenstrain_name = solid_swelling_eigenstrain
    block = pellet
    swelling_name = 'solid_swelling'
  []
  [metal_fuel_thermal]
    type = UPuZrThermal
    block = pellet
    X_Zr = 0.225
    X_Pu = 0.0
    spheat_model = savage
    thcond_model = lanl
    porosity = porosity
  []
  [fuel_density]
    type = StrainAdjustedDensity
    block = pellet
    strain_free_density = 15800.0
  []
  [fission_gas_behavior]
    type = UPuZrFissionGasRelease
    block = pellet
    critical_porosity = 0.24
    fractional_fgr_initial = 0.8
    fractional_fgr_post = 1.0
    fission_rate = fission_rate
  []
  [clad_elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 1.88e11
    poissons_ratio = 0.236
    block = clad
  []
  [clad_stress]
    type = ComputeMultipleInelasticStress
    tangent_operator = nonlinear
    inelastic_models = 'clad_ht9creep'
    block = clad
  []
  [fast_flux]
    type = FastNeutronFlux
    block = clad
    factor = 2.47e19
  []
  [clad_ht9creep]
    type = HT9CreepUpdate
    block = clad
  []
  [thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = clad
    thermal_expansion_coeff = 1.2e-5
    stress_free_temperature = 295.0
    eigenstrain_name = clad_thermal_eigenstrain
  []
  [clad_thermal]
    type = HT9Thermal
    block = clad
  []
  [clad_density]
    type = StrainAdjustedDensity
    block = clad
    strain_free_density = 7874.0
  []
  [longHT9_failure]
    type = HT9FailureClad
    boundary = '1 2 3'
    method = cdf_long
    hoop_stress = stress_zz # Since 2D-RZ
  []
[]

[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 -ksp_gmres_restart'
  petsc_options_value = 'lu       superlu_dist                  51'
  line_search = 'none'

  l_max_its = 60
  l_tol = 8e-3
  nl_max_its = 40
  nl_rel_tol = 5e-4
  nl_abs_tol = 1e-7

  end_time = 5e3
  dtmin = 10
  dtmax = 5e5

  [Quadrature]
    order = fifth
    side_order = seventh
  []
  [TimeStepper]
    type = IterationAdaptiveDT
    timestep_limiting_postprocessor = creep_timestep
    dt = 1e3
    time_t = '0   1e4'
    time_dt = '1e2 1e2'
    iteration_window = 4
    optimal_iterations = 10
  []
[]

[Postprocessors]
  [approx_FCT]
    type = AverageNodalVariableValue
    boundary = 12
    variable = temp
  []
  [max_approx_FCT]
    type = TimeExtremeValue
    value_type = max
    postprocessor = approx_FCT
  []
  [ave_FST]
    type = SideAverageValue
    boundary = 10
    variable = temp
  []
  [max_ave_FST]
    type = TimeExtremeValue
    value_type = max
    postprocessor = ave_FST
  []
  [ave_CIT]
    type = SideAverageValue
    boundary = 5
    variable = temp
  []
  [max_ave_CIT]
    type = TimeExtremeValue
    value_type = max
    postprocessor = ave_CIT
  []
  [avg_clad_temp]
    type = ElementAverageValue
    variable = temp
    block = clad
  []
  [peak_clad_temp]
    type = ElementExtremeValue
    variable = temp
    value_type = max
    block = clad
  []
  [peak_fuel_temp]
    type = ElementExtremeValue
    variable = temp
    value_type = max
    block = pellet
  []
  [max_hydro]
    type = ElementExtremeValue
    variable = hydrostatic_stress
    value_type = max
    block = pellet
  []
  [min_hydro]
    type = ElementExtremeValue
    variable = hydrostatic_stress
    value_type = min
    block = pellet
  []
  [peak_porosity]
    type = ElementExtremeValue
    variable = porosity
    value_type = max
    block = pellet
  []
  [clad_inner_vol]
    type = InternalVolume
    boundary = 7
  []
  [clad_fuel_gap]
    type = NodalExtremeValue
    variable = penetration
    boundary = 10
  []
  [max_cont_press]
    type = NodalExtremeValue
    variable = contact_pressure
    boundary = 10
  []
  [flux_from_clad]
    type = SideDiffusiveFluxIntegral
    variable = temp
    boundary = 5
    diffusivity = thermal_conductivity
  []
  [flux_from_fuel]
    type = SideDiffusiveFluxIntegral
    variable = temp
    boundary = 10
    diffusivity = thermal_conductivity
  []
  [rod_total_power]
    type = ElementIntegralPower
    variable = temp
    use_material_fission_rate = true
    fission_rate_material = fission_rate
    block = pellet
  []
  [LHGR_W_per_cm]
    type = FunctionValuePostprocessor
    function = power_history
    scale_factor = 0.01
  []
  [max_cdf]
    type = ElementExtremeValue
    value_type = max
    variable = cumulative_damage_index
  []
  [creep_timestep]
    type = MaterialTimeStepPostprocessor
    block = pellet
  []
  [hydrostatic_stress]
    type = ElementAverageValue
    variable = hydrostatic_stress
    execute_on = 'initial timestep_end'
    block = pellet
  []
  [solid_swelling]
    type = ElementAverageValue
    variable = solid_swell
    block = pellet
  []
  [gas_swelling]
    type = ElementAverageValue
    variable = gas_swell
    block = pellet
  []
  [volumetric_strain]
    type = ElementAverageValue
    variable = volumetric_strain
    block = pellet
  []
  [fission_rate]
    type = ElementAverageValue
    variable = fission_rate
    block = pellet
  []
  [porosity]
    type = ElementAverageValue
    variable = porosity
    block = pellet
  []
  [gaseous_porosity]
    type = ElementAverageValue
    variable = gaseous_porosity
    block = pellet
  []
  [max_clad_hoop_creep]
    type = ElementExtremeValue
    value_type = max
    block = clad
    variable = hoop_creep_strain
  []
  [max_clad_creep_strain_mag]
    type = ElementExtremeValue
    value_type = max
    block = clad
    variable = effective_creep_strain
  []
  [max_fuel_radial_strain]
    type = ElementExtremeValue
    value_type = max
    block = pellet
    variable = strain_xx
  []
  [max_fuel_axial_strain]
    type = ElementExtremeValue
    value_type = max
    block = pellet
    variable = strain_yy
  []
  [max_fuel_elongation]
    type = NodalExtremeValue
    variable = disp_y
    boundary = 'all_pellet_exterior'
  []
  [max_fuel_radial_disp]
    type = NodalExtremeValue
    variable = disp_x
    boundary = 'all_pellet_exterior'
  []
  [engr_strain_fuel_radial]
    type = FunctionValuePostprocessor
    function = engr_radial_strain_fuel
  []
  [engr_strain_fuel_axial]
    type = FunctionValuePostprocessor
    function = engr_axial_strain_fuel
  []
  [max_clad_elongation]
    type = NodalExtremeValue
    variable = disp_y
    boundary = 'clad_outside_top clad_outside_right'
  []
[]

[VectorPostprocessors]
  [clad_total_hoop_strain]
    type = LineValueSampler
    variable = total_hoop_strain
    start_point = '2.90e-3 2.55e-3 0.0'
    end_point = '2.90e-3 0.05 0.0'
    num_points = 300
    sort_by = y
    outputs = 'vec1'
  []
[]

[StandardMetallicFuelRodOutputs]
  initial_pressure = 0.084e6 # Pa
  fuel_pellet_blocks = 'pellet'
  plenum_boundary_name = 'inside_surfaces'
[]

[PerformanceMetricOutputs]
[]

[Outputs]
  time_step_interval =  1
  color = true
  exodus = true
  perf_graph = true
  csv = true
  sync_times = '1e3 5e3'
  [console]
    type = Console
    max_rows = 25
    time_step_interval =  1
    output_linear = true
  []
  [chkfile]
    type = CSV
    file_base = x441_mini_fuel_rod_chkfile
    show = 'max_approx_FCT max_ave_FST max_ave_CIT average_burnup fission_gas_released_percentage max_clad_hoop_creep max_clad_creep_strain_mag max_fuel_elongation max_clad_elongation max_total_hoop_strain'
    execute_on = 'FINAL'
  []
  [vec1]
    type = CSV
    file_base = x441_mini_fuel_rod_vec1
    execute_on = 'FINAL'
  []
[]

[Debug]
  show_var_residual = 'disp_x disp_y temp'
  show_var_residual_norms = true
[]

(test/tests/standard_lwr_outputs_action/pellet_only.i)

initial_fuel_density = 10431.0

[GlobalParams]
  displacements = 'disp_x disp_y'
  order = SECOND
  density = ${initial_fuel_density}
  energy_per_fission = 3.20435313e-11
  temperature = temp
[]

[Mesh]
  coord_type = RZ
  patch_size = 10
  patch_update_strategy = auto
  [mesh]
    type = FileMeshGenerator
    file = pellet_only.e
  []
[]

[Variables]
  [temp]
    initial_condition = 293
  []
[]

[Functions]
  [power_profile]
    type = PiecewiseLinear
    x = '0 100 1e8'
    y = '0 20000 20000'
  []
  [axial_peaking_factors]
    type = ParsedFunction
    expression = 1
  []
  [q]
    type = CompositeFunction
    functions = 'power_profile axial_peaking_factors'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [pellets]
    block = fuel_pellet
    add_variables = true
    strain = FINITE
    eigenstrain_names = fuel_thermal_strain
    decomposition_method = EigenSolution
  []
[]

[Kernels]
  [heat]
    type = HeatConduction
    variable = temp
  []
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temp
  []
  [heat_source]
    type = NeutronHeatSource
    variable = temp
    burnup_function = burnup
  []
[]

[Burnup]
  [burnup]
    block = fuel_pellet
    rod_ave_lin_pow = power_profile
    axial_power_profile = axial_peaking_factors
    num_radial = 80
    num_axial = 20
    a_upper = 0.01496
    a_lower = 0.00226
    fuel_inner_radius = 0.0
    fuel_outer_radius = 0.005305
    fuel_volume_ratio = 1
    RPF = RPF
  []
[]

[BCs]
  [fuel_wall_temp]
    type = DirichletBC
    preset = false
    variable = temp
    boundary = '10'
    value = 673
  []
  [no_x_all]
    type = DirichletBC
    variable = disp_x
    boundary = 12
    value = 0.0
  []
  [no_y_fuel_bottom]
    type = DirichletBC
    variable = disp_y
    boundary = '1020'
    value = 0.0
  []
[]

[Materials]
  [fuel_thermal]
    type = UO2Thermal
    temperature = temp
    burnup_function = burnup
    thermal_conductivity_model = NFIR
    initial_porosity = 0.05
  []
  [fuel_elasticity_tensor]
    type = UO2ElasticityTensor
    block = fuel_pellet
  []
  [fuel_elastic_stress]
    type = ComputeFiniteStrainElasticStress
    block = fuel_pellet
  []
  [fuel_thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = fuel_pellet
    thermal_expansion_coeff = 10.0e-6
    temperature = temp
    stress_free_temperature = 293
    eigenstrain_name = fuel_thermal_strain
  []
  [fuel_density]
    type = StrainAdjustedDensity
    strain_free_density = ${initial_fuel_density}
  []
  [fission_gas_release]
    type = UO2Sifgrs
    temperature = temp
    burnup_function = burnup
  []
[]

[Executioner]
  type = Transient
  solve_type = 'PJFNK'
  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 = 15
  nl_abs_tol = 1e-10

  dtmax = 200
  dtmin = 200
  end_time = 200
[]

[StandardLWRFuelRodOutputs]
  temperature = temp
  rod_component = fuel
  fuel_pellet_blocks = fuel_pellet
[]

[Outputs]
  perf_graph = true
  exodus = false
  color = false
  csv = true
[]

(test/tests/standard_lwr_outputs_action/clad_nonstandard_names.i)

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  coord_type = RZ
  [mesh]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 2
    ny = 10
    xmin = 1
    xmax = 2
    ymax = 5
  []
[]

[Variables]
  [temperature]
    initial_condition = 300.
  []
[]

[Functions]
  [temperature_func]
    type = PiecewiseBilinear
    x = '0.0 0.025'
    y = '0.0 300.0 300.0'
    z = '1000.0 1400.0 1400.0 2000.0 1400.0 1400.0'
    axis = 1 # (0,1,2) => (x,y,z)
  []
  [power_history]
    type = PiecewiseLinear
    x = '0.000000 10800'
    y = '0.000000 16404.200000' #LHR5
    scale_factor = 1
  []
  [axial_peaking_factors]
    type = PiecewiseBilinear
    x = '0.00324 3.77797'
    y = '0.000000 10800'
    z = '1.0 1.0 1.0 1.0'
    axis = 1
    scale_factor = 1
  []
  [inner_pressure_func]
    type = PiecewiseLinear
    x = '0.     1000.  2000.  '
    y = '1.e+05 1.e+05 7.4e+05'
  []
  [pressure_ramp]
    type = PiecewiseLinear
    scale_factor = 1
    x = '0 10800.0'
    y = '0.00651 1.0'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [clad]
    add_variables = true
    strain = FINITE
  []
[]

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

[BCs]
  [no_x_all]
    type = DirichletBC
    variable = disp_x
    boundary = 'left'
    value = 0.0
  []
  [no_y_clad_bottom]
    type = DirichletBC
    variable = disp_y
    boundary = 'bottom'
    value = 0.0
  []
  [Pressure]
    [coolantPressure]
      boundary = 'right'
      factor = 15.5e6
      function = pressure_ramp
    []
  []
[]

[CoolantChannel]
  [convective_clad_surface]
    boundary = 'right'
    variable = temperature
    inlet_temperature = 580
    inlet_pressure    = 15.5e6
    inlet_massflux    = 3800
    rod_diameter      = 0.948e-2
    rod_pitch         = 1.26e-2
    linear_heat_rate  = power_history
    axial_power_profile = axial_peaking_factors
  []
[]

[Materials]
  [thermal]
    type = ZryThermal
    temperature = temperature
  []
  [clad_elasticity_tensor]
    type = ZryElasticityTensor
  []
  [clad_stress]
    type = ComputeFiniteStrainElasticStress
  []
  [density]
    type = StrainAdjustedDensity
    strain_free_density = 6550
  []
[]

[Executioner]
  type = Transient
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -sub_pc_type'
  petsc_options_value = 'asm       lu'
  line_search = 'none'
  l_max_its = 100
  l_tol = 1e-8
  nl_max_its = 15
  nl_rel_tol = 1.e-8
  nl_abs_tol = 1.e-10
  start_time = 0.
  end_time = 50.
  dtmax = 10.0
  dtmin = 10
[]

[StandardLWRFuelRodOutputs]
  rod_component = clad
  cladding_blocks = 0
  plenum_boundary_name = 'left'
  external_clad_boundary_name = 'right'
[]

[Outputs]
  perf_graph = true
  csv = true
[]

Description
Example Input File Syntax
Input Parameters
Associated Actions
Available Actions