Boundary condition for convective heat flux where temperature and heat transfer coefficient are given by auxiliary variables. More...

#include <CoupledConvectiveHeatFluxBC.h>

Inheritance diagram for CoupledConvectiveHeatFluxBC:

[legend]

Public Member Functions
	CoupledConvectiveHeatFluxBC (const InputParameters &parameters)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

Protected Attributes
unsigned int	_n_components
	The number of components. More...

std::vector< const VariableValue * >	_T_infinity
	Far-field temperatue fields for each component. More...

std::vector< const VariableValue * >	_htc
	Convective heat transfer coefficient. More...

std::vector< const VariableValue * >	_alpha
	Volume fraction of individual phase. More...

Detailed Description

Boundary condition for convective heat flux where temperature and heat transfer coefficient are given by auxiliary variables.

Typically used in multi-app coupling scenario. It is possible to couple in a vector variable where each entry corresponds to a "phase".

Definition at line 24 of file CoupledConvectiveHeatFluxBC.h.

Constructor & Destructor Documentation

◆ CoupledConvectiveHeatFluxBC()

CoupledConvectiveHeatFluxBC::CoupledConvectiveHeatFluxBC ( const InputParameters & parameters )

Definition at line 30 of file CoupledConvectiveHeatFluxBC.C.

   : IntegratedBC(parameters), _n_components(coupledComponents("T_infinity"))
 {
   if (coupledComponents("alpha") != _n_components)
     paramError(
         "alpha",
         "The number of coupled components does not match the number of `T_infinity` components.");
   if (coupledComponents("htc") != _n_components)
     paramError(
         "htc",
         "The number of coupled components does not match the number of `T_infinity` components.");
  
   _htc.resize(_n_components);
   _T_infinity.resize(_n_components);
   _alpha.resize(_n_components);
   for (std::size_t c = 0; c < _n_components; c++)
   {
     _htc[c] = &coupledValue("htc", c);
     _T_infinity[c] = &coupledValue("T_infinity", c);
     _alpha[c] = &coupledValue("alpha", c);
   }
 }

Member Function Documentation

◆ computeQpJacobian()

Real CoupledConvectiveHeatFluxBC::computeQpJacobian ( )

protectedvirtual

Definition at line 63 of file CoupledConvectiveHeatFluxBC.C.

 {
   Real dq = 0;
   for (std::size_t c = 0; c < _n_components; c++)
     dq += (*_alpha[c])[_qp] * (*_htc[c])[_qp] * _phi[_j][_qp];
   return _test[_i][_qp] * dq;
 }

◆ computeQpResidual()

Real CoupledConvectiveHeatFluxBC::computeQpResidual ( )

protectedvirtual

Definition at line 54 of file CoupledConvectiveHeatFluxBC.C.

 {
   Real q = 0;
   for (std::size_t c = 0; c < _n_components; c++)
     q += (*_alpha[c])[_qp] * (*_htc[c])[_qp] * (_u[_qp] - (*_T_infinity[c])[_qp]);
   return _test[_i][_qp] * q;
 }

◆ validParams()

InputParameters CoupledConvectiveHeatFluxBC::validParams ( )

static

Definition at line 17 of file CoupledConvectiveHeatFluxBC.C.

 {
   InputParameters params = IntegratedBC::validParams();
   params.addClassDescription(
       "Convective heat transfer boundary condition with temperature and heat "
       "transfer coefficent given by auxiliary variables.");
   params.addCoupledVar("alpha", 1., "Volume fraction of components");
   params.addRequiredCoupledVar("T_infinity", "Field holding far-field temperature");
   params.addRequiredCoupledVar("htc", "Heat transfer coefficient");
  
   return params;
 }