MercuryDPM  Trunk
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
ThermalParticle.cc
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 
27 #include "InteractionHandler.h"
31 #include "ParticleHandler.h"
32 #include "DPMBase.h"
33 
39 {
40  temperature_ = 0;
41  //temperatureDependentDensity_
42 }
43 
53  : BaseParticle(p)
54 {
57 }
58 
64 = default;
65 
72 {
73  return new ThermalParticle(*this);
74 }
75 
82 void ThermalParticle::write(std::ostream& os) const
83 {
85  os << " temperature " << temperature_;
86 }
87 
92 std::string ThermalParticle::getName() const
93 {
94  return "ThermalParticle";
95 }
96 
97 //todo Does mass and interaction radius change when a liquid film is added?
98 
107 void ThermalParticle::read(std::istream& is)
108 {
109  BaseParticle::read(is);
110  std::string dummy;
111  is >> dummy >> temperature_;
112 }
113 
115 {
116  return temperature_;
117 }
118 
120 {
121  temperature_ = temperature;
122 }
123 
125 {
126  temperature_ += temperature;
127 }
128 
130 {
132  {
133  temperature_ = timeDependentTemperature_(getHandler()->getDPMBase()->getTime());
134  }
136  {
138  radius_ = getRadius() * cbrt(getMass() / (getVolume() * density));
139  }
140 }
141 
142 const std::function<double(double)>& ThermalParticle::getTimeDependentTemperature() const
143 {
145 }
146 
147 void ThermalParticle::setTimeDependentTemperature(const std::function<double(double)>& timeDependentTemperature)
148 {
149  timeDependentTemperature_ = timeDependentTemperature;
150  temperature_ = timeDependentTemperature(0);
151  logger(INFO, "Setting initial temperature to %", temperature_);
152 }
void read(std::istream &is) override
Particle read function, which accepts an std::istream as input.
ThermalParticle * copy() const override
Particle copy method. It calls to copy constructor of this Particle, useful for polymorfism.
void actionsAfterTimeStep() override
std::string getName() const override
virtual Mdouble getVolume() const
Get Particle volume function, which required a reference to the Species vector. It returns the volume...
Logger< MERCURY_LOGLEVEL > logger("MercuryKernel")
double Mdouble
Definition: GeneralDefine.h:34
const std::function< double(double)> & getTemperatureDependentDensity() const
const ParticleSpecies * getSpecies() const
Returns a pointer to the species of this BaseInteractable.
~ThermalParticle() override
Particle destructor, needs to be implemented and checked if it removes tangential spring information...
Mdouble getTemperature() const
ParticleHandler * getHandler() const
Returns pointer to the particle's ParticleHandler.
ThermalParticle()
Basic Particle constructor, creates a particle at (0,0,0) with radius, mass and inertia equal to 1...
const std::function< double(double)> & getTemperatureDependentDensity() const
Mdouble getMass() const
Returns the particle's mass.
Definition: BaseParticle.h:322
void write(std::ostream &os) const override
Particle print function, which accepts an std::ostream as input.
void addTemperature(Mdouble temperature)
const std::function< double(double)> & getTimeDependentTemperature() const
Mdouble getRadius() const
Returns the particle's radius.
Definition: BaseParticle.h:345
std::function< double(double temperature)> timeDependentTemperature_
void setTemperature(Mdouble temperature)
Mdouble radius_
Definition: BaseParticle.h:652
void read(std::istream &is) override
void write(std::ostream &os) const override
void setTimeDependentTemperature(const std::function< double(double)> &timeDependentTemperature)