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

InputParameters
CNSFVHLLC::validParams()
{
  return CNSFVHLLCBase::validParams();
}

CNSFVHLLC::CNSFVHLLC(const InputParameters & params) : CNSFVHLLCBase(params) {}

ADReal
CNSFVHLLC::computeQpResidual()
{
  _normal_speed_elem = _normal * _vel_elem[_qp];
  _normal_speed_neighbor = _normal * _vel_neighbor[_qp];
  auto wave_speeds = waveSpeed(hllcData(), _normal);
  _SL = std::move(wave_speeds[0]);
  _SM = std::move(wave_speeds[1]);
  _SR = std::move(wave_speeds[2]);
  if (_SL >= 0)
    return fluxElem();
  else if (_SR <= 0)
    return fluxNeighbor();
  else
  {
    if (_SM >= 0)
    {
      ADReal f = _rho_elem[_qp] * (_SL - _normal_speed_elem) / (_SL - _SM);
      return fluxElem() + _SL * (f * hllcElem() - conservedVariableElem());
    }
    else
    {
      ADReal f = _rho_neighbor[_qp] * (_SR - _normal_speed_neighbor) / (_SR - _SM);
      return fluxNeighbor() + _SR * (f * hllcNeighbor() - conservedVariableNeighbor());
    }
  }
  mooseError("Should never get here");
}