MeltableInteraction.h

Go to the documentation of this file.

// This file is part of the MercuryDPM project (https://www.mercurydpm.org).
// Copyright (c), The MercuryDPM Developers Team. All rights reserved.
// License: BSD 3-Clause License; see the LICENSE file in the root directory.
 
#ifndef MERCURY_MELTABLEINTERACTION_H
#define MERCURY_MELTABLEINTERACTION_H
 
#include "Interactions/BaseInteraction.h"
#include "ParticleHandler.h"
#include "InteractionHandler.h"
 
class BaseParticle;
class MeltableParticle;
class BaseInteractable;
class MeltableNormalSpecies;
 
class MeltableInteraction : public virtual BaseInteraction {
public:
    MeltableInteraction(BaseInteractable *P, BaseInteractable *I, unsigned timeStamp)
            : BaseInteraction(P, I, timeStamp) {}
 
    MeltableInteraction(const MeltableInteraction &p) = default;
 
    MeltableInteraction() = default;
 
    void read(std::istream &is) override;
 
    void write(std::ostream &os) const override;
 
    void computeNormalForce();
 
    void actionsAfterTimeStep() override;
 
    Mdouble getElasticEnergy() const override;
 
    std::string getBaseName() const {
        return "Meltable";
    }
 
    const MeltableNormalSpecies* getMeltableNormalSpecies() const;
 
    Mdouble getBondingOverlap() const
    {
        return bondingOverlap_;
    }
 
    void setBondingOverlap(Mdouble bondingOverlap)
    {
        bondingOverlap_ = bondingOverlap;
    }
 
    void addBondingOverlap(Mdouble bondingOverlap)
    {
        bondingOverlap_ += bondingOverlap;
    }
 
    Mdouble getNeckRadius() const {
        return std::sqrt(neckRadius2_);
    };
 
    void addNeckRadius2(Mdouble diffNeckRadius2) {
        neckRadius2_ = std::max(neckRadius2_+diffNeckRadius2,0.0);
    }
 
    Mdouble getOverlapGrowthRate() const  {
        return Vec3D::dot(getI()->getVelocity()-getP()->getVelocity(),getNormal());
    }
 
private:
    Mdouble bondingOverlap_ = 0.0;
 
    Mdouble neckRadius2_ = constants::NaN;
};
 
#endif