TALK=T;RUN(1,1)
PHOTON USE
P
gr ou z 1
msg PIPE WITH SHIFTED THERMAL BRIDGES
msg Temperature distribution at the pipe outlet
con temp z 1 fil;.00001
msg Press Enter to continue
pause;con off;red
msg Angle between Velocity vector and Z-axis
con angz z 1 fil;.00001
msg Press E to end
ENDUSE
DISPLAY
The thermal bridge is the internal longitudinal fin with
the hight equals the pipe diameter. It acts as an extra wall
providing the extra heat transfer by its both sides. In this
case the advanced fin configuration is examined. It provides
the same surface area as single thermal bridge in impeller-
like manner of surface interruption and shifting it around
axes. The user can choose the bridge type, the flow regime,
Reynolds number and type of the solution procedure in
interactive mode. The PLANT is used to calculate the angle
between the velocity vector and Z-axis to asses the validity
of the parabolic option.
ENDDIS
PLANTBEGIN
* Angle between velocity vector and Z-axis
ANGZ=ACOS(W1/sqrt(U1**2+V1**2+W1**2))
PLANTEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(BRIDGED PIPE HEAT TRANSFER:114 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 14
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 10
XULAST =1.57
XFRAC(1)=0.032413 ;XFRAC(2)=0.105311
XFRAC(3)=0.20981 ;XFRAC(4)=0.342155
XFRAC(5)=0.5 ;XFRAC(6)=0.657845
XFRAC(7)=0.790189 ;XFRAC(8)=0.894689
XFRAC(9)=0.967587 ;XFRAC(10)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 10
YVLAST =4.1605E-03
YFRAC(1)=0.032413 ;YFRAC(2)=0.105311
YFRAC(3)=0.20981 ;YFRAC(4)=0.342155
YFRAC(5)=0.5 ;YFRAC(6)=0.657845
YFRAC(7)=0.790189 ;YFRAC(8)=0.894689
YFRAC(9)=0.967587 ;YFRAC(10)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 100
ZWLAST =0.3
ZFRAC(1)=3.241313E-03 ;ZFRAC(2)=0.010531
ZFRAC(3)=0.020981 ;ZFRAC(4)=0.034216
ZFRAC(5)=0.05 ;ZFRAC(6)=0.065784
ZFRAC(7)=0.079019 ;ZFRAC(8)=0.089469
ZFRAC(9)=0.096759 ;ZFRAC(10)=0.1
ZFRAC(11)=0.103241 ;ZFRAC(12)=0.110531
ZFRAC(13)=0.120981 ;ZFRAC(14)=0.134216
ZFRAC(15)=0.15 ;ZFRAC(16)=0.165784
ZFRAC(17)=0.179019 ;ZFRAC(18)=0.189469
ZFRAC(19)=0.196759 ;ZFRAC(20)=0.2
ZFRAC(21)=0.203241 ;ZFRAC(22)=0.210531
ZFRAC(23)=0.220981 ;ZFRAC(24)=0.234216
ZFRAC(25)=0.25 ;ZFRAC(26)=0.265785
ZFRAC(27)=0.279019 ;ZFRAC(28)=0.289469
ZFRAC(29)=0.296759 ;ZFRAC(30)=0.3
ZFRAC(31)=0.303241 ;ZFRAC(32)=0.310531
ZFRAC(33)=0.320981 ;ZFRAC(34)=0.334215
ZFRAC(35)=0.35 ;ZFRAC(36)=0.365784
ZFRAC(37)=0.379019 ;ZFRAC(38)=0.389469
ZFRAC(39)=0.396759 ;ZFRAC(40)=0.4
ZFRAC(41)=0.403241 ;ZFRAC(42)=0.410531
ZFRAC(43)=0.420981 ;ZFRAC(44)=0.434216
ZFRAC(45)=0.45 ;ZFRAC(46)=0.465784
ZFRAC(47)=0.479019 ;ZFRAC(48)=0.489469
ZFRAC(49)=0.496759 ;ZFRAC(50)=0.5
ZFRAC(51)=0.503241 ;ZFRAC(52)=0.510531
ZFRAC(53)=0.520981 ;ZFRAC(54)=0.534216
ZFRAC(55)=0.55 ;ZFRAC(56)=0.565785
ZFRAC(57)=0.579019 ;ZFRAC(58)=0.589469
ZFRAC(59)=0.596759 ;ZFRAC(60)=0.6
ZFRAC(61)=0.603241 ;ZFRAC(62)=0.610531
ZFRAC(63)=0.620981 ;ZFRAC(64)=0.634216
ZFRAC(65)=0.65 ;ZFRAC(66)=0.665785
ZFRAC(67)=0.679019 ;ZFRAC(68)=0.689469
ZFRAC(69)=0.696759 ;ZFRAC(70)=0.7
ZFRAC(71)=0.703241 ;ZFRAC(72)=0.710531
ZFRAC(73)=0.720981 ;ZFRAC(74)=0.734215
ZFRAC(75)=0.75 ;ZFRAC(76)=0.765785
ZFRAC(77)=0.779019 ;ZFRAC(78)=0.789469
ZFRAC(79)=0.796759 ;ZFRAC(80)=0.8
ZFRAC(81)=0.803241 ;ZFRAC(82)=0.810531
ZFRAC(83)=0.820981 ;ZFRAC(84)=0.834216
ZFRAC(85)=0.85 ;ZFRAC(86)=0.865785
ZFRAC(87)=0.879019 ;ZFRAC(88)=0.889469
ZFRAC(89)=0.896759 ;ZFRAC(90)=0.9
ZFRAC(91)=0.901995 ;ZFRAC(92)=0.906483
ZFRAC(93)=0.912915 ;ZFRAC(94)=0.921062
ZFRAC(95)=0.930779 ;ZFRAC(96)=0.941962
ZFRAC(97)=0.954534 ;ZFRAC(98)=0.968431
ZFRAC(99)=0.983601 ;ZFRAC(100)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(7)=W1
NAME(12)=KE ;NAME(13)=EP
NAME(14)=TEMP ;NAME(149)=ENUT
NAME(150)=ANGZ
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,N,N,N,Y)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(TEMP,Y,Y,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
SOLUTN(ANGZ,Y,N,N,N,N,Y)
VIST = 149
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(U1,Y,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(TEMP,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWENT = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1.22 ;TMP1 =0. ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =1.465E-05 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(W1)=1. ;PRNDTL(KE)=1.
PRNDTL(EP)=1. ;PRNDTL(TEMP)=0.736
PRT(U1)=1. ;PRT(V1)=1.
PRT(W1)=1. ;PRT(KE)=1.
PRT(EP)=1.314 ;PRT(TEMP)=0.86
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
FIINIT(V1)=1.0E-10 ;FIINIT(W1)=0.880303
FIINIT(KE)=3.487198E-03 ;FIINIT(EP)=0.090358
FIINIT(TEMP)=0.5 ;FIINIT(ENUT)=1.0E-10
FIINIT(ANGZ)=1.0E-10
No PATCHes yet used for this Group
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(INL ,LOW , 1, 10, 1, 10, 1, 1, 1, 1)
COVAL(INL ,P1 , FIXFLU ,1.073969 )
COVAL(INL ,U1 ,0. ,0. )
COVAL(INL ,V1 ,0. ,0. )
COVAL(INL ,W1 ,0. ,0.880303 )
COVAL(INL ,KE ,0. ,3.487198E-03 )
COVAL(INL ,EP ,0. ,0.090358 )
COVAL(INL ,TEMP,0. ,0. )
PATCH(EXIT ,HIGH , 1, 10, 1, 10, 100, 100, 1, 1)
COVAL(EXIT ,P1 ,1000. ,0. )
COVAL(EXIT ,U1 ,0. ,0. )
COVAL(EXIT ,V1 ,0. ,0. )
COVAL(EXIT ,W1 ,0. ,0. )
COVAL(EXIT ,KE ,0. , SAME )
COVAL(EXIT ,EP ,0. , SAME )
COVAL(EXIT ,TEMP,0. , SAME )
PATCH(NWAL ,NWALL , 1, 10, 10, 10, 1, 100, 1, 1)
COVAL(NWAL ,U1 , GRND2 ,0. )
COVAL(NWAL ,W1 , GRND2 ,0. )
COVAL(NWAL ,KE , GRND2 , GRND2 )
COVAL(NWAL ,EP , GRND2 , GRND2 )
COVAL(NWAL ,TEMP, GRND2 ,1. )
PATCH(1BRIDGE ,WWALL , 1, 1, 1, 10, 1, 100, 1, 1)
COVAL(1BRIDGE ,V1 , GRND2 ,0. )
COVAL(1BRIDGE ,W1 , GRND2 ,0. )
COVAL(1BRIDGE ,KE , GRND2 , GRND2 )
COVAL(1BRIDGE ,EP , GRND2 , GRND2 )
COVAL(1BRIDGE ,TEMP, GRND2 ,1. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 250 ;ISWC1 = 1
LITHYD = 15 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=10
LITER(V1)=10 ;LITER(W1)=10
LITER(KE)=20 ;LITER(EP)=20
LITER(TEMP)=20
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
ENDIT(TEMP)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,0.3)
RELAX(U1,FALSDT,1.0E-03)
RELAX(V1,FALSDT,1.0E-03)
RELAX(W1,FALSDT,1.0E-03)
RELAX(KE,FALSDT,1.0E-05)
RELAX(EP,FALSDT,1.0E-05)
RELAX(TEMP,FALSDT,1.0E+05)
RELAX(ENUT,LINRLX,0.5)
RELAX(ANGZ,LINRLX,1.)
KELIN = 0
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(TEMP)=1.0E+10 ;VARMIN(TEMP)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
VARMAX(ANGZ)=1.0E+10 ;VARMIN(ANGZ)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
NAMSAT =MOSG
GENK = T
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
RG( 1) =1.460066E-05
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(U1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(TEMP,Y,N,Y,Y,Y,Y)
OUTPUT(ENUT,Y,N,Y,N,N,N)
OUTPUT(ANGZ,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 3 ;IYMON = 3 ;IZMON = 25
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP