ViewFactorObjectSurfaceRadiation computes radiative heat transfer between side sets and the view factors are computed by a ViewFactor object. More...

#include <ViewFactorObjectSurfaceRadiation.h>

Inheritance diagram for ViewFactorObjectSurfaceRadiation:

[legend]

Public Types
enum	RAD_BND_TYPE { VARIABLE_TEMPERATURE = 0, FIXED_TEMPERATURE = 4, ADIABATIC = 8 }
	Define enum for boundary type. More...

Public Member Functions
	ViewFactorObjectSurfaceRadiation (const InputParameters &parameters)

virtual void	execute () override

virtual void	initialize () override

virtual void	finalize () override


Real	getSurfaceIrradiation (BoundaryID id) const
	public interface of this UserObject More...

Real	getSurfaceHeatFluxDensity (BoundaryID id) const

Real	getSurfaceTemperature (BoundaryID id) const

Real	getSurfaceRadiosity (BoundaryID id) const

Real	getSurfaceEmissivity (BoundaryID id) const

Real	getViewFactor (BoundaryID from_id, BoundaryID to_id) const

std::set< BoundaryID >	getSurfaceIDs () const

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual std::vector< std::vector< Real > >	setViewFactors () override
	a purely virtual function that defines where view factors come from More...

virtual void	threadJoin (const UserObject &y) override

Protected Attributes
const Real	_sigma_stefan_boltzmann
	Stefan-Boltzmann constant. More...

unsigned int	_n_sides
	number of active boundary ids More...

const VariableValue &	_temperature
	the coupled temperature variable More...

const std::vector< Real >	_emissivity
	constant emissivity for each boundary More...

std::vector< const Function * >	_fixed_side_temperature
	side id to index map, side ids can have holes or be out of order More...

std::vector< Real >	_radiosity
	the radiosity of each surface More...

std::vector< Real >	_heat_flux_density
	the heat flux density qdot More...

std::vector< Real >	_side_temperature
	the average temperature: this could be important for adiabatic walls More...

std::vector< enum RAD_BND_TYPE >	_side_type
	the type of the side, allows lookup index -> type More...

std::map< BoundaryID, unsigned int >	_side_id_index
	side id to index map, side ids can have holes or be out of order More...

std::vector< Real >	_areas
	the area by participating side set More...

std::vector< Real >	_beta
	the average value of sigma * eps * T^4 More...

std::vector< Real >	_surface_irradiation
	the irradiation into each surface More...

std::map< unsigned int, unsigned int >	_fixed_side_id_index
	side id to index map for isothermal boundaries, side ids can have holes or be out of order More...

std::set< unsigned int >	_adiabatic_side_ids
	the set of adiabatic boundaries More...

std::vector< std::vector< Real > >	_view_factors
	the view factors which are set by setViewFactors by derived classes More...

Detailed Description

ViewFactorObjectSurfaceRadiation computes radiative heat transfer between side sets and the view factors are computed by a ViewFactor object.

Definition at line 25 of file ViewFactorObjectSurfaceRadiation.h.

Member Enumeration Documentation

◆ RAD_BND_TYPE

enum GrayLambertSurfaceRadiationBase::RAD_BND_TYPE

inherited

Define enum for boundary type.

Enumerator
VARIABLE_TEMPERATURE
FIXED_TEMPERATURE
ADIABATIC

Definition at line 37 of file GrayLambertSurfaceRadiationBase.h.

   {
     VARIABLE_TEMPERATURE = 0,
     FIXED_TEMPERATURE = 4,
     ADIABATIC = 8
   };

Constructor & Destructor Documentation

◆ ViewFactorObjectSurfaceRadiation()

ViewFactorObjectSurfaceRadiation::ViewFactorObjectSurfaceRadiation ( const InputParameters & parameters )

Definition at line 29 of file ViewFactorObjectSurfaceRadiation.C.

   : GrayLambertSurfaceRadiationBase(parameters)
 {
 }

Member Function Documentation

◆ execute()

void GrayLambertSurfaceRadiationBase::execute ( )

overridevirtualinherited

Definition at line 142 of file GrayLambertSurfaceRadiationBase.C.

 {
   mooseAssert(_side_id_index.find(_current_boundary_id) != _side_id_index.end(),
               "Current boundary id not in _side_id_index.");
   unsigned int index = _side_id_index.find(_current_boundary_id)->second;
  
   for (unsigned int qp = 0; qp < _qrule->n_points(); qp++)
   {
     _areas[index] += _JxW[qp] * _coord[qp];
  
     Real temp = 0;
     if (_side_type[index] == ADIABATIC)
       continue;
     else if (_side_type[index] == VARIABLE_TEMPERATURE)
       temp = _temperature[qp];
     else if (_side_type[index] == FIXED_TEMPERATURE)
     {
       unsigned int iso_index = _fixed_side_id_index.find(_current_boundary_id)->second;
       temp = _fixed_side_temperature[iso_index]->value(_t, _q_point[qp]);
     }
  
     _beta[index] += _JxW[qp] * _coord[qp] * _sigma_stefan_boltzmann * _emissivity[index] *
                     MathUtils::pow(temp, 4);
     _side_temperature[index] += _JxW[qp] * _coord[qp] * temp;
   }
 }

◆ finalize()

void GrayLambertSurfaceRadiationBase::finalize ( )

overridevirtualinherited

Definition at line 186 of file GrayLambertSurfaceRadiationBase.C.

 {
   // need to do some parallel communiction here
   gatherSum(_areas);
   gatherSum(_beta);
   gatherSum(_side_temperature);
  
   // first compute averages from the totals
   for (unsigned int j = 0; j < _n_sides; ++j)
   {
     _beta[j] /= _areas[j];
     _side_temperature[j] /= _areas[j];
   }
  
   // matrix and rhs vector for the view factor calculation
   DenseMatrix<Real> matrix(_n_sides, _n_sides);
   DenseVector<Real> rhs(_n_sides);
   DenseVector<Real> radiosity(_n_sides);
   for (unsigned int i = 0; i < _n_sides; ++i)
   {
     rhs(i) = _beta[i];
     matrix(i, i) = 1;
     for (unsigned int j = 0; j < _n_sides; ++j)
     {
       if (_side_type[i] == ADIABATIC)
         matrix(i, j) -= _view_factors[i][j];
       else
         matrix(i, j) -= (1 - _emissivity[i]) * _view_factors[i][j];
     }
   }
  
   // compute the radiosityes
   matrix.lu_solve(rhs, radiosity);
  
   // store the radiosity, temperatures and heat flux density for each surface
   for (unsigned int i = 0; i < _n_sides; ++i)
   {
     _radiosity[i] = radiosity(i);
  
     // _heat_flux_density is obtained from a somewhat cumbersome relation
     // but it has the advantage that we do not divide by 1 - emissivity
     // which blows up for black bodies
     _heat_flux_density[i] = radiosity(i);
     for (unsigned int j = 0; j < _n_sides; ++j)
       _heat_flux_density[i] -= _view_factors[i][j] * radiosity(j);
  
     if (_side_type[i] == ADIABATIC)
       _side_temperature[i] =
           std::pow((radiosity(i) + (1 - _emissivity[i]) / _emissivity[i] * _heat_flux_density[i]) /
                        _sigma_stefan_boltzmann,
                    0.25);
  
     // compute the surface irradiation into i from the radiosities
     _surface_irradiation[i] = 0;
     for (unsigned int j = 0; j < _n_sides; ++j)
       _surface_irradiation[i] += _view_factors[i][j] * radiosity(j);
   }
 }

◆ getSurfaceEmissivity()

Real GrayLambertSurfaceRadiationBase::getSurfaceEmissivity ( BoundaryID id ) const

inherited

Definition at line 301 of file GrayLambertSurfaceRadiationBase.C.

 {
   if (_side_id_index.find(id) == _side_id_index.end())
     return 1;
   return _emissivity[_side_id_index.find(id)->second];
 }

Referenced by GrayLambertNeumannBC::computeQpJacobian(), GrayLambertNeumannBC::computeQpResidual(), and SurfaceRadiationVectorPostprocessor::execute().

◆ getSurfaceHeatFluxDensity()

Real GrayLambertSurfaceRadiationBase::getSurfaceHeatFluxDensity ( BoundaryID id ) const

inherited

Definition at line 277 of file GrayLambertSurfaceRadiationBase.C.

 {
   if (_side_id_index.find(id) == _side_id_index.end())
     return 0;
   return _heat_flux_density[_side_id_index.find(id)->second];
 }

