Bifurcations in spiral tip dynamics induced by natural convection in the Belousov–Zhabotinsky reaction

Marangoni- vs. buoyancy-driven flows: competition for spatio-temporal oscillations in A + B → C systems

The emergence of self-organized behaviors such as spatio-temporal oscillations is well-known for complex reactions involving nonlinear chemical or thermal feedback. Recently, it was shown that local oscillations of the chemical species concentration can be induced under isothermal batch conditions for simple bimolecular A + B → C reactions, provided they are actively coupled with hydrodynamics. When two reactants A and B, initially separated in space, react upon diffusive contact, damped spatio-temporal oscillations could develop when the surface tension increases sufficiently in the reaction zone. Additionally, if the density decreases, the coupling of both surface tension- and buoyancy-driven contributions to the flow can further sustain this oscillatory instability. Here, we investigate the opposite case of a reaction inducing a localized decrease in surface tension and an increase in density in the reacting zones. In this case, the competition arising from the two antagonistic flows is needed to create oscillatory dynamics, i.e., no oscillations are observed for pure chemically driven Marangoni flows. We study numerically these scenarios in a 2-dimensional system and show how they are controlled by the following key parameters: (i) ΔM and ΔR governing the surface tension and density variation during the reaction, respectively, (ii) the layer thickness of the system, and (iii) its lateral length. This work is a further step toward inducing and controlling chemical oscillations in simple reactions.

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.

Chemo-hydrodynamic pulsations in simple batch A + B → C systems

Spatio-temporal oscillations can be induced under batch conditions with ubiquitous bimolecular reactions in the absence of any nonlinear chemical feedback, thanks to an active interplay between the chemical process and chemically driven hydrodynamic flows. When two reactants A and B, initially separated in space, react upon diffusive contact, they can power convective flows by inducing a localized variation of surface tension and density at the mixing interface. These flows feedback with the reaction-diffusion dynamics, bearing damped or sustained spatio-temporal oscillations of the concentrations and flow field. By means of numerical simulations, we detail the mechanism underlying these chemohydrodynamic oscillations and classify the main dynamical scenarios in the relevant space drawn by parameters ΔM and ΔR, which rule the surface tension- and buoyancy-driven contributions to convection, respectively. The reactor height is found to play a critical role in the control of the dynamics. The analysis reveals the intimate nature of these oscillatory phenomena and the hierarchy among the different phenomena at play: oscillations are essentially hydrodynamic and the chemical process features the localized trigger for Marangoni flows unstable toward oscillatory instabilities. The characteristic size of Marangoni convective rolls mainly determines the critical conditions and properties of the oscillations, which can be further tuned or suppressed by the buoyancy competition. We finally discuss the possible experimental implementation of such a class of chemo-hydrodynamic oscillator and its implications in fundamental and applied terms.

Dissipative structures: From reaction-diffusion to chemo-hydrodynamic patterns

The interplay of reaction and diffusion processes can trigger localized spatiotemporal patterns when two solutions containing separate reactants A and B of an oscillating reaction are put in contact. Using the Brusselator, a classical model for chemical oscillations, we show numerically that localized waves and Turing patterns as well as reaction-diffusion (RD) patterns due to an interaction between these two kinds of modes can develop in time around the reactive contact zone depending on the initial concentration of reactants and diffusion coefficients of the intermediate species locally produced. We further explore the possible hydrodynamic destabilization of an initially buoyantly stable stratification of such an A + B → oscillator system, when the chemical reaction provides a buoyant periodic forcing via localized density changes. Guided by the properties of the underlying RD dynamics, we predict new chemo-hydrodynamic instabilities on the basis of the dynamic density profiles which are here varying with the concentration of one of the intermediate species of the oscillator. Nonlinear simulations of the related reaction-diffusion-convection equations show how the active coupling between the localized oscillatory kinetics and buoyancy-driven convection can induce pulsatile convective fingering and pulsatile plumes as well as rising or sinking Turing spots, depending on the initial concentration of the reactants and their contribution to the density.

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.

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.

Segmented waves in a reaction-diffusion-convection system

The interaction of traveling waves, with both Marangoni and buoyancy driven flows, can generate an extraordinary rich array of patterns ranging from stationary structures to chaotic waves. However, the inherent complexity of reaction-diffusion-convection (RDC) systems makes the explanation of the patterning mechanisms very difficult, both numerically and experimentally. In this paper, we describe the appearance of segmented waves in a shallow layer of an excitable Belousov-Zhabotinsky solution. The segmentation process was found to be dependent both on the depth of the solution and on the excitability of the reaction. We caught the essential features of the system through a RDC model, where the chemical waves were coupled both with surface and bulk fluid motions and we found that by varying the excitability of the reaction, and in turn the wavelength of the chemical fronts, it is possible to create a sort of hydrodynamic resonance structures (corridors), which are responsible for the segmentation process.

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.