On the classification of buoyancy-driven chemo-hydrodynamic instabilities of chemical fronts

Thermal and compositional driven convection in thin reaction fronts

Chemical reaction fronts separate regions of reacted and unreacted substances as they propagate in liquids. These fronts may induce density gradients due to different chemical compositions and temperatures across the front. In this paper, we investigate buoyancy-induced convection driven by both types of gradients. We consider a thin front approximation where the normal front velocity depends only on the front curvature. This model applies for small curvature fronts independent of the specific type of chemical reaction. For density changes due only to heat variations near the front, we find that convection can take place for either upward or downward propagating fronts if density gradients are above a threshold. Convection can set in even if the fluid with lower density is above the higher density fluid. Our model consists of Navier-Stokes equations coupled to the front propagation equation. We carry out a linear stability analysis to determine the parameters for the onset of convection. We study the nonlinear front propagation for liquids confined in narrow two-dimensional domains. Convection leads to fronts of steady shape, propagating with constant velocities.

Front instabilities in the presence of convection due to thermal and compositional gradients

Reaction fronts separate fluids of different densities due to thermal and compositional gradients that may lead to convection. The stability of convectionless flat fronts propagating in the vertical direction depends not only on fluid properties but also in the dynamics of a front evolution equation. In this work, we analyze fronts described by the Kuramoto-Sivashinsky (KS) equation coupled to hydrodynamics. Without density gradients, the KS equation has a flat front solution that is unstable to perturbations of long wavelengths. Buoyancy enhances this instability if a fluid of lower density is underneath a denser fluid. In the reverse situation, with the denser fluid underneath, the front can be stabilized with appropriate thermal and compositional gradients. However, in this situation, a different instability develops for large enough thermal gradients. We also solve numerically the nonlinear KS equation coupled to the Navier-Stokes equations to analyze the front propagation in two-dimensional rectangular domains. As convection takes place, the reaction front curves, increasing its velocity.

Interfacial Mass Transfer in Trichloroethylene/Surfactants/ Water Systems: Implications for Remediation Strategies

The fate of dense non-aqueous phase liquids (DNAPLs) in the environment and the consequential remediation problems have been intensively studied over the last 50 years. However, a scarce literature is present about the mass transfer at the DNAPL/water interface. In this paper, we present a fast method for the evaluation of the mass transfer performance of a surfactant that can easily be employed to support an effective choice for the so-called enhanced remediation strategies. We developed a lab-scale experimental system modelled by means of simple ordinary differential equations to calculate the mass transfer coefficient (K) of trichloroethylene, chosen as representative DNAPL, in the presence and in the absence of two ethoxylated alcohols belonging to the general class of Synperonic surfactants. Our findings revealed that it exists an optimal surfactant concentration range, where K increases up to 40% with respect to pure water.

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.

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.

Thermal effects on the diffusive layer convection instability of an exothermic acid-base reaction front

A buoyancy-driven hydrodynamic instability appearing when an aqueous acid solution of HCl overlies a denser alkaline aqueous solution of NaOH in a vertically oriented Hele-Shaw cell is studied both experimentally and theoretically. The peculiarity of this reactive convection pattern is its asymmetry with regard to the initial contact line between the two solutions as convective plumes develop in the acidic solution only. We investigate here by a linear stability analysis (LSA) of a reaction-diffusion-convection model of a simple A+B→C reaction the relative role of solutal versus thermal effects in the origin and location of this instability. We show that heat effects are much weaker than concentration-related ones such that the heat of reaction only plays a minor role on the dynamics. Computation of density profiles and of the stability analysis eigenfunctions confirm that the convective motions result from a diffusive layer convection mechanism whereby a locally unstable density stratification develops in the upper acidic layer because of the difference in the diffusion coefficients of the chemical species. The growth rate and wavelength of the pattern are determined experimentally as a function of the Brinkman parameter of the problem and compare favorably with the theoretical predictions of both LSA and nonlinear simulations.

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.

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.

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.

Hot spots in density fingering of exothermic autocatalytic chemical fronts

Measurements of two-dimensional (2D) temperature fields are performed by an interferometric method during density fingering of the autocatalytic chlorite-tetrathionate reaction in a Hele-Shaw cell. These measures confirm that, because of heat losses through the glass walls of the reactor, the temperature profile across the front is a pulse rather than a front. Moreover, the full 2D temperature field shows the presence in the reactive zone of hot spots where the temperature exceeds the maximum temperature measured in a stable planar front. We investigate here experimentally the increase of temperature in the hot spots when the composition of the reactants is varied to increase the exothermicity of the reaction. We back up these experimental observations by nonlinear simulations of a reaction-diffusion-convection model which show that the maximum temperature reached in the system depends on the intensity of convection.

Spatiotemporal chaos in the dynamics of buoyantly and diffusively unstable chemical fronts

Nonlinear dynamics resulting from the interplay between diffusive and buoyancy-driven Rayleigh-Taylor (RT) instabilities of autocatalytic traveling fronts are analyzed numerically for various values of the relevant parameters. These are the Rayleigh numbers of the reactant A and autocatalytic product B solutions as well as the ratio D=D(B)/D(A) between the diffusion coefficients of the two key chemical species. The interplay between the coarsening dynamics characteristic of the RT instability and the constant short wavelength modulation of the diffusive instability can lead in some regimes to complex dynamics dominated by irregular succession of birth and death of fingers. By using spectral entropy measurements, we characterize the transition between order and spatial disorder in this system. The analysis of the power spectrum and autocorrelation function, moreover, identifies similarities between the various spatial patterns. The contribution of the diffusive instability to the complex dynamics is discussed.

