dallavalle_drag_law.h

Go to the documentation of this file.

// Author: Joaquin Gonzalez-Usua (jgonzalez@cimne.upc.edu)
// Date: April 2021
 
#if !defined(SDEM_DALLAVALLE_DRAG_LAW_H_INCLUDED)
#define SDEM_DALLAVALLE_DRAG_LAW_H_INCLUDED
#include <algorithm>
#include "stokes_drag_law.h"
 
namespace Kratos {
 
    class KRATOS_API(SWIMMING_DEM_APPLICATION) DallavalleDragLaw : public StokesDragLaw {
 
    public:
        typedef Node NodeType;
        KRATOS_CLASS_POINTER_DEFINITION(DallavalleDragLaw);
 
        DallavalleDragLaw(): StokesDragLaw(){}
 
        DallavalleDragLaw(Parameters& r_parameters): StokesDragLaw(r_parameters){}
 
        ~DallavalleDragLaw(){}
 
        DragLaw::Pointer Clone() const override;
 
        void Initialize(const ProcessInfo& r_process_info) override;
 
        std::string GetTypeOfLaw() override;
 
        void ComputeForce(SphericParticle* p_particle,
                          const double reynolds_number,
                          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) override;
 
        double CalculateEquivalentDiameter(SphericParticle* p_particle);
 
        double GetParticleMassFraction(SphericParticle* p_particle);
 
        double CalculateWeightingSum(SphericParticle* p_particle, const double& equivalent_diameter);
 
    private:
 
        friend class Serializer;
 
        virtual void save(Serializer& rSerializer) const override {
            KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, DragLaw)
        }
 
        virtual void load(Serializer& rSerializer) override {
            KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, DragLaw)
        }
 
    }; //class DallavalleDragLaw
 
} // Namespace Kratos
 
#endif /* SDEM_DALLAVALLE_DRAG_LAW_H_INCLUDED  defined */