TALK=T;RUN(1,1) PHOTON USE AUTOPLOT file PHI 5 cl;d 1 h1;d 1 ha;col3 1;blb4 2;redr msg temperature profile; pressto continue pause cl;d 1 radx;d 1 ra;col3 1;blb4 2;redr msg radiation-flux profile; press e to end pause;end ENDUSE DISPLAY The problem considered is that of radiative heat transfer in a 1d plane-parallel slab in radiative equilibrium. The slab is of thickness L and the gray medium may absorb, emit and scatter radiation. At x=L the wall is a diffuse emitter and reflecter kept at a fixed temperature. At the x=0 the net radiative heat flux Qrad is specified at the wall. ENDDIS ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(X-D Radiative Equilibrium In Slab ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 14 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 50 XULAST =1. XFRAC(1)=0.02 ;XFRAC(6)=0.12 XFRAC(11)=0.22 ;XFRAC(16)=0.32 XFRAC(21)=0.42 ;XFRAC(26)=0.52 XFRAC(31)=0.62 ;XFRAC(36)=0.72 XFRAC(41)=0.82 ;XFRAC(46)=0.92 ************************************************************ Group 4. Y-Direction Grid Spacing NY = 1 YVLAST =1. YFRAC(1)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 1 ZWLAST =1. ZFRAC(1)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(14)=H1 ;NAME(147)=RA NAME(148)=HA ;NAME(149)=EMPO NAME(150)=RADX * Y in SOLUTN argument list denotes: * 1-stored 2-solved 3-whole-field * 4-point-by-point 5-explicit 6-harmonic averaging SOLUTN(H1,Y,Y,N,N,N,Y) SOLUTN(RA,Y,N,N,N,N,Y) SOLUTN(HA,Y,N,N,N,N,Y) SOLUTN(EMPO,Y,N,N,N,N,Y) SOLUTN(RADX,Y,Y,N,N,N,N) ************************************************************ 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(H1,N,N,N,N,Y,Y) TERMS(RADX,N,N,Y,N,Y,N) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F 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 = GRND2 ;EL1 =0. TSURR =0. ;TEMP0 =0. ;PRESS0 =0. DVO1DT =0. ;DRH1DP =0. TMP1A =0. ;TMP1B =1. ;TMP1C =0. TMP2A =0. ;TMP2B =0. ;TMP2C =0. EMISS =0. ;SCATT =0. RADIA =0.5 ;RADIB =0.5 ENUL =1.0E-05 ;ENUT =0. PRNDTL(RADX)=-GRND1 PRT(RADX)=1.0E+10 CP1 =1. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(H1)=1000. ;FIINIT(RA)=1.0E-10 FIINIT(HA)=1.0E-10 ;FIINIT(EMPO)=1.0E-10 FIINIT(RADX)=5.6697E+04 PATCH(IN1 ,INIVAL, 1, 1, 1, 1, 1, 1, 1, 1) INIT(IN1 ,RA ,0. ,1.69524E+05 ) INIT(IN1 ,HA ,0. ,1314.975952 ) PATCH(IN2 ,INIVAL, 2, 2, 1, 1, 1, 1, 1, 1) INIT(IN2 ,RA ,0. ,1.683901E+05 ) INIT(IN2 ,HA ,0. ,1312.771362 ) PATCH(IN3 ,INIVAL, 3, 3, 1, 1, 1, 1, 1, 1) INIT(IN3 ,RA ,0. ,1.672562E+05 ) INIT(IN3 ,HA ,0. ,1310.555786 ) PATCH(IN4 ,INIVAL, 4, 4, 1, 1, 1, 1, 1, 1) INIT(IN4 ,RA ,0. ,1.661222E+05 ) INIT(IN4 ,HA ,0. ,1308.328857 ) PATCH(IN5 ,INIVAL, 5, 5, 1, 1, 1, 1, 1, 1) INIT(IN5 ,RA ,0. ,1.649883E+05 ) INIT(IN5 ,HA ,0. ,1306.090454 ) PATCH(IN6 ,INIVAL, 6, 6, 1, 1, 1, 1, 1, 1) INIT(IN6 ,RA ,0. ,1.638543E+05 ) INIT(IN6 ,HA ,0. ,1303.840454 ) PATCH(IN7 ,INIVAL, 7, 7, 1, 1, 1, 1, 1, 1) INIT(IN7 ,RA ,0. ,1.627204E+05 ) INIT(IN7 ,HA ,0. ,1301.578857 ) PATCH(IN8 ,INIVAL, 8, 8, 1, 1, 1, 1, 1, 1) INIT(IN8 ,RA ,0. ,1.615865E+05 ) INIT(IN8 ,HA ,0. ,1299.305298 ) PATCH(IN9 ,INIVAL, 9, 9, 1, 1, 1, 1, 1, 1) INIT(IN9 ,RA ,0. ,1.604525E+05 ) INIT(IN9 ,HA ,0. ,1297.019775 ) PATCH(IN10 ,INIVAL, 10, 10, 1, 1, 1, 1, 1, 1) INIT(IN10 ,RA ,0. ,1.593186E+05 ) INIT(IN10 ,HA ,0. ,1294.722168 ) PATCH(IN11 ,INIVAL, 11, 11, 1, 1, 1, 1, 1, 1) INIT(IN11 ,RA ,0. ,1.581846E+05 ) INIT(IN11 ,HA ,0. ,1292.412231 ) PATCH(IN12 ,INIVAL, 12, 12, 1, 1, 1, 1, 1, 1) INIT(IN12 ,RA ,0. ,1.570507E+05 ) INIT(IN12 ,HA ,0. ,1290.089844 ) PATCH(IN13 ,INIVAL, 13, 13, 1, 1, 1, 1, 1, 1) INIT(IN13 ,RA ,0. ,1.559168E+05 ) INIT(IN13 ,HA ,0. ,1287.754761 ) PATCH(IN14 ,INIVAL, 14, 14, 1, 1, 1, 1, 1, 1) INIT(IN14 ,RA ,0. ,1.547828E+05 ) INIT(IN14 ,HA ,0. ,1285.406982 ) PATCH(IN15 ,INIVAL, 15, 15, 1, 1, 1, 1, 1, 1) INIT(IN15 ,RA ,0. ,1.536489E+05 ) INIT(IN15 ,HA ,0. ,1283.046265 ) PATCH(IN16 ,INIVAL, 16, 16, 1, 1, 1, 1, 1, 1) INIT(IN16 ,RA ,0. ,1.525149E+05 ) INIT(IN16 ,HA ,0. ,1280.672485 ) PATCH(IN17 ,INIVAL, 17, 17, 1, 1, 1, 1, 1, 1) INIT(IN17 ,RA ,0. ,1.51381E+05 ) INIT(IN17 ,HA ,0. ,1278.2854 ) PATCH(IN18 ,INIVAL, 18, 18, 1, 1, 1, 1, 1, 1) INIT(IN18 ,RA ,0. ,1.50247E+05 ) INIT(IN18 ,HA ,0. ,1275.884766 ) PATCH(IN19 ,INIVAL, 19, 19, 1, 1, 1, 1, 1, 1) INIT(IN19 ,RA ,0. ,1.491131E+05 ) INIT(IN19 ,HA ,0. ,1273.470703 ) PATCH(IN20 ,INIVAL, 20, 20, 1, 1, 1, 1, 1, 1) INIT(IN20 ,RA ,0. ,1.479792E+05 ) INIT(IN20 ,HA ,0. ,1271.042725 ) PATCH(IN21 ,INIVAL, 21, 21, 1, 1, 1, 1, 1, 1) INIT(IN21 ,RA ,0. ,1.468452E+05 ) INIT(IN21 ,HA ,0. ,1268.600708 ) PATCH(IN22 ,INIVAL, 22, 22, 1, 1, 1, 1, 1, 1) INIT(IN22 ,RA ,0. ,1.457113E+05 ) INIT(IN22 ,HA ,0. ,1266.144531 ) PATCH(IN23 ,INIVAL, 23, 23, 1, 1, 1, 1, 1, 1) INIT(IN23 ,RA ,0. ,1.445773E+05 ) INIT(IN23 ,HA ,0. ,1263.67395 ) PATCH(IN24 ,INIVAL, 24, 24, 1, 1, 1, 1, 1, 1) INIT(IN24 ,RA ,0. ,1.434434E+05 ) INIT(IN24 ,HA ,0. ,1261.188843 ) PATCH(IN25 ,INIVAL, 25, 25, 1, 1, 1, 1, 1, 1) INIT(IN25 ,RA ,0. ,1.423095E+05 ) INIT(IN25 ,HA ,0. ,1258.688965 ) PATCH(IN26 ,INIVAL, 26, 26, 1, 1, 1, 1, 1, 1) INIT(IN26 ,RA ,0. ,1.411755E+05 ) INIT(IN26 ,HA ,0. ,1256.174072 ) PATCH(IN27 ,INIVAL, 27, 27, 1, 1, 1, 1, 1, 1) INIT(IN27 ,RA ,0. ,1.400416E+05 ) INIT(IN27 ,HA ,0. ,1253.644043 ) PATCH(IN28 ,INIVAL, 28, 28, 1, 1, 1, 1, 1, 1) INIT(IN28 ,RA ,0. ,1.389077E+05 ) INIT(IN28 ,HA ,0. ,1251.098633 ) PATCH(IN29 ,INIVAL, 29, 29, 1, 1, 1, 1, 1, 1) INIT(IN29 ,RA ,0. ,1.377737E+05 ) INIT(IN29 ,HA ,0. ,1248.537476 ) PATCH(IN30 ,INIVAL, 30, 30, 1, 1, 1, 1, 1, 1) INIT(IN30 ,RA ,0. ,1.366398E+05 ) INIT(IN30 ,HA ,0. ,1245.960449 ) PATCH(IN31 ,INIVAL, 31, 31, 1, 1, 1, 1, 1, 1) INIT(IN31 ,RA ,0. ,1.355058E+05 ) INIT(IN31 ,HA ,0. ,1243.367432 ) PATCH(IN32 ,INIVAL, 32, 32, 1, 1, 1, 1, 1, 1) INIT(IN32 ,RA ,0. ,1.343719E+05 ) INIT(IN32 ,HA ,0. ,1240.757935 ) PATCH(IN33 ,INIVAL, 33, 33, 1, 1, 1, 1, 1, 1) INIT(IN33 ,RA ,0. ,1.33238E+05 ) INIT(IN33 ,HA ,0. ,1238.13208 ) PATCH(IN34 ,INIVAL, 34, 34, 1, 1, 1, 1, 1, 1) INIT(IN34 ,RA ,0. ,1.32104E+05 ) INIT(IN34 ,HA ,0. ,1235.489258 ) PATCH(IN35 ,INIVAL, 35, 35, 1, 1, 1, 1, 1, 1) INIT(IN35 ,RA ,0. ,1.309701E+05 ) INIT(IN35 ,HA ,0. ,1232.829468 ) PATCH(IN36 ,INIVAL, 36, 36, 1, 1, 1, 1, 1, 1) INIT(IN36 ,RA ,0. ,1.298361E+05 ) INIT(IN36 ,HA ,0. ,1230.152222 ) PATCH(IN37 ,INIVAL, 37, 37, 1, 1, 1, 1, 1, 1) INIT(IN37 ,RA ,0. ,1.287022E+05 ) INIT(IN37 ,HA ,0. ,1227.45752 ) PATCH(IN38 ,INIVAL, 38, 38, 1, 1, 1, 1, 1, 1) INIT(IN38 ,RA ,0. ,1.275682E+05 ) INIT(IN38 ,HA ,0. ,1224.744873 ) PATCH(IN39 ,INIVAL, 39, 39, 1, 1, 1, 1, 1, 1) INIT(IN39 ,RA ,0. ,1.264343E+05 ) INIT(IN39 ,HA ,0. ,1222.01416 ) PATCH(IN40 ,INIVAL, 40, 40, 1, 1, 1, 1, 1, 1) INIT(IN40 ,RA ,0. ,1.253004E+05 ) INIT(IN40 ,HA ,0. ,1219.264893 ) PATCH(IN41 ,INIVAL, 41, 41, 1, 1, 1, 1, 1, 1) INIT(IN41 ,RA ,0. ,1.241664E+05 ) INIT(IN41 ,HA ,0. ,1216.496948 ) PATCH(IN42 ,INIVAL, 42, 42, 1, 1, 1, 1, 1, 1) INIT(IN42 ,RA ,0. ,1.230325E+05 ) INIT(IN42 ,HA ,0. ,1213.709961 ) PATCH(IN43 ,INIVAL, 43, 43, 1, 1, 1, 1, 1, 1) INIT(IN43 ,RA ,0. ,1.218985E+05 ) INIT(IN43 ,HA ,0. ,1210.903687 ) PATCH(IN44 ,INIVAL, 44, 44, 1, 1, 1, 1, 1, 1) INIT(IN44 ,RA ,0. ,1.207646E+05 ) INIT(IN44 ,HA ,0. ,1208.077759 ) PATCH(IN45 ,INIVAL, 45, 45, 1, 1, 1, 1, 1, 1) INIT(IN45 ,RA ,0. ,1.196307E+05 ) INIT(IN45 ,HA ,0. ,1205.231934 ) PATCH(IN46 ,INIVAL, 46, 46, 1, 1, 1, 1, 1, 1) INIT(IN46 ,RA ,0. ,1.184967E+05 ) INIT(IN46 ,HA ,0. ,1202.365723 ) PATCH(IN47 ,INIVAL, 47, 47, 1, 1, 1, 1, 1, 1) INIT(IN47 ,RA ,0. ,1.173628E+05 ) INIT(IN47 ,HA ,0. ,1199.47876 ) PATCH(IN48 ,INIVAL, 48, 48, 1, 1, 1, 1, 1, 1) INIT(IN48 ,RA ,0. ,1.162288E+05 ) INIT(IN48 ,HA ,0. ,1196.570923 ) PATCH(IN49 ,INIVAL, 49, 49, 1, 1, 1, 1, 1, 1) INIT(IN49 ,RA ,0. ,1.150949E+05 ) INIT(IN49 ,HA ,0. ,1193.641724 ) PATCH(IN50 ,INIVAL, 50, 50, 1, 1, 1, 1, 1, 1) INIT(IN50 ,RA ,0. ,1.13961E+05 ) INIT(IN50 ,HA ,0. ,1190.690796 ) 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(RADISO ,FREEVL, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(RADISO ,H1 , GRND1 , GRND1 ) COVAL(RADISO ,RADX, GRND1 , GRND1 ) PATCH(WALLRA ,WEST , 1, 1, 1, 1, 1, 1, 1, 1) COVAL(WALLRA ,RADX, FIXFLU ,5.6697E+04 ) PATCH(WALLRB ,EAST , 50, 50, 1, 1, 1, 1, 1, 1) COVAL(WALLRB ,RADX,1. ,5.6697E+04 ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 50 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations LITER(H1)=20 ;LITER(RADX)=2 ENDIT(H1)=1.0E-03 ;ENDIT(RADX)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(H1,FALSDT,8.818809E-04) RELAX(RA,LINRLX,1.) RELAX(HA,LINRLX,1.) RELAX(EMPO,LINRLX,1.) RELAX(RADX,FALSDT,1.0E+09) OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10 VARMAX(RA)=1.0E+10 ;VARMIN(RA)=-1.0E+10 VARMAX(HA)=1.0E+10 ;VARMIN(HA)=-1.0E+10 VARMAX(EMPO)=1.0E+10 ;VARMIN(EMPO)=-1.0E+10 VARMAX(RADX)=1.0E+10 ;VARMIN(RADX)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND PARSOL = F ISG62 = 1 SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ 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(H1,Y,N,N,Y,Y,Y) OUTPUT(RA,Y,N,Y,N,N,N) OUTPUT(HA,Y,N,Y,N,N,N) OUTPUT(EMPO,Y,N,Y,N,N,N) OUTPUT(RADX,Y,N,Y,Y,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 25 ;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 = 5 ;IXPRF = 1 ;IXPRL = 10000 IPLTF = 1 ;IPLTL = -1 ;NPLT = 5 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.4 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20 PATCH(XWISE ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(XWISE ,H1 ,0. ,0. ) PLOT(XWISE ,RADX,0. ,0. ) ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM STOP