一种高速飞行器翼舵前缘转捩预测装置及方法与流程

本发明涉及一种高速飞行器翼舵前缘转捩预测装置及方法，属于转捩预测技术领域。

背景技术：

高速飞行器(马赫数大于5)在飞行试验过程中，飞行器的翼舵前缘会出现层流和湍流两种流态，流态的不同会导致前缘出现高热流区域，湍流流态比层流流态的前缘热流高3-5倍，属于飞行器热防护设计的重点区域。准确的预测飞行器后掠尖前缘的转捩是飞行器设计的关键，决定飞行试验成败。

目前，高速飞行器的翼舵前缘均采用带有后掠角的尖前缘(前缘半径在2-10mm之间)，目前常用的翼舵前缘转捩预测方法为计算机数值模拟。计算机数值模拟，通过对流动转捩过程进行分析、抽象，建立流动数学模型，采用计算机模拟计算转捩。其优点是能够模拟流动的转捩过程，获得飞行器表面精确的转捩位置。缺点是所建立的数学模型需要经过试验验证，计算量大，无法应用于工程实际。

相对于飞行器头部前缘的转捩预测，飞行器翼舵前缘由于受到飞行器机体干扰，流动更加复杂。对于飞行器翼舵前缘的转捩预测，采用传统的计算机数值模拟方法，计算量大，成本高，无法高效快速的适用于工程计算，而且预测结果需要经过试验验证。目前工程上无可用的快速预测飞行器翼舵前缘的方法。

技术实现要素：

本发明的目的在于克服现有技术不足，提供一种降低计算量、精度高、能应用于工程的高速飞行器翼舵前缘转捩预测装置及方法。

本发明的技术解决方案：一种高速飞行器翼舵前缘转捩预测装置，包括斜激波波后流动参数确定模块、转捩判定准则系数确定模块和转捩预测模块，所述的斜激波波后流动参数确定模块包括激波干扰判断模块、有激波干扰斜激波波后流动参数确定模块和无激波干扰斜激波波后流动参数确定模块；

所述的激波干扰判断模块依据高速飞行器的飞行工况(包括马赫数、雷诺数、攻角等)和飞行器前缘结构(前缘后掠角和前缘半径)，采用数值模拟方法计算飞行器的无粘流场，判别头部激波与翼舵前缘斜激波是否相交，若相交则判定为有激波干扰，若不相交则判定为无激波干扰；

有激波干扰斜激波波后流动参数确定模块从飞行器的无粘流场获得激波相交位置处沿前缘后掠方向到翼舵前缘顶部的头部激波后的流动参数(密度、速度和粘性系数)，再根据头部激波后的流动参数和翼舵前缘几何参数，采用工程计算方法(无粘流动的斜激波方程计算，详见《气体动力学》，童秉纲)得到翼舵前缘斜波波后参数，翼舵前缘斜波波后参数包括斜激波波后的密度ρe、斜激波波后的速度ue和斜激波波后的粘性系数μe；

无激波干扰斜激波波后流动参数确定模块从飞行器的无粘流场获得头部激波波后的流动参数(密度、速度、粘性系数)，再根据头部激波波后的流动参数和翼舵前缘几何参数，采用工程计算方法得到翼舵前缘后掠斜激波波后参数，翼舵前缘斜波波后参数包括斜激波波后的密度ρe、斜激波波后的速度ue和斜激波波后的粘性系数μe；

所述的转捩判定准则系数确定模块利用公式得到转捩判定准则系数red；

所述的转捩预测模块将转捩判定准则系数red与预设的翼舵前缘转捩判定准则阈值redδ进行比较，若red≥redδ，则翼舵前缘发生转捩，为湍流流态，若red＜redδ，则翼舵前缘不发生转捩，为层流流态。

本发明还包括翼舵前缘转捩判定准则阈值确定模块，所述的翼舵前缘转捩判定准则阈值确定模块包括转捩后掠角确定模块、转捩斜激波波后参数提取模块、风洞试验转捩判定准则系数确定模块和转捩判定准则阈值对比判定模块，

所述的转捩后掠角确定模块根据若干头部前缘一致、翼舵后掠角不同、翼舵前缘半径一致的风洞模型在不同试验工况下风洞试验得到的每个试验工况下热流与后掠角的关系，再根据在每个试验工况下采用数值模拟方法得到的全层流和全湍流曲线，得到每个试验工况下发生转捩的后掠角；

所述的转捩斜激波波后参数提取模块提取与每个试验工况下转捩后掠角最接近的风洞模型进行的风洞试验，得到每个试验工况下翼舵前缘斜激波波后流动参数；

所述的风洞试验转捩判定准则系数确定模块利用公式得到每个试验工况下的转捩判定准则系数；

所述的转捩判定准则阈值对比判定模块将每个试验工况下的转捩判定准则系数进行对比，选取所有工况下的判定准则系数的最小值，作为翼舵前缘转捩判定准则阈值redδ。

