fluid_particles.h

Go to the documentation of this file.

/* -------------------------------------------------------------------------*
 *                              SPHinXsys                                   *
 * --------------------------------------------------------------------------*
 * SPHinXsys (pronunciation: s'finksis) is an acronym from Smoothed Particle    *
 * Hydrodynamics for industrial compleX systems. It provides C++ APIs for   *
 * physical accurate simulation and aims to model coupled industrial dynamic *
 * systems including fluid, solid, multi-body dynamics and beyond with SPH  *
 * (smoothed particle hydrodynamics), a meshless computational method using *
 * particle discretization.                                                 *
 *                                                                          *
 * SPHinXsys is partially funded by German Research Foundation              *
 * (Deutsche Forschungsgemeinschaft) DFG HU1527/6-1, HU1527/10-1                *
 * and HU1527/12-1.                                                         *
 *                                                                           *
 * Portions copyright (c) 2017-2020 Technical University of Munich and      *
 * the authors' affiliations.                                               *
 *                                                                           *
 * Licensed under the Apache License, Version 2.0 (the "License"); you may   *
 * not use this file except in compliance with the License. You may obtain a *
 * copy of the License at http://www.apache.org/licenses/LICENSE-2.0.        *
 *                                                                           *
 * --------------------------------------------------------------------------*/
#ifndef FLUID_PARTICLES_H
#define FLUID_PARTICLES_H

#include "base_particles.hpp"
#include "fluid_particles_variable.h"

#include "particle_generator_lattice.h"

namespace SPH
{
    class Fluid;
    class Oldroyd_B_Fluid;
    class CompressibleFluid;

    class FluidParticles : public BaseParticles
    {
    public:
        StdLargeVec<Real> p_;                
        StdLargeVec<Real> drho_dt_;          
        StdLargeVec<Real> rho_sum_;          
        StdLargeVec<int> surface_indicator_; 
        FluidParticles(SPHBody &sph_body, Fluid *fluid);
        virtual ~FluidParticles(){};

        virtual void initializeOtherVariables() override;
        virtual FluidParticles *ThisObjectPtr() override { return this; };
    };

    class ViscoelasticFluidParticles : public FluidParticles
    {
    public:
        StdLargeVec<Matd> tau_;     
        StdLargeVec<Matd> dtau_dt_; 
        ViscoelasticFluidParticles(SPHBody &sph_body, Oldroyd_B_Fluid *oldroyd_b_fluid);
        virtual ~ViscoelasticFluidParticles(){};

        virtual void initializeOtherVariables() override;
        virtual ViscoelasticFluidParticles *ThisObjectPtr() override { return this; };
    };

    class CompressibleFluidParticles : public FluidParticles
    {
    public:
        StdLargeVec<Vecd> mom_;     
        StdLargeVec<Vecd> dmom_dt_; 
        StdLargeVec<Vecd> dmom_dt_prior_;
        StdLargeVec<Real> E_;     
        StdLargeVec<Real> dE_dt_; 
        StdLargeVec<Real> dE_dt_prior_;

        CompressibleFluidParticles(SPHBody &sph_body, CompressibleFluid *compressible_fluid);
        virtual ~CompressibleFluidParticles(){};

        virtual void initializeOtherVariables() override;
        virtual CompressibleFluidParticles *ThisObjectPtr() override { return this; };
    };

    class WeaklyCompressibleFluidParticles : public FluidParticles
    {
    public:
        StdLargeVec<Real> dmass_dt_;      
        StdLargeVec<Vecd> mom_;           
        StdLargeVec<Vecd> dmom_dt_;       
        StdLargeVec<Vecd> dmom_dt_prior_; 
        WeaklyCompressibleFluidParticles(SPHBody &sph_body, Fluid *fluid);
        virtual ~WeaklyCompressibleFluidParticles(){};

        virtual void initializeOtherVariables() override;
        virtual WeaklyCompressibleFluidParticles *ThisObjectPtr() override { return this; };
    };
}
#endif // FLUID_PARTICLES_H