b2heat_material.H

//------------------------------------------------------------------------
// b2heat_material.H --
//
//
// written by Mathias Doreille
//            Neda Ebrahimi Pour <neda.ebrahimipour@dlr.de>
//
// (c) 2004-2016 SMR Engineering & Development SA
//               2502 Bienne, Switzerland
//
// (c) 2023 Deutsches Zentrum für Luft- und Raumfahrt (DLR) e.V.
//          Linder Höhe, 51147 Köln
//
// All Rights Reserved.  Proprietary source code.  The contents of
// this file may not be disclosed to third parties, copied or
// duplicated in any form, in whole or in part, without the prior
// written permission of SMR.
//------------------------------------------------------------------------
 
#ifndef B2HEAT_MATERIAL_H_
#define B2HEAT_MATERIAL_H_
 
#include "model/b2element.H"
#include "utils/b2linear_algebra.H"
 
namespace b2000 {
 
class HeatMaterial : virtual public ElementProperty {
public:
    virtual int number_of_layer() const { return 1; }
 
    virtual void laminate(
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          double& laminate_begin, double& laminate_end) {
        laminate_begin = -1;
        laminate_end = 1;
    }
 
    virtual void layer(
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          const int layer_id, double& layer_begin, double& layer_end) {
        layer_begin = -1;
        layer_end = 1;
    }
 
    enum LinearType { yes, linear_on_sub_domain, nonlinear };
 
    virtual LinearType linear(const int layer_id = -1) const = 0;
 
    virtual bool isotropic(const int layer_id = -1) const = 0;
 
    typedef ObjectTypeIncomplete<HeatMaterial, ElementProperty::type_t> type_t;
    static type_t type;
};
 
class HeatMaterial3D : public HeatMaterial {
public:
    virtual void get_conduction_and_capacity_heat(
          Model* model, const Element* element, const double coordinates[3],
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          const double element_coordinate[3], const int layer_id, const double covariant_base[3][3],
          const double temperature, const double temperature_dot, const double grad_temperature[3],
          const double time, const double linear, double heat_conduction[3],
          double d_heat_conduction_d_temperature[3], double d_heat_conduction_d_grad_temperature[6],
          double& heat_capacity, double& d_heat_capacity_d_temperature,
          double& d_heat_capacity_d_temperature_dot, GradientContainer* gradient_container,
          SolverHints* solver_hints) = 0;
 
    virtual void get_radiation_and_convection_heat(
          Model* model, const Element* element, const double coordinates[3],
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          const double element_coordinate[], const double temperature, const double time,
          const bool linear, double& heat, double& d_heat_d_temperature,
          GradientContainer* gradient_container, SolverHints* solver_hints) = 0;
 
    virtual bool get_radiation_properties(
          Model* model, const Element* element, const double coordinates[3],
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          const double element_coordinate[], const double temperature, double& heat_radiation_unit,
          double& receptivity) = 0;
};
 
class HeatMaterial2D : public HeatMaterial {
public:
    virtual double get_thickness(
          const Element* element,
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation) const = 0;
 
    virtual void get_conduction_and_capacity_heat(
          Model* model, const Element* element, const double coordinates[3],
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          const double element_coordinate[2], const int layer_id, const double covariant_base[2][2],
          const double temperature, const double temperature_dot, const double grad_temperature[2],
          const double time, const double linear, double heat_conduction[2],
          double d_heat_conduction_d_temperature[2], double d_heat_conduction_d_grad_temperature[3],
          double& heat_capacity, double& d_heat_capacity_d_temperature,
          double& d_heat_capacity_d_temperature_dot, GradientContainer* gradient_container,
          SolverHints* solver_hints) = 0;
 
    virtual void get_radiation_and_convection_heat(
          Model* model, const Element* element, const double coordinates[3],
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          const double element_coordinate[], const double temperature, const double time,
          const bool linear, double& heat, double& d_heat_d_temperature,
          GradientContainer* gradient_container, SolverHints* solver_hints) = 0;
 
    virtual bool get_radiation_properties(
          Model* model, const Element* element, const double coordinates[3],
          const b2linalg::Vector<double, b2linalg::Vdense_constref> nodes_interpolation,
          const double element_coordinate[], const double temperature, double& heat_radiation_unit,
          double& receptivity) = 0;
};
 
}  // namespace b2000
 
#endif  // B2HEAT_MATERIAL_H_