#include <NSMassInviscidFlux.h>

Inheritance diagram for NSMassInviscidFlux:

Public Member Functions
	NSMassInviscidFlux (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 &	_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 NSMassInviscidFlux.h.

Constructor & Destructor Documentation

◆ NSMassInviscidFlux()

NSMassInviscidFlux::NSMassInviscidFlux ( const InputParameters & parameters )

Definition at line 23 of file NSMassInviscidFlux.C.

23 : NSKernel(parameters) {}

Member Function Documentation

◆ computeQpJacobian()

Real NSMassInviscidFlux::computeQpJacobian ( )

protectedvirtual

Definition at line 35 of file NSMassInviscidFlux.C.

 {
   // This seems weird at first glance, but remember we have to differentiate
   // wrt the *conserved* variables
   //
   // [ U_0 ] = [ rho       ]
   // [ U_1 ] = [ rho * u_1 ]
   // [ U_2 ] = [ rho * u_2 ]
   // [ U_3 ] = [ rho * u_3 ]
   // [ U_4 ] = [ rho * E   ]
   //
   // and the inviscid mass flux residual, in terms of these variables, is:
   //
   // f(U) = ( U_k * dphi_i/dx_k ), summation over k=1,2,3
   //
   // ie. does not depend on U_0, the on-diagonal Jacobian component.
   return 0.0;
 }

◆ computeQpOffDiagJacobian()

Real NSMassInviscidFlux::computeQpOffDiagJacobian ( unsigned int jvar )

protectedvirtual

Definition at line 55 of file NSMassInviscidFlux.C.

 {
   if (isNSVariable(jvar))
   {
     // Map jvar into the variable m for our problem, regardless of
     // how Moose has numbered things.
     unsigned int m = mapVarNumber(jvar);
  
     switch (m)
     {
       // Don't handle the on-diagonal case here
       // case 0: // density
       case 1:
       case 2:
       case 3: // momentums
         return -_phi[_j][_qp] * _grad_test[_i][_qp](m - 1);
  
       case 4: // energy
         return 0.0;
  
       default:
         mooseError("Should not get here!");
         break;
     }
   }
   else
     return 0.0;
 }

◆ computeQpResidual()

Real NSMassInviscidFlux::computeQpResidual ( )

protectedvirtual

Definition at line 26 of file NSMassInviscidFlux.C.

 {
   const RealVectorValue mom(_rho_u[_qp], _rho_v[_qp], _rho_w[_qp]);
  
   // -(rho*U) * grad(phi), negative sign comes from integration-by-parts
   return -(mom * _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(), computeQpOffDiagJacobian(), NSEnergyInviscidFlux::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;
 }