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
* GROUP 1. Run identifiers and other preliminaries.
TEXT(BRIDGED PIPE HEAT TRANSFER:114
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
REAL(WINL,EPS,RE,DIAM,FRIC);EPS=0.001
CHAR(TURLAM,ELLPAR,BRIDGES)
MESG( **** EFFICIENCY OF HEAT TRANSFER AUGMENTATION DEVICES ****
MESG( **** Case study: PIPE WITH THERMAL BRIDGES ***
MESG( **** User defined parametric study variables ***
MESG(
MESG(CONTINOUS BRIDGE ?(Y/N)
READVDU(BRIDGES,CHAR,Y)
CASE :BRIDGES: OF
WHEN Y,1
MESG(CONTINOUS BRIDGE
WHEN N,1
MESG(SHIFTED BRIDGE
ENDCASE
MESG(TURBULENT FLOW ?(Y/N)
READVDU(TURLAM,CHAR,Y)
CASE :TURLAM: OF
WHEN Y,1
MESG(TURBULENT FLOW
FRIC=0.018
MESG(INLET TURBULENCE PARAMETER (FRIC) ? (:FRIC: at RE=1.E5)
READVDU(FRIC,REAL,FRIC)
MESG(FRIC = :FRIC:
WHEN N,1
MESG(LAMINAR FLOW
ENDCASE
RE=500
MESG(Re= Win*Diam/ENUL ? (:RE:)
READVDU(RE,REAL,RE)
MESG(RE = :RE:
MESG(ELLIPTIC PROCEDURE ?(Y/N)
READVDU(ELLPAR,CHAR,Y)
CASE :ELLPAR: OF
WHEN Y,1
MESG(ELLIPTIC FLOW
WHEN N,1
MESG(PARABOLIC FLOW
ENDCASE
MESG(
MESG( **** CONCLUDING NOTE:
MESG( **** To specify the rest of your input data edit Q1 file
* GROUP 3. x-direction grid specification.
CARTES=F
NREGX=1;REGEXT(X,1.5700E+00)
IREGX=1;GRDPWR(X,-10,1.5700E+00,1.7)
* GROUP 4. y-direction grid specification.
DIAM=2.*4.1605E-3
NREGY=1;REGEXT(Y,1)
IREGY=1;GRDPWR(Y,-10,DIAM/2.,1.7)
* GROUP 5. z-direction grid specification.
CASE :ELLPAR: OF
WHEN N,1
PARAB=T
ENDCASE
*
NREGZ=10;REGEXT(Z,3.0000E-02)
IREGZ=1;GRDPWR(Z,-10,3.0000E-02,1.7000E+00)
IREGZ=2;GRDPWR(Z,-10,3.0000E-02,1.7000E+00)
IREGZ=3;GRDPWR(Z,-10,3.0000E-02,1.7000E+00)
IREGZ=4;GRDPWR(Z,-10,3.0000E-02,1.7000E+00)
IREGZ=5;GRDPWR(Z,-10,3.0000E-02,1.7000E+00)
IREGZ=6;GRDPWR(Z,-10,3.0000E-02,1.7000E+00)
IREGZ=7;GRDPWR(Z,-10,3.0000E-02,1.7000E+00)
IREGZ=8;GRDPWR(Z,-10,3.0000E-02,1.7000E+00)
IREGZ=9;GRDPWR(Z,-10,3.0000E-02,1.7000E+00)
IREGZ=10;GRDPWR(Z,10,3.0000E-02,1.7000E+00)
* GROUP 7. Variables (including porosities) named,
stored & solved.
SOLVE(P1,U1,V1,W1,H1)
NAME(H1)=TEMP
* Activate the angle storage
STORE(ANGZ)
CASE :TURLAM: OF
WHEN Y,1
TURMOD(KEMODL)
STORE(ENUT)
ENDCASE
* GROUP 8. Terms (in differential equations) and devices.
TERMS(TEMP,N,Y,Y,Y,Y,Y)
TERMS(ANGZ,N,N,N,N,Y,N)
* GROUP 9. Properties of the medium (or media).
RHO1=1.2200E+00
PRNDTL(TEMP)=7.3600E-01; PRT(TEMP)=0.86
ENUL=1.4650E-05
* GROUP 11. Initialization of fields of variables,
* porosities, etc.
WINL=RE*ENUL/DIAM
FIINIT(W1)=WINL;FIINIT(TEMP)=0.5
CASE :TURLAM: OF
WHEN Y,1
REAL(TKEIN,EPIN,GMIXL)
** TKEIN=0.25*WIN*WIN*FRIC where FRIC=0.018 at REYNO=1.E5
TKEIN=0.25*WINL*WINL*FRIC;FIINIT(KE)=TKEIN
**EPIN=0.1643*KIN**1.5/LMIX where LMIX=0.09*RADIUS
GMIXL=0.09*YVLAST;EPIN=0.1643*TKEIN**1.5/GMIXL
FIINIT(EP)=EPIN
ENDCASE
* GROUP 13. Boundary conditions and special sources
INLET(INL,LOW,#1,#NREGX,#1,#NREGY,#1,#1,#1,#NREGT)
VALUE(INL,P1,WINL*RHO1)
VALUE(INL,U1,0)
VALUE(INL,V1,0)
VALUE(INL,W1,WINL)
VALUE(INL,TEMP,0.)
CASE :TURLAM: OF
WHEN Y,1
VALUE(INL,KE,TKEIN)
VALUE(INL,EP,EPIN)
ENDCASE
CASE :ELLPAR: OF
WHEN Y,1
* OUTLET boundary condition, name EXIT
OUTLET(EXIT,HIGH,#1,#NREGX,#1,#NREGY,#10,#10,#1,#NREGT)
ENDCASE
* WALL boundary condition, name NWAL
PATCH(NWAL,NWALL,#1,#NREGX,#1,#1,#1,#10,#1,#NREGT)
COVAL(NWAL,U1,LOGLAW,0.0)
COVAL(NWAL,W1,LOGLAW,0.0)
CASE :TURLAM: OF
WHEN Y,1
COVAL(NWAL,KE,LOGLAW,LOGLAW)
COVAL(NWAL,EP,LOGLAW,LOGLAW)
COVAL(NWAL,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(NWAL,TEMP,1/PRNDTL(TEMP),1)
ENDCASE
CASE :BRIDGES: OF
WHEN Y,1
* PLATE boundary condition, name 1BRIDGE
PATCH(1BRIDGE,WWALL,#1,#1,#1,#1,#1,#10,1,1)
COVAL(1BRIDGE,V1,LOGLAW,0.)
COVAL(1BRIDGE,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(1BRIDGE,KE,LOGLAW,LOGLAW)
COVAL(1BRIDGE,EP,LOGLAW,LOGLAW)
COVAL(1BRIDGE,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(1BRIDGE,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
WHEN N,1
* PLATE boundary condition, name PH1
PATCH(PH1,WWALL,#1,#1,#1,#1,#1,#1,1,1)
COVAL(PH1,V1,LOGLAW,0.)
COVAL(PH1,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(PH1,KE,LOGLAW,LOGLAW)
COVAL(PH1,EP,LOGLAW,LOGLAW)
COVAL(PH1,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(PH1,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
* PLATE boundary condition, name PV1
PATCH(PV1,EWALL,#1,#1,#1,#1,#2,#2,1,1)
COVAL(PV1,V1,LOGLAW,0.)
COVAL(PV1,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(PV1,KE,LOGLAW,LOGLAW)
COVAL(PV1,EP,LOGLAW,LOGLAW)
COVAL(PV1,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(PV1,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
* PLATE boundary condition, name PH2
PATCH(PH2,WWALL,#1,#1,#1,#1,#3,#3,1,1)
COVAL(PH2,V1,LOGLAW,0.)
COVAL(PH2,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(PH2,KE,LOGLAW,LOGLAW)
COVAL(PH2,EP,LOGLAW,LOGLAW)
COVAL(PH2,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(PH2,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
* PLATE boundary condition, name PV2
PATCH(PV2,EWALL,#1,#1,#1,#1,#4,#4,1,1)
COVAL(PV2,V1,LOGLAW,0.)
COVAL(PV2,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(PV2,KE,LOGLAW,LOGLAW)
COVAL(PV2,EP,LOGLAW,LOGLAW)
COVAL(PV2,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(PV2,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
* PLATE boundary condition, name PH3
PATCH(PH3,WWALL,#1,#1,#1,#1,#5,#5,1,1)
COVAL(PH3,V1,LOGLAW,0.)
COVAL(PH3,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(PH3,KE,LOGLAW,LOGLAW)
COVAL(PH3,EP,LOGLAW,LOGLAW)
COVAL(PH3,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(PH3,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
* PLATE boundary condition, name PV3
PATCH(PV3,EWALL,#1,#1,#1,#1,#6,#6,1,1)
COVAL(PV3,V1,LOGLAW,0.)
COVAL(PV3,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(PV3,KE,LOGLAW,LOGLAW)
COVAL(PV3,EP,LOGLAW,LOGLAW)
COVAL(PV3,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(PV3,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
* PLATE boundary condition, name PH4
PATCH(PH4,WWALL,#1,#1,#1,#1,#7,#7,1,1)
COVAL(PH4,V1,LOGLAW,0.)
COVAL(PH4,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(PH4,KE,LOGLAW,LOGLAW)
COVAL(PH4,EP,LOGLAW,LOGLAW)
COVAL(PH4,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(PH4,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
* PLATE boundary condition, name PV4
PATCH(PV4,EWALL,#1,#1,#1,#1,#8,#8,1,1)
COVAL(PV4,V1,LOGLAW,0.)
COVAL(PV4,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(PV4,KE,LOGLAW,LOGLAW)
COVAL(PV4,EP,LOGLAW,LOGLAW)
COVAL(PV4,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(PV4,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
* PLATE boundary condition, name PH5
PATCH(PH5,WWALL,#1,#1,#1,#1,#9,#9,1,1)
COVAL(PH5,V1,LOGLAW,0.)
COVAL(PH5,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(PH5,KE,LOGLAW,LOGLAW)
COVAL(PH5,EP,LOGLAW,LOGLAW)
COVAL(PH5,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(PH5,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
* PLATE boundary condition, name PV5
PATCH(PV5,EWALL,#1,#1,#1,#1,#10,#10,1,1)
COVAL(PV5,V1,LOGLAW,0.)
COVAL(PV5,W1,LOGLAW,0.)
CASE :TURLAM: OF
WHEN Y,1
COVAL(PV5,KE,LOGLAW,LOGLAW)
COVAL(PV5,EP,LOGLAW,LOGLAW)
COVAL(PV5,TEMP,LOGLAW,1.)
WHEN N,1
COVAL(PV5,TEMP,1./PRNDTL(TEMP),1.)
ENDCASE
ENDCASE
* GROUP 15. Termination criteria for sweeps and
LITHYD=15;LSWEEP=250
* GROUP 16. Termination criteria for inner iterations.
REAL(FAREA);FAREA=3.14159*YVLAST**2/4.
RG(1)=RHO1*WINL*FAREA
RESREF(P1)=RG(1)*EPS
RESREF(W1)=RHO1*WINL**2*FAREA*EPS
RESREF(U1)=0.01*RESREF(W1)
RESREF(V1)=0.01*RESREF(W1)
RESREF(TEMP)=RHO1*.001*WINL*FAREA*EPS
CASE :TURLAM: OF
WHEN Y,1
RESREF(KE)=RHO1*TKEIN*WINL*FAREA*EPS
RESREF(EP)=RHO1*EPIN*WINL*FAREA*EPS
ENDCASE
* GROUP 17. Under-relaxation and related devices.
RELAX(P1,LINRLX,0.3)
RELAX(U1,FALSDT,1.E-3)
RELAX(V1,FALSDT,1.E-3)
RELAX(W1,FALSDT,1.E-3)
RELAX(TEMP,FALSDT,1.E5)
CASE :TURLAM: OF
WHEN Y,1
RELAX(KE,FALSDT,0.00001)
RELAX(EP,FALSDT,0.00001)
relax(enut,linrlx,0.5)
ENDCASE
* GROUP 18. Limits on variables values or increments
* to them.
* GROUP 19. Data communicated by SATELLITE to GROUND
NAMSAT=MOSG
* GROUP 21. Preliminary printout
PLANTBEGIN
* Angle between velocity vector and Z-axis
ANGZ=ACOS(W1/sqrt(U1**2+V1**2+W1**2))
PLANTEND
* GROUP 22. Location of spot-value & frequency of
* residual printout.
IXMON=3; IYMON=3; IZMON=25
* GROUP 23. Variable-by-variable field printout and plot
* and/or tabulation of spot-values and residuals.
tstswp=-1
dmpstk=t