#include <NSEnergyInviscidFlux.h>

Inheritance diagram for NSEnergyInviscidFlux:

Public Member Functions
	NSEnergyInviscidFlux (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 VariableValue &	_enthalpy

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

Detailed Description

Definition at line 20 of file NSEnergyInviscidFlux.h.

Constructor & Destructor Documentation

◆ NSEnergyInviscidFlux()

NSEnergyInviscidFlux::NSEnergyInviscidFlux ( const InputParameters & parameters )

Definition at line 29 of file NSEnergyInviscidFlux.C.

   : NSKernel(parameters), _enthalpy(coupledValue(NS::enthalpy))
 {
 }

Member Function Documentation

◆ computeQpJacobian()

Real NSEnergyInviscidFlux::computeQpJacobian ( )

protectedvirtual

Definition at line 55 of file NSEnergyInviscidFlux.C.

 {
   // Derivative of this kernel wrt rho*E
   const RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
   // Ratio of specific heats
   const Real gam = _fp.gamma();
  
   // -gamma * phi_j * (U*grad(phi_i))
   return -gam * _phi[_j][_qp] * (vel * _grad_test[_i][_qp]);
 }

◆ computeQpOffDiagJacobian()

Real NSEnergyInviscidFlux::computeQpOffDiagJacobian ( unsigned int jvar )

protectedvirtual

Definition at line 68 of file NSEnergyInviscidFlux.C.

 {
   if (isNSVariable(jvar))
   {
     RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
     Real V2 = vel.norm_sq();
  
     // Ratio of specific heats
     const Real gam = _fp.gamma();
  
     // Derivative wrt density
     if (jvar == _rho_var_number)
       return -((0.5 * (gam - 1) * V2 - _enthalpy[_qp]) * _phi[_j][_qp] *
                (vel * _grad_test[_i][_qp]));
  
     // Derivatives wrt momentums
     else if ((jvar == _rhou_var_number) || (jvar == _rhov_var_number) || (jvar == _rhow_var_number))
     {
       // Map jvar into jlocal = {0,1,2}, regardless of how Moose has numbered things.
       unsigned jlocal = 0;
  
       if (jvar == _rhov_var_number)
         jlocal = 1;
       else if (jvar == _rhow_var_number)
         jlocal = 2;
  
       // Scale the velocity vector by the scalar (1-gamma)*vel(jlocal)
       vel *= (1.0 - gam) * vel(jlocal);
  
       // Add in the enthalpy in the jlocal'th entry
       vel(jlocal) += _enthalpy[_qp];
  
       // Return -1 * (vel * grad(phi_i)) * phi_j
       return -(vel * _grad_test[_i][_qp]) * _phi[_j][_qp];
     }
  
     else
     {
       std::ostringstream oss;
       oss << "Invalid jvar=" << jvar << " requested!\n"
           << "Did not match:\n"
           << " _rho_var_number =" << _rho_var_number << "\n"
           << " _rhou_var_number=" << _rhou_var_number << "\n"
           << " _rhov_var_number=" << _rhov_var_number << "\n"
           << " _rhow_var_number=" << _rhow_var_number << std::endl;
       mooseError(oss.str());
     }
   }
   else
     return 0.0;
 }

◆ computeQpResidual()

Real NSEnergyInviscidFlux::computeQpResidual ( )

protectedvirtual

Definition at line 35 of file NSEnergyInviscidFlux.C.

 {
   // H = total enthalpy = E + P/rho
   // => rho * u * H = rho * u ( E + P/rho)
   //                =       u ( rho*E + P)
  
   // velocity vector
   RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
   // Multiply vector U by the scalar value (rho*E + P) to get rho * U * H
   // vel *= (_u[_qp] + _pressure[_qp]);
  
   // Multiply velocity vector by the scalar (rho * H)
   vel *= (_rho[_qp] * _enthalpy[_qp]);
  
   // Return -1 * vel * grad(phi_i)
   return -(vel * _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.

 {
   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;
 }

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

◆ mapVarNumber()

unsigned NSKernel::mapVarNumber ( unsigned var )

protectedinherited

Definition at line 88 of file NSKernel.C.

 {
   // 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;
 }