HeatStructure2DCoupler

This component couples two 2D heat structures via a heat transfer coefficient.

Usage

This component has the following restrictions:

The coupled heat structures must be 2D heat structures.
The coupled heat structures must be of the same type.
Only one boundary name may be provided in each of the "primary_boundary" and "secondary_boundary" parameters.
The meshes along the coupled boundaries must be coincident, i.e., each node on each side must be at an identical location as a node on the other side.

Input Parameters

heat_transfer_coefficientHeat transfer coefficient function [W/(m^2-K)]
C++ Type:FunctionName
Controllable:No
Description:Heat transfer coefficient function [W/(m^2-K)]
primary_boundaryThe boundary of the first heat structure to couple
C++ Type:BoundaryName
Controllable:No
Description:The boundary of the first heat structure to couple
primary_heat_structureThe first heat structure to couple
C++ Type:std::string
Controllable:No
Description:The first heat structure to couple
secondary_boundaryThe boundary of the second heat structure to couple
C++ Type:BoundaryName
Controllable:No
Description:The boundary of the second heat structure to couple
secondary_heat_structureThe second heat structure to couple
C++ Type:std::string
Controllable:No
Description:The second heat structure to couple

Required 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:No
Description:Set the enabled status of the MooseObject.

Advanced Parameters

Formulation