一种高速飞行器翼舵前缘转捩预测方法，通过以下步骤实现：

第一步，确定高速飞行器翼舵前缘的斜激波波后流动参数，斜激波波后流动参数包括斜激波波后的密度ρe、斜激波波后的速度ue和斜激波波后的粘性系数μe；

a1.1、判断头部激波与翼舵前缘斜激波是否相交，若相交则进行步骤a1.2，若不相交则进行步骤a1.3；

依据高速飞行器的飞行工况(包括马赫数、雷诺数、攻角等)和飞行器前缘结构(前缘后掠角和前缘半径)，采用数值模拟方法计算飞行器的无粘流场，判别头部激波与翼舵前缘斜激波是否相交，若相交则判定为有激波干扰，若不相交则判定为无激波干扰。

数值模拟方法为本领域公知技术，可以采用现有的商用软件实现。

a1.2、从飞行器的无粘流场获得激波相交位置处沿前缘后掠方向到翼舵前缘顶部的头部激波后的流动参数，再根据头部激波后的流动参数和翼舵前缘几何参数，采用工程计算方法得到翼舵前缘斜波波后参数，翼舵前缘斜波波后参数包括斜激波波后的密度ρe、斜激波波后的速度ue和斜激波波后的粘性系数μe，转入第二步；

a1.3、从飞行器的无粘流场获得头部激波波后的流动参数，再根据头部激波波后的流动参数和翼舵前缘几何参数，采用工程计算方法得到翼舵前缘后掠斜激波波后参数，翼舵前缘斜波波后参数包括斜激波波后的密度ρe、斜激波波后的速度ue和斜激波波后的粘性系数μe，转入第二步；

第二步，根据第一步得到的斜激波波后流动参数，利用公式(1)得到高速飞行器翼舵前缘的转捩判定准则系数，

其中red为转捩判定准则系数，d为前缘直径；

第三步，转捩预测，

将第二步得到的转捩判定准则系数red与翼舵前缘转捩判定准则阈值redδ进行比较，若满足red≥redδ，则翼舵前缘发生转捩，为湍流流态，若不满足red≥redδ，则翼舵前缘不发生转捩，为层流流态。

所述的翼舵前缘转捩判定准则阈值redδ通过以下步骤确定，

a3.1、制作与高速飞行器头部结构一致、翼舵前缘半径一致、不同翼舵前缘后掠角的1：1的风洞模型；

后掠角的范围为0～90°，从中选取若干后掠角，包括极值，后掠角取值越多，风洞模型就越多，精度越高，风洞模型的设计与制作为本领域公知技术。

a3.2、根据高速飞行器飞行轨迹参数，确定风洞试验参数的范围，风洞试验参数包括马赫数、攻角和雷诺数，确定一系列涵盖风洞试验参数范围的试验工况，在每一个试验工况条件下对步骤a3.1制作的一系列风洞模型开展风洞试验，得到每个试验工况下热流与后掠角的关系，再根据在每个试验工况下采用数值模拟方法得到的全层流和全湍流曲线，得到每个试验工况下发生转捩的后掠角，提取与每个试验工况下转捩后掠角最接近的风洞模型进行的风洞试验，得到每个试验工况下头部前缘斜激波波后参数；

风洞试验和数值模拟方法为本领域公知技术。

a3.3、利用公式(1)得到不同工况试验状态下的转捩判定准则系数；

a3.4、对步骤a3.3得到的转捩判定准则系数进行对比，选取所有工况下的判定准则系数的最小值，获得翼舵前缘转捩判定准则阈值redδ。

本发明与现有技术相比的有益效果：

(1)本发明通过确定更为精确转捩判定准则，使转捩判定简单、精确，适用于工程应用；

(2)本发明转捩判定准则结合了飞行数据、数值模拟和风洞试验，获得的转捩判定流态准确，可以应用于后续飞行器设计，可信度高。

附图说明

图1为本发明原理框图；

图2为本发明流程图；

图3为本发明实施例中某工况下热流与后掠角的关系图。

具体实施方式

下面结合具体实例及附图对本发明进行详细说明。

实施例1

本发明如图1所示，提供了一种转捩预测装置，包括斜激波波后流动参数确定模块、转捩判定准则系数确定模块、翼舵前缘转捩判定准则阈值确定模块和转捩预测模块。斜激波波后流动参数确定模块包括激波干扰判断模块、有激波干扰斜激波波后流动参数确定模块和无激波干扰斜激波波后流动参数确定模块。翼舵前缘转捩判定准则阈值确定模块包括转捩后掠角确定模块、转捩斜激波波后参数提取模块、风洞试验转捩判定准则系数确定模块和转捩判定准则阈值对比判定模块。

激波干扰判断模块依据高速飞行器的飞行工况，采用数值模拟方法计算飞行器的无粘流场，判别头部激波与翼舵前缘斜激波是否相交，若相交则判定为有激波干扰，若不相交则判定为无激波干扰；本实例计算发现头部激波与翼舵前缘激波不相交。

