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

Frontal polymerization in thin layers: Hydrodynamic effects and asymptotic dynamics

Buoyancy-driven convection currents arise from temperature gradients in thermal frontal polymerization (FP) when the spatially localized polymerization reaction travels perpendicularly to the gravity field. We propose a theoretical study of the system dynamics under adiabatic conditions. The polymer and the reactant mixture are considered to be in the same liquid phase, but the viscosity can increase with the degree of polymerization. We find that the reaction zone propagates as a hot spot-like pattern with a broken symmetry in both the vertical and horizontal directions. Furthermore, the system can reach an asymptotic dynamics characterized by a front with a steady shape that propagates at constant speed with a steady vortex surrounding it. As the strength of the vortex is increased, either by decreasing the reactants' viscosity or by increasing the layer's thickness, we observe a transition between (i) a passive regime predicted by pure reaction-diffusion and hydrodynamic models and (ii) an active chemo-hydrodynamic regime where such models separately break down. In the active regime (ii), the front speed decreases as convection intensifies. By means of a scaling analysis, we explain how hydrodynamic currents might lower the velocity of a polymerization wave. As the viscosity of the polymer is enlarged, the flow is shifted ahead of the reaction zone and becomes more symmetrical with respect to the middle of the system, as recently observed in solid-liquid FP experiments [Y. Gao et al., Phys. Rev. Lett. 130, 028101 (2023) and Y. Gao et al., Int. J. Heat Mass Transf. 240, 126622 (2025)].

Propagating fronts in fluids with solutal feedback

We numerically study the propagation of reacting fronts in a shallow and horizontal layer of fluid with solutal feedback and in the presence of a thermally driven flow field of counterrotating convection rolls. We solve the Boussinesq equations along with a reaction-convection-diffusion equation for the concentration field where the products of the nonlinear autocatalytic reaction are less dense than the reactants. For small values of the solutal Rayleigh number the characteristic fluid velocity scales linearly, and the front velocity and mixing length scale quadratically, with increasing solutal Rayleigh number. For small solutal Rayleigh numbers the front geometry is described by a curve that is nearly antisymmetric about the horizontal midplane. For large values of the solutal Rayleigh number the characteristic fluid velocity, the front velocity, and the mixing length exhibit square-root scaling and the front shape collapses onto an asymmetric self-similar curve. In the presence of counterrotating convection rolls, the mixing length decreases while the front velocity increases. The complexity of the front geometry increases when both the solutal and convective contributions are significant and the dynamics can exhibit chemical oscillations in time for certain parameter values. Last, we discuss the spatiotemporal features of the complex fronts that form over a range of solutal and thermal driving.

Interaction of Pure Marangoni Convection with a Propagating Reactive Interface under Microgravity

A reactive interface in the form of an autocatalytic reaction front propagating in a bulk phase can generate a dynamic contact line upon reaching the free surface when a surface tension gradient builds up due to the change in chemical composition. Experiments in microgravity evidence the existence of a self-organized autonomous and localized coupling of a pure Marangoni flow along the surface with the reaction in the bulk. This dynamics results from the advancement of the contact line at the surface that acts as a moving source of the reaction, leading to the reorientation of the front propagation. Microgravity conditions allow one to isolate the transition regime during which the surface propagation is enhanced, whereas diffusion remains the main mode of transport in the bulk with negligible convective mixing, a regime typically concealed on Earth because of buoyancy-driven convection.

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.

Chlorine dioxide-induced and Congo red-inhibited Marangoni effect on the chlorite-trithionate reaction front

