Convective chemical-wave propagation in the Belousov-Zhabotinsky reaction

Buoyancy-Driven Chemohydrodynamic Patterns in A + B → Oscillator Two-Layer Stratifications

Chemohydrodynamic patterns due to the interplay of buoyancy-driven instabilities and reaction-diffusion patterns are studied experimentally in a vertical quasi-two-dimensional reactor in which two solutions A and B containing separate reactants of the oscillating Belousov-Zhabotinsky system are placed in contact along a horizontal contact line where excitable or oscillating dynamics can develop. Different types of buoyancy-driven instabilities are selectively induced in the reactive zone depending on the initial density jump between the two layers, controlled here by the bromate salt concentration. Starting from a less dense solution above a denser one, two possible differential diffusion instabilities are triggered depending on whether the fast diffusing sulfuric acid is in the upper or lower solution. Specifically, when the solution containing malonic acid and sulfuric acid is stratified above the one containing the slow-diffusing bromate salt, a diffusive layer convection (DLC) instability is observed with localized convective rolls around the interface. In that case, the reaction-diffusion wave patterns remain localized above the initial contact line, scarcely affected by the flow. If, on the contrary, sulfuric acid diffuses upward because it is initially dissolved in the lower layer, then a double-diffusion (DD) convective mode develops. This triggers fingers across the interface that mix the reactants such that oscillatory dynamics and rippled waves develop throughout the whole reactor. If the denser solution is put on top of the other one, then a fast developing Rayleigh-Taylor (RT) instability induces fast mixing of all reactants such that classical reaction-diffusion waves develop later on in the convectively mixed solutions.

Mechanism and Phenomenology of an Oscillating Chemical Reaction

Chemical reactions, which are far from equilibrium, are capable of displaying oscillations in species concentrations and hence in colour, electrode potential, pH and/or temperature. The oscillations arise from the interplay between positive and negative kinetic feedback. Mechanisms for such reactions are presented, along with the rich phenomenology that these systems exhibit, from complex oscillations and chemical waves, to stationary concentration patterns. This review will focus on the Belousov-Zhabotinksy reaction but reference to other reactions will be made where appropriate.

Surface tension- and buoyancy-driven flows across horizontally propagating chemical fronts

Chemical reactions can interplay with hydrodynamic flows to generate various complex phenomena. Because of their relevance in many research areas, chemically-induced hydrodynamic flows have attracted increasing attention in the last decades. In this context, we propose to give a review of the past and recent theoretical and experimental works which have considered the interaction of such flows with chemical fronts, i.e. reactive interfaces, formed between miscible solutions. We focus in particular on the influence of surface tension- (Marangoni) and buoyancy-driven flows on the dynamics of chemical fronts propagating horizontally in the gravity field.

Control of chemical chaos through medium viscosity in a batch ferroin-catalysed Belousov–Zhabotinsky reaction

A macroscopic parameter, such as medium viscosity, can be used to fine tune chemical chaos in a reaction–diffusion–convection system.

One-Directional Fluidic Flow Induced by Chemical Wave Propagation in a Microchannel

A one-directional flow induced by chemical wave propagation was investigated to understand the origin of its dynamic flow. A cylindrical injection port was connected with a straight propagation channel; the chemical wave was initiated at the injection port. Chemical waves propagated with a constant velocity irrespective of the channel width, indicating that the dynamics of the chemical waves were governed by a geometry-independent interplay between the chemical reaction and diffusion. In contrast, the velocity of the one-directional flow was dependent on the channel width. Furthermore, enlargement of the injection port volume increased the flow velocity and volume flux. These results imply that the one-directional flow in the microchannel is due to a hydrodynamic effect induced in the injection port. Spectroscopic analysis of a pH indicator revealed the simultaneous behavior between the pH increase near the injection port and the one-directional flow. Hence, we can conclude that the one-directional flow in the microchannel with chemical wave propagation was caused by a proton consumption reaction in the injection port, probably through liquid volume expansion by the reaction products and the reaction heat. It is a characteristic feature of the present system that the hydrodynamic flow started from the chemical wave initiation point and not the propagation wavefront, as observed for previous systems.

Propagation Behaviors of an Acid Wavefront Through a Microchannel Junction

