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
************************************************************
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)
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

************************************************************
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