CN111090907B - 一种飞行试验转捩判断方法 - Google Patents

一种飞行试验转捩判断方法 Download PDF

Info

Publication number
CN111090907B
CN111090907B CN201911043334.6A CN201911043334A CN111090907B CN 111090907 B CN111090907 B CN 111090907B CN 201911043334 A CN201911043334 A CN 201911043334A CN 111090907 B CN111090907 B CN 111090907B
Authority
CN
China
Prior art keywords
heat flow
transition
measuring point
phi
state
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
CN201911043334.6A
Other languages
English (en)
Other versions
CN111090907A (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.)
China Academy of Aerospace Aerodynamics CAAA
Original Assignee
China Academy of Aerospace Aerodynamics CAAA
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 China Academy of Aerospace Aerodynamics CAAA filed Critical China Academy of Aerospace Aerodynamics CAAA
Priority to CN201911043334.6A priority Critical patent/CN111090907B/zh
Publication of CN111090907A publication Critical patent/CN111090907A/zh
Application granted granted Critical
Publication of CN111090907B publication Critical patent/CN111090907B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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

  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Abstract

本发明给出了一种飞行试验转捩判断方法。该方法通过采用数值计算的层流状态热流数据，对飞行试验的热流进行无量纲化，根据该无量纲热流的大小，对飞行器各测点的转捩状态进行判断。该方法不仅可以剔除来流参数、攻角变化等对热流的影响，而且可以对背风面转捩进行判断，更加直观地显示转捩阵面推进过程。

Description

