contactmodelthermal.h

#pragma once
// contactmodelthermal.h
 
#include "contactmodel.h"
#include "contactmodel/interface/icontactmodelthermal.h"
 
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 string &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
 
// EoF