Steady dynamics of exothermic chemical wave fronts in van der Waals fluids

We study the steady dynamics of an exothermic Fisher-Kolmogorov-Petrovsky-Piskunov chemical wave front traveling in a one-dimensional van der Waals fluid. The propagating wave is initiated by a nonuniformity in reactant concentration contrary to usual combustion ignition processes. The heat release and activation energy of the reaction play the role of control parameters. We recently proved that the propagation of an exothermic chemical wave front in a perfect gas displays a forbidden interval of stationary wave front speeds [G. Dumazer, M. Leda, B. Nowakowski, and A. Lemarchand, Phys. Rev. E 78, 016309 (2008)]. We examine how this result is modified for nonideal fluids and determine the effect of the van der Waals parameters and fluid density on the bifurcation between diffusion flames and Chapman-Jouguet detonation waves as heat release increases. Analytical predictions are confirmed by the numerical solution of the hydrodynamic equations including reaction kinetics.

Measurement of the temperature profile of an exothermic autocatalytic reaction front

Autocatalytic reactions may propagate as solitary waves, namely, at a constant front velocity and with a stationary concentration profile, resulting from a balance between molecular diffusion and chemical reaction. When the reaction is exothermic, a thermal wave is linked to the chemical front. As the thermal diffusivity is nearly two orders of magnitude larger than the molecular one, the temperature profile spreads over length scales (mm) two orders of magnitude larger than the concentration one. Using an infrared camera, we measure the temperature profiles for a chlorite-tetrathionate autocatalytic reaction. The profiles are compared quantitatively to lattice Bhatnagar-Gross-Krook (BGK) numerical simulations. Our analysis also accounts for the lack of observation of the thermal wave for the iodate arsenous acid reaction.

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. 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.

Miscible viscous fingering induced by a simple A+B-->C chemical reaction

Viscous fingering (VF) is a hydrodynamical instability that occurs in porous media when a less viscous fluid displaces a more viscous one. We investigate here numerically how such an instability can be triggered by a simple A+B-->C reaction when a solution of one reactant is displacing linearly a miscible solution of another reactant of same viscosity producing a more viscous product C at the interface. The properties of the fingering pattern observed in the zone where the less viscous reactant pushes more viscous products are studied as a function of the relevant parameters of the problem. These are the Damköhler number, the viscosity contrast between reactants and product, the ratio of initial concentrations of A and B , and the diffusion coefficients of each species. Our study shows that the fingering pattern can in some cases be different whether A displaces B or vice versa, enlightening recent experimental observations of such asymmetries in micellar systems. In particular, we show that in this asymmetric case, VF is more intense when the invading chemical solution is either the less concentrated one for equal diffusivities or when it contains the slower diffusing reactant for fixed equimolar initial concentrations.

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.

Hot spots revealed by simultaneous experimental measurement of the two-dimensional concentration and temperature fields of an exothermic chemical front during finger-pattern formation

A noninvasive optical technique combining digital interferometry in transmission and transparency measurement of concentration is developed to analyze spatiotemporal dynamics of physicochemical systems. This technique allows one to measure simultaneously the two-dimensional (2D) dynamics of concentration and temperature fields in both reactive and nonreactive systems contained inside a transparent cell. When used to experimentally analyze buoyancy-driven fingering of an exothermic autocatalytic chemical front, this method reveals in the 2D temperature field the presence of hot spots where the temperature locally exceeds the adiabatic one.

Transition between an exothermic chemical wave front and a generic flame

We consider the two classes of exothermic chemical wave fronts, propagating toward a stable or an unstable steady state. The hydrodynamic equations for stream velocity, temperature, and concentrations are solved numerically for increasing values of the reaction heat. For a critical value of the heat release, we find a transition between a chemical front, whose speed depends on the chemical dynamics, and a generic flame, whose speed is entirely determined by heat release. We derive an analytical expression of the flame speed from the invariants of the hydrodynamic equations. This result substantiates macroscopic approaches widely used in combustion, in which the chemical models include only simplified reaction mechanisms.

Convective chemical fronts in a Poiseuille flow

Autocatalytic reaction fronts propagating in a Poiseuille flow present a change of speed and curvature depending on the strength of the flow and on the direction of front propagation. These chemical fronts separate reacted and unreacted fluids of different densities, consequently convection will always be present due to the horizontal density gradient of the curved front. In this paper, we find the change of speed caused by gravity for fronts propagating in vertical tubes under a Poiseuille flow. For small density differences, we find axisymmetric fronts. Our theory predicts a transition to nonaxisymmetric fronts as the distance between the walls is increased. The transition depends on the average speed of the Poiseuille flow.

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)].

Front fingering and complex dynamics driven by the interaction of buoyancy and diffusive instabilities

Traveling fronts can become transversally unstable either because of a diffusive instability arising when the key variables diffuse at sufficiently different rates or because of a buoyancy-driven Rayleigh-Taylor mechanism when the density jump across the front is statically unfavorable. The interaction between such diffusive and buoyancy instabilities of fronts is analyzed theoretically for a simple model system. Linear stability analysis and nonlinear simulations show that their interplay changes considerably the stability properties with regard to the pure Rayleigh-Taylor or diffusive instabilities of fronts. In particular, an instability scenario can arise which triggers convection around statically stable fronts as a result of differential diffusion. Moreover, spatiotemporal chaos can be observed when both buoyancy and diffusive effects cooperate to destabilize the front. Experimental conditions to test our predictions are suggested.