Mixing and reaction fronts in laminar flows

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.

Influence of fast advective flows on pattern formation of Dictyostelium discoideum

We report experimental and numerical results on pattern formation of self-organizing Dictyostelium discoideum cells in a microfluidic setup under a constant buffer flow. The external flow advects the signaling molecule cyclic adenosine monophosphate (cAMP) downstream, while the chemotactic cells attached to the solid substrate are not transported with the flow. At high flow velocities, elongated cAMP waves are formed that cover the whole length of the channel and propagate both parallel and perpendicular to the flow direction. While the wave period and transverse propagation velocity are constant, parallel wave velocity and the wave width increase linearly with the imposed flow. We also observe that the acquired wave shape is highly dependent on the wave generation site and the strength of the imposed flow. We compared the wave shape and velocity with numerical simulations performed using a reaction-diffusion model and found excellent agreement. These results are expected to play an important role in understanding the process of pattern formation and aggregation of D. discoideum that may experience fluid flows in its natural habitat.

Avalanches dynamics in reaction fronts in disordered flows

We report on numerical studies of avalanches of an autocatalytic reaction front in a porous medium. The front propagation is controlled by an adverse flow resulting in upstream, static, or downstream regimes. In an earlier study focusing on front shape, we identified three different universality classes associated with this system by following the front dynamics experimentally and numerically. Here, using numerical simulations in the vicinity of the second-order transition, we identify an avalanche dynamics characterized by power-law distributions of avalanche sizes, durations, and lateral extensions. The related exponents agree well with the quenched-Kardar-Parisi-Zhang theory, which describes the front dynamics. However, the geometry of the propagating front differs slightly from that of the theoretical one. We show that this discrepancy can be understood in terms of the nonquasistatic correction induced by the finite front velocity.

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.

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.

Propagation velocities of chemical reaction fronts advected by Poiseuille flow

Poiseuille flow between parallel plates advects chemical reaction fronts, distorting them and altering their propagation velocities. Analytical solutions of the cubic reaction-diffusion-advection equation resolve the chemical concentration for narrow gaps, wide gaps, and small-amplitude flow. Numerical solutions supply a general description for fluid flow in the direction of propagation of the chemical reaction front, and for flow in the opposite direction. Empirical relations for the velocity agree with numerical solutions to within a few percent, and agree exactly with the analytical limits. Applications to nonlinear fingering are discussed.

Diffusion-limited patterning of molecules in nanofluidic channels

Diffusion-limited patterning (DLP) is a new technique that enables patterning of labile molecular species in solution phase onto surfaces that are not easily accessible. This technique is self-aligning and is simple to implement for patterning multiple species. We demonstrated DLP by patterning alternating bands of fluorescently labeled and unlabeled streptavidin in biotin-functionalized nanofluidic channels with spatial resolution better than 1 microm. The methodology of DLP also enables experimental measurement of a unique parameter that relates molecular surface grafting density, concentration, diffusivity, and channel geometry.