PCNSFVHLLC.h

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

#pragma once

#include "FVFluxKernel.h"

class HLLCUserObject;
class SinglePhaseFluidProperties;

class PCNSFVHLLC : public FVFluxKernel
{
public:
  static InputParameters validParams();
  PCNSFVHLLC(const InputParameters & params);

  static std::array<ADReal, 3> waveSpeed(const ADReal & rho_elem,
                                         const ADRealVectorValue & vel_elem,
                                         const ADReal & e_elem,
                                         Real eps_elem,
                                         const ADReal & rho_neighbor,
                                         const ADRealVectorValue & vel_neighbor,
                                         const ADReal & e_neighbor,
                                         Real eps_neighbor,
                                         const SinglePhaseFluidProperties & fluid,
                                         const ADRealVectorValue & normal);

protected:
  virtual ADReal computeQpResidual() override;

  virtual ADReal fluxElem() = 0;
  virtual ADReal fluxNeighbor() = 0;

  virtual ADReal hllcElem() = 0;
  virtual ADReal hllcNeighbor() = 0;

  virtual ADReal conservedVariableElem() = 0;
  virtual ADReal conservedVariableNeighbor() = 0;

  ADReal _SL;
  ADReal _SM;
  ADReal _SR;

  ADReal _normal_speed_elem;
  ADReal _normal_speed_neighbor;

  const SinglePhaseFluidProperties & _fluid;

  const ADMaterialProperty<Real> & _specific_internal_energy_elem;
  const ADMaterialProperty<Real> & _specific_internal_energy_neighbor;

  const ADMaterialProperty<Real> & _rho_et_elem;
  const ADMaterialProperty<Real> & _rho_et_neighbor;

  const ADMaterialProperty<RealVectorValue> & _vel_elem;
  const ADMaterialProperty<RealVectorValue> & _vel_neighbor;

  const ADMaterialProperty<Real> & _speed_elem;
  const ADMaterialProperty<Real> & _speed_neighbor;

  const ADMaterialProperty<Real> & _rho_elem;
  const ADMaterialProperty<Real> & _rho_neighbor;

  const ADMaterialProperty<Real> & _pressure_elem;
  const ADMaterialProperty<Real> & _pressure_neighbor;

  const MaterialProperty<Real> & _eps_elem;
  const MaterialProperty<Real> & _eps_neighbor;
};