基于直角网格法的单个和阵列布置下柔性水翼绕流数值模拟

辛建建; 陈振雷; 石凡; 石伏龙

doi:10.7498/aps.69.20191711

摘要

研究柔性水翼在不可压缩流体中的水动力特性, 对于船舵和减摇鳍等海洋结构物的设计和性能优化具有重要意义. 本文将自主开发的径向基函数虚拟网格法求解器扩展到模拟绕单个或多个柔性水翼的不可压缩流动问题. 数值模型基于虚拟网格有限差分法考虑浸入边界对流场的影响, 引入紧支径向基函数(compact supported radial basis function, CSRBF)以物面Lagrangian质点追踪复杂的柔性动边界. 基于该方法, 首先模拟了均匀流中主动拍动的柔性水翼, 升阻力系数良好的网格收敛性结果验证了本文方法的精度和可靠性. 并研究了柔性水翼在不同振荡频率下的水动力特性, 阐述了柔性水翼的推力生成机制. 然后模拟了绕阵列布置柔性水翼的流动现象, 研究了不同间距和不同振荡频率下水翼表面的升阻力系数变化规律和尾涡特性, 观察到紧密布置的柔性水翼在高频振荡下推力系数存在显著的放大效应, 同时推力为零时的临界频率提前.

关键词:

Abstract

Studies on hydrodynamic characteristics of viscous incompressible flows around flexible hydrofoils are of practical importance for the design and performance optimization of marine structures such as ship rudders and stabilizing fins. The aim of this paper is to extend a radial basis function based ghost cell method to simulate flows around single or multiple flexible moving hydrofoils in array arrangement. The numerical model is based on a ghost cell finite difference method for considering the influence of the immersed boundaries on the flow. Also, a compact supported radial basis function (CSRBF) is introduced to track the complex flexible boundary with some controlling points of the body surface. Based on the present method, the uniform flow around a flexible hydrofoil swimming like a fish is simulated. Good grid convergence of drag and lift coefficients demonstrates the accuracy and reliability of the present method. Also, the hydrodynamics patterns of the flexible hydrofoil under different oscillation frequencies are studied. Further, the thrust generation mechanism of the hydrofoil is explained. Afterwards, flows around the undulating hydrofoils in array arrangement are simulated. The force coefficients and wake patterns under different distances and oscillation frequencies are investigated. It is observed that the thrust coefficients of the hydrofoils under narrow arrangement and high oscillation frequencies have significant amplification effects. In addition, the critical frequency at zero thrust reduces.

Keywords:

PACS:

作者及机构信息

1.
宁波大学海运学院, 宁波　315211

2.
长沙理工大学水利工程学院, 长沙　410114

通信作者: 陈振雷, chenzhenlei@nbu.edu.cn

基金项目: 国家级-国家青年科学基金项目(51909124)

Authors and contacts

1.
Faculty of Maritime and Transportation, Ningbo University, Ningbo 315211, China

2.
School of Hydraulic Engineering, Changsha University of Science and Technology, Changsha 410114, China

Corresponding author: Chen Zhen-Lei, chenzhenlei@nbu.edu.cn

文章全文

参考文献