Hydrodynamic flows can exert multiple effects on an exothermal autocatalytic reaction, such as buoyancy and the Marangoni convection, which can change the structure and velocity of chemical waves. Here we report that in the chlorite-trithionate reaction, the production and consumption of chlorine dioxide can induce and inhibit Marangoni flow, respectively, leading to different chemo-hydrodynamic patterns. The horizontal propagation of a reaction-diffusion-convection front was investigated with the upper surface open to the air. The Marangoni convection, induced by gaseous chlorine dioxide on the surface, produced from chlorite disproportionation after the proton autocatalysis, has the same effect as the heat convection. When the Marangoni effect is removed by the reaction of chlorine dioxide with the Congo red (CR) indicator, an oscillatory propagation of the front tip is observed under suitable conditions. Replacing CR with bromophenol blue (BPB) distinctly enhanced the floating, resulting in multiple vortexes, owing to the coexistence between BPB and chlorine dioxide. Using the incompressible Navier-Stokes equations coupled with reaction-diffusion and heat conduction equations, we numerically obtain various experimental scenarios of front instability for the exothermic autocatalytic reaction coupled with buoyancy-driven convection and Marangoni convection.

Convective flow reversal in self-powered enzyme micropumps

Surface-bound enzymes can act as pumps that drive large-scale fluid flows in the presence of their substrates or promoters. Thus, enzymatic catalysis can be harnessed for “on demand” pumping in nano- and microfluidic devices powered by an intrinsic energy source. The mechanisms controlling the pumping have not, however, been completely elucidated. Herein, we combine theory and experiments to demonstrate a previously unreported spatiotemporal variation in pumping behavior in urease-based pumps and uncover the mechanisms behind these dynamics. We developed a theoretical model for the transduction of chemical energy into mechanical fluid flow in these systems, capturing buoyancy effects due to the solution containing nonuniform concentrations of substrate and product. We find that the qualitative features of the flow depend on the ratios of diffusivities δ=D(P)/D(S) and expansion coefficients β=β(P)/β(S) of the reaction substrate (S) and product (P). If δ>1 and δ>β (or if δ<1 and δ<β ), an unexpected phenomenon arises: the flow direction reverses with time and distance from the pump. Our experimental results are in qualitative agreement with the model and show that both the speed and direction of fluid pumping (i) depend on the enzyme activity and coverage, (ii) vary with the distance from the pump, and (iii) evolve with time. These findings permit the rational design of enzymatic pumps that accurately control the direction and speed of fluid flow without external power sources, enabling effective, self-powered fluidic devices.

Three-dimensional convection-driven fronts of the exothermic chlorite-tetrathionate reaction

Horizontally propagating autocatalytic reaction fronts in fluids are often accompanied by convective motion in the presence of gravity. We experimentally and numerically investigate the stable complex three-dimensional pattern arising in the exothermic chlorite-tetrathionate reaction as a result of the antagonistic thermal and solutal contribution to the density change. By particle image velocimetry measurements, we construct the flow field that stabilizes the front structure. The calculations applied for incompressible fluids using the empirical rate-law model reproduce the experimental observations with good agreement.

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.

Relation between shape of liquid-gas interface and evolution of buoyantly unstable three-dimensional chemical fronts

Buoyantly unstable 3D chemical fronts were seen traveling through an iodate-arsenous acid reaction solution. The experiments were performed in channel reactors with rectangular cross sections, where the top of the reaction solution was in contact with air. A concave or convex meniscus was pinned to reactor lateral walls. Influence of the meniscus shape on front development was investigated. For the concave meniscus, an asymptotic shape of fronts holding negative curvature was observed. On the other hand, fronts propagating in the solution with the convex meniscus kept only positive curvature. Those fronts were also a bit faster than fronts propagating in the solution with the concave meniscus. A relation between the meniscus shape, flow distribution, velocity, and shape is discussed.

CHEMO-hydrodynamic coupling between forced advection in porous media and self-sustained chemical waves

