The Net Radiation Method System

Assumed geometry

The GrayDiffuseRadiation syntax assumes that radiation is exchanged in a cavity that is filled by a vacuum or a medium that is transparent to thermal radiation. The cavity must be enclosed by the boundaries specified in parameter "boundary". If the cavity is not enclosed, the result of the computation may be incorrect without MOOSE issuing a warning or error.

By default, MOOSE's raytracing module is used for calculating view factors. The method is discussed in detail in RayTracingViewFactor. The parameter "view_factor_calculator" can be set to "analytical" to force the use of UnobstructedPlanarViewFactor; at this point this is maintained as backward compatibility option and it may be removed in the future.

The "ray_tracing_face_order" parameter is important because it controls the accuracy of the view factor computation. The raytracing computation uses the QGRID quadrature available in libMesh which places quadrature points at uniform distances on the from and to boundary faces. Increasing the "ray_tracing_face_order" increases the accuracy of the view factors. Note, that mesh refinement also enhances the accuracy of view factors, but at the expense of more elements being crossed during the raytracing procedure. Uniformly increasing the number of quadrature points using "ray_tracing_face_order" is essentially equivalent to selectively refining the mesh on the faces in radiative transfer.

Available Actions

• Heat Transfer App
• RadiationTransferActionThis action sets up the net radiation calculation between specified sidesets.

Example Input syntax

[Problem<<<{"href": "../Problem/index.html"}>>>]
  kernel_coverage_check = false
[]

[Mesh<<<{"href": "../Mesh/index.html"}>>>]
  type = MeshGeneratorMesh

  [./cmg]
    type = CartesianMeshGenerator<<<{"description": "This CartesianMeshGenerator creates a non-uniform Cartesian mesh.", "href": "../../source/meshgenerators/CartesianMeshGenerator.html"}>>>
    dim<<<{"description": "The dimension of the mesh to be generated"}>>> = 2
    dx<<<{"description": "Intervals in the X direction"}>>> = '1 1.3 1.9'
    ix<<<{"description": "Number of grids in all intervals in the X direction (default to all one)"}>>> = '3 3 3'
    dy<<<{"description": "Intervals in the Y direction (required when dim>1 otherwise ignored)"}>>> = '2 1.2 0.9'
    iy<<<{"description": "Number of grids in all intervals in the Y direction (default to all one)"}>>> = '3 3 3'
    subdomain_id<<<{"description": "Block IDs (default to all zero)"}>>> = '0 1 0
                    4 5 2
                    0 3 0'
  [../]

  [./inner_bottom]
    type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = cmg
    primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = 1
    paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = 5
    new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = 'inner_bottom'
  [../]

  [./inner_left]
    type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = inner_bottom
    primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = 4
    paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = 5
    new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = 'inner_left'
  [../]

  [./inner_right]
    type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = inner_left
    primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = 2
    paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = 5
    new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = 'inner_right'
  [../]

  [./inner_top]
    type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = inner_right
    primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = 3
    paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = 5
    new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = 'inner_top'
  [../]

  [./rename]
    type = RenameBlockGenerator<<<{"description": "Changes the block IDs and/or block names for a given set of blocks defined by either block ID or block name. The changes are independent of ordering. The merging of blocks is supported.", "href": "../../source/meshgenerators/RenameBlockGenerator.html"}>>>
    old_block<<<{"description": "Elements with these block ID(s)/name(s) will be given the new block information specified in 'new_block'"}>>> = '1 2 3 4'
    new_block<<<{"description": "The new block ID(s)/name(s) to be given by the elements defined in 'old_block'."}>>> = '0 0 0 0'
    input<<<{"description": "The mesh we want to modify"}>>> = inner_top
  [../]
[]

[Variables<<<{"href": "../Variables/index.html"}>>>]
  [./temperature]
    block = 0
  [../]
[]

[Kernels<<<{"href": "../Kernels/index.html"}>>>]
  [./heat_conduction]
    type = HeatConduction<<<{"description": "Diffusive heat conduction term $-\\nabla\\cdot(k\\nabla T)$ of the thermal energy conservation equation", "href": "../../source/kernels/HeatConduction.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 0
    diffusion_coefficient<<<{"description": "Property name of the diffusion coefficient"}>>> = 5
  [../]
[]

[GrayDiffuseRadiation<<<{"href": "index.html"}>>>]
  [./cavity]
    boundary<<<{"description": "The boundaries that participate in the radiative exchange."}>>> = '4 5 6 7'
    emissivity<<<{"description": "Emissivities for each boundary."}>>> = '0.9 0.8 0.4 1'
    n_patches<<<{"description": "Number of radiation patches per sideset."}>>> = '2 2 2 3'
    partitioners<<<{"description": "Specifies a mesh partitioner to use when preparing the radiation patches."}>>> = 'centroid centroid centroid centroid'
    centroid_partitioner_directions<<<{"description": "Specifies the sort direction if using the centroid partitioner. Available options: x, y, z, radial"}>>> = 'x y y x'
    temperature<<<{"description": "The coupled temperature variable."}>>> = temperature
    adiabatic_boundary<<<{"description": "The adiabatic boundaries that participate in the radiative exchange."}>>> = '7'
    fixed_temperature_boundary<<<{"description": "The fixed temperature boundaries that participate in the radiative exchange."}>>> = '4'
    fixed_boundary_temperatures<<<{"description": "The temperatures of the fixed boundary."}>>> = '1200'
    view_factor_calculator<<<{"description": "The view factor calculator being used."}>>> = analytical
  [../]
[]

[BCs<<<{"href": "../BCs/index.html"}>>>]
  [./left]
    type = DirichletBC<<<{"description": "Imposes the essential boundary condition $u=g$, where $g$ is a constant, controllable value.", "href": "../../source/bcs/DirichletBC.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = left
    value<<<{"description": "Value of the BC"}>>> = 600
  [../]

  [./right]
    type = DirichletBC<<<{"description": "Imposes the essential boundary condition $u=g$, where $g$ is a constant, controllable value.", "href": "../../source/bcs/DirichletBC.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = right
    value<<<{"description": "Value of the BC"}>>> = 300
  [../]
[]

[Postprocessors<<<{"href": "../Postprocessors/index.html"}>>>]
  [./average_T_inner_right]
    type = SideAverageValue<<<{"description": "Computes the average value of a variable on a sideset. Note that this cannot be used on the centerline of an axisymmetric model.", "href": "../../source/postprocessors/SideAverageValue.html"}>>>
    variable<<<{"description": "The name of the variable which this postprocessor integrates"}>>> = temperature
    boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = inner_right
  [../]
[]

[Executioner<<<{"href": "../Executioner/index.html"}>>>]
  type = Steady
[]

[Outputs<<<{"href": "../Outputs/index.html"}>>>]
  exodus<<<{"description": "Output the results using the default settings for Exodus output."}>>> = true
[]

Input Parameters