The inviscid flux (convective + pressure terms) for the momentum conservation equations. More...

#include <NSMomentumInviscidFlux.h>

Inheritance diagram for NSMomentumInviscidFlux:

Public Member Functions
	NSMomentumInviscidFlux (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 &	_pressure

const unsigned int	_component

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 (unsigned int m)

Detailed Description

The inviscid flux (convective + pressure terms) for the momentum conservation equations.

Definition at line 24 of file NSMomentumInviscidFlux.h.

Constructor & Destructor Documentation

◆ NSMomentumInviscidFlux()

NSMomentumInviscidFlux::NSMomentumInviscidFlux ( const InputParameters & parameters )

Definition at line 33 of file NSMomentumInviscidFlux.C.

   : NSKernel(parameters),
     _pressure(coupledValue(NS::pressure)),
     _component(getParam<unsigned int>("component"))
 {
 }

Member Function Documentation

◆ computeJacobianHelper()

Real NSMomentumInviscidFlux::computeJacobianHelper ( unsigned int m )

private

Definition at line 74 of file NSMomentumInviscidFlux.C.

 {
   const RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
   // Ratio of specific heats
   const Real gam = _fp.gamma();
  
   switch (m)
   {
     case 0: // density
     {
       const Real V2 = vel.norm_sq();
       return vel(_component) * (vel * _grad_test[_i][_qp]) -
              0.5 * (gam - 1.0) * V2 * _grad_test[_i][_qp](_component);
     }
  
     case 1:
     case 2:
     case 3: // momentums
     {
       // Map m into m_local = {0,1,2}
       unsigned int m_local = m - 1;
  
       // Kronecker delta
       const Real delta_kl = (_component == m_local ? 1. : 0.);
  
       return -1.0 * (vel(_component) * _grad_test[_i][_qp](m_local) +
                      delta_kl * (vel * _grad_test[_i][_qp]) +
                      (1. - gam) * vel(m_local) * _grad_test[_i][_qp](_component)) *
              _phi[_j][_qp];
     }
  
     case 4: // energy
       return -1.0 * (gam - 1.0) * _phi[_j][_qp] * _grad_test[_i][_qp](_component);
   }
  
   mooseError("Shouldn't get here!");
 }

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ computeQpJacobian()

Real NSMomentumInviscidFlux::computeQpJacobian ( )

protectedvirtual

Definition at line 58 of file NSMomentumInviscidFlux.C.

 {
   // The on-diagonal entry corresponds to variable number _component+1.
   return computeJacobianHelper(_component + 1);
 }

◆ computeQpOffDiagJacobian()

Real NSMomentumInviscidFlux::computeQpOffDiagJacobian ( unsigned int jvar )

protectedvirtual

Definition at line 65 of file NSMomentumInviscidFlux.C.

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

◆ computeQpResidual()

Real NSMomentumInviscidFlux::computeQpResidual ( )

protectedvirtual

Definition at line 41 of file NSMomentumInviscidFlux.C.

 {
   // For _component = k,
  
   // (rho*u) * u_k = (rho*u_k) * u <- we write it this way
   RealVectorValue vec(_u[_qp] * _u_vel[_qp],  // (U_k) * u_1
                       _u[_qp] * _v_vel[_qp],  // (U_k) * u_2
                       _u[_qp] * _w_vel[_qp]); // (U_k) * u_3
  
   // (rho*u_k) * u + e_k * P [ e_k = unit vector in k-direction ]
   vec(_component) += _pressure[_qp];
  
   // -((rho*u_k) * u + e_k * P) * grad(phi)
   return -(vec * _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(), NSEnergyInviscidFlux::computeQpOffDiagJacobian(), NSMomentumInviscidFluxWithGradP::computeQpOffDiagJacobian(), 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;
 }