CladdingOxideCompositeHeatConduction

Computes thermal conductivity.

Description

The kernel CladdingOxideCompositeHeatConduction computes an effective thermal conductivity based on the thermal conductivity of the underlying metal, the thermal conductivity of an oxide layer, and the oxide thickness. Note that the oxide layer is not explicitly modeled.

If a given element is wholly contained in either the metal or the oxide based on the location of the interface, it uses the thermal conductivity for that material. For elements where the metal/oxide interface is within the element, the effective thermal conductivity is computed as an anisotropic quantity. In the radial direction, the effective thermal conductivity is the series combination of the metal and oxide thermal conductivities based on the element length in the radial direction. Correspondingly, the effective thermal conductivity in the axial and azimuthal directions is the parallel combination of the metal and oxide thermal conductivities.

[Kernels<<<{"href": "../../syntax/Kernels/index.html"}>>>]
  [heat]
    type = CladdingOxideCompositeHeatConduction<<<{"description": "Computes thermal conductivity.", "href": "CladdingOxideCompositeHeatConduction.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    oxide_thickness<<<{"description": "coupled oxide layer thickness"}>>> = oxide_thickness
    clad_outer_radius<<<{"description": "clad outer radius (m)"}>>> = 0.01
    pilling_bedworth_ratio<<<{"description": "Pilling-Bedworth ratio. 1.56 for ZrO2"}>>> = 1.56
  []
[]

Input Parameters

clad_outer_radiusclad outer radius (m)
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:clad outer radius (m)
pilling_bedworth_ratioPilling-Bedworth ratio. 1.56 for ZrO2
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Pilling-Bedworth ratio. 1.56 for ZrO2
variableThe name of the variable that this residual object operates on
C++ Type:NonlinearVariableName
Unit:(no unit assumed)
Controllable:No
Description:The name of the variable that this residual object operates on

Required Parameters

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
displacementsThe displacements
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:The displacements
input_oxide_thickness0input oxide layer thickness
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:input oxide layer thickness
oxide_thicknesscoupled oxide layer thickness
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:coupled oxide layer thickness

Optional Parameters

absolute_value_vector_tagsThe tags for the vectors this residual object should fill with the absolute value of the residual contribution
C++ Type:std::vector<TagName>
Controllable:No
Description:The tags for the vectors this residual object should fill with the absolute value of the residual contribution
extra_matrix_tagsThe extra tags for the matrices this Kernel should fill
C++ Type:std::vector<TagName>
Controllable:No
Description:The extra tags for the matrices this Kernel should fill
extra_vector_tagsThe extra tags for the vectors this Kernel should fill
C++ Type:std::vector<TagName>
Controllable:No
Description:The extra tags for the vectors this Kernel should fill
matrix_tagssystemThe tag for the matrices this Kernel should fill
Default:system
C++ Type:MultiMooseEnum
Options:nontime, system
Controllable:No
Description:The tag for the matrices this Kernel should fill
vector_tagsnontimeThe tag for the vectors this Kernel should fill
Default:nontime
C++ Type:MultiMooseEnum
Options:nontime, time
Controllable:No
Description:The tag for the vectors this Kernel should fill

Contribution To Tagged Field Data 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.
diag_save_inThe name of auxiliary variables to save this Kernel's diagonal Jacobian contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
C++ Type:std::vector<AuxVariableName>
Unit:(no unit assumed)
Controllable:No
Description:The name of auxiliary variables to save this Kernel's diagonal Jacobian contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
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
save_inThe name of auxiliary variables to save this Kernel's residual contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
C++ Type:std::vector<AuxVariableName>
Unit:(no unit assumed)
Controllable:No
Description:The name of auxiliary variables to save this Kernel's residual contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
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_meshTrueWhether 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:True
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/thermalZry/thermal_test_zry.i)

(test/tests/thermalZry/thermal_test_zry.i)

#--------------------------------------------------------------------------------
#
#   - Geometry:
#      Plate 0.01x0.01
#
#
#
#   - B.C.                          ---------------
#                                   |      |      |
#                                   |      |      |
#      heat flux = 60,000 W/m^2     | mtl  | ox   |   T = 500 K
#                                   |      |      |
#                                   |      |      |
#                                   ---------------
#
#  The oxide thickness is 0 at time 0 and 1, 0.00156 at time 2, and 0.00324
#    at time 3.
#  After dividing by the Pilling-Bedworth ratio, 1.56, the metal/oxide
#    interface is at the far right of the domain at time 0 and 1, at the
#    boundary between the last two elements at time 2, and in the middle of
#    the second to last element at time 3.
#  k = 10 for the metal.  k = 2 for the oxide.
#  flux = deltaT*k/length
#  At time 1, deltaT = 60000*.01/10 = 60.  500 + 60 = 560, the temperature at
#    the left side of the domain.
#  At time 2, k_eff = .01*10*2/(10*.001+2*.009) = 7.142857...
#    60000*.01/k_eff = 84.  500 + 84 = 584.
#  At time 3, k_eff = .01*10*2/(10*.0015+2*.0085) = 6.25.
#    60000*.01/k_eff = 96.  500 + 96 = 596.
#  At time 3, the code computes a mixed thermal conductivity in the second to
#    last element.
#--------------------------------------------------------------------------------

[Mesh]
  [gen]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 10
    ny = 1
    xmax = 0.01
    ymax = 0.01
  []
[]

[Variables]
  [temperature]
    order = FIRST
    family = LAGRANGE
    initial_condition = 500
  []
[]

