本发明属于气动热与热防护领域，具体涉及一种采用尺度效应修正湍流无量纲热流的方法。
背景技术：
1、飞行器气动热环境预测的精准度直接决定着高速飞行器热防护系统的可靠性与先进性，若气动热环境预测偏低，则影响高速飞行器的安全，若气动热环境预测偏高，则热防护系统的冗余质量将挤占高速飞行器宝贵的有效载荷。
2、当前，气动热风洞试验依然是获取高速飞行器表面热环境的主要手段，但由于风洞尺寸的限制，在开展气动热风洞试验时需要对模型进行缩比，小模型上获取的表面热环境与全尺寸模型的表面热环境存在差异。在湍流情况下，不仅测得的绝对热流不能直接应用于飞行器表面气动热环境设计，而且以各自驻点归一化的无量纲热流也随着尺度的变化有所不同，这就给湍流气动热数据应用带来很大的困难。
3、当前，亟需发展一种采用尺度效应修正湍流无量纲热流的方法。
技术实现思路
1、本发明所要解决的技术问题是提供一种采用尺度效应修正湍流无量纲热流的方法，用以克服现有技术的缺陷。
2、尺度效应对湍流气动热的影响，可以通过边界层方程来推导证明。在湍流情况下，计算热流的公式包括walker公式、zoby公式和平板参考焓公式。在 re ∞＜107 1/m时，walker公式、zoby公式和平板参考焓公式推导得到的尺度效应影响因子是一致的，本发明的采用尺度效应修正湍流无量纲热流的方法选择平板参考焓公式计算热流。
3、本发明的采用尺度效应修正湍流无量纲热流的方法，包括以下步骤：
4、s10.根据飞行器尺度和风洞模拟能力，确定模型缩比；其中，下标1表示小模型，下标2表示大模型； l 1为小模型总长， l 2为大模型总长；
5、s20.通过气动热风洞试验，获得湍流流态对应的小尺度模型表面的热流信息，热流信息包括驻点热流、物面热流；
6、s30.采用驻点热流对物面热流进行归一化，计算小模型表面每个特征点的无量纲热流；
7、s40.引入尺度效应关系式，计算出全流场内每个特征点处的修正因子；
8、a取0.3；如果忽略边界层厚度变化的影响，需要考虑气流压缩性的影响，当单位雷诺数 re ∞≥107 1/m时，a的取值变更为0.357；
9、a的取值对结果的影响为±10%；
10、s50.利用，求解大模型表面的无量纲热流；
11、s60.将大模型表面的无量纲热流再乘以大模型驻点热流 q s2，获得大模型物面热流 q w2。
12、进一步地，s30中所述的归一化，包括以下步骤：
13、s31.采用fay-riddell公式计算驻点热流：
14、；
15、其中，为驻点热流，为普朗特数，为物面密度，为物面粘性系数，为驻点密度，为驻点粘性系数，为驻点焓，为物面焓，为边界层外缘速度，为从头部驻点至物面当地计算位置的流线长；
16、s32.对于单位雷诺数 re ∞＜107 1/m的情况，采用层流流态平板参考焓公式计算物面热流；
17、；
18、其中，为物面热流，为以x作为特征尺度的雷诺数，为参考密度，为边界层外缘密度，为参考粘性系数，为边界层外缘粘性系数；
19、忽略边界层厚度变化的影响，利用驻点热流进行无量纲化，则有：
20、；
21、其中，为边界层外缘马赫数，为驻点压力，为驻点温度，表示模型的等效头部半径；
22、假设来流参数一致，并以大模型总长 l 2和小模型总长 l 1进行归一化，则大模型和小模型上的无量纲热流分布为：
23、；
24、其中，a=0.3；
25、s33.对于单位雷诺数 re ∞≥107 1/m的情况，采用湍流流态平板参考焓公式计算湍流流态下的热流；
26、；
27、其中，为湍流流态下的热流；
28、同样的，假设来流参数一致，并以大模型总长 l 2和小模型总长 l 1进行归一化，则大模型和小模型上的无量纲热流分布为：
29、；
30、其中，a=0.357。
31、本发明的采用尺度效应修正湍流无量纲热流的方法，通过引入不同缩比模型的总长比值，对气动热风洞试验中获得的湍流无量纲热流进行修正，将气动热风洞试验中小模型上获得的湍流无量纲热流外推到全尺寸模型，进行全尺寸模型热流分布预测。本发明的采用尺度效应修正湍流无量纲热流的方法揭示了湍流无量纲热流随模型尺度的变化规律，解决了气动热风洞试验数据应用于多个轨道点的设计问题，提高了气动热风洞试验数据应用精准度。
1.一种采用尺度效应修正湍流无量纲热流的方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的采用尺度效应修正湍流无量纲热流的方法，其特征在于，s30中所述的归一化，包括以下步骤：
