NSEnergyWeakStagnationBC.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

#include "NSEnergyWeakStagnationBC.h"

// FluidProperties includes
#include "IdealGasFluidProperties.h"

registerMooseObject("NavierStokesApp", NSEnergyWeakStagnationBC);

InputParameters
NSEnergyWeakStagnationBC::validParams()
{
  InputParameters params = NSWeakStagnationBaseBC::validParams();
  params.addClassDescription("The inviscid energy BC term with specified normal flow.");
  return params;
}

NSEnergyWeakStagnationBC::NSEnergyWeakStagnationBC(const InputParameters & parameters)
  : NSWeakStagnationBaseBC(parameters)
{
}

Real
NSEnergyWeakStagnationBC::computeQpResidual()
{
  // Compute stagnation values
  Real T_s = 0.0, p_s = 0.0, rho_s = 0.0;
  staticValues(T_s, p_s, rho_s);

  // And get velocity magnitude, squared
  Real velmag2 = this->velmag2();

  // Compute static total energy, E_s
  Real E_s = _fp.cv() * T_s + 0.5 * velmag2;

  // Compute the product rho_s * H_s (H_s = static enthalpy)
  Real rhoH_s = rho_s * E_s + p_s;

  // rho_s * H_s * |u| * (s.n) * phi_i
  return rhoH_s * std::sqrt(velmag2) * this->sdotn() * _test[_i][_qp];
}

Real
NSEnergyWeakStagnationBC::computeQpJacobian()
{
  // TODO
  return 0.0;
}

Real
NSEnergyWeakStagnationBC::computeQpOffDiagJacobian(unsigned /*jvar*/)
{
  // TODO
  return 0.0;
}