复杂构型水滴收集率的拉格朗日快速算法研究

曾涛; 王昭力; 熊华杰; 陈宇; 周志宏

doi:10.13224/j.cnki.jasp.20220539

复杂构型水滴收集率的拉格朗日快速算法研究

doi: 10.13224/j.cnki.jasp.20220539

曾涛^{1, 2,},
王昭力¹,
熊华杰¹,
陈宇^{1, 2},
周志宏^{1, 2,*, ,}

1.
四川大学建筑与环境学院力学系，成都 610000
2.
宜宾产业技术研究院，四川宜宾 644000

基金项目: 国家自然科学基金（12072213）；国家科技重大专项（J2019-Ⅲ-0010-0054）；国家数值风洞工程（NNW2019-JT01-023）

详细信息

作者简介:
曾涛（1997-），男，硕士生，研究方向为空气动力学、飞机/风力机结冰

通讯作者:
周志宏（1981-），教授，博士，研究方向为低速空气动力学、飞机/风力机结冰。E-mail：zhouzhihong@scu.edu.cn

中图分类号: V211
计量
- 文章访问数: 399
- HTML浏览量: 238
- PDF量: 33
- 被引次数: 0
出版历程
- 收稿日期: 2022-07-24
- 网络出版日期: 2025-04-10

Research on a fast Lagrangian algorithm for water droplet collection efficiency in complex configurations

ZENG Tao^{1, 2
,},
WANG Zhaoli¹,
XIONG Huajie¹,
CHEN Yu^{1, 2},
ZHOU Zhihong^{1, 2,*
, ,}

1.
Department of Mechanics，College of Architecture and Environment，Sichuan University，Chengdu 610000，China
2.
Yibin Institute of Industrial Technology，Yibin Sichuan 644000，China

摘要

摘要:
针对飞机结冰数值模拟中拉格朗日方法计算复杂构型效率低、适应性差的问题，提出一种适用于复杂构型的水滴收集率的拉格朗日快速算法。基于蒙特卡洛粒子统计方法计算水滴收集率以提高算法对于复杂构型的适应性，采用分块追踪结合壁面距离的自适应分析方法确定适当覆盖所有撞击部件的初始水滴释放阵列，提出采用结合沿程水滴网格阵列加密技术和基于误差控制的变步长加速技术提高计算效率，分析并获得了计算效率与加速方法中各控制变量的关系。计算结果表明：所发展的算法能够精确、高效地获得三维水滴收集率，适用于复杂构型，为飞机水滴收集率计算提供了新方法和思路，可以为飞机结冰特性研究、防除冰系统设计等提供技术参。通过调控加密比例系数n（n>1）和加密释放面位置控制系数k（0<k<1），可在有效提升计算效率的同时将计算精度误差3%的范围内。
- 蒙特卡洛算法 /
- 水滴收集率 /
- 拉格朗日 /
- 误差控制 /
- 飞机结冰
Abstract:
In order to solve the problems of low efficiency and poor adaptability of Lagrange method in the numerical simulation of aircraft icing, a fast Lagrangian algorithm for water droplet collection efficiency in complex configurations was proposed. Monte Carlo particle release statistical method was used to calculate the water droplet collection efficiency, which can significantly improve the adaptability of the algorithm to complex configurations. By using the adaptive analysis method of grid block tracking and wall impact judgment, the range of initial water droplet release array covering all impacting components can be determined. A method combining the encryption technology of grid array of water droplets along the way and the variable step acceleration technology based on error control was proposed, which can improve the calculation efficiency. Based on this method, the relationship between calculation efficiency and control variables in acceleration method was analyzed and obtained. The calculation results showed that the developed algorithm can obtain the three-dimensional water droplet collection efficiency accurately and efficiently, and can be applied to complex configurations, providing a new method and idea for calculating the collection rate of aircraft water droplets. It can also provide technical parameters for the study of aircraft icing characteristics and the design of anti-icing system. By adjusting the encryption proportional coefficient n (n>1) and the encryption release surface position control coefficient k (0<k<1), it can effectively improve the calculation efficiency and reduce the calculation accuracy error within 3%.
- Monte Carlo algorithm /
- droplet collection efficiency /
- Lagrangian /
- error control /
- aircraft icing

HTML

图 1 三维水滴收集率计算示意图

Figure 1. Calculation of water droplet collection efficiency in three-dimensional configuration

下载: 全尺寸图片幻灯片

图 2 蒙特卡洛方法近似计算水滴收集率

Figure 2. Monte-Carlo approximate calculation of droplet collection efficiency

下载: 全尺寸图片幻灯片

图 3 水滴粒子初始释放单元自适应识别

Figure 3. Adaptive identification of initial release units of water droplet particles

下载: 全尺寸图片幻灯片

