eddy.h

Go to the documentation of this file.

 
#pragma once
 
#include <vector>
 
class domain;
 
using namespace std;
 
 
class eddy {
 
    public:
 
 
        domain              *domn;              
        domain              *eddl;              
 
        double              eddySize;           
        double              leftEdge;           
        double              rightEdge;          
        double              invTauEddy;         
        double              Pa;                 
        bool                LperiodicEddy;      
        vector<double>      cCoef;              
        vector<double>      bCoef;              
        vector<double>      K;                  
        vector<double>      dxc;                
        vector<double>      pos0;               
 
        double              esdp1;              
        double              esdp2;
        double              esdp3;
        double              esdp4;
 
        vector<double> cca,ccb,ccc,ccd;                 
 
 
        void   sampleEddySize();
        void   sampleEddyPosition();
        void   tripMap(domain *line,    const int iS, int iE, const double C, const bool LsplitAtEddy=false);
        bool   eddyTau(const double Z_value, const double C);
        void   computeEddyAcceptanceProb(const double dtSample);
        void   applyVelocityKernels(domain *line, const int iS, const int iE);
 
    private:
 
        void   fillKernel();
        void   fillKernel_planarAnalytic();
        double eddyFavreAvgVelocity(const vector<double> &dxc);
        void   set_kernel_coefficients();
 
 
    public:
 
        eddy(){}
        void init(domain *p_domn, domain *p_eddl);
        ~eddy(){}
 
};