[AuxVariables]
  [oxide_thickness]
  []
[]


[Functions]
  [oxide_thickness]
    type = PiecewiseLinear
    x = '0 1 2       3'
    y = '0 0 0.00156 0.00234'
  []
[]

[Kernels]
  [heat]
    type = CladdingOxideCompositeHeatConduction
    variable = temperature
    oxide_thickness = oxide_thickness
    clad_outer_radius = 0.01
    pilling_bedworth_ratio = 1.56
  []
[]

[AuxKernels]
  [oxide_thickness]
    type = FunctionAux
    variable = oxide_thickness
    function = oxide_thickness
  []
[]

[BCs]
  [temperature_bc_1]
    type = NeumannBC
    variable = temperature
    boundary = 3
    value = 60000
  []
  [temperature_bc_3]
    type = DirichletBC
    boundary = 1
    value = 500
    variable = temperature
  []
[]

[Materials]
  [metal]
    type = GenericConstantMaterial
    prop_names = thermal_conductivity
    prop_values = 10
  []
  [oxide]
    type = GenericConstantMaterial
    prop_names = oxide_thermal_conductivity
    prop_values = 2
  []
[]

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options       = '-snes_ksp_ew'
  petsc_options_iname = '-ksp_gmres_restart'
  petsc_options_value = '101'

  line_search = 'none'

  l_max_its  = 100
  nl_max_its = 100
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-6
  l_tol      = 1e-5
  start_time = 0.0
  end_time   = 3
  dt         = 1.0
[]

[Postprocessors]
  [left_temperature]
    type = NodalVariableValue
    variable = temperature
    nodeid = 0
    execute_on = 'initial timestep_end'
  []
  [right_temperature]
    type = NodalVariableValue
    variable = temperature
    nodeid = 20
    execute_on = 'initial timestep_end'
  []
[]

[Outputs]
  [out]
    type = Exodus
  []
[]

(test/tests/thermalZry/thermal_test_zry.i)

#--------------------------------------------------------------------------------
#
#   - Geometry:
#      Plate 0.01x0.01
#
#
#
#   - B.C.                          ---------------
#                                   |      |      |
#                                   |      |      |
#      heat flux = 60,000 W/m^2     | mtl  | ox   |   T = 500 K
#                                   |      |      |
#                                   |      |      |
#                                   ---------------
#
#  The oxide thickness is 0 at time 0 and 1, 0.00156 at time 2, and 0.00324
#    at time 3.
#  After dividing by the Pilling-Bedworth ratio, 1.56, the metal/oxide
#    interface is at the far right of the domain at time 0 and 1, at the
#    boundary between the last two elements at time 2, and in the middle of
#    the second to last element at time 3.
#  k = 10 for the metal.  k = 2 for the oxide.
#  flux = deltaT*k/length
#  At time 1, deltaT = 60000*.01/10 = 60.  500 + 60 = 560, the temperature at
#    the left side of the domain.
#  At time 2, k_eff = .01*10*2/(10*.001+2*.009) = 7.142857...
#    60000*.01/k_eff = 84.  500 + 84 = 584.
#  At time 3, k_eff = .01*10*2/(10*.0015+2*.0085) = 6.25.
#    60000*.01/k_eff = 96.  500 + 96 = 596.
#  At time 3, the code computes a mixed thermal conductivity in the second to
#    last element.
#--------------------------------------------------------------------------------

[Mesh]
  [gen]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 10
    ny = 1
    xmax = 0.01
    ymax = 0.01
  []
[]

[Variables]
  [temperature]
    order = FIRST
    family = LAGRANGE
    initial_condition = 500
  []
[]

[AuxVariables]
  [oxide_thickness]
  []
[]


[Functions]
  [oxide_thickness]
    type = PiecewiseLinear
    x = '0 1 2       3'
    y = '0 0 0.00156 0.00234'
  []
[]

[Kernels]
  [heat]
    type = CladdingOxideCompositeHeatConduction
    variable = temperature
    oxide_thickness = oxide_thickness
    clad_outer_radius = 0.01
    pilling_bedworth_ratio = 1.56
  []
[]

[AuxKernels]
  [oxide_thickness]
    type = FunctionAux
    variable = oxide_thickness
    function = oxide_thickness
  []
[]

[BCs]
  [temperature_bc_1]
    type = NeumannBC
    variable = temperature
    boundary = 3
    value = 60000
  []
  [temperature_bc_3]
    type = DirichletBC
    boundary = 1
    value = 500
    variable = temperature
  []
[]

[Materials]
  [metal]
    type = GenericConstantMaterial
    prop_names = thermal_conductivity
    prop_values = 10
  []
  [oxide]
    type = GenericConstantMaterial
    prop_names = oxide_thermal_conductivity
    prop_values = 2
  []
[]

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options       = '-snes_ksp_ew'
  petsc_options_iname = '-ksp_gmres_restart'
  petsc_options_value = '101'

  line_search = 'none'

  l_max_its  = 100
  nl_max_its = 100
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-6
  l_tol      = 1e-5
  start_time = 0.0
  end_time   = 3
  dt         = 1.0
[]

[Postprocessors]
  [left_temperature]
    type = NodalVariableValue
    variable = temperature
    nodeid = 0
    execute_on = 'initial timestep_end'
  []
  [right_temperature]
    type = NodalVariableValue
    variable = temperature
    nodeid = 20
    execute_on = 'initial timestep_end'
  []
[]

[Outputs]
  [out]
    type = Exodus
  []
[]

Description
Input Parameters
Input Files