domain.h

Go to the documentation of this file.

 
#pragma once
 
#include "dv.h"
#include "domaincase.h"
#include "inputoutput.h"
#include "param.h"
#include "streams.h"
#include "micromixer.h"
#include "eddy.h"
#include "meshManager.h"
#include "solver.h"
#include "randomGenerator.h"
#include "cantera/thermo/IdealGasPhase.h"
#include "cantera/transport.h"
#include <vector>
#include <string>
#include <map>
 
using namespace std;
using namespace Cantera;
 
 
class domain {
 
    public:
 
 
        domain                  *domn;     
 
        int                     ngrd;      
        int                     ngrdf;     
 
        vector<dv*>             v;         
 
        dv*                     pos;       
        dv*                     posf;      
        dv*                     rho;
        dv*                     dvisc;
        dv*                     uvel;
        dv*                     vvel;
        dv*                     wvel;
        dv*                     sdiff;
        dv*                     sca;
        dv*                     phase;
        dv*                     enth;
        dv*                     temp;
        dv*                     mixf;
        dv*                     chi;
        dv*                     hr;
        dv*                     aDL;
        vector<dv*>::iterator   ysp;       
        vector<dv*>::iterator   svar;      
        vector<dv*>::iterator   eta;       
 
        map<string,dv*>         varMap;
 
        IdealGasPhase           *gas;        
        Transport               *tran;       
        streams                 *strm;       
        inputoutput             *io;         
        param                   *pram;       
        micromixer              *mimx;       
        eddy                    *ed;         
        domain                  *eddl;       
        solver                  *solv;       
        meshManager             *mesher;     
 
        randomGenerator         *rand;
 
        int                     nTrans;      
 
        domaincase              *domc;       
 
 
 
        int    domainPositionToIndex(double position, const bool LowSide, int dbg);
        void   setDomainFromRegion(const int i1, const int i2);
        double cyclePeriodicDomain(const int icycle);
        void   backCyclePeriodicDomain(const double backCycleDistance);
        double Ldomain();
 
    private:
 
        void initEddyDomain();
 
 
 
    public:
 
        void init(inputoutput     *p_io,
                  meshManager     *p_mesher,
                  streams         *p_strm,
                  IdealGasPhase   *p_gas,
                  Transport       *p_tran,
                  micromixer      *p_mimx,
                  eddy            *p_ed,
                  domain          *p_eddl,
                  solver          *p_solv,
                  randomGenerator *p_rand,
                  bool             LisEddyDomain=false);
        domain(domain *p_domn, param *p_pram);
        virtual ~domain() {
            for(int k=0; k<v.size(); k++)
                delete v.at(k);
            delete domc;     
        }
 
};