contactmodelfluid.h

#pragma once
// contactmodelfluid.h
 
#include "contactmodel/interface/icontactmodelfluid.h"
#include "contactmodel.h"
 
#ifndef __GNUC__
#pragma warning(push)
#pragma warning(disable : 4251)
#endif
namespace itascaxd {
 class IContactFluid;
 class IContactMechanical;
}
 
namespace cmodelsxd {
 using namespace itasca;
 using namespace itascaxd;
 
 class ContactModelFluidState : public ContactModelState { //DPfuture: remove?
 public:
 ContactModelFluidState()
 : gapInc_(0.0), dp_geom_(0.0) {};
 
 ~ContactModelFluidState() override{};
 
 virtual const IContactFluid* getFluidContact() const = 0;
 virtual const IContactMechanical* getMechanicalContact() const = 0;
 const IProgram* getProgram() const override = 0;
 
 double gapInc_; // gap increment in contact since last coupling step
 double dp_geom_; // strain-induced pressure increment
 };
 
 class CONTACTMODEL_EXPORT ContactModelFluid : public IContactModelFluid,
 public ContactModel {
 public:
 ContactModelFluid();
 ~ContactModelFluid() override;
 
 static const char* getPluginPrefix() { return "contactmodelfluid"; }
 static const char* getPluginDirectory() { return "contactmodels/fluid"; }
 IContactModel* getContactModel() override { return this; }
 const IContactModel* getContactModel() const override { return this; }
 
 virtual bool endPropertyUpdated(const QString& name, const IContactFluid* c) = 0;
 virtual bool validate( ContactModelFluidState* state ) = 0;
 virtual bool checkActivity( ContactModelFluidState* state ) const = 0;
 double getActivityDistance() const override { return 0.0; }
 
 virtual double getMassDensity() const {
 throw Exception( "Fluid model %1 does not have property mass density.", getName() );
 }
 virtual double getMobilityCoef() const { 
 throw Exception( "Fluid model %1 does not have property mobility coefficient.", getName() );
 }
 virtual double getSpecificStorage() const {
 throw Exception( "Fluid model %1 does not have property specific storage.", getName() );
 }
 
 // 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 {}
 };
}
#ifndef __GNUC__
#pragma warning(pop)
#endif