-
针对流固耦合问题,发展了基于浸入边界-多松弛时间格子玻尔兹曼通量求解法(immersed boundary method multi-relaxation-time lattice Boltzmann flux solver,IB-MRT-LBFS)的弱耦合算法.依据多尺度Chapman-Enskog展开,建立不可压宏观方程状态变量和通量与格子玻尔兹曼方程中粒子密度分布函数之间的关系;采用强制浸入边界法处理流固界面使固壁表面满足无滑移边界条件,根据修正的速度求解动量方程力源项;结构运动方程采用四阶龙格-库塔法求解.格子模型与浸入边界法的引入使流固耦合计算可以在笛卡尔网格下进行,无需生成贴体网格及运用动网格技术,简化了计算过程.数值模拟了单圆柱横向涡激振动、单圆柱及串列双圆柱双自由度涡激振动问题.结果表明,IB-MRT-LBFS能够准确预测圆柱涡激振动的锁定区间、振动响应、受力情况以及捕捉尾流场结构形态,验证了该算法在求解流固耦合问题的有效性和可行性.
-
关键词:
- 流固耦合 /
- 浸入边界法 /
- 格子玻尔兹曼通量求解法 /
- 涡激振动
This paper performs a newly developed method, which combines the immersed boundary method (IBM) with multi-relaxation-time lattice Boltzmann flux solver (MRT-LBFS), for solving fluid-structure interaction problems. Finite volume discretization is used to solve the macroscopic governing equations with the flow variables defined at cell centers. Based on the multi-scale Chapman-Enskog expansion analysis, LBFS builds a relationship between the variables and fluxes in incompressible Navier-Stokes equations and density distribution functions in lattice Boltzmann equation. In order to ensure no-slip boundary condition, boundary condition-enforced immersed boundary method is used to treat the fluid-structure interface. The restoring force can be resolved by making a velocity correction in the flow field. The four-stage RungeKutta scheme is used to solve the motion equation of structure. Using the lattice model and immersed boundary method, fluid-structure coupling calculation can be implemented in a Cartesian grid, without generating the body-fitted mesh and using moving mesh technique. Therefore, the computational process is considerably simplified. In order to verify the validity and feasibility of IB-MRT-LBFS to solve fluid-structure interaction problems, both one-and two-degree of freedom vortex-induced vibrations (VIV) of a circular cylinder and two-degree of freedom VIV of two cylinders in a tandem arrangement are simulated by this proposal method. For a VIV cylinder system, the transverse vibration response is much stronger than the axial response. When the vibration occurs in the range of lock-in regime, the shedding vortex frequency of the wake is close to natural frequency of the cylinder so that resonance appears, consequently causing larger amplitude. For two VIV cylinders in a tandem arrangement, the dynamic behavior of each cylinder is significantly different from that of a single cylinder. The gap spacing between the two cylinder centers is a significant parameter which effects vibration characteristics and the spacing is fixed in the simulations of two tandem cylinders. With the effects of upstream cylinder wake, the axial and transverse amplitudes of downstream cylinder obviously increase with adding the reduced velocity. The downstream cylinder is delayed, coming into lock-in regime, and the range of lock-in regime is expanded under the effects of the wake of the upstream cylinder. As the reduced velocity is relatively large, the vibration response of the upstream cylinder is close to a single cylinder and the vibration response of the downstream cylinder is more intense than the upstream cylinder. Compared with the existing literature results, our result illustrates that IB-MRT-LBFS owns the ability to correctly predict the lock-in regime, dynamic response and the forces of vortex-induced vibrations of cylinders. And this method can accurately capture the wake structures.-
Keywords:
- fluid-structure interaction /
- immersed boundary method /
- lattice Boltzmann flux solver /
- vortex-induced vibration
[1] Xing J T, Zhou S, Cui E J 1997 Adv. Mech. 27 19 (in Chinese) [邢景棠, 周盛, 崔尔杰 1997 力学进展 27 19]
[2] Qian R J, Dong S L, Yuan X F (in Chinese) [钱若军, 董石麟, 袁行飞 2008 空间结构 14 3]
[3] Guo P, Liu J, Wu W H 2013 Chin. J. Theor. Appl. Mech. 45 283 (in Chinese) [郭攀, 刘君, 武文华 2013 力学学报 45 283]
[4] Zhou D, He T, Tu J H (in Chinese) [周岱, 何涛, 涂佳黄 2012 力学学报 44 494]
[5] Zhong G H, Liang A, Sun X F 2007 J. Eng. Thermophys. 28 399 (in Chinese) [钟国华, 梁岸, 孙晓峰 2007 工程热 28 399]
[6] Liu Q Y 2012 M. S. Dissertation (Nanjing: Nanjing University of Aeronautics and Astronautics) (in Chinese) [刘齐迎 2012 硕士学位论文(南京: 南京航空航天大学]
[7] Luo H X, Dai H, Ferreira D S, Paulo J S A, Yin B 2012 Comput. Fluids 56 61
[8] Feng Z, Michaelides E 2004 J. Comput. Phys. 195 602
[9] Chen Y, Cai Q D, Xia Z H, Wang M, Chen S Y 2013 Phys. Rev. E 88 013303
[10] Wang W Q, Zhang G W, Yan Y (in Chinese) [王文全, 张国威, 闫妍 2017 北京理工大学学报 37 151]
[11] Wang W Q, Su S Q, Yan Y (in Chinese) [王文全, 苏仕琪, 闫妍 2015 计算力学学报 32 560]
[12] Ming P J, Zhang W P 2009 Chin. J. Aeronaut. 22 480
[13] Ming P J, Zhang W P, Lu X Q, Zhu M G (in Chinese) [明平剑, 张文平, 卢熙群, 朱明刚 2010 水动力研究与进展 25 321]
[14] Li S Y, Cheng Y G, Zhang C Z 2016 J. Huazhong Univ. Sci. Tech. (Natural Science Edition) 44 122 (in Chinese) [李师尧, 程永光, 张春泽 2016 华中科技大学学报 (自然科学版) 44 122]
[15] Shu C, Wang Y, Teo C J, Wu J 2014 Adv. Appl. Math. Mech. 6 436
[16] Wang Y, Shu C, Teo C J, Wu J 2015 J. Fluids Struct. 54 440
[17] Wang Y, Shu C, Yang L M, Sun Y 2017 Int. J. Numer. Meth. Fluids 83 331
[18] Suzuki K, Inamuro T 2011 Comput. Fluids 49 173
[19] Ahn H T, Kallinderis Y 2006 J. Comput. Phys. 219 671
[20] Borazjani I, Ge L, Sotiropoulos F 2008 J. Comput. Phys. 227 7587
[21] Jiang R J, Lin J Z, Chen Z L 2013 Phys. Rev. E 88 023009
[22] Wang C L, Tang H, Duan F, Yu S C M 2016 J. Fluids Struct. 60 160
[23] Han Z L, Zhou D, Tu J H 2014 J. Eng. Mech 140 04014059
[24] Prasanth A K, Mittal S 2008 J. Fluids Mech. 594 463
[25] Bao Y, Huang C, Zhou D, Tu J H, Han Z L 2012 J.Fluids Struct. 35 50
[26] Yu K R, Etienne S Scolan, Y M, Hay A, Fontaine E, Pelletier D 2016 J. Fluids Struct. 60 37
-
[1] Xing J T, Zhou S, Cui E J 1997 Adv. Mech. 27 19 (in Chinese) [邢景棠, 周盛, 崔尔杰 1997 力学进展 27 19]
[2] Qian R J, Dong S L, Yuan X F (in Chinese) [钱若军, 董石麟, 袁行飞 2008 空间结构 14 3]
[3] Guo P, Liu J, Wu W H 2013 Chin. J. Theor. Appl. Mech. 45 283 (in Chinese) [郭攀, 刘君, 武文华 2013 力学学报 45 283]
[4] Zhou D, He T, Tu J H (in Chinese) [周岱, 何涛, 涂佳黄 2012 力学学报 44 494]
[5] Zhong G H, Liang A, Sun X F 2007 J. Eng. Thermophys. 28 399 (in Chinese) [钟国华, 梁岸, 孙晓峰 2007 工程热 28 399]
[6] Liu Q Y 2012 M. S. Dissertation (Nanjing: Nanjing University of Aeronautics and Astronautics) (in Chinese) [刘齐迎 2012 硕士学位论文(南京: 南京航空航天大学]
[7] Luo H X, Dai H, Ferreira D S, Paulo J S A, Yin B 2012 Comput. Fluids 56 61
[8] Feng Z, Michaelides E 2004 J. Comput. Phys. 195 602
[9] Chen Y, Cai Q D, Xia Z H, Wang M, Chen S Y 2013 Phys. Rev. E 88 013303
[10] Wang W Q, Zhang G W, Yan Y (in Chinese) [王文全, 张国威, 闫妍 2017 北京理工大学学报 37 151]
[11] Wang W Q, Su S Q, Yan Y (in Chinese) [王文全, 苏仕琪, 闫妍 2015 计算力学学报 32 560]
[12] Ming P J, Zhang W P 2009 Chin. J. Aeronaut. 22 480
[13] Ming P J, Zhang W P, Lu X Q, Zhu M G (in Chinese) [明平剑, 张文平, 卢熙群, 朱明刚 2010 水动力研究与进展 25 321]
[14] Li S Y, Cheng Y G, Zhang C Z 2016 J. Huazhong Univ. Sci. Tech. (Natural Science Edition) 44 122 (in Chinese) [李师尧, 程永光, 张春泽 2016 华中科技大学学报 (自然科学版) 44 122]
[15] Shu C, Wang Y, Teo C J, Wu J 2014 Adv. Appl. Math. Mech. 6 436
[16] Wang Y, Shu C, Teo C J, Wu J 2015 J. Fluids Struct. 54 440
[17] Wang Y, Shu C, Yang L M, Sun Y 2017 Int. J. Numer. Meth. Fluids 83 331
[18] Suzuki K, Inamuro T 2011 Comput. Fluids 49 173
[19] Ahn H T, Kallinderis Y 2006 J. Comput. Phys. 219 671
[20] Borazjani I, Ge L, Sotiropoulos F 2008 J. Comput. Phys. 227 7587
[21] Jiang R J, Lin J Z, Chen Z L 2013 Phys. Rev. E 88 023009
[22] Wang C L, Tang H, Duan F, Yu S C M 2016 J. Fluids Struct. 60 160
[23] Han Z L, Zhou D, Tu J H 2014 J. Eng. Mech 140 04014059
[24] Prasanth A K, Mittal S 2008 J. Fluids Mech. 594 463
[25] Bao Y, Huang C, Zhou D, Tu J H, Han Z L 2012 J.Fluids Struct. 35 50
[26] Yu K R, Etienne S Scolan, Y M, Hay A, Fontaine E, Pelletier D 2016 J. Fluids Struct. 60 37
计量
- 文章访问数: 7716
- PDF下载量: 390
- 被引次数: 0