inertial.hpp

#ifndef ABLATELIBRARY_INERTIAL_HPP
#define ABLATELIBRARY_INERTIAL_HPP
 
#include <array>
#include "process.hpp"
 
namespace ablate::particles::processes {
 
class Inertial : public Process {
   private:
    const PetscReal fluidDensity;
 
    const PetscReal fluidViscosity;
 
    const std::array<PetscReal, 3> gravityField;
 
    const std::string eulerianVelocityField;
 
   public:
    explicit Inertial(const std::shared_ptr<parameters::Parameters>& parameters, const std::string& eulerianVelocityField = {});
 
    void ComputeRHS(PetscReal time, accessors::SwarmAccessor& swarmAccessor, accessors::RhsAccessor& rhsAccessor, accessors::EulerianAccessor& eulerianAccessor) override;
};
 
}  // namespace ablate::particles::processes
#endif  // ABLATELIBRARY_INERTIAL_HPP