base flow and fluctuations?

Emmanuel Germaine, 2009-02-10 07:47

I was working with the channel flow code V2.1.2. Last week, I have installed the latest version available on the wiki. 1.3.3. It seems that this version does not generate the same results than the previous that I used. If I save the velocity field using :

FlowField u(Nx,Ny,Nz,3,Lx,Lz,a,b);
...
 for(int nz=0;nz<u.Nz();++nz){
   for(int ny=0;ny<u.Ny();++ny){
     for(int nx=0;nx<u.Nx();++nx){
       fileName <<  x(nx) << ' ' << y(ny) << ' ' << z(nz) <<' ' << ' ';
       fileName <<  u(nx,ny,nz,0) << ' ' << u(nx,ny,nz,1) << ' '<< u(nx,ny,nz,2) <<  endl ;
     }
   }
 }

The version v1.2.2 saved an instantaneous field. The new version saves a velocity field that looks centered on zero (there is no mean). I would like to work with the instantaneous velocity field instead of its fluctuating part. Do you know if there is a way to do this ?

thanks!

~~DISCUSSION~~