CN111553119B - 非均匀来流下高超声速飞行器三维曲面激波反设计方法 - Google Patents
非均匀来流下高超声速飞行器三维曲面激波反设计方法 Download PDFInfo
- Publication number
- CN111553119B CN111553119B CN202010210844.4A CN202010210844A CN111553119B CN 111553119 B CN111553119 B CN 111553119B CN 202010210844 A CN202010210844 A CN 202010210844A CN 111553119 B CN111553119 B CN 111553119B
- Authority
- CN
- China
- Prior art keywords
- shock wave
- dimensional
- infinitesimal
- curved surface
- flow
- 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
- 230000035939 shock Effects 0.000 title claims abstract description 97
- 238000013461 design Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 239000002243 precursor Substances 0.000 description 9
- 230000001419 dependent effect Effects 0.000 description 6
- 230000008878 coupling Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 5
- 238000005859 coupling reaction Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
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
- 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
- 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
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Other Liquid Machine Or Engine Such As Wave Power Use (AREA)
Abstract
本发明公开了一种非均匀来流条件下的三维曲面激波反设计方法,将其命名为微元密切轴对称流场求解方法,所述方法沿目标激波面的周向和流向构建一系列微元密切面,在每个微元面内进行三维向二维流动的等效转换,然后将局部各微元面内的二维求解结果返回至三维空间,最终实现任意非均匀来流条件下的三维曲面激波反设计。本发明为飞行器前体、进气道及二者的气动一体化设计提供了全新的解决思路,在吸气式高超声速推进领域具有重要应用前景。
Description
技术领域
本发明属于高超声速飞行器设计技术领域,具体涉及一种非均匀来流下高超声速飞行器三维曲面激波反设计方法。
背景技术
在过去的几十年间,高超声速技术取得了极大发展,并不断推动着高超声速飞行器的相关研究与研制工作,高超声速飞行器前体、进气道及二者的气动一体化设计是其中的重要研究方向。
乘波理论因其显著的高升阻比和高效预压缩特性,长期以来被公认是高超声速飞行器设计领域的一种理想概念。乘波理论中,前缘激波与飞行器边缘在精心设计下能够完美贴合,飞行器犹如“乘坐”在激波面之上,故而得名。在利用乘波概念进行高超声速飞行器前体、进气道设计时,首先需要进行无粘基准流场构建,其中关键在于对前缘激波的设计和控制,然后从激波面上选定的前缘线进行流线追踪,即可得到具备乘波特性的气动压缩面。
在乘波设计理念的指导下,很自然地逐渐发展形成了两种设计思路:(1)给定激波生成体,通过数值计算得到激波面和波后流场;(2)给定激波形状,通过求解反问题得到压缩面和波后流场。两种思路均能得到满足要求的基准流场,但各有长短。
第一种思路在设计时能够得到更多复杂构型下的超声速基准流场,无需拘泥于采用二维平面激波、二维轴对称激波等简单几何形状,但其生成的激波可控性较差。第二种思路作为一种反向设计,能够对所乘激波形状进行精细控制,但实现难度较大,尤其是三维空间中,到目前为止仅在某些特殊约束下(如激波形状为三维密切轴对称曲面,来流条件为均匀水平气流等)实现了对曲面激波反问题的求解,更多一般意义上的三维曲面激波反设计问题仍处于初步的探索阶段,造成高超声速飞行器前体与进气道的气动一体化设计灵活性受限。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种非均匀来流下高超声速飞行器三维曲面激波反设计方法,方便高超声速飞行器前体与进气道在流向布局下的气动一体化设计。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种用于求解非均匀来流下三维曲面激波反问题的微元密切轴对称方法,所述方法沿目标激波面(即三维激波曲面)的周向和流向构建一系列微元密切面,在每个微元密切面内进行三维向二维流动的等效转换,然后将局部各微元密切面内的二维求解结果返回至三维空间,最终实现任意非均匀来流条件下的三维曲面激波反设计。
所述沿目标激波面周向和流向构建一系列微元密切面的方法为:在整个三维激波曲面周向和流向的每一个离散点上,将波前和波后速度方向确定的微元密切面作为该激波点附近的微元密切面。
所述在每个微元面内进行三维向二维流动等效转换的方法为:首先计算周向相邻两个微元密切平面的交线,以此交线为对称轴,将该微元密切面内的流动等效为绕该微元对称轴的轴对称二维流动;然后运用本领域的成熟技术二维特征线法,沿流向各微元密切面依次求解,得到微元密切平面内的激波波后二维流场和激波依赖区二维壁面。
所述将局部各微元面内的二维求解结果返回值三维空间的方法为:将每个微元密切平面内的激波波后二维流场和激波依赖区二维壁面按照各微元密切面与原三维空间的位置关系进行叠加变换,最终形成空间中三维曲面激波的波后流场与激波依赖区壁面。
具体步骤包括以下流程:
(1)根据Rankine-Hugoniot方程,由非均匀来流条件和三维曲面激波形状计算激波波后参数;
(2)按照激波前后速度方向,在三维曲面激波面周向和流向各激波点附近构造微元密切面;
(3)计算周向相邻两个微元密切面的交线,作为微元密切面上轴对称二维流动的微元对称轴;
(4)利用轴对称二维特征线法依次求解沿流向各微元密切轴对称面上的波后流动;
(5)若所有微元密切面上的解均存在,分别将各微元密切面上的轴对称二维流场和二维流线按照空间位置关系进行叠加,组成空间三维流场和三维曲面激波的依赖区压缩面;否则,说明目标三维曲面激波在给定的非均匀来流条件下不存在,调整输入条件,重复上述过程。
进一步的,步骤(1)中,所述激波波后参数包括气流压力、密度、速度大小和方向。
进一步的,步骤(2)中,在整个三维激波曲面周向和流向的每一个离散点上,将波前和波后速度方向确定的微元面作为过该激波点的微元密切面。此处意为该点前、后两个速度方向(即两个交于该点的直线)依据两条相交直线确定一个平面的理论,确定出经过该点的微元密切面。
进一步的,步骤(5)中,所述调整输入条件包括三维曲面激波形状或非均匀来流参数。
有益效果:本发明提供的非均匀来流下高超声速飞行器三维曲面激波反设计方法,能够在非均匀来流下,通过预先给定所需的三维曲面激波形状,得到该三维曲面激波的波后流场与激波依赖区壁面。该方法显著增强了乘波体、进气道所需基准流场的设计灵活性,能够为高超声速飞行器的前体、进气道及二者在流向布局下的一体化设计提供解决途径。
附图说明
图1为微元密切轴对称方法流程图。
图2为微元密切轴对称方法示意图。
图3为弯曲流面上的特征线网格。
图4为一体化外锥与内锥耦合基准流场。
图5为应用本发明设计得到的内外锥耦合流场数值模拟结果(三维视图)。
图6为应用本发明设计得到的内外锥耦合流场数值模拟结果(正视图)。
其中:1为表示三维曲面激波的网格线,2为微元对称轴,3为原基准轴,4为一系列流向排布的微元密切面拼接而成的流面,5为三维曲面激波的依赖区壁面,6为微元密切面,7为左行特征线,8为流线,9为外锥流场对称轴,10为内锥激波前缘线,11为内外锥激波曲面交线,12为内锥激波曲面,13为被包裹在外锥激波以内的内锥激波曲面,14为内锥流场对称轴,15为外锥激波曲面,16为外锥压缩面,17为内锥压缩面,18为内锥流场壁面流线,19为外锥流场壁面流线。
具体实施方式
下面结合附图对本发明作更进一步的说明。
本发明公开了一种非均匀来流下高超声速飞行器三维曲面激波反设计方法,可用于前体、进气道及二者的气动一体化设计,如图1所示,该方法具体包括以下流程:
(1)根据Rankine-Hugoniot方程,由非均匀来流条件和三维曲面激波形状计算激波波后参数;
(2)按照激波前后速度方向,在激波面周向和流向各点(如P11、P21、P12、P22等)附近构造微元密切面,如图2所示;
(3)计算周向相邻两个微元密切面的交线,作为微元密切面上轴对称二维流动的微元对称轴,如图2所示;
(4)利用轴对称二维特征线法依次求解沿流向各微元密切轴对称面上的波后流动,如图3所示;
(5)若所有微元密切面上的解均存在,分别将各微元密切面上的轴对称二维流场和二维流线按照空间位置关系进行叠加,最终组成空间三维流场和三维曲面激波的依赖区压缩面;否则,说明目标三维曲面激波在给定的非均匀来流条件下不存在,调整输入条件,重复上述过程。
实施例
下面结合一个实施例对本发明作进一步说明。
如图4所示,该实施例为一种一体化外锥与内锥耦合流场,其中外锥流场为常见的锥形流,外锥压缩面16的半锥角为10°,内锥流场布置在外锥流场之内,轴对称的内锥激波曲面12的母线方程为五次拟合函数y=A0+A1x+A2x2+A3x3+A4x4+A5x5。由于外锥激波曲面15波后是非均匀流动,内锥激波曲面12波前参数非均匀,因此是一种典型的非均匀来流下三维曲面激波的反设计问题。该实施例的具体步骤为:第一步,根据锥形流的控制方程Taylor-Maccoll方程,计算外锥激波曲面15的波后流场。
第二步,确定内外锥流场的相对位置关系,包括内外锥流场进口的流向距离和内锥流场对称轴14和外锥流场对称轴9的横向距离。
第三步,计算一体化外锥与内锥耦合流场中的内外锥激波曲面交线11,同时得到被包裹在外锥激波曲面15以内的内锥激波曲面13。
第四步,根据Rankine-Hugoniot方程,计算被包裹在外锥激波以内的内锥激波曲面13的波后流动参数。
第五步,沿被包裹在外锥激波以内的内锥激波曲面13的波后速度方向,从其前缘线上离散的各个周向点开始,向后推进生成激波曲面上离散的各条流向点集。
第六步,在被包裹在外锥激波以内的内锥激波曲面13上的每一个离散点处,按照激波前后速度方向构建通过该点的微元密切平面。
第七步,计算周向相邻两个微元密切面6的交线,作为该微元密切面上轴对称二维流动的微元对称轴2。
第八步,运用二维特征线法(本领域成熟技术),沿流向在各微元密切面6内部求解二维激波波后流场和激波依赖区壁面。
第九步,将各个微元密切面上的微元轴对称二维流动按照各自空间位置组成整体三维流动,得到目标内被包裹在外锥激波以内的锥激波曲面13的波后流场和内锥压缩面17。
采用计算流体力学(CFD)模拟了该实施例的设计结果,如图5和图6所示,在外锥激波曲面15的波后非均匀流场中,得到了被包裹在外锥激波以内的内锥激波曲面13,外锥流场壁面流线19和内锥流场壁面流线18分别呈现出典型的外转与内转特性。同时,横截面上的内锥激波形状为标准的圆弧形,对称面上的形状与给定的母线形状一致,说明应用本发明提出的设计方法能够在非均匀来流条件下生成与目标形状完全相同的三维曲面激波。进一步地,以高超飞行器前体激波的波后流动作为所述设计方法的非均匀来流条件,以进气道前缘激波作为所述设计方法的三维曲面激波,即可实现高超声速飞行器流向布局的前体/进气道气动一体化设计。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (3)
1.一种非均匀来流条件下高超声速飞行器三维曲面激波反设计方法,其特征在于:包括以下步骤:
(1)根据Rankine-Hugoniot方程,由非均匀来流条件和三维曲面激波形状计算激波波后参数;
(2)按照激波前后速度方向,在三维曲面激波面周向和流向各激波点附近构造微元密切面,具体指的是在整个三维激波曲面周向和流向的每一个离散点上,将波前和波后速度方向确定的微元面作为过该激波点的微元密切面;
(3)计算周向相邻两个微元密切面的交线,作为微元密切面上轴对称二维流动的微元对称轴;
(4)利用轴对称二维特征线法依次求解沿流向各微元密切轴对称面上的波后流动;
(5)若所有微元密切面上的解均存在,分别将各微元密切面上的轴对称二维流场和二维流线按照空间位置关系进行叠加,组成空间三维流场和三维曲面激波的依赖区压缩面;否则,说明目标三维曲面激波在给定的非均匀来流条件下不存在,调整输入条件,重复上述步骤(1)-(5)。
2.根据权利要求1所述的一种非均匀来流条件下高超声速飞行器三维曲面激波反设计方法,其特征在于:步骤(1)中,所述激波波后参数包括气流压力、密度、速度大小和方向。
3.根据权利要求1所述的一种非均匀来流条件下高超声速飞行器三维曲面激波反设计方法,其特征在于:步骤(5)中,所述调整输入条件包括三维曲面激波形状或非均匀来流参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010210844.4A CN111553119B (zh) | 2020-03-24 | 2020-03-24 | 非均匀来流下高超声速飞行器三维曲面激波反设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010210844.4A CN111553119B (zh) | 2020-03-24 | 2020-03-24 | 非均匀来流下高超声速飞行器三维曲面激波反设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111553119A CN111553119A (zh) | 2020-08-18 |
CN111553119B true CN111553119B (zh) | 2024-02-13 |
Family
ID=72005549
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010210844.4A Active CN111553119B (zh) | 2020-03-24 | 2020-03-24 | 非均匀来流下高超声速飞行器三维曲面激波反设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111553119B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101392685A (zh) * | 2008-10-29 | 2009-03-25 | 南京航空航天大学 | 基于任意激波形状的内乘波式高超声速进气道及设计方法 |
CN105173116A (zh) * | 2015-09-25 | 2015-12-23 | 北京航空航天大学 | 高超声速飞行器密切曲面乘波体设计方法 |
CN105697150A (zh) * | 2016-03-07 | 2016-06-22 | 厦门大学 | 一种兼顾锥导与密切乘波的鼓包进气道设计方法 |
CN105947230A (zh) * | 2016-05-24 | 2016-09-21 | 中国人民解放军63820部队吸气式高超声速技术研究中心 | 一种乘波体和进气道一体化构型的设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7207520B2 (en) * | 2005-05-31 | 2007-04-24 | Lockheed Martin Corporation | System, method, and apparatus for designing streamline traced, mixed compression inlets for aircraft engines |
-
2020
- 2020-03-24 CN CN202010210844.4A patent/CN111553119B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101392685A (zh) * | 2008-10-29 | 2009-03-25 | 南京航空航天大学 | 基于任意激波形状的内乘波式高超声速进气道及设计方法 |
CN105173116A (zh) * | 2015-09-25 | 2015-12-23 | 北京航空航天大学 | 高超声速飞行器密切曲面乘波体设计方法 |
CN105697150A (zh) * | 2016-03-07 | 2016-06-22 | 厦门大学 | 一种兼顾锥导与密切乘波的鼓包进气道设计方法 |
CN105947230A (zh) * | 2016-05-24 | 2016-09-21 | 中国人民解放军63820部队吸气式高超声速技术研究中心 | 一种乘波体和进气道一体化构型的设计方法 |
Non-Patent Citations (4)
Title |
---|
三维曲面激波反问题的参考平面解法;赵玉新;蓝庆生;赵一龙;;推进技术(11);全文 * |
基于内乘波概念的三维变截面高超声速进气道;尤延铖;梁德旺;;中国科学(E辑:技术科学)(08);全文 * |
基于特征线追踪的气动反设计;刘红阳;赵玉新;;推进技术(02);全文 * |
高超声速三维内收缩式进气道/乘波前体一体化设计研究评述;尤延铖;梁德旺;郭荣伟;黄国平;;力学进展(05);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111553119A (zh) | 2020-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107963236B (zh) | 基于密切锥理论的定平面乘波体设计方法 | |
CN104210672B (zh) | 高超声速乘波机身与进气道一体化设计方法 | |
CN112340014A (zh) | 内外流解耦的双乘波高速吸气式飞行器及其生成方法 | |
CN105667812A (zh) | 高超声速飞行器前体、进气道和机翼乘波一体化设计方法 | |
CN112035952B (zh) | 一种模拟飞行器外流的引射喷管实验装置的设计方法 | |
CN105329462A (zh) | 基于可变壁面压力分布规律的吻切流场乘波前体设计方法 | |
CN110589010B (zh) | 高超声速大装载空间乘波体设计方法 | |
CN113153529B (zh) | 基于双入射弯曲激波的宽速域进气道设计方法 | |
CN110450963B (zh) | 高超声速飞行器机体与内转式进气道一体化设计方法及系统 | |
CN111553976B (zh) | 一种激波与压力分布同时可控的三维内转进气道设计方法 | |
CN111003196B (zh) | 一种全乘波飞行器及其设计方法和系统 | |
CN201983920U (zh) | 一种超声速风洞用分段式喷管 | |
CN103950544A (zh) | 通用亚声速扩压器设计方法 | |
CN107016199B (zh) | 一种无激波边界层排移鼓包的设计方法 | |
CN111553119B (zh) | 非均匀来流下高超声速飞行器三维曲面激波反设计方法 | |
CN103500265A (zh) | 一种风洞收缩段曲线的确定方法 | |
CN105539863B (zh) | 高超声速飞行器前体、进气道及支板一体化气动布局方法 | |
CN103823921B (zh) | 一种大涵道比发动机分开式喷管设计方法 | |
CN109815549B (zh) | 一种单对超声速流向旋涡发生装置的设计方法 | |
CN110925091B (zh) | 用于高超内转式进气道的非轴对称广义内锥基准流场设计方法 | |
CN107628266B (zh) | 一种带边界层排移的轴对称预压缩前体的设计方法 | |
CN115199434A (zh) | 一种降低水下排气噪声喷管的设计方法及收缩式波瓣喷管 | |
CN107423481B (zh) | 基于渗透边界反问题方法的全三维内转式进气道设计方法 | |
CN115056998A (zh) | 一种锥导乘波前体纵向分段分级压缩设计方法 | |
CN114261530A (zh) | 最小阻力锥导乘波体和三维内转式进气道一体化设计方法 |
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 |