图 4 起始截面S_left、加密截面S_k与近壁截面S_right位置关系

Figure 4. Positional relationship between S_left of initial section, S_k of compact section and S_right of near-wall section

下载: 全尺寸图片幻灯片

图 5 加密计算后部分轨迹的分布效果图

Figure 5. Distribution effect diagram of partial trajectories after encrypted calculation

下载: 全尺寸图片幻灯片

图 6 三维圆柱网格模型

Figure 6. Grid model of three-dimensional cylinder

下载: 全尺寸图片幻灯片

图 7 三维圆柱表面水滴收集率云图

Figure 7. Cloud chart of water droplet collection efficiency on the surface of a 3D cylinder

下载: 全尺寸图片幻灯片

图 8 圆柱水滴收集率计算结果与试验结果对比（V=80 m/s, D_eq=16 μm, α=0°）

Figure 8. Comparison of calculation results and experimental results for cylindrical water droplet collection efficiency （V=80 m/s, D_eq=16 μm, α=0°）

下载: 全尺寸图片幻灯片

图 9 三维NACA0012机翼网格模型

Figure 9. Grid model of 3D NACA0012 airfoil

下载: 全尺寸图片幻灯片

图 10 三维NACA0012机翼水滴轨迹分布与收集率云图

Figure 10. Droplet trajectory distribution and counter of collection efficiency of the 3D NACA0012 airfoil

下载: 全尺寸图片幻灯片

图 11 三维NACA0012机翼水滴收集率与Lewice计算结果对比（V=102.8 m/s, D_eq=40 μm, α=4°）

Figure 11. Comparison of water droplet collection efficiency of 3D NACA0012 airfoil with Lewice’s calculation results （V=102.8 m/s, D_eq=40 μm, α=4°）

下载: 全尺寸图片幻灯片

图 12 计算效率与加密位置的关系

Figure 12. Relationship between computational efficiency and the location of encryption

下载: 全尺寸图片幻灯片

图 13 计算精度与加密位置的关系

Figure 13. Relationship between computational precision and the location of encryption

下载: 全尺寸图片幻灯片

图 14 计算精度和加密比例的关系

Figure 14. Relationship between computational precision and encryption ratio

下载: 全尺寸图片幻灯片

图 15 计算精度和水滴粒子数量的关系

Figure 15. Relationship between computational precision and the number of water droplets

下载: 全尺寸图片幻灯片

图 16 着陆构型水滴收集率云图

Figure 16. Cloud chart of water drop collection efficiency of landing configuration

下载: 全尺寸图片幻灯片

图 17 着陆构型截面上各部件的位置关系图

Figure 17. Position relationship diagram of various components on landing configuration section

下载: 全尺寸图片幻灯片

图 18 飞机各部件截面收集率系数分布图（V=90 m/s, D_eq=40 μm, α=0°）

Figure 18. Distribution diagram of sectional collection coefficient of aircraft components （V=90 m/s, D_eq=40 μm, α=0°）

下载: 全尺寸图片幻灯片

参考文献(24)

