中国科学技术大学学报 ›› 2018, Vol. 48 ›› Issue (5): 385-391.DOI: 10.3969/j.issn.0253-2778.2018.05.007

• 论著 • 上一篇    下一篇

黏弹性流体中柔性丝线拍动的数值模拟研究

卢杨,阙夏,刘难生   

  1. 1.中国科学技术大学工程科学学院近代力学系,安徽合肥 230027;2.合肥工业大学计算机与信息学院,安徽合肥 230009
  • 收稿日期:2017-02-14 修回日期:2017-05-30 接受日期:2017-05-30 出版日期:2018-05-31 发布日期:2017-05-30
  • 通讯作者: 刘难生
  • 作者简介:卢杨,女,1992年生,硕士.研究方向:黏弹性流体流固耦合数值模拟.E-mail: luyang18@mail.ustc.edu.cn
  • 基金资助:
    国家自然科学基金(11472268, 11402261)资助

Numerical investigation on flapping of a flexible filament in a viscoelastic fluid flow

LU Yang, QUE Xia, LIU Nansheng   

  1. 1. Department of Modern Mechanics, School Engineering Science, University of Science and Technology of China, Hefei 230027, China; 2. School of Computer and Information Science, Hefei University of Technology, Hefei 230009, China
  • Received:2017-02-14 Revised:2017-05-30 Accepted:2017-05-30 Online:2018-05-31 Published:2017-05-30

摘要: 采用格子波尔兹曼方法结合有限差分方法求解黏弹性流体的运动方程,非线性有限元方法求解柔性丝线的运动方程,浸没边界法处理流固耦合,数值求解了柔性丝线在均匀来流中因流固耦合作用产生的自由拍动问题.通过与牛顿流情况的对比,重点考察了流体弹性效应对丝线拍动特性的影响.计算结果发现,对于流体弹性效应较弱的情况(We<20),柔性丝线从稳定(静止)模态过渡为周期性拍动模态的临界质量比随着We数的增加显著增大;而对于弹性效应较强的情况(We>20),临界质量比随着We 数的增加渐近地趋于一常值.另外还发现,当质量比给定时,拍动丝线的阻力系数、拍动幅值和拍动频率均随We 的增加而减小.以上发现表明,流体弹性效应的增强对丝线拍动和尾迹流动失稳脱涡具有明显的抑制作用.

关键词: 非牛顿流体, 流固耦合, 柔性丝线

Abstract: The flapping of a flexible filament in a uniform incoming flow of viscoelastic fluid of polymeric solution was numerically investigated. This paper aims at examining the effects of fluid elasticity on the filament flapping behavior by a comparison with its counterpart in the Newtonian flow. Specifically, the FENE-MCR model was used as the constitutive equation for the viscoelastic fluid. The simulations of fluid flow were based on the lattice Boltzmann method to solve the Navier-Stokes (N-S) equations, the finite difference method to solve the constitutive equation of the polymer stress, the finite element method to solve the motion equation of the filament, and a penalty immersed boundary method to deal with the fluid-structure interaction. From the numerical results obtained, it was found that for a relatively weak fluid elasticity of We<20, the critical mass ratio of the filament to the fluid increases dramatically for the transition from a steady state to a periodically flapping state. However, for a stronger elasticity of We>20, increasing We number leads to a tendency of the critical mass ratio to be an approximate value of ~0.43. Moreover, it was demonstrated that increasing We number is also commensurate with the decrease of the time-averaged drag coefficient, the flapping amplitude, and the flapping frequency of the flexible filament. The above-mentioned facts indicate that enhancing the fluid elasticity has the increasing suppressing effects on the filament flapping behaviors.

Key words: viscoelastic flow, fluid-structure interaction, flexible filament

中图分类号: