arbitrarySource.hpp

#ifndef ABLATELIBRARY_BOUNDARYARBITRARYSOURCE_HPP
#define ABLATELIBRARY_BOUNDARYARBITRARYSOURCE_HPP
 
#include "boundarySolver/boundaryProcess.hpp"
#include "eos/transport/transportModel.hpp"
#include "finiteVolume/processes/navierStokesTransport.hpp"
#include "finiteVolume/processes/pressureGradientScaling.hpp"
namespace ablate::boundarySolver::physics {
 
class ArbitrarySource : public BoundaryProcess {
   private:
    const std::map<std::string, std::shared_ptr<ablate::mathFunctions::MathFunction>> functions;
 
    const BoundarySolver::BoundarySourceType boundarySourceType;
 
    static PetscErrorCode ArbitrarySourceFunction(PetscInt dim, const ablate::boundarySolver::BoundarySolver::BoundaryFVFaceGeom *fg, const PetscFVCellGeom *boundaryCell, const PetscInt uOff[],
                                                  const PetscScalar *boundaryValues, const PetscScalar *stencilValues[], const PetscInt aOff[], const PetscScalar *auxValues,
                                                  const PetscScalar *stencilAuxValues[], PetscInt stencilSize, const PetscInt stencil[], const PetscScalar stencilWeights[], const PetscInt sOff[],
                                                  PetscScalar source[], void *ctx);
 
    PetscReal currentTime = 0.0;
 
   public:
    explicit ArbitrarySource(std::map<std::string, std::shared_ptr<ablate::mathFunctions::MathFunction>> functions, BoundarySolver::BoundarySourceType boundarySourceType);
 
    void Setup(ablate::boundarySolver::BoundarySolver &bSolver) override;
};
 
}  // namespace ablate::boundarySolver::physics
 
#endif  // ABLATELIBRARY_SUBLIMATION_HPP