Autocatalytic reaction fronts between two reacting species in the absence of fluid flow, propagate as solitary waves. The coupling between autocatalytic reaction front and forced simple hydrodynamic flows leads to stationary fronts whose velocity and shape depend on the underlying flow field. We address the issue of the chemico-hydrodynamic coupling between forced advection in porous media and self-sustained chemical waves. Towards that purpose, we perform experiments over a wide range of flow velocities with the well characterized iodate arsenious acid and chlorite-tetrathionate autocatalytic reactions in transparent packed beads porous media. The characteristics of these porous media such as their porosity, tortuosity, and hydrodynamics dispersion are determined. In a pack of beads, the characteristic pore size and the velocity field correlation length are of the order of the bead size. In order to address these two length scales separately, we perform lattice Boltzmann numerical simulations in a stochastic porous medium, which takes into account the log-normal permeability distribution and the spatial correlation of the permeability field. In both experiments and numerical simulations, we observe stationary fronts propagating at a constant velocity with an almost constant front width. Experiments without flow in packed bead porous media with different bead sizes show that the front propagation depends on the tortuous nature of diffusion in the pore space. We observe microscopic effects when the pores are of the size of the chemical front width. We address both supportive co-current and adverse flows with respect to the direction of propagation of the chemical reaction. For supportive flows, experiments and simulations allow observation of two flow regimes. For adverse flow, we observe upstream and downstream front motion as well as static front behaviors over a wide range of flow rates. In order to understand better these observed static state fronts, flow experiments around a single obstacle were used to delineate the range of steady state behavior. A model using the "eikonal thin front limit" explains the observed steady states.

Horizontally propagating three-dimensional chemo-hydrodynamic patterns in the chlorite-tetrathionate reaction

Planar reaction fronts resulting from the coupling of exothermic autocatalytic reactions and transport processes can be deformed by convection in the presence of gravity field. We have experimentally investigated how buoyancy affects the spatiotemporal pattern formation at various solution thicknesses in three-dimensional medium. In the chlorite-tetrathionate reaction, a stable structure propagating horizontally with constant velocity and geometry develops when appropriately thick solutions are studied. Both the horizontal and the vertical projections of the resulting three-dimensional structures are quantitatively characterized: the smooth leading edge of the front is independent of the solution thickness and the structured trailing edge ends in a center cusp with a constant angle.

The heads and tails of buoyant autocatalytic balls

Buoyancy produced by autocatalytic reaction fronts can produce fluid flows that advect the front position, giving rise to interesting feedback between chemical and hydrodynamic effects. In this paper, we numerically investigate the evolution of autocatalytic iodate-arsenous acid reaction fronts initialized in spherical configurations. Deformation of these "autocatalytic balls" is driven by buoyancy produced by the reaction. In our simulations, we have found that depending on the initial ball radius, the reaction front will develop in one of three different ways. In an intermediate range of ball size, the flow can evolve much like an autocatalytic plume: the ball develops a reacting head and tail that is akin to the head and conduit of an autocatalytic plume. In the limit of large autocatalytic balls, however, growth of a reacting tail is suppressed and the resemblance to plumes disappears. Conversely, very small balls of product solution fail to initiate sustained fronts and eventually disappear.

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.

Buoyancy-driven convection may switch between reactive states in three-dimensional chemical waves

Traveling waves in an extended reactor, whose width cannot be neglected, represent a three-dimensional (3D) reaction-diffusion-convection system. We investigate the effects of buoyancy-driven convection in such a setting. The 3D waves traveled through horizontal layers of the iodate-arsenous acid (IAA) reaction solution containing excess of arsenous acid. The depth of the reaction solution was the examined parameter. An increase in the intensity of buoyancy-driven flow caused an increase of the traveling wave velocities. Convection distorted the front of the chemical waves. For layers deeper than h>13 mm, heat release became smaller than heat production causing the emergence of Rayleigh-Bénard convection cells. At the interface, a dependency of wave shape on solution depth was observed. For h<7 mm, the waves adopted a stable V-like shape, while for h>13 mm a parabolic shape dominated. For 7<h<13 mm, both shapes were realized with the same probability. Finally, an intermittent switch between stoichiometric regimes is observed as an unexpected effect of the buoyancy-driven convection. The switch is expressed by iodine enrichment in the product. Hence, the experiments demonstrate that the buoyancy-driven convective flow can cause long-lived, but nevertheless transient, changes in the chemical composition by inducing a local transition between different regimes of the IAA reaction.

Numerical simulations of a buoyant autocatalytic reaction front in tilted Hele-Shaw cells

