TALK=T;RUN(1,1) DISPLAY provided by l001.htm Stirred reactor with a 1D population distribution, and reactedness, ranging from zero to 1, as the population- distinguishing attribute. ___________ It is supposed that two streams of fluid | | enter a reactor which is sufficiently A ====> | well-stirred for space-wise differences | stirred | of conditions to be negligible, but not | \\\|/// C===> sufficiently for micro-mixing to be | reactor | complete. B ====> | |___________| The two streams have the same elemental composition; but one may be more reacted than the other. The flow is steady, and the total mass flow rate per unit volume is 1 kg/s m**3 . The chemical reaction rate is supposed to depend upon reactedness in accordance with the formula: rate (in kg/s m**3) = const * (1 - R) / R ** ratexp where R = reactedness, ratexp = a prescribed constant const = conrea / integral [{(1 - R) / R ** ratexp} dR] conrea = a second prescribed constant. The micro-mixing rate (as defined in MFM theory) is taken as conmix kg/s m**3 times the product of the mass fractions of interacting fluids, in accordance with the "promiscuous Mendelian" hypothesis, where conmix is again a prescribed constant. The task is to compute: (1) the fluid-population distribution in the reactor, and (2) the population-average reactedness (3) the root-mean-square of the reactedness fluctuations These values depend upon 7 dimensionless parameters, namely:- (a) the prescribed constants: RATEXP, CONREA & CONMIX; (b) the relative proportions of the two mass flow rates, namely FLOWA and FLOWB (= 1 - FLOWA); (c) the reactednesses of the two streams, RA and RB; and on the number of fluids in terms of which the population distribution is to be described, NFLR. In the following Q1, the values employed are as follows: RATEXP CONREA CONMIX FLOWA RA RB NFLR 5.0 10.0 10.0 0.5 0.0 1.0 25 The opportunity is provided to change NFLR interactively in the satellite, and CONREA and CONMIX during the EARTH run. Cases 1-1 and 102 allow choice of FLOWA and RATEXP. ENDDIS DISPLAY Notes: (1) four macros are provided in order reduce the repetiveness of MFM q1s. (2) they are named mfm, mfm1, mfm2 and mfm3, and are normally called in that order at the top of the q1 file, which then appears as: #mfm settings of nflr, nflf, etc #mfm1 settings of #mfm2 #mfm3 (3) the character variable mfm is declared and set in the always-loaded core-library macro 014.htm, which, in its turn, declares and sets the character variables mfm1, mfm2, mfm3. (4) the values are: mfm = $L004 mfm1 = $L003 mfm2 = $L002 mfm3 = $L001 ENDDIS ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(WSR, 1D Population; Flowa=0.5 ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 101 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 1 XULAST =1. XFRAC(1)=1. ************************************************************ 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(1)=P1 ;NAME(123)=F1 NAME(124)=F2 ;NAME(125)=F3 NAME(126)=F4 ;NAME(127)=F5 NAME(128)=F6 ;NAME(129)=F7 NAME(130)=F8 ;NAME(131)=F9 NAME(132)=F10 ;NAME(133)=F11 NAME(134)=F12 ;NAME(135)=F13 NAME(136)=F14 ;NAME(137)=F15 NAME(138)=F16 ;NAME(139)=F17 NAME(140)=F18 ;NAME(141)=F19 NAME(142)=F20 ;NAME(143)=F21 NAME(144)=F22 ;NAME(145)=F23 NAME(146)=F24 ;NAME(147)=F25 NAME(148)=MNSQ ;NAME(149)=AVER NAME(150)=RATE * 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(F1,Y,Y,N,N,N,Y) SOLUTN(F2,Y,Y,N,N,N,Y) SOLUTN(F3,Y,Y,N,N,N,Y) SOLUTN(F4,Y,Y,N,N,N,Y) SOLUTN(F5,Y,Y,N,N,N,Y) SOLUTN(F6,Y,Y,N,N,N,Y) SOLUTN(F7,Y,Y,N,N,N,Y) SOLUTN(F8,Y,Y,N,N,N,Y) SOLUTN(F9,Y,Y,N,N,N,Y) SOLUTN(F10,Y,Y,N,N,N,Y) SOLUTN(F11,Y,Y,N,N,N,Y) SOLUTN(F12,Y,Y,N,N,N,Y) SOLUTN(F13,Y,Y,N,N,N,Y) SOLUTN(F14,Y,Y,N,N,N,Y) SOLUTN(F15,Y,Y,N,N,N,Y) SOLUTN(F16,Y,Y,N,N,N,Y) SOLUTN(F17,Y,Y,N,N,N,Y) SOLUTN(F18,Y,Y,N,N,N,Y) SOLUTN(F19,Y,Y,N,N,N,Y) SOLUTN(F20,Y,Y,N,N,N,Y) SOLUTN(F21,Y,Y,N,N,N,Y) SOLUTN(F22,Y,Y,N,N,N,Y) SOLUTN(F23,Y,Y,N,N,N,Y) SOLUTN(F24,Y,Y,N,N,N,Y) SOLUTN(F25,Y,Y,N,N,N,Y) SOLUTN(MNSQ,Y,N,N,N,N,Y) SOLUTN(AVER,Y,N,N,N,N,Y) SOLUTN(RATE,Y,N,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(F1,N,Y,Y,Y,N,Y) TERMS(F2,N,Y,Y,Y,Y,Y) TERMS(F3,N,Y,Y,Y,N,Y) TERMS(F4,N,Y,Y,Y,Y,Y) TERMS(F5,N,Y,Y,Y,N,Y) TERMS(F6,N,Y,Y,Y,Y,Y) TERMS(F7,N,Y,Y,Y,N,Y) TERMS(F8,N,Y,Y,Y,Y,Y) TERMS(F9,N,Y,Y,Y,N,Y) TERMS(F10,N,Y,Y,Y,Y,Y) TERMS(F11,N,Y,Y,Y,N,Y) TERMS(F12,N,Y,Y,Y,Y,Y) TERMS(F13,N,Y,Y,Y,N,Y) TERMS(F14,N,Y,Y,Y,Y,Y) TERMS(F15,N,Y,Y,Y,N,Y) TERMS(F16,N,Y,Y,Y,Y,Y) TERMS(F17,N,Y,Y,Y,N,Y) TERMS(F18,N,Y,Y,Y,Y,Y) TERMS(F19,N,Y,Y,Y,N,Y) TERMS(F20,N,Y,Y,Y,Y,Y) TERMS(F21,N,Y,Y,Y,N,Y) TERMS(F22,N,Y,Y,Y,Y,Y) TERMS(F23,N,Y,Y,Y,N,Y) TERMS(F24,N,Y,Y,Y,Y,Y) TERMS(F25,N,Y,Y,Y,N,Y) 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 =0. EL1 =0. TSURR =0. ;TEMP0 =0. PRESS0 =0. DVO1DT =0. ;DRH1DP =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. ENUL =1.0E-05 ;ENUT =0. PRNDTL(F1)=1. ;PRNDTL(F2)=1. PRNDTL(F3)=1. ;PRNDTL(F4)=1. PRNDTL(F5)=1. ;PRNDTL(F6)=1. PRNDTL(F7)=1. ;PRNDTL(F8)=1. PRNDTL(F9)=1. ;PRNDTL(F10)=1. PRNDTL(F11)=1. ;PRNDTL(F12)=1. PRNDTL(F13)=1. ;PRNDTL(F14)=1. PRNDTL(F15)=1. ;PRNDTL(F16)=1. PRNDTL(F17)=1. ;PRNDTL(F18)=1. PRNDTL(F19)=1. ;PRNDTL(F20)=1. PRNDTL(F21)=1. ;PRNDTL(F22)=1. PRNDTL(F23)=1. ;PRNDTL(F24)=1. PRNDTL(F25)=1. PRT(F1)=1. ;PRT(F2)=1. PRT(F3)=1. ;PRT(F4)=1. PRT(F5)=1. ;PRT(F6)=1. PRT(F7)=1. ;PRT(F8)=1. PRT(F9)=1. ;PRT(F10)=1. PRT(F11)=1. ;PRT(F12)=1. PRT(F13)=1. ;PRT(F14)=1. PRT(F15)=1. ;PRT(F16)=1. PRT(F17)=1. ;PRT(F18)=1. PRT(F19)=1. ;PRT(F20)=1. PRT(F21)=1. ;PRT(F22)=1. PRT(F23)=1. ;PRT(F24)=1. PRT(F25)=1. CP1 =1. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.0E-10 ;FIINIT(F1)=1.0E-10 FIINIT(F2)=1.0E-10 ;FIINIT(F3)=1.0E-10 FIINIT(F4)=1.0E-10 ;FIINIT(F5)=1.0E-10 FIINIT(F6)=1.0E-10 ;FIINIT(F7)=1.0E-10 FIINIT(F8)=1.0E-10 ;FIINIT(F9)=1.0E-10 FIINIT(F10)=1.0E-10 ;FIINIT(F11)=1.0E-10 FIINIT(F12)=1.0E-10 ;FIINIT(F13)=1.0E-10 FIINIT(F14)=1.0E-10 ;FIINIT(F15)=1.0E-10 FIINIT(F16)=1.0E-10 ;FIINIT(F17)=1.0E-10 FIINIT(F18)=1.0E-10 ;FIINIT(F19)=1.0E-10 FIINIT(F20)=1.0E-10 ;FIINIT(F21)=1.0E-10 FIINIT(F22)=1.0E-10 ;FIINIT(F23)=1.0E-10 FIINIT(F24)=1.0E-10 ;FIINIT(F25)=1. FIINIT(MNSQ)=1.0E-10 ;FIINIT(AVER)=1.0E-10 FIINIT(RATE)=1. 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(INFLOWA ,CELL , 1, 1, 1, 1, 1, 1, 1, 1) COVAL(INFLOWA ,P1 , FIXFLU ,0.5 ) COVAL(INFLOWA ,F1 ,0. ,1. ) PATCH(INFLOWB ,CELL , 1, 1, 1, 1, 1, 1, 1, 1) COVAL(INFLOWB ,P1 , FIXFLU ,0.5 ) COVAL(INFLOWB ,F25 ,0. ,1. ) PATCH(OUTFLOW ,VOLUME, 1, 1, 1, 1, 1, 1, 1, 1) COVAL(OUTFLOW ,P1 ,1. ,0. ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 400 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-20 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(F1)=20 LITER(F2)=20 ;LITER(F3)=20 LITER(F4)=20 ;LITER(F5)=20 LITER(F6)=20 ;LITER(F7)=20 LITER(F8)=20 ;LITER(F9)=20 LITER(F10)=20 ;LITER(F11)=20 LITER(F12)=20 ;LITER(F13)=20 LITER(F14)=20 ;LITER(F15)=20 LITER(F16)=20 ;LITER(F17)=20 LITER(F18)=20 ;LITER(F19)=20 LITER(F20)=20 ;LITER(F21)=20 LITER(F22)=20 ;LITER(F23)=20 LITER(F24)=20 ;LITER(F25)=20 ENDIT(P1)=1.0E-03 ;ENDIT(F1)=1.0E-03 ENDIT(F2)=1.0E-03 ;ENDIT(F3)=1.0E-03 ENDIT(F4)=1.0E-03 ;ENDIT(F5)=1.0E-03 ENDIT(F6)=1.0E-03 ;ENDIT(F7)=1.0E-03 ENDIT(F8)=1.0E-03 ;ENDIT(F9)=1.0E-03 ENDIT(F10)=1.0E-03 ;ENDIT(F11)=1.0E-03 ENDIT(F12)=1.0E-03 ;ENDIT(F13)=1.0E-03 ENDIT(F14)=1.0E-03 ;ENDIT(F15)=1.0E-03 ENDIT(F16)=1.0E-03 ;ENDIT(F17)=1.0E-03 ENDIT(F18)=1.0E-03 ;ENDIT(F19)=1.0E-03 ENDIT(F20)=1.0E-03 ;ENDIT(F21)=1.0E-03 ENDIT(F22)=1.0E-03 ;ENDIT(F23)=1.0E-03 ENDIT(F24)=1.0E-03 ;ENDIT(F25)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(F1,LINRLX,0.5) RELAX(F2,LINRLX,0.5) RELAX(F3,LINRLX,0.5) RELAX(F4,LINRLX,0.5) RELAX(F5,LINRLX,0.5) RELAX(F6,LINRLX,0.5) RELAX(F7,LINRLX,0.5) RELAX(F8,LINRLX,0.5) RELAX(F9,LINRLX,0.5) RELAX(F10,LINRLX,0.5) RELAX(F11,LINRLX,0.5) RELAX(F12,LINRLX,0.5) RELAX(F13,LINRLX,0.5) RELAX(F14,LINRLX,0.5) RELAX(F15,LINRLX,0.5) RELAX(F16,LINRLX,0.5) RELAX(F17,LINRLX,0.5) RELAX(F18,LINRLX,0.5) RELAX(F19,LINRLX,0.5) RELAX(F20,LINRLX,0.5) RELAX(F21,LINRLX,0.5) RELAX(F22,LINRLX,0.5) RELAX(F23,LINRLX,0.5) RELAX(F24,LINRLX,0.5) RELAX(F25,LINRLX,0.5) RELAX(MNSQ,LINRLX,1.) RELAX(AVER,LINRLX,1.) RELAX(RATE,LINRLX,1.) OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10 VARMAX(F1)=1. ;VARMIN(F1)=0. VARMAX(F2)=1. ;VARMIN(F2)=0. VARMAX(F3)=1. ;VARMIN(F3)=0. VARMAX(F4)=1. ;VARMIN(F4)=0. VARMAX(F5)=1. ;VARMIN(F5)=0. VARMAX(F6)=1. ;VARMIN(F6)=0. VARMAX(F7)=1. ;VARMIN(F7)=0. VARMAX(F8)=1. ;VARMIN(F8)=0. VARMAX(F9)=1. ;VARMIN(F9)=0. VARMAX(F10)=1. ;VARMIN(F10)=0. VARMAX(F11)=1. ;VARMIN(F11)=0. VARMAX(F12)=1. ;VARMIN(F12)=0. VARMAX(F13)=1. ;VARMIN(F13)=0. VARMAX(F14)=1. ;VARMIN(F14)=0. VARMAX(F15)=1. ;VARMIN(F15)=0. VARMAX(F16)=1. ;VARMIN(F16)=0. VARMAX(F17)=1. ;VARMIN(F17)=0. VARMAX(F18)=1. ;VARMIN(F18)=0. VARMAX(F19)=1. ;VARMIN(F19)=0. VARMAX(F20)=1. ;VARMIN(F20)=0. VARMAX(F21)=1. ;VARMIN(F21)=0. VARMAX(F22)=1. ;VARMIN(F22)=0. VARMAX(F23)=1. ;VARMIN(F23)=0. VARMAX(F24)=1. ;VARMIN(F24)=0. VARMAX(F25)=1. ;VARMIN(F25)=0. VARMAX(MNSQ)=1. ;VARMIN(MNSQ)=0. VARMAX(AVER)=1. ;VARMIN(AVER)=0. VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND PARSOL = F ISG62 = 1 SPEDAT(SET,MFM,MFMMOD,C,MFM) SPEDAT(SET,MFM,NFLUIDS,I,25) SPEDAT(SET,MFM,NFLR,I,25) SPEDAT(SET,MFM,NFLF,I,1) SPEDAT(SET,MFM,CONREA,R,10.) SPEDAT(SET,MFM,CONMIX,R,10.) SPEDAT(SET,MFM,RATEXP,R,5.) SPEDAT(SET,MFM,HIGHF,R,1.) SPEDAT(SET,MFM,LOWF,R,0.) SPEDAT(SET,MFM,VISCON,R,0.) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=0. ;EX(F1)=0. EX(F2)=0. ;EX(F3)=0. EX(F4)=0. ;EX(F5)=0. EX(F6)=0. ;EX(F7)=0. EX(F8)=0. ;EX(F9)=0. EX(F10)=0. ;EX(F11)=0. EX(F12)=0. ;EX(F13)=0. EX(F14)=0. ;EX(F15)=0. EX(F16)=0. ;EX(F17)=0. EX(F18)=0. ;EX(F19)=0. EX(F20)=0. ;EX(F21)=0. EX(F22)=0. ;EX(F23)=0. EX(F24)=0. ;EX(F25)=0. EX(MNSQ)=0. ;EX(AVER)=0. EX(RATE)=0. ************************************************************ 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(F1,Y,N,Y,Y,Y,Y) OUTPUT(F2,Y,N,Y,Y,Y,Y) OUTPUT(F3,Y,N,Y,Y,Y,Y) OUTPUT(F4,Y,N,Y,Y,Y,Y) OUTPUT(F5,Y,N,Y,Y,Y,Y) OUTPUT(F6,Y,N,Y,Y,Y,Y) OUTPUT(F7,Y,N,Y,Y,Y,Y) OUTPUT(F8,Y,N,Y,Y,Y,Y) OUTPUT(F9,Y,N,Y,Y,Y,Y) OUTPUT(F10,Y,N,Y,Y,Y,Y) OUTPUT(F11,Y,N,Y,Y,Y,Y) OUTPUT(F12,Y,N,Y,Y,Y,Y) OUTPUT(F13,Y,N,Y,Y,Y,Y) OUTPUT(F14,Y,N,Y,Y,Y,Y) OUTPUT(F15,Y,N,Y,Y,Y,Y) OUTPUT(F16,Y,N,Y,Y,Y,Y) OUTPUT(F17,Y,N,Y,Y,Y,Y) OUTPUT(F18,Y,N,Y,Y,Y,Y) OUTPUT(F19,Y,N,Y,Y,Y,Y) OUTPUT(F20,Y,N,Y,Y,Y,Y) OUTPUT(F21,Y,N,Y,Y,Y,Y) OUTPUT(F22,Y,N,Y,Y,Y,Y) OUTPUT(F23,Y,N,Y,Y,Y,Y) OUTPUT(F24,Y,N,Y,Y,Y,Y) OUTPUT(F25,Y,N,Y,Y,Y,Y) OUTPUT(MNSQ,Y,N,Y,N,N,N) OUTPUT(AVER,Y,N,Y,N,N,N) OUTPUT(RATE,Y,N,Y,N,N,N) ************************************************************ 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 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 0 ;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