power_law_hydrodynamic_interaction_law.h

Go to the documentation of this file.

#ifndef KRATOS_SDEM_POWER_LAW_FLUID_HYDRODYNAMIC_INTERACTION_LAW_H
#define KRATOS_SDEM_POWER_LAW_FLUID_HYDRODYNAMIC_INTERACTION_LAW_H
 
#include "includes/define.h"
#include "includes/serializer.h"
#include "includes/model_part.h"
#include "containers/flags.h"
 
#include "hydrodynamic_interaction_law.h"
 
namespace Kratos {
 
class KRATOS_API(SWIMMING_DEM_APPLICATION) PowerLawFluidHydrodynamicInteractionLaw : public HydrodynamicInteractionLaw {
 
public:
    typedef Node NodeType;
 
    // Pointer types for PowerLawFluidHydrodynamicInteractionLaw
    KRATOS_CLASS_POINTER_DEFINITION(PowerLawFluidHydrodynamicInteractionLaw);
 
    PowerLawFluidHydrodynamicInteractionLaw()
        : HydrodynamicInteractionLaw(){}
 
    PowerLawFluidHydrodynamicInteractionLaw(Properties::Pointer pProp, Parameters& r_hydrodynamic_parameters)
        : HydrodynamicInteractionLaw(pProp, r_hydrodynamic_parameters){}
 
    PowerLawFluidHydrodynamicInteractionLaw(const PowerLawFluidHydrodynamicInteractionLaw &rPowerLawFluidHydrodynamicInteractionLaw)
        : HydrodynamicInteractionLaw(rPowerLawFluidHydrodynamicInteractionLaw){}
 
    void Initialize(const ProcessInfo& r_process_info) override;
 
    std::string GetTypeOfLaw() override;
 
 
    virtual ~PowerLawFluidHydrodynamicInteractionLaw();
 
    PowerLawFluidHydrodynamicInteractionLaw::Pointer Clone() const;
 
    double ComputeShahParticleReynoldsNumber(const double particle_radius,
                                             const double fluid_density,
                                             const double consistency_index,
                                             const double flow_behavior_index,
                                             const double modulus_of_minus_slip_velocity);
 
    void ComputeDragForce(SphericParticle* p_particle,
                          double particle_radius,
                          double fluid_density,
                          double fluid_kinematic_viscosity,
                          array_1d<double, 3>& minus_slip_velocity,
                          array_1d<double, 3>& drag_force,
                          const ProcessInfo& r_current_process_info);
 
private:
 
    friend class Serializer;
 
    virtual void save(Serializer& rSerializer) const override {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, Flags)
 
    }
 
    virtual void load(Serializer& rSerializer) override {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, Flags)
    }
 
}; // Class PowerLawFluidHydrodynamicInteractionLaw : public MainCL
 
KRATOS_DEFINE_APPLICATION_VARIABLE(SWIMMING_DEM_APPLICATION, PowerLawFluidHydrodynamicInteractionLaw::Pointer, SDEM_POWER_LAW_FLUID_HYDRODYNAMIC_INTERACTION_LAW_POINTER)
 
 
} // Namespace Kratos
 
#endif // KRATOS_SDEM_POWER_LAW_FLUID_HYDRODYNAMIC_INTERACTION_LAW_H