cellular_flow_partial_derivatives.h

Go to the documentation of this file.

#if !defined(KRATOS_CELLULAR_FLOW_PARTIAL_DERIVATIVES_H)
#define KRATOS_CELLULAR_FLOW_PARTIAL_DERIVATIVES_H
 
/* System includes */
#include <limits>
#include <iostream>
#include <iomanip>
 
/* External includes */
#ifdef _OPENMP
#include <omp.h>
#endif
 
/* Project includes */
#include "includes/define.h"
#include "includes/variables.h"
#include "includes/model_part.h"
 
namespace Kratos
{
class CellularFlowPartialDerivatives
{
public:
 
KRATOS_CLASS_POINTER_DEFINITION(CellularFlowPartialDerivatives);
 
 
CellularFlowPartialDerivatives(): mL(1.0), mU(0.0), mK(2.72), mOmega(Globals::Pi)
{
    mOneOverL = 1.0 / mL;
    mOmegaUOverL = mOmega * mU / mL;
}
 
CellularFlowPartialDerivatives(const double half_wavelength, const double max_flow_speed, const double oscillation_relative_amplitude, const double oscillation_angular_frequency)
                 :mL(half_wavelength), mU(max_flow_speed), mK(oscillation_relative_amplitude), mOmega(oscillation_angular_frequency)
{
    mOneOverL = 1.0 / mL;
    mOmegaUOverL = mOmega * mU / mL;
}
 
 
virtual ~CellularFlowPartialDerivatives(){}
 
 
//***************************************************************************************************************
//***************************************************************************************************************
 
void UpdateCoordinates(const double time, const array_1d<double, 3>& coor);
 
// Values
 
double U0();
 
double U1();
 
double U2();
 
// First-order derivatives
 
double U0DT();
double U0D0();
double U0D1();
double U0D2();
 
double U1DT();
double U1D0();
double U1D1();
double U1D2();
 
double U2DT();
double U2D0();
double U2D1();
double U2D2();
 
// Second-order derivatives
 
double U0DTDT();
double U0DTD0();
double U0DTD1();
double U0DTD2();
double U0D0D0();
double U0D0D1();
double U0D0D2();
double U0D1D1();
double U0D1D2();
double U0D2D2();
 
double U1DTDT();
double U1DTD0();
double U1DTD1();
double U1DTD2();
virtual double U1D0D0();
virtual double U1D0D1();
virtual double U1D0D2();
virtual double U1D1D1();
virtual double U1D1D2();
virtual double U1D2D2();
 
virtual double U2DTDT();
virtual double U2DTD0();
virtual double U2DTD1();
virtual double U2DTD2();
virtual double U2D0D0();
virtual double U2D0D1();
virtual double U2D0D2();
virtual double U2D1D1();
virtual double U2D1D2();
virtual double U2D2D2();
 
//***************************************************************************************************************
//***************************************************************************************************************
 
 
 
 
 
virtual std::string Info() const
{
    return "";
}
 
 
virtual void PrintInfo(std::ostream& rOStream) const
{
}
 
 
virtual void PrintData(std::ostream& rOStream) const
{
}
 
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
 
double mL;
double mU;
double mK;
double mOneOverL;
double mOmegaUOverL;
double mOmega;
double mSinOmegaT;
double mCosOmegaT;
double mSinPiX0;
double mCosPiX0;
double mSinPiX1;
double mCosPiX1;
 
 
 
 
 
 
 
 
CellularFlowPartialDerivatives & operator=(CellularFlowPartialDerivatives const& rOther);
 
 
 
}; // Class CellularFlowPartialDerivatives
 
 
 
 
 
} // namespace Kratos.
 
#endif // KRATOS_CELLULAR_FLOW_PARTIAL_DERIVATIVES_H defined