This Dirichlet condition imposes the condition p_0 = p_0_desired, where p_0 is the stagnation pressure, defined as: p_0 = p * (1 + (gam-1)/2 * M^2)^(gam/(gam-1)) More...

#include <NSStagnationPressureBC.h>

Inheritance diagram for NSStagnationPressureBC:

[legend]

Public Member Functions
	NSStagnationPressureBC (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual ()

Protected Attributes
const VariableValue &	_pressure

const Real	_desired_stagnation_pressure

const VariableValue &	_mach

const IdealGasFluidProperties &	_fp

Detailed Description

This Dirichlet condition imposes the condition p_0 = p_0_desired, where p_0 is the stagnation pressure, defined as: p_0 = p * (1 + (gam-1)/2 * M^2)^(gam/(gam-1))

Definition at line 26 of file NSStagnationPressureBC.h.

Constructor & Destructor Documentation

◆ NSStagnationPressureBC()

NSStagnationPressureBC::NSStagnationPressureBC ( const InputParameters & parameters )

Definition at line 31 of file NSStagnationPressureBC.C.

   : NSStagnationBC(parameters),
     _pressure(coupledValue(NS::pressure)),
     _desired_stagnation_pressure(getParam<Real>("desired_stagnation_pressure"))
 {
 }

Member Function Documentation

◆ computeQpResidual()

Real NSStagnationPressureBC::computeQpResidual ( )

protectedvirtual

Definition at line 39 of file NSStagnationPressureBC.C.

 {
   // p_0 = p*(1 + 0.5*(gam-1)*M^2)^(gam/(gam-1))
   const Real computed_stagnation_pressure =
       _pressure[_qp] * std::pow(1. + 0.5 * (_fp.gamma() - 1.) * _mach[_qp] * _mach[_qp],
                                 _fp.gamma() / (_fp.gamma() - 1.));
  
   // Return the difference between the current solution's stagnation pressure
   // and the desired.  The Dirichlet condition asserts that these should be equal.
   return computed_stagnation_pressure - _desired_stagnation_pressure;
 }

Member Data Documentation

◆ _desired_stagnation_pressure

const Real NSStagnationPressureBC::_desired_stagnation_pressure

protected

Definition at line 41 of file NSStagnationPressureBC.h.

Referenced by computeQpResidual().

◆ _fp

const IdealGasFluidProperties& NSStagnationBC::_fp

protectedinherited

Definition at line 36 of file NSStagnationBC.h.

Referenced by NSStagnationTemperatureBC::computeQpResidual(), and computeQpResidual().

◆ _mach

const VariableValue& NSStagnationBC::_mach

protectedinherited

Definition at line 33 of file NSStagnationBC.h.

Referenced by NSStagnationTemperatureBC::computeQpResidual(), and computeQpResidual().

◆ _pressure

const VariableValue& NSStagnationPressureBC::_pressure