Bursting phenomena in a simplified Oregonator flow system model

John Rinzel, William C. Troy

Research output: Contribution to journalArticle

Abstract

We investigate a model for the Belousov-Zhabotinskii reaction in a continuous flow, stirred tank reactor. The model consists of a system of three ordinary differential equations derived from a more complicated five-dimensional model introduced by Janz, Vanecek, and Field [J. Chem. Phys. 73, 3132 (1980)]. It is basically a simplified, irreversible, batch reactor Oregonator with flow terms and with instantaneous dependence of the stoichiometric parameter f on the amount of brominated organic substance. Over a range of physical parameters (e.g., lower flow rates) the system exhibits "bursts" of oscillations. As in some experiments [e.g.,. Marek and Svobodova, Biophys. Chem. 3 263, (1975)], the observer sees several spikes followed by an interval of quiescence (IQ) which is subsequently followed by a resumption of the spikes, etc. As Janz et al. found, bursting occurs as f sweeps slowly back and forth through values for which the batch reactor subsystem exhibits bistability and hysteresis. The CSTR solution alternates between the batch model's stable periodic solution, during the oscillatory phase, and the batch model's stable steady state of low oxidation, during the IQ. Analytically and numerically we describe these subsystem solutions as functions of f; bistability arises as a subcritical Hopf bifurcation. Also, we estimate quantitatively features of the CSTR bursts, such as IQ duration, and their dependence on parameter values. For other parameter ranges we find qualitatively different bursting phenomena. In one example, the CSTR system is excitable; there is a stable steady state and the response to an adequate perturbation is an excitation burst of several pulses and then a return to the steady state. At considerably higher flow rates, there are repetitive, single-spike bursts with IQ's of high oxidation; such patterns resemble those calculated by Showalter, Noyes, and Bar-Eli [J. Chem. Phys. 69, 2514 (1978)] and some responses observed experimentally by Schmitz, Graziani, and Hudson [J. Chem. Phys. 67, 3040 (1977)]

Original languageEnglish (US)
Pages (from-to)1775-1789
Number of pages15
JournalThe Journal of chemical physics
Volume76
Issue number4
StatePublished - 1982

Fingerprint

bursts
spikes
reactors
Batch reactors
intervals
flow velocity
Flow rate
Oxidation
oxidation
Hopf bifurcation
Ordinary differential equations
Hysteresis
differential equations
hysteresis
perturbation
oscillations
estimates
pulses
excitation
Experiments

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Bursting phenomena in a simplified Oregonator flow system model. / Rinzel, John; Troy, William C.

In: The Journal of chemical physics, Vol. 76, No. 4, 1982, p. 1775-1789.

Research output: Contribution to journalArticle

@article{9edf35b820bc47619a17eace55d1ff38,
title = "Bursting phenomena in a simplified Oregonator flow system model",
abstract = "We investigate a model for the Belousov-Zhabotinskii reaction in a continuous flow, stirred tank reactor. The model consists of a system of three ordinary differential equations derived from a more complicated five-dimensional model introduced by Janz, Vanecek, and Field [J. Chem. Phys. 73, 3132 (1980)]. It is basically a simplified, irreversible, batch reactor Oregonator with flow terms and with instantaneous dependence of the stoichiometric parameter f on the amount of brominated organic substance. Over a range of physical parameters (e.g., lower flow rates) the system exhibits {"}bursts{"} of oscillations. As in some experiments [e.g.,. Marek and Svobodova, Biophys. Chem. 3 263, (1975)], the observer sees several spikes followed by an interval of quiescence (IQ) which is subsequently followed by a resumption of the spikes, etc. As Janz et al. found, bursting occurs as f sweeps slowly back and forth through values for which the batch reactor subsystem exhibits bistability and hysteresis. The CSTR solution alternates between the batch model's stable periodic solution, during the oscillatory phase, and the batch model's stable steady state of low oxidation, during the IQ. Analytically and numerically we describe these subsystem solutions as functions of f; bistability arises as a subcritical Hopf bifurcation. Also, we estimate quantitatively features of the CSTR bursts, such as IQ duration, and their dependence on parameter values. For other parameter ranges we find qualitatively different bursting phenomena. In one example, the CSTR system is excitable; there is a stable steady state and the response to an adequate perturbation is an excitation burst of several pulses and then a return to the steady state. At considerably higher flow rates, there are repetitive, single-spike bursts with IQ's of high oxidation; such patterns resemble those calculated by Showalter, Noyes, and Bar-Eli [J. Chem. Phys. 69, 2514 (1978)] and some responses observed experimentally by Schmitz, Graziani, and Hudson [J. Chem. Phys. 67, 3040 (1977)]",
author = "John Rinzel and Troy, {William C.}",
year = "1982",
language = "English (US)",
volume = "76",
pages = "1775--1789",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "4",

}

