This Dirichlet condition imposes the condition T_0 = T_0_desired, where T_0 is the stagnation temperature, defined as: T_0 = T * (1 + (gam-1)/2 * M^2) More...

#include <NSStagnationTemperatureBC.h>

Inheritance diagram for NSStagnationTemperatureBC:

[legend]

Public Member Functions
	NSStagnationTemperatureBC (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual ()

Protected Attributes
const VariableValue &	_temperature

const Real	_desired_stagnation_temperature

const VariableValue &	_mach

const IdealGasFluidProperties &	_fp

Detailed Description

This Dirichlet condition imposes the condition T_0 = T_0_desired, where T_0 is the stagnation temperature, defined as: T_0 = T * (1 + (gam-1)/2 * M^2)

Definition at line 26 of file NSStagnationTemperatureBC.h.

Constructor & Destructor Documentation

◆ NSStagnationTemperatureBC()

NSStagnationTemperatureBC::NSStagnationTemperatureBC ( const InputParameters & parameters )

Definition at line 30 of file NSStagnationTemperatureBC.C.

   : NSStagnationBC(parameters),
     _temperature(coupledValue(NS::temperature)),
     _desired_stagnation_temperature(getParam<Real>("desired_stagnation_temperature"))
 {
 }

Member Function Documentation

◆ computeQpResidual()

Real NSStagnationTemperatureBC::computeQpResidual ( )

protectedvirtual

Definition at line 38 of file NSStagnationTemperatureBC.C.

 {
   // T_0 = T*(1 + 0.5*(gam-1)*M^2)
   Real computed_stagnation_temperature =
       _temperature[_qp] * (1. + 0.5 * (_fp.gamma() - 1.) * _mach[_qp] * _mach[_qp]);
 
   // Return the difference between the current solution's stagnation temperature
   // and the desired.  The Dirichlet condition asserts that these should be equal.
   return computed_stagnation_temperature - _desired_stagnation_temperature;
 }