constant.hpp

#ifndef ABLATELIBRARY_RADIATIONPROPERTIESCONSTANT_H
#define ABLATELIBRARY_RADIATIONPROPERTIESCONSTANT_H
 
#include "radiation/radiation.hpp"
#include "radiationProperties.hpp"
 
namespace ablate::eos::radiationProperties {
 
class Constant : public RadiationModel {
   private:
    struct FunctionContext {
        const PetscReal absorptivity;
        const PetscReal emissivity;
    };
 
    PetscReal absorptivityIn;
    PetscReal emissivityIn;
 
    static inline constexpr PetscReal GetRefractiveIndex() { return 1; }
 
    static PetscErrorCode ConstantAbsorptionTemperatureFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
 
    static PetscErrorCode ConstantEmissionTemperatureFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
 
   public:
    explicit Constant(double absorptivity, double emissivity);
    explicit Constant(const Constant&) = delete;
    void operator=(const Constant&) = delete;
 
    [[nodiscard]] ThermodynamicTemperatureFunction GetRadiationPropertiesTemperatureFunction(RadiationProperty property, const std::vector<domain::Field>& fields) const override;
};
 
}  // namespace ablate::eos::radiationProperties
 
#endif  // ABLATELIBRARY_CONSTANT_H