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;
};