TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2010, Date 03/03/10
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
Whole-object heat sources which vary with cross-stream position
are applied by In-Form to an array of steel pipes held at right
angles to a cold air stream.
Use is illustrated of both of the two methods of applying sources
to facetted objects, whereby method 1 uses the statement
(SOURCE of TEM1 at PIPEx is 10*YG with WHOL)
and method 2 uses the statement
> OBJ, NAME, PIPEx
> OBJ, INFSRC_TEM1, 10*YG with WHOL
The same dual approach is used for inital values, which are
specified for some of the objects via:
(INITIAL of TEM1 at PIPEx is INTEM)
and for others via
> OBJ, INFINI_TEM1, INTEM
The results are the same whichever method is used.
In-Form is also used to calculate the average Nusselt number and
to print it in a special INFOROUT file.
The air flow is laminar. The Reynolds number equals 26.
The Prandtl number equals 0.71
Nusselt number calculated by the formula for flow at right angles
to a cylinder according to the textbook formula:
Nu = 0.3*Re^0.6*Pr^1/3
equals 1.9.
ENDDIS
PHOTON USE
p;;;
vi z
con tem1 z 1 fi;0.0001
enduse
************************************************************
IRUNN = 1 ;LIBREF = 147
************************************************************
Group 1. Run Title
TEXT(Heat exchange about a array of pipes. )
************************************************************
Echo save-block settings for Group 1
SAVE1BEGIN
Declarations and settings
REAL(INTEM,DIAM,DOMPRPS)
INTEM=18.0 ! inlet temperature
DIAM=0.04 ! rod diameter
DOMPRPS=0.0 ! domain material index, constant-property air
SAVE1END
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,100,100,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(144)=DEN1 ;NAME(145)=ENUL
NAME(146)=KOND ;NAME(147)=CP1
NAME(148)=PRAN ;NAME(149)=PRPS
NAME(150) =TEM1
* Solved variables list
SOLVE(P1,U1,V1,TEM1)
* Stored variables list
STORE(PRPS,PRAN,CP1,KOND,ENUL,DEN1)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0 =1.0E+05 ;TEMP0 =273.
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
SETPRPS(1, 0)
ENUT =0.
DVO1DT =3.41E-03
PRNDTL(TEM1)=-0.0258
************************************************************
Echo save-block settings for Group 9
save9begin
ENUT =0.
SETPRPS(1,DOMPRPS) ! set pil properties as those of material
! with index DOMPRPS
save9end
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS)=-1. ;FIINIT(TEM1)=18.
No PATCHes used for this Group
INIADD = F
************************************************************
Echo save-block settings for Group 11
SAVE11BEGIN
FIINIT(TEM1) = INTEM
!!!! note that the some init statements are inactive because
they do not start in the first or second column.
Their effects are provided by the infini_ lines
in the VR-Object-related settings below.
(INITIAL of TEM1 at PIPE1 is INTEM)
(INITIAL of TEM1 at PIPE2 is INTEM)
(INITIAL of TEM1 at PIPE3 is INTEM)
(INITIAL of TEM1 at PIPE4 is INTEM)
(INITIAL of TEM1 at PIPE5 is INTEM)
(INITIAL of TEM1 at PIPE6 is INTEM)
(INITIAL of TEM1 at PIPE7 is INTEM)
SAVE11END
************************************************************
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
************************************************************
Echo save-block settings for Group 13
SAVE13BEGIN
Echo InForm settings for Group 13
!!!! note that the some source statements are inactive because
they do not start in the first or second column.
Their effects are provided by the infsrc_ lines
in the VR-Object-related settings below.
Heat source
(SOURCE of TEM1 at PIPE1 is 10*YG with WHOL)
(SOURCE of TEM1 at PIPE2 is 10*YG with WHOL)
(SOURCE of TEM1 at PIPE3 is 10*YG with WHOL)
(SOURCE of TEM1 at PIPE4 is 10*YG with WHOL)
(SOURCE of TEM1 at PIPE5 is 10*YG with WHOL)
(SOURCE of TEM1 at PIPE6 is 10*YG with WHOL)
(SOURCE of TEM1 at PIPE7 is 10*YG with WHOL)
Next In-Form statements are used for calculation
average Nusselt number
** Pipe length
REAL(PLEN); PLEN=0.1
** Reynolds number
REAL(REYNO,WIN); WIN=0.01; 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(NUS0); NUS0=0.3*REYNO**0.6*PRANDT**(1/3)
NUS0
** Summa of temperatures in solid
(MAKE STSL is 0.)
(STORE1 of STSL is SUM(TEM1) with swpfin!if(prps.eq.111))
** Number of solid cells
(MAKE SNSL is 0.)
(STORE1 of SNSL is SUM(1) with swpfin!if(prps.eq.111))
** Average solid temperature
(MAKE AVTS is 0.)
(STORE1 of AVTS is STSL/SNSL with swpfin)
** The heat flux supplied to one average pipe
This is linear function of Y coordinate: 10*Y
The value of heat flux for first pipe is 10.*0.05
where 0.05 is y coordinate of centre of PIPE1
REAL(HTFL); HTFL=10.*0.05
** Surface of one average pipe
REAL(PARE,PI); PI=3.1415927; PARE=PI*DIAM*PLEN
** Nusselt number for average temteratures
(MAKE NUSS is 0.)
(STORE1 of NUSS is :HTFL:*:DIAM:/(:COND:*(AVTS-:INTEM:)*:PARE:) wit$
h swpfin)
** Print-out into INFOROUT file
(PRINT of Solid_Temp_Summ is STSL)
(PRINT of Sol_Number_Cell is SNSL)
(PRINT of Averag_Sol_Temp is AVTS)
(PRINT of Average_Nusselt is NUSS)
STORE(CP1,KOND,ENUL,DEN1) ! provision for printing Prandtl
(STORED VAR PRAN IS CP1*DEN1*ENUL/KOND with imat<100)
SAVE13END
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1000
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(TEM1)=200
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
CONWIZ=T ! make active the convergence wizard
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=0.01106 ;EX(U1)=0.01004
EX(V1)=1.97E-03 ;EX(DEN1)=572.700012
EX(ENUL)=0.022 ;EX(KOND)=3.175
EX(CP1)=966. ;EX(PRAN)=0.6628
EX(PRPS)=9.063 ;EX(TEM1)=43.93
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 66 ;IYMON = 52 ;IZMON = 1
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.,0.,1.)
GVIEW(UP,0.,1.,0.)
> DOM, SIZE, 3.100000E-01, 3.100000E-01, 1.000000E-01
> DOM, MONIT, 2.010526E-01, 1.550000E-01, 5.000000E-02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 10, 1.000000E+00
> GRID, RSET_X_2, 13, 1.000000E+00
> GRID, RSET_X_3, 7, 1.000000E+00
> GRID, RSET_X_4, 13, 1.000000E+00
> GRID, RSET_X_5, 57, 1.000000E+00
> GRID, RSET_Y_1, 10, 1.000000E+00
> GRID, RSET_Y_2, 11, 1.000000E+00
> GRID, RSET_Y_3, 2, 1.000000E+00
> GRID, RSET_Y_4, 10, 1.000000E+00
> GRID, RSET_Y_5, 2, 1.000000E+00
> GRID, RSET_Y_6, 10, 1.000000E+00
> GRID, RSET_Y_7, 2, 1.000000E+00
> GRID, RSET_Y_8, 9, 1.000000E+00
> GRID, RSET_Y_9, 2, 1.000000E+00
> GRID, RSET_Y_10, 9, 1.000000E+00
> GRID, RSET_Y_11, 2, 1.000000E+00
> GRID, RSET_Y_12, 9, 1.000000E+00
> GRID, RSET_Y_13, 2, 1.000000E+00
> GRID, RSET_Y_14, 11, 1.000000E+00
> GRID, RSET_Y_15, 9, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> OBJ, NAME, INLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 3.100000E-01, 1.000000E-01
> OBJ, GEOMETRY, cube3t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.
> OBJ, VELOCITY, 1.0E-02 ,0. ,0.
> OBJ, TEMPERATURE, 18.
> OBJ, NAME, OUTLET
> OBJ, POSITION, 3.100000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 3.100000E-01, 1.000000E-01
> OBJ, GEOMETRY, cube12t
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.
> OBJ, NAME, PIPE1
> 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, TYPE, BLOCKAGE
> OBJ, MATERIAL, 111,STEEL at 27 deg c (C = 1%)
> OBJ, NAME, PIPE2
> OBJ, POSITION, 3.000000E-02, 1.000000E-01, 0.000000E+00
> OBJ, SIZE, 4.000000E-02, 4.000000E-02, 1.000000E-01
> OBJ, GEOMETRY, cylinder
> OBJ, TYPE, BLOCKAGE
> OBJ, INFSRC_TEM1, 10*YG with WHOL
> OBJ, INFINI_TEM1, INTEM
> OBJ, MATERIAL, 111,STEEL at 27 deg c (C = 1%)
> OBJ, NAME, PIPE3
> OBJ, POSITION, 3.000000E-02, 1.700000E-01, 0.000000E+00
> OBJ, SIZE, 4.000000E-02, 4.000000E-02, 1.000000E-01
> OBJ, GEOMETRY, cylinder
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 111,STEEL at 27 deg c (C = 1%)
> OBJ, NAME, PIPE4
> OBJ, POSITION, 3.000000E-02, 2.400000E-01, 0.000000E+00
> OBJ, SIZE, 4.000000E-02, 4.000000E-02, 1.000000E-01
> OBJ, GEOMETRY, cylinder
> OBJ, TYPE, BLOCKAGE
> OBJ, INFSRC_TEM1, 10*YG with WHOL
> OBJ, INFINI_TEM1, INTEM
> OBJ, MATERIAL, 111,STEEL at 27 deg c (C = 1%)
> OBJ, NAME, PIPE5
> OBJ, POSITION, 9.000000E-02, 6.500000E-02, 0.000000E+00
> OBJ, SIZE, 4.000000E-02, 4.000000E-02, 1.000000E-01
> OBJ, GEOMETRY, cylinder
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 111,STEEL at 27 deg c (C = 1%)
> OBJ, NAME, PIPE6
> OBJ, POSITION, 9.000000E-02, 1.350000E-01, 0.000000E+00
> OBJ, SIZE, 4.000000E-02, 4.000000E-02, 1.000000E-01
> OBJ, GEOMETRY, cylinder
> OBJ, TYPE, BLOCKAGE
> OBJ, INFSRC_TEM1, 10*YG with WHOL
> OBJ, INFINI_TEM1, INTEM
> OBJ, MATERIAL, 111,STEEL at 27 deg c (C = 1%)
> OBJ, NAME, PIPE7
> OBJ, POSITION, 9.000000E-02, 2.050000E-01, 0.000000E+00
> OBJ, SIZE, 4.000000E-02, 4.000000E-02, 1.000000E-01
> OBJ, GEOMETRY, cylinder
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 111,STEEL at 27 deg c (C = 1%)
STOP