Research Papers: Heat and Mass Transfer

Modified Outlet Boundary Condition Schemes for Large Density Ratio Lattice Boltzmann Models

[+] Author and Article Information
Long Li

School of Mechanical Engineering,
Shanghai Jiao Tong University,
Shanghai 200240, China
e-mail: lilong315@sjtu.edu.cn

Xiaodong Jia

School of Mechanical Engineering,
Shanghai Jiao Tong University,
Shanghai 200240, China
e-mail: jxd_sjtu@126.com

Yongwen Liu

School of Mechanical Engineering,
Shanghai Jiao Tong University,
Shanghai 200240, China
e-mail: ywliu@sjtu.edu.cn

1Corresponding author.

Presented at the 2016 ASME 5th Micro/Nanoscale Heat & Mass Transfer International Conference. Paper No. MNHMT2016-6374.

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received June 15, 2016; final manuscript received January 29, 2017; published online March 7, 2017. Assoc. Editor: Robert D. Tzou.

J. Heat Transfer 139(5), 052003 (Mar 07, 2017) (8 pages) Paper No: HT-16-1391; doi: 10.1115/1.4036001 History: Received June 15, 2016; Revised January 29, 2017

Outlet boundary conditions (OBCs) and their numerical descriptions are critical to computational fluid dynamics (CFD) since they have significant influence on the numerical accuracy and stability. They present significant challenges to the two-phase lattice Boltzmann (LB) method, especially in the limit of large density ratio. In this study, three commonly used OBCs: convection boundary condition (CBC), Neumann boundary condition (NBC), and extrapolation boundary condition (EBC), are investigated and improved on basis of two LB models for large density ratios (single and double distribution function models). The existing numerical schemes for the OBCs are not directly applicable to the LB models because of the deviation of the momentum balance at the outlet boundary. The deviation becomes substantial at a large density ratio. Thus, in this work, modified OBC schemes are proposed to make the OBCs suitable for the two-phase LB models by adding an independent equation to obtain the outlet velocity. Numerical tests on droplet flowing in a channel are performed to evaluate the performance of the modified OBC schemes. Results indicate that the modified OBC schemes may be extended to tackle large density ratio situations. The modified NBC and EBC schemes are only suitable for the LB model with single distribution function. Three modified CBC schemes exhibit optimum performance for both single and double distribution function LB models which can be implemented for large density ratios.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


