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. )
 ************************************************************
  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