TALK=T;RUN(1,1) DISPLAY The problem considered is 1d fully-developed turbulent flow and heat and mass transfer in a pipe at a Reynolds number of 5.E4. The turbulence is simulated by use of the Reynolds stress transport model (RSTM), and optionally by use of the k-eps model. The energy equation may be solved by means of the H1 or TEM1 equation, and two scalar variables are solved which should produce identical solutions to the energy equation. For heat and scalar transport with the RSTM, one of the following models may be employed: a simple gradient-diffusion model (IRSMSM=0); a generalised gradient-diffusion model (IRSMSM=1); or a full transport model (IRSMSM=2). The laminar Prandtl number is 3.0. ENDDIS ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(RSTM_1D DEVELOPED PIPE FLOW :T602) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 14 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = F NX = 1 XULAST =0.1 XFRAC(1)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 30 YVLAST =0.5 YFRAC(1)=0.063391 ;YFRAC(2)=0.110369 YFRAC(3)=0.152659 ;YFRAC(4)=0.192164 YFRAC(5)=0.229721 ;YFRAC(6)=0.265794 YFRAC(7)=0.300679 ;YFRAC(8)=0.334577 YFRAC(9)=0.367636 ;YFRAC(10)=0.399967 YFRAC(11)=0.431657 ;YFRAC(12)=0.462775 YFRAC(13)=0.493377 ;YFRAC(14)=0.523512 YFRAC(15)=0.553219 ;YFRAC(16)=0.582533 YFRAC(17)=0.611482 ;YFRAC(18)=0.640092 YFRAC(19)=0.668386 ;YFRAC(20)=0.696383 YFRAC(21)=0.724102 ;YFRAC(22)=0.751558 YFRAC(23)=0.778765 ;YFRAC(24)=0.805737 YFRAC(25)=0.832485 ;YFRAC(26)=0.85902 YFRAC(27)=0.885351 ;YFRAC(28)=0.911488 YFRAC(29)=0.937439 ;YFRAC(30)=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(5)=V1 ;NAME(7)=W1 NAME(12)=KE ;NAME(13)=EP NAME(124)=PRPS ;NAME(125)=DSDY NAME(126)=VTRS ;NAME(127)=VSC1 NAME(128)=VSC2 ;NAME(129)=U2RS NAME(130)=V2RS ;NAME(131)=W2RS NAME(132)=VWRS ;NAME(133)=LEN1 NAME(134)=ENUT ;NAME(135)=W2DK NAME(136)=V2DK ;NAME(137)=U2DK NAME(138)=VWDK ;NAME(139)=FWAL NAME(140)=EPDK ;NAME(141)=PK NAME(142)=DVW ;NAME(143)=PU2 NAME(144)=PV2 ;NAME(145)=PW2 NAME(146)=PVW ;NAME(147)=DWDY NAME(148)=SC2 ;NAME(149)=SC1 NAME(150)=TEM1 * Y in SOLUTN argument list denotes: * 1-stored 2-solved 3-whole-field * 4-point-by-point 5-explicit 6-harmonic averaging SOLUTN(V1,Y,N,N,N,N,Y) SOLUTN(W1,Y,Y,N,N,N,Y) SOLUTN(KE,Y,N,N,N,N,N) SOLUTN(EP,Y,Y,N,N,N,N) SOLUTN(PRPS,Y,N,N,N,N,Y) SOLUTN(DSDY,Y,N,N,N,N,Y) SOLUTN(VTRS,Y,Y,N,N,N,N) SOLUTN(VSC1,Y,Y,N,N,N,N) SOLUTN(VSC2,Y,Y,N,N,N,N) SOLUTN(U2RS,Y,Y,N,N,N,N) SOLUTN(V2RS,Y,Y,N,N,N,N) SOLUTN(W2RS,Y,Y,N,N,N,N) SOLUTN(VWRS,Y,Y,N,N,N,N) SOLUTN(LEN1,Y,N,N,N,N,N) SOLUTN(ENUT,Y,N,N,N,N,N) SOLUTN(W2DK,Y,N,N,N,N,Y) SOLUTN(V2DK,Y,N,N,N,N,Y) SOLUTN(U2DK,Y,N,N,N,N,Y) SOLUTN(VWDK,Y,N,N,N,N,Y) SOLUTN(FWAL,Y,N,N,N,N,Y) SOLUTN(EPDK,Y,N,N,N,N,Y) SOLUTN(PK,Y,N,N,N,N,Y) SOLUTN(DVW,Y,N,N,N,N,Y) SOLUTN(PU2,Y,N,N,N,N,Y) SOLUTN(PV2,Y,N,N,N,N,Y) SOLUTN(PW2,Y,N,N,N,N,Y) SOLUTN(PVW,Y,N,N,N,N,Y) SOLUTN(DWDY,Y,N,N,N,N,Y) SOLUTN(SC2,Y,Y,N,N,N,Y) SOLUTN(SC1,Y,Y,N,N,N,Y) SOLUTN(TEM1,Y,Y,N,N,N,Y) VIST = 134 LEN1 = 133 PRPS = 124 ************************************************************ 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(W1,Y,N,Y,Y,Y,Y) TERMS(EP,Y,N,Y,Y,N,Y) TERMS(VTRS,N,N,Y,Y,Y,Y) TERMS(VSC1,N,N,Y,Y,N,Y) TERMS(VSC2,N,N,Y,Y,Y,Y) TERMS(U2RS,N,N,Y,Y,N,Y) TERMS(V2RS,N,N,Y,Y,Y,Y) TERMS(W2RS,N,N,Y,Y,N,Y) TERMS(VWRS,N,N,Y,Y,Y,Y) TERMS(SC2,N,N,Y,Y,Y,Y) TERMS(SC1,N,N,Y,Y,N,Y) TERMS(TEM1,N,N,Y,Y,Y,Y) DIFCUT =0. ;ZDIFAC =1. GALA = F ;ADDDIF = F NEWENT = F UDIFF = T UCONNE = T USOURC = 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. ;TMP1 =0. EL1 = GRND 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 =2.0E-05 ;ENUT = GRND ENUTA =0. ;ENUTB =0. ENUTC =0. IENUTA = 0 PRNDTL(W1)=1. ;PRNDTL(EP)=1. PRNDTL(VTRS)=1. ;PRNDTL(VSC1)=1. PRNDTL(VSC2)=1. ;PRNDTL(U2RS)=1. PRNDTL(V2RS)=1. ;PRNDTL(W2RS)=1. PRNDTL(VWRS)=1. ;PRNDTL(SC2)=3. PRNDTL(SC1)=3. ;PRNDTL(TEM1)=CONDFILE PRT(W1)=1. ;PRT(EP)=1.314 PRT(VTRS)=1. ;PRT(VSC1)=1. PRT(VSC2)=1. ;PRT(U2RS)=1. PRT(V2RS)=1. ;PRT(W2RS)=1. PRT(VWRS)=1. ;PRT(SC2)=1. PRT(SC1)=1. ;PRT(TEM1)=1. CP1 =1. ;CP2 =1. * List of user-defined materials to be read by EARTH MATFLG=T;IMAT=1 * Name *Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr. *36 1. 2.E-5 1.0 6.667E-6 0 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(V1)=0. ;FIINIT(W1)=1. FIINIT(KE)=5.232592E-03 ;FIINIT(EP)=1.381975E-03 FIINIT(PRPS)=36. ;FIINIT(DSDY)=1.0E-10 FIINIT(VTRS)=1.0E-10 ;FIINIT(VSC1)=1.0E-10 FIINIT(VSC2)=1.0E-10 ;FIINIT(U2RS)=3.488394E-03 FIINIT(V2RS)=3.488394E-03 ;FIINIT(W2RS)=3.488394E-03 FIINIT(VWRS)=1.569778E-03 ;FIINIT(LEN1)=1.0E-10 FIINIT(ENUT)=1.0E-10 ;FIINIT(W2DK)=1.0E-10 FIINIT(V2DK)=1.0E-10 ;FIINIT(U2DK)=1.0E-10 FIINIT(VWDK)=1.0E-10 ;FIINIT(FWAL)=1.0E-10 FIINIT(EPDK)=1.0E-10 ;FIINIT(PK)=1.0E-10 FIINIT(DVW)=1.0E-10 ;FIINIT(PU2)=1.0E-10 FIINIT(PV2)=1.0E-10 ;FIINIT(PW2)=1.0E-10 FIINIT(PVW)=1.0E-10 ;FIINIT(DWDY)=1.0E-10 FIINIT(SC2)=5. ;FIINIT(SC1)=5. FIINIT(TEM1)=5. 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(WAL1 ,NWALL , 1, 1, 30, 30, 1, 1, 1, 1) COVAL(WAL1 ,W1 , GRND2 ,0. ) COVAL(WAL1 ,EP , GRND2 , GRND2 ) COVAL(WAL1 ,SC2 , GRND2 ,10. ) COVAL(WAL1 ,SC1 , GRND2 ,10. ) COVAL(WAL1 ,TEM1, GRND2 ,10. ) PATCH(PFOR ,VOLUME, 1, 1, 1, 30, 1, 1, 1, 1) COVAL(PFOR ,W1 , FIXFLU ,0.010465 ) PATCH(WDMPN03 ,PHASEM, 1, 1, 1, 30, 1, 1, 1, 1) COVAL(WDMPN03 ,VTRS, FIXFLU , GRND1 ) COVAL(WDMPN03 ,VSC1, FIXFLU , GRND1 ) COVAL(WDMPN03 ,VSC2, FIXFLU , GRND1 ) COVAL(WDMPN03 ,U2RS, FIXFLU , GRND1 ) COVAL(WDMPN03 ,V2RS, FIXFLU , GRND1 ) COVAL(WDMPN03 ,W2RS, FIXFLU , GRND1 ) COVAL(WDMPN03 ,VWRS, FIXFLU , GRND1 ) PATCH(SMPLS ,SOUTH , 1, 1, 1, 1, 1, 1, 1, 1) COVAL(SMPLS ,VTRS, GRND1 ,0. ) COVAL(SMPLS ,VSC1, GRND1 ,0. ) COVAL(SMPLS ,VSC2, GRND1 ,0. ) COVAL(SMPLS ,VWRS, GRND1 ,0. ) PATCH(FDFCWT ,PHASEM, 1, 1, 1, 30, 1, 1, 1, 1) COVAL(FDFCWT ,SC2 ,0.027333 ,10. ) COVAL(FDFCWT ,SC1 ,0.027333 ,10. ) COVAL(FDFCWT ,TEM1,0.027333 ,10. ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 300 ;ISWC1 = 1 LITHYD = 8 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations LITER(W1)=10 ;LITER(EP)=20 LITER(VTRS)=20 ;LITER(VSC1)=20 LITER(VSC2)=20 ;LITER(U2RS)=20 LITER(V2RS)=20 ;LITER(W2RS)=20 LITER(VWRS)=20 ;LITER(SC2)=20 LITER(SC1)=20 ;LITER(TEM1)=20 ENDIT(W1)=1.0E-03 ;ENDIT(EP)=1.0E-03 ENDIT(VTRS)=1.0E-03 ;ENDIT(VSC1)=1.0E-03 ENDIT(VSC2)=1.0E-03 ;ENDIT(U2RS)=1.0E-03 ENDIT(V2RS)=1.0E-03 ;ENDIT(W2RS)=1.0E-03 ENDIT(VWRS)=1.0E-03 ;ENDIT(SC2)=1.0E-03 ENDIT(SC1)=1.0E-03 ;ENDIT(TEM1)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(V1,LINRLX,1.) RELAX(W1,FALSDT,0.1) RELAX(KE,LINRLX,1.) RELAX(EP,FALSDT,0.1) RELAX(PRPS,LINRLX,1.) RELAX(DSDY,LINRLX,1.) RELAX(VTRS,FALSDT,0.05) RELAX(VSC1,FALSDT,0.05) RELAX(VSC2,FALSDT,0.05) RELAX(U2RS,FALSDT,0.05) RELAX(V2RS,FALSDT,0.05) RELAX(W2RS,FALSDT,0.05) RELAX(VWRS,FALSDT,0.05) RELAX(LEN1,LINRLX,1.) RELAX(ENUT,LINRLX,1.) RELAX(W2DK,LINRLX,1.) RELAX(V2DK,LINRLX,1.) RELAX(U2DK,LINRLX,1.) RELAX(VWDK,LINRLX,1.) RELAX(FWAL,LINRLX,1.) RELAX(EPDK,LINRLX,1.) RELAX(PK,LINRLX,1.) RELAX(DVW,LINRLX,1.) RELAX(PU2,LINRLX,1.) RELAX(PV2,LINRLX,1.) RELAX(PW2,LINRLX,1.) RELAX(PVW,LINRLX,1.) RELAX(DWDY,LINRLX,1.) RELAX(SC2,FALSDT,0.5) RELAX(SC1,FALSDT,0.5) RELAX(TEM1,FALSDT,0.5) OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits 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(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(DSDY)=1.0E+10 ;VARMIN(DSDY)=-1.0E+10 VARMAX(VTRS)=1.0E+10 ;VARMIN(VTRS)=-1.0E+10 VARMAX(VSC1)=1.0E+10 ;VARMIN(VSC1)=-1.0E+10 VARMAX(VSC2)=1.0E+10 ;VARMIN(VSC2)=-1.0E+10 VARMAX(U2RS)=1.0E+10 ;VARMIN(U2RS)=-1.0E+10 VARMAX(V2RS)=1.0E+10 ;VARMIN(V2RS)=-1.0E+10 VARMAX(W2RS)=1.0E+10 ;VARMIN(W2RS)=-1.0E+10 VARMAX(VWRS)=1.0E+10 ;VARMIN(VWRS)=-1.0E+10 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10 VARMAX(W2DK)=1.0E+10 ;VARMIN(W2DK)=-1.0E+10 VARMAX(V2DK)=1.0E+10 ;VARMIN(V2DK)=-1.0E+10 VARMAX(U2DK)=1.0E+10 ;VARMIN(U2DK)=-1.0E+10 VARMAX(VWDK)=1.0E+10 ;VARMIN(VWDK)=-1.0E+10 VARMAX(FWAL)=1.0E+10 ;VARMIN(FWAL)=-1.0E+10 VARMAX(EPDK)=1.0E+10 ;VARMIN(EPDK)=-1.0E+10 VARMAX(PK)=1.0E+10 ;VARMIN(PK)=-1.0E+10 VARMAX(DVW)=1.0E+10 ;VARMIN(DVW)=-1.0E+10 VARMAX(PU2)=1.0E+10 ;VARMIN(PU2)=-1.0E+10 VARMAX(PV2)=1.0E+10 ;VARMIN(PV2)=-1.0E+10 VARMAX(PW2)=1.0E+10 ;VARMIN(PW2)=-1.0E+10 VARMAX(PVW)=1.0E+10 ;VARMIN(PVW)=-1.0E+10 VARMAX(DWDY)=1.0E+10 ;VARMIN(DWDY)=-1.0E+10 VARMAX(SC2)=1.0E+10 ;VARMIN(SC2)=-1.0E+10 VARMAX(SC1)=1.0E+10 ;VARMIN(SC1)=-1.0E+10 VARMAX(TEM1)=1.0E+10 ;VARMIN(TEM1)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND GENK = T RSTM = T YPLS = T FDFSOL = T PARSOL = F IRSMSM = 2 ISG62 = 1 CSG10 ='q1' SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,MATERIAL,36,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(V1,Y,N,Y,N,N,N) OUTPUT(W1,Y,N,Y,Y,Y,Y) OUTPUT(KE,Y,N,Y,N,N,N) OUTPUT(EP,Y,N,Y,Y,Y,Y) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(DSDY,Y,N,Y,N,N,N) OUTPUT(VTRS,Y,N,Y,Y,Y,Y) OUTPUT(VSC1,Y,N,Y,Y,Y,Y) OUTPUT(VSC2,Y,N,Y,Y,Y,Y) OUTPUT(U2RS,Y,N,Y,Y,Y,Y) OUTPUT(V2RS,Y,N,Y,Y,Y,Y) OUTPUT(W2RS,Y,N,Y,Y,Y,Y) OUTPUT(VWRS,Y,N,Y,Y,Y,Y) OUTPUT(LEN1,Y,N,Y,N,N,N) OUTPUT(ENUT,Y,N,Y,N,N,N) OUTPUT(W2DK,Y,N,Y,N,N,N) OUTPUT(V2DK,Y,N,Y,N,N,N) OUTPUT(U2DK,Y,N,Y,N,N,N) OUTPUT(VWDK,Y,N,Y,N,N,N) OUTPUT(FWAL,Y,N,Y,N,N,N) OUTPUT(EPDK,Y,N,Y,N,N,N) OUTPUT(PK,Y,N,Y,N,N,N) OUTPUT(DVW,Y,N,Y,N,N,N) OUTPUT(PU2,Y,N,Y,N,N,N) OUTPUT(PV2,Y,N,Y,N,N,N) OUTPUT(PW2,Y,N,Y,N,N,N) OUTPUT(PVW,Y,N,Y,N,N,N) OUTPUT(DWDY,Y,N,Y,N,N,N) OUTPUT(SC2,Y,N,Y,Y,Y,Y) OUTPUT(SC1,Y,N,Y,Y,Y,Y) OUTPUT(TEM1,Y,N,Y,Y,Y,Y) WALPRN = T ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 29 ;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 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 10000 IPLTF = 1 ;IPLTL = -1 ;NPLT = 20 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