Concentration-dependent diffusion instability in reactive miscible fluids

Nonlinear development of convective patterns driven by a neutralization reaction in immiscible two-layer systems

This article provides the results of a theoretical and experimental study of buoyancy-driven instabilities triggered by a neutralization reaction in an immiscible two-layer system placed in a vertical Hele-Shaw cell. Flow patterns are predicted by a reaction-induced buoyancy number [Formula: see text], which we define as the ratio of densities of the reaction zone and the lower layer. In experiments, we observed the development of cellular convection ([Formula: see text]), the fingering process with an aligned line of fingertips at a slightly denser reaction zone ([Formula: see text]) and the typical Rayleigh-Taylor convection for [Formula: see text]. A mathematical model includes a set of reaction-diffusion-convection equations written in the Hele-Shaw approximation. The model's novelty is that it accounts for the water produced during the reaction, a commonly neglected effect. The persisting regularity of the fingering during the collapse of the reaction zone is explained by the dynamic release of water, which compensates for the heavy fluid falling and stabilizes the pattern. Finally, we present a stability map on the plane of the initial concentrations of solutions. Good agreement between the experimental data and theoretical results is observed. This article is part of the theme issue 'New trends in pattern formation and nonlinear dynamics of extended systems'.

Switching Modes of Mixing Due to an Adjustable Gap in a Continuous-Flow Microreactor

Microreactors are an important development in chemical engineering since the pharmaceutical industry needs flexible production rather than a large amount of product yield. The size of the microreactor may be so small that it requires the development of non-mechanical methods for reagent mixing. In this paper, we propose the design of a continuous-flow microreactor in the form of a narrow cell with a variable gap. By tuning the gap width in time and space, one can control the reaction rate and regulate the product yield. We show that the governing equation for the fluid flow can be reduced to the Darcy equation with permeability varying in space and time. As a test reaction, we consider the neutralization of nitric acid with sodium hydroxide resulting in the solutal convection in the presence of gravity. We show numerically that the prototyping spatially-distributed relief of the reactor walls can successfully separate the incoming and outgoing flows of reagents, control the mixing intensity, increase or decrease the product yield. We demonstrate also the dynamic control of the reactor efficiency via real-time local changes in the gap width.

Spatial analog of the two-frequency torus breakup in a nonlinear system of reactive miscible fluids

We present a theoretical study on pattern formation occurring in miscible fluids reacting by a second-order reaction A+B→C in a vertical Hele-Shaw cell under constant gravity. We have recently reported that the concentration-dependent diffusion of species coupled with a frontal neutralization reaction can produce a multilayer system where low-density depleted zones could be embedded between the denser layers. This leads to the excitation of chemoconvective modes spatially separated from each other by a motionless fluid. In this Rapid Communication, we show that the layers can interact via a diffusion mechanism. Since diffusively coupled instabilities initially have different wavelengths, this causes a long-wave modulation of one pattern by another. We have developed a mathematical model which includes a system of reaction-diffusion-convection equations. The linear stability of a transient base state is studied by calculating the growth rate of the Lyapunov exponent for each unstable layer. Numerical simulations supported by phase portrait reconstruction and Fourier spectra calculation have revealed that nonlinear dynamics consistently passes through (i) a perfect spatially periodic system of chemoconvective cells, (ii) a quasiperiodic system of the same cells, and (iii) a disordered fingering structure. We show that in this system, the coordinate codirected to the reaction front paradoxically plays the role of time, time itself acts as a bifurcation parameter, and a complete spatial analog of the two-frequency torus breakup is observed.

Experimental methods for visualization of hydrodynamic instability caused by the neutralization reaction in a miscible two-layer system

This paper focuses on experimental investigations of the spatio-temporal distributions of fluid velocity and temperature and the concentration of reagents and reaction products. We study concentration-dependent diffusion (CDD) convection driven by the neutralization reaction in a two-layer miscible system in a vertical Hele-Shaw cell using the original experimental complex. A comprehensive understanding of the physical mechanisms of convective motion and instabilities requires employing various experimental methods simultaneously. The proposed experimental complex provides simultaneous visualization and facilitates identification of the location of the reaction front, which is of importance to the study of its characteristics.

Shock-wave-like structures induced by an exothermic neutralization reaction in miscible fluids

We report shock-wave-like structures that are strikingly different from previously observed fingering instabilities, which occur in a two-layer system of miscible fluids reacting by a second-order reaction A+B→S in a vertical Hele-Shaw cell. While the traditional analysis expects the occurrence of a diffusion-controlled convection, we show both experimentally and theoretically that the exothermic neutralization reaction can also trigger a wave with a perfectly planar front and nearly discontinuous change in density across the front. This wave propagates fast compared with the characteristic diffusion times and separates the motionless fluid and the area with anomalously intense convective mixing. We explain its mechanism and introduce a new dimensionless parameter, which allows to predict the appearance of such a pattern in other systems. Moreover, we show that our governing equations, taken in the inviscid limit, are formally analogous to well-known shallow-water equations and adiabatic gas flow equations. Based on this analogy, we define the critical velocity for the onset of the shock wave which is found to be in the perfect agreement with the experiments.

Chemical control of dissolution-driven convection in partially miscible systems: nonlinear simulations and experiments

Numerical simulations combined with experimental results from two laboratory-scale model systems show how to control convective dissolution by chemical reactions.

Chemo-hydrodynamic patterns in porous media

Chemical reactions can interplay with hydrodynamic flows to generate chemo-hydrodynamic instabilities affecting the spatio-temporal evolution of the concentration of the chemicals. We review here such instabilities for porous media flows. We describe the influence of chemical reactions on viscous fingering, buoyancy-driven fingering in miscible systems, convective dissolution as well as precipitation patterns. Implications for environmental systems are discussed.This article is part of the themed issue 'Energy and the subsurface'.

Density profiles around A+B→C reaction-diffusion fronts in partially miscible systems: A general classification

Various spatial density profiles can develop in partially miscible stratifications when a phase A dissolves with a finite solubility into a host phase containing a dissolved reactant B. We investigate theoretically the impact of an A+B→C reaction on such density profiles in the host phase and classify them in a parameter space spanned by the ratios of relative contributions to density and diffusion coefficients of the chemical species. While the density profile is either monotonically increasing or decreasing in the nonreactive case, reactions combined with differential diffusivity can create eight different types of density profiles featuring up to two extrema in density, at the reaction front or below it. We use this framework to predict various possible hydrodynamic instability scenarios inducing buoyancy-driven convection around such reaction fronts when they propagate parallel to the gravity field.