Referenced by GrayLambertNeumannBC::computeQpResidual(), SurfaceRadiationVectorPostprocessor::execute(), and GrayLambertSurfaceRadiationPP::getValue().

◆ getSurfaceIDs()

std::set< BoundaryID > GrayLambertSurfaceRadiationBase::getSurfaceIDs ( ) const

inherited

Definition at line 260 of file GrayLambertSurfaceRadiationBase.C.

 {
   std::set<BoundaryID> surface_ids;
   for (auto & p : _side_id_index)
     surface_ids.insert(p.first);
   return surface_ids;
 }

Referenced by SurfaceRadiationVectorPostprocessor::initialize(), and ViewfactorVectorPostprocessor::initialize().

◆ getSurfaceIrradiation()

Real GrayLambertSurfaceRadiationBase::getSurfaceIrradiation ( BoundaryID id ) const

inherited

public interface of this UserObject

Definition at line 269 of file GrayLambertSurfaceRadiationBase.C.

 {
   if (_side_id_index.find(id) == _side_id_index.end())
     return 0;
   return _surface_irradiation[_side_id_index.find(id)->second];
 }

Referenced by GrayLambertNeumannBC::computeQpResidual().

◆ getSurfaceRadiosity()

Real GrayLambertSurfaceRadiationBase::getSurfaceRadiosity ( BoundaryID id ) const

inherited

Definition at line 293 of file GrayLambertSurfaceRadiationBase.C.

 {
   if (_side_id_index.find(id) == _side_id_index.end())
     return 0;
   return _radiosity[_side_id_index.find(id)->second];
 }

Referenced by SurfaceRadiationVectorPostprocessor::execute(), and GrayLambertSurfaceRadiationPP::getValue().

◆ getSurfaceTemperature()

Real GrayLambertSurfaceRadiationBase::getSurfaceTemperature ( BoundaryID id ) const

inherited

Definition at line 285 of file GrayLambertSurfaceRadiationBase.C.

 {
   if (_side_id_index.find(id) == _side_id_index.end())
     return 0;
   return _side_temperature[_side_id_index.find(id)->second];
 }

Referenced by SurfaceRadiationVectorPostprocessor::execute(), and GrayLambertSurfaceRadiationPP::getValue().

◆ getViewFactor()

Real GrayLambertSurfaceRadiationBase::getViewFactor	(	BoundaryID	from_id,
		BoundaryID	to_id
	)		const

inherited

Definition at line 309 of file GrayLambertSurfaceRadiationBase.C.

 {
   if (_side_id_index.find(from_id) == _side_id_index.end())
     return 0;
   if (_side_id_index.find(to_id) == _side_id_index.end())
     return 0;
   return _view_factors[_side_id_index.find(from_id)->second][_side_id_index.find(to_id)->second];
 }

Referenced by ViewfactorVectorPostprocessor::execute().

◆ initialize()

void GrayLambertSurfaceRadiationBase::initialize ( )

overridevirtualinherited

Reimplemented in ConstantViewFactorSurfaceRadiation.

Definition at line 170 of file GrayLambertSurfaceRadiationBase.C.

 {
   // view factors are obtained here to make sure that another object had
   // time to compute them on exec initial
   _view_factors = setViewFactors();
  
   // initialize areas, beta, side temps
   for (unsigned int j = 0; j < _n_sides; ++j)
   {
     _areas[j] = 0;
     _beta[j] = 0;
     _side_temperature[j] = 0;
   }
 }

Referenced by ConstantViewFactorSurfaceRadiation::initialize().

◆ setViewFactors()

std::vector< std::vector< Real > > ViewFactorObjectSurfaceRadiation::setViewFactors ( )

overrideprotectedvirtual

a purely virtual function that defines where view factors come from

Implements GrayLambertSurfaceRadiationBase.

Definition at line 36 of file ViewFactorObjectSurfaceRadiation.C.

 {
   const ViewFactorBase & view_factor_uo = getUserObject<ViewFactorBase>("view_factor_object_name");
   std::vector<BoundaryName> boundary_names = getParam<std::vector<BoundaryName>>("boundary");
   std::vector<std::vector<Real>> vf(_n_sides);
  
   for (unsigned int i = 0; i < _n_sides; ++i)
   {
     vf[i].resize(_n_sides);
  
     for (unsigned int j = 0; j < _n_sides; ++j)
       vf[i][j] = view_factor_uo.getViewFactor(boundary_names[i], boundary_names[j]);
   }
   return vf;
 }

