This class corresponds to the inviscid part of the "natural" boundary condition for the momentum equations, i.e. More...

#include <NSMomentumInviscidBC.h>

Inheritance diagram for NSMomentumInviscidBC:

[legend]

Public Member Functions
	NSMomentumInviscidBC (const InputParameters &parameters)

Protected Member Functions
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

Friends
template<class U >
class	NSPressureDerivs

Detailed Description

This class corresponds to the inviscid part of the "natural" boundary condition for the momentum equations, i.e.

int_{Gamma} n . (rho*uu + Ip) . v

While this kernel implements the convective and pressure term residuals and jacobians, it does not itself implement any of the computeQp* functions. For that, use one of the derived classes: 1.) NSMomentumInviscidSpecifiedPressureBC 2.) NSMomentumInviscidSpecifiedNormalFlowBC 3.) NSMomentumInviscidNoPressureImplicitFlowBC

The first kernel above would be used for a subsonic outflow BC in the Euler or Navier-Stokes equations in which one physical value (the pressure) is specified. In this case, the residual and Jacobian contrbutions of the n.(rho*u)(u.v) term are computed and added to the matrix/rhs. For the pressure term, the residual contribution due to the specified pressure is computed but there is no corresponding Jacobian entry since the value is given.

The second kernel above would be used if, instead of the pressure, the value of the vector (rho*u)(u.n) is given. This situation is not common, but a special case of it, u.n=0, is very common for implementing a free-slip boundary: in that situation the pressure is an unknown and u.n = 0 would be imposed by the kernel (though in this trivial case you could omit the term completely).

The third kernel above would be used when the pressure term has not been integrated by parts in the momentum equations, and therefore there is no pressure in the boundary term at all.

We note that other combinations are also theoretically possible (e.g. unspecified pressure and normal flow, fully-specified pressure and normal flow) however they have not yet been implemented since I'm not sure if they are physically-relevant.

Definition at line 60 of file NSMomentumInviscidBC.h.

Constructor & Destructor Documentation

◆ NSMomentumInviscidBC()

NSMomentumInviscidBC::NSMomentumInviscidBC ( const InputParameters & parameters )

Definition at line 24 of file NSMomentumInviscidBC.C.

   : NSIntegratedBC(parameters),
     _component(getParam<unsigned>("component")),
     // Object for computing deriviatives of pressure
     _pressure_derivs(*this)
 {
 }

Member Function Documentation

◆ convectiveQpJacobianHelper()

Real NSMomentumInviscidBC::convectiveQpJacobianHelper ( unsigned var_number )

protected

Definition at line 64 of file NSMomentumInviscidBC.C.

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

 {
   // 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;
 }

◆ convectiveQpResidualHelper()

Real NSMomentumInviscidBC::convectiveQpResidualHelper ( Real rhou_udotn )

protected

Definition at line 52 of file NSMomentumInviscidBC.C.

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

 {
   // 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;
 }

◆ 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 )

protected

Definition at line 45 of file NSMomentumInviscidBC.C.

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

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

◆ pressureQpResidualHelper()

Real NSMomentumInviscidBC::pressureQpResidualHelper ( Real pressure )

protected

Definition at line 33 of file NSMomentumInviscidBC.C.

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

 {
   // 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;
 }

Friends And Related Function Documentation

◆ NSPressureDerivs

template<class U >

friend class NSPressureDerivs

friend

Definition at line 76 of file NSMomentumInviscidBC.h.

Member Data Documentation

◆ _component

const unsigned NSMomentumInviscidBC::_component

protected

Definition at line 68 of file NSMomentumInviscidBC.h.

Referenced by NSMomentumInviscidSpecifiedNormalFlowBC::computeQpJacobian(), NSMomentumInviscidSpecifiedPressureBC::computeQpJacobian(), NSMomentumInviscidNoPressureImplicitFlowBC::computeQpJacobian(), NSMomentumInviscidNoPressureImplicitFlowBC::computeQpResidual(), NSMomentumInviscidSpecifiedPressureBC::computeQpResidual(), convectiveQpJacobianHelper(), pressureQpJacobianHelper(), and pressureQpResidualHelper().

◆ _dynamic_viscosity

const MaterialProperty<Real>& NSIntegratedBC::_dynamic_viscosity

protectedinherited

Definition at line 59 of file NSIntegratedBC.h.

◆ _fp

const IdealGasFluidProperties& NSIntegratedBC::_fp

protectedinherited

Definition at line 63 of file NSIntegratedBC.h.

