material_data.hpp

Go to the documentation of this file.

#pragma once
 
#include <meltpooldg/compressible_flow/eos_utils.hpp>
#include <meltpooldg/compressible_flow/material.hpp>
 
#include <memory>
 
namespace MeltPoolDG::CompressibleFlow
{
  template <int dim, typename number>
  class Material
  {
  public:
    const MaterialPhaseData<number> &data;
 
    std::shared_ptr<EOS::EquationOfStateUtils<dim, number>> eos_utils;
 
    explicit Material(const MaterialPhaseData<number> &material_data_in)
      : data(material_data_in)
    {
      eos_utils = make_eos_utils(data);
    }
 
  private:
    std::shared_ptr<EOS::EquationOfStateUtils<dim, number>>
    make_eos_utils(const MaterialPhaseData<number> &material_data)
    {
      switch (material_data.eos_type)
        {
          case EquationOfState::ideal_gas:
            return std::make_shared<EOS::IdealGas<dim, number>>(material_data);
          case EquationOfState::stiffened_gas:
            return std::make_shared<EOS::StiffenedGas<dim, number>>(material_data);
          case EquationOfState::noble_abel_stiffened_gas:
            return std::make_shared<EOS::NobleAbelStiffenedGas<dim, number>>(material_data);
          default:
            DEAL_II_NOT_IMPLEMENTED();
        }
    }
  };
} // namespace MeltPoolDG::CompressibleFlow