sootTransportModel.hpp

#ifndef ABLATELIBRARY_SOOTTRANSPORTMODEL_HPP
#define ABLATELIBRARY_SOOTTRANSPORTMODEL_HPP
#include <Kokkos_Macros.hpp>
#ifndef KOKKOS_ENABLE_CUDA
 
#include <eos/eos.hpp>
#include <memory>
#include "eos/transport/transportModel.hpp"
 
namespace ablate::eos::tChemSoot {
 
class SootTransportModel : public eos::transport::TransportModel {
   private:
    const std::string fieldName;
    const std::shared_ptr<TransportModel> transport;
 
    constexpr static PetscReal solidCarbonFactor = 1.0 / 100.0;
 
    constexpr static std::size_t offset = 0;
 
    static PetscErrorCode VectorizeSpeciesDiffusionFunction(const PetscReal conserved[], PetscReal* property, void* ctx);
 
    static PetscErrorCode AdjustSpeciesDiffusionFunction(const PetscReal conserved[], PetscReal* property, void* ctx);
 
    static PetscErrorCode VectorizeSpeciesDiffusionTemperatureFunction(const PetscReal conserved[], PetscReal T, PetscReal* property, void* ctx);
 
    static PetscErrorCode AdjustSpeciesDiffusionTemperatureFunction(const PetscReal conserved[], PetscReal T, PetscReal* property, void* ctx);
 
   public:
    explicit SootTransportModel(const std::shared_ptr<TransportModel>& transport, std::string componentName);
 
    [[nodiscard]] ThermodynamicFunction GetTransportFunction(eos::transport::TransportProperty property, const std::vector<domain::Field>& fields) const override;
 
    [[nodiscard]] ThermodynamicTemperatureFunction GetTransportTemperatureFunction(eos::transport::TransportProperty property, const std::vector<domain::Field>& fields) const override;
};
}  // namespace ablate::eos::tChemSoot
 
#endif  // ABLATELIBRARY_SOOTTRANSPORTMODEL_HPP
#endif