buoyancy.hpp

#ifndef ABLATELIBRARY_BUOYANCY_HPP
#define ABLATELIBRARY_BUOYANCY_HPP
 
#include "flowProcess.hpp"
 
namespace ablate::finiteVolume::processes {
 
class Buoyancy : public FlowProcess {
   private:
    const std::vector<PetscReal> buoyancyVector;
    PetscReal densityAvg = NAN;
 
    PetscErrorCode UpdateAverageDensity(TS flowTs, ablate::solver::Solver& flow);
 
    static PetscErrorCode ComputeBuoyancySource(PetscInt dim, PetscReal time, const PetscFVCellGeom* cg, const PetscInt uOff[], const PetscScalar u[], const PetscInt aOff[], const PetscScalar a[],
                                                PetscScalar f[], void* ctx);
 
   public:
    explicit Buoyancy(std::vector<double> buoyancyVector);
    void Setup(ablate::finiteVolume::FiniteVolumeSolver& fv) override;
};
 
}  // namespace ablate::finiteVolume::processes
#endif  // ABLATELIBRARY_BUOYANCY_HPP