Waves in reaction-diffusion systems yield a wealth of dynamic self-assembling phenomena in nature. Recent studies have been devoted to utilizing these active waves in conjunction with microscale technology. To provide a compass for controlling reaction-diffusion waves in microspaces, we have investigated the propagation behavior of one specific variety of the reaction-diffusion wave: an acid wave that utilizes an autocatalytic proton-production reaction. Furthermore, the acid wave that we have investigated occurs in a microchannel with a junction connecting circular and straight regions. The obtained results were compared with a neutralization wave that involves only a neutralization reaction. The acid wave was ignited by the addition of the appropriate amount of H2SO4 into the circular region that was filled with a substrate solution, where proton-consuming and proton-producing reactions followed a rapid neutralization reaction. At this stage, the wave penetrated and propagated into the channel region. Comparison between the acid and the neutralization waves clarified that the acid wave required a minimum threshold of H2SO4 concentration in order to be ignited and that the propagation of the acid wave was temporarily delayed because of the presence of intermediate chemical reaction steps. Furthermore, the propagation dynamics was found to be tuned through the configuration of the microchannel. The importance of microchannel configuration, especially for systems with a junction connecting different shapes, is discussed in terms of Fick's law and in terms of the proton flux from the circular to the straight regions.

Convective dynamics of traveling autocatalytic fronts in a modulated gravity field

Modulation of the gravity field, spanning from the hyper-gravity to micro-gravity of a parabolic flight, reveals the contribution of Marangoni flow in a propagating reaction front with an open air–liquid interface.

Marangoni-driven convection around exothermic autocatalytic chemical fronts in free-surface solution layers

Gradients of concentration and temperature across exothermic chemical fronts propagating in free-surface solution layers can initiate Marangoni-driven convection. We investigate here the dynamics arising from such a coupling between exothermic autocatalytic reactions, diffusion, and Marangoni-driven flows. To this end, we numerically integrate the incompressible Navier-Stokes equations coupled through the tangential stress balance to evolution equations for the concentration of the autocatalytic product and for the temperature. A solutal and a thermal Marangoni numbers measure the coupling between reaction-diffusion processes and surface-driven convection. In the case of an isothermal system, the asymptotic dynamics is characterized by a steady fluid vortex traveling at a constant speed with the front, deforming and accelerating it [L. Rongy and A. De Wit, J. Chem. Phys. 124, 164705 (2006)]. We analyze here the influence of the reaction exothermicity on the dynamics of the system in both cases of cooperative and competitive solutal and thermal effects. We show that exothermic fronts can exhibit new unsteady spatio-temporal dynamics when the solutal and thermal effects are antagonistic. The influence of the solutal and thermal Marangoni numbers, of the Lewis number (ratio of thermal diffusivity over molecular diffusivity), and of the height of the liquid layer on the spatio-temporal front evolution are investigated.

Chemical activity induces dynamical force with global structure in a reaction-diffusion-convection system

We found a rotating global structure induced by the dynamical force of local chemical activity in a thin solution layer of excitable Belousov-Zhabotinsky reaction coupled with diffusion. The surface flow and deformation associated with chemical spiral waves (wavelength about 1 mm) represents a global unidirectional structure and a global tilt in the entire Petri dish (100 mm in diameter), respectively. For these observations, we scanned the condition of hierarchal pattern selection. From this result, the bromomalonic acid has an important role to induce the rotating global structure. An interaction between a reaction-diffusion process and a surface-tension-driven effect leads to such hierarchal pattern with different scales.

Influence of thermal effects on buoyancy-driven convection around autocatalytic chemical fronts propagating horizontally

The spatiotemporal dynamics of vertical autocatalytic fronts traveling horizontally in thin solution layers closed to the air can be influenced by buoyancy-driven convection induced by density gradients across the front. We perform here a combined experimental and theoretical study of the competition between solutal and thermal effects on such convection. Experimentally, we focus on the antagonistic chlorite-tetrathionate reaction for which solutal and thermal contributions to the density jump across the front have opposite signs. We show that in isothermal conditions the heavier products sink below the lighter reactants, providing an asymptotic constant finger shape deformation of the front by convection. When thermal effects are present, the hotter products, on the contrary, climb above the reactants for strongly exothermic conditions. These various observations as well as the influence of the relative weight of the solutal and thermal effects and of the thickness of the solution layer on the dynamics are discussed in terms of a two-dimensional reaction-diffusion-convection model parametrized by a solutal R(C) and a thermal R(T) Rayleigh number.

Interaction between buoyancy and diffusion-driven instabilities of propagating autocatalytic reaction fronts. II. Nonlinear simulations

The nonlinear dynamics resulting from the interplay between diffusive and buoyancy-driven Rayleigh-Taylor (RT) instabilities of autocatalytic traveling fronts are analyzed numerically for fronts ascending or descending in the gravity field and for various values of the relevant parameters, the Rayleigh numbers R(a) and R(b) of the reactant A and autocatalytic product B, respectively, and the ratio D=D(B)/D(A) of the diffusion coefficients of the two key chemical species. The interaction between the coarsening dynamics characteristic of the RT instability and the fixed short wavelength dynamics of the diffusive instability leads in some parameter regimes to complex dynamics dominated by the irregular succession of birth and death of fingers. Large single convective fingers with a tip deformed by the short wavelength diffusive instability are also observed. If D is sufficiently small and the RT instability is active, the concentration of the slower diffusing species B can be convected to values above its fully reacted concentration. Experimental conditions that would allow the observation of the dynamics predicted here are described.