一种飞行试验转捩判断方法
技术领域
本发明涉及一种飞行试验转捩判断方法，属于空气动力学技术领域。
背景技术
边界层转捩对高超声速飞行器的推进和热防护系统设计有重要影响。但由于转捩过程受到多种因素的影响，机理复杂，无论是理论分析、数值计算或风洞试验，准确预测转捩困难重重。国内外诸多研究机构试图通过开展转捩飞行试验提供真实飞行数据，从而完善转捩理论、验证和发展转捩计算模型、确认风洞试验结果，进一步增强转捩预测能力。
在分析转捩飞行试验结果时，一般把热流曲线发生突跃的位置作为转捩起始位置。但由于“曲线突跃”没有定量的判断准则，而且随着飞行状态(高度、马赫数、攻角、侧滑角)变化，也会导致热流突然增加；此外飞行器背风面热流绝对数值小且变化平缓。以上情况均无法简单的通过曲线突跃去准确判断是否转捩。
此时可以借助高精度的层流数值计算去判断转捩。相较于转捩和湍流状态，层流数值模拟不受各种物理模型的影响，精度高、可靠性强。采用数值计算获得飞行器表面的层流热流，对飞行试验热流进行无量纲化，获得无量纲热流，以衡量飞行试验热流偏离层流状态的程度。当该无量纲参数超过一定数值时，即认为转捩发生。无量纲热流不仅可以剔除来流参数、飞行姿态等对热流的影响，而且可以对背风面转捩进行判断，直观地研究转捩阵面推进过程。
发明内容
本发明解决的技术问题是：克服现有方法在判断转捩时的主观性以及无法对背风面转捩进行判断，提供了一种飞行试验转捩判断方法。
本发明解决技术的方案是：一种飞行试验转捩判断方法，步骤如下：
(1)通过飞行试验，实测飞行器壁面不同测点处温度数据，根据实测的温度数据，确定净热流以及净热流随时间的演化曲线；
(2)根据上述飞行试验实际飞行情况重构飞行器弹道，在该弹道上选取一系列状态点，采用数值方法，得到层流状态下各状态点壁面热流分布，进而确定飞行器试验中各测点位置的计算热流；
(3)根据步骤(1)(2)的结果计算每个测点的无量纲热流φ，φ≥2时，判断该测点处于转捩或湍流状态，否则处于层流状态。
优选的，所述的无量纲热流φ通过下述方式确定：
首先，将步骤(1)中确定的净热流转换至壁面温度为T时的热流所述的T与步骤(2)中计算热流时的壁面温度相同；
然后，利用公式计算无量纲热流；
其中，为步骤(2)中计算的热流。
优选的，所述的无量纲热流φ通过下述方式确定：
首先，将步骤(2)中计算的热流转换为各测点的净热流转换过程中使用的壁面温度采用飞行试验中实测的温度数据；
然后，利用公式计算无量纲热流；
其中，为步骤(1)中计算的净热流。
优选的，所述步骤(2)中的状态点要求能够反映重构弹道的重要特征，包括弹道顶点、攻角突变前后、侧滑角突变前后，以及其它会引起热流突变的状态点。
优选的，在步骤(3)中根据得到无量纲热流φ结合对应测点的位置，生成φ的云图，直观、清晰的看到转捩阵面随着状态点的变化过程。
优选的，根据每个时刻各个测点的无量纲热流φ，用横坐标表示时间，纵坐标表示各测点，生成φ的二维表格，将φ≥2的单元格突出显示，清晰直观的显示转捩的时空变化。
优选的，步骤(2)中的数值方法采用有限体积法求解三维可压缩雷诺平均Navier-Stokes方程，空间离散选择AUSM+格式并使用Minmod限制器，时间推进选择无条件稳定的LU-SGS方法。
优选的，在步骤(3)之后执行如下步骤：
首先，根据步骤(3)的判断结果，确定转捩位置；
然后，根据数值计算中确定的每个时刻、各个测点位置的基于动量损失厚度的雷诺数Reθ和边界层外缘马赫数Mae，得到转捩位置与Reθ和Mae的关系，将该关系作为转捩准则。
优选的，通过热流辨识方法，得到测点的净热流随时间的演化曲线。
优选的，所述步骤(2)中由壁面热流分布确定测点位置的计算热流，通过Tecplot软件的数据插值功能而获得，进一步采用Python进行批处理，快速得到测点位置在系列状态点下的计算热流。
本发明与现有技术相比的有益效果是：
目前飞行器转捩的判断主要是根据热流曲线的突跃进行确定，更多依赖于主观观察，误差较大，而且无法对热流变化平缓的背风面转捩情况进行判断。采用无量纲热流，不仅可以剔除来流参数、飞行姿态等对热流的影响，而且可以对背风面转捩进行判断，更加直观地显示转捩阵面推进过程。
附图说明：
图1为本发明方法流程图。
具体实施方式
下面结合实施例对本发明作进一步阐述。
参照图1，本发明提供了一种飞行器背风面转捩的判断方法，步骤如下：
(1)在飞行器特定位置布置温度传感器，通过飞行试验，可以获得飞行器壁面不同位置处温度传感器的实测数据，即温度随时间的演化；
飞行试验结束后，首先要对测量采集系统的数据进行有效性分析，尤其是同轴热电偶的测温数据。一方面可以与理论计算数值进行比较，另一方面，如果相近位置安装了热流传感器，也可以先通过热流辨识获得热流，再与热流传感器的结果作比较。
热流辨识方法有一维模型辨识、二维轴对称模型和三维模型。下面介绍一维热流辨识，其它方法可以参照相关文献资料。一维辨识的计算公式如下：
其中，k,ρ和Cp为热电偶传感器的物性参数，tj和Tj分别是时间序列和对应的温度序列。
当测温数据的有效性得到了验证，就可以利用温度数据进行转捩分析。对每个温度传感器得到的实测温度进行热流辨识，得到这些测点的净热流随时间的演化。考察热流曲线，可以根据热流曲线的突跃，初步判断转捩时间和位置。为了方便与数值计算获得的结果进行比较，按照以下公式，换算得到300K等壁温时的热流
其中，TF是飞行试验的实测温度，QF是辨识热流，T0是相应时刻的总温。总温可以根据重构弹道获得。
采用300K等壁温热流，是为了方便与数值计算的结果直接进行比较。
(2)从飞行试验的重构弹道上选取一系列状态点，这些状态点应该能够反映重构弹道的重要特征，如攻角、侧滑角的改变等；状态点的选取要适当，一方面要包括弹道最高点、飞行器攻角、侧滑角变化前后的状态点，另一方面不能过多选取以导致庞大的计算量。
进行高精度数值计算时，可以采用有限体积法求解三维可压缩雷诺平均Navier-Stokes方程，选择合适的数值格式，即空间离散选择AUSM+格式并使用Minmod限制器，时间推进选择无条件稳定的LU-SGS方法，壁面设置为300K等壁温，计算各状态点在层流时的表面热流分布，用表示。通过插值，得到飞行器测点位置的计算热流。
(3)根据飞行试验和数值计算，由试验辨识和数值计算的300K等壁温热流，定义无量纲热流φ，即：
该无量纲热流可以衡量飞行试验热流偏离层流状态的程度。如果测点的φ接近1，说明测点处于层流状态；反之，当测点发生转捩或处于湍流状态时，无量纲热流φ会明显大于1。
综合层流计算精度和网格的影响，φ＞2时，即实际飞行热流大于层流的2倍时，则可以判断为处于转捩或湍流状态，否则处于层流状态。
(4)由各个测点的无量纲热流随时间的变化，可以判断各个测点发生转捩的时刻。在每个弹道点，使用专业绘图软件Origin，根据该时刻各测点的无量纲热流φ，得到φ的云图，就可以直观、清晰的看到转捩阵面随着弹道点的变化过程。
(5)使用微软办公软件Excel，根据每个时刻各个测点的无量纲热流φ，即用横坐标表示时间，纵坐标表示各测点，生成φ的二维表格。然后利用Excel的数据分析功能，将φ＞2的单元格突出显示，也可以清晰直观的显示转捩的时空变化。
本发明未详细说明部分属于本领域技术人员的公知常识。

