Applicability of two kinds of computational-fluid-dynamics method adopting Cahn-Hilliard (CH) and Allen-Cahn (AC)-type diffuse-interface advection equations based on a phase-field model (PFM) is examined to simulation of motions of microscopic incompressible two-phase fluid on solid surface. A capillarity-driven gas-liquid motion in rectangular channel is simulated by use of a PFM method for solving Navier-Stokes (NS) equations and a CH equation, whereas an immiscible liquid-liquid flow in a microchannel with T-junction and square cross section is simulated by use of another PFM method proposed in this study, which adopts a lattice-Boltzmann method based on fictitious particles kinematics as numerical scheme for solving NS equations and an AC equation that is modified to improve volume-of-fluid conservation. The major findings are as follows: (1) effect of capillary force on the dynamic two-phase fluid system with a high density ratio is well predicted for cross-sectional aspect ratio of the channel = 1 and 2; (2) mono-dispersed slug flow pattern transition is reproduced in good agreement with experimental observations in terms of variation in length and interval of droplets as increasing their volumetric flow rates at a constant flow rate ratio = 1. These results prove that the PFM methods can be used for analyzing two-phase fluid motions in various microfluidic devices and micro fabrication processes.
References
1.
NorthrupM. A.JensenK. F.HarrisonD. J. eds. Micro Total Analysis Systems 2003 (Transducers Research Foundation, USA, 2003).
2.
MatsumotoS.TaiY. -C., A nL/min-Scale Valveless Flow Regulator for Integrated Micro Chemical Devices, JSME Int. J. Ser. B Fluids Therm. Eng., 2004, 47 (3), 528–533.
3.
WangY.PetersonG. P., Investigation of a Novel Flat Heat Pipe, J. Heat Transfer, 2005, 127 (2), 165–170.
4.
BerthierJ.ClementzPh.RaccurtO.JaryD.ClaustreP.PeponnetC.FouilletY., Computer Aided Design of an EWOD Microdevice, Sens. Actuator A – Phys., 2006, Vol. 127 (2), 283–294.
5.
FeenstraB. J.HayesR. A.van DijkR.BoomR. G. H.WagemansM. M. H.CampsI. G. J.GiraldoA.v. d. HeijdenB., Electrowetting-based Displays: Bringing Microfluidics Alive On-screen, Proceedings of MEMS 2006, 2006, 48–53.
6.
MatsumotoS.TakadaN.MatsumotoJ., Microfabrication Process of Cellular Structures in Hollow Fiber-shaped Substrates, Transactions of the Japan Society of Mechanical Engineers Series C, 2010, 76 (768), 1911–1913.
7.
LuJ.TakagiH.MaedaR., Chip to Wafer Temporary Bonding with Self-alignment by Patterned FDTS Layer for Size-free MEMS Integration, Proceedings of IEEE Sensors 2011 Conference, 2011, 1121–1124.
8.
InoueT.OhtakiK.MurakamiS.MatsumotoS., Direct Synthesis of Hydrogen Peroxide Based on Microreactor Technology, Fuel Process. Technol., 2013, 108, 8–11.
9.
AndersonD. M.McFaddenG. B.WheelerA. A., Diffuse-interface Methods in Fluid Mechanics. Annu. Rev. Fluid Mech., 1998, 30, 139–165.
10.
NadigaB. T.ZaleskiS., Investigations of a Two-phase Fluid Model, Eur. J. Mech. B Fluids, 1996, 15 (6), 885–896.
11.
JacqminD., Calculation of Two-phase Navier–Stokes Flows Using Phase-field Modelling, J. Comput. Phys, 1999, 155 (1), 96–127.
12.
JametD.LebaiqueO.CoutrisN.DelhayeJ. M., The Second Gradient Method for the Direct Numerical Simulation of Liquid-Vapor Flows with Phase Change, J. Comput. Phys., 2001, 169 (2), 624–651.
13.
JametD.TorresD.BrackbillJ. U., On the Theory and Computation of Surface Tension: The Elimination of Parasitic Currents through Energy Conservation in the Second-Gradient Method, J. Comput. Phys., 2002, 182 (1), 262–276.
14.
BadalassiV. E.CenicerosH. D.BanerjeeS., Computation of Multiphase Systems with Phase Field Model, J. Comput. Phys., 2003, 190 (2), 371–397.
15.
KimJ., A Continuous Surface Tension Force Formulation for Diffuse-interface Models, J. Comput. Phys., 2005, 204 (2), 784–804.
16.
KimJ., Phase Field Computations for Ternary Fluid Flows, Comput. Meth. Appl. Mech. Eng., 2007, 196 (45–48), 4779–4788.
17.
van der WaalsJ. E., The Thermodynamic Theory of Capillarity under the Hypothesis of a Continuous Density Variation, Transl. RowlinsonJ. S., J. Stat. Phys., 1979, 20, 197–244.
18.
CahnJ. W.HilliardJ. E., Free Energy of a Nonuniform System. I. Interfacial Free Energy, J. Chem. Phys., 1958, 28 (2), 258–267.
19.
SeppecherP., Moving Contact Lines in the Cahn-Hilliard Theory, Int. J. Eng. Sci., 1996, 34 (9), 977–992.
20.
BriantA. J.PapatzacosP.YeomansJ. M., Lattice Boltzmann Simulations of Contact Line Motion in a Liquid-Gas System, Phil. Trans. R. Soc. Lond. A, 2002, 360, 485–495.
21.
YoshinoM.MizutaniY., Lattice Boltzmann Simulation of Liquid-Gas Flows through Solid Bodies in a Square Duct, Math. Comput. Simul., 2006, 72 (2–6), 264–269.
TakadaN.TomiyamaA., A Numerical Method for Two-Phase Flow Based on a Phase-Field Model. JSME Int. J. Ser. B Fluids Therm. Eng., 2006, 49 (3), 636–644.
24.
TakadaN.MisawaM.TomiyamaA., A Phase-Field Method for Interface-Tracking Simulation of Two-Phase Flows. Math. Comput. Simul., 2006, 72 (2–6), 220–226.
25.
TakadaN.TomiyamaA., Application of Interface-Tracking Method Based on Phase-Field Model to Numerical Analysis of Free Surface Flow, Theoretical and Applied Mechanics Japan, 2006, 55, 149–156.
26.
MatsumotoJ.TakadaN.Two-phase Flow Analysis Based on a Phase-field Model Using Orthogonal Basis Bubble Function Finite Element Method, Int. J. Comput. Fluid Dyn., 2008, 22 (8), 555–568.
27.
TakadaN.MatsumotoJ.MatsumotoS.IchikawaN., Application of a Phase-field Method to the Numerical Analysis of Motions of a Two-phase Fluid with High Density Ratio on a Solid Surface, Journal of Computational Science and Technology, 2008, 2 (2), 318–329.
28.
TakadaN.MatsumotoJ.MatsumotoS., Phase-field Model-based Simulation of Motions of a Two-phase Fluid on Solid Surface, Journal of Computational Science and Technology, 2013, 7 (2), 322–337.
29.
ChenS.DoolenG. D., Lattice Boltzmann Method for Fluid Flows, Annu. Rev. Fluid Mech., 1998, 30, 329–364.
30.
SucciS., The Lattice Boltzmann Equation for Fluid Dynamics and Beyond, Oxford at the Clarendon Press, UK, 2001.
31.
SwiftM. R.OrlandiniE.OsbornW. R.YeomansJ. M., Lattice Bolztammn Simulations of Liquid-Gas and Binary Fluid Systems, Phys. Rev. E, 1996, 54 (5), 5041–5052.
32.
TakadaN.MisawaM.TomiyamaA.HosokawaS., Simulation of Bubble Motion under Gravity by Lattice Boltzmann Method, J. Nucl. Sci. Technol., 2001, 38 (5), 330–341.
33.
InamuroT.OgataT.TajimaS.KonishiN., A Lattice Boltzmann Method for Incompressible Two-phase Flows with Large Density Differences, J. Comput. Phys., 2004, 198 (2), 628–644.
34.
YanY. Y.ZuY. Q., A Lattice Boltzmann Method for Incompressible Two-phase Flows on Partial Wetting Surface with Large Density Ratio, J. Comput. Phys., 2007, 227 (1), 763–775.
35.
TengS.ChenY.OhashiH., Lattice Boltzmann Simulation of Multiphase Fluid Flows through the Total Variation Diminishing with Artificial Compression Scheme, Int. J. Heat Fluid Flow, 2000, 21 (1), 112–121.
36.
SetaT.KonoK., Thermal Lattice Boltzmann Method for Liquid-Gas Two-phase Flows in Two Dimensions, JSME Int. J. Ser. B Fluids Therm. Eng., 2004, 47 (3), 572–583.
37.
SetaT.OkuiK., The Single Component Thermal Lattice Boltzmann Simulation of Pool Boiling in Two Dimensions, J. Therm. Sci. Technol., 2006, 1 (2), 125–137.
38.
ChorinA. J., Numerical Solution of the Navier-Stokes Equations, Math. Comput., 1968, 22 (104), 745–762.
39.
KawamuraT.KuwaharaK., Computation of High Reynolds Number Flow around a Circular Cylinder with Surface Roughness, AIAA Paper, 1984, 84–0340.
40.
ChiuP.-H.LinY.-T., A Conservative Phase Field Method for Solving Incompressible Two-phase Flows, J. Comput. Phys., 2011, 230 (1), 185–204.
41.
BibenT.MisbahC., Tumbling of Vesicles under Shear Flow within an Advected-field Approach, Phys. Rev. E., 2003, 67 (3), 031908.
42.
BeaucourtJ.BibenT.LeyatA.VerdierC., Modeling Breakup and Relaxation of Newtonian Droplets Using the Advected Phase-field Approach, Phys. Rev. E., 2007, 75 (2), 021405.
43.
OlssonE.KreissG., A Conservative Level Set Method for Two Phase Flow, J. Comput. Phys., 2005, 210 (1), 225–246.
44.
TanN.AokiT.InoueK.YoshitaniK., Numerical Simulation of Two-phase Flow Driven by Rotating Object, Transactions of the Japan Society of Mechanical Engineers Series B, 2011, 77 (781), 1699–1714.
45.
SatoY.NicenoB., A New Contact Line Treatment for a Conservative Level Set Method, J. Comput. Phys., 2012, 231 (10), 3887–3895.
46.
TanakaY.YoshinoM.HirataT., Lattice Boltzmann Simulation of Nucleate Pool Boiling in Saturated Liquid, Commun. Comput. Phys., 2011, 9 (5), 1347–1361.
47.
OrlandiniE.SwiftM. R.YeomansJ. M., A Lattice Boltzmann Model of Binary-fluid Mixtures, Europhys. Lett., 1995, 32 (6), 463–468.
48.
QianY. H.OrszagS. A., Scalings in Reaction A+B?C: Numerical Simulation by Lattice Boltzmann Models, J. Stat. Phys., 1995, 81 (1&2), 237–253.
49.
HirabayashiM.ChenY.OhashiH., The Lattice BGK Model for the Poisson Equation, JSME Int. J. Ser. B Fluids Therm. Eng., 2001, 44 (1), 45–52.
50.
IchikawaN.HosokwaK.MaedaR., Interface Motion of Capillary-driven Flow in Rectilinear Microchannel, J. Colloid Interface Sci., 2004, 280 (1), 155–164.
51.
GarsteckiP.FuerstmanM. J.StoneH. A.WhitesidesG. M., Formation of Droplets and Bubbles in a Microfluidic T-junction – Scaling and Mechanism of Break-up, Lab Chip, 2006, 6 (3), 437–446.
52.
KawaharaA.SadatomiM.NeiK.MatsuoH., Experimental Study on Bubble Velocity, Void Fraction and Pressure Drop for Gas-Liquid Two-phase Flow in a Circular Microchannel, Int. J. Heat Fluid Flow, 2009, 30 (5), 831–841.
53.
NakataM.SatoM.MatsuoY.MaedaS.HayashiS., Hollow Fibers Containing Various Display Elements: A Novel Structure for Electronic Paper, Journal of the Society for Information Display, 2006, 14 (8), 723–728.
