1. Introduction
The growing interest towards microreactors is triggered by the possibility to combine continuous flow operations with an unprecedented control of operating conditions. This unique feature opens the possibility to reduce reactants, wastes and energy consumption, enabling more sustainable processes. Indeed, it is well accepted that microfluidics represents the main route towards process intensification to improve the yield and the selectivity of chemical reactions while limiting their environmental impact [
1]. It has been estimated that around 50% of the homogeneous reactions and 30% of heterogeneous reactions, which are encountered in pharmaceutical and fine chemistry processes, may benefit of microfluidics [
2,
3,
4]. Moreover, microreactors give the chance to explore new and greener synthesis routes, as the huge surface-to-volume ratio ensures very high heat transfer rates, allowing to safely deal with very exothermic or too dangerous reactions [
5].
Despite the above fascinating characteristics, the actual application of microreactors has been hindered so far; the main reason is the difficulty to efficiently mix reactants [
6]. Indeed mixing occurs solely through diffusion due to the laminar regime, which is determined by the tiny dimensions. Hence many efforts have been devoted at enhancing mixing through active or passive techniques. The latter methods are appealing as they do not need any external energy source and are more robust than active ones. These passive microreactors rely on a special geometrical arrangement of the inlet and mixing channels aimed at promoting the mixing of the liquid streams through convection [
7].
The simplest example of a passive microreactor is the T-junction; despite such a basic configuration, its fluid dynamics is very intricate and exhibits different flow regimes, depending on the Reynolds number, i.e.,
[
8]. For very low Reynolds numbers, the two inlet streams entering the T-shaped microreactor remain completely segregated in the mixing channel (segregated regime) with low mixing due to diffusion. With increasing
, 3D U-shaped vortical structures stem from the confluence between the inlet streams and generate a double pair of counter-rotating vortical legs in the mixing channel (steady symmetric or vortex regime). Further augmenting the flow rate, the top part of the U-shaped vortical structures tilts and, as a consequence, one of the two legs becomes more intense than the other, which soon disappears in the channel (steady asymmetric or engulfment regime). This latter flow regime ensures a remarkable increase in the degree of mixing, because of the associated loss of symmetry. Then, the flow becomes unsteady and time-periodic, although being still in the laminar regime, preserving an intense asymmetry (unsteady asymmetric regime). The flow oscillations further push mixing, resulting in values of the degree of mixing that are even 30–50% larger than those of the engulfment regime [
9,
10,
11]. However and unexpectedly, a further increase of
triggers a flow, which, although still time-periodic, exhibits a large degree of symmetry, thus strongly hampering mixing (unsteady symmetric regime). Finally, the flow becomes chaotic.
The literature on the above flow regimes occurring in T-shaped microreactors is impressively extensive with both experimental and numerical investigations. It covers the effect of operating conditions [
12], the insight on mechanisms promoting flow regime transitions [
13], the effect of channel aspect ratios [
14,
15,
16], the effect of the inclination [
17,
18,
19,
20] and position of the inlet channels [
21,
22,
23], the presence of obstacles [
24] as well as influence of fluid properties [
10,
25,
26,
27]. Some works investigate the influence of inlet flow conditions, e.g., not fully developed flow [
28] or swirled flow [
29]) on flow regimes. Besides, a dimensionless number taking into account the T-mixer geometry, has been proposed in the literature to identify the onset of the different steady flow regimes [
30]. A comprehensive review on this topic has been recently provided by Camarri et al. [
31]. However, in most of the above investigations the focus is on how flow regimes affect mixing rather than the progress of a chemical reaction.
In our recent works [
32,
33], we attempted to fill this gap by providing an insight on the role that steady flow regimes play on the performance of a chemical reaction in a T-shaped microreactor. More specifically, we considered different kinetic constants to comprehend the interplay between flow and kinetics. We highlighted different behaviors depending on the flow regime: in the segregated regime, the yield
of a chemical reaction depends solely on the flow to chemical time-scale ratio, i.e.,
, where
is the Damköhler number, while in the vortex and engulfment regimes, an explicit dependency on the kinetic constant
was observed with
.
The present work, aims at extending the above investigation to unsteady flows. Firstly, we want to assess the presence of unsteady motions even for reactive fluids. Secondly, we aim at analyzing the impact of these motions on the progress of a chemical reaction, comparing the behavior to that of the steady regimes. To this purpose, experiments and numerical simulations are carried out jointly.
2. Experimental Setup
The T-shaped microreactor and the experimental setup are the same as in the work of [
32], therefore only a few details are briefly recalled herein. The inlet channels of the microreactor are characterized by a square cross-section, with width
1 mm, while the mixing channel has a rectangular cross-section with aspect ratio 2:1, i.e.,
mm. Hence, the hydraulic diameter of the mixing channel is
. In the following text, coordinates will be made non-dimensional by using the hydraulic diameter, i.e.,
,
, and
. The inlet channels length is equal
, so the flow is fully developed at the T-junction, whereas the mixing channel length is equal to
and assures the complete evolution of vortical structures. The geometry of the reactor and the reference system of coordinates are provided in
Figure 1a, while the experimental setup is shown in
Figure 1b.
The test reaction is the reduction at ambient temperature of methylene blue (MB
+), a water-soluble dye molecule, to the corresponding colorless leuco compound (LMB
+) by adding ascorbic acid (AsA). An additional reaction product is the dehydroascorbic acid (DA). Such reaction is catalyzed by hydrogen chloride (HCl).
The aqueous solution of
and HCl enters the reactor from one inlet while the aqueous solution of AsA is fed into the other inlet. The two reactant solutions were prepared as follows: 17 mg/L of
powder (by Sigma Aldrich) are dissolved in aqueous solutions of 1 mol/L of HCl (by Sigma Aldrich), thus obtaining a concentration of
equal to [
]
mol/L; as for the ascorbic acid solution, 300 mg/L of L-Ascorbic Acid (Ultrafine Vitamin C-powder by Cutatonic
®) are dissolved in deionized water, to obtain an ascorbic acid concentration equal to [AsA]
mol/L. A syringe pump (KD Scientific), equipped with two Becton Dickson plastic syringes of 60 mL, was used to introduce the two streams at the same flow rate, and thus equal inlet bulk velocity,
U. The pump is stopped after each test and then started again imposing the new flow rate; this procedure is aimed at avoiding any hysteresis phenomena in the mixing and reaction process.
In case of an excess of ascorbic acid, i.e., when [AsA] ≫ [
] as in the present case, the above reaction has a pseudo-first order kinetic [
34], that depends on the concentration of
only with a kinetic constant,
1/s (see [
32] for more details).
We selected this reaction because we can exploit the decolorization of the blue reactant to visualize the progress of the reaction at the confluence region and along the mixing channel, and thus study the effect of the flow regimes and mixing on the reaction yield. To this purpose, an upright microscope (Nikon, model Eclipse 80i), equipped with a magnifying lens of 4× (aperture N.A. = 0.13), a D-LH 12 V-100 W halogen lamp and a high-speed camera (Optomotive Velociraptor HS) with an up-mounted 0.7× lens, were used. Images were acquired with a resolution of 780 × 1504 pixels and an acquisition frequency of 450 frames/s. The image exposure time is significantly shorter than the one related to the frame rate, i.e., s. Each image collects a region of (the equivalent dimension of a pixel is 2.93 μm). Among the species involved in the reaction, is the only light-absorbing species. Therefore, the resulting intensity of the light across the specimen that is captured by the camera depends only on the MB+ concentration.
The Lambert-Beer’s law was used to convert light intensity into normalized depth-averaged
concentration for each pixel of the image and was implemented in an offline procedure enabling the post processing of the flow visualizations. Additional details on calibration, image post-processing, and experimental setup are provided in [
35].
3. Numerical Setup
The flow can be described by unsteady Navier-Stokes equations and transport/reaction equations for the chemical species
,
, AsA, HCl and DA. The energy equation is not resolved as thermal effects can be neglected. The non-dimensional form of the equation is:
where the lengths are normalized with the mixing channel hydraulic diameter
d and velocities with the inlet bulk velocity
U.
is the non-dimensional time, i.e.,
. In the above equations,
represents the velocity vector,
p is the modified non-dimensional pressure, i.e.,
where
P is the pressure, and
g the gravity acceleration, while
is the non-dimensional gravity, i.e.,
.
and
are non-dimensional density and viscosity, respectively, referred to a convenient reference state, i.e.,
and
, which in our case corresponds to pure water at 20 °C.
represents the mass fraction of the
k-th chemical species while
is its rate of production or consumption due to chemical reactions, which, for instance is:
for the methylene blue.
is the non-dimensional diffusivity referred to the water self-diffusivity
.
The characteristic non-dimensional numbers are:
where is the kinematic viscosity of water, while is the density difference between the two inlet fluids.
Both fluid density and viscosity depends on the ascorbic acid content
and hence, the correlations
from [
36] are implemented in the model, see for details [
32]. The fluid stream carrying AsA is 1.17 denser, at the inlet, than the other stream; indeed some stratification was observed to take place at
for the given size of our microreactor [
32]. Despite the fact that stratification is therefore negligible in the range of
for which unsteady flows occur, we retained the full formulation of the governing equations, including the term related to gravity.
The boundary conditions consist of: uniform velocity and concentration of the reactants at the entrance of the inlet channels; no-slip velocity at the channel walls and pressure outlet condition with ambient pressure at the outlet.
The transient solver of the finite volume code ANSYS Fluent v.19 [
37], was utilized to solve the above equations. Time advancement was handled with a second order implicit advancement method, by using a time step corresponding to a CFL number
, while a second order upwind interpolation scheme was used for spatial discretization. The PISO algorithm was employed for the pressure–velocity coupling. Such numerical setup has been successfully bench-marked against a massive parallel spectral element code in [
38] where details on the grid independence study are also given. The computational grid is fully structured with 4.7 M cells; the cells are cubical at the confluence region, while they elongate along the inlet and the outlet channels.
4. Results
Figure 2 summarizes the steady flow regimes, namely segregated, vortex and engulfment regimes, which have been already characterized in Mariotti et al. [
32]. The images in
Figure 2a–c represent the depth-averaged
concentration, estimated from the numerical simulations, in the confluence region. Such concentration is reported in non-dimensional form, being 1 the
concentration in the inlet.
The agreement with the experimental flow visualizations has already been discussed in the same paper. Clearly, the engulfment regime (
Figure 2c at
) breaks the flow symmetries, thus strongly enhancing mixing with respect to the segregated (
Figure 2a at
) and vortex regimes (
Figure 2b at
), which both exhibit a prevailing dark region, corresponding to the methylene blue aqueous solution, on the left and a light region, corresponding to the aqueous solution of ascorbic acid, on the right. The non-dimensional velocity magnitude, i.e., normalized with the bulk velocity, in the horizontal mid-plane is reported in
Figure 2d–f along with the flow streamlines. We observe that in both the segregated and vortex regimes the inlet streams flow side by side in the channel; however, a recirculation region is visible at the top of the microreactor. This recirculation region is more intense in the engulfment regime, where the fingerprint of the top part of the U-shaped vortical structures, which are tilted, can be well discerned from the streamlines.
The experimental visualization of the time-periodic flow occurring at
is reported in
Figure 3, which shows six different instants within a temporal cycle, as well as in the
Videos S1 and S2. For each experimental test, images were acquired firstly in the confluence region, at
, then in the region at
and further along the mixing channel
. Indeed, the flow visualization of the entire mixing channel was not undertaken as it would lead to a poor resolution because of the large length-to-width ratio of the channel. Therefore, we remark that the flow visualizations in the different regions, illustrated in
Figure 3, are not synchronized with each other.
The temporal evolution of the
concentration, confirms the occurrence and features of the asymmetric unsteady regime, that has been previously observed in literature with non-reactive systems [
39]). More specifically, at the T-junction (top series of images) the initial flow configuration resembles the one of the engulfment regime. Here the top part of the U-shaped vortical structures is tilted, but during the temporal evolution, the tilting angle increases while the vortical structures approach each other, merge and annihilate. Simultaneously, the co-rotating vortical structures in the mixing channel generate a blob of vorticity (see top of
Figure 3c), which is convected along the mixing channel (see top of
Figure 3d,e). Moving downstream in the mixing channel (middle series of image), we notice a high degree of mixing, being the color of the images rather uniform, however exhibiting some signatures of the unsteady flow. Such signatures attenuate while moving further downstream (bottom series).
Similar experimental snapshots are shown in
Figure 4 for
. The related video is provided in the
Supplementary Material. Here we observe how at the confluence the flow preserves a large degree of symmetry, resembling the vortex regime shown in
Figure 2b, except for the dynamic feature. Hence, this symmetric unsteady regime, which has been already observed when feeding the same fluid, i.e., water, to the microreactor, occurs also in the present case involving fluids with different density and viscosity. We notice from the
concentration maps just small oscillations of the contact region. Mixing is clearly hampered with respect to the previous regime, despite the fact that the
is almost doubled. Such segregation persists further down in the mixing channel, the bottom series of images still showing dark and white signatures of the inlet streams on the left and right, respectively.
The corresponding degree of mixing at
is reported in
Figure 5 for the two Reynolds numbers at different instants, which are made non-dimensional with respect to the cycle period, i.e.,
with
being the Strouhal number. At
we estimated
from the numerical simulations, while a slightly lower
was determined for
. The experimental values of the Strouhal number were
at
and
at
, hence showing a very satisfactory agreement with numerical ones.
The degree of mixing,
is defined as in [
26,
32,
38], i.e.,
where
is the standard deviation of the volumetric dye flow and
is the maximum value of
obtained for completely segregated streams. Thus,
varies between
, indicating a completely segregated flow, and
, corresponding to a perfectly mixed flow. In our case the dye is HCl, as this component does not participate in the chemical reaction and hence acts as a passive tracer.
For both
, the mixing degree fluctuates being affected by the flow dynamics. In the unsteady asymmetric regime, i.e., at
we notice how the signal exhibits a double peak (see
Figure 5a); one peak can be related to the passage of the vorticity blob, while the other to the instant when the vortical legs in the channel are strongest. In this regime, the averaged
at
is
. In the unsteady symmetric regime, i.e., at
, the dynamics is simpler with a nearly sinusoidal shape of the
signal (see
Figure 5b). The degree of mixing at
is
as reported also in
Table 1.
The above
values are compared to those of the steady regimes in
Figure 6. We notice how the unsteady motions boost mixing as
increases from about 40% at
(i.e., steady asymmetric regime) to 65.9% at
(i.e., unsteady asymmetric regime); however a further increase of
towards the unsteady symmetric regime is detrimental.
Figure 7 compares experimental and numerical temporal signals of the average
concentration, which is made non-dimensional with respect to the inlet
concentration, at the
,
and
cross-sections (from top to bottom) and for
(
Figure 7a) and
(
Figure 7b).
We can notice how the decreases moving along the channel as the reaction proceeds and how all signals show the fingerprint of the dynamic behavior which affects the entire mixing channel. The good agreement between experiments and numerical simulations confirms that the present numerical setup is able to well capture the evolution of the chemical reaction.
The corresponding experimental and numerical reaction yields are reported in
Figure 8a–d, respectively, along with data from the steady regimes, already discussed in [
32]. The reaction yield is evaluated at the
from the
concentration as:
with
being the
concentration averaged on the cross-section. The lines refer to the scaling of the reaction yield with the Reynolds and Damköhler numbers suggested in the aforementioned paper for the different regimes, respectively, i.e.,:
in the segregated regime (see solid line)
in the vortex and engulfment regimes (dashed line).
Here the Damköhler number, representing the flow to chemical time-scale ratio, is
and
is the non-dimensional kinetic constant, i.e.,
The agreement between the experimental and numerical yield is very satisfactory. We notice how the yield at is enhanced considerably with respect to the value we would expect in case of a steady engulfment regime. In particular at while increases up to at . More precisely, in the simulations and in experiments, thus further confirming the accuracy of the numerical setup. Such an increase of reaction yield is motivated by the periodic motions which strongly boost the mixing of reactants.
Instead, the yield drops significantly at
down to
(
in simulations and
in experiments), thus at values that are even lower than those of the segregated regime, i.e.,
for
. Indeed, despite the fact that the oscillations of the contact region in the symmetric unsteady regime improve mixing with respect to the segregated regime, the yield is penalized by the much lower residence time at
. It is worth noting that
Figure 8c,d indicate how this yield is somewhat larger than the one we would expect from the
relationship, i.e., based solely on the residence time, because of the periodic motions.
5. Conclusions
In the present work experiments and numerical simulations were used jointly to investigate the presence of unsteady flows in case of reactive system and analyze their impact on the reaction progress of a chemical reaction. Two different unsteady regimes, i.e., asymmetric and symmetric, were observed with the same features already described for non-reactive flows. The numerical simulations provided results which were closely in agreement with the experiments. We highlight that a grid adaption cannot be effectively implemented in transient simulations to improve the description of the mixing of reactants; however the unsteady regimes further promote convection and increase the mixing length, making the computational grid fine enough to capture mixing.
The unsteady motions strongly boost the progress of the chemical reaction, its yield increasing by more than 30% when moving from the steady flow at to the unsteady flow at . Such reaction enhancement is higher than the one we would estimate by considering the relationship developed for the steady regime. According to this relationship, if at in the engulfment regime we would expect at , while the unsteady motions push the yield up to .
However, the onset of the unsteady symmetric regime at higher flow rates, has a detrimental effect on the reaction progress, with , thus lower than the yields of the steady segregated regime. In fact, despite oscillations of the contact region enhance somewhat the mixing of reactants, the very low residence time hampers the reaction progress.
To our knowledge the present work provides the first evidence that unsteady flow regimes may be effectively exploited to improve significantly the yield of chemical reactions. Hence we believe that this information is useful for enabling process intensification in microreactors.