TALK=T;RUN(1,1)
   PHOTON USE
   p



   vec x 1 sh
   con temp x 1 sh;int 30
   msg             Velocity vectors and temperaturee
   msg  Press Enter to continue
   ENDUSE
  DISPLAY
  The flow  in  a  thick-walled  pipe  is  considered.  In order to
  account for the conductive heat transfer through  the  wall,  the
  velocities are fixed to zero in the solid portions of the domain.
  A constant heat flux is prescribed at the outer wall, and uniform
  temperature  and  velocity  profiles  are  prescribed at the pipe
  inlet.  The problem is similar to the one from  PHOENICS  library
  but  no  porosity  technique  but  PLANT  facilities  are used to
  simulate the conjugate heat transfer.

    The locally-defined variables are:
      WIN           Inlet velocity.
      QDOT          Heat flux at outer wall
      CP            Specific heat of water
      IYMET         First cell inside solid
      IYWAT         Last cell in the water
      CONDM         Metal conductivity
      CONDW         Water conductivity
  ENDDIS
     PLANTBEGIN
   VISL=VISC
   LAMPR(TEMP)=DIFS
     PLANTEND
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(CONJUGATE HEAT TRANSFER IN THICK-WALLED )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 123
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 1
 XULAST =1.
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 30
 YVLAST =0.024
   Method of pairs used for grid setting.
 YFRAC(1)=-20. ;YFRAC(2)=0.04166
 YFRAC(3)=10. ;YFRAC(4)=0.016666
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 20
 ZWLAST =0.4
 ZFRAC(1)=2.5E-03 ;ZFRAC(2)=1.0E-02
 ZFRAC(3)=0.0225 ;ZFRAC(4)=0.04
 ZFRAC(5)=0.0625 ;ZFRAC(6)=0.09
 ZFRAC(7)=0.1225 ;ZFRAC(8)=0.16
 ZFRAC(9)=0.2025 ;ZFRAC(10)=0.25
 ZFRAC(11)=0.3025 ;ZFRAC(12)=0.36
 ZFRAC(13)=0.4225 ;ZFRAC(14)=0.49
 ZFRAC(15)=0.5625 ;ZFRAC(16)=0.64
 ZFRAC(17)=0.7225 ;ZFRAC(18)=0.81
 ZFRAC(19)=0.9025 ;ZFRAC(20)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(5)=V1
 NAME(7)=W1 ;NAME(148)=DIFS
 NAME(149)=VISC ;NAME(150)=TEMP
    * 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(V1,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(DIFS,Y,N,N,N,N,Y)
 SOLUTN(VISC,Y,N,N,N,N,Y)
 SOLUTN(TEMP,Y,Y,N,N,N,Y)
 ************************************************************
  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(V1,Y,Y,Y,Y,Y,Y)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 TERMS(TEMP,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWENL = 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 =1000. ;TMP1 =0.
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL = GRND ;ENUT =0.
 ENULA =0. ;ENULB =0.
 ENULC =0.
 PRNDTL(V1)=1. ;PRNDTL(W1)=1.
 PRNDTL(TEMP)=-GRND
 PRT(V1)=1. ;PRT(W1)=1.
 PRT(TEMP)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(V1)=1.0E-10
 FIINIT(W1)=1.0E-02 ;FIINIT(DIFS)=1.0E-10
 FIINIT(VISC)=1.0E-10 ;FIINIT(TEMP)=1.0E-10
 
 PATCH(VISCW ,INIVAL, 1, 1, 1, 20, 1, 20, 1, 1)
 INIT(VISCW ,VISC,0. ,1.0E-06 )
 
 PATCH(VISCM ,INIVAL, 1, 1, 21, 30, 1, 20, 1, 1)
 INIT(VISCM ,VISC,0. ,1.0E+10 )
 
 PATCH(KOVRW ,INIVAL, 1, 1, 1, 20, 1, 20, 1, 1)
 INIT(KOVRW ,DIFS,0. ,1.434035E-07 )
 
 PATCH(KOVRM ,INIVAL, 1, 1, 21, 30, 1, 20, 1, 1)
 INIT(KOVRM ,DIFS,0. ,8.652008E-06 )
 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(INLET ,LOW , 1, 1, 1, 20, 1, 1, 1, 1)
 COVAL(INLET ,P1 , FIXFLU ,10. )
 COVAL(INLET ,W1 ,0. ,1.0E-02 )
 COVAL(INLET ,TEMP,0. ,0. )
 
 PATCH(SOLIDW ,CELL , 1, 1, 21, 30, 1, 20, 1, 1)
 COVAL(SOLIDW ,W1 , FIXVAL ,0. )
 
 PATCH(SOLIDV ,CELL , 1, 1, 20, 30, 1, 20, 1, 1)
 COVAL(SOLIDV ,V1 , FIXVAL ,0. )
 
 PATCH(OUTER ,NORTH , 1, 1, 30, 30, 1, 20, 1, 1)
 COVAL(OUTER ,TEMP, FIXFLU ,0.239006 )
 
 PATCH(EXIT ,HIGH , 1, 1, 1, 20, 20, 20, 1, 1)
 COVAL(EXIT ,P1 ,1000. ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 120 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(V1)=10
 LITER(W1)=10 ;LITER(TEMP)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
 ENDIT(W1)=1.0E-03 ;ENDIT(TEMP)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(V1,FALSDT,40.)
 RELAX(W1,FALSDT,40.)
 RELAX(DIFS,LINRLX,1.)
 RELAX(VISC,LINRLX,1.)
 RELAX(TEMP,FALSDT,1.0E+09)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(DIFS)=1.0E+10 ;VARMIN(DIFS)=-1.0E+10
 VARMAX(VISC)=1.0E+10 ;VARMIN(VISC)=-1.0E+10
 VARMAX(TEMP)=1.0E+10 ;VARMIN(TEMP)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 PARSOL = F
 ISG62 = 1
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=5.036E+07 ;EX(V1)=7.972E-05
 EX(W1)=7.603E-03 ;EX(DIFS)=2.98E-06
 EX(VISC)=3.333E+09 ;EX(TEMP)=1.913
 ************************************************************
  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(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(DIFS,Y,N,Y,N,N,N)
 OUTPUT(VISC,Y,N,Y,N,N,N)
 OUTPUT(TEMP,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 6 ;IZMON = 19
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 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