sinusoidal_porosity_solution_transient_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_SINUSOIDAL_POROSITY_SOLUTION_TRANSIENT_BODY_FORCE_PROCESS_H
#define KRATOS_SINUSOIDAL_POROSITY_SOLUTION_TRANSIENT_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) SinusoidalPorositySolutionTransientBodyForceProcess : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(SinusoidalPorositySolutionTransientBodyForceProcess);
 
 
    SinusoidalPorositySolutionTransientBodyForceProcess();
    SinusoidalPorositySolutionTransientBodyForceProcess(
        ModelPart& rModelPart);
 
    SinusoidalPorositySolutionTransientBodyForceProcess(
        ModelPart& rModelPart,
        Parameters& rParameters);
 
    SinusoidalPorositySolutionTransientBodyForceProcess(
        Model& rModel,
        Parameters& rParameters);
 
    ~SinusoidalPorositySolutionTransientBodyForceProcess() 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;
    double                                         mMaxGradAlpha;
    double                                           mWaveNumber;
    bool                                      mInitialConditions;
    bool                                 mAlternativeFormulation;
 
 
    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(
        double &rReynoldsNumber,
        double &viscosity);
 
    void CalculatePermeability(
        double &rDamKohlerNumber,
        double &dynamic_viscosity,
        double &permeability);
 
    void CalculateWaveNumber(
        double &mMaxGradAlpha,
        double &mReynoldsNumber,
        double &dynamic_viscosity,
        double &mDeltaAlpha,
        double &mUchar,
        double &mWaveNumber);
 
    void SetInitialBodyForceAndPorosityField();
 
    void SetBodyForceAndPorosityField();
 
    void SetFluidProperties();
 
 
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "SinusoidalPorositySolutionTransientBodyForceProcess" ;
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override {rOStream << "SinusoidalPorositySolutionTransientBodyForceProcess";}
 
    void PrintData(std::ostream& rOStream) const override {}
 
 
 
 
private:
 
 
 
 
 
 
 
 
 
 
 
 
}; // Class SinusoidalPorositySolutionTransientBodyForceProcess
 
 
 
 
 
};  // namespace Kratos.
 
#endif // KRATOS_SINUSOIDAL_POROSITY_SOLUTION_TRANSIENT_BODY_FORCE_PROCESS_H