Numerical Study of the Flow and Mass Transfer in Micromixers

Author(s):  
Nassim Ait Mouheb ◽  
Camille Solliec ◽  
Agnes Montillet ◽  
Jacques Comiti ◽  
Patrick Legentilhomme ◽  
...  

Computational fluid dynamic simulations are used to characterize the flow and the liquid mixing quality in a micromixer as a function of the Reynolds number. Two micromixers are studied in steady flow conditions; they are based on two geometries, respectively T-shaped and cross-type (⊤ and + shapes). Simulations allow, in the case of ⊤ micromixers, to chart the topology of the flow and to describe the evolution of species concentration downstream the intersection. The streamline layout and the mixing quality curves reveal the three characteristic types of flow, depending on Reynolds number: stratified, vortex and engulfment flows. Vortices appear after impingement, in the exit channel. They become asymmetrical and gain in length with an increase in Re making the flow unsteady, which induces an enhancement of the mass transfer by advection between the two liquids. In the case of cross-type micromixers, the structure of the flow is strongly three-dimensional. It is characterized by symmetrical vortices in both output channels. In the zone close to the impingement, a back flow is observed which induces strong shear stresses. The results show that the + shaped system can improve the mixing process in comparison with the micromixers having ⊤ geometry. The numerical study also allows to select the locations of the most relevant zones of study, from an experimental point of view. It will allow to choose the location of PIV planes and local non intrusive sensors, such as electrochemical microprobes, in order to experimentally investigate the flow.

Author(s):  
Nassim Ait Mouheb ◽  
Agnes Montillet ◽  
Camille Solliec ◽  
Jacques Comiti ◽  
Patrick Legentilhomme ◽  
...  

The understanding of physical phenomena such as flow behaviour and mass transfer performance is needed in order to develop appropriate micromixers for industrial or biomedical applications. In this work, CFD is used to characterize the flow and the liquid mixing quality in a micromixer as a function of the Reynolds number. Two micromixers are studied in steady flow conditions; they are based on two geometries, respectively T-shaped (⊤) and cross-type (+). Simulations allow, in the case of ⊤ micromixers, to chart the topology of the flow and to describe the evolution of species concentration downstream the crossing. The streamlines layout and the mixing quality curves reveal three characteristic types of flow previously reported in the literature, depending on Reynolds number: stratified, vortex and engulfment flows. In the case of cross-type micromixers, the structure of the flow is strongly three-dimensional and is characterized by symmetrical vortices in both output channels. The results show that the + shaped system can improve the mixing process in comparison with the micromixers having ⊤ geometry. The second part of the study is experimental. Two cells are constructed, for both geometries (T-shaped and cross) using square channels with 400 μm hydraulic diameter. In order to use particle image velocimetry (PIV), a system has been adapted to measure velocity fields for various channel plans at different channel depths. This allows observing the evolution of the flow and the vortices development along the microchannels. A second experimental technique, the electrochemical one involving microelectrodes implemented at several positions on the channel wall located near the crossing, has been used. The electrochemical method can locally characterize the formation of swirling flows. These two complementary experimental results will be analysed and a comparison with the CFD results will be performed.


Author(s):  
M. Pau ◽  
G. Paniagua

Ensuring an adequate life of high pressure turbines requires efficient cooling methods, such as rim seal flow ejection from the stator-rotor wheel space cavity interface, which prevents hot gas ingress into the rotor disk. The present work addresses the potential to improve the efficiency in transonic turbines at certain rim seal ejection rates. To understand this process a numerical study was carried out combining computational fluid dynamic simulations (CFD) and experiments on a single stage axial test turbine. The three dimensional steady CFD analysis was performed modeling the purge cavity flow ejected downstream of the stator blade row, at three flow regimes, subsonic M2 = 0.73, transonic M2 = 1.12 and supersonic M2 = 1.33. Experimental static pressure measurements were used to calibrate the computational model. The main flow field-purge flow interaction is found to be governed by the vane shock structures at the stator hub. The interaction between the vane shocks at the hub and the purge flow has been studied and quantitatively characterized as function of the purge ejection rate. The ejection of 1% of the core flow from the rim seal cavity leads to an increase of the hub static pressure of approximately 7% at the vane trailing edge. This local reduction of the stator exit Mach number decreases the trailing edge losses in the transonic regime. Finally, a numerically predicted loss breakdown is presented, focusing on the relative importance of the trailing edge losses, boundary layer losses, shock losses and mixing losses, as a function of the purge rate ejected. Contrary to the experience in subsonic turbines, results in a transonic model demonstrate that ejecting purge flow improves the vane efficiency due to the shock structures modification downstream of the stator.


