Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Reconfigurable Metasurface Antenna Based on the Liquid Metal for Flexible Scattering Fields Manipulation
Next Article in Special Issue
An Investigation of Flow Patterns and Mixing Characteristics in a Cross-Shaped Micromixer within the Laminar Regime
Previous Article in Journal
All-Dielectric Color Filter with Ultra-Narrowed Linewidth
Previous Article in Special Issue
CFD Simulations of Microreactors for the Hydrolysis of Cellobiose to Glucose by β-Glucosidase Enzyme
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study on the Effect of Flow Unsteadiness on the Yield of a Chemical Reaction in a T Micro-Reactor

by
Alessandro Mariotti
,
Matteo Antognoli
,
Chiara Galletti
*,†,
Roberto Mauri
,
Maria Vittoria Salvetti
and
Elisabetta Brunazzi
Dipartimento di Ingegneria Civile e Industriale, Università di Pisa, Largo Lazzarino 2, 56122 Pisa, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Micromachines 2021, 12(3), 242; https://doi.org/10.3390/mi12030242
Submission received: 3 February 2021 / Revised: 22 February 2021 / Accepted: 24 February 2021 / Published: 27 February 2021
(This article belongs to the Special Issue Micro-Reaction Engineering)

Abstract

:
Despite the very simple geometry and the laminar flow, T-shaped microreactors have been found to be characterized by different and complex steady and unsteady flow regimes, depending on the Reynolds number. In particular, flow unsteadiness modifies strongly the mixing process; however, little is known on how this change may affect the yield of a chemical reaction. In the present work, experiments and 3-dimensional numerical simulations are carried out jointly to analyze mixing and reaction in a T-shaped microreactor with the ultimate goal to investigate how flow unsteadiness affects the reaction yield. The onset of the unsteady asymmetric regime enhances the reaction yield by more than 30 % ; however, a strong decrease of the yield back to values typical of the vortex regime is observed when the flow undergoes a transition to the unsteady symmetric regime.

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., R 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 R e , 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 R e 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., η = f ( D a ) , where D a is the Damköhler number, while in the vortex and engulfment regimes, an explicit dependency on the kinetic constant k r was observed with η = f ( k r , D a ) .
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 W i = H = 1 mm, while the mixing channel has a rectangular cross-section with aspect ratio 2:1, i.e., W o = 2 H = 2 mm. Hence, the hydraulic diameter of the mixing channel is d = 4 H / 3 . In the following text, coordinates will be made non-dimensional by using the hydraulic diameter, i.e., X = x / d , Y = y / d , and Z = z / d . The inlet channels length is equal L i = 30 d , so the flow is fully developed at the T-junction, whereas the mixing channel length is equal to L o = 45 d 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).
MB + + AsA LMB + + DA
The aqueous solution of MB + 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 MB + powder (by Sigma Aldrich) are dissolved in aqueous solutions of 1 mol/L of HCl (by Sigma Aldrich), thus obtaining a concentration of MB + equal to [ MB + ] = 5.31 × 10 5 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] = 1.7 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] ≫ [ MB + ] as in the present case, the above reaction has a pseudo-first order kinetic [34], that depends on the concentration of MB + only with a kinetic constant, k r = 10.71 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., δ T = 10 6 s. Each image collects a region of 1.70 d × 3.3 d (the equivalent dimension of a pixel is 2.93 μm). Among the species involved in the reaction, MB + 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 MB + 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 MB + , LMB + , AsA, HCl and DA. The energy equation is not resolved as thermal effects can be neglected. The non-dimensional form of the equation is:
ρ ^ θ + · ( ρ ^ u ) = 0 ,
ρ ^ ( u θ + u · u ) = p + 1 R e · [ μ ^ ( u + u T ) ] + R i ( ρ ^ 1 ) g ^ ,
ρ ^ ( ϕ k θ + u · ϕ k ) = 1 P e · ( ρ ^ D k ^ ϕ k ) + d ω ˙ k ρ 0 U ,
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., θ = t U d . In the above equations, u represents the velocity vector, p is the modified non-dimensional pressure, i.e.,
p = ( P ρ 0 g Z ) / ρ 0 U 2
where P is the pressure, and g the gravity acceleration, while g ^ is the non-dimensional gravity, i.e., g ^ = g g . ρ ^ = ρ / ρ 0 and μ ^ = μ / μ 0 are non-dimensional density and viscosity, respectively, referred to a convenient reference state, i.e., μ 0 and ρ 0 , which in our case corresponds to pure water at 20 °C. ϕ k represents the mass fraction of the k-th chemical species while ω ˙ k is its rate of production or consumption due to chemical reactions, which, for instance is:
ω ˙ MB + = k r ρ ϕ MB +
for the methylene blue. D k ^ = D / D 0 is the non-dimensional diffusivity referred to the water self-diffusivity D 0 .
The characteristic non-dimensional numbers are:
  • the Reynolds number
    R e = ρ 0 U d μ 0
  • the Richardson number
    R i = g d Δ ρ ρ 0 U 2
  • the Peclet number
    P e = U d D 0
