The worked example is based on library case B534. This is a listing of case B534 after passing through the steps described in Chapter 2.
TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.5, Date 12/11/02
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
PHOTON USE
p
use patgeo
msg Geometry
vec x 1 sh
MSG Velocity vectors
msg
msg Press return to plot pressure contours
pause
cont p1 x 1 fil;.001
msg
msg Press return to plot streamlines
pause
clear
stream sh
x 1
y
posit
3
0.51268E+02 0.20453E+04 CR
0.38451E+02 0.19150E+04 CR
0.44859E+02 0.17847E+04 CR
0.76901E+02 0.16740E+04 CR
0.57676E+02 0.15502E+04 CR
0.57676E+02 0.14200E+04 CR
0.51268E+02 0.13548E+04 CR
0.11792E+04 0.12571E+04 t
exit
use patgeo
msg
msg Type e to End
ENDUSE
DISPLAY
The incompressible, single-phase flow of water through a
fully-open axi-symmetric ball valve is solved. The pipe-work
considered is 2 m in length; the radius at inlet is 0.15 m; it
diminishes to a minimum of 0.09 m in front of the ball; and the
maximum radius is 0.16 m.
A plug profile of axial velocity is prescribed at the inlet,
and a fixed-pressure condition is employed at outlet. Wall
friction is activated along the boundary of the ball, and along
the pipe wall. A fixed turbulent kinematic viscosity is set to
100 times the laminar value (the Reynolds number is of order
1.0E5).
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 1
************************************************************
Group 1. Run Title
TEXT(FLOW THROUGH A BALL VALVE : B534 )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,12,30)
* Set overall domain extent:
* xulast yvlast zwlast
name
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 1.000000E+00; YSI= 1.000000E+00; ZSI= 1.000000E+00
RSET(D,CHAM )
* Set objects: x0 y0 z0
* dx dy dz
name
XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 3.333333E-01
XSI= 1.000000E+00; YSI= 8.333337E-02; ZSI= 5.666666E-01
RSET(B,CMP0 )
XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
XSI= 1.000000E+00; YSI= 1.000000E+00; ZSI= 0.000000E+00
RSET(B,INLET )
XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 1.000000E+00
XSI= 1.000000E+00; YSI= 1.000000E+00; ZSI= 0.000000E+00
RSET(B,OUTLET )
XPO= 0.000000E+00; YPO= 1.000000E+00; ZPO= 0.000000E+00
XSI= 1.000000E+00; YSI= 0.000000E+00; ZSI= 1.000000E+00
RSET(B,WFUN )
XPO= 0.000000E+00; YPO= 8.333337E-02; ZPO= 3.333333E-01
XSI= 1.000000E+00; YSI= 0.000000E+00; ZSI= 5.666666E-01
RSET(B,VALVEWLL)
************************************************************
Group 6. Body-Fitted coordinates
BFC=T
READCO(grid4)
**********
NONORT = T
NCRT = 1
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(143) =REST ; NAME(144) =MOMZ
NAME(145) =MOMY ; NAME(146) =NPOR
NAME(147) =VPOR ; NAME(148) =WCRT
NAME(149) =VCRT ; NAME(150) =UCRT
* Solved variables list
SOLVE(P1 ,V1 ,W1 )
* Stored variables list
STORE(UCRT,VCRT,WCRT,VPOR,NPOR,MOMY,MOMZ,REST)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,N)
************************************************************
Group 8. Terms & Devices
DIFCUT = 0.000000E+00
************************************************************
Group 9. Properties
RHO1 = 1.000000E+03
ENUL = 1.000000E-06
CP1 = 1.000000E+00
ENUT = 1.000000E-04
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(W1 ) = 2.000000E+00 ;FIINIT(NPOR) = 1.000000E+00
FIINIT(VPOR) = 1.000000E+00 ;FIINIT(WCRT) = 1.001000E-10
FIINIT(VCRT) = 1.001000E-10 ;FIINIT(UCRT) = 1.001000E-10
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
INLET (INLET ,LOW ,2,0,0,0,0,0,1,1)
VALUE (INLET ,P1 , 2.000000E+03)
VALUE (INLET ,W1 , 2.000000E+00)
PATCH (GENPAT ,CELL ,0,0,0,0,0,0,1,1)
COVAL (GENPAT ,V1 , FIXFLU , GRND )
COVAL (GENPAT ,W1 , FIXFLU , GRND )
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 200
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 2.000000E-01)
RELAX(V1 ,FALSDT, 3.333333E-03)
RELAX(W1 ,FALSDT, 3.333333E-03)
RELAX(MOMZ,LINRLX, 7.000000E-01)
RELAX(MOMY,LINRLX, 7.000000E-01)
************************************************************
Group 18. Limits
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
L($G001
GENTR = T
*------------------------------------------------------
* GENTRA GROUP 1: Particle physics
*------------------------------------------------------
* Particle type - 30
GPTYPE = 30
* Gravity components in GENTRA Cartesian system
GGRAX = 0.000000E+00 ;GGRAY = 0.000000E+00
GGRAZ = 9.800000E+00
* Buoyancy forces
GBUOYA = F ;GSURPR = F
* Stochastic model of turbulence
GSTOCH = F
* Data for isothermal particles
GDRAG = GRND1
*-------------------------------------------------------
* GENTRA GROUP 2: Boundary conditions for particles
*-------------------------------------------------------
* Inlet-data file-name
GINFIL ='Q1'
<GENTRA-INLET-DATA>
*__YP__ZP__VP__WP__DI_LDEN_MDOT_(NUM)
0.01 0.0 0 1 0.001 500.0 1.0E-5
0.04 0.0 0 2 0.0001 1000.0 1.0E-5
0.07 0.0 0 3 0.0015 1000.0 1.0E-5
0.10 0.0 0 0.5 0.002 1000.0 1.0E-5
0.13 0.0 0 1 0.001 1500.0 1.0E-5
<END-GENTRA-INLET>
* Wall treatment, and rest coefficient if appropriate
GWALLC = 3
GWREST = 7.500000E-01
* Porosity threshold
GPOROS = 0.000000E+00
*-------------------------------------------------------
* GENTRA GROUP 3: Numerical controls
*-------------------------------------------------------
* 1st GENTRA sweep; frequency of calls
GSWEP1 = 190 ;GSWEPF = 1
* Maximum Lagrangian time-step; time step size multplier
GDTMAX = 1.000000E+00 ;GRTFRC = 7.000000E-01
* Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS = 5 ;GSTEMX = -100
GTIMMX = 1.000000E+01
*-------------------------------------------------------
* GENTRA GROUP 4: Output controls
*-------------------------------------------------------
* Restart-file, history-file and frequency for output
GRSFIL ='NONE'
GHFILE ='GHIS'
GOUTFR = 1
* The identifier of the individual history and
trajectory files
GH1STC ='NONE'
GT1STC ='T'
* The first, last particles and the interval for
writing history and trajectory files
NGWSTR = 1 ;NGWEND = 5 ;NGWINT = 1
GSWOUT = 1
* Particle number for residence-time calculation
GRESTI = 2
L(GENSET
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 2 ;IZMON = 20
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NYPRIN = 2
NZPRIN = 5
NPLT = 2
ISWPRF = 1 ;ISWPRL = 100000
PATCH (DOMAIN ,CONTUR,1,1,1,12,1,30,1,1)
PLOT(DOMAIN ,P1 , 0.000000E+00, 1.500000E+01)
PATCH (INNER ,PROFIL,1,1,2,2,1,30,1,1)
PLOT(INNER ,W1 , 0.000000E+00, 0.000000E+00)
PATCH (OUTER ,PROFIL,1,1,12,12,1,30,1,1)
PLOT(OUTER ,W1 , 0.000000E+00, 0.000000E+00)
PATCH (FRONT ,PROFIL,1,1,2,12,10,10,1,1)
PLOT(FRONT ,P1 , 0.000000E+00, 0.000000E+00)
PLOT(FRONT ,V1 , 0.000000E+00, 0.000000E+00)
PLOT(FRONT ,W1 , 0.000000E+00, 0.000000E+00)
PATCH (BACK ,PROFIL,1,1,2,12,28,28,1,1)
PLOT(BACK ,P1 , 0.000000E+00, 0.000000E+00)
PLOT(BACK ,V1 , 0.000000E+00, 0.000000E+00)
PLOT(BACK ,W1 , 0.000000E+00, 0.000000E+00)
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,-1.000000E+00,0.000000E+00,0.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 1.000000E+00, 1.200000E+01, 3.000000E+01
> DOM, MONIT, 1.000000E+00, 2.000000E+00, 2.000000E+01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, CMP0
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.000000E+01
> OBJ, SIZE, 1.000000E+00, 1.000000E+00, 1.700000E+01
> OBJ, CLIPART, CMP0
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 199
> OBJ, NAME, INLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 1.200000E+01, 0.000000E+00
> OBJ, CLIPART, INLET
> OBJ, TYPE, USER_DEFINED
> OBJ, NAME, OUTLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 3.000000E+01
> OBJ, SIZE, 1.000000E+00, 1.200000E+01, 0.000000E+00
> OBJ, CLIPART, OUTLET
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, -1.026000E+04
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 2.000000E+00
> OBJ, NAME, WFUN
> OBJ, POSITION, 0.000000E+00, 1.200000E+01, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 0.000000E+00, 3.000000E+01
> OBJ, CLIPART, WFUN
> OBJ, TYPE, PLATE
> OBJ, NAME, VALVEWLL
> OBJ, POSITION, 0.000000E+00, 1.000000E+00, 1.000000E+01
> OBJ, SIZE, 1.000000E+00, 0.000000E+00, 1.700000E+01
> OBJ, CLIPART, VALVEWLL
> OBJ, TYPE, PLATE
> OBJ, POROSITY, 0.000000E+00
> OBJ, SIDE, HIGH
STOP