2020 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Rawad Himo ◽  
Charbel Bou-Mosleh ◽  
Charbel Habchi

Purpose Flow separation on wings, blades and vehicles can be delayed or even suppressed by the use of vortex generators (VG). Numerous studies, documented in the literature, extensively describe the performance of triangular and rectangular VG winglets. This paper aims to focus on the use of non-conventional VG shapes, more specifically an array of trapezoildal-perforated VG tabs. Design/methodology/approach In this study, computational fluid dynamic simulations are performed on an inline array of trapezoidal VG with various dimensions and inclination angles, in addition to considering perforations in the VG centers. The methodology of the present numerical study is validated with experimental data from the literature. Findings The performance and the associated flow structures of these tested non-conventional VG are compared to classical triangular winglets. For the proposed non-conventional trapezoidal VG, at the onset of stall, a 21% increase of lift over drag on the airfoil is observed. The trapezoidal VG enhancement is also witnessed during stall where the lift over drag ratio is increased by 120% for the airfoil and by 10% with respect to the triangular winglets documented in the literature. Originality/value The originality of this paper is the use of non-conventional vortex generator shape to enhance lift over drag coefficient using three-dimensional numerical simulations.


2014 ◽  
Vol 660 ◽  
pp. 674-678 ◽  
Author(s):  
Azmahani Sadikin ◽  
Nurul Akma Mohd Yunus ◽  
Kamil Abdullah ◽  
Akmal Nizam Mohammed

The unsteady three dimensional flow simulation around sphere using numerical simulation computational fluid dynamic for moderate Reynolds Number between 20 ≤ Re ≤ 500 is presented. The aim of this work is to analyze the flow regimes around sphere and flow separation. Extensive comparisons were made between the present predicted results and available experimental and numerical investigations, and showed that they are in close agreement. The results show that the vortex shedding increases with the Reynolds number. The flow separates early when Reynolds number increases, therefore the separation angle is found to be smaller when high Reynolds number is present.


2009 ◽  
Vol 13 (3) ◽  
pp. 59-67 ◽  
Author(s):  
Enrico Mollica ◽  
Eugenio Giacomazzi ◽  
Marco di

In this article a combustor burning hydrogen and air in mild regime is numerically studied by means of computational fluid dynamic simulations. All the numerical results show a good agreement with experimental data. It is seen that the flow configuration is characterized by strong exhaust gas recirculation with high air preheating temperature. As a consequence, the reaction zone is found to be characteristically broad and the temperature and concentrations fields are sufficiently homogeneous and uniform, leading to a strong abatement of nitric oxide emissions. It is also observed that the reduction of thermal gradients is achieved mainly through the extension of combustion in the whole volume of the combustion chamber, so that a flame front no longer exists ('flameless oxidation'). The effect of preheating, further dilution provided by inner recirculation and of radiation model for the present hydrogen/air mild burner are analyzed.


Author(s):  
T. Fondelli ◽  
D. Massini ◽  
A. Andreini ◽  
B. Facchini ◽  
F. Leonardi

The reduction of fluid-dynamic losses in high speed gearing systems is nowadays increasing importance in the design of innovative aircraft propulsion systems, which are particularly focused on improving the propulsive efficiency. Main sources of fluid-dynamic losses in high speed gearing systems are windage losses, inertial losses resulting by impinging oil jets used for jet lubrication and the losses related to the compression and the subsequent expansion of the fluid trapped between gears teeth. The numerical study of the latter is particularly challenging since it faces high speed multiphase flows interacting with moving surfaces, but it paramount for improving knowledge of the fluid behavior in such regions. The current work aims to analyze trapping losses in a gear pair by means of three-dimensional CFD simulations. In order to reduce the numerical effort, an approach for restricting computational domain was defined, thus only a portion of the gear pair geometry was discretized. Transient calculations of a gear pair rotating in an oil-free environment were performed, in the context of conventional eddy viscosity models. Results were compared with experimental data from the open literature in terms of transient pressure within a tooth space, achieving a good agreement. Finally, a strategy for meshing losses calculation was developed and results as a function of rotational speed were discussed.


2012 ◽  
Vol 23 (04) ◽  
pp. 1250030 ◽  
Author(s):  
FAYÇAL HAMMAMI ◽  
NADER BEN-CHEIKH ◽  
ANTONIO CAMPO ◽  
BRAHIM BEN-BEYA ◽  
TAIEB LILI