where ν 0 = μ 0 / ρ 0 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 ϕ A s A and hence, the correlations
μ = μ o f ϕ A s A
ρ = ρ o f ϕ A s A
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 R e < 60 for the given size of our microreactor [32]. Despite the fact that stratification is therefore negligible in the range of R e 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 5 , 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 MB + concentration, estimated from the numerical simulations, in the confluence region. Such concentration is reported in non-dimensional form, being 1 the MB + 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 R e = 230 ) breaks the flow symmetries, thus strongly enhancing mixing with respect to the segregated (Figure 2a at R e = 60 ) and vortex regimes (Figure 2b at R e = 100 ), 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 R e = 320 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 2.55 < Y < 0.75 , then in the region at 9.15 < Y < 5.85 and further along the mixing channel 26.65 < Y < 23.35 . 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 MB + 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 R e = 650 . 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 MB + concentration maps just small oscillations of the contact region. Mixing is clearly hampered with respect to the previous regime, despite the fact that the R e 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 Y = 25 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., τ = d U · S t with S t being the Strouhal number. At R e = 320 we estimated S t = 0.253 from the numerical simulations, while a slightly lower S t = 0.228 was determined for R e = 650 . The experimental values of the Strouhal number were S t = 0.246 at R e = 320 and S t = 0.224 at R e = 650 , hence showing a very satisfactory agreement with numerical ones.
The degree of mixing, δ m is defined as in [26,32,38], i.e.,
δ m = 1 σ b σ m a x ,
where σ b is the standard deviation of the volumetric dye flow and σ m a x is the maximum value of σ b obtained for completely segregated streams. Thus, δ m varies between δ m = 0 , indicating a completely segregated flow, and δ m = 1 , 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 R e , the mixing degree fluctuates being affected by the flow dynamics. In the unsteady asymmetric regime, i.e., at R e = 320 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 δ m at Y = 25 is δ m = 65.9 ± 1.6 % . In the unsteady symmetric regime, i.e., at R e = 650 , the dynamics is simpler with a nearly sinusoidal shape of the δ m signal (see Figure 5b). The degree of mixing at Y = 25 is δ m = 50.1 ± 2.1 % as reported also in Table 1.
The above δ m values are compared to those of the steady regimes in Figure 6. We notice how the unsteady motions boost mixing as δ m increases from about 40% at R e = 280 (i.e., steady asymmetric regime) to 65.9% at R e = 320 (i.e., unsteady asymmetric regime); however a further increase of R e towards the unsteady symmetric regime is detrimental.
Figure 7 compares experimental and numerical temporal signals of the average MB + concentration, which is made non-dimensional with respect to the inlet MB + concentration, at the Y = 2 , Y = 7.5 and Y = 25 cross-sections (from top to bottom) and for R e = 320 (Figure 7a) and R e = 650 (Figure 7b).
We can notice how the MB + 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 Y = 25 from the MB + concentration as:
η = 1 2 C B M + , Y = 25 C B M 0 +
with C B M + , Y 25 being the MB + 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)
    η D a 0.3
  • in the vortex and engulfment regimes (dashed line).
    η k r ˜ 0.1 D a 0.4
