ParticleMixture.h

Go to the documentation of this file.

// SPDX-License-Identifier: MIT
// Copyright (c) 2021-2025 MPM-Geomechanics Development Team
 
#ifndef INC_PARTICLE_PARTICLEMIXTURE_H_
#define INC_PARTICLE_PARTICLEMIXTURE_H_
 
#include <Particle/Particle.h>
 
class ParticleMixture: public Particle {
 
public:
 
    ParticleMixture(const Vector3d& position, Material* material, const Vector3d& size);
 
    virtual ~ParticleMixture();
 
    virtual void updatePorosity();
 
    virtual inline double getMassFluid() const { return this->massFluid; }
 
    virtual inline const Vector3d* getVelocityFluid() const { return &(this->velocityFluid); }
 
    virtual inline const Vector3d* getExternalForceFluid() const { return &(this->externalForceFluid); }
 
    virtual Vector3d getDragForceFluid() const;
 
    virtual inline double getPressureFluid() const { return this->pressureFluid; }
 
    virtual inline double getPorosity() const { return this->porosityMixture; }
 
    virtual inline double getInitialPorosity() const { return this->initialPorosityMixture; }
 
    virtual inline double getDensityFluid() const { return this->densityFluid; }
 
    virtual inline void addExternalForceFluid(const Vector3d& delta_external_fluid_force) { this->externalForceFluid+=delta_external_fluid_force; }
 
    virtual inline void setVelocityFluid(const Vector3d& particle_velocity_fluid) { this->velocityFluid=particle_velocity_fluid;}
 
    virtual inline void setStrainIncrementFluid(const Matrix3d& strain_increment_fluid) { this->strainIncrementFluid=strain_increment_fluid; }
 
    virtual void setMaterial(Material* material);
    
    virtual inline const Matrix3d* getStrainIncrementFluid() const { return &(this->strainIncrementFluid); }
 
    virtual void updatePressure(double dt);
 
    virtual inline double getSaturation() const { return 1.0; }
 
    virtual inline void setPressureFluid(double pressure) { this->pressureFluid=pressure; }
 
    virtual double getCurrentVolume() const;
 
private:
 
    double massFluid; 
    double densityFluid; 
    double pressureFluid; 
    double porosityMixture; 
    double initialPorosityMixture; 
    
    Vector3d velocityFluid; 
    Vector3d externalForceFluid; 
    Vector3d internalForceFluid; 
    Matrix3d strainIncrementFluid; 
};
 
#endif /* INC_PARTICLE_PARTICLEMIXTURE_H_ */