考虑表面质量引射效应的气动热计算方法与流程

1.本发明涉及气动热和热防护领域，更为具体的，涉及一种考虑表面质量引射效应的气动热计算方法。


背景技术：

2.高超声速飞行器面临着严重的气动热环境，研究人员采用各种热防护手段以突破热障。在各种热防护手段中，烧蚀防热技术已被广泛应用了数十年；而发汗冷却技术因其独特优势，被寄厚望于解决临近空间高超声速飞行器的热防护难题。上述的两种热防护技术均涉及表面质量引射：如烧蚀防热技术中的热解、蒸发/升华与各种化学反应等都包含气体从材料表面析出进入流场的现象和发汗冷却技术中的冷却剂从多孔介质表面射入流场的现象等。
3.现有的考虑表面质量引射效应的气动热计算方法，主要存在如下缺点：1. 总结工程经验的计算方法，虽然计算速度快，但不能探究表面质量引射在高超声速飞行器上对气动热环境产生的细节影响；2. 只考虑了高超声速飞行器表面质量引射现象中垂直表面射流的影响，而未考虑平行表面流动粘性带来的影响；3. 在计算过程中进行模型表面数据的处理时，仅用表面质量引射气体的法向速度来模拟表面质量引射中垂直表面射流，在壁面网格点和临近壁面第一层网格点之间压力梯度为零，不能较真实的反映表面质量引射效应等。
4.综上，现有的考虑表面质量引射效应的气动热计算方法存在考虑现象单一和预测现象不准确、不精细的问题。


技术实现要素：