Referenced by NSEnergyWeakStagnationBC::computeQpResidual(), NSEnergyInviscidBC::qpResidualHelper(), and NSWeakStagnationBaseBC::staticValues().

◆ _grad_rho

const VariableGradient& NSIntegratedBC::_grad_rho

protectedinherited

Definition at line 46 of file NSIntegratedBC.h.

Referenced by NSEnergyViscousBC::NSEnergyViscousBC().

◆ _grad_rho_E

const VariableGradient& NSIntegratedBC::_grad_rho_E

protectedinherited

Definition at line 50 of file NSIntegratedBC.h.

Referenced by NSEnergyViscousBC::NSEnergyViscousBC().

◆ _grad_rho_u

const VariableGradient& NSIntegratedBC::_grad_rho_u

protectedinherited

Definition at line 47 of file NSIntegratedBC.h.

Referenced by NSEnergyViscousBC::NSEnergyViscousBC().

◆ _grad_rho_v

const VariableGradient& NSIntegratedBC::_grad_rho_v

protectedinherited

Definition at line 48 of file NSIntegratedBC.h.

Referenced by NSEnergyViscousBC::NSEnergyViscousBC().

◆ _grad_rho_w

const VariableGradient& NSIntegratedBC::_grad_rho_w

protectedinherited

Definition at line 49 of file NSIntegratedBC.h.

Referenced by NSEnergyViscousBC::NSEnergyViscousBC().

◆ _pressure_derivs

NSPressureDerivs<NSMomentumInviscidBC> NSMomentumInviscidBC::_pressure_derivs

protected

Definition at line 72 of file NSMomentumInviscidBC.h.

Referenced by pressureQpJacobianHelper().

◆ _rho

const VariableValue& NSIntegratedBC::_rho

protectedinherited

Definition at line 40 of file NSIntegratedBC.h.

Referenced by NSEnergyViscousBC::computeQpOffDiagJacobian(), NSMomentumInviscidSpecifiedPressureBC::computeQpResidual(), NSMomentumInviscidNoPressureImplicitFlowBC::computeQpResidual(), and NSEnergyInviscidBC::qpJacobianTermA().

◆ _rho_E

const VariableValue& NSIntegratedBC::_rho_E

protectedinherited

Definition at line 44 of file NSIntegratedBC.h.

Referenced by NSEnergyInviscidBC::qpJacobianTermA(), and NSEnergyInviscidBC::qpResidualHelper().

◆ _rho_u

const VariableValue& NSIntegratedBC::_rho_u

protectedinherited

Definition at line 41 of file NSIntegratedBC.h.

Referenced by NSEnergyViscousBC::computeQpOffDiagJacobian(), and NSMassUnspecifiedNormalFlowBC::computeQpResidual().

◆ _rho_v

const VariableValue& NSIntegratedBC::_rho_v

protectedinherited

Definition at line 42 of file NSIntegratedBC.h.

Referenced by NSEnergyViscousBC::computeQpOffDiagJacobian(), and NSMassUnspecifiedNormalFlowBC::computeQpResidual().

◆ _rho_var_number

unsigned NSIntegratedBC::_rho_var_number

protectedinherited

Definition at line 52 of file NSIntegratedBC.h.

Referenced by NSIntegratedBC::isNSVariable(), and NSIntegratedBC::mapVarNumber().

◆ _rho_w

const VariableValue& NSIntegratedBC::_rho_w

protectedinherited

Definition at line 43 of file NSIntegratedBC.h.

Referenced by NSEnergyViscousBC::computeQpOffDiagJacobian(), and NSMassUnspecifiedNormalFlowBC::computeQpResidual().

◆ _rhoE_var_number

unsigned NSIntegratedBC::_rhoE_var_number

protectedinherited

Definition at line 56 of file NSIntegratedBC.h.

Referenced by NSIntegratedBC::isNSVariable(), and NSIntegratedBC::mapVarNumber().

◆ _rhou_var_number

unsigned NSIntegratedBC::_rhou_var_number

protectedinherited

Definition at line 53 of file NSIntegratedBC.h.

Referenced by NSIntegratedBC::isNSVariable(), and NSIntegratedBC::mapVarNumber().

◆ _rhov_var_number

unsigned NSIntegratedBC::_rhov_var_number

protectedinherited

Definition at line 54 of file NSIntegratedBC.h.

Referenced by NSIntegratedBC::isNSVariable(), and NSIntegratedBC::mapVarNumber().

◆ _rhow_var_number

unsigned NSIntegratedBC::_rhow_var_number

protectedinherited

