ConstantViewFactorSurfaceRadiation

Description

ConstantViewFactorSurfaceRadiation implements the exchange of heat by radiation between sidesets by the net radiation method. The net radiation method is valid if the surfaces are gray, diffuse radiators. The purpose of the ConstantViewFactorSurfaceRadiation is to provide the means to model radiative exchange for coarse-grained model. ConstantViewFactorSurfaceRadiation does not compute radiative exchange between element surfaces, but it computes radiative transfer only in an average sense between sidesets.

The ConstantViewFactorSurfaceRadiation allows coupling radiative heat transfer to regions where the heat equation is solved. The net heat transfer caused by radiative heat transfer is coupled to the temperature field by GrayLambertNeumannBC. ConstantViewFactorSurfaceRadiation also supports the definition of adiabatic and fixed temperature sidesets. The temperature variable, i.e. the variable of the heat conduction equation, does not need to be defined on adiabatic and fixed temperature boundaries. This is particularly useful in cavities, where temperature is only defined on sidesets immediately adjacent to the heat conduction domain and not on the adiabatic and isothermal walls enclosing the cavity. There are three different types of boundary conditions in ConstantViewFactorSurfaceRadiation:

VARIABLE_TEMPERATURE are sidesets where temperature is provided by the temperature variable. The heat equation is coupled to the ConstantViewFactorSurfaceRadiation through these boundaries.
FIXED_TEMPERATURE are sidesets with temperature given as a function. The difference to VARIABLE_TEMPERATURE is that we do not need to solve for temperature on FIXED_TEMPERATURE. FIXED_TEMPERATURE sidesets can for example be the outside of a cavity that is in radiative heat transfer with the domain but is kept at constant temperature by being in contact with an effective coolant.
ADIABATIC sidesets have a zero net heat-flux. This implies that the inflow of radiation has to equal the outflow of radiation. The surface temperature assumes the value at which this condition is true. ADIABATIC sidesets do not require the temperature variable to be defined.

Explanation of the input parameters

This paragraph describes the input structure of the ConstantViewFactorSurfaceRadiation object. The following parameters are defined in detail:

