Conjugate heat transfer of rotating obj


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



  vi z
  GR Y 1;GR Y M;GR Y 10 10 Z 1 1;GR Y 14 14 Z 1 1
  GR X 4 4 Y 1 9 Z 1 1
  GR X 11 11 Y 1 9 Z 1 1
  GR X 18 18 Y 1 9 Z 1 1
  GR X 7 7 Y 14 21 Z 1 1
  GR X 14 14 Y 14 21 Z 1 1
  GR X 21 21 Y 14 21 Z 1 1
  msg  Vectors
  vec z 1 sh
  msg  Hit Enter to continue
  pause;cl
  msg  Temperature contours
  con temp z 1 fil;.001
  GR Y 1;GR Y M;GR Y 10 10 Z 1 1;GR Y 14 14 Z 1 1
  GR X 4 4 Y 1 9 Z 1 1
  GR X 11 11 Y 1 9 Z 1 1;GR X 18 18 Y 1 9 Z 1 1
  GR X 7 7 Y 14 21 Z 1 1;GR X 14 14 Y 14 21 Z 1 1
  GR X 21 21 Y 14 21 Z 1 1
  msg  Hit Enter to continue
  ENDUSE
  DISPLAY
  This run   simulates  the  flow  and  and  conjugate  heat
  transfer in a cross section of concentric pipes  sectioned
  by radial fins. The fluid is sent in motion by rotation of
  the middle thick wall pipe,  the thermal  conductivity  of
  which  is  3  times higher than surrounding fluids.  It is
  heated up by volumetric heat source,  while the inner  and
  outer  walls  of  the  domain are hold under constant cold
  temperature.

  The example demonstrates how to use PLANT for  calculation
  of conjugate heat transfer in moving objects.

  PHOTON use commands are supplied.
  ENDDIS
       PLANTBEGIN
     VELAD=-3.*RG2D
    REGION() 1
    The above statement removes the rotation  velocity  from
    the solid.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
     VISL=0.001
     LAMPR(U1)=1.e5
    REGION() 1
     LAMPR(V1)=1.e5
    REGION() 1
     LAMPR(TEMP)=.003
     LAMPR(TEMP)=204.
    REGION() 1
    The above  settings  introduce   the   artificial   high
    viscosities  and  actual  conductivity  for  solid  pipe
    material marked by MARK=1.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
  ** Rotated pipe  velocity
PATCH(SS001,VOLUME,1,NX,1,NY,1,NZ,1,lstep)
   CO=1.e10
   VAL=3.0*RG2D
COVAL(SS001,U1,GRND,GRND)
   CO=1.e10
   VAL=0.0
COVAL(SS001,V1,GRND,GRND)
   CO=FIXFLU
   VAL=3.0/FIXFLU
