基于欧拉法流线的三维表面水滴撞击特性计算方法

  • 申晓斌 ,
  • 叶泽坤 ,
  • 赵静宇 ,
  • 郭洁涛 ,
  • 林贵平
展开
  • 北京航空航天大学

收稿日期: 2025-09-03

  修回日期: 2025-12-02

  网络出版日期: 2025-12-08

A New Method Based on Eulerian Droplet Streamline for Impingement Characteristics Calculation on Three Dimensional Surfaces

  • SHEN Xiao-Bin ,
  • YE Ze-Kun ,
  • ZHAO Jing-Yu ,
  • GUO Jie-Tao ,
  • LIN Gui-Ping
Expand

Received date: 2025-09-03

  Revised date: 2025-12-02

  Online published: 2025-12-08

摘要

飞机及其发动机表面的水滴撞击特性计算是其结冰分析与防除冰系统设计的首要工作。为弥补传统欧拉法与拉格朗日法的不足,建立了一种基于欧拉法流线的三维表面水滴撞击特性计算方法。其仅求解水滴动量方程来获得速度分布,不需要计算水滴连续性方程;通过三维表面网格节点的逆向积分获得水滴流线,进而计算水滴撞击特性,不需要提前确定水滴释放位置与追踪大量的运动轨迹。将所建立的方法应用于三维圆球、发动机整流帽罩及进气道唇口的水滴收集系数计算,计算效率较高,结果与文献数据都吻合较好,验证了基于欧拉法流线的三维表面水滴撞击特性计算方法的可行性与有效性。工作可以为实际飞机与发动机的结冰分析与防除冰系统设计提供参考。

本文引用格式

申晓斌 , 叶泽坤 , 赵静宇 , 郭洁涛 , 林贵平 . 基于欧拉法流线的三维表面水滴撞击特性计算方法[J]. 航空学报, 0 : 1 -0 . DOI: 10.7527/S1000-68932025.32751

Abstract

The computation of droplet impingement characteristics on the three-dimensional surface is the primary task for ice accretion analysis and anti-icing/de-icing system design for an aircraft and its engine. To overcome the shortcomings of the traditional Eulerian method and Lagrangian method, a streamline-based Eulerian method was established to obtain three-dimensional surface water droplet impingement characteristic. This method only solves the momentum equation to derive the velocity distribution, and there is no need to calculate the droplet continuity equation. Droplet streamlines are generated via backward integration from the vertices of the three-dimensional surface mesh, allowing impingement characteristics to be calculated without predetermination for droplet release locations or tracking a large number of droplet trajectories. The proposed method is applied to compute droplet collection efficiencies on a three-dimensional sphere, a spinner, and an engine inlet. The results show high computational efficiency and good agreement with the data in the literature, thereby confirming the feasibility and effectiveness of the streamline-based Eulerian method for three-dimensional droplet impingement analysis. This work can provide a reference for the ice accretion analysis and anti-icing/de-icing system design of aircraft and engines.

参考文献

