Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #pragma once
      11             : 
      12             : #include "InternalSideUserObject.h"
      13             : #include "MooseMesh.h"
      14             : 
      15             : using side_type = std::pair<const Elem *, unsigned int>;
      16             : 
      17             : // Forward Declarations
      18             : class HLLCUserObject;
      19             : class SinglePhaseFluidProperties;
      20             : 
      21             : class HLLCUserObject : public InternalSideUserObject
      22             : {
      23             : public:
      24             :   static InputParameters validParams();
      25             : 
      26             :   HLLCUserObject(const InputParameters & parameters);
      27             : 
      28             :   virtual void initialSetup() override;
      29          34 :   virtual void initialize() override {}
      30             :   virtual void execute() override;
      31          28 :   virtual void finalize() override {}
      32             :   virtual void threadJoin(const UserObject & y) override;
      33             : 
      34             :   /// accessor for the wave speed
      35             :   std::vector<ADReal> waveSpeed(const Elem * elem, unsigned int side) const;
      36             : 
      37             :   /// Query whether this processor has data for the provided element and side
      38             :   bool hasData(const Elem * elem, unsigned int side) const;
      39             : 
      40             : protected:
      41             :   /// helper function for returning the FaceInfo object for an elem/side pair
      42             :   const FaceInfo & faceInfoHelper(const Elem * elem, unsigned int side) const;
      43             : 
      44             :   /// quadrature point dummy
      45             :   unsigned int _qp = 0;
      46             : 
      47             :   /// fluid properties
      48             :   const SinglePhaseFluidProperties & _fluid;
      49             : 
      50             :   /// FV face info from MooseMesh
      51             :   const std::vector<const FaceInfo *> & _face_info;
      52             : 
      53             :   /// face info lookup allows searching for face info entry from elem/side pair
      54             :   std::map<side_type, unsigned int> _side_to_face_info;
      55             : 
      56             :   /// data structure storing the wave speeds SL, SM, SR
      57             :   std::map<side_type, std::vector<ADReal>> _wave_speed;
      58             : 
      59             :   ///@{ material properties computed by VarMat that Riemann solver needs
      60             :   const ADMaterialProperty<RealVectorValue> & _vel_elem;
      61             :   const ADMaterialProperty<RealVectorValue> & _vel_neighbor;
      62             :   const ADMaterialProperty<Real> & _speed_elem;
      63             :   const ADMaterialProperty<Real> & _speed_neighbor;
      64             :   const ADMaterialProperty<Real> & _pressure_elem;
      65             :   const ADMaterialProperty<Real> & _pressure_neighbor;
      66             :   const ADMaterialProperty<Real> & _rho_elem;
      67             :   const ADMaterialProperty<Real> & _rho_neighbor;
      68             :   const ADMaterialProperty<Real> & _specific_internal_energy_elem;
      69             :   const ADMaterialProperty<Real> & _specific_internal_energy_neighbor;
      70             :   const MaterialProperty<Real> * const _eps_elem;
      71             :   const MaterialProperty<Real> * const _eps_neighbor;
      72             :   ///@}
      73             : };