Here the Damköhler number, representing the flow to chemical time-scale ratio, is
D a = d k r / U = k r ˜ / R e
and k r ˜ is the non-dimensional kinetic constant, i.e.,
k r ˜ = k r d 2 ν
The agreement between the experimental and numerical yield is very satisfactory. We notice how the yield at R e = 320 is enhanced considerably with respect to the value we would expect in case of a steady engulfment regime. In particular η 28 % at R e = 280 while increases up to η 35 % at R e = 320 . More precisely, η = 34.9 ± 1.6 % in the simulations and η = 35.4 ± 2.5 % 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 R e = 650 down to η 11 % ( η = 11.4 ± 2.7 % in simulations and η = 10.7 ± 4.4 % in experiments), thus at values that are even lower than those of the segregated regime, i.e., η 15 % for R e = 60 . 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 R e = 650 . It is worth noting that Figure 8c,d indicate how this yield is somewhat larger than the one we would expect from the η D a 0.3 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 R e = 280 to the unsteady flow at R e = 320 . Such reaction enhancement is higher than the one we would estimate by considering the relationship η k r ˜ 0.1 D a 0.4 developed for the steady regime. According to this relationship, if η 28 % at R e = 280 in the engulfment regime we would expect η 30 % at R e = 320 , while the unsteady motions push the yield up to η = 35 % .
However, the onset of the unsteady symmetric regime at higher flow rates, has a detrimental effect on the reaction progress, with η 11 % , 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.

Supplementary Materials

The following videos are available online at https://www.mdpi.com/2072-666X/12/3/242/s1, Video S1: Experimental depth-averaged MB + concentration at R e = 320 .; Video S2: Experimental depth-averaged MB + concentration at R e = 650 .

Author Contributions

Conceptualization and methodology, A.M., C.G., E.B., R.M. and M.V.S.; experimental investigation, A.M., M.A. and E.B.; numerical modeling, A.M., M.A. and C.G.; validation, R.M. and M.V.S.; writing—original draft preparation, C.G., A.M. and E.B.; writing—review and editing, R.M. and M.V.S. All authors have read and agreed to the published version of the manuscript.

Funding

Support by the University of Pisa through the “Progetti di Ricerca di Ateneo PRA 2017–2018” funding program.

Acknowledgments

The authors are grateful to Simone Egidi, Giacomo Flori, Sara Tomasi Masoni and Cesare Merello.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following symbols and abbreviations are used in this manuscript:
Δ ρ density difference between the two inlet fluids, kg m 3
δ m degree of mixing, -
δ T image exposure time, s
η reaction yield, -
θ non-dimensional time, -
μ dynamic viscosity, kg m 1 s 1
μ ^ non-dimensional kinematic viscosity, -
μ 0 dynamic viscosity of pure water at 20 °C, kg m 1 s 1
ν 0 kinematic viscosity of pure water at 20 °C, m 2 s 1
ρ density, kg m 3
ρ ^ non-dimensional density, -
ρ 0 density of pure water at 20 °C, kg m 3
σ b standard deviation of the volumetric dye flow, -
σ m a x maximum value of the standard deviation of the volumetric dye flow, -
τ cycle period, s
ϕ k mass fraction of the k-th chemical species, -
ω ˙ k rate of production or consumption of ϕ k due to chemical reactions, kg m 3 s 1
AsA ascorbic acid
C B M + methylene blue concentration, mol m 3
CFLCourant–Friedrichs–Lewy number, -
dmixing channel hydraulic diameter, m
D diffusivity, m 2 s 1
D k ^ non-dimensional diffusivity, -
D 0 water self-diffusivity, m 2 s 1
D a Damköhler number, -
DA dehydroascorbic acid
ggravity acceleration, m s 2
g ^ non-dimensional gravity, -
Hchannel height, m
HCl hydrogen chloride
k r kinetic constant, s 1
k r ˜ non-dimensional kinetic constant, -
L i inlet channel length, m
LMB + leuco compound
L o mixing channel length, m
MB + methylene blue
N.A.numeric aperture, -
pmodified non-dimensional pressure, -
Ppressure, Pa
P e Peclet number, -
PISOPressure implicit with splitting of operator algorithm
R e Reynolds number, -
R i Richardson number, -
S t Strouhal number, -
ttime, s
Ubulk velocity, m s 1
u velocity vector, m s 1
W i inlet channels width, m
W o mixing channel width, m
xx-coordinate, m
Xnon-dimensional x-coordinate, -
yy-coordinate, m
Ynon-dimensional y-coordinate, -
zz-coordinate, m
Znon-dimensional z-coordinate, -

