CN108304600B - 一种高超声速飞行器转捩位置预测方法 - Google Patents
一种高超声速飞行器转捩位置预测方法 Download PDFInfo
- Publication number
- CN108304600B CN108304600B CN201710673355.0A CN201710673355A CN108304600B CN 108304600 B CN108304600 B CN 108304600B CN 201710673355 A CN201710673355 A CN 201710673355A CN 108304600 B CN108304600 B CN 108304600B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- transition
- flow
- hypersonic aircraft
- transition position
- 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
Links
Images
Classifications
-
- 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
-
- 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)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供一种高超声速飞行器转捩位置预测方法,先计算得到飞行条件下飞行器表面的层流流场,并将流场插值到沿飞行器表面法向的正交化网格上;沿正交化网格下的流场的势流方向对第二模态和更低频的扰动波的增长率进行积分得到飞行器表面的N值分布;根据触发转捩的N0值确定飞行器表面转捩的位置。该方法从横流稳定性的角度去分析和预测转捩的发生,改进了传统e‑N方法,考虑了横流稳定性,对横流模态的扰动波进行了积分,同时考虑了飞行器壁面温度对转捩的影响,获得了一种基于稳定性分析的转捩预测方法,可以较为准确的预测高超声速飞行器表面转捩位置。
Description
技术领域
本发明涉及飞行器转捩预测技术领域,具体涉及一种高超声速飞行器边界层的转捩位置预测方法,尤其涉及一种基于边界层稳定性理论的应用于升力体/乘波体高超声速飞行器边界层的转捩预测方法。
背景技术
高超声速飞行器(下文简称为飞行器)在飞行试验过程中,飞行器表面存在层流和湍流两种流态。流态的不同导致飞行器升力、阻力、表面热流、发动机性能和进气道启动性能的不同。本发明中定义从层流到湍流的流动过程为转捩。准确的预测飞行器表面的转捩位置是飞行器设计的关键,决定飞行试验成败。
预测飞行器表面的转捩位置,目前常采用的方法有:基于线性稳定性理论的半经验的传统e-N方法、求解抛物化稳定性方程的方法(PSE)、直接数值模拟(DNS)、湍流转捩模式和工程转捩准则,其中,PSE方法无法应用于工程实际;DNS和湍流转捩模式计算量巨大,无法应用于工程实际;工程转捩准则没有理论依据,无法模拟转捩过程,给出精确的转捩位置。因此,在工程应用中,基于稳定性理论的e-N方法应用最为广泛。在航空工业领域,e-N方法通常被认为是预测转捩位置最有效的方法。
传统e-N方法:边界层内存在各种频率的小扰动,当它们向下游传播时,分别在不同的位置幅值开始增长;各种频率的扰动波从其开始增长的位置起,可用线性稳定性理论计算它们向下游方向累计的增长倍数,增长倍数首先达到某一特定值的扰动波引发转捩。具体表示为:
其中,A为扰动的幅值,A0为扰动开始增长时的幅值。
通常情况下,e-N方法中的增长率通过使用线性稳定性理论计算得到,其中引入了平行流假设,即假设基本流沿流向的展向是平行的。这种假设的理论依据是基本流在流向和展向的尺度远大于扰动波在流向和展向上的长度尺度,因此,大多数情况下,这个假设是合理的。但对于复杂外形高超声速飞行器的流场,横流明显、三维性较强,流向和展向变化比较大,因此使用线性稳定性给出的增长率通常存在较大的误差,可靠性较差。而由于存在横流,仅仅使用流向不稳定性(传统e-N方法)来分析和预测复杂外形飞行器是不合适的。
发明内容
本发明的目的在于克服现有技术不足,提供了一种高超声速飞行器转捩位置预测方法,该预测方法基于一种改进的e-N方法,考虑了由于外形复杂引起的三维横流的影响,改进了传统e-N方法的不足,适用于复杂外形高超声速飞行器转捩位置的预测。
本发明的技术解决方案:
一种高超声速飞行器转捩位置预测方法,通过以下步骤实现:
计算目标飞行器在典型状态工况下飞行器表面的层流流场,其中,壁面温度根据沿弹道气动热和飞行器结构计算得到;
该步骤中,层流流场的计算为本领域公知的技术,可采用商业软件或工具获得;
将计算得到的层流流场插值到飞行器表面沿法向的正交网格上;
以正交化的流场作为基本流,分析中性曲线,得到考虑横流影响的中性曲线;
基于得到的中性曲线,选取幅值增长率大于0的第二模态和横流模态扰动波,采用e-N方法对所述的扰动波沿势流方向积分,并取所计算的所有扰动波的N值的包络线作为飞行器表面的N值分布;
根据经过试验验证或经验得到的触发转捩的N0值,确定飞行器表面的转捩位置。
本发明实施例提供的一种基于改进的e-N稳定性分析的复杂外形高超声速飞行器转捩预测方法,先计算得到飞行条件下飞行器表面的层流流场,并将流场插值到沿飞行器表面法向的正交化网格上;沿正交化网格下的流场的势流方向对第二模态和更低频的扰动波的增长率进行积分得到飞行器表面的N值分布;根据触发转捩的N0值确定飞行器表面转捩的位置。该方法从横流稳定性的角度去分析和预测转捩的发生,改进了传统e-N方法,考虑了横流稳定性,对横流模态的扰动波进行了积分,同时考虑了飞行器壁面温度对转捩的影响,获得了一种基于稳定性分析的转捩预测方法,可以较为准确的预测高超声速飞行器表面转捩位置。
附图说明
所包括的附图用来提供对本发明实施例的进一步的理解,其构成了说明书的一部分,用于例示本发明的实施例,并与文字描述一起来阐释本发明的原理。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的高超声速飞行器转捩位置预测方法的一种实施例的流程图;
图2为本发明提供的高超声速飞行器转捩位置预测方法中考虑和不考虑横流的中性曲线对比图;
图3为本发明提供的高超声速飞行器转捩位置预测方法中不同频率扰动波沿势流方向的N值分布;
图4为本发明提供的高超声速飞行器转捩位置预测方法中N值分布示意图;
图5为采用本发明的方法得到的转捩位置示意图;
图6为试验结果出示的转捩位置示意图。
其中:图5和6中黑线表示转捩位置,N0值=10。
具体实施方式
下面将结合附图对本发明的具体实施例进行详细说明。在下面的描述中,出于解释而非限制性的目的,阐述了具体细节,以帮助全面地理解本发明。然而,对本领域技术人员来说显而易见的是,也可以在脱离了这些具体细节的其它实施例中实践本发明。
在此需要说明的是,为了避免因不必要的细节而模糊了本发明,在附图中仅仅示出了与根据本发明的方案密切相关的设备结构和/或处理步骤,而省略了与本发明关系不大的其他细节。
本发明实施例提供了一种高超声速飞行器转捩位置预测方法,如图1所示,包括下述步骤:
1、弹道分析;
针对目标飞行器的飞行弹道,通过分析统计弹道的各个参数的变化范围,获得马赫数、高度、攻角、壁面温度的变化范围,得到建立数据库所需要的样本,此步骤为本领域公知的步骤;
2、建立N值分布数据库,包括:
2.1获得飞行器样本工况的层流流场;
根据步骤1所获得的样本工况,采用工程计算软件Fluent计算得到飞行器样本工况的层流流场,计算过程中,以飞行弹道条件下的气动热环境和真实的飞行器结构为基础,通过结构传热分析得到飞行器表面的温度,以此温度作为Fluent计算的壁面条件;
其中,通过结构传热分析得到飞行器表面的温度可采用本领域常规技术手段获得,具体见《数值传热学》,杨世铭,高等教育出版社;
2.2步骤2.1得到的层流流场转换到飞行器表面沿法向的正交网格上,得到正交化的层流流场;
2.3以得到的正交化的层流流场作为基本流,采用线性稳定性理论分析中性曲线特征,得到考虑横流后的中性曲线,并将该中性曲线中受横流影响最大的第一模态作为横流模态;
该步骤中,采用线性稳定性理论分析中性曲线特征可采用本领域公知的手段,具体见《流动稳定性》,周恒等,国防工业出版社;
举例说明,参见图2,通常情况下,采用线性稳定性理论分析中性曲线时,基本流中只考虑x和y两个方向的分量,强制z方向分量为0;但在考虑横流后,基本流中z方向的分量不再强制为0。如图2所示,给出了不考虑横流和考虑横流时中性曲线的比较。图中横坐标为扰动波展向波数β,纵坐标为扰动波频率f,云图颜色的变化表示扰动波的幅值增长率(-αi)的大小,满足-αi=0条件的点的集合称为中性曲线。从图中可以看出,相比与第二模态(Mack模态),低频扰动模态受横流的影响很大,横流使得第一模态失稳区域扩大并且向更低频区域移动。可以将此时受横流影响大的第一模态,看作为横流模态。
2.4根据步骤一得到的中性曲线,选取幅值增长率-αi>0的第二模态扰动波和更低频的横流模态波,采用e-N方法沿势流方向积分,得到不同频率扰动波沿势流方向的N值分布,如图3所示;
该步骤中,所述N值的具体计算公式为e-N方法中的常规公式,如下所示:
其中,A为扰动的幅值,A0为扰动开始增长时的幅值,一般取0.0001,x0表示某一频率的扰动波使得增长率αi=0的流向位置,x表示流向位置;
2.5统计步骤2.4得到的沿势流方向的N值分布,得到沿势流方向的N值的包络线(不同x位置处最大N值的集合),形成飞行器表面的N值分布,如图4所示;
步骤3转捩位置判定;
基于所述的飞行器表面的N值分布,结合根据经过试验验证或经验得到的触发转捩的N0值=10;确定飞行器表面的转捩位置,见图5中黑色曲线表示N0=10的位置,即飞行器转捩位置。
其中,对飞行器的任意一条弹道,可采用插值获得飞行器表面沿弹道的N值。
采用考虑横流模态的改进e-N方法,对目标飞行器进行了边界层转捩预测。结果显示流动三维性强的区域提前转捩发展为湍流,考虑横流模态改进后的e-N方法显示该部位发生转捩,并与试验结果转捩位置(图6)比较,说明了考虑横流模态的改进e-N方法的有效性,而传统e-N方法未能预测出转捩。
如上针对一种实施例描述和/或示出的特征可以以相同或类似的方式在一个或更多个其它实施例中使用,和/或与其它实施例中的特征相结合或替代其它实施例中的特征使用。
应该强调,术语“包括/包含”在本文使用时指特征、整件、步骤或组件的存在,但并不排除一个或更多个其它特征、整件、步骤、组件或其组合的存在或附加。
本发明以上的装置和方法可以由硬件实现,也可以由硬件结合软件实现。本发明涉及这样的计算机可读程序,当该程序被逻辑部件所执行时,能够使该逻辑部件实现上文所述的装置或构成部件,或使该逻辑部件实现上文所述的各种方法或步骤。本发明还涉及用于存储以上程序的存储介质,如硬盘、磁盘、光盘、DVD、flash存储器等。
这些实施例的许多特征和优点根据该详细描述是清楚的,因此所附权利要求旨在覆盖这些实施例的落入其真实精神和范围内的所有这些特征和优点。此外,由于本领域的技术人员容易想到很多修改和改变,因此不是要将本发明的实施例限于所例示和描述的精确结构和操作,而是可以涵盖落入其范围内的所有合适修改和等同物。
本发明未详细说明部分为本领域技术人员公知技术。
Claims (5)
1.一种高超声速飞行器转捩位置预测方法,其特征在于,通过以下步骤实现:
计算目标飞行器在典型状态工况下飞行器表面的层流流场;
将计算得到的层流流场插值到飞行器表面沿法向的正交网格上;
以正交化的流场作为基本流,分析中性曲线,得到考虑横流影响的中性曲线,其中,将该中性曲线中受横流影响最大的第一模态作为横流模态;
基于得到的中性曲线,选取幅值增长率大于0的第二模态和横流模态扰动波,采用e-N方法对所述的扰动波沿势流方向积分,并取所计算的所有扰动波的N值的包络线作为飞行器表面的N值分布;
根据触发转捩的N0值,确定飞行器表面的转捩位置。
2.根据权利要求1所述的一种高超声速飞行器转捩位置预测方法,其特征在于,所述计算目标飞行器在典型状态工况下飞行器表面的层流流场时,所需要的壁面温度根据沿弹道气动热和飞行器结构计算得到。
3.根据权利要求1所述的一种高超声速飞行器转捩位置预测方法,其特征在于,所述的N0值根据经过试验验证或经验得到。
4.根据权利要求3所述的一种高超声速飞行器转捩位置预测方法,其特征在于,所述的N0值为8~10。
5.根据权利要求1-4中任一项权利要求所述的一种高超声速飞行器转捩位置预测方法,其特征在于,所述的高超声速飞行器为升力体/乘波体外形的高超声速飞行器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710673355.0A CN108304600B (zh) | 2017-08-09 | 2017-08-09 | 一种高超声速飞行器转捩位置预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710673355.0A CN108304600B (zh) | 2017-08-09 | 2017-08-09 | 一种高超声速飞行器转捩位置预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108304600A CN108304600A (zh) | 2018-07-20 |
CN108304600B true CN108304600B (zh) | 2020-03-31 |
Family
ID=62872634
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710673355.0A Active CN108304600B (zh) | 2017-08-09 | 2017-08-09 | 一种高超声速飞行器转捩位置预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108304600B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108304601B (zh) * | 2017-08-09 | 2019-12-31 | 北京空天技术研究所 | 一种高超声速飞行器边界层转捩的判断方法 |
CN109583067B (zh) * | 2018-11-22 | 2020-04-07 | 北京空天技术研究所 | 基于温度平衡的高速飞行器转捩位置测量传感器设计方法 |
CN112231847B (zh) * | 2020-11-04 | 2024-04-02 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 转捩位置确定方法、装置、电子设备及存储介质 |
CN112613250B (zh) * | 2020-12-29 | 2021-12-10 | 中国航天空气动力技术研究院 | 火星进入器表面流动转捩位置预测方法 |
CN112861447B (zh) * | 2021-02-09 | 2022-03-04 | 天津大学 | 一种基于稳定性理论的水下回转体首部线型优选设计方法 |
CN113221350B (zh) * | 2021-05-10 | 2022-02-18 | 天津大学 | 基于全局稳定性分析的高超声速飞行器转捩预测方法 |
CN113486440B (zh) * | 2021-05-25 | 2023-07-14 | 北京临近空间飞行器系统工程研究所 | 基于高频压力传感器测量高速边界层扰动波的布置方法 |
CN113468655B (zh) * | 2021-05-25 | 2023-05-09 | 北京临近空间飞行器系统工程研究所 | 基于飞行试验数据的高速边界层转捩判据标定方法和系统 |
CN113505542B (zh) * | 2021-06-11 | 2023-03-21 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种高超声速飞行器后掠机翼的湍流转捩模型构建方法 |
CN114061889A (zh) * | 2021-12-14 | 2022-02-18 | 南京航空航天大学 | 一种用于第二模态波直接显示的超声同频观测装置及方法 |
CN114117648B (zh) * | 2022-01-24 | 2022-04-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 同时抑制Mack模态和横流失稳的组合被动式控制结构 |
CN114528665B (zh) * | 2022-02-24 | 2022-12-09 | 天津大学 | 一种壁面加热/冷却的水下平板边界层转捩位置的预测方法 |
CN115293070B (zh) * | 2022-09-28 | 2022-12-23 | 中国人民解放军国防科技大学 | 基于一般时间根方尺度的四方程转捩模型的数值方法 |
CN115659522B (zh) * | 2022-12-27 | 2023-03-28 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种飞行器转捩位置预测方法、装置、设备及介质 |
CN115809513B (zh) * | 2023-02-08 | 2023-05-26 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种强迫转捩-俯仰震荡数值模拟方法 |
CN115859482B (zh) * | 2023-02-16 | 2023-04-28 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种飞行器流场稳态基本流快速计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101148195A (zh) * | 2006-09-18 | 2008-03-26 | 刘维国 | 新型翼刀及其强化辅具——翼梢喷气翼 |
CN103954425A (zh) * | 2014-04-30 | 2014-07-30 | 北京大学 | 高超声速静风洞喷管设计方法及该喷管转捩位置确定方法 |
CN106021665A (zh) * | 2016-05-10 | 2016-10-12 | 中国航天空气动力技术研究院 | 高超声速飞行器多批次气动数据快速预测方法 |
CN106547954A (zh) * | 2016-10-17 | 2017-03-29 | 北京航空航天大学 | 一种低雷诺数错层翼的翼型优化设计方法 |
CN106596033A (zh) * | 2016-11-04 | 2017-04-26 | 中国航天空气动力技术研究院 | 一种临近空间螺旋桨地面试验测试系统及测试方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8538738B2 (en) * | 2011-03-22 | 2013-09-17 | Aerion Corporation | Predicting transition from laminar to turbulent flow over a surface |
CN105258904B (zh) * | 2015-11-02 | 2017-12-22 | 中国航天空气动力技术研究院 | 一种旋转导弹锥形运动稳定性风洞试验方法 |
CN106124626B (zh) * | 2016-06-30 | 2018-11-20 | 中国航天空气动力技术研究院 | 延迟高超声速边界层转捩材料的研究系统及实验研究方法 |
-
2017
- 2017-08-09 CN CN201710673355.0A patent/CN108304600B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101148195A (zh) * | 2006-09-18 | 2008-03-26 | 刘维国 | 新型翼刀及其强化辅具——翼梢喷气翼 |
CN103954425A (zh) * | 2014-04-30 | 2014-07-30 | 北京大学 | 高超声速静风洞喷管设计方法及该喷管转捩位置确定方法 |
CN106021665A (zh) * | 2016-05-10 | 2016-10-12 | 中国航天空气动力技术研究院 | 高超声速飞行器多批次气动数据快速预测方法 |
CN106547954A (zh) * | 2016-10-17 | 2017-03-29 | 北京航空航天大学 | 一种低雷诺数错层翼的翼型优化设计方法 |
CN106596033A (zh) * | 2016-11-04 | 2017-04-26 | 中国航天空气动力技术研究院 | 一种临近空间螺旋桨地面试验测试系统及测试方法 |
Non-Patent Citations (2)
Title |
---|
高超声速三维边界层e_N方法的应用及扰动演化的预测;于高通;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20170715;全文 * |
高超音速圆锥边界层的转捩预测及e-N方法的改进;苏彩虹;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20090715;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108304600A (zh) | 2018-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108304600B (zh) | 一种高超声速飞行器转捩位置预测方法 | |
CN108304601B (zh) | 一种高超声速飞行器边界层转捩的判断方法 | |
JP5747343B2 (ja) | 超音速航空機の実機相当レイノルズ数における自然層流翼の設計方法 | |
CN113221350B (zh) | 基于全局稳定性分析的高超声速飞行器转捩预测方法 | |
Kostić et al. | CFD modeling of supersonic airflow generated by 2D nozzle with and without an obstacle at the exit section | |
CN113657051B (zh) | 一种水下航行体表面层流区壁面脉动压力频谱预测方法 | |
CN105631125A (zh) | 基于降阶模型的气动-热-结构耦合分析方法 | |
Chwalowski et al. | Fun3d analyses in support of the second aeroelastic prediction workshop | |
Aftosmis et al. | Applications of a Cartesian mesh boundary-layer approach for complex configurations | |
Goyal et al. | Dynamic effects in transition from regular to Mach reflection in steady supersonic flows | |
CN115879216A (zh) | 一种内流道强波系干扰控制下的流场重构设计方法 | |
Yamashita et al. | Full-field sonic boom simulation in real atmosphere | |
Hadidoolabi et al. | Computational investigation of vortex structure and breakdown over a delta wing at supersonic pitching maneuver | |
Li et al. | The simulation of wraparound fins’ aerodynamic characteristics | |
Bui | CFD Analysis of the Nozzle Jet Plume Effects on Sonic Boom Signature | |
Katzenmeier et al. | Correction technique for quality improvement of doublet lattice unsteady loads by introducing CFD small disturbance aerodynamics | |
Hindman et al. | Q3D-wing: validation of a modern lifting-line method for application to clean and iced wings | |
Matsushima et al. | Investigation into far-field drag calculation methods and pseudo total enthalpy generation for airplane CFD | |
Serdarevic-Kadic et al. | Effects of Base Shape to Drag at Transonic and Supersonic Speeds by CFD | |
Kuzmin et al. | Transonic aerofoils admitting anomalous behaviour of lift coefficient | |
Karpuzcu et al. | On Linear Stability of Compression Corner Flows Obtained by Kinetic Theory | |
Ullah et al. | Effect of Reduced Frequency and Mach Number on Damping Characteristics of SACCON UCAV | |
Bui | Computational Fluid Dynamics Analysis of Nozzle Plume Effects on Sonic Boom Signature | |
Hoshi et al. | Numerical simulations of transonic flutter on a three-dimensional wing | |
Carter et al. | USM3D Predictions of Supersonic Nozzle Flow |
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 |