全部

力学与实践 ›› 2014, Vol. 36 ›› Issue (4): 383-395.doi: 10.6052/1000-0879-14-168

• 专题综述 • 上一篇    下一篇

非牛顿流体流动的格子Boltzmann方法研究进展

李勇, 柳文琴   

  1. 华东理工大学机械与动力工程学院, 上海200237
  • 收稿日期:2014-04-28 修回日期:2014-05-19 出版日期:2014-08-08 发布日期:2014-11-13
  • 通讯作者: 李勇,博士,副教授.E-mail:yong_li@ecust.edu.cn E-mail:yong_li@ecust.edu.cn
  • 基金资助:

    国家自然科学基金资助项目(11372111).

THE RESEARCH PROGRESS OF LATTICE BOLTZMANN METHOD IN NON-NEWTONIAN FLUID FLOW

LI Yong, LIU Wenqin   

  1. School of Mechanical and Power Engineering, East China University of Science and Technology, Shanghai 200237, China
  • Received:2014-04-28 Revised:2014-05-19 Online:2014-08-08 Published:2014-11-13

摘要:

格子玻尔兹曼方法(lattice Boltzmann method,LBM)能够直接计算局部剪切速率并可以达到二次精度,因此在非牛顿流动数值模拟中展现出一定优势。尽管已证实LBM 对于非牛顿流动的适用性,但是LBM 需要通过即时调节BGK(Bhatnagar-Gross-Krook)碰撞项中的松弛时间来实时反映黏度改变,当松弛时间接近1/2 时,迭代会出现数值不稳定现象。该文对LBM 在非牛顿流体研究中的进展进行了总结,介绍了增加数值稳定性的方法并对结果的精度进行了比较,在此基础上对LBM 在非牛顿研究中的进一步发展进行了展望。

关键词:

格子玻尔兹曼方法|广义牛顿流体|黏弹性流体|分布函数

Abstract:

LBM (lattice Boltzmann method) allows straightforward calculations of the local shear rate to the second-order accuracy, so it has some advantages in the simulation of the non-Newtonian behavior of flows. Although the suitability of LBM for the non-Newtonian flow has been well demonstrated, the algorithms are mainly through instantaneously regulating the local relaxation time in the BGK (Bhatnagar-Gross-Krook) approximated collision term to match its corresponding viscosity for reflecting the effect of the local rheology variation.This approach, however, might cause some numerical instabilities when the relation time is close to 1/2.In this paper, the advances of LBM applied to generalized Newtonian flows are reviewed.The methods to mitigate numerical instabilities are described and the results are compared.Several ways of extending LBM to viscoelastic flows in the literature are introduced.Finally the prospective developments of LBM for non-Newtonian flows are suggested.

Key words:

lattice Boltzmann method|generalized Newtonian flow|viscoelastic fluid|distribution function

中图分类号: 

  • O373