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