Claims (8)

1.一种飞行试验转捩判断方法，其特征在于，步骤如下：
(1)通过飞行试验，实测飞行器壁面不同测点处温度数据，根据实测的温度数据，确定净热流以及净热流随时间的演化曲线；
(2)根据上述飞行试验实际飞行情况重构飞行器弹道，在该弹道上选取一系列状态点，采用数值方法，得到层流状态下各状态点壁面热流分布，进而确定飞行器试验中各测点位置的计算热流；
(3)根据步骤(1)(2)的结果计算每个测点的无量纲热流φ，φ≥2时，判断该测点处于转捩或湍流状态，否则处于层流状态；
所述的无量纲热流φ通过以下任一方式确定：
1)、首先，将步骤(1)中确定的净热流转换至壁面温度为T时的热流所述的T与步骤(2)中计算热流时的壁面温度相同；
然后，利用公式计算无量纲热流；
其中，为步骤(2)中计算的热流；
2)、首先，将步骤(2)中计算的热流转换为各测点的净热流转换过程中使用的壁面温度采用飞行试验中实测的温度数据；
然后，利用公式计算无量纲热流；
其中，为步骤(1)中的净热流。
2.根据权利要求1所述的方法，其特征在于：所述步骤(2)中的状态点要求能够反映重构弹道的重要特征，包括弹道顶点、攻角突变前后、侧滑角突变前后，以及其它会引起热流突变的状态点。
3.根据权利要求1所述的方法，其特征在于：在步骤(3)中根据得到无量纲热流φ结合对应测点的位置，生成φ的云图，显示转捩阵面随着状态点的变化过程。
4.根据权利要求1所述的方法，其特征在于：根据每个时刻各个测点的无量纲热流φ，用横坐标表示时间，纵坐标表示各测点，生成φ的二维表格，将φ≥2的单元格突出显示，显示转捩的时空变化。
5.根据权利要求1所述的方法，其特征在于：步骤(2)中的数值方法采用有限体积法求解三维可压缩雷诺平均Navier-Stokes方程，空间离散选择AUSM+格式并使用Minmod限制器，时间推进选择无条件稳定的LU-SGS方法。
6.根据权利要求1所述的方法，其特征在于：在步骤(3)之后执行如下步骤：
首先，根据步骤(3)的判断结果，确定转捩位置；
然后，根据数值计算中确定的每个时刻、各个测点位置的基于动量损失厚度的雷诺数Reθ和边界层外缘马赫数Mae，得到转捩位置与Reθ和Mae的关系，将该关系作为转捩准则。
7.根据权利要求1所述的方法，其特征在于：通过热流辨识方法，得到测点的净热流随时间的演化曲线。
8.根据权利要求1所述的方法，其特征在于：所述步骤(2)中由壁面热流分布确定测点位置的计算热流，通过Tecplot软件的数据插值功能而获得，进一步采用Python进行批处理，快速得到测点位置在系列状态点下的计算热流。
CN201911043334.6A 2019-10-30 2019-10-30 一种飞行试验转捩判断方法 Active CN111090907B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911043334.6A CN111090907B (zh) 2019-10-30 2019-10-30 一种飞行试验转捩判断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911043334.6A CN111090907B (zh) 2019-10-30 2019-10-30 一种飞行试验转捩判断方法

Publications (2)

Publication Number Publication Date
CN111090907A CN111090907A (zh) 2020-05-01
CN111090907B true CN111090907B (zh) 2023-09-29

Family

ID=70393184

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911043334.6A Active CN111090907B (zh) 2019-10-30 2019-10-30 一种飞行试验转捩判断方法

Country Status (1)

