apparent_capacity.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <deal.II/base/vectorization.h>
 
#include <meltpooldg/core/material_data.hpp>
 
namespace MeltPoolDG::Heat
{
  template <typename number>
  struct QLQParameters
  {
    QLQParameters(const MaterialData<number> &material)
      : dT(material.liquidus_temperature - material.solidus_temperature)
      , T_peak(0.5 * (material.solidus_temperature + material.liquidus_temperature))
      , C_peak(2.0 * material.latent_heat_of_fusion / dT)
      , delta(dT / 6.0)
    {}
    const number dT;
 
    const number T_peak;
 
    const number C_peak;
 
    const number delta;
  };
 
  template <typename number>
  class ApparentCapacity
  {
  public:
    ApparentCapacity(const MaterialData<number> &material_data_in);
 
    dealii::VectorizedArray<number>
    evaluate(const dealii::VectorizedArray<number> &solid_fraction) const;
 
    number
    evaluate(const number &solid_fraction) const;
 
    dealii::VectorizedArray<number>
    compute_temperature_derivative(const dealii::VectorizedArray<number> &solid_fraction) const;
 
    number
    compute_temperature_derivative(const number &solid_fraction) const;
 
  private:
    const MaterialData<number> &material;
 
    const number inv_mushy_zone_radius;
 
    const number peak;
 
    const ApparentCapacityType apparent_capacity_type;
 
    const QLQParameters<number> qlq;
 
    number
    compute_temperature_derivative_qlq_from_solid_fraction(const number &solid_fraction) const;
 
    number
    evaluate_qlq_from_solid_fraction(const number &solid_fraction) const;
  };
 
} // namespace MeltPoolDG::Heat