[1]	Zabihi M, Lari K, Amiri H 2017 J. Mech. Sci. Technol. 31 3539 Google Scholar
[2]	Yang J M 2016 J. Hydrodyn. Ser. B 28 713 Google Scholar
[3]	王力, 田方宝 2018 中国科学: 物理学力学天文学 48 094703 Google Scholar Wang L, Tian F B 2018 Sci. China Phys. Mech. Astron. 48 094703 Google Scholar
[4]	Al-Marouf M, Samtaney R 2017 J. Comput. Phys. 337 339 Google Scholar
[5]	Huang W X, Chang C B, Sung H J 2011 J. Comput. Phys. 230 5061 Google Scholar
[6]	Tullio M D D, Pascazio G 2016 J. Comput. Phys. 325 201 Google Scholar
[7]	吴晓笛, 刘华坪, 陈浮 2017 22 224702 Google Scholar Wu X D, Liu H P, Chen F 2017 Acta Phys. Sin. 22 224702 Google Scholar
[8]	Yeo K S, Ang S J, Shu C 2010 Comput. Fluids 39 403 Google Scholar
[9]	Tian F B, Wang W, Wu J, Sui Y 2016 Comput. Fluids 124 1 Google Scholar
[10]	Bergmann M, Iollo A, Mittal R 2014 Bioinspiration Biomimetics 9 046001 Google Scholar
[11]	Khalid M S U, Akhtar I, Dong H 2016 J. Fluids Struct. 66 19 Google Scholar
[12]	Khalid M S U, Akhtar I, Imtiaz H, Dong H, Wu B 2018 Ocean Eng. 157 108 Google Scholar
[13]	Mittal R, Dong H, Bozkurttas M, Najjar F M, Vargas A, Loebbecke A V 2008 J. Comput. Phys. 227 4825 Google Scholar
[14]	王亮, 吴锤结 2011 力学学报 43 18 Google Scholar Wang L, Wu C J 2011 Chin. J. Theor. Appl. Mech. 43 18 Google Scholar
[15]	王文全 2014 计算力学学报 31 646 Google Scholar Wang W Q 2014 Chin. J. Comput. Mech. 31 646 Google Scholar
[16]	辛建建, 石伏龙, 金秋 2017 66 186 Xin J J, Shi F L, Jin Q 2017 Acta Phys. Sin. 66 186
[17]	Xin J J, Shi F L, Jin Q, Lin C 2018 Comput. Fluids 176 210 Google Scholar
[18]	Leer B V 1979 J. Comput. Phys. 32 101 Google Scholar
[19]	Van d V H A, Kohn R V, Iserles A, Ciarlet P G, Wright M H 2006 Iterative Krylov Methods for Large Linear Systems (Cambridge: Cambridge University Press) pp133–135
[20]	Liu G R, Gu Y T 2005 An Introduction to Meshfree Methods and Their Programming (Netherlands: Springer) p73
[21]	Sui Y, Chew Y T, Roy P, Low H T 2007 Int. J. Numer. Methods Fluids 53 1727 Google Scholar

图 1 虚拟网格法的二维插值格式　(a)无虚拟网格; (b)一个虚拟网格; (c)两个虚拟网格

Fig. 1. Two-dimensional interpolation stencil of the ghost cell method: (a) No ghost cells; (b) one ghost cell; (c) two ghost cells.

下载: 全尺寸图片幻灯片

图 2 静止或运动边界流动求解算法的流程图

Fig. 2. Flowchart of numerical algorithm for stationary or moving boundary flows.

下载: 全尺寸图片幻灯片

图 3 五个频率下的力系数比较　(a) 平均阻力系数; (b) 均方根升力系数

Fig. 3. Comparison of force coefficients under five oscillation frequencies: (a) Average drag coefficient; (b) root mean square lift coefficient.

下载: 全尺寸图片幻灯片

图 4 三个频率下升阻力系数随无因次时间的变化趋势　(a) 阻力系数; (b) 升力系数

Fig. 4. Variation of drag and lift coefficients with the dimensionless time under three oscillation frequencies: (a) Drag coefficient; (b) lift coefficient.

下载: 全尺寸图片幻灯片

