The inviscid energy BC term with specified normal flow. More...

#include <NSEnergyInviscidSpecifiedNormalFlowBC.h>

Inheritance diagram for NSEnergyInviscidSpecifiedNormalFlowBC:

Public Member Functions
	NSEnergyInviscidSpecifiedNormalFlowBC (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

virtual Real	computeQpOffDiagJacobian (unsigned jvar)

Real	qpResidualHelper (Real pressure, Real un)

Real	qpResidualHelper (Real rho, RealVectorValue u, Real pressure)

Real	qpJacobianTermA (unsigned var_number, Real pressure)

Real	qpJacobianTermB (unsigned var_number, Real un)

Real	qpJacobianTermC (unsigned var_number, Real un)

bool	isNSVariable (unsigned var)

unsigned	mapVarNumber (unsigned var)

Protected Attributes
const VariableValue &	_pressure

const Real	_un

const VariableValue &	_temperature

NSPressureDerivs< NSEnergyInviscidBC >	_pressure_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 (unsigned var_number)

Detailed Description

The inviscid energy BC term with specified normal flow.

Definition at line 23 of file NSEnergyInviscidSpecifiedNormalFlowBC.h.

Constructor & Destructor Documentation

◆ NSEnergyInviscidSpecifiedNormalFlowBC()

NSEnergyInviscidSpecifiedNormalFlowBC::NSEnergyInviscidSpecifiedNormalFlowBC ( const InputParameters & parameters )

Definition at line 26 of file NSEnergyInviscidSpecifiedNormalFlowBC.C.

   : NSEnergyInviscidBC(parameters), _pressure(coupledValue(NS::pressure)), _un(getParam<Real>("un"))
 {
 }

Member Function Documentation

◆ computeJacobianHelper()

Real NSEnergyInviscidSpecifiedNormalFlowBC::computeJacobianHelper ( unsigned var_number )

private

Definition at line 54 of file NSEnergyInviscidSpecifiedNormalFlowBC.C.

 {
   // For specified u.n, term "A" is zero, see base class for details.
   return qpJacobianTermB(var_number, _un) + qpJacobianTermC(var_number, _un);
 }

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ computeQpJacobian()

Real NSEnergyInviscidSpecifiedNormalFlowBC::computeQpJacobian ( )

protectedvirtual

Definition at line 39 of file NSEnergyInviscidSpecifiedNormalFlowBC.C.

 {
   return computeJacobianHelper(/*on-diagonal variable is energy=*/4);
 }

◆ computeQpOffDiagJacobian()

Real NSEnergyInviscidSpecifiedNormalFlowBC::computeQpOffDiagJacobian ( unsigned jvar )

protectedvirtual

Definition at line 45 of file NSEnergyInviscidSpecifiedNormalFlowBC.C.

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

◆ computeQpResidual()

Real NSEnergyInviscidSpecifiedNormalFlowBC::computeQpResidual ( )

protectedvirtual

Definition at line 33 of file NSEnergyInviscidSpecifiedNormalFlowBC.C.

 {
   return qpResidualHelper(_pressure[_qp], _un);
 }

◆ isNSVariable()

bool NSIntegratedBC::isNSVariable ( unsigned var )

protectedinherited

Definition at line 80 of file NSIntegratedBC.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;
 }

◆ mapVarNumber()

unsigned NSIntegratedBC::mapVarNumber ( unsigned var )

protectedinherited

Definition at line 90 of file NSIntegratedBC.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 int mapped_var_number;
  
   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;
 }

◆ qpJacobianTermA()

Real NSEnergyInviscidBC::qpJacobianTermA	(	unsigned	var_number,
		Real	pressure
	)

protectedinherited

Definition at line 55 of file NSEnergyInviscidBC.C.

 {
   Real result = 0.0;
  
   switch (var_number)
   {
     case 0: // density
     {
       // Velocity vector object
       RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
       result = -(vel * _normals[_qp]);
       break;
     }
  
     case 1:
     case 2:
     case 3: // momentums
       result = _normals[_qp](var_number - 1);
       break;
  
     case 4: // energy
       result = 0.;
       break;
  
     default:
       mooseError("Shouldn't get here!");
       break;
   }
  
   // Notice the division by _rho[_qp] here.  This comes from taking the
   // derivative wrt to either density or momentum.
   return (_rho_E[_qp] + pressure) / _rho[_qp] * result * _phi[_j][_qp] * _test[_i][_qp];
 }

Referenced by NSEnergyInviscidSpecifiedPressureBC::computeJacobianHelper(), and NSEnergyInviscidUnspecifiedBC::computeJacobianHelper().

◆ qpJacobianTermB()

Real NSEnergyInviscidBC::qpJacobianTermB	(	unsigned	var_number,
		Real	un
	)

protectedinherited

Definition at line 92 of file NSEnergyInviscidBC.C.

 {
   Real result = 0.0;
   switch (var_number)
   {
     case 0: // density
     case 1:
     case 2:
     case 3: // momentums
     {
       result = 0.;
       break;
     }
  
     case 4: // energy
     {
       result = _phi[_j][_qp] * un * _test[_i][_qp];
       break;
     }
  
     default:
       mooseError("Shouldn't get here!");
       break;
   }
  
   return result;
 }

Referenced by NSEnergyInviscidSpecifiedPressureBC::computeJacobianHelper(), NSEnergyInviscidUnspecifiedBC::computeJacobianHelper(), NSEnergyInviscidSpecifiedBC::computeJacobianHelper(), and computeJacobianHelper().

◆ qpJacobianTermC()

Real NSEnergyInviscidBC::qpJacobianTermC	(	unsigned	var_number,
		Real	un
	)

protectedinherited

Definition at line 122 of file NSEnergyInviscidBC.C.

 {
   return _pressure_derivs.get_grad(var_number) * _phi[_j][_qp] * un * _test[_i][_qp];
 }

Referenced by NSEnergyInviscidUnspecifiedBC::computeJacobianHelper(), and computeJacobianHelper().

◆ qpResidualHelper() [1/2]

Real NSEnergyInviscidBC::qpResidualHelper	(	Real	pressure,
		Real	un
	)

protectedinherited

Definition at line 37 of file NSEnergyInviscidBC.C.

 {
   return (_rho_E[_qp] + pressure) * un * _test[_i][_qp];
 }

Referenced by NSEnergyInviscidSpecifiedBC::computeQpResidual(), NSEnergyInviscidUnspecifiedBC::computeQpResidual(), computeQpResidual(), NSEnergyInviscidSpecifiedPressureBC::computeQpResidual(), and NSEnergyInviscidSpecifiedDensityAndVelocityBC::computeQpResidual().

◆ qpResidualHelper() [2/2]

Real NSEnergyInviscidBC::qpResidualHelper	(	Real	rho,
		RealVectorValue	u,
		Real	pressure
	)

protectedinherited

Definition at line 43 of file NSEnergyInviscidBC.C.

 {
   // return (rho*(cv*_temperature[_qp] + 0.5*u.norm_sq()) + pressure) * (u*_normals[_qp]) *
   // _test[_i][_qp];
   // We can also expand pressure in terms of rho... does this make a difference?
   // Then we don't use the input pressure value.
   return rho * (_fp.gamma() * _fp.cv() * _temperature[_qp] + 0.5 * u.norm_sq()) *
          (u * _normals[_qp]) * _test[_i][_qp];
 }