SinterNeck.hpp

Go to the documentation of this file.

// This file is part of GranOO, a workbench for DEM simulation.
//
// Author(s)    : - Damien Andre       IRCER/UNILIM, Limoges France
//                  <damien.andre@unilim.fr>
//                - Jean-luc Charles   Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <jean-luc.charles@ensam.eu>
//                - Jeremie Girardot   Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <jeremie.girardot@ensam.eu>
//                - Cedric Hubert      LAMIH/UPHF, Valenciennes France
//                  <cedric.hubert@uphf.fr>
//                - Ivan Iordanoff     Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <ivan.iordanoff@ensam.eu>
//
// Copyright (C) 2008-2019 D. Andre, JL. Charles, J. Girardot, C. Hubert, I. Iordanoff
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program.  If not, see <http://www.gnu.org/licenses/>.
 
 
#ifndef _libDEM_SinterNeck_hpp_
#define _libDEM_SinterNeck_hpp_
 
#include <string>
#include <vector>
 
 
 
 
#include "GranOO3/DEM/Bond.hpp"
#include "GranOO3/DEM/Element.hpp"
#include "GranOO3/Core/SetOf.hpp"
 
 
#include "GranOO3/Geom/Point.hpp"
 
namespace GranOO3
{
  namespace DEM
  {
 
 
    class SinterNeck : public Bond,
               public Core::Register< SinterNeck >,
               public Physic::CriticalTimeStep
    {
 
      GRANOO_CLASS(DEM::SinterNeck, Bond);
 
    public:
      //CONSTRUCTORS & DESTRUCTORS
      SinterNeck(Element& el1, Element& el2);
      virtual ~SinterNeck();
 
      //USEFULL
      
      virtual double get_linear_stiffness() const;
      virtual double get_angular_stiffness() const;
      virtual double compute_critical_time_step() const;
      virtual void draw();
      virtual std::string info() const;
 
      // COMMON METHOD
      void init_all()     ; // initialize all parameters
      void update_T()     ; // update temperature
      void update_n()     ; // update normal vector
      void update_r()     ; // update equivalent radius
      void update_h()     ; // update neck thickness 
      void update_vn()    ; // update normal relative velocity
      void update_a()     ; // update neck radius (only AFTER computing reaction forces !)
      virtual void update_delta_b() = 0; // update effective grain boundary diffusion coef 
      void force_neck_radius_to_zero(); // force a to be null
      
      //ACCESSORS
      GRANOO_ACCESS(T      , double, _T       ); // temperature
      GRANOO_ACCESS(a      , double, _a       ); // neck radius
      GRANOO_ACCESS(a0     , double, _a0      ); // initial neck radius
      GRANOO_ACCESS(h      , double, _h       ); // neck thickness
      GRANOO_ACCESS(h0     , double, _h0      ); // initial neck thickness
      GRANOO_ACCESS(r      , double, _r       ); // equivalent radius
      GRANOO_ACCESS(gamma_s, double, _gamma_s ); // surface energy
      GRANOO_ACCESS(omega  , double, _omega   ); // atomic volume 
      GRANOO_ACCESS(delta_b, double, _delta_b ); // effective grain boundary diffusion coef 
      GRANOO_ACCESS(vn     , double, _vn      ); // normal relative velocity
      GRANOO_ACCESS_GET(k  , double, k        ); // boltzmann constant
      GRANOO_ACCESS_GET(R  , double, R        ); // gas constant
      GRANOO_ACCESS(Fnv    , double, _Fnv     ); // normal viscous force (repulsive)
      GRANOO_ACCESS(Fna    , double, _Fna     ); // normal sintering force (attractive)
      GRANOO_ACCESS(n, Geom::Vector, _n       ); // normal direction
      
    private:
      SinterNeck() = delete;
      SinterNeck(const SinterNeck& frame) = delete;
      SinterNeck & operator=(const SinterNeck &) = delete;
 
    private:
      //SERIALIZATION
      friend class boost::serialization::access;
      template<class Archive> void save(Archive&, const unsigned int ) const;
      template<class Archive> void load(Archive&, const unsigned int);
      BOOST_SERIALIZATION_SPLIT_MEMBER();
 
    protected:
      double _T = 0.; // current temperature
      
      double _a  = 0.; // neck radius       (a_s: Martin,       a: Rojek)
      double _a0 = 0.; // initial neck radius (new)
      double _h  = 0.; // neck thickness    (h  : Martin,       δ: Rojek)
      double _r  = 0.; // equivalent radius (R* : Martin, \bar{r}: Rojek)
      double _h0 = 0.; // initial neck thickness (new)
      
      double _vn = 0.; // normal relative velocity     (dh/dt: Martin, v_rn: Rojek)
      
      double _gamma_s = 1.72    ; // surface energy (𝛾s: Martin,  𝛾s: Rojek)
      double _omega   = 1.18e-29; // atomic volume  ( Ω: Martin,    Ω: Rojek)
      
      double _delta_b = 0.; // effective grain boundary diffusion coef (Δb : Martin, Db: Rojek)
      
      const double k = 1.3806503e-23; // boltzmann constant
      const double R = 8.314;         // gas constant
      
      double _Fnv = 0.;   // normal viscous force (repulsive)
      double _Fna = 0.;   // normal sintering force (attractive)
 
      Geom::Vector _n = Geom::Vector(0.,0.,0.); // normal direction
 
    };
 
 
    template<class Archive> void
    SinterNeck::save(Archive& ar, const unsigned int v) const {
      ar << boost::serialization::base_object<Bond>(*this);
      ar << _T;
 
      ar << _a;
      if (v >= 1)
     ar << _a0;
 
      ar << _h;
      if (v >= 1)
     ar << _h0;
 
      ar << _r;
      ar << _gamma_s;
      ar << _omega;
      ar << _delta_b;
      ar << _vn;
      ar << _Fnv;
      ar << _Fna;
    }
 
    template<class Archive> void
    SinterNeck::load(Archive& ar, const unsigned int v) {
      ar >> boost::serialization::base_object<Bond>(*this);
      ar >> _T;
      
      ar >> _a;
      if (v >= 1)
    ar >> _a0;
 
      ar >> _h;
      if (v >= 1)
    ar >> _h0;
 
      ar >> _r;
      ar >> _gamma_s;
      ar >> _omega;
      ar >> _delta_b;
      ar >> _vn;
      ar >> _Fnv;
      ar >> _Fna;
    }
 
    inline void
    SinterNeck::draw() {
#ifndef SERVER
      Bond::draw();
#endif
    }
 
  }
}
 
 
#include <boost/serialization/version.hpp>
BOOST_CLASS_VERSION(GranOO3::DEM::SinterNeck, 1)
 
 
namespace GranOO3
{
  GRANOO_CLASS_DECLARE_TPL(DEM::SinterNeck);
}
 
 
#endif