Definition at line 55 of file NSIntegratedBC.h.

Referenced by NSIntegratedBC::isNSVariable(), and NSIntegratedBC::mapVarNumber().

◆ _u_vel

const VariableValue& NSIntegratedBC::_u_vel

protectedinherited

Definition at line 36 of file NSIntegratedBC.h.

Referenced by NSEnergyInviscidSpecifiedPressureBC::computeJacobianHelper(), NSEnergyInviscidUnspecifiedBC::computeJacobianHelper(), NSEnergyInviscidUnspecifiedBC::computeQpResidual(), NSEnergyInviscidSpecifiedPressureBC::computeQpResidual(), NSMomentumInviscidNoPressureImplicitFlowBC::computeQpResidual(), NSMomentumInviscidSpecifiedPressureBC::computeQpResidual(), NSPenalizedNormalFlowBC::computeQpResidual(), NSEnergyViscousBC::computeQpResidual(), convectiveQpJacobianHelper(), NSEnergyInviscidBC::qpJacobianTermA(), and NSWeakStagnationBaseBC::velmag2().

◆ _v_vel

const VariableValue& NSIntegratedBC::_v_vel

protectedinherited

Definition at line 37 of file NSIntegratedBC.h.

Referenced by NSEnergyInviscidSpecifiedPressureBC::computeJacobianHelper(), NSEnergyInviscidUnspecifiedBC::computeJacobianHelper(), NSEnergyInviscidUnspecifiedBC::computeQpResidual(), NSEnergyInviscidSpecifiedPressureBC::computeQpResidual(), NSMomentumInviscidNoPressureImplicitFlowBC::computeQpResidual(), NSMomentumInviscidSpecifiedPressureBC::computeQpResidual(), NSPenalizedNormalFlowBC::computeQpResidual(), NSEnergyViscousBC::computeQpResidual(), convectiveQpJacobianHelper(), NSEnergyInviscidBC::qpJacobianTermA(), and NSWeakStagnationBaseBC::velmag2().

◆ _viscous_stress_tensor

const MaterialProperty<RealTensorValue>& NSIntegratedBC::_viscous_stress_tensor

protectedinherited

Definition at line 60 of file NSIntegratedBC.h.

Referenced by NSEnergyViscousBC::computeQpOffDiagJacobian(), NSMomentumViscousBC::computeQpResidual(), and NSEnergyViscousBC::computeQpResidual().

◆ _w_vel

const VariableValue& NSIntegratedBC::_w_vel

protectedinherited

Definition at line 38 of file NSIntegratedBC.h.

Referenced by NSEnergyInviscidSpecifiedPressureBC::computeJacobianHelper(), NSEnergyInviscidUnspecifiedBC::computeJacobianHelper(), NSEnergyInviscidUnspecifiedBC::computeQpResidual(), NSEnergyInviscidSpecifiedPressureBC::computeQpResidual(), NSMomentumInviscidNoPressureImplicitFlowBC::computeQpResidual(), NSMomentumInviscidSpecifiedPressureBC::computeQpResidual(), NSPenalizedNormalFlowBC::computeQpResidual(), NSEnergyViscousBC::computeQpResidual(), convectiveQpJacobianHelper(), NSEnergyInviscidBC::qpJacobianTermA(), and NSWeakStagnationBaseBC::velmag2().

The documentation for this class was generated from the following files:

navier_stokes/include/bcs/NSMomentumInviscidBC.h
navier_stokes/src/bcs/NSMomentumInviscidBC.C

Public Member Functions

Protected Member Functions

Protected Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ NSMomentumInviscidBC()

Member Function Documentation

◆ convectiveQpJacobianHelper()

◆ convectiveQpResidualHelper()

◆ isNSVariable()

◆ mapVarNumber()

◆ pressureQpJacobianHelper()

◆ pressureQpResidualHelper()

Friends And Related Function Documentation

◆ NSPressureDerivs

Member Data Documentation

◆ _component

◆ _dynamic_viscosity

◆ _fp

◆ _grad_rho

◆ _grad_rho_E

◆ _grad_rho_u

◆ _grad_rho_v

◆ _grad_rho_w

◆ _pressure_derivs

◆ _rho

◆ _rho_E

◆ _rho_u

◆ _rho_v

◆ _rho_var_number

◆ _rho_w

◆ _rhoE_var_number

◆ _rhou_var_number

◆ _rhov_var_number

◆ _rhow_var_number

◆ _u_vel

◆ _v_vel

◆ _viscous_stress_tensor

◆ _w_vel