## Appendix E; Listing of the Q1 File for the Example

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