DG kernel for interfacing diffusion between two variables on adjacent blocks. More...

#include <SideSetHeatTransferKernel.h>

Inheritance diagram for SideSetHeatTransferKernel:

Public Member Functions
	SideSetHeatTransferKernel (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual (Moose::DGResidualType type) override

virtual Real	computeQpJacobian (Moose::DGJacobianType type) override

Protected Attributes
const MaterialProperty< Real > &	_cond
	Conductivity of gap divided by effective gap width. More...

const VariableValue *	_Tbulk_var
	Bulk temperature of gap. More...

const MaterialProperty< Real > *	_Tbulk_mat

const MaterialProperty< Real > &	_hp
	Convective heat transfer coefficient (master face) More...

const MaterialProperty< Real > &	_hm
	Convective heat transfer coefficient (neighbor face) More...

const MaterialProperty< Real > &	_eps_p
	Master face effective emissivity \epsilon^+\sigma(1-\rho^-)/(1-\rho^+\rho^-) More...

const MaterialProperty< Real > &	_eps_m
	Neighbor face effective emissivity \epsilon^-\sigma(1-\rho^+)/(1-\rho^+\rho^-) More...

Detailed Description

DG kernel for interfacing diffusion between two variables on adjacent blocks.

Definition at line 23 of file SideSetHeatTransferKernel.h.

Constructor & Destructor Documentation

◆ SideSetHeatTransferKernel()

SideSetHeatTransferKernel::SideSetHeatTransferKernel ( const InputParameters & parameters )

Definition at line 51 of file SideSetHeatTransferKernel.C.

   : InterfaceKernel(parameters),
     _cond(getMaterialProperty<Real>("conductance")),
     _Tbulk_var(isParamValid("Tbulk_var") ? &coupledValue("Tbulk_var") : nullptr),
     _Tbulk_mat(_Tbulk_var ? nullptr : &getMaterialProperty<Real>("Tbulk_mat")),
     _hp(getMaterialProperty<Real>("h_master")),
     _hm(getMaterialProperty<Real>("h_neighbor")),
     _eps_p(getMaterialProperty<Real>("emissivity_eff_master")),
     _eps_m(getMaterialProperty<Real>("emissivity_eff_neighbor"))
 {
   if (parameters.isParamSetByUser("Tbulk_mat") && _Tbulk_var)
     paramError("Tbulk_var", "Both Tbulk_mat and Tbulk_var set by user, cannot use both.");
  
   if (_var.number() == _neighbor_var.number() && _var.isNodal())
     mooseError(
         "Variable and neighbor variable are the same, but they are not elemental variables.");
 }

Member Function Documentation

◆ computeQpJacobian()

Real SideSetHeatTransferKernel::computeQpJacobian ( Moose::DGJacobianType type )

overrideprotectedvirtual

Definition at line 125 of file SideSetHeatTransferKernel.C.

 {
   Real jac = 0;
  
   if (_cond[_qp] != 0.0) // Conduction
   {
     switch (type)
     {
       case Moose::ElementElement:
         jac += _cond[_qp] * _phi[_j][_qp] * _test[_i][_qp];
         break;
  
       case Moose::NeighborNeighbor:
         jac += _cond[_qp] * _phi_neighbor[_j][_qp] * _test_neighbor[_i][_qp];
         break;
  
       case Moose::ElementNeighbor:
         jac -= _cond[_qp] * _phi_neighbor[_j][_qp] * _test[_i][_qp];
         break;
  
       case Moose::NeighborElement:
         jac -= _cond[_qp] * _phi[_j][_qp] * _test_neighbor[_i][_qp];
         break;
     }
   }
  
   if (_hp[_qp] != 0.0 && _hm[_qp] != 0.0) // Convection
   {
     switch (type)
     {
       case Moose::ElementElement:
         jac += _hp[_qp] * _phi[_j][_qp] * _test[_i][_qp];
         break;
  
       case Moose::NeighborNeighbor:
         jac += _hm[_qp] * _phi_neighbor[_j][_qp] * _test_neighbor[_i][_qp];
         break;
  
       case Moose::NeighborElement:
       case Moose::ElementNeighbor:
         break;
     }
   }
  
   if (_eps_p[_qp] != 0.0 && _eps_m[_qp] != 0.0) // Radiation
   {
     switch (type)
     {
       case Moose::ElementElement:
         jac += 4.0 * _eps_p[_qp] * (_u[_qp] * _u[_qp] * _u[_qp]) * _phi[_j][_qp] * _test[_i][_qp];
         break;
  
       case Moose::NeighborNeighbor:
         jac += 4.0 * _eps_m[_qp] *
                (_neighbor_value[_qp] * _neighbor_value[_qp] * _neighbor_value[_qp]) *
                _phi_neighbor[_j][_qp] * _test_neighbor[_i][_qp];
         break;
  
       case Moose::ElementNeighbor:
         jac -= 4.0 * _eps_m[_qp] *
                (_neighbor_value[_qp] * _neighbor_value[_qp] * _neighbor_value[_qp]) *
                _phi_neighbor[_j][_qp] * _test[_i][_qp];
         break;
  
       case Moose::NeighborElement:
         jac -= 4.0 * _eps_p[_qp] * (_u[_qp] * _u[_qp] * _u[_qp]) * _phi[_j][_qp] *
                _test_neighbor[_i][_qp];
         break;
     }
   }
  
   return jac;
 }

◆ computeQpResidual()

Real SideSetHeatTransferKernel::computeQpResidual ( Moose::DGResidualType type )

overrideprotectedvirtual

Definition at line 70 of file SideSetHeatTransferKernel.C.

 {
   Real r = 0;
  
   if (_cond[_qp] != 0.0) // Conduction
   {
     Real jump = _u[_qp] - _neighbor_value[_qp];
     switch (type)
     {
       case Moose::Element:
         r += _cond[_qp] * jump * _test[_i][_qp];
         break;
  
       case Moose::Neighbor:
         r -= _cond[_qp] * jump * _test_neighbor[_i][_qp];
         break;
     }
   }
  
   if (_hp[_qp] != 0.0 && _hm[_qp] != 0.0) // Convection
   {
     Real Tb = (_Tbulk_var ? (*_Tbulk_var)[_qp] : (*_Tbulk_mat)[_qp]);
     switch (type)
     {
       case Moose::Element:
         r += _hp[_qp] * (_u[_qp] - Tb) * _test[_i][_qp];
         break;
  
       case Moose::Neighbor:
         r += _hm[_qp] * (_neighbor_value[_qp] - Tb) * _test_neighbor[_i][_qp];
         break;
     }
   }
  
   if (_eps_p[_qp] != 0.0 && _eps_m[_qp] != 0.0) // Radiation
   {
     Real Rp = _eps_p[_qp] * (_u[_qp] * _u[_qp] * _u[_qp] * _u[_qp]);
     Real Rm = _eps_m[_qp] * (_neighbor_value[_qp] * _neighbor_value[_qp] * _neighbor_value[_qp] *
                              _neighbor_value[_qp]);
     switch (type)
     {
       case Moose::Element:
         r += (Rp - Rm) * _test[_i][_qp];
         break;
  
       case Moose::Neighbor:
         r += (Rm - Rp) * _test_neighbor[_i][_qp];
         break;
     }
   }
  
   return r;
 }