sublimationModel.hpp

#ifndef ABLATELIBRARY_SUBLIMATIONMODEL_HPP
#define ABLATELIBRARY_SUBLIMATIONMODEL_HPP
 
#include <petsc.h>
#include "boundarySolver/boundarySolver.hpp"
 
namespace ablate::boundarySolver::physics::subModels {
 
class SublimationModel : public io::Serializable {
   private:
    // static name of this model
    inline const static std::string sublimationModelId = "SublimationModel";
 
   public:
    struct SurfaceState {
        PetscReal massFlux;
 
        PetscReal temperature;
 
        PetscReal regressionRate;
    };
 
    virtual bool RequiresUpdate() { return false; };
 
    virtual void Initialize(ablate::boundarySolver::BoundarySolver &bSolver){};
 
    virtual PetscErrorCode Update(PetscInt faceId, PetscReal dt, PetscReal heatFluxToSurface, PetscReal &temperature) { return PETSC_SUCCESS; };
 
    virtual PetscErrorCode Compute(PetscInt faceId, PetscReal heatFluxToSurface, SurfaceState &) = 0;
 
    ~SublimationModel() override = default;
 
    [[nodiscard]] SerializerType Serialize() const override { return io::Serializable::SerializerType::none; }
 
    [[nodiscard]] const std::string &GetId() const override { return sublimationModelId; }
 
    PetscErrorCode Save(PetscViewer viewer, PetscInt sequenceNumber, PetscReal time) override { return PETSC_SUCCESS; };
 
    PetscErrorCode Restore(PetscViewer viewer, PetscInt sequenceNumber, PetscReal time) override { return PETSC_SUCCESS; }
};
}  // namespace ablate::boundarySolver::physics::subModels
 
#endif  // ABLATELIBRARY_SUBLIMATIONMODEL_HPP