MercuryDPM  Trunk
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
ParhamiMcMeekingSinterSpecies.h
Go to the documentation of this file.
1 //Copyright (c) 2013-2020, The MercuryDPM Developers Team. All rights reserved.
2 //For the list of developers, see <http://www.MercuryDPM.org/Team>.
3 //
4 //Redistribution and use in source and binary forms, with or without
5 //modification, are permitted provided that the following conditions are met:
6 // * Redistributions of source code must retain the above copyright
7 // notice, this list of conditions and the following disclaimer.
8 // * Redistributions in binary form must reproduce the above copyright
9 // notice, this list of conditions and the following disclaimer in the
10 // documentation and/or other materials provided with the distribution.
11 // * Neither the name MercuryDPM nor the
12 // names of its contributors may be used to endorse or promote products
13 // derived from this software without specific prior written permission.
14 //
15 //THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
16 //ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
17 //WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
18 //DISCLAIMED. IN NO EVENT SHALL THE MERCURYDPM DEVELOPERS TEAM BE LIABLE FOR ANY
19 //DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
20 //(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
21 //LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
22 //ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
23 //(INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
24 //SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
25 
26 #ifndef PARHAMIMCMEEKINGSINTERSPECIES_H
27 #define PARHAMIMCMEEKINGSINTERSPECIES_H
28 
30 #include "Math/ExtendedMath.h"
32 
38 {
39 public:
42 
45 
48 
51 
53  void read(std::istream& is);
54 
56  void write(std::ostream& os) const;
57 
59  std::string getBaseName() const;
60 
63 
64 //adhesion-specific functions
65 
66  //setters and getters
67  void set(Mdouble alpha, Mdouble beta, Mdouble atomicVolume /*Omega*/, Mdouble surfaceEnergy /*gamma_s*/,
68  Mdouble thicknessDiffusion /*deltaB*D0B*/, Mdouble activationEnergy /*QB*/, Mdouble temperature /*T*/,
69  Mdouble pseudoSlidingFrictionCoefficient /*\etaPart*/)
70  {
71  alpha_ = alpha;
72  beta_ = beta;
73  atomicVolume_ = atomicVolume;
74  surfaceEnergy_ = surfaceEnergy;
75  thicknessDiffusion_ = thicknessDiffusion;
76  activationEnergy_ = activationEnergy;
77  temperature_ = temperature;
78  pseudoSlidingFrictionCoefficient_ = pseudoSlidingFrictionCoefficient;
79 
80  Mdouble boltzmannConstant /*k_B*/ = 1.38064852e-23;
81  Mdouble gasConstant /*R_g*/ = 8.314459848;
82  Mdouble thicknessDiffusionVacancy /*DB*/ =
84  std::cout << thicknessDiffusionVacancy << "|" << thicknessDiffusion_ << std::endl;
85  Mdouble diffusionParameter /*DeltaB*/ =
86  atomicVolume_ / boltzmannConstant / temperature_ * thicknessDiffusionVacancy;
87  viscosityCoefficient_ = constants::pi / (2.0 * beta * diffusionParameter);
90  pseudoSlidingFrictionCoefficient_ * constants::pi / (2.0 * beta * diffusionParameter);
91  }
92 
94  { return viscosityCoefficient_; }
95 
97  { return adhesionCoefficient_; }
98 
100  { return slidingFrictionCoefficient_; }
101 
102 private:
103 
108  Mdouble surfaceEnergy_; /*gamma_s*/
109  Mdouble thicknessDiffusion_; /*deltaB*D0B*/
113 
116 
119 
122 };
123 
124 #endif
void mix(ParhamiMcMeekingSinterSpecies *S, ParhamiMcMeekingSinterSpecies *T)
creates default values for mixed species
ParhamiMcMeekingSinterInteraction InteractionType
The correct Interaction type for this AdhesiveForceSpecies.
double Mdouble
Definition: GeneralDefine.h:34
Mdouble exp(Mdouble Exponent)
Definition: ExtendedMath.cc:84
Mdouble beta(Mdouble z, Mdouble w)
This is the beta function, returns the approximation based on cmath's implementation of ln(gamma) ...
Mdouble viscosityCoefficient_
viscous force is viscosityCoefficient_*contactRadius^4*normalRelativeVelocity
Mdouble adhesionCoefficient_
adhesion force is adhesionCoefficient_*radius
ParhamiMcMeekingSinterSpecies()
The default constructor.
void set(Mdouble alpha, Mdouble beta, Mdouble atomicVolume, Mdouble surfaceEnergy, Mdouble thicknessDiffusion, Mdouble activationEnergy, Mdouble temperature, Mdouble pseudoSlidingFrictionCoefficient)
ParhamiMcMeekingSinterSpecies contains the parameters used to describe a linear reversible short-rang...
void read(std::istream &is)
Reads the species properties from an input stream.
void write(std::ostream &os) const
Writes the species properties to an output stream.
const Mdouble pi
Definition: ExtendedMath.h:45
Mdouble slidingFrictionCoefficient_
tangential force is slidingFrictionCoefficient_*contactRadius^2*radius*tangentialRelativeVelocity ...
~ParhamiMcMeekingSinterSpecies()
The default constructor.
Mdouble alpha_
viscous force is adhesionCoefficient_*temperature*contactRadius^4*normalRelativeVelocity ...
std::string getBaseName() const
Used in Species::getName to obtain a unique name for each Species.