```
DISPLAY
Stirred reactor, with 1D (reactedness population), transient.

The is initially full of reacted fluid, measured by F:NFLUIDS: .

Completely-unreacted fluid, measured by F1, enters in the A
stream; no fluid at all enters in the B stream.

The final steady state is of interest; and, in particular,
whether this consists almost completely of low-numbered fluids,
F1, F2 etc, in which case the reactor is said to have been
'extinguished', or whether significant mass fractions of
higher-numbered fluids are to be found.

The final value of the variable AVER  (i.e. average reactedness)
is the simplest indicator. When extinction occurs, it tends
asymptotically to zero; otherwise it maintains a finite value of
between (say) 35 and 100 %.

Its precise value depends upon the 'physical' constants CONMIX,
CONREA, and RATEXP; and also on the 'numerical-grid-related'
quantities, LSTEP and NFLR.
ENDDIS
#mfm
Set the micro-mixing and reaction-rate constants
CONMIX=100.0; CONREA=1.0; RATEXP=5.0
NFLR=100   ! set the number of reactedness subdivisions
Make consequential settings
#mfm1
FLOWA=1.0  ! inflow rate of A stream
FLOWB=0.0  ! inflow rate of B stream
RA=0.0     ! reactedness of A stream
RB=1.0     ! reactedness of B stream
Make further consequential settings
#mfm2
GRDPWR(T,LSTEP,TLAST,1.0)

FIINIT(F:NFLUIDS:)=1.0
re-set final patch argument to LSTEP
PATCH(INFLOWA,CELL,1,1,1,1,1,1,1,LSTEP)
PATCH(INFLOWB,CELL,1,1,1,1,1,1,1,LSTEP)
PATCH(OUTFLOW,VOLUME,1,1,1,1,1,1,1,LSTEP)

REAL(DUM);DUM = 1.0  +  RB * (NFLUIDS-1)
INTEGER(IB);IB=DUM
FLUIDB=F:IB:
FLUIDB
COVAL(INFLOWB,:FLUIDB:,ONLYMS,1.)
conmix
conrea
LSWEEP=25

TEXT(AWSR transient, conmix=:conmix: conrea=:conrea:

SPEDAT(SET,MFM,CONMIX,R,:CONMIX:)   ! transmit data to gxmfm
SPEDAT(SET,MFM,CONREA,R,:CONREA:)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)

SOLVE(RCTD)                         ! reactedness of ideal WSR
PATCH(WHOLE,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP) ! source via In-Form
REAL(CON)
CON=CONREA*(1. + RATEXP) * (2. + RATEXP)
(source of RCTD at WHOLE is :CON:*(1-RCTD)*RCTD^:RATEXP: with line)
FIINIT(RCTD)=1.0                    ! initialise so as to ignite
VARMAX(RCTD)=1.0 ; VARMIN(RCTD)=0.0 ! set limits
RELAX(RCTD,LINRLX,0.1)              ! relax

PATCH(PROFILES,PROFIL,1,1,1,1,1,1,1,LSTEP) ! use line-printer plot
INTEGER(NFL)                               ! for time-wise profiles
NFL=:NFLUIDS:                              ! for some of the fluid
COVAL(PROFILES,F:NFL:,0.0,1.0)             ! mass fractions and of
COVAL(PROFILES,F:NFL-2:,0.0,0.1)
COVAL(PROFILES,F:NFL-4:,0.0,0.1)
COVAL(PROFILES,F:NFL-6:,0.0,0.1)
COVAL(PROFILES,AVER,0.0,1.0)               ! average reactedness &
COVAL(PROFILES,RCTD,0.0,1.0)               ! IWSR reactedness
DISTIL=T
EX(P1  )=   1.000E+00; EX(RCTD)=   9.735E-01; EX(F1  )=   2.001E-02
EX(F2  )=   5.298E-04; EX(F3  )=   5.134E-04; EX(F4  )=   5.163E-04
EX(F5  )=   5.226E-04; EX(F6  )=   5.305E-04; EX(F7  )=   5.397E-04
EX(F8  )=   5.498E-04; EX(F9  )=   5.607E-04; EX(F10 )=   5.723E-04
EX(F11 )=   5.845E-04; EX(F12 )=   5.973E-04; EX(F13 )=   6.107E-04
EX(F14 )=   6.247E-04; EX(F15 )=   6.392E-04; EX(F16 )=   6.544E-04
EX(F17 )=   6.701E-04; EX(F18 )=   6.865E-04; EX(F19 )=   7.034E-04
EX(F20 )=   7.210E-04; EX(F21 )=   7.391E-04; EX(F22 )=   7.579E-04
EX(F23 )=   7.774E-04; EX(F24 )=   7.975E-04; EX(F25 )=   8.183E-04
EX(F26 )=   8.398E-04; EX(F27 )=   8.620E-04; EX(F28 )=   8.849E-04
EX(F29 )=   9.085E-04; EX(F30 )=   9.329E-04; EX(F31 )=   9.580E-04
EX(F32 )=   9.840E-04; EX(F33 )=   1.011E-03; EX(F34 )=   1.038E-03
EX(F35 )=   1.067E-03; EX(F36 )=   1.096E-03; EX(F37 )=   1.126E-03
EX(F38 )=   1.157E-03; EX(F39 )=   1.189E-03; EX(F40 )=   1.222E-03
EX(F41 )=   1.256E-03; EX(F42 )=   1.292E-03; EX(F43 )=   1.328E-03
EX(F44 )=   1.365E-03; EX(F45 )=   1.404E-03; EX(F46 )=   1.444E-03
EX(F47 )=   1.485E-03; EX(F48 )=   1.528E-03; EX(F49 )=   1.572E-03
EX(F50 )=   1.618E-03; EX(F51 )=   1.666E-03; EX(F52 )=   1.715E-03
EX(F53 )=   1.767E-03; EX(F54 )=   1.820E-03; EX(F55 )=   1.876E-03
EX(F56 )=   1.934E-03; EX(F57 )=   1.995E-03; EX(F58 )=   2.059E-03
EX(F59 )=   2.126E-03; EX(F60 )=   2.196E-03; EX(F61 )=   2.269E-03
EX(F62 )=   2.347E-03; EX(F63 )=   2.429E-03; EX(F64 )=   2.516E-03
EX(F65 )=   2.608E-03; EX(F66 )=   2.706E-03; EX(F67 )=   2.810E-03
EX(F68 )=   2.921E-03; EX(F69 )=   3.040E-03; EX(F70 )=   3.167E-03
EX(F71 )=   3.305E-03; EX(F72 )=   3.453E-03; EX(F73 )=   3.613E-03
EX(F74 )=   3.788E-03; EX(F75 )=   3.978E-03; EX(F76 )=   4.187E-03
EX(F77 )=   4.416E-03; EX(F78 )=   4.669E-03; EX(F79 )=   4.950E-03
EX(F80 )=   5.263E-03; EX(F81 )=   5.614E-03; EX(F82 )=   6.011E-03
EX(F83 )=   6.462E-03; EX(F84 )=   6.978E-03; EX(F85 )=   7.575E-03
EX(F86 )=   8.272E-03; EX(F87 )=   9.093E-03; EX(F88 )=   1.007E-02
EX(F89 )=   1.126E-02; EX(F90 )=   1.273E-02; EX(F91 )=   1.458E-02
EX(F92 )=   1.696E-02; EX(F93 )=   2.012E-02; EX(F94 )=   2.451E-02
EX(F95 )=   3.090E-02; EX(F96 )=   4.093E-02; EX(F97 )=   5.854E-02
EX(F98 )=   9.591E-02; EX(F99 )=   1.887E-01; EX(F100)=   2.724E-01
EX(MNSQ)=   2.086E-01; EX(AVER)=   8.947E-01; EX(RATE)=   1.000E+00
```