The heat conduction equation is the following: $var element = document.getElementById("moose-equation-86be707d-bae5-4848-8a25-e8b9b01860bc");katex.render(" \\rho c_p \\pd{T}{t} - \\nabla \\cdot (k \\nabla T) = q''' \\eqc", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ where

$var element = document.getElementById("moose-equation-b4983533-06ad-4ddd-ac66-554caada555d");katex.render("\\rho", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is density,
$var element = document.getElementById("moose-equation-852a9dd1-6861-4656-bae5-1aada190caa7");katex.render("c_p", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is specific heat capacity,
$var element = document.getElementById("moose-equation-fbc7b08e-6545-4edd-88ec-a8ef547ab236");katex.render("k", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is thermal conductivity,
$var element = document.getElementById("moose-equation-a862980e-f972-460b-9f1f-e583546c3b7c");katex.render("T", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is temperature, and
$var element = document.getElementById("moose-equation-c5105d43-eb61-4386-8553-ea9e44714635");katex.render("q'''", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is a volumetric heat source.

Multiplying by a test function $var element = document.getElementById("moose-equation-3379b034-9f38-4815-a048-621d899ec8d4");katex.render("\\phi_i", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ and integrating by parts over the domain $var element = document.getElementById("moose-equation-5855e511-bcce-46c1-9e48-9259f68fb6e9");katex.render("\\Omega", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ gives $var element = document.getElementById("moose-equation-7e5a0f37-7230-49be-bba2-b4c7f4c58aac");katex.render(" \\pr{\\rho c_p \\pd{T}{t}, \\phi_i}_\\Omega + \\pr{k \\nabla T, \\nabla\\phi_i}_\\Omega - \\left\\langle k \\nabla T, \\phi_i\\mathbf{n}\\right\\rangle_{\\partial\\Omega} = \\pr{q''', \\phi_i}_\\Omega \\eqc", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ where $var element = document.getElementById("moose-equation-922cf81b-5682-4296-85c8-55bd05950c6b");katex.render("\\partial\\Omega", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the boundary of the domain $var element = document.getElementById("moose-equation-7ec7c9c2-9de9-4010-8464-aae10cf221c1");katex.render("\\Omega", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ .

For Neumann boundary conditions on the boundary $var element = document.getElementById("moose-equation-b0b411b4-ac46-4b9b-a4a4-fd2a455607cf");katex.render("\\Gamma", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , $var element = document.getElementById("moose-equation-2783a863-6b3b-4f01-86df-a93d5c5c2aeb");katex.render("k \\nabla T \\cdot \\mathbf{n}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is replaced with a known incoming heat flux function $var element = document.getElementById("moose-equation-cab89c44-c7e5-45a7-bc92-b125448dcb09");katex.render("q_b", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ :

var element = document.getElementById("moose-equation-bdff8977-ef97-4234-a59e-255c8beabd63");katex.render("k \\nabla T \\cdot \\mathbf{n} = q_b \\qquad \\mathbf{x} \\in \\Gamma \\eqp", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});

For the heat structure $var element = document.getElementById("moose-equation-bb1d4e81-75e7-401d-9860-f5bb742b61ed");katex.render("k", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , the incoming boundary heat flux $var element = document.getElementById("moose-equation-bdfbadc6-43d5-45d3-99d3-277d1dc6e11c");katex.render("q_b", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is computed as

$var element = document.getElementById("moose-equation-6e510d43-5d2f-416c-afba-b62161c79614");katex.render(" q_b = \\mathcal{H} (T_j - T_k) \\eqc", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ where

$var element = document.getElementById("moose-equation-30f39fa9-cf2c-463e-a912-2c922168cfb3");katex.render("\\mathcal{H}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the heat transfer coefficient,
$var element = document.getElementById("moose-equation-be67c4a3-248c-4891-89ac-8371a634108e");katex.render("T_k", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the surface temperature of the heat structure $var element = document.getElementById("moose-equation-63b12285-b896-4978-a8c1-9bccfc855bdb");katex.render("k", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , and
$var element = document.getElementById("moose-equation-9ce4ee32-13a0-48bb-a194-ccc3f4470e3d");katex.render("T_j", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the surface temperature of the coupled heat structure $var element = document.getElementById("moose-equation-1c778da1-090f-4709-bce0-b82a61ed8054");katex.render("j", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ .

Input Files

(modules/thermal_hydraulics/test/tests/components/heat_structure_2d_coupler/heat_structure_2d_coupler.i)

(modules/thermal_hydraulics/test/tests/components/heat_structure_2d_coupler/heat_structure_2d_coupler.i)

[HeatStructureMaterials]
  [hs_mat]
    type = SolidMaterialProperties
    k = 15
    cp = 500
    rho = 8000
  []
[]

[Components]
  [hs1]
    type = HeatStructureCylindrical
    position = '-0.5 0 0'
    orientation = '1 0 0'
    length = 0.5
    n_elems = 5

    names = 'region1'
    widths = '0.5'
    n_part_elems = '5'
    materials = 'hs_mat'

    initial_T = 500
  []

  [hs2]
    type = HeatStructureCylindrical
    position = '0 0 0'
    orientation = '1 0 0'
    length = '0.5 0.5'
    n_elems = '5 5'
    axial_region_names = 'axregion1 axregion2'

    names = 'region1 region2'
    widths = '0.5 0.2'
    n_part_elems = '5 3'
    materials = 'hs_mat hs_mat'

    initial_T = 300
  []

  [hs3]
    type = HeatStructureCylindrical
    position = '0.5 0 0'
    orientation = '1 0 0'
    length = 0.5
    n_elems = 5

    names = 'region1'
    widths = '0.5'
    n_part_elems = '5'
    materials = 'hs_mat'

    initial_T = 500
  []

  [hs_coupling_1_2]
    type = HeatStructure2DCoupler
    primary_heat_structure = hs2
    secondary_heat_structure = hs1
    primary_boundary = hs2:region1:start
    secondary_boundary = hs1:end
    heat_transfer_coefficient = 1000
  []

  [hs_coupling_2_3]
    type = HeatStructure2DCoupler
    primary_heat_structure = hs2
    secondary_heat_structure = hs3
    primary_boundary = hs2:axregion2:outer
    secondary_boundary = hs3:inner
    heat_transfer_coefficient = 500
  []
[]

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

[Postprocessors]
  [E_tot]
    type = ADHeatStructureEnergyRZ
    block = 'hs1:region1 hs2:region1 hs2:region2 hs3:region1'
    execute_on = 'INITIAL TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  scheme = 'bdf2'

  start_time = 0
  dt = 1000
  num_steps = 10
  abort_on_solve_fail = true

  solve_type = 'NEWTON'
  nl_rel_tol = 1e-6
  nl_abs_tol = 1e-6
  nl_max_its = 30

  l_tol = 1e-4
  l_max_its = 300
[]

[Outputs]
  file_base = 'cylindrical'
  exodus = true
[]

Usage
Input Parameters
Formulation
Input Files

Install MOOSE

New Users

Examples and Tutorials

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools

INL Applications and Remote Access