This class is responsible for computing residuals and Jacobian terms for the k * grad(T) * grad(phi) term in the Navier-Stokes energy equation. More...

#include <NSEnergyThermalFlux.h>

Inheritance diagram for NSEnergyThermalFlux:

[legend]

Public Member Functions
	NSEnergyThermalFlux (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

virtual Real	computeQpOffDiagJacobian (unsigned int jvar)

bool	isNSVariable (unsigned var)
	Helper functions for mapping Moose variable numberings into the "canonical" numbering for the compressible NS equations. More...

unsigned	mapVarNumber (unsigned var)

Protected Attributes
const VariableGradient &	_grad_temp

const MaterialProperty< Real > &	_thermal_conductivity

NSTemperatureDerivs< NSEnergyThermalFlux >	_temp_derivs

const VariableValue &	_u_vel

const VariableValue &	_v_vel

const VariableValue &	_w_vel

const VariableValue &	_rho

const VariableValue &	_rho_u

const VariableValue &	_rho_v

const VariableValue &	_rho_w

const VariableValue &	_rho_E

const VariableGradient &	_grad_rho

const VariableGradient &	_grad_rho_u

const VariableGradient &	_grad_rho_v

const VariableGradient &	_grad_rho_w

const VariableGradient &	_grad_rho_E

unsigned	_rho_var_number

unsigned	_rhou_var_number

unsigned	_rhov_var_number

unsigned	_rhow_var_number

unsigned	_rhoE_var_number

const MaterialProperty< Real > &	_dynamic_viscosity

const MaterialProperty< RealTensorValue > &	_viscous_stress_tensor

const IdealGasFluidProperties &	_fp

Private Member Functions
Real	computeJacobianHelper_value (unsigned var_number)

Private Attributes
std::vector< const VariableGradient * >	_gradU

Friends
template<class U >
class	NSTemperatureDerivs

Detailed Description

This class is responsible for computing residuals and Jacobian terms for the k * grad(T) * grad(phi) term in the Navier-Stokes energy equation.

Definition at line 27 of file NSEnergyThermalFlux.h.

Constructor & Destructor Documentation

◆ NSEnergyThermalFlux()

NSEnergyThermalFlux::NSEnergyThermalFlux ( const InputParameters & parameters )

Definition at line 28 of file NSEnergyThermalFlux.C.

   : NSKernel(parameters),
     _grad_temp(coupledGradient(NS::temperature)),
     _thermal_conductivity(getMaterialProperty<Real>("thermal_conductivity")),
     // Temperature derivative computing object
     _temp_derivs(*this)
 {
   // Store pointers to all variable gradients in a single vector.
   _gradU.resize(5);
   _gradU[0] = &_grad_rho;
   _gradU[1] = &_grad_rho_u;
   _gradU[2] = &_grad_rho_v;
   _gradU[3] = &_grad_rho_w;
   _gradU[4] = &_grad_rho_E;
 }

Member Function Documentation

◆ computeJacobianHelper_value()

Real NSEnergyThermalFlux::computeJacobianHelper_value ( unsigned var_number )

private

Definition at line 69 of file NSEnergyThermalFlux.C.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

 {
   // The value to return
   Real result = 0.0;
 
   // I used "ell" here as the loop counter since it matches the
   // "\ell" used in my LaTeX notes.
   for (unsigned int ell = 0; ell < 3; ++ell)
   {
     // Accumulate the first dot product term
     Real intermediate_result = _temp_derivs.get_grad(var_number) * _grad_phi[_j][_qp](ell);
 
     // Now accumulate the Hessian term
     Real hess_term = 0.0;
     for (unsigned n = 0; n < 5; ++n)
     {
       // hess_term += get_hess(m,n) * gradU[n](ell); // ideally... but you can't have a
       // vector<VariableGradient&> :-(
       hess_term += _temp_derivs.get_hess(var_number, n) *
                    (*_gradU[n])[_qp](ell); // dereference pointer to get value
     }
 
     // Accumulate the second dot product term
     intermediate_result += hess_term * _phi[_j][_qp];
 
     // Hit intermediate_result with the test function, accumulate in the final value
     result += intermediate_result * _grad_test[_i][_qp](ell);
   }
 
   // Return result, don't forget to multiply by "k"!
   return _thermal_conductivity[_qp] * result;
 }

◆ computeQpJacobian()

Real NSEnergyThermalFlux::computeQpJacobian ( )

protectedvirtual

Definition at line 52 of file NSEnergyThermalFlux.C.

 {
   // The "on-diagonal" Jacobian for the energy equation
   // corresponds to variable number 4.
   return computeJacobianHelper_value(/*var_number=*/4);
 }

◆ computeQpOffDiagJacobian()

Real NSEnergyThermalFlux::computeQpOffDiagJacobian ( unsigned int jvar )

protectedvirtual

Definition at line 60 of file NSEnergyThermalFlux.C.

 {
   if (isNSVariable(jvar))
     return computeJacobianHelper_value(mapVarNumber(jvar));
   else
     return 0.0;
 }

◆ computeQpResidual()

Real NSEnergyThermalFlux::computeQpResidual ( )

protectedvirtual

Definition at line 45 of file NSEnergyThermalFlux.C.

 {
   // k * grad(T) * grad(phi)
   return _thermal_conductivity[_qp] * (_grad_temp[_qp] * _grad_test[_i][_qp]);
 }

◆ isNSVariable()

bool NSKernel::isNSVariable ( unsigned var )

protectedinherited

Helper functions for mapping Moose variable numberings into the "canonical" numbering for the compressible NS equations.

Definition at line 78 of file NSKernel.C.

Referenced by NSSUPGEnergy::computeJacobianHelper(), NSSUPGMass::computeJacobianHelper(), NSSUPGMomentum::computeJacobianHelper(), NSMassInviscidFlux::computeQpOffDiagJacobian(), NSEnergyInviscidFlux::computeQpOffDiagJacobian(), NSMomentumInviscidFluxWithGradP::computeQpOffDiagJacobian(), NSMomentumInviscidFlux::computeQpOffDiagJacobian(), NSEnergyViscousFlux::computeQpOffDiagJacobian(), NSMomentumViscousFlux::computeQpOffDiagJacobian(), and computeQpOffDiagJacobian().

 {
   if (var == _rho_var_number || var == _rhou_var_number || var == _rhov_var_number ||
       var == _rhow_var_number || var == _rhoE_var_number)
     return true;
   else
     return false;
 }

◆ mapVarNumber()

unsigned NSKernel::mapVarNumber ( unsigned var )

protectedinherited

Definition at line 88 of file NSKernel.C.

Referenced by NSSUPGEnergy::computeJacobianHelper(), NSSUPGMass::computeJacobianHelper(), NSSUPGMomentum::computeJacobianHelper(), NSMassInviscidFlux::computeQpOffDiagJacobian(), NSMomentumInviscidFluxWithGradP::computeQpOffDiagJacobian(), NSEnergyViscousFlux::computeQpOffDiagJacobian(), NSMomentumInviscidFlux::computeQpOffDiagJacobian(), NSMomentumViscousFlux::computeQpOffDiagJacobian(), and computeQpOffDiagJacobian().

 {
   // Convert the Moose numbering to:
   // 0 for rho
   // 1 for rho*u
   // 2 for rho*v
   // 3 for rho*w
   // 4 for rho*e
   // regardless of the problem dimension, etc.
   unsigned mapped_var_number = 99;
 
   if (var == _rho_var_number)
     mapped_var_number = 0;
   else if (var == _rhou_var_number)
     mapped_var_number = 1;
   else if (var == _rhov_var_number)
     mapped_var_number = 2;
   else if (var == _rhow_var_number)
     mapped_var_number = 3;
   else if (var == _rhoE_var_number)
     mapped_var_number = 4;
   else
     mooseError("Invalid var!");
 
   return mapped_var_number;
 }