#include <SideSetHeatTransferMaterial.h>

Inheritance diagram for SideSetHeatTransferMaterial:

Public Member Functions
	SideSetHeatTransferMaterial (const InputParameters &parameters)

Protected Member Functions
virtual void	computeQpProperties () override

Protected Attributes
const Function &	_kgap
	Material conductivity. More...

const VariableValue *	_Tk
	Variable to compute conductivity. More...

const Function &	_dgap
	Gap width. More...

const Function &	_Tb
	Gap bulk temperature. More...

const Function &	_hp
	Heat transfer coefficients. More...

const Function &	_hm

const Function &	_eps_p
	Surface emissivities. More...

const Function &	_eps_m

const Function *	_rho_p
	Surface reflectivities. More...

const Function *	_rho_m

MaterialProperty< Real > &	_cond
	Material property handles. More...

MaterialProperty< Real > &	_Tbulk

MaterialProperty< Real > &	_h_master

MaterialProperty< Real > &	_h_neighbor

MaterialProperty< Real > &	_emmissivity_eff_master

MaterialProperty< Real > &	_emmissivity_eff_neighbor

const Real	_sigma
	Stefan-Boltzmann constant in W/cm^2-K^4. More...

Detailed Description

Definition at line 23 of file SideSetHeatTransferMaterial.h.

Constructor & Destructor Documentation

◆ SideSetHeatTransferMaterial()

SideSetHeatTransferMaterial::SideSetHeatTransferMaterial ( const InputParameters & parameters )

Definition at line 45 of file SideSetHeatTransferMaterial.C.

   : InterfaceMaterial(parameters),
     _kgap(isParamValid("conductivity_temperature_function")
               ? getFunction("conductivity_temperature_function")
               : getFunction("conductivity")),
     _Tk(isCoupled("gap_temperature") ? &coupledValue("gap_temperature") : nullptr),
     _dgap(getFunction("gap_length")),
     _Tb(getFunction("Tbulk")),
     _hp(getFunction("h_master")),
     _hm(getFunction("h_neighbor")),
     _eps_p(getFunction("emissivity_master")),
     _eps_m(getFunction("emissivity_neighbor")),
     _rho_p(isParamValid("reflectivity_master") ? &getFunction("reflectivity_master") : nullptr),
     _rho_m(isParamValid("reflectivity_neighbor") ? &getFunction("reflectivity_neighbor") : nullptr),
     _cond(declareProperty<Real>("gap_conductance")),
     _Tbulk(declareProperty<Real>("gap_Tbulk")),
     _h_master(declareProperty<Real>("gap_h_master")),
     _h_neighbor(declareProperty<Real>("gap_h_neighbor")),
     _emmissivity_eff_master(declareProperty<Real>("gap_emissivity_eff_master")),
     _emmissivity_eff_neighbor(declareProperty<Real>("gap_emissivity_eff_neighbor")),
     _sigma(SIGMA)
 {
   if ((parameters.isParamSetByUser("conductivity") ||
        isParamValid("conductivity_temperature_function")) &&
       !parameters.isParamSetByUser("gap_length"))
     paramError("gap_length", "gap_length not set, but conduction term requested.");
   if (parameters.isParamSetByUser("gap_length") &&
       !(parameters.isParamSetByUser("conductivity") ||
         isParamValid("conductivity_temperature_function")))
     paramError("conductivity", "conductivity not set, but conduction term requested.");
   if (isParamValid("conductivity_temperature_function") &&
       parameters.isParamSetByUser("conductivity"))
     paramError("conductivity",
                "Cannot specify both conductivity and conductivity_temperature_function.");
   if (isParamValid("conductivity_temperature_function") && !_Tk)
     paramError("gap_temperature",
                "Variable specification for temp needed if specifying a temperature dependent "
                "conductivity.");
  
   if (parameters.isParamSetByUser("h_master") && !parameters.isParamSetByUser("h_neighbor"))
     paramError("h_neighbor", "h_neighbor not set, but convection term requested.");
   if (parameters.isParamSetByUser("h_neighbor") && !parameters.isParamSetByUser("h_master"))
     paramError("h_master", "h_master not set, but convection term requested.");
  
   if (parameters.isParamSetByUser("emissivity_master") &&
       !parameters.isParamSetByUser("emissivity_neighbor"))
     paramError("emissivity_neighbor", "emissivity_neighbor not set, but radiation term requested");
   if (parameters.isParamSetByUser("emissivity_neighbor") &&
       !parameters.isParamSetByUser("emissivity_master"))
     paramError("emissivity_master", "emissivity_master not set, but radiation term requested");
 }

Member Function Documentation

◆ computeQpProperties()

void SideSetHeatTransferMaterial::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 98 of file SideSetHeatTransferMaterial.C.

 {
   // Conductance defined as k_{gap}/\delta
   Real tqp = (_Tk ? (*_Tk)[_qp] : _t);
   _cond[_qp] = _kgap.value(tqp, _q_point[_qp]) / _dgap.value(_t, _q_point[_qp]);
  
   // Convection parameters
   _Tbulk[_qp] = _Tb.value(_t, _q_point[_qp]);
   _h_master[_qp] = _hp.value(_t, _q_point[_qp]);
   _h_neighbor[_qp] = _hm.value(_t, _q_point[_qp]);
  
   // If reflectivity not provided, assume 1 - epsilon
   Real rhop = (_rho_p ? (*_rho_p).value(_t, _q_point[_qp]) : 1.0 - _eps_p.value(_t, _q_point[_qp]));
   Real rhom = (_rho_m ? (*_rho_m).value(_t, _q_point[_qp]) : 1.0 - _eps_m.value(_t, _q_point[_qp]));
  
   // Compute effective emissivity (master): \frac{\sigma\epsilon^+(1-\rho^-)}{1-rho^+\rho^-}
   _emmissivity_eff_master[_qp] = _sigma * _eps_p.value(_t, _q_point[_qp]) * (1.0 - rhom);
   if (_emmissivity_eff_master[_qp] != 0.0) // Making sure we don't devide by zero
     _emmissivity_eff_master[_qp] /= 1.0 - rhop * rhom;
  
   // Compute effective emissivity (neighbor): \frac{\sigma\epsilon^-(1-\rho^+)}{1-rho^+\rho^-}
   _emmissivity_eff_neighbor[_qp] = _sigma * _eps_m.value(_t, _q_point[_qp]) * (1.0 - rhop);
   if (_emmissivity_eff_neighbor[_qp] != 0.0) // Making sure we don't devide by zero
     _emmissivity_eff_neighbor[_qp] /= 1.0 - rhop * rhom;
 }