[1]	LYNCH F T,KHODADOUST A. Effects of ice accretions on aircraft aerodynamics[J]. Progress in Aerospace Sciences,2001,37(8): 669-767.
[2]	孙志国,朱程香,朱春玲. 飞机结冰数值仿真软件开发[J]. 计算机仿真,2012,29(4): 104-107,111. SUN Zhiguo,ZHU Chengxiang,ZHU Chunling. Development of software for aircraft icing simulation[J]. Computer Simulation,2012,29(4): 104-107,111. (in Chinese doi: 10.3969/j.issn.1006-9348.2012.04.025 SUN Zhiguo, ZHU Chengxiang, ZHU Chunling. Development of software for aircraft icing simulation[J]. Computer Simulation, 2012, 29(4): 104-107, 111. (in Chinese) doi: 10.3969/j.issn.1006-9348.2012.04.025
[3]	DONG Wei,ZHU Jianjun,LEI Guilin,et al. Numerical simulation of hot air anti-icing charateristics of an aero-engine strut: AIAA 2015-0537 [R]. Reston,Virginia: AIAA,2015.
[4]	黄平,卜雪琴,刘一鸣,等. 混合相/冰晶条件下的结冰研究综述[J]. 航空学报,2022,43(5): 025178. HUANG Ping,BU Xueqing,LIU Yiming,et al. A review of icing research under mixed phase/ice crystal conditions[J]. Acta Aeronautica et Astronautica Sinica,2022,43(5): 025178. (in Chinese HUANG Ping, BU Xueqing, LIU Yiming, et al. A review of icing research under mixed phase/ice crystal conditions[J]. Acta Aeronautica et Astronautica Sinica, 2022, 43(5): 025178. (in Chinese)
[5]	易贤,桂业伟,朱国林,等. 运输机翼型结冰的计算和实验[J]. 航空动力学报,2011,26(4): 808-813. YI Xian,GUI Yewei,ZHU Guolin,et al. Experimental and computational investigation into ice accretion on airfoil of a transport aircraft[J]. Journal of Aerospace Power,2011,26(4): 808-813. (in Chinese YI Xian, GUI Yewei, ZHU Guolin, et al. Experimental and computational investigation into ice accretion on airfoil of a transport aircraft[J]. Journal of Aerospace Power, 2011, 26(4): 808-813. (in Chinese)
[6]	桑为民,蔡旸,鲁天. 变形破碎特性对SLD结冰过程影响[J]. 航空动力学报,2017,32(7): 1537-1544. SANG Weimin,CAI Yang,LU Tian. Effect of deformation and breakup characteristic on supercooled large droplet icing process[J]. Journal of Aerospace Power,2017,32(7): 1537-1544. (in Chinese SANG Weimin, CAI Yang, LU Tian. Effect of deformation and breakup characteristic on supercooled large droplet icing process[J]. Journal of Aerospace Power, 2017, 32(7): 1537-1544. (in Chinese)
[7]	WANG Qiang,YI Xian,LIU Yu,et al. Simulation and analysis of wind turbine ice accretion under yaw condition via an Improved Multi-Shot Icing Computational Model[J]. Renewable Energy,2020,162: 1854-1873. doi: 10.1016/j.renene.2020.09.107
[8]	常士楠,苏新明,邱义芬. 三维机翼结冰模拟[J]. 航空学报,2011,32(2): 212-222. CHANG Shinan,SU Xinming,QIU Yifen. Ice accretion simulation on three dimensional wings[J]. Acta Aeronautica et Astronautica Sinica,2011,32(2): 212-222. (in Chinese CHANG Shinan, SU Xinming, QIU Yifen. Ice accretion simulation on three dimensional wings[J]. Acta Aeronautica et Astronautica Sinica, 2011, 32(2): 212-222. (in Chinese)
[9]	易贤,王开春,马洪林,等. 大型风力机结冰过程水滴收集率三维计算[J]. 空气动力学学报,2013,31(6): 745-751. YI Xian,WANG Kaichun,MA Honglin,et al. 3-D numerical simulation of droplet collection efficiency in large-scale wind turbine icing[J]. Acta Aerodynamica Sinica,2013,31(6): 745-751. (in Chinese YI Xian, WANG Kaichun, MA Honglin, et al. 3-D numerical simulation of droplet collection efficiency in large-scale wind turbine icing[J]. Acta Aerodynamica Sinica, 2013, 31(6): 745-751. (in Chinese)
[10]	申晓斌,张志强,林贵平,等. 旋转部件复杂表面水滴撞击计算[J]. 空气动力学学报,2016,34(6): 709-713. SHEN Xiaobin,ZHANG Zhiqiang,LIN Guiping,et al. Droplet impingement calculation on complex suface of rotating part[J]. Acta Aerodynamica Sinica,2016,34(6): 709-713. (in Chinese SHEN Xiaobin, ZHANG Zhiqiang, LIN Guiping, et al. Droplet impingement calculation on complex suface of rotating part[J]. Acta Aerodynamica Sinica, 2016, 34(6): 709-713. (in Chinese)
[11]	孙志国,朱春玲. 三维机翼表面水滴撞击特性计算[J]. 计算物理,2011,28(5): 677-685. SUN Zhiguo,ZHU Chunling. Calculation of water-droplet impingement on wing surface[J]. Chinese Journal of Computational Physics,2011,28(5): 677-685. (in Chinese doi: 10.3969/j.issn.1001-246X.2011.05.006 SUN Zhiguo, ZHU Chunling. Calculation of water-droplet impingement on wing surface[J]. Chinese Journal of Computational Physics, 2011, 28(5): 677-685. (in Chinese) doi: 10.3969/j.issn.1001-246X.2011.05.006
[12]	周志宏,易贤,桂业伟,等. 水滴撞击特性的高效计算方法[J]. 空气动力学学报,2014,32(5): 712-716. ZHOU Zhihong,YI Xian,GUI Yewei,et al. An efficient methpd tp simulate water drpplet trajectpry and impingement[J]. Acta Aerodynamica Sinica,2014,32(5): 712-716. (in Chinese doi: 10.7638/kqdlxxb-2012.0179 ZHOU Zhihong, YI Xian, GUI Yewei, et al. An efficient methpd tp simulate water drpplet trajectpry and impingement[J]. Acta Aerodynamica Sinica, 2014, 32(5): 712-716. (in Chinese) doi: 10.7638/kqdlxxb-2012.0179
[13]	任靖豪,王强,刘宇,等. 大型商用运输机机翼增升构型水滴撞击特性计算[J]. 空气动力学学报,2021,39(1): 52-58,72. REN Jinghao,WANG Qiang,LIU Yu,et al. Numerical simulation of droplet impingement characteristics on a high-lift configuration of a large commercial transport aircraft[J]. Acta Aerodynamica Sinica,2021,39(1): 52-58,72. (in Chinese REN Jinghao, WANG Qiang, LIU Yu, et al. Numerical simulation of droplet impingement characteristics on a high-lift configuration of a large commercial transport aircraft[J]. Acta Aerodynamica Sinica, 2021, 39(1): 52-58, 72. (in Chinese)
[14]	POTAPCZUK M,BIDWELL C. Swept wing ice accretion modeling: AIAA 1990-756 [R]. Reston,Virigina: AIAA,1990.
[15]	WIDHALM M,RONZHEIMER A,MEYER J. Lagrangian particle tracking on large unstructured three-dimensional meshes: AIAA 2008-472 [R]. Reston,Virigina: AIAA,2008.
[16]	XIE Liang,LI Peizhe,CHEN Hang,et al. Robust and efficient prediction of the collection efficiency in icing accretion simulation for 3D complex geometries using the Lagrangian approach I: an adaptive interpolation method based on the restricted radial basis functions[J]. International Journal of Heat and Mass Transfer,2020,150: 119290.
[17]	HAMED A,DAS K,BASU D. Numerical simulations of ice droplet trajectories and collection efficiency on aero-engine rotating machinery: AIAA 2005-1248 [R]. Reston,Virigina: AIAA,2005.
[18]	BIDWELL C S. Particle Trajectory and icing analysis of the E (sup 3) turbofan engine using LEWICE3D Version 3:AIAA 2011-38-0048[R]. Reston, Virigina: AIAA, 2011.
[19]	韩雅慧,柯鹏,杨春信,等. 水滴撞击特性的粒子统计法研究[J]. 航空学报,2013,34(7): 1588-1595. HAN Yahui,KE Peng,YANG Chunxin,et al. Research on particle statistic method of water droplet impingement characteristics[J]. Acta Aeronautica et Astronautica Sinica,2013,34(7): 1588-1595. (in Chinese HAN Yahui, KE Peng, YANG Chunxin, et al. Research on particle statistic method of water droplet impingement characteristics[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(7): 1588-1595. (in Chinese)
[20]	赵海波,郑楚光. 描述颗粒沉积动力学演变过程的一种随机算法[J]. 空气动力学学报,2006,24(2): 141-146. ZHAO Haibo,ZHENG Chuguang. A stochastic algorithm for the dynamic evolution process by particle deposition[J]. Acta Aerodynamica Sinica,2006,24(2): 141-146. (in Chinese doi: 10.3969/j.issn.0258-1825.2006.02.001 ZHAO Haibo, ZHENG Chuguang. A stochastic algorithm for the dynamic evolution process by particle deposition[J]. Acta Aerodynamica Sinica, 2006, 24(2): 141-146. (in Chinese) doi: 10.3969/j.issn.0258-1825.2006.02.001
[21]	郭庆,徐甘生,赵洪利. 基于蒙特卡罗发动机竞争失效的下发仿真模型[J]. 航空动力学报,2019,34(3): 616-626. GUO Qing,XU Gansheng,ZHAO Hongli. Monte Carlo-based competitive failure delivery simulation model of engine[J]. Journal of Aerospace Power,2019,34(3): 616-626. (in Chinese GUO Qing, XU Gansheng, ZHAO Hongli. Monte Carlo-based competitive failure delivery simulation model of engine[J]. Journal of Aerospace Power, 2019, 34(3): 616-626. (in Chinese)
[22]	赵钟,何磊,何先耀. 风雷(PHengLEI)通用CFD软件设计[J]. 计算机工程与科学,2020,42(2): 210-219. ZHAO Zhong,HE Lei,HE Xianyao. Design of general CFD software PHengLEI[J]. Computer Engineering & Science,2020,42(2): 210-219. (in Chinese doi: 10.3969/j.issn.1007-130X.2020.02.004 ZHAO Zhong, HE Lei, HE Xianyao. Design of general CFD software PHengLEI[J]. Computer Engineering & Science, 2020, 42(2): 210-219. (in Chinese) doi: 10.3969/j.issn.1007-130X.2020.02.004
[23]	TONG Xiaoling,LUKE E. Eulerian simulations of icing collection efficiency using a singularity diffusion model: AIAA 2005-1246 [R]. Reston,Virigina: AIAA,2005.
[24]	MORENCY F,BEAUGENDRE H,HABASHI W. FENSAP-ICE: effect of ice shapes on 3D eulerian droplet impingement: AIAA 2003-1223 [R]. Reston,Virigina: AIAA,2003.