Convective mixing induced by acid-base reactions

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

Investigation of chemoconvection in vibration fields

This study is devoted to the investigation of the chemoconvection in a two-layer miscible system caused by the neutralisation reaction proceeding in the convective-controlled (CC) regime under the influence of vertical vibrations. The CC regime without vibrational influence is characterized by the development of a density wave and vigorous convection in the upper layer, ensuring a high reaction rate and forcing the reaction front to move downwards more rapidly than in the well-known diffusive-controlled (DC) regime. It is shown that vibrations lead to some deceleration of the convection that depends both on the magnitude of the vibrational acceleration and on the initial concentrations of the reagents. Analysis of the system behaviour depending on the dimensionless parameters is carried out. It is demonstrated that the theory of thermal vibrational convection may be applied for reacting systems on quasi-steady time intervals.

Hydrodynamically-enhanced transfer of dense non-aqueous phase liquids into an aqueous reservoir

The use of surfactants represents a viable strategy to boost the removal yield of Dense Non-Aqueous Phase Liquids (DNAPLs) from groundwater and to shorten the operational timing of the remediation process. Surfactants, in general, help in reducing the interfacial tension at the DNAPL/water interface and enhance the solubility of the pollutant in the water phase through the formation of dispersed systems, such as micelles and emulsions. In this paper, we show that a suitable choice of a surfactant, in this case belonging to the bio-degradable class of ethoxylated alcohols, allows for the formation of hydrodynamic interfacial instabilities that further enhances the dissolution rate of the organic pollutant into the water phase. In a stratified configuration (denser organic phase at the bottom and lighter water phase on top), the instabilities appear as upward-pointing fingers that originate from the inversion of the local density at the interface. This inversion stems from the synergetic coupling of two effects promoted by the ethoxylated surfactant: i) the enhanced co-solubility of the DNAPL into the water (and viceversa), and (ii) the differential diffusion of the DNAPL and the surfactant in the aqueous phase. By dissolving into the DNAPL, the surfactant also reduces locally the surface tension at the liquid-liquid interface, thereby inducing transversal Marangoni flows. In our work, we carefully evaluated the effects of the concentration of different surfactants (two different ethoxylated alcohols, sodium dodecylsulphate, cetyltrimethyl ammonium bromide, N-tetradecyl-N, N-dimethylamine oxide and bis(2-ethylhexyl) sulfosuccinate sodium salt) on the onset of the instabilities in 3 different DNAPLs/water stratifications, namely chloroform, trichloroethylene and tetrachloroethylene, with a special emphasis on the trichloroethylene/water system. By means of a theoretical model and nonlinear simulations, supported by surface tension, density and diffusivity measurements, we could provide a solid explanation to the observed phenomena and we found that the type of the dispersed system, the solubility of the DNAPL into the water phase, the solubility of the surfactant in the organic phase, as well as the relative diffusion and density of the surfactant and the DNAPL in the aqueous phase, are all key parameters for the onset of the instabilities. These results can be exploited in the most common remediation techniques.

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.

Theoretical Study on Ionization of Boric Acid in Aqueous Solution by Ab Initio and DFT Methods at T=298.15 K

Objectives

The aim of this research work was to theoretically calculate the pK_a value of boric acid in aqueous solution by theoretical methods at T=298.15 K.

Materials and Methods

Boric acid has antifungal and antiviral properties. It is used in various prescription pharmaceutical products. The ab initio and density functional theory (DFT) methods were used in this research work.

Results

To explain the determined acidic dissociation constant, the various molecular conformations and solute-solvent interactions of the species of boric acid were considered. The basis set at the B3LYP/6-31+G (d) level of theory was selected for DFT calculations. We analyzed the formation of intermolecular hydrogen bonds between several species of boric acid and water molecules through Tomasi's method.

Conclusion

The result showed that there was comparable agreement between the experimentally and theoretically determined pK_a values for boric acid.

Unpredictable Dynamics of Polymeric Reacting Flow by Comparison between Pre- and Post-Reaction Fluid Properties: Hydrodynamics Involving Molecular Diagnosis via ATR-FTIR Spectroscopy

In reacting flows, changes in fluid properties induced by the chemical reaction can alter the flow dynamics. Generally, these changes in fluid properties are evaluated by comparison between their pre- and post-reaction properties. If a fluid property such as viscosity decreases between pre- and post-reaction, we expect a decrease in viscosity to occur in the reacting flow. However, this study demonstrates a reacting polymeric liquid flow where a remarkable increase in the viscoelasticity temporally occurs despite the viscosity slightly decreasing after the reaction. We elucidated the underlying reaction mechanism, which involves a structural change in the side functional group (carboxyl) in polyacrylamide at ultrahigh molecular weights ( M_w > 10⁶) with ultralow concentrations ([polymer] < 1 wt %) by using ATR-FTIR spectroscopy. This study demonstrates the existence of a reacting flow in which examination of microscopic molecular structure is required to understand the macroscopic flow dynamics. The findings will be valuable not only for industrial application such as reactor designs and rheology control but also for opening a new research area: chemically reacting flow involving the diagnosis of molecule structure.

