CNSFVHLLCMomentumBC.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 "CNSFVHLLCMomentumBC.h"

registerMooseObject("NavierStokesApp", CNSFVHLLCSpecifiedMassFluxAndTemperatureMomentumBC);
registerMooseObject("NavierStokesApp", CNSFVHLLCSpecifiedPressureMomentumBC);

template <typename T>
void
CNSFVHLLCMomentumBC<T>::addCommonParams(InputParameters & params)
{
  MooseEnum momentum_component("x=0 y=1 z=2");
  params.addRequiredParam<MooseEnum>(
      "momentum_component",
      momentum_component,
      "The component of the momentum equation that this kernel applies to.");
}

template <typename T>
InputParameters
CNSFVHLLCMomentumBC<T>::validParams()
{
  InputParameters params = T::validParams();
  addCommonParams(params);
  params.addClassDescription("Implements the momentum boundary flux portion of the free-flow HLLC "
                             "discretization given specified mass fluxes and fluid temperature");
  return params;
}

template <>
InputParameters
CNSFVHLLCMomentumBC<CNSFVHLLCSpecifiedPressureBC>::validParams()
{
  InputParameters params = CNSFVHLLCSpecifiedPressureBC::validParams();
  addCommonParams(params);
  params.addClassDescription("Implements the momentum boundary flux portion of the free-flow HLLC "
                             "discretization given specified pressure");
  return params;
}

template <typename T>
CNSFVHLLCMomentumBC<T>::CNSFVHLLCMomentumBC(const InputParameters & params)
  : T(params), _index(this->template getParam<MooseEnum>("momentum_component"))
{
}

template <typename T>
ADReal
CNSFVHLLCMomentumBC<T>::fluxElem()
{
  return this->_normal_speed_elem * this->_rho_elem[this->_qp] *
             this->_vel_elem[this->_qp](_index) +
         this->_normal(_index) * this->_pressure_elem[this->_qp];
}

template <typename T>
ADReal
CNSFVHLLCMomentumBC<T>::fluxBoundary()
{
  return this->_normal_speed_boundary * this->_rho_boundary * this->_vel_boundary(_index) +
         this->_normal(_index) * this->_pressure_boundary;
}

template <typename T>
ADReal
CNSFVHLLCMomentumBC<T>::hllcElem()
{
  auto vel_nonnormal = this->_vel_elem[this->_qp] - this->_normal_speed_elem * this->_normal;
  return this->_normal(_index) * this->_SM + vel_nonnormal(_index);

  // For some reason, the below expression doesn't give as good results as the
  // above one.
  // return this->_normal(_index) * (_SM - this->_normal_speed_elem) +
  // this->_vel_elem[this->_qp](_index);
}

template <typename T>
ADReal
CNSFVHLLCMomentumBC<T>::hllcBoundary()
{
  auto vel_nonnormal = this->_vel_boundary - this->_normal_speed_boundary * this->_normal;
  return this->_normal(_index) * this->_SM + vel_nonnormal(_index);

  // For some reason, the below expression doesn't give as good results as the
  // above one.
  // return this->_normal(_index) * (_SM - this->_normal_speed_boundary) +
  // this->_vel_elem[this->_qp](_index);
}

template <typename T>
ADReal
CNSFVHLLCMomentumBC<T>::conservedVariableElem()
{
  return this->_rho_elem[this->_qp] * this->_vel_elem[this->_qp](_index);
}

template <typename T>
ADReal
CNSFVHLLCMomentumBC<T>::conservedVariableBoundary()
{
  return this->_rho_boundary * this->_vel_boundary(_index);
}

template class CNSFVHLLCMomentumBC<CNSFVHLLCSpecifiedMassFluxAndTemperatureBC>;
template class CNSFVHLLCMomentumBC<CNSFVHLLCSpecifiedPressureBC>;