Menu

Volume-of-Fluid Simulations of Gas-Liquid-Liquid Flows in Microchannels

Volume-of-Fluid Simulations of Gas-Liquid-Liquid Flows in Microchannels

Authors: 
Buwa, V. - Presenter, Indian Institute of Technology Delhi
Rajesh, V. M. - Presenter, Shiv Nadar University




v\:* {behavior:url(#default#VML);} o\:* {behavior:url(#default#VML);} w\:* {behavior:url(#default#VML);} .shape {behavior:url(#default#VML);}




Print
Clean
Clean
false


false
false
false

EN-US
X-NONE
X-NONE




MicrosoftInternetExplorer4














DefSemiHidden="false" DefQFormat="false" DefPriority="99"
LatentStyleCount="371">


UnhideWhenUsed="true" QFormat="true" Name="heading 2"/>
UnhideWhenUsed="true" QFormat="true" Name="heading 3"/>
UnhideWhenUsed="true" QFormat="true" Name="heading 4"/>
UnhideWhenUsed="true" QFormat="true" Name="heading 5"/>
UnhideWhenUsed="true" QFormat="true" Name="heading 6"/>
UnhideWhenUsed="true" QFormat="true" Name="heading 7"/>
UnhideWhenUsed="true" QFormat="true" Name="heading 8"/>
UnhideWhenUsed="true" QFormat="true" Name="heading 9"/>
Name="index 1"/>
Name="index 2"/>
Name="index 3"/>
Name="index 4"/>
Name="index 5"/>
Name="index 6"/>
Name="index 7"/>
Name="index 8"/>
Name="index 9"/>
UnhideWhenUsed="true" Name="toc 1"/>
UnhideWhenUsed="true" Name="toc 2"/>
UnhideWhenUsed="true" Name="toc 3"/>
UnhideWhenUsed="true" Name="toc 4"/>
UnhideWhenUsed="true" Name="toc 5"/>
UnhideWhenUsed="true" Name="toc 6"/>
UnhideWhenUsed="true" Name="toc 7"/>
UnhideWhenUsed="true" Name="toc 8"/>
UnhideWhenUsed="true" Name="toc 9"/>
Name="Normal Indent"/>
Name="footnote text"/>
Name="annotation text"/>
Name="header"/>
Name="footer"/>
Name="index heading"/>
UnhideWhenUsed="true" QFormat="true" Name="caption"/>
Name="table of figures"/>
Name="envelope address"/>
Name="envelope return"/>
Name="footnote reference"/>
Name="annotation reference"/>
Name="line number"/>
Name="page number"/>
Name="endnote reference"/>
Name="endnote text"/>
Name="table of authorities"/>
Name="macro"/>
Name="toa heading"/>
Name="List"/>
Name="List Bullet"/>
Name="List Number"/>
Name="List 2"/>
Name="List 3"/>
Name="List 4"/>
Name="List 5"/>
Name="List Bullet 2"/>
Name="List Bullet 3"/>
Name="List Bullet 4"/>
Name="List Bullet 5"/>
Name="List Number 2"/>
Name="List Number 3"/>
Name="List Number 4"/>
Name="List Number 5"/>

Name="Closing"/>
Name="Signature"/>
UnhideWhenUsed="true" Name="Default Paragraph Font"/>
Name="Body Text"/>
Name="Body Text Indent"/>
Name="List Continue"/>
Name="List Continue 2"/>
Name="List Continue 3"/>
Name="List Continue 4"/>
Name="List Continue 5"/>
Name="Message Header"/>

Name="Salutation"/>
Name="Date"/>
Name="Body Text First Indent"/>
Name="Body Text First Indent 2"/>
Name="Note Heading"/>
Name="Body Text 2"/>
Name="Body Text 3"/>
Name="Body Text Indent 2"/>
Name="Body Text Indent 3"/>
Name="Block Text"/>
Name="Hyperlink"/>
Name="FollowedHyperlink"/>


Name="Document Map"/>
Name="Plain Text"/>
Name="E-mail Signature"/>
Name="HTML Top of Form"/>
Name="HTML Bottom of Form"/>
Name="Normal (Web)"/>
Name="HTML Acronym"/>
Name="HTML Address"/>
Name="HTML Cite"/>
Name="HTML Code"/>
Name="HTML Definition"/>
Name="HTML Keyboard"/>
Name="HTML Preformatted"/>
Name="HTML Sample"/>
Name="HTML Typewriter"/>
Name="HTML Variable"/>
Name="Normal Table"/>
Name="annotation subject"/>
Name="No List"/>
Name="Outline List 1"/>
Name="Outline List 2"/>
Name="Outline List 3"/>
Name="Table Simple 1"/>
Name="Table Simple 2"/>
Name="Table Simple 3"/>
Name="Table Classic 1"/>
Name="Table Classic 2"/>
Name="Table Classic 3"/>
Name="Table Classic 4"/>
Name="Table Colorful 1"/>
Name="Table Colorful 2"/>
Name="Table Colorful 3"/>
Name="Table Columns 1"/>
Name="Table Columns 2"/>
Name="Table Columns 3"/>
Name="Table Columns 4"/>
Name="Table Columns 5"/>
Name="Table Grid 1"/>
Name="Table Grid 2"/>
Name="Table Grid 3"/>
Name="Table Grid 4"/>
Name="Table Grid 5"/>
Name="Table Grid 6"/>
Name="Table Grid 7"/>
Name="Table Grid 8"/>
Name="Table List 1"/>
Name="Table List 2"/>
Name="Table List 3"/>
Name="Table List 4"/>
Name="Table List 5"/>
Name="Table List 6"/>
Name="Table List 7"/>
Name="Table List 8"/>
Name="Table 3D effects 1"/>
Name="Table 3D effects 2"/>
Name="Table 3D effects 3"/>
Name="Table Contemporary"/>
Name="Table Elegant"/>
Name="Table Professional"/>
Name="Table Subtle 1"/>
Name="Table Subtle 2"/>
Name="Table Web 1"/>
Name="Table Web 2"/>
Name="Table Web 3"/>
Name="Balloon Text"/>

Name="Table Theme"/>























Name="List Paragraph"/>

Name="Intense Quote"/>














































































Name="Subtle Emphasis"/>
Name="Intense Emphasis"/>
Name="Subtle Reference"/>
Name="Intense Reference"/>

UnhideWhenUsed="true" Name="Bibliography"/>
UnhideWhenUsed="true" QFormat="true" Name="TOC Heading"/>













Name="Grid Table 1 Light Accent 1"/>




Name="Grid Table 6 Colorful Accent 1"/>
Name="Grid Table 7 Colorful Accent 1"/>
Name="Grid Table 1 Light Accent 2"/>




Name="Grid Table 6 Colorful Accent 2"/>
Name="Grid Table 7 Colorful Accent 2"/>
Name="Grid Table 1 Light Accent 3"/>




Name="Grid Table 6 Colorful Accent 3"/>
Name="Grid Table 7 Colorful Accent 3"/>
Name="Grid Table 1 Light Accent 4"/>




Name="Grid Table 6 Colorful Accent 4"/>
Name="Grid Table 7 Colorful Accent 4"/>
Name="Grid Table 1 Light Accent 5"/>




Name="Grid Table 6 Colorful Accent 5"/>
Name="Grid Table 7 Colorful Accent 5"/>
Name="Grid Table 1 Light Accent 6"/>




Name="Grid Table 6 Colorful Accent 6"/>
Name="Grid Table 7 Colorful Accent 6"/>







Name="List Table 1 Light Accent 1"/>




Name="List Table 6 Colorful Accent 1"/>
Name="List Table 7 Colorful Accent 1"/>
Name="List Table 1 Light Accent 2"/>




Name="List Table 6 Colorful Accent 2"/>
Name="List Table 7 Colorful Accent 2"/>
Name="List Table 1 Light Accent 3"/>




Name="List Table 6 Colorful Accent 3"/>
Name="List Table 7 Colorful Accent 3"/>
Name="List Table 1 Light Accent 4"/>




Name="List Table 6 Colorful Accent 4"/>
Name="List Table 7 Colorful Accent 4"/>
Name="List Table 1 Light Accent 5"/>




Name="List Table 6 Colorful Accent 5"/>
Name="List Table 7 Colorful Accent 5"/>
Name="List Table 1 Light Accent 6"/>




Name="List Table 6 Colorful Accent 6"/>
Name="List Table 7 Colorful Accent 6"/>

/* Style Definitions */ table.MsoNormalTable {mso-style-name:"Table Normal"; mso-tstyle-rowband-size:0; mso-tstyle-colband-size:0; mso-style-noshow:yes; mso-style-priority:99; mso-style-parent:""; mso-padding-alt:0in 5.4pt 0in 5.4pt; mso-para-margin:0in; mso-para-margin-bottom:.0001pt; mso-pagination:widow-orphan; font-size:10.0pt; font-family:"Calibri","sans-serif";}







justify;line-height:normal;text-autospace:none"> mso-ansi-language:EN-IN'>Microreactors EN-IN'> are miniaturized and continuous devices, often comprised of several
channels (of different size and shape) with characteristic dimensions of the
order of 10 - 1000 μm. Microreactors
offer several advantages and the key advantage is that they offer a large
interfacial area (per unit volume) in comparison to conventional reactors.
  color:black;mso-ansi-language:EN-IN">Micro−reactors have been used to
intensify processes involving gas−liquid and liquid−liquid reactions.
However, applications of microreactors for catalytic reactions
involving gas−liquid−liquid (G/L/L) flows still remain unexplored.
An extensive literature survey showed that several experimental and numerical
studies were performed to understand two−phase gas−liquid and liquid−liquid
flows in microchannels ( EN-IN">e.g. Shao et al., 2010; Raj et al., 2010 EN-IN">). However, only a very few experimental studies were reported on
characterization of hydrodynamics of gas−liquid−liquid flows in microchannels ( 14.0pt;font-family:TimesNewRomanPSMT;color:#0070C0;mso-ansi-language:EN-IN">e.g.
Khan and Duraiswamy, 2009; Wang et al., 2010 and Yue
et al., 2014 color:black;mso-ansi-language:EN-IN">) 14.0pt;font-family:TimesNewRomanPSMT;color:#241F20;mso-ansi-language:EN-IN">.

justify;line-height:normal;text-autospace:none"> EN-IN">For development of microreactors to intensify
processes involving three−phase 14.0pt;font-family:TimesNewRomanPSMT;color:black;mso-ansi-language:EN-IN">G/L/L mso-ansi-language:EN-IN"> reactions, it is important to understand the
hydrodynamics of TimesNewRomanPSMT;color:black;mso-ansi-language:EN-IN">G/L/L mso-ansi-language:EN-IN"> flows in microchannels. In
particular, it is important to understand three−phase flow regimes and
formation of bubbles/drops or slugs of gas and liquid into another liquid.
Moreover, it is important to understand how to manipulate distributor
configurations and flow rates to produce desired flow patterns and to control
bubble/drop or slug sizes. While there exists significant literature on
numerical simulations of gas−liquid and liquid−liquid flows using
various numerical methods e.g. the Volume−of−Fluid (VOF) or the
Level Set or the Lattice−Boltzmann method; to the best of author?s
knowledge, numerical simulations of three−phase flows using interface
tracking/capturing methods have not been attempted so far. The objective of the
present contribution is to perform the VOF simulations of gas (air) and liquid
(water) in another liquid (kerosene) 14.0pt;font-family:TimesNewRomanPSMT;color:black;mso-ansi-language:EN-IN">immiscible
with the first liquid font-family:TimesNewRomanPSMT;color:#241F20;mso-ansi-language:EN-IN">flowing continuously
in a microchannel for different flow regimes and mso-ansi-language:EN-IN">understand the formation dynamics of bubbles/drops mso-ansi-language:EN-IN"> font-family:TimesNewRomanPSMT;color:black;mso-ansi-language:EN-IN">or slugs
under different flow conditions and to verify the predictions using experiments.

justify;line-height:normal;text-autospace:none">Experiments
were performed using a single T−junction microchannel
with a rectangular cross−section (H = 950 μm;
W = 1000 μm, Linlet
= 30 mm, Lchannel =210 mm) machined on a
PMMA sheet (Figure 1a). The machined plate was EN-IN'>glued with another PMMA sheet with a thin glue sheet. Air and water were
pumped from the side inlets and kerosene was pumped through the main channel
(Figure 1(a)), to generate bubbles/drops or slugs of air and water in kerosene
(continuous phase). The Sodium Dodecyl Sulfate (SDS) was used as the surfactant in the present
work with the concentrations of 0.3 and 2 wt/wt%.  The surface/interfacial tensions of air−water (air−water+SDS
(0.3 wt/wt %) = 0.035; air−water+SDS (2 wt/wt %) = 0.031) and oil-water (oil−
water+SDS (0.3 wt/wt %) = 0.005; oil−water+SDS
(2 wt/wt %) = 0.003) were measured using Tensiometer
(Kruss GmbH, Germany). A high − speed
digital camera was used (at 2000 fps) to visualize the formation dynamics of
bubble / drop / slug and the G/L/L/ flow regimes. These observed flow regimes
and formation mechanisms of bubble / drop or slugs and their lengths were
further used to verify the numerical predictions.

justify;line-height:normal;text-autospace:none"> mso-ansi-language:EN-IN'> 

justify;line-height:normal;text-autospace:none"> font-family:TimesNewRomanPSMT;mso-no-proof:yes">

justify;line-height:normal;text-autospace:none"> EN-IN'>Numerical simulations were performed using the VOF method implemented in commercial CFD solver Fluent
14 (Ansys Inc, USA) to
simulate the formation dynamics of bubbles/drops in different flow regimes for
0.3 and 2 wt/wt% SDS
mixture. The channel geometry and dimensions (see Figure 1 (b)) used in the simulations were the same as that
considered in the experiments (1000 µm x 950 µm). Simulations were performed
with air as the gas phase (ρair =
1.78 kg/m3, μair= 1.37 x
10−5Pa.s) and water+SDS
as dispersed liquid phase (ρwater =
998.2 kg/m3, μwater= 0.001
Pa.s) and kerosene as the continuous liquid phase (ρkerosene = 780 kg/m3, μkerosene= 0.00115 Pa.s).
The volumetric flow rates of phases were used in the
simulations ranges of Qoil = 3?15 ml/min; Qair= 5.24 ? 11.18 ml/min; Qwater =
0.85?5.61 ml/min. The predictions were verified using the experiments. In this
abstract few key results are reported. The volumetric flow rates of phases were
used in the simulations ranges of Qoil = 3?15 ml/min; Qair= 5.24 ? 11.18 ml/min; Qwater = 0.85?5.61 ml/min. The predictions were verified
using the experiments. In this abstract few key results are reported. A
time step of 5x10−6 s was used for the numerical simulations.
 Further details of the computational model and numerics
will be provided in the full length manuscript.

justify;line-height:normal;text-autospace:none"> mso-ansi-language:EN-IN'>Figure. font-family:"Times New Roman","serif";color:red;mso-ansi-language:EN-IN'>2a color:black;mso-ansi-language:EN-IN'>shows the predicted formation dynamics of
air slug (by iso surface of volume fraction) at 1stT mso-ansi-language:EN-IN'>−junction for Qoil=
5 ml/min for 0.3 wt/wt % of
water+SDS mixture. As shown in Figure. 2a the magnitude of shear
force exerted on the air−oil interface
at 1stT−junction was less
compared to the surface tension force which opposes the slug movement and as a
result slugs were grown and occupy the main channel cross-section and that led
to obstruct the continuous phase flow. Then the pressure was build up in the downstream
of air−oil interface at the 1stT−junction. The increase in pressure drop across the slug led
to its squeezing and the slug was finally detached at the 1stT−junction as shown at t = 21 ms in
Figure. 2a. This
regime was referred as the squeezing regime. In addition, simulations of bubble
formation in dripping formation regime were performed (results not shown here).

justify;line-height:normal;text-autospace:none"> mso-ansi-language:EN-IN'>The pressure in the dispersed and continuous phase was
identified as essential element to understand the formation dynamics of bubbles/drops
or slugs. In the present work, simulated point pressure at both 1stT mso-ansi-language:EN-IN'>−junction (see Figure.
3(a)) and 2nd
T−junction (see Figure. 3(b)) was considered to
analyse the formation dynamics of bubbles/drops or slugs. The time history of
pressures was recorded at each time step at the aforementioned points in the
solution domain for Qoil = 5 and 15 ml/min as shown in Figure. 3.  The detailed
pressure analysis (see Figure. 3 (b)) during the squeezing and dripping regime will be reported
in the full length manuscript.

normal;text-autospace:none"> mso-no-proof:yes">

justify;line-height:normal;text-autospace:none"> font-family:"Times New Roman","serif";color:black;mso-no-proof:yes'>

justify;line-height:normal;text-autospace:none"> mso-ansi-language:EN-IN'> 

justify;line-height:normal;text-autospace:none"> font-family:TimesNewRomanPSMT;mso-no-proof:yes">

justify;line-height:normal;text-autospace:none"> mso-ansi-language:EN-IN'>Figure. font-family:"Times New Roman","serif";color:red;mso-ansi-language:EN-IN'>4(a) color:black;mso-ansi-language:EN-IN'>shows the predicted formation dynamics of
a water slug at the 2 font-family:"Times New Roman","serif";mso-ansi-language:EN-IN'>nd mso-ansi-language:EN-IN'> font-family:"Times New Roman","serif";mso-ansi-language:EN-IN'>T−junction for 0.3 wt/wt
% of water+SDS mixture. Figure. 4(b) shows the time
evolution of pressure at point P2. Air slug formed at the 1st
T−junction led the detachment of water+SDS−oil
interface at the 2ndT−junction as shown in Figure. 4(a). At t =0 ms, rear end of air slug
was attached to the water-oil interface and water slug started entering into
the main channel. It should be noted that addition of 0.3 wt/wt % SDS led to reduction in the interfacial tension force
between water+SDS and oil (0.005 N/m) compared to
that between the pure water and oil (0.046 N/m). As time progressed, the
deformed water+SDS-oil interface was dragged into the
main channel as shown at t = 4 - 8 ms. Then the water+SDS-oil interface at the 2nd T−junction was squeezed by the incoming air slug as shown at
t = 14 and 16 ms and the corresponding interface was
detached at t = 18 ms as shown in Figure 4(a). The corresponding
experimental image is shown as the first image in Figure. 5(c) and the predicted
detachment mechanism and water slug shape agreed very well with the experiments.

justify;line-height:normal;text-autospace:none"> mso-ansi-language:EN-IN'>A comparison of the experimental and predicted flow
regimes for SDS concentrations in aqueous phase of 0.3 wt/wt% is shown in Figure. mso-ansi-language:EN-IN'>5 font-family:"Times New Roman","serif";color:black;mso-ansi-language:EN-IN'>. In the images shown in Figure. mso-ansi-language:EN-IN'>5 (a) ? 5(d), mso-ansi-language:EN-IN'>the first two images correspond to the experiments, in
which the first image was captured at the 2nd T mso-ansi-language:EN-IN'>−junction, and second
image was captured at a distance of 4.6 cm downstream the 2nd T−junction. The corresponding predicted flow regime is shown
in the third image of Figure. 5(a) ? 5(d). In experiments and numerical investigations, air
bubbles/slugs were formed at the 1st T−junction (not visible in the figures) and in all situations
these bubbles/slugs lead to the detachment of water drops/slugs at the 2nd
T−junction as shown in the first images
of Figure. 5(a) ? 5(d). The predicted three−phase
flow regimes (third image) showed an excellent agreement with an experimental
flow regimes (1st image) viz. the Bubble−Drop (B−D) (Figure. 5(a)) and Slug−Slug (S−S) (Figure.5 (b) ? 5(d)) flow regimes. The corresponding 3D images (of predicted
flow regimes) are shown in Figure. 6.

justify;line-height:normal;text-autospace:none"> font-family:TimesNewRomanPSMT;mso-no-proof:yes">

normal;text-autospace:none">

justify;line-height:normal;text-autospace:none"> EN-IN'>While the predicted flow regimes and detachment mechanisms agreed well
with those observed in experiments, a quantitative verification was performed
by comparing the predicted and measured lengths of bubbles/drops or slugs with the effects of Qoil (figure not shown here),
Qair (see Figure. 7) and Qwater
(see Figure. 8).
The predictions showed quantitative agreement with experiments as shown in
Figure. 7 and Figure. 8. The detailed analysis
(with the help of dimensionless numbers such as Caoil,
Weair and Wewater)
of results obtained from the present work will be reported in the full length
manuscript. The understanding of formation mechansims
and flow regimes, experimentally validated computational model is expected to
help in devising distributors and optimizing flow conditions to control
bubble/drops/slug sizes and flow patterns in microchannels.
   

normal;text-autospace:none">

justify;line-height:normal;text-autospace:none"> normal"> yes">

justify;line-height:normal;text-autospace:none">References

justify;line-height:normal;text-autospace:none"> normal"> yes">