MercuryDPM  Trunk
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
ThermalInteraction.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 THERMALINTERACTION_H
27 #define THERMALINTERACTION_H
28 
31 
32 template<class NormalForceSpecies>
34 
35 template<class NormalForceInteraction>
37 {
38 public:
40 
44  ThermalInteraction(BaseInteractable* P, BaseInteractable* I, unsigned timeStamp);
45 
50 
55 
59  virtual ~ThermalInteraction();
60 
64  void computeNormalForce();
65 
66 };
67 
68 template<class NormalForceInteraction>
70  unsigned timeStamp)
71  : BaseInteraction(P, I, timeStamp), NormalForceInteraction(P, I, timeStamp)
72 {}
73 
74 template<class NormalForceInteraction>
77 {}
78 
79 template<class NormalForceInteraction>
82 {}
83 
84 template<class NormalForceInteraction>
86 {}
87 
88 template<class NormalForceInteraction>
90 {
91  NormalForceInteraction::computeNormalForce();
92  Mdouble radius = 2.0 * NormalForceInteraction::getEffectiveRadius();
93  Mdouble contactArea = constants::pi * radius * NormalForceInteraction::getOverlap();
94  const SpeciesType* species = static_cast<const SpeciesType*>(NormalForceInteraction::getBaseSpecies()->getNormalForce());
95  ThermalParticle* pParticle = dynamic_cast<ThermalParticle*>(NormalForceInteraction::getP());
96  ThermalParticle* iParticle = dynamic_cast<ThermalParticle*>(NormalForceInteraction::getI());
97  // if both p and i are particles
98  if (pParticle && iParticle)
99  {
100  /* heat transfer rate Q=m*c*dT/dt */
101  Mdouble heatTransfer = species->getThermalConductivity()
102  * (pParticle->getTemperature() - iParticle->getTemperature())
103  * contactArea / NormalForceInteraction::getDistance();
104  /* m*dT = Q/c*dt */
105  Mdouble mdT = heatTransfer / species->getHeatCapacity()
106  * NormalForceInteraction::getHandler()->getDPMBase()->getTimeStep();
107  pParticle->addTemperature(-mdT * pParticle->getInvMass());
108  iParticle->addTemperature(mdT * iParticle->getInvMass());
109  }
110 }
111 
112 #endif
void computeNormalForce()
Computes the normal forces due to linear plastic visco elastic interaction.
double Mdouble
Definition: GeneralDefine.h:34
Mdouble getTemperature() const
Mdouble getInvMass() const override
Returns the inverse of the particle's mass.
Definition: BaseParticle.h:293
Stores information about interactions between two interactable objects; often particles but could be ...
const Mdouble pi
Definition: ExtendedMath.h:45
void addTemperature(Mdouble temperature)
ThermalInteraction()
Default Constructor.
ThermalSpecies< typename NormalForceInteraction::SpeciesType > SpeciesType
Mdouble getHeatCapacity() const
Allows heatCapacity_ to be accessed.
Defines the basic properties that a interactable object can have.
virtual ~ThermalInteraction()
Destructor.
Mdouble getThermalConductivity() const
Allows heatCapacity_ to be accessed.