surfaceForce.hpp

#ifndef ABLATELIBRARY_FINITEVOLUME_CHEMISTRY_HPP
#define ABLATELIBRARY_FINITEVOLUME_CHEMISTRY_HPP
 
#include <petsc.h>
#include <memory>
#include <vector>
#include "domain/range.hpp"
#include "finiteVolume/fluxCalculator/fluxCalculator.hpp"
#include "flowProcess.hpp"
#include "process.hpp"
#include "solver/solver.hpp"
#include "twoPhaseEulerAdvection.hpp"
 
namespace ablate::finiteVolume::processes {
 
class SurfaceForce : public Process {
    PetscReal sigma;
 
   private:
    struct VertexStencil {
        std::vector<PetscInt> stencil;
        PetscInt stencilSize;
        PetscInt vertexId;
        std::vector<PetscScalar> gradientWeights;
        std::vector<PetscScalar> stencilCoord;
    };
    DM dmData;
 
   public:
    // Hold a list of VortexStencils
    std::vector<VertexStencil> vertexStencils;
 
    explicit SurfaceForce(PetscReal sigma);
 
    ~SurfaceForce() override;
 
    void Setup(ablate::finiteVolume::FiniteVolumeSolver &flow) override;
    static PetscErrorCode ComputeSource(const FiniteVolumeSolver &solver, DM dm, PetscReal time, Vec locX, Vec locFVec, void *ctx);
};
}  // namespace ablate::finiteVolume::processes
#endif