L220

TALK=T;RUN(1,1) L220

  
  Display
  A turbulent combustible mixture flows steadily from left to right
  towards and into a hot porous plug, which is represented by 
  fixing the mass fraction of the fully-reacted gas to unity there.
  
  A distorted distance scale is used, so that the whole flame lies
  between x=0 and x=1.
  
  The RATE value is fixed arbitrarily, as is also the diffusivity 
  viscosity.
  
  The case can be used for testing the influence of the fineness
  of the population grid on the reaction-rate distribution.
  
  ENDDIS
  Declarations: ***************************************************
#mfm
  Data: ***********************************************************
   the coupling and chemical rate constants: conmix, conrea, ratexp
CONMIX=1; CONREA=0.025; RATEXP=5
NFLR=40; NFLF=1; 
NFLUIDS=NFLR*NFLF
TEXT(1D Flame Propgn; Conmix=:Conmix: Conrea=:Conrea:
TITLE

GRDPWR(X,40,1,1)

MFMMOD=XFL1D
STORE(REAC)
#mfm1
  ! 
#mfm2

FIINIT(RATE)=1.0; ENUL=1 ! set micro-mixing rate and diffusivity to
                         ! unity
PATCH(COLD,CELL,1,1,1,1,1,1,1,LSTEP)  ! upstream boundary
PATCH(HOT,CELL,NX,NX,1,1,1,1,1,LSTEP) ! downstream boundary
DO II = 1,NFLUIDS                     ! settings for all fluids
 RELAX(F:II:,LINRLX,0.5)                ! under-relaxation
 FIINIT(F:II:)=0.0                      ! initial values
 COVAL(COLD,F:II:,FIXVAL,0.0)           ! upstream boundary values
 COVAL(HOT,F:II:,FIXVAL,0.0)            ! downstream boundary values
ENDDO
COVAL(COLD,F1,FIXVAL,1.0)             ! fluid 1 upstream
COVAL(HOT,F:NFLUIDS:,1.E3,1.0)        ! last fluid downstream
 
LSWEEP=5000
TSTSWP=-1
IXMON=NX/2
distil=t

   line-printer profile plots
IPROF=1
PATCH(AVES,PROFIL,1,NX,1,1,1,1,1,1)
COVAL(AVES,REAC,0,0)

PATCH(EXT,PROFIL,1,NX,1,1,1,1,1,1)    ! profiles of extreme fluids
COVAL(EXT,F1,0,1); COVAL(EXT,F:NFLUIDS:,0,0)

PATCH(INT,PROFIL,1,NX,1,1,1,1,1,1)    ! profiles of some 
DO II=2, :NFLUIDS:-1, 1+:NFLUIDS:/5  ! intermediate fluids
 COVAL(INT,F:II:,0,.03)
ENDDO
 EX(F1  )=   4.197E-01; EX(F2  )=   1.603E-02; EX(F3  )=   7.670E-03
 EX(F4  )=   6.074E-03; EX(F5  )=   5.199E-03; EX(F6  )=   4.627E-03
 EX(F7  )=   4.217E-03; EX(F8  )=   3.905E-03; EX(F9  )=   3.656E-03
 EX(F10 )=   3.452E-03; EX(F11 )=   3.281E-03; EX(F12 )=   3.133E-03
 EX(F13 )=   3.003E-03; EX(F14 )=   2.888E-03; EX(F15 )=   2.785E-03
 EX(F16 )=   2.692E-03; EX(F17 )=   2.606E-03; EX(F18 )=   2.529E-03
 EX(F19 )=   2.457E-03; EX(F20 )=   2.392E-03; EX(F21 )=   2.333E-03
 EX(F22 )=   2.279E-03; EX(F23 )=   2.231E-03; EX(F24 )=   2.188E-03
 EX(F25 )=   2.152E-03; EX(F26 )=   2.122E-03; EX(F27 )=   2.098E-03
 EX(F28 )=   2.082E-03; EX(F29 )=   2.075E-03; EX(F30 )=   2.077E-03
 EX(F31 )=   2.090E-03; EX(F32 )=   2.117E-03; EX(F33 )=   2.162E-03
 EX(F34 )=   2.229E-03; EX(F35 )=   2.328E-03; EX(F36 )=   2.475E-03
 EX(F37 )=   2.701E-03; EX(F38 )=   3.078E-03; EX(F39 )=   4.043E-03
 EX(F40 )=   4.528E-01; EX(REAC)=   2.263E-03; EX(RATE)=   1.000E+00