Momentum equation boundary condition used when pressure is not integrated by parts, i.e. More...

#include <NSMomentumInviscidNoPressureImplicitFlowBC.h>

Inheritance diagram for NSMomentumInviscidNoPressureImplicitFlowBC:

Public Member Functions
	NSMomentumInviscidNoPressureImplicitFlowBC (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

virtual Real	computeQpOffDiagJacobian (unsigned jvar)

Real	pressureQpResidualHelper (Real pressure)

Real	pressureQpJacobianHelper (unsigned var_number)

Real	convectiveQpResidualHelper (Real rhou_udotn)

Real	convectiveQpJacobianHelper (unsigned var_number)

bool	isNSVariable (unsigned var)

unsigned	mapVarNumber (unsigned var)

Protected Attributes
const unsigned	_component

NSPressureDerivs< NSMomentumInviscidBC >	_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

Momentum equation boundary condition used when pressure is not integrated by parts, i.e.

when there is "no pressure" term on the boundary. In this case the flow is not specified either, so that the entire term is handled implicitly.

Definition at line 26 of file NSMomentumInviscidNoPressureImplicitFlowBC.h.

Constructor & Destructor Documentation

◆ NSMomentumInviscidNoPressureImplicitFlowBC()

NSMomentumInviscidNoPressureImplicitFlowBC::NSMomentumInviscidNoPressureImplicitFlowBC ( const InputParameters & parameters )

Definition at line 24 of file NSMomentumInviscidNoPressureImplicitFlowBC.C.

   : NSMomentumInviscidBC(parameters)
 {
 }

Member Function Documentation

◆ computeQpJacobian()

Real NSMomentumInviscidNoPressureImplicitFlowBC::computeQpJacobian ( )

protectedvirtual

Definition at line 46 of file NSMomentumInviscidNoPressureImplicitFlowBC.C.

 {
   // There is no Jacobian for the pressure term when the pressure is specified,
   // so all we have left is the convective part.  The on-diagonal variable number
   // is _component+1
   return convectiveQpJacobianHelper(_component + 1);
 }

◆ computeQpOffDiagJacobian()

Real NSMomentumInviscidNoPressureImplicitFlowBC::computeQpOffDiagJacobian ( unsigned jvar )

protectedvirtual

Definition at line 55 of file NSMomentumInviscidNoPressureImplicitFlowBC.C.

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

◆ computeQpResidual()

Real NSMomentumInviscidNoPressureImplicitFlowBC::computeQpResidual ( )

protectedvirtual

Definition at line 31 of file NSMomentumInviscidNoPressureImplicitFlowBC.C.

 {
   // Velocity vector object
   RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
   // Velocity vector dotted with normal
   Real u_dot_n = vel * _normals[_qp];
  
   // The current value of the vector (rho*u)(u.n)
   RealVectorValue rhou_udotn = u_dot_n * _rho[_qp] * vel;
  
   return convectiveQpResidualHelper(rhou_udotn(_component));
 }

◆ convectiveQpJacobianHelper()

Real NSMomentumInviscidBC::convectiveQpJacobianHelper ( unsigned var_number )

protectedinherited

Definition at line 64 of file NSMomentumInviscidBC.C.

 {
   // Velocity vector object
   RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
   // Variable to store convective contribution to boundary integral.
   Real conv_term = 0.0;
  
   // Inviscid components
   switch (var_number)
   {
     case 0: // density
       // Note: the minus sign here is correct, it comes from differentiating wrt U_0
       // (rho) which is in the denominator.
       conv_term = -vel(_component) * (vel * _normals[_qp]) * _phi[_j][_qp] * _test[_i][_qp];
       break;
  
     case 1:
     case 2:
     case 3: // momentums
       if (var_number - 1 == _component)
         // See Eqn. (68) from the notes for the inviscid boundary terms
         conv_term = ((vel * _normals[_qp]) + vel(_component) * _normals[_qp](_component)) *
                     _phi[_j][_qp] * _test[_i][_qp];
       else
         // off-diagonal
         conv_term =
             vel(_component) * _normals[_qp](var_number - 1) * _phi[_j][_qp] * _test[_i][_qp];
       break;
  
     case 4: // energy
       // No derivative wrt energy
       conv_term = 0.0;
       break;
  
     default:
       mooseError("Shouldn't get here!");
       break;
   }
  
   // Return the result.  We could return it directly from the switch statement, but this is
   // convenient for printing...
   return conv_term;
 }

Referenced by computeQpJacobian(), NSMomentumInviscidSpecifiedPressureBC::computeQpJacobian(), NSMomentumInviscidSpecifiedPressureBC::computeQpOffDiagJacobian(), and computeQpOffDiagJacobian().

◆ convectiveQpResidualHelper()

Real NSMomentumInviscidBC::convectiveQpResidualHelper ( Real rhou_udotn )

protectedinherited

Definition at line 52 of file NSMomentumInviscidBC.C.

 {
   // n . (rho*uu) . v = rho*(u.n)*(u.v) = (rho*u)(u.n) . v
  
   // The "inviscid" contribution: (rho*u)(u.n) . v
   Real conv_term = rhou_udotn * _test[_i][_qp];
  
   // Return value, or print it first if debugging...
   return conv_term;
 }

Referenced by NSMomentumInviscidSpecifiedNormalFlowBC::computeQpResidual(), computeQpResidual(), and NSMomentumInviscidSpecifiedPressureBC::computeQpResidual().

◆ 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;
 }

◆ pressureQpJacobianHelper()

Real NSMomentumInviscidBC::pressureQpJacobianHelper ( unsigned var_number )

protectedinherited

Definition at line 45 of file NSMomentumInviscidBC.C.

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

Referenced by NSMomentumInviscidSpecifiedNormalFlowBC::computeQpJacobian(), and NSMomentumInviscidSpecifiedNormalFlowBC::computeQpOffDiagJacobian().

◆ pressureQpResidualHelper()

Real NSMomentumInviscidBC::pressureQpResidualHelper ( Real pressure )

protectedinherited

Definition at line 33 of file NSMomentumInviscidBC.C.

 {
   // n . (Ip) . v
  
   // The pressure contribution: p * n(component) * phi_i
   Real press_term = pressure * _normals[_qp](_component) * _test[_i][_qp];
  
   // Return value, or print it first if debugging...
   return press_term;
 }