TALK=T;RUN( 1, 1)

  DISPLAY

  DELTA-WINGLETS BEHIND A CYLINDER

  PARSOL is used here aimed at the demonstration of the method
  for the simulation of turbulent 3D flow over a cylinder with
  attached Delta-Wing Pair.

  The demonstration case considered consists of a half of cylinder
  placed in an air flow   within a rectangular duct. The delta-
  winglet is mounted on the wall behind a cylinder.   The position
  of the wingle tip is a half of cylinder diameter relative to the
  center of cylinder.The right triangle winglet with a hight equal
  to the duct width has an angle   of attack of 45 degree.

  The pressure difference between the pressure and the suction
  sides of the winglet is   anticipated to generate longitudinal
  vortex which supresses the flow separation on   cylinder and
  develops well pronounced secondary flow.

  The task is to calculate the flow and pressure fields along with
  all related field distributions.


  ENDDIS

 ************************************************************
   Q1 created by VDI menu, Version 3.3, Date 10/02/00
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN   =       1 ;LIBREF =     290
 ************************************************************
  Group 1. Run Title
 TEXT(Delta-winglets behind a cylinder)
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,60,60,10,7.620000E-05)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast
            name
 XSI= 1.200000E+00; YSI= 5.000000E-01; ZSI= 2.000000E-01
 RSET(D,CHAM    )
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(147) =SMOK ; NAME(148) =EPKE
 NAME(149) =LEN1 ; NAME(150) =ENUT
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  ,W1  ,SMOK)
    * Stored variables list
 STORE(ENUT,LEN1,EPKE)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,N)
 TURMOD(KEMODL)

 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  Group 9. Properties
 SETPRPS(1,  0)
 RHO1    = 1.189000E+00
 PRESS0  = 1.000000E+05
 CP1     = 1.005000E+03
 ENUL    = 1.544000E-05
 DVO1DT  = 3.410000E-03
 PRT   (EP  ) =  1.314000E+00
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(U1  ) =  2.000000E+01
   No PATCHes used for this Group


 INIADD  =    F
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
   No PATCHes used for this Group

 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =     500
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER (U1  ) =   20 ;LITER (V1  ) =   20
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 3.000000E-01)
 RELAX(U1  ,FALSDT, 1.000000E-03)
 RELAX(V1  ,FALSDT, 1.000000E-02)
 RELAX(W1  ,FALSDT, 1.000000E-02)
 RELAX(KE  ,LINRLX, 3.000000E-01)
 RELAX(EP  ,LINRLX, 3.000000E-01)
 RELAX(SMOK,LINRLX, 1.000000E+00)
 KELIN   =       3
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 GENK    =    T
 ASAP    =    T
 LSG3    =    T
 CSG3    ='FGEM'
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(W1  ,Y,Y,Y,N,Y,Y)
 OUTPUT(SMOK,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =       4 ;IYMON  =      32 ;IZMON  =       5
 NPRMON  =     100
 NPRMNT  =       1
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =  100000
 NXPRIN  =       1
 NYPRIN  =       1
 IPLTF   =       1 ;IPLTL  =     100
 NPLT    =       1
 ISWPRF  =       1 ;ISWPRL =  100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 NOWIPE  =    T

 GVIEW(P,-3.389963E-01,3.455267E-01,8.750392E-01)
 GVIEW(UP,1.248200E-01,9.384090E-01,-3.221937E-01)

> DOM,    SIZE,        1.200000E+00, 5.000000E-01, 2.000000E-01
> DOM,    MONIT,       7.000000E-02, 2.625000E-01, 9.000000E-02
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    5.000000E-03
> DOM,    VARIABLES,  SMOK
> DOM,    RELAX,       1.000000E+00

> OBJ1,   NAME,        OUTLET
> OBJ1,   POSITION,    1.200000E+00, 0.000000E+00, 0.000000E+00
> OBJ1,   SIZE,        0.000000E+00, 5.000000E-01, 2.000004E-01
> OBJ1,   CLIPART,     cube12t
> OBJ1,   ROTATION24,        1
> OBJ1,   VISIBLE,    -1.000000E+00
> OBJ1,   GRID,        2.000000E+00
> OBJ1,   TYPE,        OUTLET
> OBJ1,   PRESSURE,      0.000000E+00
> OBJ1,   TEMPERATURE,  -1.026000E+04
> OBJ1,   COEFFICIENT,   1.000000E+03
> OBJ1,   TURBULENCE,   -1.026000E+04,-1.026000E+04
> OBJ1,   SMOKE,         0.000000E+00

> OBJ2,   NAME,        IN3
> OBJ2,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ2,   SIZE,        0.000000E+00, 5.000000E-01, 2.000004E-01
> OBJ2,   CLIPART,     cube3t
> OBJ2,   ROTATION24,        1
> OBJ2,   GRID,        2.000000E+00
> OBJ2,   TYPE,        INLET
> OBJ2,   PRESSURE,      1.189000E+00
> OBJ2,   VELOCITY,      2.000000E+01, 0.000000E+00, 0.000000E+00
> OBJ2,   TEMPERATURE,   0.000000E+00
> OBJ2,   SMOKE,         0.000000E+00
> OBJ2,   TURB-INTENS,   5.000000E+00

> OBJ3,   NAME,        B12
> OBJ3,   POSITION,    8.320003E-01, 1.060000E-01, 0.000000E+00
> OBJ3,   SIZE,        1.080001E-01, 1.500000E-02, 2.000000E-01
> OBJ3,   ARBORIEN,    3.000000E+01, 0.000000E+00, 0.000000E+00
> OBJ3,   CLIPART,     wedge1
> OBJ3,   ROTATION24,        2
> OBJ3,   GRID,        2.000000E+00
> OBJ3,   TYPE,        BLOCKAGE
> OBJ3,   MATERIAL,      103
> OBJ3,   FIXED_TMP,     0.000000E+00, 0.000000E+00

> OBJ4,   NAME,        B5
> OBJ4,   POSITION,    4.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ4,   SIZE,        4.000000E-01, 2.000000E-01, 2.000000E-01
> OBJ4,   CLIPART,     ahcyl
> OBJ4,   ROTATION24,        4
> OBJ4,   GRID,        2.000000E+00
> OBJ4,   TYPE,        BLOCKAGE
> OBJ4,   MATERIAL,      198
STOP
distil=T
EX(P1  )=   1.031E+02
EX(U1  )=   2.012E+01
EX(V1  )=   2.710E+00
EX(W1  )=   3.723E-01
EX(KE  )=   2.485E+00
EX(EP  )=   2.429E+02
EX(PRPS)=   1.974E+01
EX(EPKE)=   3.875E+01
EX(LEN1)=   8.973E-03
EX(ENUT)=   6.369E-03