无激波干扰斜激波波后流动参数确定模块从飞行器的无粘流场获得头部激波波后的流动参数，再根据头部激波波后的流动参数和翼舵前缘几何参数，采用工程计算方法得到翼舵前缘后掠斜激波波后参数，翼舵前缘斜波波后参数包括斜激波波后的密度ρe＝＝0.653kg/m³、斜激波波后的速度ue＝769m/s和斜激波波后的粘性系数μe＝0.00000961pa·s。

转捩判定准则系数确定模块利用公式得到转捩判定准则系数red＝209000，本实例中前缘半径为5mm。

转捩后掠角确定模块根据若干头部前缘一致、翼舵后掠角不同、翼舵前缘半径一致的风洞模型在不同试验工况下风洞试验得到的每个试验工况下热流与后掠角的关系，再根据在每个试验工况下采用数值模拟方法得到的全层流和全湍流曲线，得到每个试验工况下发生转捩的后掠角。转捩斜激波波后参数提取模块提取与每个试验工况下转捩后掠角最接近的风洞模型进行的风洞试验，得到每个试验工况下翼舵前缘斜激波波后流动参数。风洞试验转捩判定准则系数确定模块利用公式得到每个试验工况下的转捩判定准则系数。转捩判定准则阈值对比判定模块将每个试验工况下的转捩判定准则系数进行对比，选取所有工况下的判定准则系数的最小值，作为翼舵前缘转捩判定准则阈值redδ。本实例中翼舵前缘转捩判定准则阈值redδ＝170000。

转捩预测模块将转捩判定准则系数red与预设的翼舵前缘转捩判定准则阈值redδ＝170000进行比较，red≥redδ，则翼舵前缘发生转捩，为湍流流态。

实施例2

本发明如图2所示，还提供了一种转捩预测方法，通过以下步骤实现：

1、确定高速飞行器翼舵前缘的斜激波波后流动参数，斜激波波后流动参数包括斜激波波后的密度ρe、斜激波波后的速度ue和斜激波波后的粘性系数μe。

(1)判断头部激波与翼舵前缘斜激波是否相交

依据高速飞行器的飞行工况，采用数值模拟方法计算飞行器的无粘流场，头部激波与翼舵前缘斜激波相交，判定为有激波干扰。

(2)从飞行器的无粘流场获得激波相交位置处处沿前缘后掠方向到翼舵前缘顶部的头部激波后的流动参数，再根据头部激波后的流动参数和翼舵前缘几何参数，采用工程计算方法得到翼舵前缘斜波波后参数，翼舵前缘斜波波后参数包括斜激波波后的密度ρe＝＝1.47kg/m³、斜激波波后的速度ue＝280m/s和斜激波波后的粘性系数μe＝0.0000257pa·s。

2、利用公式(1)得到高速飞行器翼舵后掠尖前缘的转捩判定准则系数red＝72500，本实例中前缘半径为5mm。

3、确定翼舵前缘转捩判定准则阈值redδ

(1)制作与高速飞行器头部结构一致、翼舵前缘半径一致、不同翼舵前缘后掠角的1：1的风洞模型。

后掠角范围0～90°，根据每5°或10°选取一个角度，制作0°、10°、20°、25°、30°、35°、40°、45°、50°、60°、70°、80°、90°共13个风洞模型。

(2)根据高速飞行器飞行轨迹参数，确定风洞试验参数范围，本实例中马赫数为4～8，攻角-20°～20°，雷诺数为(0.45～3.6)×10⁷，从中确定一系列试验工况，每个工况下对所有的风洞模型进行风洞试验。例如，采用工况马赫数5.0、雷诺数2.3×10⁷，攻角0°进行风洞试验，得到如图3所示的热流与后掠角的关系图，再根据在此工况下采用数值模拟方法得到的全层流和全湍流曲线，从中得到在此工况条件下转捩后掠角为20°。提取后掠角为20°的风洞模型在该工况下的风洞试验，得到飞行器翼舵前缘斜激波波后参数。同理，得到其他工况条件下，飞行器翼舵前缘斜激波波后参数。

(3)利用公式(1)得到不同工况试验状态下的转捩判定准则系数。

上述工况马赫数5.0、雷诺数2.3×10⁷，攻角0°，后掠角20°风洞模型的转捩判定准则系数为1.26e+05，其他工况条件下转捩判定准则系数分别为1.32e+05、1.12e+05、1.55e+05、1.73e+05、2e+04、1.45e+05、1.81e+05…。

(4)对步骤(3)得到的转捩判定准则系数进行对比，选取最小值作为翼舵前缘转捩判定准则阈值redδ＝20000。

4、转捩预测

对比步骤2得到的red和步骤3得到的redδ，red＞redδ，则翼舵前缘发生转捩，为湍流流态。

本发明未详细说明部分为本领域技术人员公知技术。