pressureGradientScaling.hpp

#ifndef ABLATELIBRARY_PRESSUREGRADIENTSCALING_HPP
#define ABLATELIBRARY_PRESSUREGRADIENTSCALING_HPP
 
#include <memory>
#include "eos/eos.hpp"
#include "finiteVolume/finiteVolumeSolver.hpp"
#include "flowProcess.hpp"
#include "monitors/logs/log.hpp"
 
namespace ablate::finiteVolume::processes {
class PressureGradientScaling : public FlowProcess, public io::Serializable {
   private:
    const std::string id = "PressureGradientScaling";
 
    const std::shared_ptr<eos::EOS> eos;
 
    const PetscReal maxMachAllowed = 0.7;
 
    const PetscReal maxAlphaChange = 0.10;
 
    const PetscReal maxAlphaAllowed;
 
    const PetscReal maxDeltaPressureFac;
 
    const PetscReal domainLength;
 
    const std::shared_ptr<ablate::monitors::logs::Log> log;
 
    PetscReal maxMach = 0.0;
    PetscReal alpha;
 
   public:
    PressureGradientScaling(std::shared_ptr<eos::EOS> eos, double alphaInit, double domainLength, double maxAlphaAllowed = {}, double maxDeltaPressureFac = {},
                            std::shared_ptr<ablate::monitors::logs::Log> = {});
 
    PetscErrorCode UpdatePreconditioner(TS flowTs, ablate::solver::Solver& flow);
 
    void Setup(ablate::finiteVolume::FiniteVolumeSolver& fv) override;
 
    // Alpha accessor
    inline const PetscReal& GetAlpha() { return alpha; }
 
    // Alpha accessor
    inline const PetscReal& GetMaxMach() { return maxMach; }
 
    const std::string& GetId() const override { return id; }
 
    PetscErrorCode Save(PetscViewer viewer, PetscInt sequenceNumber, PetscReal time) override;
 
    PetscErrorCode Restore(PetscViewer viewer, PetscInt sequenceNumber, PetscReal time) override;
};
}  // namespace ablate::finiteVolume::processes
#endif  // ABLATELIBRARY_PRESSUREGRADIENTSCALING_HPP