In this work, a numerical study devoted to the two-dimensional and three-dimensional flow of a viscous, incompressible fluid inside a lid-driven cavity is undertaking. All transport equations are solved using the finite volume formulation on a staggered grid system and multi-grid acceleration. Quantitative aspects of two and three-dimensional flows in a lid-driven cavity for Reynolds number Re = 1000 show good agreement with benchmark results. An analysis of the flow evolution demonstrates that, with increments in Re beyond a certain critical value Rec, the steady flow becomes unstable and bifurcates into unsteady flow. It is observed that the transition from steadiness to unsteadiness follows the classical Hopf bifurcation. The time-dependent velocity distribution is studied in detail and the critical Reynolds number is localized for both 2D and 3D cases. Benchmark solutions for 2D and 3D lid-driven cavity flows are performed for Re = 1500 and 6000.


2009 ◽  
Vol 131 (5) ◽  
Author(s):  
Alfred von Loebbecke ◽  
Rajat Mittal ◽  
Frank Fish ◽  
Russell Mark

Three-dimensional fully unsteady computational fluid dynamic simulations of five Olympic-level swimmers performing the underwater dolphin kick are used to estimate the swimmer’s propulsive efficiencies. These estimates are compared with those of a cetacean performing the dolphin kick. The geometries of the swimmers and the cetacean are based on laser and CT scans, respectively, and the stroke kinematics is based on underwater video footage. The simulations indicate that the propulsive efficiency for human swimmers varies over a relatively wide range from about 11% to 29%. The efficiency of the cetacean is found to be about 56%, which is significantly higher than the human swimmers. The computed efficiency is found not to correlate with either the slender body theory or with the Strouhal number.


Author(s):  
Abdulrahman Alenezi ◽  
Abdulrahman Almutairi ◽  
Hamad Alhajeri ◽  
Abdulaziz Gamil ◽  
Faisal Alshammari

Abstract A detailed heat transfer numerical study of a three-dimensional impinging jet on a roughened isothermal surface is presented and is investigated from flow physics vantage point under the influence of different parameters. The effects of the Reynolds number, roughness location, and roughness dimension on the flow physics and heat transfer parameters are studied. Additionally, the relations between average heat transfer coefficient (AHTC) and flow physics including pressure, wall shear and flow vortices with thermodynamic nonequilibrium are offered. This paper studies the effect of varying both location and dimension of the roughness element which took the shape of square cross-sectional continuous ribs to deliver a favorable trade-off between total pressure loss and heat transfer rate. The roughness element was tested for three different radial locations (R/D) = 1, 1.5, and 2 and at each location its height (i.e., width) (e) was changed from 0.25 to 1 mm in incremental steps of 0.25. The study used a jet angle (α) of 90 deg, jet-to-target distance (H/D = 6), and Re ranges from 10,000 to 50,000, where H is the vertical distance between the target plate and jet exit. The results show that the AHTC can be significantly affected by changing the geometry and dimensions of the roughness element. This variation can be either an augmentation of, or decrease in, the (HTC) when compared with the baseline case. An enhancement of 12.9% in the AHTC was achieved by using optimal location and dimensions of the roughness element at specific Reynolds number. However, a diminution between 10% and 30% in (AHTC) was attained by the use of rib height e = 1 mm at Re = 50k. The variation of both rib location and height showed better contribution in increasing heat transfer for low-range Reynolds numbers.


Author(s):  
Azita Soleymani ◽  
Eveliina Takasuo ◽  
Piroz Zamankhan ◽  
William Polashenski

Results are presented from a numerical study examining the flow of a viscous, incompressible fluid through random packing of nonoverlapping spheres at moderate Reynolds numbers (based on pore permeability and interstitial fluid velocity), spanning a wide range of flow conditions for porous media. By using a laminar model including inertial terms and assuming rough walls, numerical solutions of the Navier-Stokes equations in three-dimensional porous packed beds resulted in dimensionless pressure drops in excellent agreement with those reported in a previous study (Fand et al., 1987). This observation suggests that no transition to turbulence could occur in the range of Reynolds number studied. For flows in the Forchheimer regime, numerical results are presented of the lateral dispersivity of solute continuously injected into a three-dimensional bounded granular bed at moderate Peclet numbers. Lateral fluid dispersion coefficients are calculated by comparing the concentration profiles obtained from numerical and analytical methods. Comparing the present numerical results with data available in the literature, no evidence has been found to support the speculations by others for a transition from laminar to turbulent regimes in porous media at a critical Reynolds number.


Sign in / Sign up

Export Citation Format

Share Document