Lateral movements in Rayleigh-Taylor instabilities due to frontiers. Numerical analysis

Numerical simulations were performed for Rayleigh-Taylor (RT) hydrodynamic instabilities when a frontier is present. The frontier formed by the interface between two fluids prevents the free movement of the fingers created by the instability. As a consequence, transversal movements at the rear of the fingers are observed in this area. These movements produce collapse of the fingers (two or more fingers join in one finger) or oscillations in the case that there is no collapse. The transversal velocity of the fingers, the amplitude of the oscillations, and the wave number of the RT instabilities as a function of the Rayleigh number (Ra) were studied near the frontier. We verified numerically that in classical RT instabilities, without a frontier, these lateral movements do not occur; only with a physical frontier, the transversal displacements of the fingers appear. The transverse displacement velocity and the initial wave number increase with Ra. This leads to the collapse of the fingers, diminishing the wave number of the instabilities at the interface. Instead, no significant changes in the amplitude of the oscillations are observed modifying Ra. The numerical results are independent of the type or origin of the frontier (gas-liquid, liquid-liquid, or solid-liquid). The numerical results are in good agreement with the experimental results reported by Binda et al. [Chaos 28, 013107 (2018)]. Based on these results, it was possible to determine the cause of the transverse displacements, which had not been explained until now.

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.

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

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.

Cross-diffusion-driven hydrodynamic instabilities in a double-layer system: General classification and nonlinear simulations

Cross diffusion, whereby a flux of a given species entrains the diffusive transport of another species, can trigger buoyancy-driven hydrodynamic instabilities at the interface of initially stable stratifications. Starting from a simple three-component case, we introduce a theoretical framework to classify cross-diffusion-induced hydrodynamic phenomena in two-layer stratifications under the action of the gravitational field. A cross-diffusion-convection (CDC) model is derived by coupling the fickian diffusion formalism to Stokes equations. In order to isolate the effect of cross-diffusion in the convective destabilization of a double-layer system, we impose a starting concentration jump of one species in the bottom layer while the other one is homogeneously distributed over the spatial domain. This initial configuration avoids the concurrence of classic Rayleigh-Taylor or differential-diffusion convective instabilities, and it also allows us to activate selectively the cross-diffusion feedback by which the heterogeneously distributed species influences the diffusive transport of the other species. We identify two types of hydrodynamic modes [the negative cross-diffusion-driven convection (NCC) and the positive cross-diffusion-driven convection (PCC)], corresponding to the sign of this operational cross-diffusion term. By studying the space-time density profiles along the gravitational axis we obtain analytical conditions for the onset of convection in terms of two important parameters only: the operational cross-diffusivity and the buoyancy ratio, giving the relative contribution of the two species to the global density. The general classification of the NCC and PCC scenarios in such parameter space is supported by numerical simulations of the fully nonlinear CDC problem. The resulting convective patterns compare favorably with recent experimental results found in microemulsion systems.

Concentration-dependent diffusion instability in reactive miscible fluids

We report on chemoconvective pattern formation phenomena observed in a two-layer system of miscible fluids filling a vertical Hele-Shaw cell. We show both experimentally and theoretically that the concentration-dependent diffusion coupled with frontal acid-base neutralization can give rise to the formation of a local unstable zone low in density, resulting in a perfectly regular cell-type convective pattern. The described effect gives an example of yet another powerful mechanism which allows the reaction-diffusion processes to govern the flow of reacting fluids under gravity conditions.

Buoyancy-driven instabilities around miscible A+B→C reaction fronts: a general classification

Upon contact between miscible solutions of reactants A and B along a horizontal interface in the gravity field, various buoyancy-driven instabilities can develop when an A+B→C reaction takes place and the density varies with the concentrations of the various chemicals. To classify the possible convective instability scenarios, we analyze the spatial dependence of the large time asymptotic density profiles as a function of the key parameters of the problem, which are the ratios of diffusion coefficients and of solutal expansion coefficients of species A, B, and C. We find that 62 different density profiles can develop in the reactive problem, whereas only 6 of them can be obtained in the nonreactive one.

Nonmonotonic Rayleigh-Taylor instabilities driven by gas-liquid CO2 chemisorption

Density variations induced by gas absorption in reactive aqueous solutions often trigger buoyancy-induced motions, generally in the form of plumes monotonically sinking into the bulk liquid and enhancing the absorption rate. Here, we contrast two types of CO2-absorbing alkaline solutions, studying their dynamics inside a vertical Hele-Shaw cell by interferometry. While the first one indeed behaves as expected, the second one leads to a quite unusual oscillatory (phase-slipping) dynamics of convective plumes, which moreover does not lead to a significant transfer enhancement. Thanks to a simplified model of momentum and species transport, we show that this particular dynamics is related to a nonmonotonic density stratification, resulting in a stagnant layer close to the interface. Conditions for this to occur are highlighted in terms of the ratios of species' diffusivities and their contribution to density, a classification deemed to be useful for optimizing chemisorption (e.g., for CO2 capture or sequestration) processes.

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.