NSStagnationTemperatureBC.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

// Navier-Stokes includes
#include "NSStagnationTemperatureBC.h"
#include "NS.h"

// FluidProperties includes
#include "IdealGasFluidProperties.h"

registerMooseObject("NavierStokesApp", NSStagnationTemperatureBC);

InputParameters
NSStagnationTemperatureBC::validParams()
{
  InputParameters params = NSStagnationBC::validParams();
  params.addClassDescription("This Dirichlet condition imposes the condition T_0 = T_0_desired.");
  params.addRequiredCoupledVar(NS::temperature, "temperature");
  params.addRequiredParam<Real>("desired_stagnation_temperature", "");
  return params;
}

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

Real
NSStagnationTemperatureBC::computeQpResidual()
{
  // 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;
}