1 / 43

Chaos and Control in Combustion

Chaos and Control in Combustion. Steve Scott School of Chemistry University of Leeds. Outline. Review of H 2 and CO combustion Use of flow reactors Oscillatory ignition Mechanistic comments Complex oscillations Chaos Control of Chaos. The H 2 + O 2 reaction.

Audrey
Download Presentation

Chaos and Control in Combustion

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Chaos and Controlin Combustion Steve Scott School of Chemistry University of Leeds

  2. Outline • Review of H2 and CO combustion • Use of flow reactors • Oscillatory ignition • Mechanistic comments • Complex oscillations • Chaos • Control of Chaos

  3. The H2 + O2 reaction The classic example of a branched chain reaction simplest combustion reaction etc.

  4. H2 + O2 branching cycle H2 OH + O H2O + H rds H + O2 H + OH H2 H2 H2O + H Overall: H + 3H2 + O2 3H + 2H2O rb = 2 kb [H][O2]

  5. Mechanism at 2nd limit balance between chain branching and gas-phase (termolecular) termination H + O2 3 H rb = 2kb[H][O2] H + O2 + M ® HO2 + M rt = kt[H][O2][M]

  6. Then: where is the net branching factor. f < 0: evolve to low steady state f> 0: exponential growth

  7. Condition for limit Critical condition is  = 0 2 kb = kt [M]

  8. Studies in flow reactors • Continuous-flow, well-stirred reactor (CSTR) • Also shows p-Ta ignition limits • Study in vicinity of 2nd limit

  9. p-Ta diagram for H2 + O2 in CSTR tres = 8 s

  10. Oscillatory ignition

  11. How does oscillation vary with experimental operating conditions?

  12. “Limit cycles” Oscillation in time corresponds to “lapping” on limit cycle

  13. Extinction at low Ta tres = 2 s

  14. “SNIPER” bifurcation

  15. More complex behaviour different oscillations at same operating conditions: birhythmicity

  16. Mixed-mode oscillations H2-rich systems

  17. Why do oscillations occur? • Need to consider “third body efficiencies” remember ignition limit condition 2 kb = kt [M] this assumes all species have same ability to stabilise HO2-species in fact, different species have different efficiencies: aO2 ~ 0.3, aH2O ~ 6 so: overall efficiency of reacting mixture changes with composition

  18. Allow for this in following way: In ignition region: f> 0, based on reactant composition. After “ignition”, composition now has H2 and O2 replaced by H2O, so overall efficiency is increased, such that for this compositionf < 0. H2O outflow and H2+O2 inflow causes f to increase again – next ignition can develop.

  19. Explains: oscillatory nature and importance of flow; period varies with Ta – through kb; upper Ta limit to oscillatory region (f > 0 even for “ignited composition”; extinction of oscillations at ignition limit. Doesn’t explain: complex oscillations. Need to include: a few more reactions + temperature effects

  20. CO + O2 in closed vessels • shows p-Ta ignition limit • chemiluminescent reaction (CO2*) “glow” • can get “steady glow” and “oscillatory glow” – the lighthouse effect (Ashmore & Norrish, Linnett) • very sensitive to trace quantities of H-containing species

  21. CO + O2 in a CSTR • p-T ignition limit diagram shows region of “oscillatory ignition”

  22. Complex oscillations

  23. Record data under steady operation

  24. Next-maximum map example chaotic trace next- maximum Map

  25. parameter lower boundary upper boundary value used Temperaturea (K) 786 ( 2) 791 ( 2) 789 O2flowb (sccm) 4.0 ( 0.1) 9.0 ( 0.15) 5.6 CO flowc (sccm) 6.9 ( 0.5) 7.4 ( 0.2) 7.14 Extent of chaotic region for system with p = 19 mmHg. sccm = standard cubic centimetre per minute; awith = 5.6 sccm and fCO = 7.14 sccm; bwith T = 789 K and fCO = 7.14 sccm; cwith T = 789 K and = 5.6 sccm.

  26. A quick guide to maps xn+1 = Axn (1 – xn) 1 < A < 4

  27. iteration of the map

  28. Perturbing the map fixed point shifts

  29. targeting the fixed point Ott, Grebogi, Yorke 1990; Petrov, Peng, Showalter 1991 need to determine : location of fixed point of unperturbed system slope of map in vicinity of fixed point shift in fixed point as system is perturbed

  30. experimental strategy From the experimental time series: • collect enough data to plot the map • fit the data to get the fixed point and the slope in its region • perturb one of the experimental parameters • determine the new map – fit to find shift in fixed point

  31. control constant Can calculate a “control constant” g where m is the slope of the map and dxF/df is the rate of change of the fixed point with some experimental parameter Note: m and dxF/dfcan be measured experimentally

  32. Calculate appropriate perturbation If we observe system and it comes “near to” the fixed point of the map : Dx = x-xF Can calculate the appropriate perturbation to the operating conditions

  33. Exploiting the map Chaos control Map varies with the exptl conditions

  34. Control of Chaos by suitable, very small amplitude dynamic perturbations can control chaos perturbations determined from Experiment Davies et al., J. Phys. Chem. A: 16/11/00

  35. some unexpected features control transient time depends on how long perturbation is applied for

  36. optimal control occursfor perturbation applied foronly 25% of oscillatory period

  37. Conclusions • Oscillations, including complex oscillations and even chaotic evolution, arise naturally in chemical reactions as a consequence of “normal” mechanisms with “feedback” • Chaos occurs for a range of experimental conditions. • Chaotic systems can be “controlled” using simple experimental strategies • These need no information regarding the chemical mechanisms and we can determine all the parameters necessary from experiments even if only one signal can be measured

  38. Acknowledgements Barry Johnson Matt Davies, Mark Tinsley, Peter Halford-Maw Istvan Kiss, Vilmos Gaspar (Debrecen) British Council – Hungarian Academy ESF Scientific Programme REACTOR

More Related