5.本发明的目的在于克服现有技术的不足，提供考虑表面质量引射效应的气动热计算方法，可以更精准地预测高超声速飞行器表面质量引射效应对气动热环境的影响。
6.本发明的目的是通过以下方案实现的：一种考虑表面质量引射效应的气动热计算方法，包括步骤：s1，获取高超声速飞行器的几何外型；s2，针对获取的高超声速飞行器的几何外型进行网格划分；s3，获取高超声速来流速度数据、高超声速来流温度数据、高超声速来流密度数据和高超声速来流压力数据，并输入表面质量引射气体质量流率数据和表面质量引射气体温度数据到计算机处理器中；s4，在计算机处理器中，基于步骤s3中的数据和利用表面质量引射气体常数数据来计算表面质量引射气体密度数据、表面质量引射气体速度数据、表面质量引射气体压力数据和表面质量引射气体温度数据；s5，在计算机处理器中，将步骤s1~s4中的数据均带入n
‑
s方程进行计算得到高超
声速飞行器壁面热流数据，通过所述壁面热流数据表达高超声速飞行器的气动热环境，然后基于计算机显示单元进行呈现所述气动热环境。
7.进一步地，在步骤s2中，包括如下子步骤：针对高超声速飞行器外型结合工况按由点到线、由线到面、由面到体的顺序逐步绘制高超声速飞行器的几何外型对应的网格，然后完成网格划分。
8.进一步地，在步骤s4中，包括子步骤：s40，输入表面质量引射气体质量流率、表面质量引射气体温度 和表面质量引射气体常数 ，并利用步骤s3中获取的高超声速来流速度数据、高超声速来流温度数据、高超声速来流密度数据和高超声速来流压力数据，通过流场计算获得临近壁面第一层网格点上的压力 、临近壁面第一层网格点上的密度 、临近壁面第一层网格点上的法向速度 ；s41，按照如下公式计算表面质量引射气体密度数据：s42，利用步骤s2计算得到的表面质量引射气体密度数据 和步骤s3输入的表面质量引射气体质量流率 ，按照如下公式计算表面质量引射的垂直射流气体的法向速度 ： ；s43，利用得到的表面质量引射气体密度数据 和步骤s3输入的表面质量引射气体温度 和表面质量引射气体常数 ，按照如下公式计算得到表面质量引射气体压力： 。
9.进一步地，所述高超声速飞行器包括航天器返回舱。
10.进一步地，所述高超声速飞行器包括钝头高超声速再入导弹。
11.进一步地，所述高超声速飞行器包括临近空间高超声速飞行器。
12.进一步地，所述高超声速飞行器包括火星及其他行星进入器。
13.进一步地，在步骤s3中，根据已知的高超声速飞行器弹道即不同时刻高超声速飞行器所处的高度数据和高超声速飞行器马赫数数据，查气体参数表得到高超声速来流声速数据、高超声速来流温度数据、高超声速来流密度数据和高超声速来流压力数据；根据已知的高超声速马赫数数据和查表得到的高超声速来流声速数据求得高超声速飞行器来流速度数据， ， 为高超声速飞行器来流速度、 为高超声速飞行器马赫数、 为高超声速来流声速。
14.进一步地，在步骤s4中，所述表面质量引射气体常数由气体的物性参数表中查得。
15.本发明的有益效果包括：本发明实施例中提供了一种考虑表面质量引射效应的气动热计算方法，该方法较现有技术更加贴近表面质量引射效应。具体的，该方法既考虑了高超声速飞行器上含质量引射区域的壁面垂直射流带来的影响又考虑了平行于壁面流动的气体具有非常大的粘性带来的影响，按以上表面质量引射的特点采用了垂直壁面方向压力梯度不为零的动量守恒方程计算处理，相较于现有技术给出的用法向速度模拟壁面垂直射流的方案更加贴近物理现象，实现更精准地预测高超声速飞行器表面质量引射效应对气动热环境的影响。
16.本发明实施例的计算方法输入端为表面质量引射气体质量流率、表面质量引射气体温度和表面质量引射气体常数数据，表面质量引射气体质量流率在实际的工程实践和实验中更易直接获取，可以更精准地预测高超声速飞行器表面质量引射效应对气动热环境、流动特点等的影响及其影响规律。
17.本发明实施例的计算方法通过更精准、真实的模拟表面质量引射效应，可以为热防护技术中的烧蚀技术和发汗冷却技术的流固热耦合计算研究打下良好的基础。
附图说明
18.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
19.图1为本发明实施例中的计算方法的实现中的方向设置示意图；其中， 是法向方向， ， 是切面上的两个方向；图2为本发明实施例中的考虑表面质量引射效应的气动热计算方法中的模型表面数据处理示意图；图3为本发明实施例中的考虑表面质量引射效应的气动热计算方法中的模型表面数据处理流程图；图4为现有理论计算和实验模型采用的模型图；其中，该模型为半锥角为5
°
的圆锥，其中，从模型头部沿表面的前3.75英寸l为无表面质量引射区域（无阴影部分），即正常固壁面；从3.75英寸处开始再向后16.25英寸，即一直到模型尾端，为表面质量引射区域（阴影部分）；图5为现有理论计算结果和实验数据与采用本发明考虑表面质量引射效应的气动热计算方法进行计算得到的壁面热流结果对比图；其中，横坐标为无量纲化距离（沿模型表面长度s/参考长度l=95.25mm。纵坐标为无量纲化的壁面热流（壁面热流 /参考热流 ）。图5右上角图示中capter代表在capter计算平台中采用本发明的考虑表面质量引射效应的气动热计算方法进行计算的结果，theory代表用理论方法计算得到的结果，experiment代表实验得到的结果。m=0代表注入率为0%，m=0.04代表注入率为0.04%，m=0.08代表注入率为0.08%，对应各个结果的线和符号如图5右上角图示显示。注入率：表面质量引射气体质量流率比来流质量流率；图6为本发明方法实施例的步骤流程图。
具体实施方式
20.本说明书中所有实施例公开的所有特征，或隐含公开的所有方法或过程中的步骤，除了互相排斥的特征和/或步骤以外，均可以以任何方式组合和/或扩展、替换。
21.实施例1如图1~图6所示，一种考虑表面质量引射效应的气动热计算方法，包括步骤：s1，获取高超声速飞行器的几何外型；s2，针对获取的高超声速飞行器的几何外型进行网格划分；s3，获取高超声速来流速度数据、高超声速来流温度数据、高超声速来流密度数据和高超声速来流压力数据，并输入表面质量引射气体质量流率数据和表面质量引射气体温度数据到计算机处理器中；s4，在计算机处理器中，基于步骤s3中的数据和利用表面质量引射气体常数数据来计算表面质量引射气体密度数据、表面质量引射气体速度数据、表面质量引射气体压力数据和表面质量引射气体温度数据；s5，在计算机处理器中，将步骤s1~s4中的数据均带入n
‑
s方程进行计算得到高超声速飞行器壁面热流数据，通过所述壁面热流数据表达高超声速飞行器的气动热环境，然后基于计算机显示单元进行呈现所述气动热环境。
22.实施例2在实施例1的基础上，在步骤s2中，包括如下子步骤：针对高超声速飞行器外型结合工况按由点到线、由线到面、由面到体的顺序逐步绘制高超声速飞行器的几何外型对应的网格，然后完成网格划分。
23.实施例3在实施例1的基础上，在步骤s4中，包括子步骤：s40，输入表面质量引射气体质量流率 、表面质量引射气体温度 和表面质量引射气体常数 ，并利用步骤s3中获取的高超声速来流速度数据、高超声速来流温度数据、高超声速来流密度数据和高超声速来流压力数据，通过流场计算获得临近壁面第一层网格点上的压力 、临近壁面第一层网格点上的密度、临近壁面第一层网格点上的法向速度 ；s41，按照如下公式计算表面质量引射气体密度数据 ：s42，利用步骤s2计算得到的表面质量引射气体密度数据 和步骤s3输入的表面质量引射气体质量流率 ，按照如下公式计算表面质量引射的垂直射流气体的法向速度 ： ；s43，利用得到的表面质量引射气体密度数据 和步骤s3输入的表面质量引射
气体温度 和表面质量引射气体常数 ，按照如下公式计算得到表面质量引射气体压力：。
24.在本发明基于实施例1的其他实施方式中，需要说明的是，所述高超声速飞行器包括航天器返回舱。
25.在本发明基于实施例1的其他实施方式中，需要说明的是，所述高超声速飞行器包括钝头高超声速再入导弹。
26.在本发明基于实施例1的其他实施方式中，需要说明的是，所述高超声速飞行器包括临近空间高超声速飞行器。
27.在本发明基于实施例1的其他实施方式中，需要说明的是，所述高超声速飞行器包括火星及其他行星进入器。
28.在本发明基于实施例1的其他实施方式中，需要说明的是，在步骤s3中，根据已知的高超声速飞行器弹道即不同时刻高超声速飞行器所处的高度数据和高超声速飞行器马赫数数据，查气体参数表得到高超声速来流声速数据、高超声速来流温度数据、高超声速来流密度数据和高超声速来流压力数据；根据已知的高超声速马赫数数据和查表得到的高超声速来流声速数据求得高超声速飞行器来流速度数据， ， 为高超声速飞行器来流速度、 为高超声速飞行器马赫数、 为高超声速来流声速。
29.在本发明基于实施例1的其他实施方式中，需要说明的是，在步骤s4中，所述表面质量引射气体常数由气体的物性参数表中查得。
30.将采用本发明实施例考虑表面质量引射效应的气动热计算方法进行计算得到的结果与marvin, j. g.等人（marvin j, akin, c m, combined effects of mass addition and nose bluntness on boundary
‑
layer transition, aiaa journal, vol. 8, no. 5, 1970, pp. 857
‑
863.ss）先前得到的理论计算解和实验数据进行了对比评估，该理论解是由marvin, j. g.等人研究的通过一种考虑表面质量引射的非相似层流边界层方程计算得到的，该实验数据是marvin, j. g.等人在ames 3.5
‑
ft 高超声速风洞中测得的。计算和实验模型如图4所示，计算和实验的工况如表1所示。图5展示了采用该方法计算得到的壁面热流结果和marvin, j. g.等人计算和实验得到的壁面热流结果的比较。计算时按条件采用对应的边界条件，其中无表面质量引射的壁面区域设置为粘性壁面边界条件，等温壁壁温 ；表面质量引射区域采用本发明实施例所示考虑表面质量引射效应的气动热计算方法中的模型表面数据处理方法，其中表面质量引射气体温度；表面质量引射气体常数按理论计算与实验对应条件设置为空气气体常数；表面质量引射气体质量率由无量纲数注入率得来，注入率等于 ,其中为表面质量引射气体密度， 为表面质量引射气体法向速度， 为来流气体密度， 为来流气体速度，本实施例中对应实验和计算选择了注入率为0%（即无表面质量引射粘性壁面）、0.04%和0.08%三个条件。由图5中可以看到结果大致相同，且趋势一致。
31.表1 计算和实验的工况数据
ꢀ
（来流空气马赫数） （来流空气温度） （来流空气密度）
ꢀꢀ
（来流空气压力）（来流空气压力）（来流空气压力）（来流空气压力）
本发明功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，在一台计算机设备（可以是个人计算机，服务器，或者网络设备等）以及相应的软件中执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、或者光盘等各种可以存储程序代码的介质，进行测试或者实际的数据在程序实现中存在于只读存储器（random access memory，ram）、随机存取存储器（random access memory，ram）等。