porosity_solution_and_sinusoidal_body_force_process.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Joaquin Gonzalez-Usua
//
//
 
#ifndef KRATOS_POROSITY_SOLUTION_AND_SINUSOIDAL_BODY_FORCE_PROCESS_H
#define KRATOS_POROSITY_SOLUTION_AND_SINUSOIDAL_BODY_FORCE_PROCESS_H
 
// System includes
#include <string>
#include <iostream>
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "processes/process.h"
#include "includes/kratos_parameters.h"
 
// Application includes
 
 
namespace Kratos
{
 
 
 
 
 
 
class KRATOS_API(SWIMMING_DEM_APPLICATION) PorositySolutionAndSinusoidalBodyForceProcess : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(PorositySolutionAndSinusoidalBodyForceProcess);
 
 
    PorositySolutionAndSinusoidalBodyForceProcess();
    PorositySolutionAndSinusoidalBodyForceProcess(
        ModelPart& rModelPart);
 
    PorositySolutionAndSinusoidalBodyForceProcess(
        ModelPart& rModelPart,
        Parameters& rParameters);
 
    PorositySolutionAndSinusoidalBodyForceProcess(
        Model& rModel,
        Parameters& rParameters);
 
    ~PorositySolutionAndSinusoidalBodyForceProcess() override {}
 
 
 
 
    void Execute() override;
 
    void ExecuteInitialize() override;
 
    void ExecuteBeforeSolutionLoop() override;
 
    void ExecuteInitializeSolutionStep() override;
 
    void ExecuteFinalizeSolutionStep() override;
 
 
    void CheckDefaultsAndProcessSettings(Parameters &rParameters);
 
    const Parameters GetDefaultParameters() const override;
 
    void CalculateKinematicViscosity();
 
    void CalculatePermeability(double &dynamic_viscosity);
 
    void SetInitialBodyForceAndPorosityField();
 
    void SetFluidProperties();
 
    bool IsInsideEllipticalSupport(
        const double x1,
        const double x2,
        const double c,
        const double R);
 
 
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "PorositySolutionAndSinusoidalBodyForceProcess" ;
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override {rOStream << "PorositySolutionAndSinusoidalBodyForceProcess";}
 
    void PrintData(std::ostream& rOStream) const override {}
 
    ModelPart&                                       mrModelPart;
    double                                              mDensity;
    double                                            mViscosity;
    double                                         mPermeability;
    double                                                mUchar;
    double                                           mDeltaAlpha;
    double                                               mLength;
    double                                   mMaxSqueezeFraction;
    double                                                mOmega;
    double                                     mSqueezeAmplitude;
    double                                              mNSafety;
    double                                             mX1Origin;
    double                                             mX2Origin;
    double                                       mReynoldsNumber;
    double                                      mDamKohlerNumber;
    bool                                      mInitialConditions;
    bool                                 mAlternativeFormulation;
 
 
private:
 
 
 
 
 
 
 
 
 
 
 
 
}; // Class PorositySolutionAndSinusoidalBodyForceProcess
 
 
 
 
 
};  // namespace Kratos.
 
#endif // KRATOS_POROSITY_SOLUTION_AND_SINUSOIDAL_BODY_FORCE_PROCESS_H