HeatFluidCoupledParticle.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 HeatFluidCoupledParticle_H
#define HeatFluidCoupledParticle_H
 
#include "ThermalParticle.h"
#include "LiquidFilmParticle.h"
 
template<class Particle>
class HeatFluidCoupled : public Thermal<LiquidFilm<Particle>>
{
public:
    HeatFluidCoupled() = default;
 
    HeatFluidCoupled(const HeatFluidCoupled& p) = default;
 
    ~HeatFluidCoupled() override = default;
 
    HeatFluidCoupled* copy() const override
    {
        return new HeatFluidCoupled<Particle>(*this);
    }
 
    std::string getName() const override
    {
        return "HeatFluidCoupled" + Particle::getName();
    }
    unsigned getNumberOfFieldsVTK() const override { return 5; }
 
    std::string getTypeVTK(unsigned) const override { return "Float32"; }
 
    std::string getNameVTK(unsigned i) const override {
        if (i==0)
            return "fullLiquidVolume";
        else if (i==1)
            return "liquidFilmVolume";
        else if (i==2)
            return "liquidBridgeVolume";
        else if (i==3)
            return "totalEvaporatedLiquidVolume";
        else /* i=4 */
            return "temperature";
    }
 
    std::vector<Mdouble> getFieldVTK(unsigned i) const override {
        if (i==0)
            return { this->getFullLiquidVolume() };
        else if (i==1)
            return { this->liquidVolume_ };
        else if (i==2)
            return { this->getLiquidBridgeVolume() };
        else if (i==3)
            return { this->totalEvaporatedLiquidVolume_ };
        else /* i=4 */
            return { this->temperature_ };
    }
 
    void actionsAfterTimeStep() override {
        this->getSpecies()->actionsAfterTimeStep(this);
    }
};
 
typedef HeatFluidCoupled<SphericalParticle> HeatFluidCoupledParticle;
#endif