Charles N. Baroud,Fridolin Okkels,Laure Menetrier,Patrick Tabeling
Abstract:We confront, quantitatively, the theoretical description of the reaction-diffusion of a second order reaction to experiment. The reaction at work is \ca/CaGreen, and the reactor is a T-shaped microchannel, 10 $\mu$m deep, 200 $\mu$m wide, and 2 cm long. The experimental measurements are compared with the two-dimensional numerical simulation of the reaction-diffusion equations. We find good agreement between theory and experiment. From this study, one may propose a method of measurement of various quantities, such as the kinetic rate of the reaction, in conditions yet inaccessible to conventional methods.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **How to verify the effectiveness of the two - dimensional reaction - diffusion equation in describing the chemical reaction kinetics in microfluidic reactors through quantitative comparison between experiments and theories, and propose a new method for measuring the chemical reaction rate constant**.
Specifically, the authors focus on the second - order irreversible reaction between calcium ions (Ca²⁺) and the fluorescent probe CaGreen. They use a T - shaped microchannel as a reactor, measure the concentration distribution of reaction products by fluorescence microscopy, and compare it with the results of numerical simulations. In this way, they hope to verify the accuracy of the theoretical model and propose a method for measuring the chemical reaction rate constant, especially on the sub - millisecond time scale that is difficult to achieve with traditional methods.
### Main problems and research objectives:
1. **Verify the theoretical model**: Verify through experiments whether the two - dimensional reaction - diffusion equation can accurately describe the chemical reaction process in microfluidic reactors.
2. **Measure the reaction rate constant**: Propose an experimental method based on microfluidic reactors for measuring the rate constant of chemical reactions, especially for fast reactions on the sub - millisecond time scale.
3. **Understand the morphology of the reaction zone**: Study the morphological characteristics of the reaction zone (such as the position of the maximum concentration, width, etc.), and explore the relationship between these characteristics and the initial concentration of reactants and the diffusion coefficient.
### Research background:
The reaction - diffusion equation has existed for more than a century, but it was not until recently that there were analytical solutions for step - initial conditions in an infinite domain. These analytical solutions reveal that the position of the maximum production rate in the reaction zone grows with time \( t^{1/2} \), while the width of the reaction zone grows with time \( t^{1/6} \). However, in actual experiments, it is usually affected by convection, so it is necessary to conduct research in small capillaries without large - scale flow. In addition, the traditional stopped - flow technique cannot measure the reaction rate on the sub - millisecond time scale.
### Experimental design:
The authors used a T - shaped microchannel reactor with a channel depth of 10 μm, a width of 200 μm, and a length of 2 cm. Two reactants (Ca²⁺ and CaGreen) are introduced from the two inlets of the T - shaped channel respectively and react in the main channel. The concentration distribution of reaction products is measured by fluorescence microscopy and compared with the numerical simulation results.
### Key formulas:
The reaction - diffusion equation can be expressed as:
\[
U \frac{\partial A}{\partial X} = D_a \frac{\partial^2 A}{\partial Y^2} - kAB
\]
\[
U \frac{\partial B}{\partial X} = D_b \frac{\partial^2 B}{\partial Y^2} - kAB
\]
where \( A \) and \( B \) are the concentrations of reactants A and B respectively, \( D_a \) and \( D_b \) are their diffusion coefficients, and \( k \) is the chemical reaction rate constant.
For the reaction product C, its evolution equation is:
\[
U \frac{\partial C}{\partial X} = D_c \frac{\partial^2 C}{\partial Y^2} + kAB
\]
By introducing dimensionless parameters, these equations can be simplified:
\[
\ell^2 = \sqrt{\frac{D_a D_b}{k \sqrt{A_0 B_0}}}
\]
The finally obtained dimensionless equations are:
\[
\frac{\partial a}{\partial x} = \chi \frac{\partial^2 a}{\partial y^2} - ab
\]
\[
\frac{\partial b}{\partial x} = \frac{1}{\chi} \frac{\partial^2 b}{\partial y^2} - \beta ab
\]
\[
\frac{\partial c}{\partial x} = \frac{1}{\chi} \frac{\partial^2 c}{\partial y^2} + \beta ab
\]