radiative_transport_data.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/parameter_handler.h>
 
#include <meltpooldg/core/finite_element_data.hpp>
#include <meltpooldg/linear_algebra/linear_solver_data.hpp>
#include <meltpooldg/radiative_transport/pseudo_rte_data.hpp>
#include <meltpooldg/utilities/enum.hpp>
 
 
namespace MeltPoolDG::RadiativeTransport
{
  // Choose the formulation of the absorptivity coefficient
  BETTER_ENUM(AbsorptivityType,
              char,
              // absorptivity coefficient as material constants
              constant,
              // absorptivity coefficient as a function of the heaviside gradient, and applying a
              // MacAulay bracket
              gradient_based)
 
  BETTER_ENUM(RTEPredictorType,
              char,
              // no predictor specified; use old value as initial guess
              none,
              // use pseude time stepping operation as a predictor
              pseudo_time_stepping)
 
 
  template <typename number>
  struct RadiativeTransportData
  {
    RadiativeTransportData();
 
    FiniteElementData        fe;
    LinearSolverData<number> linear_solver;
    int                      verbosity_level                    = -1;
    RTEPredictorType         predictor_type                     = RTEPredictorType::none;
    AbsorptivityType         absorptivity_type                  = AbsorptivityType::gradient_based;
    number                   avoid_singular_matrix_absorptivity = 1e-16;
 
    struct AbsorptivityConstantData
    {
      number absorptivity_gas    = 0.1;
      number absorptivity_liquid = 0.9;
 
      void
      add_parameters(dealii::ParameterHandler &prm);
 
    } absorptivity_constant_data;
 
    struct AbsorptivityGradientBasedData
    {
      number avoid_div_zero_constant = 1e-16;
 
      void
      add_parameters(dealii::ParameterHandler &prm);
 
    } absorptivity_gradient_based_data;
 
    PseudoTimeSteppingData<number> pseudo_time_stepping;
 
    void
    add_parameters(dealii::ParameterHandler &prm);
 
    void
    post(const FiniteElementData &base_fe_data, const unsigned int base_verbosity_level);
 
    void
    check_input_parameters(const FiniteElementData &base_fe_data) const;
  };
} // namespace MeltPoolDG::RadiativeTransport