contactmodelthermal.h

#pragma once
// contactmodelthermal.h
 
#include "contactmodel.h"
#include "contactmodel/interface/icontactmodelthermal.h"
 
#ifndef __GNUC__
#pragma warning(push)
#pragma warning(disable:4251)
#endif
namespace itascaxd
{
 class IContactThermal;
 class IContactMechanical;
}
 
namespace cmodelsxd
{
 using namespace itasca;
 using namespace itascaxd;
 
 class ContactModelThermalState : public ContactModelState
 {
 public:
 ContactModelThermalState() : power_(0.0), length_(0.0),
 end1Temperature_(0.0), end2Temperature_(0.0), tempInc_(0.0), gapInc_(0.0)
 { };
 
 ~ContactModelThermalState() override { };
 
 virtual const IContactThermal *getThermalContact() const=0;
 virtual const IContactMechanical *getMechanicalContact() const=0;
 const IProgram * getProgram() const override =0;
 
 double power_; // current total power (from 1 to 2) at this contact (computed by the contact model)
 double length_; // length of the pipe
 double end1Temperature_; // temperature in reservoir 1
 double end2Temperature_; // temperature in reservoir 2
 double tempInc_; // temperature increment in contact
 double gapInc_; // gap increment in contact since last coupling step
 };
 
 class CONTACTMODEL_EXPORT ContactModelThermal : public IContactModelThermal
 ,public ContactModel
 {
 public:
 ContactModelThermal();
 ~ContactModelThermal() override;
 
 static const char * getPluginPrefix() { return "contactmodelthermal"; }
 static const char * getPluginDirectory() { return "contactmodels/thermal"; }
 
 // IContactModelThermal functions.
 
 IContactModel * getContactModel() override {return this;}
 const IContactModel * getContactModel() const override {return this;}
 
 // Returns true if the property was actually updated in the contact model.
 virtual bool endPropertyUpdated(const QString &name,const IContactThermal *c)=0;
 
 // For the contact model state
 virtual bool validate(ContactModelThermalState *state,const double &timestep)=0;
 virtual bool checkActivity(ContactModelThermalState *d) const=0;
 
 // Returns Activity Distance (defaults to 0.0) -- distance between objects that contact is considered active.
 double getActivityDistance() const override {return 0.0;}
 
 // Updates power from end1 to end2.
 // Returns true if contact is active
 virtual bool updatePower(ContactModelThermalState *state,const double &timestep)=0;
 
 // Returns timestep limits (defaults to (0.0,infinity))
 // These are further timestep restrictions that the contact model can place on the simulation, over and above those implied by stiffness.
 virtual DVect2 getTimestepLimits(ContactModelThermalState *,const double &) const { return DVect2(0.0,limits<double>::max());}
 
 // Returns the Resistance per unit length
 virtual double getEffectiveResistance() const {return 0.0;}
 
 // For contact specific plotting
 void getSphereList(const IContact *,std::vector<DVect> *,std::vector<double> *,
 std::vector<double> *) override { }
#ifdef THREED
 void getDiskList(const IContact *,std::vector<DVect> *,std::vector<DVect> *,
 std::vector<double> *,std::vector<double> *) override { }
#endif
 void getCylinderList(const IContact *,std::vector<DVect> *,std::vector<DVect> *,
 std::vector<double> *,std::vector<double> *,
 std::vector<double> *) override { }
 };
} // namespace cmodel
 
#ifndef __GNUC__
#pragma warning(pop)
#endif
 
// EoF