Initialization.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.
 
#include "Mercury3D.h"
#include "cmath"
 
class Initialization : public Mercury3D {
protected:
    double radius = 50e-6;
    double density = 1050;
    double elasticModulus = 23.6e8;
    double poissonRatio = 0.4;
    double latentHeat = 56.4e3;// /(1-pow(1-relMoltenLayerThicknessMax,3));
    double solidHeatCapacity = 1200;
    double liquidHeatCapacity = 1.2*solidHeatCapacity;
    double temperatureInterval = 20;
    double meltingTemperature = 451.15;
    double thermalConductivity = 0.12;
    double thermalConvection = 150;
    double emmisivity = 0.9;
    double ambientTemperature = meltingTemperature-2.0*temperatureInterval;
    double surfaceTension = 34.4e-3;
    double activationEnergy = 32.5e3;
    double mass = density*4./3.*constants::pi*pow(radius,3);
    double dissipation = 0.9;
    double referenceViscosity = 1e-2;
    double gravity = -9.81;
    double timeStep = 19e-9;
 
    void setScaleFactor (double scaleFactor) {
        this->scaleFactor = scaleFactor;
        density *= scaleFactor;
        surfaceTension *= sqrt(scaleFactor); //0.035 * 1e-3 * sqrt(scaleMass);// * 2e-3;
        referenceViscosity *= sqrt(scaleFactor); //2e-4 * sqrt(scaleMass);
        gravity /= scaleFactor;
        solidHeatCapacity /= scaleFactor;
        liquidHeatCapacity /= scaleFactor;
        timeStep *= sqrt(scaleFactor);
    }
 
    double getScaleFactor() {
        return scaleFactor;
    }
 
private:
    double scaleFactor = 1.0;
};