arrheniusSublimation.hpp

#ifndef ABLATELIBRARY_ARRHENIUSSUBLIMATION_HPP
#define ABLATELIBRARY_ARRHENIUSSUBLIMATION_HPP
 
#include <map>
#include <memory>
#include "oneDimensionHeatTransfer.hpp"
#include "solver/cellSolver.hpp"
#include "solver/timeStepper.hpp"
#include "sublimationModel.hpp"
 
namespace ablate::boundarySolver::physics::subModels {
 
class ArrheniusSublimation : public SublimationModel {
   private:
    std::map<PetscInt, std::shared_ptr<OneDimensionHeatTransfer>> oneDimensionHeatTransfer;
 
    const std::shared_ptr<ablate::parameters::Parameters> properties;
 
    const std::shared_ptr<ablate::mathFunctions::MathFunction> initialization;
 
    const std::shared_ptr<ablate::parameters::Parameters> options;
 
    const PetscReal latentHeatOfFusion;
 
    const PetscReal solidDensity;
 
    const PetscReal preExponentialFactor;
    const PetscReal activationEnergy;
    const PetscReal parameterB;
 
    constexpr inline static PetscReal ugc = 8.31446261815324;  // J/K-mol
 
    [[nodiscard]] PetscReal ComputeMassFluxRate(PetscReal temperature) const;
 
   public:
    ArrheniusSublimation(const std::shared_ptr<ablate::parameters::Parameters> &properties, const std::shared_ptr<ablate::mathFunctions::MathFunction> &initialization,
                         const std::shared_ptr<ablate::parameters::Parameters> &options = {});
 
    void Initialize(ablate::boundarySolver::BoundarySolver &bSolver) override;
 
    bool RequiresUpdate() override { return true; };
 
    PetscErrorCode Update(PetscInt faceId, PetscReal dt, PetscReal heatFluxToSurface, PetscReal &temperature) override;
 
    PetscErrorCode Compute(PetscInt faceId, PetscReal heatFluxToSurface, SurfaceState &) override;
 
    [[nodiscard]] SerializerType Serialize() const override { return SerializerType::serial; }
 
    PetscErrorCode Save(PetscViewer viewer, PetscInt sequenceNumber, PetscReal time) override;
 
    PetscErrorCode Restore(PetscViewer viewer, PetscInt sequenceNumber, PetscReal time) override;
};
 
}  // namespace ablate::boundarySolver::physics::subModels
#endif  // ABLATELIBRARY_SUBLIMATIONTEMPERATURE_HPP