This class computes the residual and Jacobian contributions for the incompressible Navier-Stokes temperature (energy) equation. More...

#include <INSTemperature.h>

Inheritance diagram for INSTemperature:

[legend]

Public Member Functions
	INSTemperature (const InputParameters &parameters)

virtual	~INSTemperature ()

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

virtual Real	computeQpOffDiagJacobian (unsigned jvar)

Protected Attributes
const VariableValue &	_u_vel

const VariableValue &	_v_vel

const VariableValue &	_w_vel

unsigned	_u_vel_var_number

unsigned	_v_vel_var_number

unsigned	_w_vel_var_number

const MaterialProperty< Real > &	_rho

const MaterialProperty< Real > &	_k

const MaterialProperty< Real > &	_cp

Detailed Description

This class computes the residual and Jacobian contributions for the incompressible Navier-Stokes temperature (energy) equation.

Definition at line 24 of file INSTemperature.h.

Constructor & Destructor Documentation

◆ INSTemperature()

INSTemperature::INSTemperature ( const InputParameters & parameters )

Definition at line 36 of file INSTemperature.C.

   : Kernel(parameters),
  
     // Coupled variables
     _u_vel(coupledValue("u")),
     _v_vel(_mesh.dimension() >= 2 ? coupledValue("v") : _zero),
     _w_vel(_mesh.dimension() == 3 ? coupledValue("w") : _zero),
  
     // Variable numberings
     _u_vel_var_number(coupled("u")),
     _v_vel_var_number(_mesh.dimension() >= 2 ? coupled("v") : libMesh::invalid_uint),
     _w_vel_var_number(_mesh.dimension() == 3 ? coupled("w") : libMesh::invalid_uint),
  
     // Material Properties
     _rho(getMaterialProperty<Real>("rho_name")),
     _k(getMaterialProperty<Real>("k_name")),
     _cp(getMaterialProperty<Real>("cp_name"))
 {
 }

◆ ~INSTemperature()

virtual INSTemperature::~INSTemperature ( )

inlinevirtual

Definition at line 29 of file INSTemperature.h.

29 {}

Member Function Documentation

◆ computeQpJacobian()

Real INSTemperature::computeQpJacobian ( )

protectedvirtual

Definition at line 73 of file INSTemperature.C.

 {
   RealVectorValue U(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
   Real convective_part = _rho[_qp] * _cp[_qp] * (U * _grad_phi[_j][_qp]) * _test[_i][_qp];
   Real conduction_part = _k[_qp] * (_grad_phi[_j][_qp] * _grad_test[_i][_qp]);
  
   return convective_part + conduction_part;
 }

◆ computeQpOffDiagJacobian()

Real INSTemperature::computeQpOffDiagJacobian ( unsigned jvar )

protectedvirtual

Definition at line 84 of file INSTemperature.C.

 {
   if (jvar == _u_vel_var_number)
   {
     Real convective_part = _rho[_qp] * _cp[_qp] * _phi[_j][_qp] * _grad_u[_qp](0) * _test[_i][_qp];
     return convective_part;
   }
  
   else if (jvar == _v_vel_var_number)
   {
     Real convective_part = _rho[_qp] * _cp[_qp] * _phi[_j][_qp] * _grad_u[_qp](1) * _test[_i][_qp];
     return convective_part;
   }
  
   else if (jvar == _w_vel_var_number)
   {
     Real convective_part = _rho[_qp] * _cp[_qp] * _phi[_j][_qp] * _grad_u[_qp](2) * _test[_i][_qp];
     return convective_part;
   }
   else
     return 0;
 }

◆ computeQpResidual()

Real INSTemperature::computeQpResidual ( )

protectedvirtual

Definition at line 57 of file INSTemperature.C.

 {
   // The convection part, rho * cp u.grad(T) * v.
   // Note: _u is the temperature variable, _grad_u is its gradient.
   Real convective_part = _rho[_qp] * _cp[_qp] *
                          (_u_vel[_qp] * _grad_u[_qp](0) + _v_vel[_qp] * _grad_u[_qp](1) +
                           _w_vel[_qp] * _grad_u[_qp](2)) *
                          _test[_i][_qp];
  
   // Thermal conduction part, k * grad(T) * grad(v)
   Real conduction_part = _k[_qp] * _grad_u[_qp] * _grad_test[_i][_qp];
  
   return convective_part + conduction_part;
 }