temperature: the user must provide the name of the temperature variable as the temperature parameter. temperature is the variable that the heat conduction equation solves for.
boundary contains the names of _all_ sidesets that participate in the radiative heat exchange.
fixed_temperature_boundary should list all sidesets that are fixed temperature boundaries. fixed_temperature_boundary must be a subset of the boundary array.
fixed_boundary_temperatures contains function names specifying the temperatures on the fixed temperature boundaries. fixed_boundary_temperatures must be the same length as fixed_temperature_boundary.
adiabatic_boundary should list all sidesets that are adiabatic boundaries. adiabatic_boundary must be a subset of the boundary array.
view_factors lists the view factors $var element = document.getElementById("moose-equation-10b95dac-8509-4fea-a937-5ba92ac925b6");katex.render("F_{i,j}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ from sideset $var element = document.getElementById("moose-equation-ab8e1f84-09c3-4c9a-937d-9ab95372748f");katex.render("i", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ to sideset $var element = document.getElementById("moose-equation-5be1df7d-a11c-43d1-b5e8-e42328941cb1");katex.render("j", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ . The ordering of the view_factors follows the ordering in boundary. View factors must be added as square arrays. This may appear redundant because missing view factors can be computed using reciprocity, but this makes it impossible to check that each row of the view factor matrix sums to 1. This is absolutely necessary for conserving energy and ensuring stability of a coupled radiative transfer, heat conduction calculation. Row sums that are less than 5 percent off, are corrected to sum to 1. If the row sum is off by more than this threshold, an error is thrown. This is an example of the format of view_factors for three surfaces:
$var element = document.getElementById("moose-equation-98e64717-81a2-49c1-ad91-826f7bee3bc3");katex.render(" '~~F_{1,1} F_{1,2} F_{1,3} ; F_{2,1} F_{2,2} F_{2,3} ; F_{3,1} F_{3,2} F_{3,3} ~~'", element, {displayMode:true,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$

Example Input syntax

[Mesh]
  type = GeneratedMesh
  dim = 2
  xmin = 0
  xmax = 1
  ymin = 0
  ymax = 2
  nx = 1
  ny = 1
[]

[Problem]
  kernel_coverage_check = false
[]

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

[UserObjects]
  [./gray_lambert]
    type = ConstantViewFactorSurfaceRadiation
    boundary = 'bottom top left right'
    fixed_temperature_boundary = 'bottom top'
    fixed_boundary_temperatures = '550 300'
    adiabatic_boundary = 'right left'
    emissivity = '1 0.75 0.75 0.75'
    temperature = temperature
    view_factors = '0 0.123 0.6928 0.1841;
                    0.123 0 0.1841 0.6928;
                    0.2771 0.0736 0.4458 0.2035;
                    0.0736 0.2771 0.2035 0.4458'
  [../]
[]

[VectorPostprocessors]
  [./lambert_vpp]
    type = SurfaceRadiationVectorPostprocessor
    surface_radiation_object_name = gray_lambert
    information = 'temperature emissivity radiosity heat_flux_density'
  [../]

  [./view_factors]
    type = ViewfactorVectorPostprocessor
    surface_radiation_object_name = gray_lambert
  [../]
[]

[Postprocessors]
  [./heat_flux_density_bottom]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = HEAT_FLUX_DENSITY
    boundary = bottom
  [../]

  [./temperature_left]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = TEMPERATURE
    boundary = left
  [../]

  [./temperature_right]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = TEMPERATURE
    boundary = right
  [../]

  [./brightness_top]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = RADIOSITY
    boundary = top
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

/opt/civet/build_0/moose/modules/heat_conduction/test/tests/gray_lambert_radiator/gray_lambert_cavity.i

[Mesh]
  type = GeneratedMesh
  dim = 2
  xmin = 0
  xmax = 1
  ymin = 0
  ymax = 2
  nx = 1
  ny = 1
[]

[Problem]
  kernel_coverage_check = false
[]

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

[UserObjects]
  [./gray_lambert]
    type = ConstantViewFactorSurfaceRadiation
    boundary = 'bottom top left right'
    fixed_temperature_boundary = 'bottom top'
    fixed_boundary_temperatures = '550 300'
    adiabatic_boundary = 'right left'
    emissivity = '1 0.75 0.75 0.75'
    temperature = temperature
    view_factors = '0 0.123 0.6928 0.1841;
                    0.123 0 0.1841 0.6928;
                    0.2771 0.0736 0.4458 0.2035;
                    0.0736 0.2771 0.2035 0.4458'
  [../]
[]

[VectorPostprocessors]
  [./lambert_vpp]
    type = SurfaceRadiationVectorPostprocessor
    surface_radiation_object_name = gray_lambert
    information = 'temperature emissivity radiosity heat_flux_density'
  [../]

  [./view_factors]
    type = ViewfactorVectorPostprocessor
    surface_radiation_object_name = gray_lambert
  [../]
[]

[Postprocessors]
  [./heat_flux_density_bottom]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = HEAT_FLUX_DENSITY
    boundary = bottom
  [../]

  [./temperature_left]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = TEMPERATURE
    boundary = left
  [../]

  [./temperature_right]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = TEMPERATURE
    boundary = right
  [../]

  [./brightness_top]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = RADIOSITY
    boundary = top
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

Input Parameters

boundaryThe list of boundary IDs from the mesh where this boundary condition applies
C++ Type:std::vector
Options:
Description:The list of boundary IDs from the mesh where this boundary condition applies
emissivityEmissivities for each boundary.
C++ Type:std::vector
Options:
Description:Emissivities for each boundary.
temperatureThe coupled temperature variable.
C++ Type:std::vector
Options:
Description:The coupled temperature variable.
view_factorsThe view factors from sideset i to sideset j.
C++ Type:std::vector>
Options:
Description:The view factors from sideset i to sideset j.

Required Parameters

adiabatic_boundaryThe list of boundary IDs from the mesh that are adiabatic.
C++ Type:std::vector
Options:
Description:The list of boundary IDs from the mesh that are adiabatic.
execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:NONE INITIAL LINEAR NONLINEAR TIMESTEP_END TIMESTEP_BEGIN FINAL CUSTOM
Description:The list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM.
fixed_boundary_temperaturesThe temperatures of the fixed boundary.
C++ Type:std::vector
Options:
Description:The temperatures of the fixed boundary.
fixed_temperature_boundaryThe list of boundary IDs from the mesh with fixed temperatures.
C++ Type:std::vector
Options:
Description:The list of boundary IDs from the mesh with fixed temperatures.
stefan_boltzmann_constant5.67037e-08The Stefan-Boltzmann constant. Default value is in units of [W / m^2 K^4].
Default:5.67037e-08
C++ Type:double
Options:
Description:The Stefan-Boltzmann constant. Default value is in units of [W / m^2 K^4].

Optional Parameters

allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Options:
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector
Options:
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Options:
Description:Set the enabled status of the MooseObject.
use_displaced_meshFalseWhether 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:False
C++ Type:bool
Options:
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

Input Files

modules/heat_conduction/test/tests/gray_lambert_radiator/coupled_heat_conduction.i
modules/heat_conduction/test/tests/gray_lambert_radiator/gray_lambert_cavity.i

modules/heat_conduction/test/tests/gray_lambert_radiator/coupled_heat_conduction.i

[Problem]
  kernel_coverage_check = false
[]

[Mesh]
  type = MeshGeneratorMesh

  [./cartesian]
    type = CartesianMeshGenerator
    dim = 2
    dx = '1 1 1'
    ix = '2 2 2'
    dy = '5'
    iy = '10'
    subdomain_id = '1 2 3'
  [../]

  [./break_sides]
    type = BreakBoundaryOnSubdomainGenerator
    boundaries = 'bottom top'
    input = cartesian
  [../]

  [./left_interior]
    type = SideSetsBetweenSubdomainsGenerator
    master_block = 1
    paired_block = 2
    new_boundary = left_interior
    input = break_sides
  [../]

  [./right_interior]
    type = SideSetsBetweenSubdomainsGenerator
    master_block = 3
    paired_block = 2
    new_boundary = right_interior
    input = left_interior
  [../]
  [./rename]
    type = RenameBlockGenerator
    input = right_interior
    old_block_id = '1 2 3'
    new_block_id = '1 4 3'
  [../]
[]


[Variables]
  [./temperature]
    initial_condition = 300
    block = '1 3'
  [../]
[]

[Kernels]
  [./heat_conduction]
    type = HeatConduction
    variable = temperature
    diffusion_coefficient = 1
    block = '1 3'
  [../]
[]

[UserObjects]
  [./cavity_radiation]
    type = ConstantViewFactorSurfaceRadiation
    boundary = 'left_interior right_interior bottom_to_2 top_to_2'
    temperature = temperature
    emissivity = '0.8 0.8 0.8 0.8'
    adiabatic_boundary = 'bottom_to_2 top_to_2'
    # these view factors are made up to exactly balance energy
    # transfer through the cavity
    view_factors = '0    0.8 0.1 0.1;
                    0.8  0   0.1 0.1;
                    0.45 0.45  0 0.1;
                    0.45 0.45 0.1  0'
    execute_on = 'INITIAL LINEAR TIMESTEP_END'
  [../]
[]

[BCs]
  [./bottom_left]
    type = DirichletBC
    preset = false
    variable = temperature
    boundary = bottom_to_1
    value = 1500
  [../]

  [./top_right]
    type = DirichletBC
    preset = false
    variable = temperature
    boundary = top_to_3
    value = 300
  [../]

  [./radiation]
    type = GrayLambertNeumannBC
    variable = temperature
    reconstruct_emission = false
    surface_radiation_object_name = cavity_radiation
    boundary = 'left_interior right_interior'
  [../]
[]

[Postprocessors]
  [./qdot_left]
    type = GrayLambertSurfaceRadiationPP
    boundary = left_interior
    surface_radiation_object_name = cavity_radiation
    return_type = HEAT_FLUX_DENSITY
  [../]

  [./qdot_right]
    type = GrayLambertSurfaceRadiationPP
    boundary = right_interior
    surface_radiation_object_name = cavity_radiation
    return_type = HEAT_FLUX_DENSITY
  [../]

  [./qdot_top]
    type = GrayLambertSurfaceRadiationPP
    boundary = top_to_2
    surface_radiation_object_name = cavity_radiation
    return_type = HEAT_FLUX_DENSITY
  [../]

  [./qdot_bottom]
    type = GrayLambertSurfaceRadiationPP
    boundary = bottom_to_2
    surface_radiation_object_name = cavity_radiation
    return_type = HEAT_FLUX_DENSITY
  [../]
[]

[Executioner]
  type = Steady
[]

[Outputs]
  exodus = true
[]

modules/heat_conduction/test/tests/gray_lambert_radiator/gray_lambert_cavity.i

[Mesh]
  type = GeneratedMesh
  dim = 2
  xmin = 0
  xmax = 1
  ymin = 0
  ymax = 2
  nx = 1
  ny = 1
[]

[Problem]
  kernel_coverage_check = false
[]

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

[UserObjects]
  [./gray_lambert]
    type = ConstantViewFactorSurfaceRadiation
    boundary = 'bottom top left right'
    fixed_temperature_boundary = 'bottom top'
    fixed_boundary_temperatures = '550 300'
    adiabatic_boundary = 'right left'
    emissivity = '1 0.75 0.75 0.75'
    temperature = temperature
    view_factors = '0 0.123 0.6928 0.1841;
                    0.123 0 0.1841 0.6928;
                    0.2771 0.0736 0.4458 0.2035;
                    0.0736 0.2771 0.2035 0.4458'
  [../]
[]

[VectorPostprocessors]
  [./lambert_vpp]
    type = SurfaceRadiationVectorPostprocessor
    surface_radiation_object_name = gray_lambert
    information = 'temperature emissivity radiosity heat_flux_density'
  [../]

  [./view_factors]
    type = ViewfactorVectorPostprocessor
    surface_radiation_object_name = gray_lambert
  [../]
[]

[Postprocessors]
  [./heat_flux_density_bottom]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = HEAT_FLUX_DENSITY
    boundary = bottom
  [../]

  [./temperature_left]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = TEMPERATURE
    boundary = left
  [../]

  [./temperature_right]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = TEMPERATURE
    boundary = right
  [../]

  [./brightness_top]
    type = GrayLambertSurfaceRadiationPP
    surface_radiation_object_name = gray_lambert
    return_type = RADIOSITY
    boundary = top
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

Description
Explanation of the input parameters
Example Input syntax
Input Parameters
Input Files

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools