TALK=T;RUN( 1, 1)
  DISPLAY                                                                       
                                                                                
  Simulation of laminar steady flow around an inclined tube
  in a uniform stream. The tube is presented in the form of
  two coaxial cylinders. The external cylinder has the material
  with number 198 and the internal cylinder has the material
  of domain.
  
  The grid is cartesian and rather coarse (20*20*20), but the flow 
  follows the curved surface smoothly because the SPARSOL technique
  is used.
  
  The fluid is atmospheric air, treated as incompressible.
  
  The variables computed are:
  * p1, the pressure,
  * u1, the x-direction velocity,
  * v1, the y-direction velocity, and
  * w1, the z-direction velocity.
  
  This file contains a macro of commands which cause the Viewer 
  (when the macro button is pressed) to display contours 
  automatically.
                                                                                
  ENDDIS                                                                        
  VRV USE
  * Start of frame
  OBJECT PAINT TYPE BLOCKAGE ON
  OBJECT WIREFRAME TYPE BLOCKAGE ON
  PAUSE
  * Start of frame
  OBJECT PAINT TYPE BLOCKAGE OFF
  OBJECT WIREFRAME TYPE BLOCKAGE ON
  VARIABLE  Pressure; CON ON
  msg Pressure contours
  PAUSE
  * Start of frame
  VARIABLE  Velocity; VEC ON
  msg Pressure contours superimposed by velocity vectors
  PAUSE
  ENDUSE
 ************************************************************
   Q1 created by VDI menu, Version 2014, Date 23/01/14
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN = 1 ;LIBREF = 802
 ************************************************************
  Group 1. Run Title
 TEXT(Flow around an inclined tube)
   save1begin
           Initial data of problem
           -----------------------
          delaration of parameters
REAL(DIAMEX,DIAMIN,ANGL,THICK,ZSIZE,ANGLY,ANGLZ) ! delaration
REAL(POSX,POSY,POSZ,DMSIZE,ZDMSZ,VELIN)          ! of parameters
                   ! settings
DIAMEX=2.0         ! diameter of external cylinder
DIAMIN=1.6         ! diameter of internal cylinder
THICK=0.5*(DIAMEX-DIAMIN) ! TUBE wall thickness
ZSIZE=3.0          ! TUBE length
ANGLY=10.          ! TUBE rotation angle about Y axis
ANGLZ=20.          ! TUBE rotation angle about Z axis
POSX=0.8           ! X start position of TUBE
POSY=1.5           ! Y start position of TUBE
POSZ=3.0           ! Z start position of TUBE
DMSIZE=4.0         ! domain size in X and Y direction
ZDMSZ=10.0         ! domain size in Z direction
VELIN=1.0          ! inlent velocity in Z direction
  save1end
 ************************************************************
  Group 2. Transience
 STEADY = T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,NX,NY,NZ)
XULAST=DMSIZE      ! PIL variables
YVLAST=DMSIZE
ZWLAST=ZDMSZ
NX=20
NY=20
NZ=20
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 NAME(150) =PRPS
    * Solved variables list
 SOLVE(P1,U1,V1,W1)
    * Stored variables list
 STORE(PRPS)
    * Additional solver options
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 SOLUTN(V1,Y,Y,Y,N,N,Y)
 SOLUTN(W1,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  Group 9. Properties
 RHO1 =1.
 ENUL =0.05
 CP1 =1005.
 ENUT =0.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(W1)=1. ;FIINIT(PRPS)=-1.
   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 = 150
 RESFAC =1.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,1.)
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 PARSOL = F ! swich off Parsol
 ISG62 = 1  ! make active SParsol
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=0.3522 ;EX(U1)=0.
 EX(V1)=0.07498 ;EX(W1)=1.022
 EX(PRPS)=5.237
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 10 ;IYMON = 10 ;IZMON = 10
 NPRMON = 100000
 NPRMNT = 1
 TSTSWP = -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000
 ISWPRF = 1 ;ISWPRL = 100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,-0.941132,0.281294,-0.187468)
 GVIEW(UP,0.246784,0.950729,0.187649)
 GVIEW(DEPTH,999.)
 GVIEW(NEARPLANE,0.02)
 GVIEW(VDIS,5.258137)
 GVIEW(CENTRE,1.,2.,5.)
 
> DOM,    SIZE,        XULAST      , YVLAST      , ZWLAST
> DOM,    MONIT,       XULAST*0.5  , YVLAST*0.5  , ZDMSZ*0.5
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
  > GRID,   RSET_X_1,     NX, 1.000000E+00
  > GRID,   RSET_Y_1,     NY, 1.000000E+00
  > GRID,   RSET_Z_1,     NZ, 1.000000E+00
 
> OBJ,    NAME,        INLET
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        XULAST      , YVLAST      , 0.000000E+00
> OBJ,    GEOMETRY,    cube3t
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,    0.
> OBJ,    VELOCITY,    0. ,0. ,VELIN
 
> OBJ,    NAME,        OUTLET
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, ZWLAST
> OBJ,    SIZE,        XULAST      , YVLAST      , 0.000000E+00
> OBJ,    GEOMETRY,    cube12t
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,    0.
> OBJ,    COEFFICIENT, 1000.
 
> OBJ,    NAME,        EXTERNAL
> OBJ,    POSITION,    POSX        , POSY        , POSZ
> OBJ,    SIZE,        DIAMEX      , DIAMEX      , ZSIZE
> OBJ,    ROT-ANGLE,   0.000000E+00, ANGLY       , ANGLZ
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cylinder
> OBJ,    GRID,        N,N,N
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        INTERNAL
> OBJ,    POSITION,    POSX+THICK  , POSY+THICK  , POSZ
> OBJ,    SIZE,        DIAMIN      , DIAMIN      , ZSIZE
> OBJ,    ROT-ANGLE,   0.000000E+00, ANGLY       , ANGLZ
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cylinder
> OBJ,    VISIBLE,     NO
> OBJ,    GRID,        N,N,N
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    DOMAIN
STOP