He, Y. , Li, Q. , Wang, Y. , and Tang, G. , 2009, “ Lattice Boltzmann Method and Its Applications in Engineering Thermophysics,” Chin. Sci. Bull., 54(22), pp. 4117–4134. [CrossRef]
Aidun, C. K. , and Clausen, J. R. , 2010, “ Lattice-Boltzmann Method for Complex Flows,” Annu. Rev. Fluid Mech., 42(1), pp. 439–472. [CrossRef]
Zhang, J. , 2010, “ Lattice Boltzmann Method for Microfluidics: Models and Applications,” Microfluid. Nanofluid., 10(1), pp. 1–28. [CrossRef]
Chen, L. , Kang, Q. , Mu, Y. , He, Y.-L. , and Tao, W.-Q. , 2014, “ A Critical Review of the Pseudopotential Multiphase Lattice Boltzmann Model: Methods and Applications,” Int. J. Heat Mass Transfer, 76, pp. 210–236. [CrossRef]
Li, Q. , Luo, K. H. , Kang, Q. J. , He, Y. L. , Chen, Q. , and Liu, Q. , 2016, “ Lattice Boltzmann Methods for Multiphase Flow and Phase-Change Heat Transfer,” Prog. Energy Combust. Sci., 52, pp. 62–105. [CrossRef]
He, X. , Chen, S. , and Zhang, R. , 1999, “ A Lattice Boltzmann Scheme for Incompressible Multiphase Flow and Its Application in Simulation of Rayleigh–Taylor Instability,” J. Comput. Phys., 152(2), pp. 642–663. [CrossRef]
Xing, X. Q. , Butler, D. L. , and Yang, C. , 2007, “ Lattice Boltzmann-Based Single-Phase Method for Free Surface Tracking of Droplet Motions,” Int. J. Numer. Methods Fluids, 53(2), pp. 333–351. [CrossRef]
Fakhari, A. , and Lee, T. , 2013, “ Multiple-Relaxation-Time Lattice Boltzmann Method for Immiscible Fluids at High Reynolds Numbers,” Phys. Rev. E, 87(2), p. 023304. [CrossRef]
Sun, K. , Jia, M. , and Wang, T. , 2013, “ Numerical Investigation of Head-On Droplet Collision With Lattice Boltzmann Method,” Int. J. Heat Mass Transfer, 58(1–2), pp. 260–275. [CrossRef]
Varmazyar, M. , and Bazargan, M. , 2011, “ Modeling of Free Convection Heat Transfer to a Supercritical Fluid in a Square Enclosure by the Lattice Boltzmann Method,” ASME J. Heat Transfer, 133(2), p. 022501. [CrossRef]
Li, L. , Mei, R. , and Klausner, J. F. , 2013, “ Heat Transfer Evaluation on Curved Boundaries in Thermal Lattice Boltzmann Equation Method,” ASME J. Heat Transfer, 136(1), p. 012403. [CrossRef]
Liu, X. , and Cheng, P. , 2013, “ Lattice Boltzmann Simulation for Dropwise Condensation of Vapor Along Vertical Hydrophobic Flat Plates,” Int. J. Heat Mass Transfer, 64, pp. 1041–1052. [CrossRef]
Li, Q. , Luo, K. H. , and Li, X. J. , 2013, “ Lattice Boltzmann Modeling of Multiphase Flows at Large Density Ratio With an Improved Pseudopotential Model,” Phys. Rev. E, 87(5), p. 053301. [CrossRef]
Liang, H. , Shi, B. C. , Guo, Z. L. , and Chai, Z. H. , 2014, “ Phase-Field-Based Multiple-Relaxation-Time Lattice Boltzmann Model for Incompressible Multiphase Flows,” Phys. Rev. E, 89(5), p. 053320. [CrossRef]
Xu, A. , Zhao, T. S. , An, L. , and Shi, L. , 2015, “ A Three-Dimensional Pseudo-Potential-Based Lattice Boltzmann Model for Multiphase Flows With Large Density Ratio and Variable Surface Tension,” Int. J. Heat Fluid Flow, 56, pp. 261–271. [CrossRef]
Cristea, A. , and Sofonea, V. , 2003, “ Reduction of Spurious Velocity in Finite Difference Lattice Boltzmann Models for Liquid–Vapor Systems,” Int. J. Mod. Phys. C, 14(09), pp. 1251–1266. [CrossRef]
Yuan, P. , and Schaefer, L. , 2006, “ Equations of State in a Lattice Boltzmann Model,” Phys. Fluids, 18(4), p. 042101. [CrossRef]
Lee, T. , and Lin, C.-L. , 2005, “ A Stable Discretization of the Lattice Boltzmann Equation for Simulation of Incompressible Two-Phase Flows at High Density Ratio,” J. Comput. Phys., 206(1), pp. 16–47. [CrossRef]
Lee, T. , and Fischer, P. F. , 2006, “ Eliminating Parasitic Currents in the Lattice Boltzmann Equation Method for Nonideal Gases,” Phys. Rev. E, 74(4), p. 046709. [CrossRef]
Lee, T. , 2009, “ Effects of Incompressibility on the Elimination of Parasitic Currents in the Lattice Boltzmann Equation Method for Binary Fluids,” Comput. Math. Appl., 58(5), pp. 987–994. [CrossRef]
Lee, T. , and Liu, L. , 2010, “ Lattice Boltzmann Simulations of Micron-Scale Drop Impact on Dry Surfaces,” J. Comput. Phys., 229(20), pp. 8045–8063. [CrossRef]
Li, L. , Jia, X. , Liu, Y. , and Su, M. , 2016, “ Simulation of Double Droplets Impact on Liquid Film by a Simplified Lattice Boltzmann Model,” Appl. Therm. Eng., 98, pp. 656–669. [CrossRef]
Connington, K. , and Lee, T. , 2013, “ Lattice Boltzmann Simulations of Forced Wetting Transitions of Drops on Superhydrophobic Surfaces,” J. Comput. Phys., 250, pp. 601–615. [CrossRef]
Baroudi, L. , Kawaji, M. , and Lee, T. , 2014, “ Effects of Initial Conditions on the Simulation of Inertial Coalescence of Two Drops,” Comput. Math. Appl., 67(2), pp. 282–289. [CrossRef]
Chen, S. , Martinez, D. , and Mei, R. , 1996, “ On Boundary Conditions in Lattice Boltzmann Methods,” Phys. Fluids, 8(9), pp. 2527–2536. [CrossRef]
Guo, Z. , Zheng, C. , and Shi, B. , 2002, “ An Extrapolation Method for Boundary Conditions in Lattice Boltzmann Method,” Phys. Fluids, 14(6), pp. 2007–2010. [CrossRef]
Junk, M. , and Yang, Z. , 2005, “ Asymptotic Analysis of Lattice Boltzmann Boundary Conditions,” J. Stat. Phys., 121(1–2), pp. 3–35. [CrossRef]
Yang, Z. , 2013, “ Lattice Boltzmann Outflow Treatments: Convective Conditions and Others,” Comput. Math. Appl., 65(2), pp. 160–171. [CrossRef]
Lou, Q. , Guo, Z. , and Shi, B. , 2013, “ Evaluation of Outflow Boundary Conditions for Two-Phase Lattice Boltzmann Equation,” Phys. Rev. E, 87(6), p. 063301. [CrossRef]


Grahic Jump Location
Fig. 1

Schematic illustration of the nodes around the outlet boundary

Grahic Jump Location
Fig. 2

Schematic illustration of the flow process

Grahic Jump Location
Fig. 3

Droplet shape and velocity fields for different OBC schemes at T∗=5.0,5.5, and 5.75. The contour lines from the exterior to the interior represent ρ=0.1,0.5, and 0.9. Re=24, We=133, and ρr=100.

Grahic Jump Location
Fig. 6

Droplet shape and velocity fields for three MCBC schemes at T∗=26.25,  27.5, and 28.75. The contour lines fromexterior to interior represent ρ=0.1,0.5,and 0.9. Re=72, We=133, and ρr=1000.

Grahic Jump Location
Fig. 5

Instantaneous relative error curves for the modified OBC schemes. Re=24, We=133, and ρr=100.

Grahic Jump Location
Fig. 4

Outlet velocity profiles along the height direction for the modified OBC schemes. Re=24, We=133, and ρr=100.

Grahic Jump Location
Fig. 8

Shape of the droplet for different Weber numbers at T∗=20. (a) We=133, (b) We=400, and (c) We=1333.

Grahic Jump Location
Fig. 7

Instantaneous relative error curves for the three MCBC schemes. Re=72, We=133, and ρr=1000.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In