[1] WAGDI G H. Handbook of numerical simulation of in-flight icing [M]. Cham, Switzerland: Springer, 2024: 1-27.
[2] 赵宾宾, 张恒, 李杰. 翼型结冰状态复杂分离流动数值模拟综述 [J]. 航空学报, 2023, 44(1): 627211.
ZHAO B B, ZHANG H, LI J. Review of numerical simulation on complex separated flow of iced airfoil [J]. Acta Aeronautica et Astronautica Sinica, 2023, 44(1): 627211 (in Chinese).
[3] 郑梅, 冯丽娟, 秦娜, 等. 短舱防冰系统三维内外流耦合计算方法 [J]. 航空学报, 2023, 44(1): 627425.
ZHENG M, FENG L J, QIN N, et al. 3D computational method for conjugate heat transfer between internal and external flow of nacelle anti-icing system [J]. Acta Aeronautica et Astronautica Sinica, 2023, 44(1): 627425 (in Chinese).
[4] 李浩然,段玉宇,张宇飞,等.结冰模拟软件AERO-ICE中的关键数值方法 [J]. 航空学报, 2021, 42(S1): 726371.
LI H R, DUAN Y Y, ZHANG Y F, et al. Numerical method of ice-accretion software AERO-ICE [J]. Acta Aeronautica et Astronautica Sinica, 2021, 42(S1): 726371 (in Chinese).
[5] 杨倩, 郭晓峰, 李芹, 等. 基于POD 和代理模型的热气防冰性能预测方法 [J]. 航空学报, 2023, 44(1): 626992.
YANG Q, GUO X F, LI Q, et al. Hot air anti-icing performance estimation method based on POD and surrogate model [J]. Acta Aeronautica et Astronautica Sinica, 2023, 44(1): 626992 (in Chinese).
[6] 石达志, 桑为民, 李世杰, 等. 基于LBM 的结冰表面水滴撞击特性数值模拟 [J]. 航空学报, 2023, 44(S2): 729192.
SHI D Z, SANG W M, LI S J, et al. Numerical simulation of water droplet impact characteristics on icing surfaces based on LBM [J]. Acta Aeronautica et Astronautica Sinica, 2023, 44(S2): 729192 (in Chinese).
[7] LAURENDEAU E, BOURGAULT-COTE S, OZCER I, et al. Summary from the 1st AIAA Ice Prediction Workshop [C]// AIAA AVIATION 2022 Forum. Chicago: 2022.
[8] IULIANO E, BRANDI V, MINGIONE G, et al. Water impingement prediction on multi-element airfoils by means of eulerian and lagrangian approach with viscous and inviscid air flow [C]// 44th AIAA Aerospace Sciences Meeting and Exhibit. Reno, Nevada, 2006.
[9] WRIGHT W B. User manual for the NASA Glenn ice accretion code LEWICE version 2.2.2 [R]. NASA-CR-2002-211793, 2002.
[10] SHEN X B, XIAO C H, NING Y J, et al. Research on the methods for obtaining droplet impingement characteristics in the lagrangian framework [J]. Aerospace, 2024, 11(3):172.
[11] BOURGAULT Y, BOUTANIOS Z, HABASHI G W. Three-dimensional eulerian approach to droplet impingement simulation using FENSAP - ICE, part 1: model, algorithm, and validation [J]. Journal of Aircraft, 2000, 37(1): 95-103.
[12] WIROGO S, SRIRAMBHATLA S. An eulerian method to calculate the collection efficiency on two and three dimensional bodies [C]// 1st Aerospace Sciences Meeting and Exhibit. Reno: AIAA, 2003.
[13] WANG C, CHANG S, WU H, Lagrangian approach for simulating supercooled large droplets’ impingement effect [J], Journal of Aircraft, 2015, 52(2): 524-537.
[14] LU S, CHEN W, ZHANG D, et al. Investigation on phase transition and collection characteristics of non-spherical ice crystals with eulerian and lagrangian methods [J]. Aerospace, 2024, 11(4): 299.
[15] XIE L, LI P Z, CHEN H, 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.
[16] SENGUPTA B, ESMAEILIFAR E, ARAGHIZADEH M, et al. Rotor-fuselage-intake aerodynamics and icing using vortex and eulerian–lagrangian computational fluid dynamics methods [J], AIAA Journal, 2025, 63(3): 909-930.
[17] ZHU C X, TAO M J, ZHAO N, et al. Study of droplet shadow zone of aircraft wing with diffusion effects [J]. AIAA Journal, 2019, 57(8): 3339-3348.
[18] 陈维建. 飞机机翼结冰的数值模拟研究 [D]. 南京: 南京航空航天大学, 2007: 59-96.
CHEN W J. Numerical simulation of aircraft wing icing [D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2007 (in Chinese).
[19] WU J, XU Q Y, WU F, et al. Droplet collection efficiency regularity of NACA0012 airfoil based on the eulerian method [J]. Aerospace, 2023, 10(5): 412.
[20] TONG X L, LUKE E. Eulerian simulations of icing collection efficiency using a singularity diffusion model [J]. Engineering Applications of Computational Fluid Mechanics, 2010, 4(4): 483-495.
[21] LIU Y, QU J, YI X, et al. A monte carlo lagrangian droplet solver with backpropagation neural network for aircraft icing simulation [J]. Transactions of Nanjing University of Aeronautics and Astronautics,2023, 40(5): 566-577.
[22] 曾涛, 王昭力, 熊华杰, 等. 复杂构型水滴收集率的拉格朗日快速算法研究 [J]. 航空动力学报, 2025, 40(7): 20220539.
ZENG T, WANG Z L, XIONG H J, et al. Research on a fast Lagrangian algorithm for water droplet collection efficiency in complex configurations [J]. Journal of Aerospace Power, 2025, 40(7): 20220539 (in Chinese).
[23] 任靖豪, 王强, 李维浩, 等. 基于梯度下降的水滴收集率计算方法 [J]. 航空学报, 2023, 44(4): 126381.
REN J H, WANG Q, LI W H, et al. A prediction algorithm of collection efficiency based on gradient descent method [J]. Acta Aeronautica et Astronautica Sinica, 2023, 44(4): 126381(in Chinese).
[24] BELLOSTA T, BALDAN G, SIRIANNI G, et al. Lagrangian and eulerian algorithms for water droplets in in-flight ice accretion [J]. Journal of Computational and Applied Mathematics, 2023, 429: 115230.
[25] WANG S, LOTH E. Droplet impact efficiency on aerodynamic surfaces with a globally eulerian/locally lagrangian method [J]. Journal of Aircraft, 2017, 54(1): 104-113.
[26] ZAYNI M K, BLANCHET M, LAURENDEAU E. Lagrangian particle tracking for ice accretion applications [C]// AIAA AVIATION FORUM AND ASCEND 2024. Las Vegas, Nevada:2024.
[27] 李静, 刘振侠, 胡剑平. 三维S形进气道沿程结冰参数分析方法 [J]. 科学技术与工程, 2018, 18(21): 141-145.
LI J, LIU Z X, HU J P. Analysis method of icing parameters along three-dimensional S-shaped inlet [J]. Science Technology and Engineering, 2018, 18(21): 141-145 (in Chinese).
[28] 吴佩佩, 晏涛, 马赛强, 等. 发动机唇口电热防冰系统性能仿真 [J]. 航空动力学报, 2020, 35(10): 2025-2063.
WU P P, YAN T, MA S Q, et al. Numerical simulation of performance of engine lip electro-thermal anti-icing system [J]. Journal of Aerospace Power, 2020, 35(10): 2025-2063 (in Chinese).
[29] TORMEN D, ZANON A, GENNARO M D. Ice protection system design for the next generation civil tiltrotor engine intake [C]// SAE International Conference on Icing of Aircraft, Engines, and Structures. Austria, 2023.
[30] JUNG S K, MYONG R S. A second-order positivity-preserving finite volume upwind scheme for air-mixed droplet flow in atmospheric icing [J]. Computers & Fluids, 2013, 86:459-469.
[31] BLANCHET M, BOURGAULT-COTE S, LAURENDEAU E. Conservative hyperbolic droplet solver for aircraft icing [C]// AIAA AVIATION 2022 Forum, Chicago: 2022.
[32] 陈希, 招启军. 考虑遮蔽区影响的旋翼三维水滴撞击特性计算新方法[J]. 航空学报, 2017, 38(6): 120745.
CHEN X, ZHAO Q J. New method for predicting 3-d water droplet impingement on rotor considering influence of shadow zone. Acta Aeronautica et Astronautica Sinica, 2017, 38(6): 120745 (in Chinese).
[33] SOTOMAYOR-ZAKHAROV D, BANSMER S. Finite - volume Eulerian solver for simulation of particle - laden flows for icing applications [J]. Computers & Fluids, 2021, 228: 105009.
[34] BOURGAULT Y, HABASHI W G, DOMPIERRE J, et al. A finite element method study of Eulerian droplets impingement models [J]. International Journal for Numerical Methods in Fluids, 1999, 29(4): 429-449.
[35] LU Y J, SUN Q Q, CHOI K S, et al. Joint effects of virtual surfaces on anti-icing and drag reduction [J]. AIAA Journal, 2025, 63(4): 1502-1511.
[36] GAO X, QIU B, WANG Z, et al. Influence of spinner shape on droplet impact over rotating spinners [J]. Aerospace, 2023, 10(1): 68.
[37] BIDWELL C, STANLEY MOHLER Y R. Collection efficiency and ice accretion calculations for a sphere, a swept MS(1)-317 wing, a swept NACA-0012 wing tip, an axisymmetric inlet, and a Boeing 737-300 inlet [C]// 33rd Aerospace Sciences Meeting and Exhibit. Reston:AIAA, 1995.
[38] JIA W, ZHANG F. Numerical investigation of supercooled large droplets impingement characteristics of the rotating spinner [J]. International Journal of Aerospace Engineering, 2024, 2024: 1683744.
[39] ZANON A, PAGé J, TORMEN D, et al. 1st AIAA ice prediction workshop: AIT numerical simulation results [C]// AIAA AVIATION 2022 Forum, Chicago: 2022.
[40] IULIANO E, MINGIONE G, et al. An Eulerian Approach to Three-Dimensional Droplet Impingement Simulation in Icing Environment[C]// AIAA Atmospheric and Space Environments Conference, Toronto: 2010.
文章导航

/