Post processing toolbox for 2d and 3d DNS and experimental results of turbulent reacting flows

Documentation and tutorial - Example and guide to use the tools

Example for the investigation of a 3 dimensional non-premixed, turbulent H2/air flame.

• First define and locate the flamefront:

  • e.g. along the stoichiometric isolevel of the mixture fraction f at f=0.5
    between the upper boundary Z=0.8 and the lower boundary Z=0.1:
    [og,ug]=<a title="flamefront3" href="index.php?site=lss&id=324&lang=de">flamefront3</a> (X,Y,Z,f,0.5,'abs',0.8,0.1)

  • and calculate the flame surface area
    [Area]=<a title="surfarea" href="index.php?site=lss&id=326&lang=de">surfarea</a> (X,Y,Z,f,0.5)

• Second compute or extract important variables along the flamefront

  • e.g. extract the flame index along flamefront defined by f=0.5;
    [inn,variso]=<a title="flameindex3" href="index.php?site=lss&id=396&lang=de">flameindex3</a> (X,Y,Z,Y_H2,Y_O2,f,'abs','yes')
  • compute strainrate along flamefront:
    [str3ff]=<a title="strnlev3" href="index.php?site=lss&id=438&lang=de">strnlev3</a> (X,Y,Z,strt,f,0.5,'abs')

• Then draw linear cuts through flame front between boundaries defined above

  • extract variables for example heatr along these cuts
    [XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX,XX1,YY1,ZZ1,VV1,VV2,VVG]=<a title="lincut3" href="index.php?site=lss&id=380&lang=de">lincut3</a> (X,Y,Z,0.0062,0.004887,0.0.0048,-0.4936,-0.8846,0.2071,f,heatr,0.0001,0.5,'abs',0.8,0.1);
  • compute integral value of heat release along the cut
    intvar=<a title="intcutlin3" href="index.php?site=lss&id=382&lang=de">intcutlin3</a> (X,Y,Z,0.0062,0.004887,0.0048,-0.4936,-0.8846,0.2071,f,heatr,0.0001,0.5,'abs',0.8,0.1)
  • compute flamethickness along linear cuts, use batch modus
    [flthick3,mathick,mithick,meanthick]=<a title="flamethicklin3b2" href="index.php?site=lss&id=336&lang=de">flamethicklin3b2</a> (X,Y,Z,f,0.5,'abs',0.8,0.1)

• Then draw nonlinear cuts through the flame front
  [XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX,XX1,YY1,ZZ1,VV1,VV2,VVG]=<a title="nonlincut3" href="index.php?site=lss&id=388&lang=de">nonlincut3</a> (X,Y,0.0062,0.004882,0.0048,f,Y_H2O2,0.0001,0.5,'abs',0.8,0.1)

  • extract variables for example the mass fraction of the radical H2O2 along the cuts
• Then do some statistical investigation
  • plot correlation between variables extracted along a linear cut, e.g. heat release and mass fraction of the radical OH:
    [cc]=<a title="correlationgen3" href="index.php?site=lss&id=424&lang=de">correlationgen3</a>(heatr,Y_OH)
  • compute mean, minimum, maximum and rms values of a variable e.g. heatrelease conditioned on f
    [cc,m,rms,mi,ma]=<a title="condmean3filt" href="index.php?site=lss&id=418&lang=de">condmean3filt</a> (X,Y,Z,f,heatr,0.01)
  • Investigate the velocity field. Plot amount and direction of the velocity vectors as cones.
    []=<a title="coneplot3" href="index.php?site=lss&id=452&lang=de">coneplot3</a> (x_velocity,y_velocity,z_velocity,X,Y,Z)