radiationLoss.hpp

#ifndef ABLATELIBRARY_RADIATIONLOSS_HPP
#define ABLATELIBRARY_RADIATIONLOSS_HPP
 
#include "eos/radiationProperties/radiationProperties.hpp"
#include "finiteVolume/compressibleFlowFields.hpp"
#include "process.hpp"
#include "utilities/constants.hpp"
 
namespace ablate::finiteVolume::processes {
class RadiationLoss : public Process {
    const std::map<std::string, std::shared_ptr<ablate::mathFunctions::MathFunction>> functions;
 
    static PetscErrorCode ComputeRadiationLoss(PetscInt dim, PetscReal time, const PetscFVCellGeom* cg, const PetscInt uOff[], const PetscScalar u[], const PetscInt aOff[], const PetscScalar a[],
                                               PetscScalar f[], void* ctx);
 
   public:
    explicit RadiationLoss(std::shared_ptr<eos::radiationProperties::RadiationModel> radiationModelIn, double tInfinityIn = 300);
 
    void Setup(ablate::finiteVolume::FiniteVolumeSolver& fvmSolver) override;
 
    inline std::shared_ptr<eos::radiationProperties::RadiationModel> GetRadiationModel() { return radiationModel; }
 
    static inline PetscReal GetIntensity(PetscReal tInfinity, PetscReal temperature, PetscReal kappa) {
        tInfinity = PetscMax(tInfinity, (temperature - 500));
 
        // Compute the losses
        PetscReal netIntensity = -4.0 * ablate::utilities::Constants::sbc * temperature * temperature * temperature * temperature;
 
        netIntensity += 4.0 * ablate::utilities::Constants::sbc * tInfinity * tInfinity * tInfinity * tInfinity;
 
        // scale by kappa
        netIntensity *= kappa;
 
        return PetscAbsReal(netIntensity) > ablate::utilities::Constants::large ? ablate::utilities::Constants::large * PetscSignReal(netIntensity) : netIntensity;
    }
 
    const std::shared_ptr<eos::radiationProperties::RadiationModel> radiationModel;
 
    eos::ThermodynamicTemperatureFunction absorptivityFunction;
 
    PetscReal tInfinity;
};
 
}  // namespace ablate::finiteVolume::processes
 
#endif  // ABLATELIBRARY_RADIATIONLOSS_HPP