SemiAnalyticalIncompressibleFluidModel.h

Go to the documentation of this file.

#pragma once
#include "SemiAnalyticalIncompressibilityModule.h"
#include "Module/TopologyModule.h"
#include "Module/GroupModule.h"
 
#include "TriangularMeshConstraint.h"

 
namespace dyno
{
    template <typename TDataType>
    class PointSetToPointSet;
    template <typename TDataType>
    class ParticleIntegrator;
    template <typename TDataType>
    class NeighborPointQuery;
    template <typename TDataType>
    class IterativeDensitySolver;
    template<typename TDataType>
    class NeighborTriangleQuery;
    template <typename TDataType>
    class TriangularMeshConstraint;
    template <typename TDataType>
    class SurfaceTension;
    template <typename TDataType>
    class ImplicitViscosity;
    template <typename TDataType>
    class Helmholtz;
    template <typename>
    class PointSetToPointSet;
    typedef typename TopologyModule::Triangle Triangle;
 
    class ForceModule;
    class ConstraintModule;
    class Attribute;
 
    template <typename TDataType>
    class SemiAnalyticalIncompressibleFluidModel : public GroupModule
    {
        DECLARE_TCLASS(SemiAnalyticalIncompressibleFluidModel, TDataType)
    public:
        typedef typename TDataType::Real  Real;
        typedef typename TDataType::Coord Coord;
 
        SemiAnalyticalIncompressibleFluidModel();
 
    public:
        DEF_INSTANCE_IN(TriangleSet<TDataType>, TriangleSet, "");

        void updateImpl() override;

        void setSmoothingLength(Real len)
        {
            m_smoothing_length.setValue(len);
        }

        void setRestDensity(Real rho)
        {
            m_restRho = rho;
        }
 
    public:
        FVar<Real> m_smoothing_length;  //searching distance for particles
 
        FVar<Real> max_vel;               //no use
        FVar<Real> var_smoothing_length;  //no use
 
        DeviceArrayField<Real> m_particle_mass;  //mass of particles
 
        DeviceArrayField<Coord> m_particle_position;  //particle positions
        DeviceArrayField<Coord> m_particle_velocity;  //particle velocities
 
        DeviceArrayField<Attribute> m_particle_attribute;  //particle attributes, used to juedge if a particle is a fluid particle
 
        DeviceArrayField<Real>     m_triangle_vertex_mass;  //mass of triangle vertex
        DeviceArrayField<Coord>    m_triangle_vertex;       //current positions of triangle vertexs
        DeviceArrayField<Coord>    m_triangle_vertex_old;   //positions of triangle vertexs at last time step, used to update triangle velocities
        DeviceArrayField<Triangle> m_triangle_index;        //triangle vertex's indexes
 
        DeviceArrayField<Coord> m_particle_force_density;  //force density of fluid particles
        DeviceArrayField<Coord> m_vertex_force_density;    //no use
        DeviceArrayField<Coord> m_vn;                      //no use
 
        DeviceArrayField<int> m_flip;
        Reduction<Real>*      pReduce;
 
        DeviceArrayField<Coord> m_velocity_mod;  //velocity norm of each particle
 
    private:
        int  m_pNum;
        Real m_restRho;
        int  first = 1;
 
        //std::shared_ptr<ConstraintModule> m_surfaceTensionSolver;
        std::shared_ptr<ConstraintModule> m_viscositySolver;  // no use
 
        std::shared_ptr<ConstraintModule> m_incompressibilitySolver;  // no use
 
        std::shared_ptr<SemiAnalyticalIncompressibilityModule<TDataType>> m_pbdModule;  
 
        std::shared_ptr<TriangularMeshConstraint<TDataType>> m_meshCollision;  
 
        std::shared_ptr<ImplicitViscosity<TDataType>>  m_visModule;             
        std::shared_ptr<SurfaceTension<TDataType>>     m_surfaceTensionSolver;  
        std::shared_ptr<Helmholtz<TDataType>>          m_Helmholtz;             
        std::shared_ptr<PointSetToPointSet<TDataType>> m_mapping;               //no use
        std::shared_ptr<ParticleIntegrator<TDataType>> m_integrator;            
        std::shared_ptr<NeighborPointQuery<TDataType>>      m_nbrQueryPoint;         
 
        std::shared_ptr<NeighborTriangleQuery<TDataType>> m_nbrQueryTri;       
    };
}  // namespace dyno