Interaction between buoyancy and diffusion-driven instabilities of propagating autocatalytic reaction fronts. I. Linear stability analysis

The interaction between buoyancy-driven and diffusion-driven instabilities that can develop along a propagating reaction front is discussed for a system based on an autocatalytic reaction. Twelve different cases are possible depending on whether the front is ascending or descending in the gravity field, whether the reactant is heavier or lighter than the products, and whether the reactant diffuses faster, slower, or at the same rate as the product. A linear stability analysis (LSA) is undertaken, in which dispersion curves (plots of the growth rate sigma against wave number k) are derived for representative cases as well as an asymptotic analysis for small wave numbers. The results from the LSA indicate that, when the initial reactant is denser than the reaction products, upward propagating fronts remain unstable with the diffusion-driven instability enhancing this instability. Buoyantly stable downward propagating fronts become unstable when the system is also diffusionally unstable. When the initial reactant is lighter than the reaction products, any diffusionally unstable upward propagating front is stabilized by small buoyancy effects. A diffusional instability enhances the buoyant instability of a downward propagating front with there being a very strong interaction between these effects in this case.

Buoyancy-driven convection around chemical fronts traveling in covered horizontal solution layers

Density differences across an autocatalytic chemical front traveling horizontally in covered thin layers of solution trigger hydrodynamic flows which can alter the concentration profile. We theoretically investigate the spatiotemporal evolution and asymptotic dynamics resulting from such an interplay between isothermal chemical reactions, diffusion, and buoyancy-driven convection. The studied model couples the reaction-diffusion-convection evolution equation for the concentration of an autocatalytic species to the incompressible Stokes equations ruling the evolution of the flow velocity in a two-dimensional geometry. The dimensionless parameter of the problem is a solutal Rayleigh number constructed upon the characteristic reaction-diffusion length scale. We show numerically that the asymptotic dynamics is one steady vortex surrounding, deforming, and accelerating the chemical front. This chemohydrodynamic structure propagating at a constant speed is quite different from the one obtained in the case of a pure hydrodynamic flow resulting from the contact between two solutions of different density or from the pure reaction-diffusion planar traveling front. The dynamics is symmetric with regard to the middle of the layer thickness for positive and negative Rayleigh numbers corresponding to products, respectively, lighter or heavier than the reactants. A parametric study shows that the intensity of the flow, the propagation speed, and the deformation of the front are increasing functions of the Rayleigh number and of the layer thickness. In particular, the asymptotic mixing length and reaction-diffusion-convection speed both scale as square root Ra for Ra>5. The velocity and concentration fields in the asymptotic dynamics are also found to exhibit self-similar properties with Ra. A comparison of the dynamics in the case of a monostable versus bistable kinetics is provided. Good agreement is obtained with experimental data on the speed of iodate-arsenous acid fronts propagating in horizontal capillaries. We furthermore compare the buoyancy-driven dynamics studied here to Marangoni-driven deformation of traveling chemical fronts in solution open to the air in the absence of gravity previously studied in the same geometry [L. Rongy and A. De Wit, J. Chem. Phys. 124, 164705 (2006)].

Active media under rotational forcing

The bubble-free Belousov-Zhabotinsky reaction has been used to study the effects of centrifugal forces on autowave propagation. The reaction parameters were chosen such that the system oscillates naturally creating target waves. In the present study, the system was forced to rotate with a constant velocity around a central axis. In studying the effects of such a forcing on the system, we focused on target dynamics. The system reacts to this forcing in different ways, the most spectacular being a dramatic increase in the period of the target, the effect growing stronger as we move away from the center of rotation. A numerical study was carried out using the two-variable Oregonator model, modified to include convective effects through the diffusion coefficient. The numerical results showed a good qualitative agreement with those of the experiments.

Steady Marangoni flow traveling with chemical fronts

