MeltableParticle.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 MeltableParticle_H
#define MeltableParticle_H
 
#include "ThermalParticle.h"
#include "Species/NormalForceSpecies/MeltableNormalSpecies.h"
 
class MeltableParticle : public ThermalParticle {
public:
    MeltableParticle() = default;
 
    MeltableParticle(const MeltableParticle &p) = default;
 
    MeltableParticle *copy() const override;
 
    //void write(std::ostream &os) const override;
 
    std::string getName() const override;
 
    Mdouble getParticleProjectedArea();
 
    const MeltableNormalSpecies* getMeltableSpecies() const {
        return meltableSpecies_;
    }
 
    void actionsBeforeTimeStep() override;
 
    void actionsAfterTimeStep() override;
 
    unsigned getNumberOfFieldsVTK() const override { return 3; }
 
    std::string getTypeVTK(unsigned i) const override { return "Float32"; }
 
    std::string getNameVTK(unsigned i) const override;
 
    std::vector<Mdouble> getFieldVTK(unsigned i) const override;
 
    void setSpecies(const ParticleSpecies* species) override {
        BaseParticle::setSpecies(species);
        meltableSpecies_ = dynamic_cast<const MeltableNormalSpecies*>(species);
    }
 
    void addHeat(double heat) {
        heat_ += heat;
    }
 
    double getMeltRate(double solidRadius) const {
        return heat_/(4.0*constants::pi*solidRadius*solidRadius*getSpecies()->getDensity()*getMeltableSpecies()->getEffectiveLatentHeat());
    }
 
    double getSolidRadius() const {
        return meltableSpecies_->getRelativeSolidRadius(temperature_)*getRadius();
    }
 
    double getMoltenLayerThickness() const {
        return getRadius() - getSolidRadius();
    }
 
private:
    double heat_ = 0;
    const MeltableNormalSpecies* meltableSpecies_ = nullptr;
};
 
#endif