radiationFieldMonitor.hpp

#ifndef ABLATELIBRARY_RADIATIONFIELDMONITOR_H
#define ABLATELIBRARY_RADIATIONFIELDMONITOR_H
 
#include "eos/eos.hpp"
#include "eos/radiationProperties/radiationProperties.hpp"
#include "fieldMonitor.hpp"
#include "finiteVolume/compressibleFlowFields.hpp"
#include "io/interval/fixedInterval.hpp"
#include "utilities/constants.hpp"
 
namespace ablate::monitors {
 
class RadiationFieldMonitor : public FieldMonitor {
    enum FieldPlacements { intensity, absorption, fieldsStart };
    using tp = ablate::eos::ThermodynamicProperty;
    using fLoc = ablate::domain::FieldLocation;
 
   public:
    RadiationFieldMonitor(const std::shared_ptr<ablate::eos::EOS> eosIn, std::shared_ptr<eos::radiationProperties::RadiationModel> radiationModelIn,
                          std::shared_ptr<io::interval::Interval> intervalIn);
 
    void Register(std::shared_ptr<ablate::solver::Solver> solverIn) override;
    PetscErrorCode Save(PetscViewer viewer, PetscInt sequenceNumber, PetscReal time) override;
    static PetscErrorCode MonitorRadiation(TS ts, PetscInt step, PetscReal crtime, Vec u, void* ctx);
    PetscMonitorFunction GetPetscFunction() override { return MonitorRadiation; }
 
   private:
    const std::shared_ptr<ablate::eos::EOS> eos;  
 
    const std::shared_ptr<eos::radiationProperties::RadiationModel> radiationModel;  
 
    eos::ThermodynamicTemperatureFunction absorptivityFunction;  
 
    const std::shared_ptr<io::interval::Interval> interval;  
    PetscInt step;
 
    // Create suffix vector
    std::vector<std::string> fieldNames{"radiationIntensity", "absorptionCoefficient"};
};
}  // namespace ablate::monitors
#endif  // ABLATELIBRARY_RADIATIONFIELDMONITOR_H