DislocationMobilityBase.h

Go to the documentation of this file.

/* This file is part of MoDELib, the Mechanics Of Defects Evolution Library.
 *
 *
 * MoDELib is distributed without any warranty under the
 * GNU General Public License (GPL) v2 <http://www.gnu.org/licenses/>.
 */
 
#ifndef _model_DislocationMobilityBase_h_
#define _model_DislocationMobilityBase_h_
 
#include <iostream>
#include <random>
#include <cmath>
#include <vector>
#include <chrono>
#include <assert.h>
#include <Eigen/Dense>
#include <Eigen/Geometry>
 // defines std::cout
#include <TerminalColors.h> // defines mode::cout
#include <PolycrystallineMaterialBase.h>
#include <StaticID.h>
#include <DefectiveCrystal/DefectiveCrystalTraitsIO.h>
 
namespace model
{
    
    struct StochasticForceGenerator
    {
        typedef std::default_random_engine T;
        std::vector<T> generators;
        std::vector<std::normal_distribution<double>> distributions;
        const int stochasticForceSeed;
        const int seed;
 
        StochasticForceGenerator(const DefectiveCrystalTraitsIO& traitsIO);
        double stochasticVelocity(const double& kB,
                               const double& T,
                               const double& B,
                               const double& L,
                               const double& dt);
        
    };
    
    struct DislocationMobilityBase : public StaticID<DislocationMobilityBase>
    {
        typedef Eigen::Matrix<double,3,3> MatrixDim;
        typedef Eigen::Matrix<double,3,1> VectorDim;
        
        static constexpr double kB_SI=1.38064852e-23; // Boltzmann constant in SI units [J/K]
        const std::string name;
 
        
        DislocationMobilityBase(const std::string& name_in) ;
        
        virtual ~DislocationMobilityBase(){};
        
        virtual double velocity(const MatrixDim& S,
                                const VectorDim& b,
                                const VectorDim& n,
                                const double& T,
                                const double& dL,
                                const double& dt,
                                const std::shared_ptr<StochasticForceGenerator>& sfg) const =0 ;
        
        double velocity(const MatrixDim& S,
                        const VectorDim& b,
                        const VectorDim& n,
                        const double& T) const ;
        
    };
    
}
#endif