References

  1. Rossetti, I.; Compagnoni, M. Chemical reaction engineering, process design and scale-up issues at the frontier of synthesis: Flow chemistry. Chem. Eng. J. 2016, 296, 56–70. [Google Scholar] [CrossRef]
  2. Roberge, D.M.; Ducry, L.; Bieler, N.; Cretton, P.; Zimmermann, B. Microreactor technology: A revolution for the fine chemical and pharmaceutical industries? Chem. Eng. Technol. 2005, 28, 318–323. [Google Scholar] [CrossRef]
  3. Pommella, A.; Tomaiuolo, G.; Chartoire, A.; Caserta, S.; Toscano, G.; Nolan, S.P.; Guido, S. Palladium-N-heterocyclic carbene (NHC) catalyzed C–N bond formation in a continuous flow microreactor. Effect of process parameters and comparison with batch operation. Chem. Eng. J. 2013, 223, 578–583. [Google Scholar] [CrossRef]
  4. Rossetti, I. Continuous flow (micro-)reactors for heterogeneously catalyzed reactions: Main design and modelling issues. Catal. Today 2018, 308, 20–31. [Google Scholar] [CrossRef]
  5. Russo, D.; Tomaiuolo, G.; Andreozzi, R.; Guido, S.; Lapkin, A.A.; Di Somma, I. Heterogeneous benzaldehyde nitration in batch and continuous flow microreactor. Chem. Eng. J. 2019, 377, 120346. [Google Scholar] [CrossRef]
  6. Gani, R.; Bałdyga, J.; Biscans, B.; Brunazzi, E.; Charpentier, J.C.; Drioli, E.; Feise, H.; Furlong, A.; Van Geem, K.M.; de Hemptinne, J.C.; et al. A multi-layered view of chemical and biochemical engineering. Chem. Eng. Res. Des. 2020, 155, A133–A145. [Google Scholar] [CrossRef]
  7. Lee, S.H.; Kang, P.K. Three-Dimensional Vortex-Induced Reaction Hot Spots at Flow Intersections. Phys. Rev. Lett. 2020, 124, 144501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Wong, S.H.; Ward, M.C.; Wharton, C.W. Micro T-mixer as a rapid mixing micromixer. Sens. Actuators B Chem. 2004, 100, 359–379. [Google Scholar] [CrossRef]
  9. Dreher, S.; Kockmann, N.; Woias, P. Characterization of laminar transient flow regimes and mixing in T-shaped micromixers. Heat Transf. Eng. 2009, 30, 91–100. [Google Scholar] [CrossRef]
  10. Galletti, C.; Brunazzi, E.; Mauri, R. Unsteady mixing of binary liquid mixtures with composition-dependent viscosity. Chem. Eng. Sci. 2017, 164, 333–343. [Google Scholar] [CrossRef]
  11. Minakov, A.; Rudyak, V.; Gavrilov, A.; Dekterev, A. Mixing in a T-shaped micromixer at moderate Reynolds numbers. Thermophys. Aeromech. 2012, 19, 385–395. [Google Scholar] [CrossRef]
  12. Kockmann, N.; Kiefer, T.; Engler, M.; Woias, P. Convective mixing and chemical reactions in microchannels with high flow rates. Sens. Actuators B Chem. 2006, 117, 495–508. [Google Scholar] [CrossRef]
  13. Fani, A.; Camarri, S.; Salvetti, M.V. Investigation of the steady engulfment regime in a three-dimensional T-mixer. Phys. Fluids 2013, 25, 064102. [Google Scholar] [CrossRef] [Green Version]
  14. Poole, R.J.; Alfateh, M.; Gauntlett, A.P. Bifurcation in a T-channel junction: Effects of aspect ratio and shear-thinning. Chem. Eng. Sci. 2013, 104, 839–848. [Google Scholar] [CrossRef]
  15. Lobasov, A.S.; Minakov, A.V.; Kuznetsov, V.V.; Rudyak, V.Y.; Shebeleva, A.A. Investigation of mixing efficiency and pressure drop in T-shaped micromixers. Chem. Eng. Process. Process Intensif. 2018, 134, 105–114. [Google Scholar] [CrossRef] [Green Version]
  16. Mariotti, A.; Lanzetta, M.; Dini, G.; Rossi, A.; Brunazzi, E.; Mauri, R.; Galletti, C. Influence of cross-sectional geometry on mixing in a T-shaped micro-junction. Chem. Eng. Trans. 2019, 74, 955–960. [Google Scholar]
  17. Mariotti, A.; Galletti, C.; Brunazzi, E.; Salvetti, M.V. Steady flow regimes and mixing performance in arrow-shaped micro-mixers. Phys. Rev. Fluids 2019, 4, 034201. [Google Scholar] [CrossRef]
  18. Mariotti, A.; Galletti, C.; Brunazzi, E.; Salvetti, M.V. Unsteady flow regimes in arrow-shaped micro-mixers with different tilting angles. Phys. Fluids 2021, 33, 012008. [Google Scholar] [CrossRef]
  19. You, B.J.; Choi, Y.; Im, S.G. Influence of adjusting the inlet channel confluence angle on mixing behaviour in inertial microfluidic mixers. Microfluid. Nanofluid. 2017, 21, 121. [Google Scholar] [CrossRef]
  20. Babu, H.; Satu, S.; Haderlein, M.; Peukert, W.; Verma, N. Numerical investigation of flow patterns and concentration profiles in Y-mixers. Chem. Eng. Technol. 2016, 39, 1963. [Google Scholar] [CrossRef]
  21. Ansari, M.A.; Kim, K.Y.; Kim, S.M. Numerical and Experimental Study on Mixing Performances of Simple and Vortex Micro T-Mixers. Micromachines 2018, 9, 204. [Google Scholar] [CrossRef] [Green Version]
  22. Chan, S.T.; Ault, J.T.; Haward, S.J.; Meiburg, E.; Shen, A.Q. Coupling of vortex breakdown and stability in a swirling flow. Phys. Rev. Fluids 2019, 4, 084701. [Google Scholar] [CrossRef] [Green Version]
  23. Liu, Z.; Lu, Y.; Wang, J.; Luo, G. Mixing characterization and scaling-up analysis of asymmetrical T-shaped micromixer: Experiment and CFD simulation. Chem. Eng. J. 2012, 181–182, 597–606. [Google Scholar] [CrossRef]
  24. Wu, C.; Lai, B. Numerical study of T-shaped micromixers with vortex-inducing obstacles in the inlet channels. Micromachines 2020, 11, 1122. [Google Scholar] [CrossRef]
  25. Lobasov, A.S.; Minakov, A.V. Analyzing mixing quality in a T-shaped micromixer for different fluids properties through numerical simulation. Chem. Eng. Process. Process Intensif. 2018, 124, 11–23. [Google Scholar] [CrossRef] [Green Version]
  26. Galletti, C.; Arcolini, G.; Brunazzi, E.; Mauri, R. Mixing of binary fluids with composition-dependent viscosity in a T-shaped micro-device. Chem. Eng. Sci. 2015, 123, 300–310. [Google Scholar] [CrossRef] [Green Version]
  27. Orsi, G.; Roudgar, M.; Brunazzi, E.; Galletti, C.; Mauri, R. Water-ethanol mixing in T-shaped microdevices. Chem. Eng. Sci. 2013, 95, 174–183. [Google Scholar] [CrossRef]
  28. Galletti, C.; Roudgar, M.; Brunazzi, E.; Mauri, R. Effect of inlet conditions on the engulfment pattern in a T-shaped micro-mixer. Chem. Eng. J. 2012, 185-186, 300–313. [Google Scholar] [CrossRef]
  29. Dundi, T.M.; Raju, V.R.K.; Chandramohan, V.P. Numerical evaluation of swirl effect on liquid mixing in a passive T-micromixer. Aust. J. Mech. Eng. 2019, 0, 1–15. [Google Scholar] [CrossRef]
  30. Soleymani, A.; Yousefi, H.; Turunen, I. Dimensionless number for identification of flow patterns inside a T-micromixer. Chem. Eng. Sci. 2008, 63, 5291–5297. [Google Scholar] [CrossRef]
  31. Camarri, S.; Mariotti, A.; Galletti, C.; Brunazzi, E.; Mauri, R.; Salvetti, M.V. An Overview of Flow Features and Mixing in Micro T and Arrow Mixers. Ind. Eng. Chem. Res. 2020, 59, 3669–3686. [Google Scholar] [CrossRef]
  32. Mariotti, A.; Antognoli, M.; Galletti, C.; Mauri, R.; Salvetti, M.V.; Brunazzi, E. The role of flow features and chemical kinetics on the reaction yield in a T-shaped micro-reactor. Chem. Eng. J. 2020, 396, 125223. [Google Scholar] [CrossRef]
  33. Mariotti, A.; Galletti, C.; Mauri, R.; Salvetti, M.V.; Brunazzi, E. Effect of stratification on the mixing and reaction yield in a T-shaped micro-mixer. Phys. Rev. Fluids 2021, 004200, 1–29. [Google Scholar]
  34. Mowry, S.; Ogren, P.J. Kinetics of Methylene Blue Reduction by Ascorbic Acid. J. Chem. Educ. 1999, 76, 970–973. [Google Scholar] [CrossRef]
  35. Mariotti, A.; Galletti, C.; Mauri, R.; Salvetti, M.V.; Brunazzi, E. Steady and unsteady regimes in a T-shaped micro-mixer: Synergic experimental and numerical investigation. Chem. Eng. J. 2018, 341, 414–431. [Google Scholar] [CrossRef]
  36. Nishikata, E.; Ishii, T.; Ohta, T. Viscosities of aqueous hydrochloric acid solutions, and densities and viscosities of aqueous hydroiodic acid solutions. J. Chem. Eng. Data 1981, 26, 254–256. [Google Scholar] [CrossRef]
  37. ANSYS Inc. ANSYS Fluent Theory Guide, 2019R2; ANSYS Inc.: Canonsburg, PA, USA, 2019. [Google Scholar]
  38. Galletti, C.; Mariotti, A.; Siconolfi, L.; Mauri, R.; Brunazzi, E. Numerical investigation of flow regimes in T-shaped micromixers: Benchmark between finite volume and spectral element methods. Can. J. Chem. Eng. 2019, 97, 528–541. [Google Scholar] [CrossRef]
  39. Mariotti, A.; Galletti, C.; Salvetti, M.V.; Brunazzi, E. Unsteady Flow Regimes in a T-Shaped Micromixer: Mixing and Characteristic Frequencies. Ind. Eng. Chem. Res. 2019, 58, 13340–13356. [Google Scholar] [CrossRef]
