The inviscid energy BC term with specified density and velocity components. More...

#include <NSEnergyInviscidSpecifiedDensityAndVelocityBC.h>

Inheritance diagram for NSEnergyInviscidSpecifiedDensityAndVelocityBC:

Public Member Functions
	NSEnergyInviscidSpecifiedDensityAndVelocityBC (const InputParameters &parameters)

virtual	~NSEnergyInviscidSpecifiedDensityAndVelocityBC ()

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

Real	_specified_density

Real	_specified_u

Real	_specified_v

Real	_specified_w

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

Detailed Description

The inviscid energy BC term with specified density and velocity components.

This was experimental code and did not really work out, do not use!

Definition at line 24 of file NSEnergyInviscidSpecifiedDensityAndVelocityBC.h.

Constructor & Destructor Documentation

◆ NSEnergyInviscidSpecifiedDensityAndVelocityBC()

NSEnergyInviscidSpecifiedDensityAndVelocityBC::NSEnergyInviscidSpecifiedDensityAndVelocityBC ( const InputParameters & parameters )

Definition at line 40 of file NSEnergyInviscidSpecifiedDensityAndVelocityBC.C.

   : NSEnergyInviscidBC(parameters),
     _pressure(coupledValue(NS::pressure)),
     _specified_density(getParam<Real>("specified_density")),
     _specified_u(getParam<Real>("specified_u")),
     _specified_v(getParam<Real>("specified_v")),
     _specified_w(getParam<Real>("specified_w"))
 {
 }

◆ ~NSEnergyInviscidSpecifiedDensityAndVelocityBC()

virtual NSEnergyInviscidSpecifiedDensityAndVelocityBC::~NSEnergyInviscidSpecifiedDensityAndVelocityBC ( )

inlinevirtual

Definition at line 29 of file NSEnergyInviscidSpecifiedDensityAndVelocityBC.h.

29 {}

Member Function Documentation

◆ computeQpJacobian()

Real NSEnergyInviscidSpecifiedDensityAndVelocityBC::computeQpJacobian ( )

protectedvirtual

Definition at line 60 of file NSEnergyInviscidSpecifiedDensityAndVelocityBC.C.

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

◆ computeQpOffDiagJacobian()

Real NSEnergyInviscidSpecifiedDensityAndVelocityBC::computeQpOffDiagJacobian ( unsigned jvar )

protectedvirtual

Definition at line 68 of file NSEnergyInviscidSpecifiedDensityAndVelocityBC.C.

 {
   // TODO
   // return computeJacobianHelper(mapVarNumber(jvar));
   return 0.;
 }

◆ computeQpResidual()

Real NSEnergyInviscidSpecifiedDensityAndVelocityBC::computeQpResidual ( )

protectedvirtual

Definition at line 52 of file NSEnergyInviscidSpecifiedDensityAndVelocityBC.C.

 {
   return qpResidualHelper(_specified_density,
                           RealVectorValue(_specified_u, _specified_v, _specified_w),
                           _pressure[_qp]);
 }

◆ 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 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 NSEnergyInviscidUnspecifiedBC::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(), NSEnergyInviscidUnspecifiedBC::computeQpResidual(), NSEnergyInviscidSpecifiedPressureBC::computeQpResidual(), NSEnergyInviscidSpecifiedNormalFlowBC::computeQpResidual(), and 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];
 }