chemistry.hpp

#ifndef ABLATELIBRARY_FINITEVOLUME_CHEMISTRY_HPP
#define ABLATELIBRARY_FINITEVOLUME_CHEMISTRY_HPP
 
#include <memory>
#include "eos/chemistryModel.hpp"
#include "process.hpp"
 
namespace ablate::finiteVolume::processes {
 
class Chemistry : public Process, public ablate::utilities::Loggable<Chemistry> {
   private:
    const std::shared_ptr<ablate::eos::ChemistryModel> chemistryModel;
 
    std::shared_ptr<ablate::eos::ChemistryModel::SourceCalculator> sourceCalculator;
 
    PetscErrorCode ChemistryPreStage(TS flowTs, ablate::solver::Solver &flow, PetscReal stagetime);
 
    static PetscErrorCode AddChemistrySourceToFlow(const FiniteVolumeSolver &solver, DM dm, PetscReal time, Vec locX, Vec fVec, void *ctx);
 
   public:
    explicit Chemistry(std::shared_ptr<ablate::eos::ChemistryModel>);
 
    void Setup(ablate::finiteVolume::FiniteVolumeSolver &flow) override;
 
    void Initialize(ablate::finiteVolume::FiniteVolumeSolver &flow) override;
 
    void AddChemistrySourceToFlow(const FiniteVolumeSolver &solver, Vec locX, Vec locFVec);
};
}  // namespace ablate::finiteVolume::processes
#endif