material.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/parameter_handler.h>
 
#include <meltpooldg/utilities/enum.hpp>
 
#include <array>
#include <cstddef>
#include <limits>
#include <memory>
 
namespace MeltPoolDG::CompressibleFlow
{
  BETTER_ENUM(EquationOfState,
              char,
              ideal_gas,
              stiffened_gas,
              noble_abel_stiffened_gas,
              weakly_compressible)
 
  
  template <typename number>
  struct EOSData
  {
    number stiffening_pressure = std::numeric_limits<number>::max();
 
    number covolume = std::numeric_limits<number>::max();
 
    number heat_bound = std::numeric_limits<number>::min();
 
    struct WeaklyCompressibleSpecificData
    {
      number reference_state_pressure = std::numeric_limits<number>::max();
 
      number reference_state_density = std::numeric_limits<number>::max();
 
      number reference_state_sound_speed = std::numeric_limits<number>::max();
    } weakly_compressible;
 
    void
    add_parameters(dealii::ParameterHandler &prm);
 
    void
    post(const EquationOfState eos_type) const;
  };
 
  template <typename number>
  struct MaterialSpeciesData
  {
    number specific_isobaric_heat = 1000.0;
 
    number specific_isochoric_heat = 1000.0;
 
    number dynamic_viscosity = 1. / 1600.;
 
    number gamma = 1.4;
 
    number specific_gas_constant = 287.1;
 
    number thermal_conductivity = std::numeric_limits<number>::max();
 
    EOSData<number> eos_data;
 
    number molar_mass = std::numeric_limits<number>::max();
 
    // Name of the species (for output purposes)
    std::string name = "species";
 
    void
    add_parameters(dealii::ParameterHandler &prm, const unsigned species_number);
 
  private:
    void
    add_parameters(dealii::ParameterHandler &prm);
  };
 
  template <typename number>
  struct MaterialPhaseData
  {
    constexpr static std::size_t n_max_species = 2;
 
    constexpr static std::size_t n_max_diffusion_coefficients =
      n_max_species * (n_max_species - 1) / 2;
 
    EquationOfState eos_type = EquationOfState::ideal_gas;
 
    unsigned int number_of_species = 1;
 
    std::array<MaterialSpeciesData<number>, n_max_species> species_data;
 
    number reference_density = 1.0;
 
    number reference_dynamic_viscosity = 1. / 1600.;
 
    number          specific_isobaric_heat;
    number          specific_isochoric_heat;
    number          dynamic_viscosity;
    number          gamma;
    number          specific_gas_constant;
    number          thermal_conductivity;
    EOSData<number> eos_data;
 
    number
    get_diffusion_coefficient(const unsigned int i, const unsigned int j) const;
 
    void
    add_parameters(dealii::ParameterHandler &prm, const bool is_gas_phase);
 
    void
    post(const bool is_gas);
 
  private:
    struct SpeciesInteraction
    {
      unsigned int species_i             = 1;
      unsigned int species_j             = 2;
      number       diffusion_coefficient = 0.0;
    };
 
    std::array<SpeciesInteraction, n_max_diffusion_coefficients> species_interactions_input{};
 
    std::array<number, n_max_diffusion_coefficients> species_interactions{};
 
    std::size_t
    get_species_interaction_index(const unsigned int i, const unsigned int j) const;
  };
} // namespace MeltPoolDG::CompressibleFlow