ethier_flow_field.h

Go to the documentation of this file.

#if !defined(KRATOS_ETHIER_FLOW_FIELD_H)
#define KRATOS_ETHIER_FLOW_FIELD_H
 
// /* External includes */
 
// System includes
 
// Project includes
#include "includes/variables.h"
 
/* System includes */
#include <limits>
#include <iostream>
#include <iomanip>
 
/* External includes */
#ifdef _OPENMP
#include <omp.h>
#endif
 
/* Project includes */
#include "includes/define.h"
#include "utilities/openmp_utils.h"
#include "real_functions.h"
#include "velocity_field.h"
 
 
namespace Kratos
{
class KRATOS_API(SWIMMING_DEM_APPLICATION) EthierFlowField : public VelocityField
{
public:
 
KRATOS_CLASS_POINTER_DEFINITION(EthierFlowField);
 
 
EthierFlowField():VelocityField(), mA(0.25 * Globals::Pi), mD(0.5 * Globals::Pi)
{
    unsigned int number_of_threads = ParallelUtilities::GetNumThreads();
    ResizeVectorsForParallelism(number_of_threads);
}
 
EthierFlowField(const double a, const double b)
                 :VelocityField(), mA(a), mD(b)
{
    unsigned int number_of_threads = ParallelUtilities::GetNumThreads();
    ResizeVectorsForParallelism(number_of_threads);
}
 
 
 
virtual ~EthierFlowField(){}
 
 
//***************************************************************************************************************
//***************************************************************************************************************
void ResizeVectorsForParallelism(const int n_threads) override;
 
void UpdateCoordinates(const double time, const array_1d<double, 3>& coor, const int i_thread = 0) override;
 
void UpdateCoordinates(const double time, const DenseVector<double>& coor, const int i_thread = 0) override;
 
void LockCoordinates(const int i_thread = 0) override;
 
void UnlockCoordinates(const int i_thread = 0) override;
//***************************************************************************************************************
//***************************************************************************************************************
 
 
 
 
 
virtual std::string Info() const override
{
    return "";
}
 
 
virtual void PrintInfo(std::ostream& rOStream) const override
{
}
 
 
virtual void PrintData(std::ostream& rOStream) const override
{
}
 
 
 
 
protected:
 
 
 
 
 
 
// Values
 
double U0(const int i_thread = 0) override;
double U1(const int i_thread = 0) override;
double U2(const int i_thread = 0) override;
 
// First-order derivatives
 
double U0DT(const int i_thread = 0) override;
double U0D0(const int i_thread = 0) override;
double U0D1(const int i_thread = 0) override;
double U0D2(const int i_thread = 0) override;
 
double U1DT(const int i_thread = 0) override;
double U1D0(const int i_thread = 0) override;
double U1D1(const int i_thread = 0) override;
double U1D2(const int i_thread = 0) override;
 
double U2DT(const int i_thread = 0) override;
double U2D0(const int i_thread = 0) override;
double U2D1(const int i_thread = 0) override;
double U2D2(const int i_thread = 0) override;
 
// Second-order derivatives
 
double U0DTDT(const int i_thread = 0) override;
double U0DTD0(const int i_thread = 0) override;
double U0DTD1(const int i_thread = 0) override;
double U0DTD2(const int i_thread = 0) override;
double U0D0D0(const int i_thread = 0) override;
double U0D0D1(const int i_thread = 0) override;
double U0D0D2(const int i_thread = 0) override;
double U0D1D1(const int i_thread = 0) override;
double U0D1D2(const int i_thread = 0) override;
double U0D2D2(const int i_thread = 0) override;
 
double U1DTDT(const int i_thread = 0) override;
double U1DTD0(const int i_thread = 0) override;
double U1DTD1(const int i_thread = 0) override;
double U1DTD2(const int i_thread = 0) override;
double U1D0D0(const int i_thread = 0) override;
double U1D0D1(const int i_thread = 0) override;
double U1D0D2(const int i_thread = 0) override;
double U1D1D1(const int i_thread = 0) override;
double U1D1D2(const int i_thread = 0) override;
double U1D2D2(const int i_thread = 0) override;
 
double U2DTDT(const int i_thread = 0) override;
double U2DTD0(const int i_thread = 0) override;
double U2DTD1(const int i_thread = 0) override;
double U2DTD2(const int i_thread = 0) override;
double U2D0D0(const int i_thread = 0) override;
double U2D0D1(const int i_thread = 0) override;
double U2D0D2(const int i_thread = 0) override;
double U2D1D1(const int i_thread = 0) override;
double U2D1D2(const int i_thread = 0) override;
double U2D2D2(const int i_thread = 0) override;
 
 
 
 
 
 
 
 
 
private:
 
 
 
double mA;
double mD;
std::vector<int> mCoordinatesAreUpToDate;
std::vector<double> mExpD2T;
std::vector<double> mExpAX;
std::vector<double> mExpAZ;
std::vector<double> mExpAY;
std::vector<double> mSinAXDY;
std::vector<double> mCosAXDY;
std::vector<double> mSinAYDZ;
std::vector<double> mCosAYDZ;
std::vector<double> mSinAZDX;
std::vector<double> mCosAZDX;
 
 
 
 
 
 
 
 
 
 
EthierFlowField & operator=(EthierFlowField const& rOther);
 
 
 
}; // Class EthierFlowField
 
 
 
 
 
} // namespace Kratos.
 
#endif // KRATOS_ETHIER_FLOW_FIELD_H defined