When autocatalytic chemical fronts propagate in thin layers of solution in contact with air, they can induce capillary flows due to surface tension gradients across the front (Marangoni flows). We investigate here such an interplay between autocatalytic reactions, diffusion, and Marangoni effects with a theoretical model coupling the incompressible Navier-Stokes equations to a conservation equation for the autocatalytic product concentration in the absence of gravity and for isothermal conditions. The boundary condition at the open liquid/air interface takes the surface activity of this product into account and introduces the solutal Marangoni number M representing the intensity of the coupling between hydrodynamics and reaction-diffusion processes. Positive and negative Marangoni numbers correspond, respectively, to the cases where the product decreases or increases surface tension behind the front. We show that, in both cases, such coupled systems reach an asymptotic dynamics characterized by a steady fluid vortex traveling at a constant speed with the front and deforming it, with, however, an asymmetry between the results for positive and negative M. A parametric study shows that increased propagation speed, front deformation, and possible transient oscillating dynamics occur when the absolute value of M is increased.

Nuclear Magnetic Resonance Studies of Convection in the 1,4-Cyclohexanedione−Bromate−Acid Reaction

The manifestation and development of convection during pattern formation in the 1,4-cyclohexanedione-acid-bromate reaction was investigated using pulsed gradient spin-echo nuclear magnetic resonance (PGSE NMR) experiments. An apparatus was devised that enabled convection to be probed inside an NMR spectrometer and prevented hydrodynamic motion arising from extraneous sources, such as poor mixing or temperature gradients imposed by the experimental setup. PGSE experiments were performed concurrently with magnetic resonance imaging (MRI) experiments to show that convection arose spontaneously from inhomogeneities associated with the chemical patterns. Quantitative data on diffusion coefficients and hydrodynamic velocities are reported.

Density Changes Accompanying Wave Propagation in the Cerium-Catalyzed Belousov−Zhabotinsky Reaction

Refractive index measurement using an interferometric imaging system and observation of chemical wave shapes were carried out during chemical wave propagation of a cerium-catalyzed Belousov-Zhabotinsky (BZ) reaction. Densities increased as chemical waves propagated in samples without NaBr, and decreased in samples with NaBr. Concentration changes of malonic acid, bromomalonic acid, and BrO3- were estimated from Raman spectral measurements in a stirred batch BZ reaction, and these also exhibited differences between samples with and without NaBr. It is proposed that a reaction subset yielding low molecular weight carboxylic acids is predominant in samples with NaBr, whereas a pathway leading to dibromoacetic acid or tribromoacetic acid production is the major process in samples without NaBr.

Spiral wave meandering induced by fluid convection in an excitable medium

An isothermal reaction-diffusion system is considered in a two-dimensional fluid medium within a gravitational field. Inhomogeneities in the concentration field of the species give rise to a fluid flow due to buoyancy forces. A two-dimensional reaction-diffusion-convection model of an excitable medium is presented. The influence of hydrodynamics on spiral wave dynamics is systematically studied. A kinematic model is also introduced to better understand the mechanisms involved here.

Convection in chemical fronts with quadratic and cubic autocatalysis

Convection in chemical fronts enhances the speed and determines the curvature of the front. Convection is due to density gradients across the front. Fronts propagating in narrow vertical tubes do not exhibit convection, while convection develops in tubes of larger diameter. The transition to convection is determined not only by the tube diameter, but also by the type of chemical reaction. We determine the transition to convection for chemical fronts with quadratic and cubic autocatalysis. We show that quadratic fronts are more stable to convection than cubic fronts. We compare these results to a thin front approximation based on an eikonal relation. In contrast to the thin front approximation, reaction-diffusion models show a transition to convection that depends on the ratio between the kinematic viscosity and the molecular diffusivity. (c) 2002 American Institute of Physics.

Convective structures in a two-layer gel-liquid excitable medium

The onset of convection due to wave propagation is investigated in the framework of the Belousov-Zhabotinsky reaction. Numerical calculations are based on a three variable Oregonator model coupled with the Navier-Stokes hydrodynamic equations under the Boussineq approximation for a system consisting of two layers, a liquid and a gel, both in close contact through an interface where chemical concentration exchange is allowed. The influence on the formation of convective rolls associated to wave front propagation is studied in terms of the exchange rate through the interface, the liquid layer width, and the coupling strength between the fluid flow and chemical dynamics. Waves are initiated on the surface of the gel and this perturbation is allowed to propagate into the liquid initiating either two counterrotating convective cells (at both sides of the front) or a disordered pattern.

Electric-field-induced front deformation of Belousov-Zhabotinsky waves

Experimental evidence is given of deformations in the vertical profiles of ferroin-catalyzed Belousov-Zhabotinsky waves propagating in thin horizontal cuvettes under an imposed dc electric field. While no deformations are seen in a zero field, for low negative field, or for any positive field, a pronounced S-shaped deformation does occur when a wave is exposed to a negative field above some critical magnitude. The observed phenomena are discussed on the basis of the convective flows that are assumed to increase in the negative field resulting from changes in the longitudinal profile of the wave.