University of Surrey

Test tubes in the lab Research in the ATI Dance Research

A scaling analysis of the effects of natural convection, when Sal'nikov's reaction: P → A → B occurs, together with diffusion and heat transfer in a batch reactor

Campbell, AN, Cardoso, SSS and Hayhurst, AN (2006) A scaling analysis of the effects of natural convection, when Sal'nikov's reaction: P → A → B occurs, together with diffusion and heat transfer in a batch reactor Chemical Engineering Research and Design, 84 (7). pp. 553-561.

very final ms update.pdf - Accepted version Manuscript
Available under License : See the attached licence file.

Download (878kB)
Text (licence)

Download (33kB)


Sal’nikov's chemical reaction in its simplest form consists of two consecutive first-order steps, producing a product B from a precursor P via an active intermediate A, in P→A→B. The first step is assumed here to be thermoneutral, with zero activation energy, whilst the second step is exothermic and has a positive activation energy. These properties make this mechanism one of the simplest to display thermokinetic oscillations, as seen in e.g., cool flames or a batch reactor. We first consider a pure gas, P, undergoing Sal’nikov's reaction in a closed spherical vessel, whose walls are held at a constant temperature. Natural convection becomes significant once the temperature is high enough for the Rayleigh number (Ra) to reach ∼ 103. The subsequent behaviour of the system depends on the interaction between convection, diffusion of heat and mass, and chemical kinetics. By examining the governing equations, we develop and evaluate scales for the characteristic velocity, the concentration of the intermediate A and the temperature rise during the progress of the reaction, for the two extreme cases when transport is dominated, in turn, by diffusion and then by natural convection. These scales depend on the characteristic timescales for the interacting phenomena of chemical reaction, diffusion and natural convection. Typically, the characteristic velocity in a relatively small reactor of radius 0.27 m is as large as 0.3 m s−1, when the temperature rise is ≈ 100 K near the centre of the vessel. Our theoretical predictions are verified by full numerical simulations. Oscillations of both the temperature and the concentration of the intermediate, A, are considered; the accompanying flow field proves to be toroidal, with the fluid ascending close to the reactor's axis, but descending adjacent to its walls. In addition, the effects of such process variables as the initial temperature of the batch reactor and its contents, the pressure and also the size of the reactor are all assessed, together with a consideration of what happens when the reaction proceeds in the liquid phase. In this case, because of the different physical properties of a liquid and a gas, natural convection is more intense than in the gas-phase and is quite likely to lead to turbulence and good mixing.

Item Type: Article
Divisions : Faculty of Engineering and Physical Sciences > Chemical and Process Engineering
Authors :
Campbell, AN
Cardoso, SSS
Hayhurst, AN
Date : 2006
DOI : 10.1205/cherd.05213
Related URLs :
Additional Information : NOTICE: this is the author’s version of a work that was accepted for publication in Chemical Engineering Research and Design. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Chemical Engineering Research and Design, 84(7), July 2006, DOI 10.1205/cherd.05213.
Depositing User : Symplectic Elements
Date Deposited : 16 May 2013 10:52
Last Modified : 31 Oct 2017 14:42

Actions (login required)

View Item View Item


Downloads per month over past year

Information about this web site

© The University of Surrey, Guildford, Surrey, GU2 7XH, United Kingdom.
+44 (0)1483 300800