TY - JOUR

T1 - Bursting phenomena in a simplified Oregonator flow system model

AU - Rinzel, John

AU - Troy, William C.

PY - 1982

Y1 - 1982

N2 - We investigate a model for the Belousov-Zhabotinskii reaction in a continuous flow, stirred tank reactor. The model consists of a system of three ordinary differential equations derived from a more complicated five-dimensional model introduced by Janz, Vanecek, and Field [J. Chem. Phys. 73, 3132 (1980)]. It is basically a simplified, irreversible, batch reactor Oregonator with flow terms and with instantaneous dependence of the stoichiometric parameter f on the amount of brominated organic substance. Over a range of physical parameters (e.g., lower flow rates) the system exhibits "bursts" of oscillations. As in some experiments [e.g.,. Marek and Svobodova, Biophys. Chem. 3 263, (1975)], the observer sees several spikes followed by an interval of quiescence (IQ) which is subsequently followed by a resumption of the spikes, etc. As Janz et al. found, bursting occurs as f sweeps slowly back and forth through values for which the batch reactor subsystem exhibits bistability and hysteresis. The CSTR solution alternates between the batch model's stable periodic solution, during the oscillatory phase, and the batch model's stable steady state of low oxidation, during the IQ. Analytically and numerically we describe these subsystem solutions as functions of f; bistability arises as a subcritical Hopf bifurcation. Also, we estimate quantitatively features of the CSTR bursts, such as IQ duration, and their dependence on parameter values. For other parameter ranges we find qualitatively different bursting phenomena. In one example, the CSTR system is excitable; there is a stable steady state and the response to an adequate perturbation is an excitation burst of several pulses and then a return to the steady state. At considerably higher flow rates, there are repetitive, single-spike bursts with IQ's of high oxidation; such patterns resemble those calculated by Showalter, Noyes, and Bar-Eli [J. Chem. Phys. 69, 2514 (1978)] and some responses observed experimentally by Schmitz, Graziani, and Hudson [J. Chem. Phys. 67, 3040 (1977)]

AB - We investigate a model for the Belousov-Zhabotinskii reaction in a continuous flow, stirred tank reactor. The model consists of a system of three ordinary differential equations derived from a more complicated five-dimensional model introduced by Janz, Vanecek, and Field [J. Chem. Phys. 73, 3132 (1980)]. It is basically a simplified, irreversible, batch reactor Oregonator with flow terms and with instantaneous dependence of the stoichiometric parameter f on the amount of brominated organic substance. Over a range of physical parameters (e.g., lower flow rates) the system exhibits "bursts" of oscillations. As in some experiments [e.g.,. Marek and Svobodova, Biophys. Chem. 3 263, (1975)], the observer sees several spikes followed by an interval of quiescence (IQ) which is subsequently followed by a resumption of the spikes, etc. As Janz et al. found, bursting occurs as f sweeps slowly back and forth through values for which the batch reactor subsystem exhibits bistability and hysteresis. The CSTR solution alternates between the batch model's stable periodic solution, during the oscillatory phase, and the batch model's stable steady state of low oxidation, during the IQ. Analytically and numerically we describe these subsystem solutions as functions of f; bistability arises as a subcritical Hopf bifurcation. Also, we estimate quantitatively features of the CSTR bursts, such as IQ duration, and their dependence on parameter values. For other parameter ranges we find qualitatively different bursting phenomena. In one example, the CSTR system is excitable; there is a stable steady state and the response to an adequate perturbation is an excitation burst of several pulses and then a return to the steady state. At considerably higher flow rates, there are repetitive, single-spike bursts with IQ's of high oxidation; such patterns resemble those calculated by Showalter, Noyes, and Bar-Eli [J. Chem. Phys. 69, 2514 (1978)] and some responses observed experimentally by Schmitz, Graziani, and Hudson [J. Chem. Phys. 67, 3040 (1977)]

UR - http://www.scopus.com/inward/record.url?scp=0000358090&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0000358090&partnerID=8YFLogxK

M3 - Article

VL - 76

SP - 1775

EP - 1789

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 4

ER -