volumeRadiation.hpp

#ifndef ABLATELIBRARY_VOLUMERADIATION_HPP
#define ABLATELIBRARY_VOLUMERADIATION_HPP
 
#include "domain/dynamicRange.hpp"
#include "io/interval/interval.hpp"
#include "radiation.hpp"
 
namespace ablate::radiation {
 
class VolumeRadiation : public solver::CellSolver, public solver::RHSFunction {
   public:
    PetscErrorCode ComputeRHSFunction(PetscReal time, Vec locXVec, Vec locFVec) override;
 
    void Initialize() override;
    void Setup() override;
    void Register(std::shared_ptr<ablate::domain::SubDomain> subDomain) override;
 
    VolumeRadiation(const std::string& solverId1, const std::shared_ptr<io::interval::Interval>& interval, std::shared_ptr<radiation::Radiation> radiation,
                    const std::shared_ptr<parameters::Parameters>& options1, const std::shared_ptr<monitors::logs::Log>& unnamed1);
 
    ~VolumeRadiation() override;
 
    PetscErrorCode PreRHSFunction(TS ts, PetscReal time, bool initialStage, Vec locX) override;
 
   private:
    const std::shared_ptr<io::interval::Interval> interval;
    std::shared_ptr<ablate::radiation::Radiation> radiation;
    ablate::domain::DynamicRange radiationCellRange;
 
    eos::ThermodynamicTemperatureFunction absorptivityFunction;
    eos::ThermodynamicTemperatureFunction emissivityFunction;
};
}  // namespace ablate::radiation
#endif  // ABLATELIBRARY_VOLUMERADIATION_HPP