COVAL(SS001,TEMP,GRND,GRND)
    PATCH name  SS001  indicates  that it is the solid pipe,
    MARK=1,  for which the velocities are set to solid  body
    rotation. It is heated up by volumetric heat source.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
    PLANTEND
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT( Conjugate heat transfer of rotating obj)
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 614
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 21
 XULAST =6.28
 XFRAC(1)=0.047619 ;XFRAC(2)=0.095238
 XFRAC(3)=0.142857 ;XFRAC(4)=0.190476
 XFRAC(5)=0.238095 ;XFRAC(6)=0.285714
 XFRAC(7)=0.333333 ;XFRAC(8)=0.380952
 XFRAC(9)=0.428571 ;XFRAC(10)=0.47619
 XFRAC(11)=0.52381 ;XFRAC(12)=0.571429
 XFRAC(13)=0.619048 ;XFRAC(14)=0.666667
 XFRAC(15)=0.714286 ;XFRAC(16)=0.761905
 XFRAC(17)=0.809524 ;XFRAC(18)=0.857143
 XFRAC(19)=0.904762 ;XFRAC(20)=0.952381
 XFRAC(21)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 20
 YVLAST =0.25
 RINNER =0.05 ;SNALFA =0.
 YFRAC(1)=0.05 ;YFRAC(2)=0.1
 YFRAC(3)=0.15 ;YFRAC(4)=0.2
 YFRAC(5)=0.25 ;YFRAC(6)=0.3
 YFRAC(7)=0.35 ;YFRAC(8)=0.4
 YFRAC(9)=0.45 ;YFRAC(10)=0.5
 YFRAC(11)=0.55 ;YFRAC(12)=0.6
 YFRAC(13)=0.65 ;YFRAC(14)=0.7
 YFRAC(15)=0.75 ;YFRAC(16)=0.8
 YFRAC(17)=0.85 ;YFRAC(18)=0.9
 YFRAC(19)=0.95 ;YFRAC(20)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 1
 ZWLAST =1.
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
    * X-cyclic boundaries switched
 XCYIZ( 1, T )
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(148)=EPOR
 NAME(149)=MARK ;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(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(EPOR,Y,N,N,N,N,N)
 SOLUTN(MARK,Y,N,N,N,N,Y)
 SOLUTN(TEMP,Y,Y,N,N,N,Y)
 EPOR = 148 ;HPOR = 0 ;NPOR = 0 ;VPOR = 0
 ************************************************************
  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(TEMP,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWENL = T
 U1AD = GRND ;V1AD =0. ;W1AD =0.
 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. ;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. ;ENULD =0.
 ENULE =0. ;ENULF =0. ;ENULG =0.
 IENULA = 0 ;IENULB = 0
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(TEMP)=-GRND
 PRT(U1)=1. ;PRT(V1)=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(U1)=1.0E-10
 FIINIT(V1)=1.0E-10 ;FIINIT(EPOR)=1.
 FIINIT(MARK)=0. ;FIINIT(TEMP)=1.0E-10
 
 PATCH(INI1 ,INIVAL, 1, 21, 9, 12, 1, 1, 1, 1)
 INIT(INI1 ,MARK,0. ,1. )
 
 PATCH(CMP1 ,INIVAL, 2, 2, 1, 8, 1, 1, 1, 1)
 INIT(CMP1 ,U1 ,0. ,0. )
 INIT(CMP1 ,EPOR,0. ,0. )
 
 PATCH(CMP2 ,INIVAL, 9, 9, 1, 8, 1, 1, 1, 1)
 INIT(CMP2 ,U1 ,0. ,0. )
 INIT(CMP2 ,EPOR,0. ,0. )
 
 PATCH(CMP3 ,INIVAL, 16, 16, 1, 8, 1, 1, 1, 1)
 INIT(CMP3 ,U1 ,0. ,0. )
 INIT(CMP3 ,EPOR,0. ,0. )
 
 PATCH(CMP4 ,INIVAL, 5, 5, 13, 20, 1, 1, 1, 1)
 INIT(CMP4 ,U1 ,0. ,0. )
 INIT(CMP4 ,EPOR,0. ,0. )
 
 PATCH(CMP5 ,INIVAL, 12, 12, 13, 20, 1, 1, 1, 1)
 INIT(CMP5 ,U1 ,0. ,0. )
 INIT(CMP5 ,EPOR,0. ,0. )
 
 PATCH(CMP6 ,INIVAL, 19, 19, 13, 20, 1, 1, 1, 1)
 INIT(CMP6 ,U1 ,0. ,0. )
 INIT(CMP6 ,EPOR,0. ,0. )
 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(WALLS ,SWALL , 1, 21, 1, 1, 1, 1, 1, 1)
 COVAL(WALLS ,U1 ,1. ,0. )
 COVAL(WALLS ,TEMP,1. ,0. )
 
 PATCH(WALLN ,NWALL , 1, 21, 20, 20, 1, 1, 1, 1)
 COVAL(WALLN ,U1 ,1. ,0. )
 COVAL(WALLN ,TEMP,1. ,0. )
 
 PATCH(FIXPR1 ,CELL , 1, 1, 1, 1, 1, 1, 1, 1)
 COVAL(FIXPR1 ,P1 ,1. ,0. )
 COVAL(FIXPR1 ,TEMP,0. , SAME )
 
 PATCH(FIXPR2 ,CELL , 1, 1, 20, 20, 1, 1, 1, 1)
 COVAL(FIXPR2 ,P1 ,1. ,0. )
 COVAL(FIXPR2 ,TEMP,0. , SAME )
 
 PATCH(SS001 ,VOLUME, 1, 21, 1, 20, 1, 1, 1, 1)
 COVAL(SS001 ,U1 , GRND , GRND )
 COVAL(SS001 ,V1 , GRND , GRND )
 COVAL(SS001 ,TEMP, GRND , GRND )
 
 PATCH(CMP1-EW ,EWALL , 2, 2, 1, 8, 1, 1, 1, 1)
 COVAL(CMP1-EW ,V1 ,1. ,0. )
 
 PATCH(CMP1-WW ,WWALL , 3, 3, 1, 8, 1, 1, 1, 1)
 COVAL(CMP1-WW ,V1 ,1. ,0. )
 
 PATCH(CMP2-EW ,EWALL , 9, 9, 1, 8, 1, 1, 1, 1)
 COVAL(CMP2-EW ,V1 ,1. ,0. )
 
 PATCH(CMP2-WW ,WWALL , 10, 10, 1, 8, 1, 1, 1, 1)
 COVAL(CMP2-WW ,V1 ,1. ,0. )
 
 PATCH(CMP3-EW ,EWALL , 16, 16, 1, 8, 1, 1, 1, 1)
 COVAL(CMP3-EW ,V1 ,1. ,0. )
 
 PATCH(CMP3-WW ,WWALL , 17, 17, 1, 8, 1, 1, 1, 1)
 COVAL(CMP3-WW ,V1 ,1. ,0. )
 
 PATCH(CMP4-EW ,EWALL , 5, 5, 13, 20, 1, 1, 1, 1)
 COVAL(CMP4-EW ,V1 ,1. ,0. )
 
 PATCH(CMP4-WW ,WWALL , 6, 6, 13, 20, 1, 1, 1, 1)
 COVAL(CMP4-WW ,V1 ,1. ,0. )
 
 PATCH(CMP5-EW ,EWALL , 12, 12, 13, 20, 1, 1, 1, 1)
 COVAL(CMP5-EW ,V1 ,1. ,0. )
 
 PATCH(CMP5-WW ,WWALL , 13, 13, 13, 20, 1, 1, 1, 1)
 COVAL(CMP5-WW ,V1 ,1. ,0. )
 
 PATCH(CMP6-EW ,EWALL , 19, 19, 13, 20, 1, 1, 1, 1)
 COVAL(CMP6-EW ,V1 ,1. ,0. )
 
 PATCH(CMP6-WW ,WWALL , 20, 20, 13, 20, 1, 1, 1, 1)
 COVAL(CMP6-WW ,V1 ,1. ,0. )
 XCYCLE = T
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 200 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-02
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(TEMP)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(TEMP)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.5)
 RELAX(U1,FALSDT,100.)
 RELAX(V1,FALSDT,100.)
 RELAX(MARK,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(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10
 VARMAX(MARK)=1.0E+10 ;VARMIN(MARK)=-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
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=0.02117 ;EX(U1)=0.1775
 EX(V1)=0.02174 ;EX(EPOR)=0.8857
 EX(MARK)=0.2 ;EX(TEMP)=1.508
 ************************************************************
  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(EPOR,Y,N,Y,N,N,N)
 OUTPUT(MARK,Y,N,Y,N,N,N)
 OUTPUT(TEMP,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 1 ;IZMON = 1
 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
 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
 IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
 CSG1    ='P'
STOP