Layered Plenum Temperature Action System

The LayeredPlenumTemperature action simplifies the input file of simulations that estimate the temperature of the gas in the plenum by using a volume waiting average by creating the necessary, AuxVariables, AuxKernels, and Postprocessors required for the Layered Plenum Temperature Model.

Constructed Objects

The objects constructed by the LayeredPlenumTemperature action are summarized in Table 1. The Generated Name column refers to the name of the object constructed by the action. Depending on the functionality of the object the name could represent an aux variable or postprocessor.

Table 1: Correspondence Among Action Functionality and Moose/BISON Objects for the LayeredPlenumTemperature Action

Functionality	Generated Name	Created Classes	Associated Parameters
Penetration distance aux variable	`pt_distance`	Radius	`inner_radius_surfaces`: fuel surfaces along inner radius
			`point1`: Point one for defining an axis
			`point2`: Point two for defining an axis
		PlenumTemperatureDistance	`inner_surfaces`: The inner surfaces (e.g. fuel not on inner radius)
			`outer_surfaces1: The outer surfaces (e.g. cladding)
Gap value aux variable	`gap_value`	GapValueDefaultSameValue	`temp`: The temperature variable name
			`inner_surfaces`: The inner surfaces (e.g. fuel not on inner radius)
			`outer_surfaces`: The outer surfaces (e.g. cladding)
			`tangential_tolerance`: Tangential distance to extend edges of contact surfaces
			`warnings`: Whether to output warning messages concerning nodes not being found
Layered plenum temperature postprocessor	`user_defined_by_subblock`	LayeredPlenumTemperaturePostprocessor	`fuel_pin_geometry`: Name of the UserObject that reads the pin geometry from the mesh
			`boundary`: The boundary ID or name over which the plenum temperature should be measured
			`temp`: The temperature variable name
			`use_current_slice_heights`: A flag for using the current or initial slice heights
Cladding layered average out-of-plane strain	`out_of_plane_strain_cladding`	LayeredAverage	`out_of_plane_strain` : variable holding the out-of-plane strain
Fuel layered average out-of-plane strain	`out_of_plane_strain_fuel`	LayeredAverage	`out_of_plane_strain` : variable holding the out-of-plane strain
Layered plenum temperature userobject	`current_layer_heights`	LayeredPlenumTemperatureUserObject	`fuel_pin_geometry` : the userobject containing information read from the mesh

Example Input File Syntax

An example of using the LayeredPlenumTemperature action is given by:

[LayeredPlenumTemperature<<<{"href": "index.html"}>>>]
  [plenum_temp]
    boundary<<<{"description": "The boundary ID or name over which the plenum temperature should be measured."}>>> = 5
    fuel_pin_geometry<<<{"description": "Name of the UserObject that reads the pin geometry from the mesh."}>>> = pin_geometry
    out_of_plane_strain<<<{"description": "The out-of-plane strain nodal variable."}>>> = strain_yy_0
    inner_surfaces<<<{"description": "The inner surfaces (e.g. fuel not on inner radius)."}>>> = '5'
    outer_surfaces<<<{"description": "The outer surfaces (e.g. cladding)."}>>> = '10'
    temperature<<<{"description": "Coupled temperature (K)"}>>> = temp
  []
[]

Input Parameters

fuel_pin_geometryName of the UserObject that reads the pin geometry from the mesh.
C++ Type:UserObjectName
Controllable:No
Description:Name of the UserObject that reads the pin geometry from the mesh.
out_of_plane_strainThe out-of-plane strain nodal variable.
C++ Type:VariableName
Unit:(no unit assumed)
Controllable:No
Description:The out-of-plane strain nodal variable.
temperatureCoupled temperature (K)
C++ Type:NonlinearVariableName
Unit:(no unit assumed)
Controllable:No
Description:Coupled temperature (K)

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
axial_directionyThe axial direction of the model.
Default:y
C++ Type:MooseEnum
Options:x, y, z
Controllable:No
Description:The axial direction of the model.
boundaryThe boundary ID or name over which the plenum temperature should be measured.
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The boundary ID or name over which the plenum temperature should be measured.
clad_blocksThe list of names of the block (subdomain) for the cladding.
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of names of the block (subdomain) for the cladding.
execute_onINITIAL LINEARThe 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:INITIAL LINEAR
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.
familyLAGRANGESpecifies the family of FE shape functions to use for this variable.
Default:LAGRANGE
C++ Type:std::string
Controllable:No
Description:Specifies the family of FE shape functions to use for this variable.
fuel_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.
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.
inner_radius_surfacesThe fuel surfaces along the inner radius.
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The fuel surfaces along the inner radius.
inner_surfacesThe inner surfaces (e.g. fuel not on inner radius).
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The inner surfaces (e.g. fuel not on inner radius).
orderFIRSTSpecifies the order of the FE shape function to use for this variable.
Default:FIRST
C++ Type:MooseEnum
Options:CONSTANT, FIRST, SECOND, THIRD, FORTH
Controllable:No
Description:Specifies the order of the FE shape function to use for this variable.
outer_surfacesThe outer surfaces (e.g. cladding).
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The outer surfaces (e.g. cladding).
point1Point one for defining an axis.
C++ Type:libMesh::Point
Controllable:No
Description:Point one for defining an axis.
point2Point two for defining an axis.
C++ Type:libMesh::Point
Controllable:No
Description:Point two for defining an axis.
tangential_toleranceTangential distance to extend edges of contact surfaces.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Tangential distance to extend edges of contact surfaces.
use_current_slice_heightsFalseWhether or not to use the current slice height in the calculation of the plenum temperature. When set to true the results are more accurate, but the calculation is more difficult to converge.
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not to use the current slice height in the calculation of the plenum temperature. When set to true the results are more accurate, but the calculation is more difficult to converge.
warningsFalseWhether to output warning messages concerning nodes not being found.
Default:False
C++ Type:bool
Controllable:No
Description:Whether to output warning messages concerning nodes not being found.

Optional Parameters

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

Advanced Parameters

Associated Actions

Available Actions

Bison App
LayeredPlenumTemperatureActionCreates the necessary AuxVariables, AuxKernels, Postprocessors and UserObjects to calculate the temperature of plenum for Layered1D or Layered2D geometries.

(test/tests/layered_1D/layered_plenum_temperature.i)

#
# Mesh contains one fuel and clad slice of height 1 and a plenum
# slice of height 9.  The gap between the fuel and clad is set to 1.
# The fuel and clad surface temperature is set to 90.
# The plenum temperature is set to 50.
#
# The layered plenum temperature, which estimates the average temperature in the
# gap is calculated as:
#
# T_p = (50 * 9 * 2  + 90 * 1 * 1) / (9 * 2 + 1 * 1) = 52.10526
#

[GlobalParams]
  order = FIRST
  family = LAGRANGE
[]

[Mesh]
  coord_type = RZ
  [layered1D_mesh]
    type = Layered1DMeshGenerator
    uniform_slice_heights = true
    fuel_height = 1.0
    slices_per_block = 1
    include_fuel = true
    include_plenum = true
    plenum_height = 9.0
    pellet_bottom_coor = 0
    pellet_outer_radius = 1.0
    clad_gap_width = 1.0
    clad_thickness = 1.0
    elem_type = EDGE2
  []
[]

[Functions]
  [temp1]
    type = PiecewiseLinear
    x = '0 1'
    y = '0 90'
  []
  [temp2]
    type = PiecewiseLinear
    x = '0 1'
    y = '0 50'
  []
[]

[Variables]
  [temp]
    initial_condition = 0
  []
[]

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

[AuxVariables]
  [strain_yy_0]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[BCs]
  [temp1]
    type = FunctionDirichletBC
    boundary = '10002 10003'
    variable = temp
    function = temp1
  []
  [temp2]
    type = FunctionDirichletBC
    boundary = 10023
    variable = temp
    function = temp2
  []
[]

[UserObjects]
  [pin_geometry]
    type = Layered1DFuelPinGeometry
    mesh_generator = layered1D_mesh
  []
[]

[LayeredPlenumTemperature]
  [plenum_temp]
    boundary = 5
    fuel_pin_geometry = pin_geometry
    out_of_plane_strain = strain_yy_0
    inner_surfaces = '5'
    outer_surfaces = '10'
    temperature = temp
  []
[]

[Materials]
  [heat1]
    type = HeatConductionMaterial
    block = 'fuel clad'
    specific_heat = 1.0
    thermal_conductivity = 100000000.0
  []
  [density]
    type = ParsedMaterial
    block = 'fuel clad'
    property_name = density
    expression = 1.0
  []
[]

[Executioner]
  type = Transient

  solve_type = 'PJFNK'

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

  line_search = 'none'

  nl_abs_tol = 1e-10
  nl_rel_tol = 1e-12

  l_tol = 1e-8
  l_max_its = 100

  start_time = 0.0
  dt = 1.0
  end_time = 1.0
[]

[Outputs]
  csv = true
[]

Constructed Objects
Example Input File Syntax
Input Parameters
Associated Actions
Available Actions