mps.hpp

Go to the documentation of this file.

#pragma once
 
#include "bucket.hpp"
#include "common.hpp"
#include "domain.hpp"
#include "input.hpp"
#include "neighbor_searcher.hpp"
#include "particles.hpp"
#include "pressure_calculator/interface.hpp"
#include "refvalues.hpp"
#include "settings.hpp"
#include "surface_detector/interface.hpp"
 
#include <Eigen/Dense>
#include <Eigen/Sparse>
#include <memory>
#include <vector>
 
class MPS {
public:
    Settings settings;                   
    Eigen::Vector3d gravity;             
    RefValues refValuesForNumberDensity; 
    RefValues refValuesForLaplacian;     
    RefValues refValuesForGradient;      
    Particles particles;                 
    Domain domain{};                     
 
    std::unique_ptr<PressureCalculator::Interface> pressureCalculator; 
 
    double courant{}; 
 
    MPS() = default;
 
    MPS(const Input& input,
        Eigen::Vector3d& gravity,
        std::unique_ptr<PressureCalculator::Interface>&& pressureCalculator,
        std::unique_ptr<SurfaceDetector::Interface>&& surfaceDetector);
 
    void stepForward();
 
private:
    NeighborSearcher neighborSearcher;                           
    std::unique_ptr<SurfaceDetector::Interface> surfaceDetector; 
 
    void calGravity();
 
    void calViscosity(const double& re);
 
    void moveParticle();
 
    void collision();
 
    void calNumberDensity(const double& re);
 
    void setBoundaryCondition();
 
    bool isFreeSurface(const Particle& pi);
 
    bool isParticleDistributionBiased(const Particle& pi);
 
    void setMinimumPressure(const double& re);
 
    void calPressureGradient(const double& re);
 
    void moveParticleUsingPressureGradient();
 
    void calCourant();
};