Figure 1. (a) Sketch of the T-reactor geometry and reference system; (b) experimental set-up.
Figure 1. (a) Sketch of the T-reactor geometry and reference system; (b) experimental set-up.
Micromachines 12 00242 g001
Figure 2. Depth-averaged MB + concentration from numerical simulations in: (a) the segregated regime, (b) the vortex regime, and (c) the engulfment regime. Non-dimensional velocity magnitude and flow streamlines in the horizontal mid-plane in: (d) the segregated regime, (e) the vortex regime, and (f) the engulfment regime.
Figure 2. Depth-averaged MB + concentration from numerical simulations in: (a) the segregated regime, (b) the vortex regime, and (c) the engulfment regime. Non-dimensional velocity magnitude and flow streamlines in the horizontal mid-plane in: (d) the segregated regime, (e) the vortex regime, and (f) the engulfment regime.
Micromachines 12 00242 g002
Figure 3. Experimental depth-averaged MB + concentration at R e = 320 , for instants: (a) t / τ = 0 , (b) t / τ = 0.17 , (c) t / τ = 0.33 , (d) t / τ = 0.5 , (e) t / τ = 0.67 , (f) t / τ = 0.83 .
Figure 3. Experimental depth-averaged MB + concentration at R e = 320 , for instants: (a) t / τ = 0 , (b) t / τ = 0.17 , (c) t / τ = 0.33 , (d) t / τ = 0.5 , (e) t / τ = 0.67 , (f) t / τ = 0.83 .
Micromachines 12 00242 g003
Figure 4. Experimental depth-averaged MB + concentration at R e = 650 for instants: (a) t / τ = 0 , (b) t / τ = 0.25 , (c) t / τ = 0.50 , (d) t / τ = 0.75 .
Figure 4. Experimental depth-averaged MB + concentration at R e = 650 for instants: (a) t / τ = 0 , (b) t / τ = 0.25 , (c) t / τ = 0.50 , (d) t / τ = 0.75 .
Micromachines 12 00242 g004
Figure 5. Numerical time signal of the mixing degree at Y = 25 for (a) R e = 320 and (b) R e = 650 .
Figure 5. Numerical time signal of the mixing degree at Y = 25 for (a) R e = 320 and (b) R e = 650 .
Micromachines 12 00242 g005
Figure 6. Mixing degree at Y = 25 as a function of Reynolds number for steady and unsteady regimes.
Figure 6. Mixing degree at Y = 25 as a function of Reynolds number for steady and unsteady regimes.
Micromachines 12 00242 g006
Figure 7. Experimental and numerical time signals of the MB + concentration at for (a) R e = 320 and (b) R e = 650 at cross- sections (from top to bottom): Y = 2 , Y = 7.5 , and Y = 25 .
Figure 7. Experimental and numerical time signals of the MB + concentration at for (a) R e = 320 and (b) R e = 650 at cross- sections (from top to bottom): Y = 2 , Y = 7.5 , and Y = 25 .
Micromachines 12 00242 g007
Figure 8. Reaction yield at Y = 25 in experiments (a,c) and numerical simulations (b,d) as a function of Reynolds number (a,b) and Damköhler number (c,d).
Figure 8. Reaction yield at Y = 25 in experiments (a,c) and numerical simulations (b,d) as a function of Reynolds number (a,b) and Damköhler number (c,d).
Micromachines 12 00242 g008
Table 1. Numerical values of the mixing degree and reaction yield at Y = 25 .
Table 1. Numerical values of the mixing degree and reaction yield at Y = 25 .
Re δ m % η %
320 65.9 ± 1.6 34.9 ± 1.6
650 50.1 ± 2.1 11.4 ± 2.7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mariotti, A.; Antognoli, M.; Galletti, C.; Mauri, R.; Salvetti, M.V.; Brunazzi, E. A Study on the Effect of Flow Unsteadiness on the Yield of a Chemical Reaction in a T Micro-Reactor. Micromachines 2021, 12, 242. https://doi.org/10.3390/mi12030242

AMA Style

Mariotti A, Antognoli M, Galletti C, Mauri R, Salvetti MV, Brunazzi E. A Study on the Effect of Flow Unsteadiness on the Yield of a Chemical Reaction in a T Micro-Reactor. Micromachines. 2021; 12(3):242. https://doi.org/10.3390/mi12030242

Chicago/Turabian Style

Mariotti, Alessandro, Matteo Antognoli, Chiara Galletti, Roberto Mauri, Maria Vittoria Salvetti, and Elisabetta Brunazzi. 2021. "A Study on the Effect of Flow Unsteadiness on the Yield of a Chemical Reaction in a T Micro-Reactor" Micromachines 12, no. 3: 242. https://doi.org/10.3390/mi12030242

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop