MercuryDPM  Alpha
 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-2014, 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 #include "Species/BaseSpecies.h"
29 #include "Math/ExtendedMath.h"
31 
37 {
38 public:
41 
44 
47 
50 
52  void read(std::istream& is);
53 
55  void write(std::ostream& os) const;
56 
58  std::string getBaseName() const;
59 
62 
63 //adhesion-specific functions
64 
67 
68  //setters and getters
69  void set(Mdouble alpha, Mdouble beta, Mdouble atomicVolume /*Omega*/, Mdouble surfaceEnergy /*gamma_s*/,
70  Mdouble thicknessDiffusion /*deltaB*D0B*/, Mdouble activationEnergy /*QB*/, Mdouble temperature /*T*/,
71  Mdouble pseudoSlidingFrictionCoefficient /*\etaPart*/) {
72  alpha_ = alpha;
73  beta_ = beta;
74  atomicVolume_ = atomicVolume;
75  surfaceEnergy_ = surfaceEnergy;
76  thicknessDiffusion_ = thicknessDiffusion;
77  activationEnergy_ = activationEnergy;
78  temperature_ = temperature;
79  pseudoSlidingFrictionCoefficient_ = pseudoSlidingFrictionCoefficient;
80 
81  Mdouble boltzmannConstant /*k_B*/ = 1.38064852e-23;
82  Mdouble gasConstant /*R_g*/ = 8.314459848;
83  Mdouble thicknessDiffusionVacancy /*DB*/ = thicknessDiffusion_*exp(-activationEnergy_/gasConstant/temperature_);
84  std::cout << thicknessDiffusionVacancy << "|" << thicknessDiffusion_ << std::endl;
85  Mdouble diffusionParameter /*DeltaB*/ = atomicVolume_/boltzmannConstant/temperature_*thicknessDiffusionVacancy;
86  viscosityCoefficient_ = constants::pi/(2.0*beta*diffusionParameter);
89  }
90 
94 
95 private:
96 
101  Mdouble surfaceEnergy_; /*gamma_s*/
102  Mdouble thicknessDiffusion_; /*deltaB*D0B*/
106 
109 
112 
115 };
116 #endif
BaseSpecies is the class from which all other species are derived.
Definition: BaseSpecies.h:44
ParhamiMcMeekingSinterInteraction InteractionType
The correct Interaction type for this AdhesiveForceSpecies.
Mdouble exp(Mdouble Exponent)
Definition: ExtendedMath.cc:78
double Mdouble
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:42
void mix(ParhamiMcMeekingSinterSpecies *const S, ParhamiMcMeekingSinterSpecies *const T)
creates default values for mixed species
Mdouble getInteractionDistance() const
returns the largest separation distance at which adhesive short-range forces can occur.
Mdouble slidingFrictionCoefficient_
tangential force is slidingFrictionCoefficient_*contactRadius^2*radius*tangentialRelativeVelocity ...
virtual ~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.