The inviscid energy BC term with specified pressure. More...

#include <NSEnergyInviscidUnspecifiedBC.h>

Inheritance diagram for NSEnergyInviscidUnspecifiedBC:

Public Member Functions
	NSEnergyInviscidUnspecifiedBC (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 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 pressure.

Definition at line 23 of file NSEnergyInviscidUnspecifiedBC.h.

Constructor & Destructor Documentation

◆ NSEnergyInviscidUnspecifiedBC()

NSEnergyInviscidUnspecifiedBC::NSEnergyInviscidUnspecifiedBC ( const InputParameters & parameters )

Definition at line 25 of file NSEnergyInviscidUnspecifiedBC.C.

   : NSEnergyInviscidBC(parameters), _pressure(coupledValue(NS::pressure))
 {
 }

Member Function Documentation

◆ computeJacobianHelper()

Real NSEnergyInviscidUnspecifiedBC::computeJacobianHelper ( unsigned var_number )

private

Definition at line 58 of file NSEnergyInviscidUnspecifiedBC.C.

 {
   // Velocity vector object
   RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
   // Normal component of velocity
   Real un = vel * _normals[_qp];
  
   // When both u.n and pressure are unspecified, all 3 Jacobian terms apply.
   // See base class for details.
   return qpJacobianTermA(var_number, _pressure[_qp]) + qpJacobianTermB(var_number, un) +
          qpJacobianTermC(var_number, un);
 }

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ computeQpJacobian()

Real NSEnergyInviscidUnspecifiedBC::computeQpJacobian ( )

protectedvirtual

Definition at line 43 of file NSEnergyInviscidUnspecifiedBC.C.

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

◆ computeQpOffDiagJacobian()

Real NSEnergyInviscidUnspecifiedBC::computeQpOffDiagJacobian ( unsigned jvar )

protectedvirtual

Definition at line 49 of file NSEnergyInviscidUnspecifiedBC.C.

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

◆ computeQpResidual()

Real NSEnergyInviscidUnspecifiedBC::computeQpResidual ( )

protectedvirtual

Definition at line 31 of file NSEnergyInviscidUnspecifiedBC.C.

 {
   // Velocity vector object
   RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
   // Normal component
   Real un = vel * _normals[_qp];
  
   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 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(), computeJacobianHelper(), NSEnergyInviscidSpecifiedBC::computeJacobianHelper(), and NSEnergyInviscidSpecifiedNormalFlowBC::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 computeJacobianHelper(), and NSEnergyInviscidSpecifiedNormalFlowBC::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(), computeQpResidual(), NSEnergyInviscidSpecifiedNormalFlowBC::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];
 }