TALK=F;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.6, Date 01/04/04
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
This case simulates a heat transfer by convection about
a single heated pipe which is described by a cylinder object.
The pipe allocates heat non-uniformly in the Z direction.
The pipe is cooled by a air flow of which moves
in the X direction from left to right.
The air flow is laminar. The Reynolds number equals 26.
The Prandtl number equals 0.7
Nusselt number calculated by the formula for a external cross
flow of a round pipe, located in depth of a pipe array,
Nu = 0.3*Re^0,6*Pr^1/3
equals 1.9.
The In-Form is used for the setting of a non-uniform heating
source per unit volume as
(SOURCE of TEM1 at PIPE is 1.e+5*ZG with VOLU)
The Q1 contains PHOTON USE commands
ENDDIS
PHOTON USE
p;;;
vi x
con tem1 x 33 fi;0.0001
gr x 33
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 146
Group 1. Run Title
TEXT(Heat exchange about a single pipe )
Group 2. Transience
STEADY = T
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,80,40,5)
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(149) =PRPS ; NAME(150) =TEM1
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 ,TEM1)
* Stored variables list
STORE(PRPS)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
Group 9. Properties
PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02
* Domain material is: 0 Air at 20 deg C, 1 atm
SETPRPS(1, 0)
ENUT = 0.000000E+00
DVO1DT = 3.410000E-03
PRNDTL(TEM1) = -2.580000E-02
** Reynolds number
REAL(REYNO,WIN,DIAM); WIN=0.01; DIAM=0.04; REYNO=WIN*DIAM/ENUL
REYNO
** Prandtl number
REAL(PRANDT,COND); COND=2.58E-02; PRANDT=CP1*ENUL*RHO1/COND
PRANDT
** Nusselt number by 0.3*REYNO**0.6*PRANDT**(1/3) formula
REAL(NUSN); NUSN=0.3*REYNO**0.6*PRANDT**(1/3)
NUSN
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS) = -1.000000E+00
Group 13. Boundary & Special Sources
Echo InForm settings for Group 13
INFORM13BEGIN
(SOURCE of TEM1 at PIPE is 1.e+5*ZG with VOLU)
INFORM13END
Group 15. Terminate Sweeps
LSWEEP = 300
RESFAC = 1.000000E-03
Group 16. Terminate Iterations
LITER (TEM1) = 600
ENDIT (TEM1) = 1.000000E-06
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
Group 20. Preliminary Printout
DISTIL = T
EX (P1 ) = 4.457000E-05 ;EX (U1 ) = 1.004000E-02
EX (V1 ) = 9.967000E-04 ;EX (W1 ) = 7.503000E-09
EX (PRPS) = 6.638000E+00 ;EX (TEM1) = 2.468000E+01
Group 22. Monitor Print-Out
IXMON = 29 ;IYMON = 20 ;IZMON = 6
TSTSWP = -1
GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 2.000000E-01, 1.000000E-01, 1.000000E-01
> DOM, MONIT, 7.125000E-02, 4.875000E-02, 5.500000E-02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, INLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E-01, 1.000000E-01
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 1.000000E-02, 0.000000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 1.800000E+01
> OBJ, NAME, OUTLET
> OBJ, POSITION, 2.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E-01, 1.000000E-01
> OBJ, GEOMETRY, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, NAME, PIPE
> OBJ, POSITION, 3.000000E-02, 3.000000E-02, 0.000000E+00
> OBJ, SIZE, 4.000000E-02, 4.000000E-02, 1.000000E-01
> OBJ, GEOMETRY, cylinder
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 111,STEEL at 27 deg c (C = 1%)
STOP