compressibleFlowFields.hpp

#ifndef ABLATELIBRARY_COMPRESSIBLEFLOWFIELDS_HPP
#define ABLATELIBRARY_COMPRESSIBLEFLOWFIELDS_HPP
 
#include <domain/region.hpp>
#include <memory>
#include <string>
#include <vector>
#include "domain/fieldDescriptor.hpp"
#include "eos/eos.hpp"
#include "parameters/mapParameters.hpp"
 
namespace ablate::finiteVolume {
 
class CompressibleFlowFields : public domain::FieldDescriptor {
   public:
    typedef enum { RHO, RHOE, RHOU, RHOV, RHOW } EulerComponents;
 
    inline const static std::string EULER_FIELD = "euler";
 
    inline const static std::string CONSERVED = "density";
 
    inline const static std::string YI_FIELD = eos::EOS::YI;
    inline const static std::string DENSITY_YI_FIELD = CONSERVED + YI_FIELD;
 
    inline const static std::string PROGRESS_FIELD = eos::EOS::PROGRESS;
    inline const static std::string DENSITY_PROGRESS_FIELD = CONSERVED + PROGRESS_FIELD;
 
    inline const static std::string EV_TAG = "EV";
 
    enum class ValidRange { POSITIVE, FULL, BOUND, MINUSONETOONE };
    inline const static std::string PositiveRange = "positive";
    inline const static std::string FullRange = "full";
    inline const static std::string BoundRange = "bound";
    inline const static std::string MinusOneToOneRange = "minusOneToOneRange";
 
    inline const static std::string EV_FIELD = "EV";
    inline const static std::string DENSITY_EV_FIELD = CONSERVED + EV_FIELD;
 
    inline const static std::string TEMPERATURE_FIELD = "temperature";
    inline const static std::string VELOCITY_FIELD = "velocity";
    inline const static std::string PRESSURE_FIELD = "pressure";
 
   protected:
    const std::shared_ptr<eos::EOS> eos;
    const std::shared_ptr<domain::Region> region;
    const std::shared_ptr<parameters::Parameters> conservedFieldOptions;
    const std::shared_ptr<parameters::Parameters> auxFieldOptions = ablate::parameters::MapParameters::Create({{"petscfv_type", "leastsquares"}, {"petsclimiter_type", "none"}});
 
   public:
    explicit CompressibleFlowFields(std::shared_ptr<eos::EOS> eos, std::shared_ptr<domain::Region> region = {}, std::shared_ptr<parameters::Parameters> conservedFieldParameters = {});
 
    std::vector<std::shared_ptr<domain::FieldDescription>> GetFields() override;
};
 
std::istream& operator>>(std::istream& is, CompressibleFlowFields::ValidRange& v);
 
}  // namespace ablate::finiteVolume
 
#endif  // ABLATELIBRARY_COMPRESSIBLEFLOWFIELDS_HPP