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