图 5 $\omega$ = 0.15时在两个时刻的局部瞬时涡等值线　(a) $\theta =$ ${\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$

Fig. 5. Local instantaneous vortex contours for $\omega$ = 0.15 at two moments: (a) $\theta = {\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$ .

下载: 全尺寸图片幻灯片

图 7 $\omega$ = 0.75时在两个时刻的局部瞬时涡等值线　(a) $\theta =$ ${\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$

Fig. 7. Local instantaneous vortex contours for $\omega$ = 0.75 at two moments: (a) $\theta = {\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$ .

下载: 全尺寸图片幻灯片

图 6 $\omega$ = 0.45时在两个时刻的局部瞬时涡等值线　(a) $\theta =$ ${\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$

Fig. 6. Local instantaneous vortex contours for $\omega$ = 0.45 at two moments: (a) $\theta = {\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$ .

下载: 全尺寸图片幻灯片

图 8 阵列布置水翼游动的几何模型

Fig. 8. Geometrical model of the undulating hydrofoil in array arrangement.

下载: 全尺寸图片幻灯片

图 9 三个间距下群游中间位置水翼和单水翼在五个频率下的力系数比较　(a)平均阻力系数; (b)均方根升力系数

Fig. 9. Comparison of force coefficients between the central hydrofoil in array arrangement for three gap distances and the single hydrofoil under five frequencies: (a) Average drag coefficient; (b) root mean square lift coefficient.

下载: 全尺寸图片幻灯片

图 10 群游中间位置水翼在三个频率下升阻力系数随无因次时间的变化趋势　(a)阻力系数; (b)升力系数

Fig. 10. Variation of drag and lift coefficients on the central hydrofoil in array arrangement with the dimensionless time under three oscillation frequencies: (a) Drag coefficient; (b) lift coefficient.

下载: 全尺寸图片幻灯片

图 11 $\omega$ = 0.15时在两时刻的水翼群游局部涡场图　 (a) $\theta = {\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$

Fig. 11. Local instantaneous vortex contours of undulating hydrofoils in array arrangement for $\omega$ = 0.15 at two moments: (a) $\theta = {\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$ .

下载: 全尺寸图片幻灯片

图 12 $\omega$ = 0.45时在两时刻的水翼群游局部涡场图　(a) $\theta =$ ${\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$

Fig. 12. Local instantaneous vortex contours of undulating hydrofoils in array arrangement for $\omega$ = 0.45 at two moments: (a) $\theta = {\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$ .

下载: 全尺寸图片幻灯片

图 13 $\omega$ = 0.75时在两时刻的水翼群游局部涡场图　(a) $\theta =$ ${\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$

Fig. 13. Local instantaneous vortex contours of undulating hydrofoils in array arrangement for $\omega$ = 0.75 at two moments: (a) $\theta = {\text{π}}/2$ ; (b) $\theta = 3{\text{π}}/2$ .

下载: 全尺寸图片幻灯片

表 1 本文力系数结果与Sui等^[21]的参考结果比较

Table 1. Comparison of the force coefficients between the present method and Sui’s method^[21].

方法平均阻力系数均方根升力系数

150 × 80 0.200 0.91
300 × 150 0.184 1.15
600 × 300 0.182 1.16
Sui等^[21] 0.180 1.20

下载: 导出CSV

map

[1]	Zabihi M, Lari K, Amiri H 2017 J. Mech. Sci. Technol. 31 3539 Google Scholar
[2]	Yang J M 2016 J. Hydrodyn. Ser. B 28 713 Google Scholar
[3]	王力, 田方宝 2018 中国科学: 物理学力学天文学 48 094703 Google Scholar Wang L, Tian F B 2018 Sci. China Phys. Mech. Astron. 48 094703 Google Scholar
[4]	Al-Marouf M, Samtaney R 2017 J. Comput. Phys. 337 339 Google Scholar
[5]	Huang W X, Chang C B, Sung H J 2011 J. Comput. Phys. 230 5061 Google Scholar
[6]	Tullio M D D, Pascazio G 2016 J. Comput. Phys. 325 201 Google Scholar
[7]	吴晓笛, 刘华坪, 陈浮 2017 22 224702 Google Scholar Wu X D, Liu H P, Chen F 2017 Acta Phys. Sin. 22 224702 Google Scholar
[8]	Yeo K S, Ang S J, Shu C 2010 Comput. Fluids 39 403 Google Scholar
[9]	Tian F B, Wang W, Wu J, Sui Y 2016 Comput. Fluids 124 1 Google Scholar
[10]	Bergmann M, Iollo A, Mittal R 2014 Bioinspiration Biomimetics 9 046001 Google Scholar
[11]	Khalid M S U, Akhtar I, Dong H 2016 J. Fluids Struct. 66 19 Google Scholar
[12]	Khalid M S U, Akhtar I, Imtiaz H, Dong H, Wu B 2018 Ocean Eng. 157 108 Google Scholar
[13]	Mittal R, Dong H, Bozkurttas M, Najjar F M, Vargas A, Loebbecke A V 2008 J. Comput. Phys. 227 4825 Google Scholar
[14]	王亮, 吴锤结 2011 力学学报 43 18 Google Scholar Wang L, Wu C J 2011 Chin. J. Theor. Appl. Mech. 43 18 Google Scholar
[15]	王文全 2014 计算力学学报 31 646 Google Scholar Wang W Q 2014 Chin. J. Comput. Mech. 31 646 Google Scholar
[16]	辛建建, 石伏龙, 金秋 2017 66 186 Xin J J, Shi F L, Jin Q 2017 Acta Phys. Sin. 66 186
[17]	Xin J J, Shi F L, Jin Q, Lin C 2018 Comput. Fluids 176 210 Google Scholar
[18]	Leer B V 1979 J. Comput. Phys. 32 101 Google Scholar
[19]	Van d V H A, Kohn R V, Iserles A, Ciarlet P G, Wright M H 2006 Iterative Krylov Methods for Large Linear Systems (Cambridge: Cambridge University Press) pp133–135
[20]	Liu G R, Gu Y T 2005 An Introduction to Meshfree Methods and Their Programming (Netherlands: Springer) p73
[21]	Sui Y, Chew Y T, Roy P, Low H T 2007 Int. J. Numer. Methods Fluids 53 1727 Google Scholar

[1]	巫梦丹, 周胜林, 叶安娜, 王敏, 张晓华, 杨朝晖. 基于中性水凝胶/取向碳纳米管阵列高电压柔性固态超级电容器. , 2019, 68(10): 108201. doi: 10.7498/aps.68.20182288
[2]	张妮, 刘丁, 冯雪亮. 直拉硅单晶生长过程中工艺参数对相变界面形态的影响. , 2018, 67(21): 218701. doi: 10.7498/aps.67.20180305
[3]	吴晓笛, 刘华坪, 陈浮. 基于浸入边界-多松弛时间格子玻尔兹曼通量求解法的流固耦合算法研究. , 2017, 66(22): 224702. doi: 10.7498/aps.66.224702
[4]	辛建建, 石伏龙, 金秋. 一种径向基函数虚拟网格法数值模拟复杂边界流动. , 2017, 66(4): 044704. doi: 10.7498/aps.66.044704
[5]	李强, 李五明. 带嵌件型腔内熔接过程的数值模拟研究. , 2016, 65(6): 064601. doi: 10.7498/aps.65.064601
[6]	魏德志, 陈福集, 郑小雪. 基于混沌理论和改进径向基函数神经网络的网络舆情预测方法. , 2015, 64(11): 110503. doi: 10.7498/aps.64.110503
[7]	杨青, 曹曙阳, 刘十一. 基于浸入式边界方法的串联双矩形柱绕流数值模拟. , 2014, 63(21): 214702. doi: 10.7498/aps.63.214702
[8]	于涛, 尹成友, 刘汉. 基于特征基函数的球面共形微带天线阵列分析. , 2014, 63(23): 230701. doi: 10.7498/aps.63.230701
[9]	刘飞飞, 魏守水, 魏长智, 任晓飞. 基于速度源修正的浸入边界-晶格玻尔兹曼法研究仿生微流体驱动模型. , 2014, 63(19): 194704. doi: 10.7498/aps.63.194704
[10]	曾喆昭. 不确定混沌系统的径向基函数神经网络反馈补偿控制. , 2013, 62(3): 030504. doi: 10.7498/aps.62.030504
[11]	陈帝伊, 柳烨, 马孝义. 基于径向基函数神经网络的混沌时间序列相空间重构双参数联合估计. , 2012, 61(10): 100501. doi: 10.7498/aps.61.100501
[12]	刘佳, 徐玲玲, 张海霖, 吕威, 朱琳, 高红, 张喜田. 一步水热法在Al掺杂ZnO纳米盘上可控自组装合成ZnO纳米棒阵列. , 2012, 61(2): 027802. doi: 10.7498/aps.61.027802
[13]	武浩, 朱拓, 孔艳, 陈卫, 杨建磊. 基于径向基函数神经网络的荧光光谱技术在菌种识别中的应用. , 2010, 59(4): 2396-2400. doi: 10.7498/aps.59.2396
[14]	白安琪, 胡迪, 丁武昌, 苏少坚, 胡炜玄, 薛春来, 樊中朝, 成步文, 俞育德, 王启明. 用PAA模板法实现硅基纳米孔阵列结构. , 2009, 58(7): 4997-5001. doi: 10.7498/aps.58.4997
[15]	戴保东, 程玉民. 势问题的径向基函数——局部边界积分方程方法. , 2007, 56(2): 597-603. doi: 10.7498/aps.56.597
[16]	司马文霞, 刘凡, 孙才新, 廖瑞金, 杨庆. 基于改进的径向基函数神经网络的铁磁谐振系统混沌控制. , 2006, 55(11): 5714-5720. doi: 10.7498/aps.55.5714
[17]	李宏斌, 刘文清, 张玉钧, 丁志群, 赵南京, 魏庆农, 王玉平, 杨立书. 基于径向基函数网络的激光诱导荧光特征光谱分离算法. , 2005, 54(9): 4451-4457. doi: 10.7498/aps.54.4451
[18]	郭会军, 刘君华. 基于径向基函数神经网络的Lorenz混沌系统滑模控制. , 2004, 53(12): 4080-4086. doi: 10.7498/aps.53.4080
[19]	刘丁, 任海鹏, 孔志强. 基于径向基函数神经网络的未知模型混沌系统控制. , 2003, 52(3): 531-535. doi: 10.7498/aps.52.531
[20]	申金媛, 刘, 常胜江, 贾佳, 张文伟, 张延, 母国光. 基于径向基函数的多目标旋转不变分类及其光电实现. , 1998, 47(12): 1968-1975. doi: 10.7498/aps.47.1968

计量

文章访问数: 7581
PDF下载量: 67

姓名
邮箱
手机号码
标题
留言内容
验证码

搜索

留言板

基于直角网格法的单个和阵列布置下柔性水翼绕流数值模拟