300 likes | 512 Views
Development of FDO Patterns in the BZ Reaction. Steve Scott University of Leeds. Acknowledgements. Jonnie Bamforth (Leeds) Rita T ó th (Debrecen) Vilmos G áspár (Debrecen) British Council/Hungarian Academy of Science ESF REACTOR programme.
E N D
Development of FDO Patterns in the BZ Reaction Steve Scott University of Leeds
Acknowledgements • Jonnie Bamforth (Leeds) • Rita Tóth (Debrecen) • Vilmos Gáspár (Debrecen) • British Council/Hungarian Academy of Science • ESF REACTOR programme
Kuznetsov, Andresen, Mosekilde, Dewel, Borckmans Flow Distributed Oscillations • patterns without differential diffusion or flow • Very simple reactor configuration: plug-flow tubular reactor fed from CSTR • reaction run under conditions so it is oscillatory in batch, but steady-state in CSTR
Simple explanation • CSTR ensures each “droplet” leaves with same “phase” • Oscillations occur in each droplet at same time after leaving CSTR and, hence, at same place in PFR
Explains: existence of stationary patterns need for “oscillatory batch” reaction BZ system with f = 0.17 cm s-1 [BrO3-] = 0.24 M, H+ = 0.15M[MA] = 0.4 M, [Ferroin] = 7 ´ 10-4 M Images taken at 2 min intervals
Doesn’t explain critical flow velocity nonlinear dependence of wavelength on flow velocity other responses observed, especially the dynamics of pattern development
Analysis • Oregonator model: Has a uniform steady state uss, vss
Perturbation: u = U + uss, v = V + vss linearised equations Seek solutions of the form
Dispersion relation Tr = j11 + j22 D = j11j22 – j12j21
Absolute to Convective Instability Look for zero group velocity, i.e. find k = k0such that gives so Setting Im(w(k0)) = 0 gives fAC
Bifurcation to Stationary Patterns Required condition is w = 0 with Im(k) = 0 Setting w = 0 yields So Im(k) = 0 gives critical flow velocity
Initial Development of Stationary Pattern • Oregonator modele = 0.25f = 1.0q = 8 10-4f = 2 0.4 time units per frame
Experimental verification BZ system with f = 0.17 cm s-1[BrO3-] = 0.2 M, H+ = 0.15M[MA] = 0.4 M, [Ferroin] = 7 ´ 10-4 M
Oregonator model as before, Pattern already established now change f from 2.0 to 4.0 Adjustment of wavelength to change in flow velocity
Nonlinear l-f response e = 0.8 e = 0.5 e = 0.25
Complex Pattern Development e = 0.25f = 1.0q = 8 10-4f = 1.5 0.4 time units per frame
space-time plot f = 1.5
more complexity f = 1.4
CDIMA reaction Patterns but unsteady
Lengyel-Epstein model • a= 0.5f = 5 0.12 time units per frame