◆ threadJoin()

void GrayLambertSurfaceRadiationBase::threadJoin ( const UserObject & y )

overrideprotectedvirtualinherited

Definition at line 246 of file GrayLambertSurfaceRadiationBase.C.

 {
   const GrayLambertSurfaceRadiationBase & pps =
       static_cast<const GrayLambertSurfaceRadiationBase &>(y);
  
   for (unsigned int j = 0; j < _n_sides; ++j)
   {
     _areas[j] += pps._areas[j];
     _side_temperature[j] += pps._side_temperature[j];
     _beta[j] += pps._beta[j];
   }
 }

◆ validParams()

InputParameters ViewFactorObjectSurfaceRadiation::validParams ( )

static

Definition at line 18 of file ViewFactorObjectSurfaceRadiation.C.

 {
   InputParameters params = GrayLambertSurfaceRadiationBase::validParams();
   params.addRequiredParam<UserObjectName>("view_factor_object_name",
                                           "Name of the ViewFactor userobjects.");
   params.addClassDescription(
       "ViewFactorObjectSurfaceRadiation computes radiative heat transfer between side sets and the "
       "view factors are computed by a ViewFactor object");
   return params;
 }

Member Data Documentation

◆ _adiabatic_side_ids

std::set<unsigned int> GrayLambertSurfaceRadiationBase::_adiabatic_side_ids

protectedinherited

the set of adiabatic boundaries

Definition at line 103 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::GrayLambertSurfaceRadiationBase().

◆ _areas

std::vector<Real> GrayLambertSurfaceRadiationBase::_areas

protectedinherited

the area by participating side set

Definition at line 91 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::execute(), GrayLambertSurfaceRadiationBase::finalize(), GrayLambertSurfaceRadiationBase::initialize(), and GrayLambertSurfaceRadiationBase::threadJoin().

◆ _beta

std::vector<Real> GrayLambertSurfaceRadiationBase::_beta

protectedinherited

the average value of sigma * eps * T^4

Definition at line 94 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::execute(), GrayLambertSurfaceRadiationBase::finalize(), GrayLambertSurfaceRadiationBase::initialize(), and GrayLambertSurfaceRadiationBase::threadJoin().

◆ _emissivity

const std::vector<Real> GrayLambertSurfaceRadiationBase::_emissivity

protectedinherited

constant emissivity for each boundary

Definition at line 70 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::execute(), GrayLambertSurfaceRadiationBase::finalize(), GrayLambertSurfaceRadiationBase::getSurfaceEmissivity(), and GrayLambertSurfaceRadiationBase::GrayLambertSurfaceRadiationBase().

◆ _fixed_side_id_index

std::map<unsigned int, unsigned int> GrayLambertSurfaceRadiationBase::_fixed_side_id_index

protectedinherited

side id to index map for isothermal boundaries, side ids can have holes or be out of order

Definition at line 100 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::execute(), and GrayLambertSurfaceRadiationBase::GrayLambertSurfaceRadiationBase().

◆ _fixed_side_temperature

std::vector<const Function *> GrayLambertSurfaceRadiationBase::_fixed_side_temperature

protectedinherited

side id to index map, side ids can have holes or be out of order

Definition at line 73 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::execute(), and GrayLambertSurfaceRadiationBase::GrayLambertSurfaceRadiationBase().

◆ _heat_flux_density

std::vector<Real> GrayLambertSurfaceRadiationBase::_heat_flux_density

protectedinherited

the heat flux density qdot

Definition at line 79 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::finalize(), and GrayLambertSurfaceRadiationBase::getSurfaceHeatFluxDensity().

◆ _n_sides

unsigned int GrayLambertSurfaceRadiationBase::_n_sides

protectedinherited

number of active boundary ids

Definition at line 64 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::finalize(), ConstantViewFactorSurfaceRadiation::initialize(), GrayLambertSurfaceRadiationBase::initialize(), setViewFactors(), ConstantViewFactorSurfaceRadiation::setViewFactors(), and GrayLambertSurfaceRadiationBase::threadJoin().

◆ _radiosity

std::vector<Real> GrayLambertSurfaceRadiationBase::_radiosity

protectedinherited