We present a numerical analysis of solutal buoyancy effects on the shape and the velocity of autocatalytic reaction fronts, propagating in thin tilted rectangular channels. We use two-dimensional (2D) lattice Bathnagar-Gross-Krook (BGK) numerical simulations of gap-averaged equations for the flow and the concentration, namely a Stokes-Darcy equation coupled with an advection-diffusion-reaction equation. We do observe stationary-shaped fronts, spanning the width of the cell and propagating along the cell axis. We show that the model accounts rather well for experiments we performed using an Iodate Arsenous Acid reaction propagating in tilted Hele-Shaw cells, hence validating our 2D modelization of a three-dimensional problem. This modelization is also able to account for results found for another chemical reaction (chlorite tetrathionate) in a horizontal cell. In particular, we show that the shape and the traveling velocity of such fronts are linked with an eikonal equation. Moreover, we show that the front velocity varies nonmonotonically with the tilt of the cell, and nonlinearly with the width of the cell.

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.

Convective instabilities in horizontally propagating vertical chemical fronts

The effect of convection on horizontally propagating fronts arising in the acid-catalyzed reaction of chlorite and tetrathionate ions is investigated experimentally in vertically oriented narrow reaction vessels of various heights. A stable finger develops which can be characterized by its mixing length defined as the standard deviation of the average front profile. The mixing length scales with the height of the reaction vessel, and the scaling law has been determined for three different chemical compositions.

Dynamics of A + B --> C reaction fronts in the presence of buoyancy-driven convection

The dynamics of A+B-->C fronts in horizontal solution layers can be influenced by buoyancy-driven convection as soon as the densities of A, B, and C are not all identical. Such convective motions can lead to front propagation even in the case of equal diffusion coefficients and initial concentration of reactants for which reaction-diffusion (RD) scalings predict a nonmoving front. We show theoretically that the dynamics in the presence of convection can in that case be predicted solely on the basis of the knowledge of the one-dimensional RD density profile across the front.

Nonlinear chemoconvection in the methylene-blue-glucose system: two-dimensional shallow layers

Interfacial hydrodynamic instabilities arise in a range of chemical systems. One mechanism for instability is the occurrence of unstable density gradients due to the accumulation of reaction products. In this paper we conduct two-dimensional nonlinear numerical simulations for a member of this class of system: the methylene-blue-glucose reaction. The result of these reactions is the oxidation of glucose to a relatively, but marginally, dense product, gluconic acid, that accumulates at oxygen permeable interfaces, such as the surface open to the atmosphere. The reaction is catalyzed by methylene-blue. We show that simulations help to disassemble the mechanisms responsible for the onset of instability and evolution of patterns, and we demonstrate that some of the results are remarkably consistent with experiments. We probe the impact of the upper oxygen boundary condition, for fixed flux, fixed concentration, or mixed boundary conditions, and find significant qualitative differences in solution behavior; structures either attract or repel one another depending on the boundary condition imposed. We suggest that measurement of the form of the boundary condition is possible via observation of oxygen penetration, and improved product yields may be obtained via proper control of boundary conditions in an engineering setting. We also investigate the dependence on parameters such as the Rayleigh number and depth. Finally, we find that pseudo-steady linear and weakly nonlinear techniques described elsewhere are useful tools for predicting the behavior of instabilities beyond their formal range of validity, as good agreement is obtained with the simulations.

Solitary Marangoni-driven convective structures in bistable chemical systems

Bistable chemical fronts can be deformed by Marangoni-driven convective flows induced by gradients of surface tension across the front. We investigate here the nonlinear dynamics of such a system by simulations of two-dimensional Navier-Stokes equations coupled to a reaction-diffusion-convection equation for a surface-active chemical species present in the bulk of the solution and involved in a bistable kinetics. We show that Marangoni flows cannot only alter the shape and speed of the front but also change the relative stability of the two stable steady states, reversing in some cases the direction of propagation of the front with regard to the pure reaction-diffusion situation. A detailed parametric study discusses the properties of the asymptotic dynamics as a function of the Marangoni number M and of a kinetic parameter d.