Country Link
CN (1) CN111090907B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111832159B (zh) * 2020-06-23 2023-08-29 北京临近空间飞行器系统工程研究所 一种基于飞行试验数据的边界层转捩阵面动态演化过程确定方法
CN113532722B (zh) * 2021-05-25 2023-04-14 北京临近空间飞行器系统工程研究所 一种基于飞行试验脉动压力数据的双谱分析转捩辨识方法
CN115809513B (zh) * 2023-02-08 2023-05-26 中国空气动力研究与发展中心计算空气动力研究所 一种强迫转捩-俯仰震荡数值模拟方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106872195A (zh) * 2017-01-06 2017-06-20 北京临近空间飞行器系统工程研究所 一种高速飞行器气动热飞行试验数据的关联分析方法
CN108195542A (zh) * 2017-12-25 2018-06-22 中国航天空气动力技术研究院 一种飞行试验测点位置的流态判读方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7251592B1 (en) * 2003-08-18 2007-07-31 United Technologies Corporation Boundary layer transition model
US8600710B2 (en) * 2011-10-05 2013-12-03 King Fahd University Of Petroleum And Minerals Method of modeling thermal problems using a non-dimensional finite element method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106872195A (zh) * 2017-01-06 2017-06-20 北京临近空间飞行器系统工程研究所 一种高速飞行器气动热飞行试验数据的关联分析方法
CN108195542A (zh) * 2017-12-25 2018-06-22 中国航天空气动力技术研究院 一种飞行试验测点位置的流态判读方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
MF-1模型飞行试验表面压力与温度测量技术研究;杨庆涛 等;空气动力学学报;第35卷(第05期);732-741 *
利用表面脉动压力探测转捩的数据处理方法研究;王猛 等;实验力学;第27卷(第06期);764-770 *
飞行试验热流辨识和边界层转捩滞后现象;国义军 等;航空学报;第38卷(第10期);90-100 *
高超声速飞行器气动热关联换算方法研究;赵金山 等;力学学报;第50卷(第05期);1235-1245 *

Also Published As

Publication number Publication date
CN111090907A (zh) 2020-05-01

Similar Documents

Publication Publication Date Title
CN111090907B (zh) 一种飞行试验转捩判断方法
Medic et al. Prediction of transition and losses in compressor cascades using large-eddy simulation
CN103901291B (zh) 一种变电设备内部绝缘缺陷的诊断方法
CN109871625A (zh) 基于Johnson变换的非高斯风压模拟方法
CN105045973A (zh) 一种变弧长自适应采样方法
CN113947035A (zh) 一种高超声速边界层转捩数据天地相关性方法
Vidanović et al. Validation of the CFD code used for determination of aerodynamic characteristics of nonstandard AGARD-B calibration model
CN113722860B (zh) 基于降阶模型的瞬态热力状态在线评估方法、设备及介质
CN110377941A (zh) 建立卫星温度场的罚盲似然克里金代理模型的方法
Bartel et al. Force calibration using errors-in-variables regression and Monte Carlo uncertainty evaluation
CN113128047A (zh) 一种航空发动机推力计算方法、存储介质及计算设备
CN115510732B (zh) 一种基于深度学习的方舱红外特性仿真快速算法
CN115455838B (zh) 一种面向时程数据的高空间分辨率流场重构方法
CN109598059B (zh) 一种基于代理模型的热防护系统优化设计方法及设计系统
CN116432548A (zh) 基于迁移学习的融合风速计算方法、装置、设备及介质
CN114462336B (zh) 一种核反应堆主管道冷却剂平均温度计算方法
JP6867072B2 (ja) 気象予測データ作成プログラム、気象予測データ作成方法、及び、移動体
CN114722695A (zh) 一种基于无量纲输入输出神经网络的fads解算系统及方法
CN114295214A (zh) 一种基于有效发射率的复杂环境下涡轮叶片辐射测温方法和装置
Shenghui et al. Sensitivity analysis of the VVP wind retrieval method for single-Doppler weather radars
CN112512267A (zh) 基于移动机器人的温度监控方法、系统、设备和存储介质
Wright Validation Process for LEWICE by Use of a Naviér-Stokes Solver
Vedam et al. Evaluation of Gradient and Curvature-Based Adaptive Mesh Refinement for Viscous Transonic Flows
Heinemann et al. Experiences with the Application of the non-hydrostatic Mesoscale Model GESIMA for assessing Wind Potential in Complex Terrain
Teng et al. Design of an optimal plan of accelerated degradation test via acceleration factor constant principle

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