CN116562192B - 一种飞机结冰冰形预测方法、装置、设备及存储介质 - Google Patents

一种飞机结冰冰形预测方法、装置、设备及存储介质 Download PDF

Info

Publication number
CN116562192B
CN116562192B CN202310824277.5A CN202310824277A CN116562192B CN 116562192 B CN116562192 B CN 116562192B CN 202310824277 A CN202310824277 A CN 202310824277A CN 116562192 B CN116562192 B CN 116562192B
Authority
CN
China
Prior art keywords
water
mass flow
control body
overflow
water mass
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202310824277.5A
Other languages
English (en)
Other versions
CN116562192A (zh
Inventor
李伟斌
崔兴达
宋超
刘钒
蓝庆生
余永刚
周铸
刘红阳
肖中云
王浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202310824277.5A priority Critical patent/CN116562192B/zh
Publication of CN116562192A publication Critical patent/CN116562192A/zh
Application granted granted Critical
Publication of CN116562192B publication Critical patent/CN116562192B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/28—Fuselage, exterior or interior
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/12—Timing analysis or timing optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Optimization (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Fluid Mechanics (AREA)
  • Transportation (AREA)
  • Manufacturing & Machinery (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Automation & Control Theory (AREA)
  • Computational Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种飞机结冰冰形预测方法、装置、设备及存储介质，应用于冰形预测领域，该方法通过在迭代计算活动水质量流量时引入溢流水流动时间步，将溢流水流动时间步和网格每一个边界的溢流速度相乘得到每一个迭代步的流出水质量流量；将相邻控制体的流出水质量流量确定为当前控制体的流入水质量流量；根据流出水质量流量和流入水质量流量迭代计算活动水质量流量直至收敛，从而确定最终的目标冻结水冻结速率以预测结冰冰形。避免了现有技术中在迭代计算活动水质量流量时，使每一次迭代每个网格的流出水完全流出到相邻的网格，同时该网格也接收来自上游的溢流水，忽视了不同网格尺寸和溢流速度的差异而导致预测冰形不精确的问题。

Description

一种飞机结冰冰形预测方法、装置、设备及存储介质
技术领域
本发明涉及冰形预测领域，特别涉及一种飞机结冰冰形预测方法、装置、设备及计算机可读存储介质。
背景技术
飞机遭遇低温云雾环境时，机翼、平尾等会发生结冰现象。飞机结冰会改变飞机的气动外形，直接影响飞机气动性能，严重时威胁飞行安全。在飞机设计阶段，数值模拟是所设计外形的结冰评估的重要手段之一，主要是进行冰形的几何形状和结冰后气动性能的计算。在飞行温度较高时（通常大于-15℃），过冷水滴撞击至飞机时，不会立即发生冻结，一部分水滴形成水膜后会发生溢流现象，对液膜溢流特性的准确预测，才能保证结冰冰形数值模拟的准确性。
现有技术在进行计算活动水质量流量的计算时，认为每一次迭代，每个网格的流出水完全流出到相邻的网格中，同时该网格也接收到了来自上游的溢流水。该做法仅以速度的方向作为溢流水流动的分配依据，忽视了网格尺寸和溢流速度的差异，即忽视了溢流水在不同网格上的流动时间差异，导致飞机结冰冰形预测误差较大，不精确。
发明内容
本发明的目的在于提供一种飞机结冰冰形预测方法、装置、设备及存储介质，应用于冰形预测领域，该方法通过在迭代计算活动水质量流量时引入溢流水流动时间步，迭代计算活动水质量流量直至收敛，从而确定最终的目标冻结水冻结速率以预测结冰冰形，避免了现有技术中在迭代计算活动水质量流量时，忽视了不同网格尺寸和溢流速度的差异而导致冰形预测不精确的问题。
为解决上述计算问题，本发明提供一种飞机结冰冰形预测方法，包括：
S101：确定溢流水流动时间步、控制体的活动水质量流量和网格边界的溢流速度；
S102：根据所述溢流速度、所述溢流水流动时间步和所述活动水质量流量计算所述控制体的流出水质量流量；
S103：根据相邻控制体的所述流出水质量流量确定所述控制体的流入水质量流量；
S104：根据所述控制体的所述流出水质量流量和所述流入水质量流量计算迭代活动水质量流量；
S105：判断所述迭代活动水质量是否收敛；若否，执行步骤S106；若是，执行步骤S107；
S106：将所述迭代活动水质量确定为所述活动水质量流量，并重新开始执行步骤S102；
S107：根据所述迭代活动水质量流量确定目标冻结速率，并根据所述目标冻结速率计算飞机预测冰形。
可选的，所述步骤S102，包括：
将所述溢流速度、所述溢流水流动时间步和所述活动水质量流量输入第一模型计算网格的每个所述网格边界的所述流出水质量流量；
将全部所述网格边界的所述流出水质量流量相加得到所述控制体的所述流出水质量流量，所述第一模型的表达式为：
；
式中，为第i个网格边界长度，/>为网格面积，/>为所述溢流水流动时间步，/>为第i个所述网格边界的所述溢流速度，/>为所述活动水质量流量，/>为第i个所述网格边界的所述流出水质量流量。
可选的，所述步骤S101，包括：
确定网格面积、网格边界长度；
将所述网格面积、所述网格边界长度和所述网格边界的所述溢流速度输入第二模型确定所述溢流水流动时间步，所述第二模型的表达式为：
；
式中，为第j个所述网格边界对应的所述网格面积，/>为第j个所述网格边界长度，/>为第j个所述网格边界的所述溢流速度，/>为所述溢流水流动时间步。
可选的，所述步骤S101，包括：
确定所述网格边界的流出方向，在预设来流条件下确定附面层外空气速度；
计算所述附面层外空气速度在所述流出方向的速度投影，并将所述速度投影确定为所述溢流速度。
可选的，所述步骤S104，包括：
在预设来流条件下计算空气流场数据，根据所述空气流场数据计算所述控制体的撞击水质量流量；
根据表面温度确定蒸发水蒸发速率，并联立求解所述控制体的质量守恒方程与能量守恒方程，得到计算冻结水冻结速率；
根据所述计算冻结水冻结速率对冻结水冻结速率进行赋值，并根据所述冻结水冻结速率更新所述表面温度和所述蒸发水蒸发速率；
根据所述蒸发水蒸发速率、所述冻结水冻结速率、所述撞击水质量流量、所述流出水质量流量和所述流入水质量流量确定所述迭代活动水质量流量。
可选的，所述控制体的所述质量守恒方程为：
；
式中，为所述撞击水质量流量，/>为所述流入水质量流量，/>为所述蒸发水蒸发速率，/>为所述流出水质量流量，/>为所述冻结水冻结速率，i为相邻网格的序号；
所述控制体的所述能量守恒方程为：
；
式中，为液态水冻结离开所述控制体的能量，/>为液态水蒸发离开所述控制体的能量，/>为与所述控制体相碰撞的水滴的能量，/>为气流摩擦对控制体表面的加热量，/>为气流与所述控制体表面的对流换热量，/>为流入所述控制体的溢流水能量，/>为流出所述控制体的所述溢流水能量。
可选的，所述根据所述计算冻结水冻结速率对冻结水冻结速率进行赋值，包括：
将所述计算冻结水冻结速率代入第三模型对所述冻结水冻结速率进行赋值，所述第三模型的表达式为：
；
式中，为所述撞击水质量流量，/>为所述流入水质量流量，/>为所述蒸发水蒸发速率，/>为所述流出水质量流量，/>为所述冻结水冻结速率，/>为所述计算冻结水冻结速率。
可选的，所述在预设来流条件下计算空气流场数据，根据所述空气流场数据计算所述控制体的撞击水质量流量，包括：
在所述预设来流条件下求解不可压N-S方程得到的外形的压力和附外层空气速度；
应用Euler方法根据所述压力和所述附外层空气速度求解水滴的输运方程，计算所述控制体的所述撞击水质量流量。
为解决上述技术问题，本发明提供一种飞机结冰冰形预测装置，包括：
数据确定模块，用于确定溢流水流动时间步、控制体的活动水质量流量和网格边界的溢流速度；
第一计算模块，用于根据所述溢流速度、所述溢流水流动时间步和所述活动水质量流量计算所述控制体的流出水质量流量；
第二计算模块，用于根据相邻控制体的所述流出水质量流量确定所述控制体的流入水质量流量；
第三计算模块，用于根据所述控制体的所述流出水质量流量和所述流入水质量流量计算迭代活动水质量流量；
判断模块，用于判断所述迭代活动水质量是否收敛；若否，执行第一执行模块；若是，执行第二执行模块；
第一执行模块，用于将所述迭代活动水质量确定为所述活动水质量流量，并重新开始执行所述第一计算模块；
第二执行模块，用于根据所述迭代活动水质量流量确定目标冻结速率，并根据所述目标冻结速率计算飞机预测冰形。
为解决上述技术问题，本发明还提供一种飞机结冰冰形预测设备，包括：
存储器，用于储存计算机程序；
处理器，用于执行所述计算机程序时实现任一项上述所述飞机结冰冰形预测方法。
为解决上述技术问题，本发明还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机可执行指令，所述计算机可执行指令被处理器执行时，实现任一项上述所述飞机结冰冰形预测方法。
可见，本发明方法通过在迭代计算活动水质量流量时引入溢流水流动时间步，将溢流水流动时间步和网格每一个边界的溢流速度相乘得到每一个迭代步的流出水质量流量；将相邻控制体的流出水质量流量确定为当前控制体的流入水质量流量；根据流出水质量流量和流入水质量流量迭代计算活动水质量流量直至收敛，从而确定最终的目标冻结水冻结速率以预测结冰冰形。避免了现有技术中在迭代计算活动水质量流量时，每一次迭代每个网格的流出水完全流出到相邻的网格，同时该网格也接收来自上游的溢流水，忽视了不同网格尺寸和溢流速度的差异而导致预测冰形不精确的问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为本发明实施例所提供的一种飞机结冰冰形预测方法的流程图；
图2为本发明实施例所提供的一种NACA0012翼型的外形图；
图3为本发明实施例所提供的一种NACA0012翼型外形的网格图；
图4为本发明实施例所提供的一种翼型中间截面第一迭代步的冻结水冻结速率分布图；
图5为本发明实施例所提供的一种翼型中间截面第一迭代步的表面温度分布图；
图6为本发明实施例所提供的一种翼型中间截面第一迭代步的蒸发水蒸发速率分布图；
图7为本发明实施例所提供的一种翼型中间截面第一迭代步的流出水质量流量分布图；
图8为本发明实施例所提供的一种翼型中间截面最终步冻结水冻结速率分布图；
图9为本发明实施例所提供的一种翼型中间截面最终步蒸发水蒸发速率分布图；
图10为本发明实施例所提供的一种翼型中间截面最终步流出水质量流量分布图；
图11为本发明实施例所提供的一种飞机结冰冰形预测装置的结构框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
以下结合图1，图1为本发明实施例所提供的一种飞机结冰冰形预测方法的流程图，该方法可以包括：
S101：确定溢流水流动时间步、控制体的活动水质量流量和网格边界的溢流速度。
本实施例中的网格为在飞机模型表面划分的网格，本实施例中的控制体是一种特定的几何体，空间中可以流动流体的体积，每一个网格上都对应一个控制体。本实施例中的溢流水流动时间步为溢流水在一个迭代步内所流动的时间，本实施例并不限定溢流水流动时间步的确定方式，一般情况下，在一个迭代步内，溢流水不能出现横跨网格的情况，所以在本实施例中可以将溢流水横跨一个网格的最小时间作为溢流水流动时间步。本实施例中的溢流水流动时间步需要通过网格边界的溢流速度参与求得，网格边界的溢流速度是溢流水沿着网格边界流出方向的速度，例如一个网格拥有四个边界，每个边界都有其对应的溢流速度。本实施例并不限定溢流速度的计算方式，在本实施例中，可以在给定的预设来流条件下计算飞机外形的流场数据，本实施例中可以在预设来流条件下求解不可压N-S方程计算得到附外层空气速度，根据流场数据计算附外层空气速度，附外层空气速度在边界流出方向上的投影即为该边界的溢流速度。
本实施例并不限定溢流水流动时间步的计算方式，一般可以先确定所有的网格面积和所有网格的网格边界长度；将网格面积及对应网格边界长度和网格边界的溢流速度输入第二模型确定溢流水流动时间步，第二模型的表达式为：
；
式中，为第j个网格边界对应的网格面积，/>为第j个网格边界长度，/>为第j个网格边界的溢流速度，/>为溢流水流动时间步。
本实施例中的j为所有网格边界的数量，本实施例并不限定第j个网格边对应的网格面积，一般一个网格边界对应两个网格面积，在本实施中可以根据实际应用进行选取。本实施例中所确定的溢流水流动时间步为横跨一个网格所需要的最小时间，在实际应用中，溢流水流动时间步也可以小于横跨一个网格所需要的最小时间，但在计算过程中迭代的次数增加。
在本实施例中，可以确定活动水质量流量，本实施例中的活动水质量流量为控制体中初始的活动水质量流量，本实施例并不限定活动水质量流量的确定方式，在本实施例中的活动水质量流量与迭代活动水质量流量的计算方式一致，区别在于用于计算活动水质量流量的流出水质量流量和流入水质量流量等数值为初始数值。具体计算方式可以为，在预设来流条件下求解不可压N-S方程得到的外形的压力和附外层空气速度，应用Euler（欧拉）方法根据压力和附外层空气速度求解水滴的输运方程，计算控制体的撞击水质量流量，撞击水质量流量为预设来流撞击进控制体的水滴质量流量。
初始化控制体的冻结水冻结速率、活动水质量流量、流入水质量流量和流出水质量流量为0，预设表面温度Ts=0℃，根据表面温度确定蒸发水蒸发速率。根据表面温度联立求解质量守恒方程和能量守恒方程求解计算冻结水冻结速率，对于每个控制体来说控制体的质量守恒方程为：
；
式中，为撞击水质量流量，/>为流入水质量流量，/>为蒸发水蒸发速率，/>为流出水质量流量，/>为冻结水冻结速率，i为相邻网格的序号；
控制体的能量守恒方程为：
；
式中，为液态水冻结离开控制体的能量，/>为液态水蒸发离开控制体的能量，/>为与控制体相碰撞的水滴的能量，/>为气流摩擦对控制体表面的加热量，/>为气流与控制体表面的对流换热量，/>为流入控制体的溢流水能量，/>为流出控制体的溢流水能量。
具体计算过程可以根据能量守恒方程求得的数值，根据/>可以求得计算冻结水冻结速率，根据计算冻结水冻结速率确定冻结水冻结速率的值，具体可以如下式所示，式中，/>为撞击水质量流量，/>为流入水质量流量，/>为蒸发水蒸发速率，/>为流出水质量流量，/>为冻结水冻结速率，/>为计算冻结水冻结速率：
；
根据赋值之后的冻结水冻结速率，联立控制体的质量守恒方程和能量守恒方程求解下一迭代步的表面温度，根据更新后的表面温度更新下一迭代步的蒸发水蒸发速率。
根据蒸发水蒸发速率、冻结水冻结速率、撞击水质量流量、流出水质量流量和流入水质量流量确定活动水质量流量，如下所示：
；
式中，为活动水质量流量，/>为撞击水质量流量，/>为流入水质量流量，/>为蒸发水蒸发速率，/>为流出水质量流量，/>为冻结水冻结速率。
S102：根据溢流速度、溢流水流动时间步和活动水质量流量计算控制体的流出水质量流量。
S103：根据相邻控制体的流出水质量流量确定控制体的流入水质量流量。
在确定初始的活动水质量流量之后，根据溢流速度、溢流水流动时间步和活动水质量流量计算控制体的流出水质量流量，具体计算方式可以将溢流速度、溢流水流动时间步和活动水质量流量输入第一模型计算网格的每个网格边界的流出水质量流量，第一模型的表达式为：
；
式中，为第i个网格边界长度，/>为网格面积，/>为溢流水流动时间步，/>为第i个网格边界的溢流速度，/>为活动水质量流量，/>为第i个网格边界的流出水质量流量。
将全部网格边界的流出水质量流量相加得到控制体的流出水质量流量。
在本实施例中，可以认为相邻控制体从与控制体共用网格边界流出的流出水质量流量为该控制体的流入水质量流量，即在本实施例中可以根据相邻控制体的流出水质量流量确定控制体的流入水质量流量。
S104：根据控制体的流出水质量流量和流入水质量流量计算迭代活动水质量流量。
本实施例中，在计算得到流出水质量流量和流入水质量流量之后，根据流出水质量流量和流入水质量流量计算迭代活动水质量流量，具体的计算方式与活动水质量流量的计算方式一致，所使用的表面温度和蒸发水蒸发速率等数值为更新后的数值。
本实施例中迭代活动水质量流量的计算过程是一个循环迭代的过程，活动水质量流量相当于初始步迭代活动水质量流量的确定，在每一个迭代步中更新下一迭代步所需数值。
S105：判断迭代活动水质量是否收敛；若否，执行步骤S106；若是，执行步骤S107。
S106：将迭代活动水质量确定为活动水质量流量，并重新开始执行步骤S102。
S107：根据迭代活动水质量流量确定目标冻结速率，并根据目标冻结速率计算飞机预测冰形。
本实施例通过判断每一迭代步的迭代活动水质量流量是否收敛来判断溢流水是否稳定，本实施例并不限定具体收敛的判定方式，可以是当每个控制体中当前迭代活动水质量流量与上一步迭代活动水质量流量的差值小于预设阈值时，判断迭代活动水质量流量收敛；也可以是每个控制体当前迭代活动水质量流量的和与上一步迭代活动水质量流量和的差值小于预设阈值时，判断迭代活动水质量流量收敛。
若判断迭代活动水质量流量收敛，可以将计算当前步迭代活动水质量流量过程中的冻结水冻结速率确定为目标冻结速率，并根据目标冻结速率计算飞机预测冰形。
若迭代活动水质量流量仍未收敛，可以继续执行迭代的过程，直至判断迭代活动水质量流量收敛。
本实施例通过在迭代计算活动水质量流量时引入溢流水流动时间步，迭代计算活动水质量流量直至收敛，从而确定最终的目标冻结水冻结速率以预测结冰冰形，避免了现有技术中在迭代计算活动水质量流量时，忽视了不同网格尺寸和溢流速度的差异而导致冰形预测不精确的问题。
以下本发明实施例所提供的一种飞机结冰冰形预测方法的具体实施例，其中，飞机计算外形为NACA0012翼型，弦长0.5334m，展长1m，如图2所示，以翼型前缘点为坐标原点构建直角坐标系，图2中x表示翼型表面点的横坐标，y表示翼型表面点的纵坐标，上表面点的纵坐标（y）为正值，下表面点纵坐标（y）为负值。来流攻角0°，速度105.4m/s，液态水含量1g/m3，中值粒径20μm，温度为-10℃。采用结构网格进行计算，外形的表面网格如图3所示，网格规模为40万。本具体实施例计算得到的数据保存在数据文件中，保存于数据文件中，该具体实施例可以包括：
在预设来流条件下求解外形表面压力、附外层空气速度和撞击水质量流量，并根据附外层空气速度在网格边界流出方向上的投影得到网格边界的溢流速度；
初始化控制体的冻结水冻结速率、活动水质量流量、流入水质量流量和流出水质量流量为0，预设表面温度Ts=0℃，根据表面温度确定蒸发水蒸发速率。根据表面温度联立求解质量守恒方程和能量守恒方程求解计算冻结水冻结速率；
根据计算冻结水冻结速率确定冻结水冻结速率的值，翼型中间截面上第一迭代步的冻结水冻结速率分布可以如图4所示，图4横坐标（y）与图2纵坐标（y）定义相同，纵坐标表示冻结水冻结速率。根据赋值之后的冻结水冻结速率，联立控制体的质量守恒方程和能量守恒方程求解下一迭代步的表面温度，翼型中间截面上第一迭代步得到的表面温度分布可以如图5所示，图5横坐标（y）与图2纵坐标（y）定义相同，纵坐标表示表面温度。根据更新后的表面温度更新下一迭代步的蒸发水蒸发速率，翼型中间截面上第一迭代步得到的蒸发水蒸发速率分布可以如图6所示，图6横坐标（y）与图2纵坐标（y）定义相同，纵坐标表示蒸发水蒸发速率；
根据蒸发水蒸发速率、冻结水冻结速率、撞击水质量流量、流出水质量流量和流入水质量流量确定活动水质量流量；根据溢流速度、溢流水流动时间步和活动水质量流量计算控制体的流出水质量流量；翼型中间截面上第一迭代步得到的流出水质量流量分布可以如图7所示，图7横坐标（y）与图2纵坐标（y）定义相同，纵坐标表示流出水质量流量；根据相邻控制体流出的流出水质量流量确定当前控制体的流入水质量流量；
根据控制体的流出水质量流量和流入水质量流量计算迭代活动水质量流量。判断迭代活动水质量是否收敛；若否，将迭代活动水质量确定为活动水质量流量，并重新迭代计算迭代活动水质量流量；若是，根据迭代活动水质量流量确定目标冻结速率，并根据目标冻结速率计算飞机预测冰形。
本具体实施例中，在迭代活动水质量流量收敛后翼型中间截面上最终的冻结水冻结速率分布可以如图8所示，图8横坐标（y）与图2纵坐标（y）定义相同，纵坐标表示冻结水冻结速率；翼型中间截面上最终的蒸发水蒸发速度分布可以如图9所示，图9横坐标（y）与图2纵坐标（y）定义相同，纵坐标表示蒸发水蒸发速率；翼型中间截面上最终的流出水质量流量分布可以如图10所示，图10横坐标（y）与图2纵坐标（y）定义相同，纵坐标表示流出水质量流量。
以下结合图11，图11为本发明实施例所提供的一种飞机结冰冰形预测装置的结构框图，该装置可以包括：
数据确定模块100，用于确定溢流水流动时间步、控制体的活动水质量流量和网格边界的溢流速度；
第一计算模块200，用于根据所述溢流速度、所述溢流水流动时间步和所述活动水质量流量计算所述控制体的流出水质量流量；
第二计算模块300，用于根据相邻控制体的所述流出水质量流量确定所述控制体的流入水质量流量；
第三计算模块400，用于根据所述控制体的所述流出水质量流量和所述流入水质量流量计算迭代活动水质量流量；
判断模块500，用于判断所述迭代活动水质量是否收敛；若否，执行第一执行模块；若是，执行第二执行模块；
第一执行模块600，用于将所述迭代活动水质量确定为所述活动水质量流量，并重新开始执行所述第一计算模块；
第二执行模块700，用于根据所述迭代活动水质量流量确定目标冻结速率，并根据所述目标冻结速率计算飞机预测冰形。
基于上述实施例，本发明通过在迭代计算活动水质量流量时引入溢流水流动时间步，迭代计算活动水质量流量直至收敛，从而确定最终的目标冻结水冻结速率以预测结冰冰形，避免了现有技术中在迭代计算活动水质量流量时，忽视了不同网格尺寸和溢流速度的差异而导致冰形预测不精确的问题。
基于上述各实施例，所述第一计算模块，可以包括：
第一模型单元，用于将所述溢流速度、所述溢流水流动时间步和所述活动水质量流量输入第一模型计算网格的每个所述网格边界的所述流出水质量流量；
求和单元，用于将全部所述网格边界的所述流出水质量流量相加得到所述控制体的所述流出水质量流量，所述第一模型的表达式为：
；
式中，为第i个网格边界长度，/>为网格面积，/>为所述溢流水流动时间步，/>为第i个所述网格边界的所述溢流速度，/>为所述活动水质量流量，/>为第i个所述网格边界的所述流出水质量流量。
基于上述各实施例，所述数据确定模块100，可以包括：
网格数据确定单元，用于确定网格面积、网格边界长度；
第二模型单元，用于将所述网格面积、所述网格边界长度和所述网格边界的所述溢流速度输入第二模型确定所述溢流水流动时间步，所述第二模型的表达式为：
；
式中，为第j个所述网格边界对应的所述网格面积，/>为第j个所述网格边界长度，/>为第j个所述网格边界的所述溢流速度，/>为所述溢流水流动时间步。
基于上述各实施例，所述数据确定模块100，可以包括：
方向确定单元，用于确定所述网格边界的流出方向，在预设来流条件下确定附面层外空气速度；
投影单元，用于计算所述附面层外空气速度在所述流出方向的速度投影，并将所述速度投影确定为所述溢流速度。
基于上述各实施例，所述第三计算模块400，可以包括：
流场计算单元，用于在预设来流条件下计算空气流场数据，根据所述空气流场数据计算所述控制体的撞击水质量流量；
冻结速率单元，用于根据表面温度确定蒸发水蒸发速率，并联立求解所述控制体的质量守恒方程与能量守恒方程，得到计算冻结水冻结速率；
更新单元，用于根据所述计算冻结水冻结速率对冻结水冻结速率进行赋值，并根据所述冻结水冻结速率更新所述表面温度和所述蒸发水蒸发速率；
计算单元，用于根据所述蒸发水蒸发速率、所述冻结水冻结速率、所述撞击水质量流量、所述流出水质量流量和所述流入水质量流量确定所述迭代活动水质量流量。
基于上述各实施例，所述控制体的所述质量守恒方程可以为：
；
式中，为所述撞击水质量流量，/>为所述流入水质量流量，/>为所述蒸发水蒸发速率，/>为所述流出水质量流量，/>为所述冻结水冻结速率，i为相邻网格的序号；
所述控制体的所述能量守恒方程为：
；
式中，为液态水冻结离开所述控制体的能量，/>为液态水蒸发离开所述控制体的能量，/>为与所述控制体相碰撞的水滴的能量，/>为气流摩擦对控制体表面的加热量，/>为气流与所述控制体表面的对流换热量，/>为流入所述控制体的溢流水能量，/>为流出所述控制体的所述溢流水能量。
基于上述各实施例，所述更新单元，可以包括：
第三模型子单元，用于将所述计算冻结水冻结速率代入第三模型对所述冻结水冻结速率进行赋值，所述第三模型的表达式为：
；
式中，为所述撞击水质量流量，/>为所述流入水质量流量，/>为所述蒸发水蒸发速率，/>为所述流出水质量流量，/>为所述冻结水冻结速率，/>为所述计算冻结水冻结速率。
基于上述各实施例，所述流场计算单元，可以包括：
求解子单元，用于在所述预设来流条件下求解不可压N-S方程得到的外形的压力和附外层空气速度；
计算子单元，用于应用Euler方法根据所述压力和所述附外层空气速度求解水滴的输运方程，计算所述控制体的所述撞击水质量流量。
基于上述实施例，本发明还提供了一种设备，该设备可以包括存储器和处理器，其中，该存储器中存有计算机程序，该处理器调用该存储器中的计算机程序时，可以实现上述实施例所提供的步骤。当然，该设备还可以包括各种必要的网络接口、电源以及其它零部件等。
本发明还提供了一种计算机可读存储介质，其上存有计算机程序，该计算机程序被执行终端或处理器执行时可以实现本发明实施例所提供的方法；该存储介质可以包括：U盘、移动硬盘、只读存储器（Read-Only Memory，ROM）、随机存取存储器（Random AccessMemory，RAM）、磁碟或者光盘等各种可以存储程序代码的介质。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其它实施例的不同之处，各个实施例之间相同或相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
最后，还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上对本发明所提供的一种飞机结冰冰形预测方法、装置、设备及存储介质进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (9)

1.一种飞机结冰冰形预测方法，其特征在于，包括：
S101：确定溢流水流动时间步、控制体的活动水质量流量和网格边界的溢流速度；
S102：根据所述溢流速度、所述溢流水流动时间步和所述活动水质量流量计算所述控制体的流出水质量流量；
S103：根据相邻控制体的所述流出水质量流量确定所述控制体的流入水质量流量；
S104：根据所述控制体的所述流出水质量流量和所述流入水质量流量计算迭代活动水质量流量；
S105：判断所述迭代活动水质量是否收敛；若否，执行步骤S106；若是，执行步骤S107；
S106：将所述迭代活动水质量确定为所述活动水质量流量，并重新开始执行步骤S102；
S107：根据所述迭代活动水质量流量确定目标冻结速率，并根据所述目标冻结速率计算飞机预测冰形；
其中，所述步骤S102，包括：
将所述溢流速度、所述溢流水流动时间步和所述活动水质量流量输入第一模型计算网格的每个所述网格边界的所述流出水质量流量；
将全部所述网格边界的所述流出水质量流量相加得到所述控制体的所述流出水质量流量，所述第一模型的表达式为：
；
式中，为第i个网格边界长度，/>为网格面积，/>为所述溢流水流动时间步，/>为第i个所述网格边界的所述溢流速度，/>为所述活动水质量流量，/>为第i个所述网格边界的所述流出水质量流量；
所述步骤S104，包括：
在预设来流条件下计算空气流场数据，根据所述空气流场数据计算所述控制体的撞击水质量流量；
根据表面温度确定蒸发水蒸发速率，并联立求解所述控制体的质量守恒方程与能量守恒方程，得到计算冻结水冻结速率；
根据所述计算冻结水冻结速率对冻结水冻结速率进行赋值，并根据所述冻结水冻结速率更新所述表面温度和所述蒸发水蒸发速率；
根据所述蒸发水蒸发速率、所述冻结水冻结速率、所述撞击水质量流量、所述流出水质量流量和所述流入水质量流量确定所述迭代活动水质量流量。
2.根据权利要求1所述的飞机结冰冰形预测方法，其特征在于，所述步骤S101，包括：
确定网格面积、网格边界长度；
将所述网格面积、所述网格边界长度和所述网格边界的所述溢流速度输入第二模型确定所述溢流水流动时间步，所述第二模型的表达式为：
；
式中，为第j个所述网格边界对应的所述网格面积，/>为第j个所述网格边界长度，为第j个所述网格边界的所述溢流速度，/>为所述溢流水流动时间步。
3.根据权利要求1所述的飞机结冰冰形预测方法，其特征在于，所述步骤S101，包括：
确定所述网格边界的流出方向，在预设来流条件下确定附面层外空气速度；
计算所述附面层外空气速度在所述流出方向的速度投影，并将所述速度投影确定为所述溢流速度。
4.根据权利要求1所述的飞机结冰冰形预测方法，其特征在于，所述控制体的所述质量守恒方程为：
；
式中，为所述撞击水质量流量，/>为所述流入水质量流量，/>为所述蒸发水蒸发速率，/>为所述流出水质量流量，/>为所述冻结水冻结速率，i为相邻网格的序号；
所述控制体的所述能量守恒方程为：
；
式中，为液态水冻结离开所述控制体的能量，/>为液态水蒸发离开所述控制体的能量，/>为与所述控制体相碰撞的水滴的能量，/>为气流摩擦对控制体表面的加热量，为气流与所述控制体表面的对流换热量，/>为流入所述控制体的溢流水能量，为流出所述控制体的所述溢流水能量。
5.根据权利要求1所述的飞机结冰冰形预测方法，其特征在于，所述根据所述计算冻结水冻结速率对冻结水冻结速率进行赋值，包括：
将所述计算冻结水冻结速率代入第三模型对所述冻结水冻结速率进行赋值，所述第三模型的表达式为：
；
式中，为所述撞击水质量流量，/>为所述流入水质量流量，/>为所述蒸发水蒸发速率，/>为所述流出水质量流量，/>为所述冻结水冻结速率，/>为所述计算冻结水冻结速率。
6.根据权利要求1所述的飞机结冰冰形预测方法，其特征在于，所述在预设来流条件下计算空气流场数据，根据所述空气流场数据计算所述控制体的撞击水质量流量，包括：
在所述预设来流条件下求解不可压N-S方程得到的外形的压力和附外层空气速度；
应用Euler方法根据所述压力和所述附外层空气速度求解水滴的输运方程，计算所述控制体的所述撞击水质量流量。
7.一种飞机结冰冰形预测装置，其特征在于，包括：
数据确定模块，用于确定溢流水流动时间步、控制体的活动水质量流量和网格边界的溢流速度；
第一计算模块，用于根据所述溢流速度、所述溢流水流动时间步和所述活动水质量流量计算所述控制体的流出水质量流量；
第二计算模块，用于将相邻控制体的所述流出水质量流量确定为所述控制体的流入水质量流量；
第三计算模块，用于根据所述流出水质量流量和所述流入水质量流量计算迭代活动水质量流量；
判断模块，用于判断所述迭代活动水质量是否收敛；若否，执行第一执行模块；若是，执行第二执行模块；
第一执行模块，用于将所述迭代活动水质量确定为所述活动水质量流量，并重新开始执行所述第一计算模块；
第二执行模块，用于根据所述迭代活动水质量流量确定目标冻结速率，并根据所述目标冻结速率计算飞机预测冰形；
其中，所述第一计算模块，包括：
第一模型单元，用于将所述溢流速度、所述溢流水流动时间步和所述活动水质量流量输入第一模型计算网格的每个所述网格边界的所述流出水质量流量；
求和单元，用于将全部所述网格边界的所述流出水质量流量相加得到所述控制体的所述流出水质量流量，所述第一模型的表达式为：
；
式中，为第i个网格边界长度，/>为网格面积，/>为所述溢流水流动时间步，/>为第i个所述网格边界的所述溢流速度，/>为所述活动水质量流量，/>为第i个所述网格边界的所述流出水质量流量；
所述第三计算模块，包括：
流场计算单元，用于在预设来流条件下计算空气流场数据，根据所述空气流场数据计算所述控制体的撞击水质量流量；
冻结速率单元，用于根据表面温度确定蒸发水蒸发速率，并联立求解所述控制体的质量守恒方程与能量守恒方程，得到计算冻结水冻结速率；
更新单元，用于根据所述计算冻结水冻结速率对冻结水冻结速率进行赋值，并根据所述冻结水冻结速率更新所述表面温度和所述蒸发水蒸发速率；
计算单元，用于根据所述蒸发水蒸发速率、所述冻结水冻结速率、所述撞击水质量流量、所述流出水质量流量和所述流入水质量流量确定所述迭代活动水质量流量。
8.一种飞机结冰冰形预测设备，其特征在于，包括：
存储器，用于储存计算机程序；
处理器，用于执行所述计算机程序时实现如权利要求1至6任一项所述飞机结冰冰形预测方法。
9.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机可执行指令，所述计算机可执行指令被处理器执行时，实现如权利要求1至6任一项所述飞机结冰冰形预测方法。
CN202310824277.5A 2023-07-06 2023-07-06 一种飞机结冰冰形预测方法、装置、设备及存储介质 Active CN116562192B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310824277.5A CN116562192B (zh) 2023-07-06 2023-07-06 一种飞机结冰冰形预测方法、装置、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310824277.5A CN116562192B (zh) 2023-07-06 2023-07-06 一种飞机结冰冰形预测方法、装置、设备及存储介质

Publications (2)

Publication Number Publication Date
CN116562192A CN116562192A (zh) 2023-08-08
CN116562192B true CN116562192B (zh) 2023-09-12

Family

ID=87503822

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310824277.5A Active CN116562192B (zh) 2023-07-06 2023-07-06 一种飞机结冰冰形预测方法、装置、设备及存储介质

Country Status (1)

Country Link
CN (1) CN116562192B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117875224B (zh) * 2024-03-12 2024-06-04 南京航空航天大学 一种非等分连续多粒径水滴耦合欧拉方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2996844A1 (en) * 2014-08-29 2016-06-16 Tzunum, Inc. Power train for a hybrid-electric aircraft
CN108460217A (zh) * 2018-03-13 2018-08-28 西北工业大学 一种非稳态三维结冰数值模拟方法
CN111291505A (zh) * 2020-05-08 2020-06-16 中国空气动力研究与发展中心低速空气动力研究所 一种基于深度置信网络的翼型结冰冰形预测方法及装置
CN111396269A (zh) * 2020-06-08 2020-07-10 中国空气动力研究与发展中心低速空气动力研究所 一种多时间步非定常结冰计算方法、系统及存储介质
CN114398844A (zh) * 2022-01-25 2022-04-26 南京航空航天大学 基于连续水膜流动的稳态防冰仿真方法
CN114462330A (zh) * 2022-01-19 2022-05-10 清华大学 飞机结冰冰形预测方法、装置、计算机设备和存储介质
CN114970308A (zh) * 2021-12-30 2022-08-30 成都流体动力创新中心 一种飞机结冰预测方法、系统及计算机程序产品
CN115081121A (zh) * 2022-08-22 2022-09-20 四川大学 一种考虑溪流现象的结冰模拟方法和存储介质
CN115238611A (zh) * 2022-08-09 2022-10-25 南京航空航天大学 一种基于多相格子玻尔兹曼通量方法的多相流模拟效率优化方法
CN115563906A (zh) * 2022-11-10 2023-01-03 南京航空航天大学 一种基于非定常欧拉两相流的多步长结冰计算方法和系统
CN115659517A (zh) * 2022-11-10 2023-01-31 南京航空航天大学 一种旋翼桨叶结冰准非定常数值模拟方法和系统

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2996844A1 (en) * 2014-08-29 2016-06-16 Tzunum, Inc. Power train for a hybrid-electric aircraft
CN108460217A (zh) * 2018-03-13 2018-08-28 西北工业大学 一种非稳态三维结冰数值模拟方法
CN111291505A (zh) * 2020-05-08 2020-06-16 中国空气动力研究与发展中心低速空气动力研究所 一种基于深度置信网络的翼型结冰冰形预测方法及装置
CN111396269A (zh) * 2020-06-08 2020-07-10 中国空气动力研究与发展中心低速空气动力研究所 一种多时间步非定常结冰计算方法、系统及存储介质
CN114970308A (zh) * 2021-12-30 2022-08-30 成都流体动力创新中心 一种飞机结冰预测方法、系统及计算机程序产品
CN114462330A (zh) * 2022-01-19 2022-05-10 清华大学 飞机结冰冰形预测方法、装置、计算机设备和存储介质
CN114398844A (zh) * 2022-01-25 2022-04-26 南京航空航天大学 基于连续水膜流动的稳态防冰仿真方法
CN115238611A (zh) * 2022-08-09 2022-10-25 南京航空航天大学 一种基于多相格子玻尔兹曼通量方法的多相流模拟效率优化方法
CN115081121A (zh) * 2022-08-22 2022-09-20 四川大学 一种考虑溪流现象的结冰模拟方法和存储介质
CN115563906A (zh) * 2022-11-10 2023-01-03 南京航空航天大学 一种基于非定常欧拉两相流的多步长结冰计算方法和系统
CN115659517A (zh) * 2022-11-10 2023-01-31 南京航空航天大学 一种旋翼桨叶结冰准非定常数值模拟方法和系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
飞机三维结冰模型及其数值求解方法;易贤 等;航空学报;第31卷(第11期);2152-2158 *

Also Published As

Publication number Publication date
CN116562192A (zh) 2023-08-08

Similar Documents

Publication Publication Date Title
CN116562192B (zh) 一种飞机结冰冰形预测方法、装置、设备及存储介质
CN112800700B (zh) 低温表面干模态结霜模拟方法、装置、电子设备和介质
CN107491608B (zh) 一种飞机编队飞行的队形参数优化方法及系统
CN105975645B (zh) 一种基于多步的含激波区域飞行器流场快速计算方法
CN109376403B (zh) 一种基于笛卡尔自适应重构技术的二维结冰数值模拟方法
CN108460217B (zh) 一种非稳态三维结冰数值模拟方法
CN105718634B (zh) 一种基于非概率区间分析模型的翼型鲁棒优化设计方法
CN112818573B (zh) 一种用于非结构网格的获取边界层非当地变量信息的方法
CN116384290B (zh) 一种考虑真实气体效应的高超声速飞行器动导数预测方法
CN113777931B (zh) 结冰翼型气动模型构造方法、装置、设备及介质
Moghtadernejad et al. SPH simulation of rivulet dynamics on surfaces with various wettabilities
CN114444216B (zh) 基于数值模拟的高空条件下飞行器姿态控制方法及系统
CN111310381A (zh) 一种三维水滴收集系数计算方法
CN115817822B (zh) 电加热防冰系统的热载荷分布设计方法
Kollar et al. Inverse design of wind turbine blade sections for operation under icing conditions
KR20140115835A (ko) 공기 흡입구의 표면 결빙 시뮬레이션 방법
CN115618498A (zh) 一种飞行器跨流域流场的预测方法、装置、设备及介质
Kashani et al. Heat and mass transfer in natural draft cooling towers
CN117408189B (zh) 高超声速边界层的转捩预测方法、装置、设备及存储介质
Xudaykulov et al. Turbulent flow out of a convex curve in a channel using the SST turbulence model
CN116738576B (zh) 一种旋翼结冰冰形预测方法、装置、设备及存储介质
Zhang et al. Lateral motion of a droplet impacting on a wettability-patterned surface: numerical and theoretical studies
WO2019186151A1 (en) Methods and apparatus for simulating liquid collection on aerodynamic components
CN116013443A (zh) 一种传热特性预测方法、装置、设备及可读存储介质
Szilder et al. Numerical simulation of ice accretion influenced by anti-icing heating system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant