phase_change_data.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/parameter_handler.h>
 
namespace MeltPoolDG::Multiphase
{
  template <typename number>
  struct PhaseChangeData
  {
    struct LiquidGas
    {
      number boiling_temperature = 3133.0;
 
      number reference_pressure = 1.e5;
 
      number latent_heat_of_vaporization = 8.84e6;
    } liquid_gas;
 
    struct SolidLiquid
    {
      bool use_darcy_damping = false;
 
      number liquidus_temperature = 2200.;
 
      number solidus_temperature = 1933.;
    } solid_liquid;
 
    void
    add_parameters(dealii::ParameterHandler &prm)
    {
      prm.enter_subsection("phase change");
      {
        prm.enter_subsection("liquid gas");
        prm.add_parameter("boiling temperature",
                          liquid_gas.boiling_temperature,
                          "Boiling temperature at given reference pressure.",
                          dealii::Patterns::Double(0., std::numeric_limits<number>::max()));
        prm.add_parameter("reference pressure",
                          liquid_gas.reference_pressure,
                          "Reference pressure for boiling temperature.",
                          dealii::Patterns::Double(0., std::numeric_limits<number>::max()));
        prm.add_parameter("latent heat of vaporization",
                          liquid_gas.latent_heat_of_vaporization,
                          "Latent heat of vaporization (J/kg).",
                          dealii::Patterns::Double(0., std::numeric_limits<number>::max()));
        prm.leave_subsection();
        prm.enter_subsection("solid liquid");
        prm.add_parameter("use darcy damping",
                          solid_liquid.use_darcy_damping,
                          "Use Darcy damping?",
                          dealii::Patterns::Bool());
        prm.add_parameter("liquidus temperature",
                          solid_liquid.liquidus_temperature,
                          "Liquidus temperature.",
                          dealii::Patterns::Double(0., std::numeric_limits<number>::max()));
        prm.add_parameter("solidus temperature",
                          solid_liquid.solidus_temperature,
                          "Solidus temperature.",
                          dealii::Patterns::Double(0., std::numeric_limits<number>::max()));
        prm.leave_subsection();
      }
      prm.leave_subsection();
    }
  };
} // namespace MeltPoolDG::Multiphase