the radiosity of each surface

Definition at line 76 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::finalize(), and GrayLambertSurfaceRadiationBase::getSurfaceRadiosity().

◆ _side_id_index

std::map<BoundaryID, unsigned int> GrayLambertSurfaceRadiationBase::_side_id_index

protectedinherited

side id to index map, side ids can have holes or be out of order

Definition at line 88 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::execute(), GrayLambertSurfaceRadiationBase::getSurfaceEmissivity(), GrayLambertSurfaceRadiationBase::getSurfaceHeatFluxDensity(), GrayLambertSurfaceRadiationBase::getSurfaceIDs(), GrayLambertSurfaceRadiationBase::getSurfaceIrradiation(), GrayLambertSurfaceRadiationBase::getSurfaceRadiosity(), GrayLambertSurfaceRadiationBase::getSurfaceTemperature(), GrayLambertSurfaceRadiationBase::getViewFactor(), and GrayLambertSurfaceRadiationBase::GrayLambertSurfaceRadiationBase().

◆ _side_temperature

std::vector<Real> GrayLambertSurfaceRadiationBase::_side_temperature

protectedinherited

the average temperature: this could be important for adiabatic walls

Definition at line 82 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::execute(), GrayLambertSurfaceRadiationBase::finalize(), GrayLambertSurfaceRadiationBase::getSurfaceTemperature(), GrayLambertSurfaceRadiationBase::initialize(), and GrayLambertSurfaceRadiationBase::threadJoin().

◆ _side_type

std::vector<enum RAD_BND_TYPE> GrayLambertSurfaceRadiationBase::_side_type

protectedinherited

the type of the side, allows lookup index -> type

Definition at line 85 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::execute(), GrayLambertSurfaceRadiationBase::finalize(), and GrayLambertSurfaceRadiationBase::GrayLambertSurfaceRadiationBase().

◆ _sigma_stefan_boltzmann

const Real GrayLambertSurfaceRadiationBase::_sigma_stefan_boltzmann

protectedinherited

Stefan-Boltzmann constant.

Definition at line 61 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::execute(), and GrayLambertSurfaceRadiationBase::finalize().

◆ _surface_irradiation

std::vector<Real> GrayLambertSurfaceRadiationBase::_surface_irradiation

protectedinherited

the irradiation into each surface

Definition at line 97 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::finalize(), and GrayLambertSurfaceRadiationBase::getSurfaceIrradiation().

◆ _temperature

const VariableValue& GrayLambertSurfaceRadiationBase::_temperature

protectedinherited

the coupled temperature variable

Definition at line 67 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::execute().

◆ _view_factors

std::vector<std::vector<Real> > GrayLambertSurfaceRadiationBase::_view_factors

protectedinherited

the view factors which are set by setViewFactors by derived classes

Definition at line 106 of file GrayLambertSurfaceRadiationBase.h.

Referenced by GrayLambertSurfaceRadiationBase::finalize(), GrayLambertSurfaceRadiationBase::getViewFactor(), ConstantViewFactorSurfaceRadiation::initialize(), and GrayLambertSurfaceRadiationBase::initialize().

The documentation for this class was generated from the following files:

heat_conduction/include/userobjects/ViewFactorObjectSurfaceRadiation.h
heat_conduction/src/userobjects/ViewFactorObjectSurfaceRadiation.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ RAD_BND_TYPE

Constructor & Destructor Documentation

◆ ViewFactorObjectSurfaceRadiation()

Member Function Documentation

◆ execute()

◆ finalize()

◆ getSurfaceEmissivity()

◆ getSurfaceHeatFluxDensity()

◆ getSurfaceIDs()

◆ getSurfaceIrradiation()

◆ getSurfaceRadiosity()

◆ getSurfaceTemperature()

◆ getViewFactor()

◆ initialize()

◆ setViewFactors()

◆ threadJoin()

◆ validParams()

Member Data Documentation

◆ _adiabatic_side_ids

◆ _areas

◆ _beta

◆ _emissivity

◆ _fixed_side_id_index

◆ _fixed_side_temperature

◆ _heat_flux_density

◆ _n_sides

◆ _radiosity

◆ _side_id_index

◆ _side_temperature

◆ _side_type

◆ _sigma_stefan_boltzmann

◆ _surface_irradiation

◆ _temperature

◆ _view_factors