本发明是一种超声速飞行器气动热环境预测领域的技术，具体是一种新型的气动热预测的智能计算方法。
背景技术：
目前高超声速飞行器是各国军事航天领域研究的重点，如第四代第五代战机，其可以实现快速精准的打击，在未来战争中具有极强的威慑力。对于高超声速飞行器的设计来讲高超声速飞行条件下飞行器的热防护问题一直是高超声速飞行器研制的关键技术。作为热防护结构设计和材料选取的设计依据，高超声速飞行器表面的气动热问题一直是高超声速技术研究的重点和关键技术。因此，对飞行器气动热环境准确预测是气动热领域的重要基础工作，对防热结构设计的合理性和高超声速飞行器的飞行安全有着十分重要的意义。
气动热预测常用的手段之一是数值计算，数值计算的过程其实就是求解n-s方程，其预测结果精度高，但是由于计算过程的复杂性，耗时时间长，耗费大量计算资源，在某些复杂飞行器设计过程中，会扩大研制周期。因此，发明一种计算速度快、精度损失很小的气动热预测新方法对飞行器设计有重要意义。
机器学习是一种数据驱动的预测方法，并在图像、自然语言处理领域取得了显著成效。其中人工神经网络在求解非线性物理方程方面具有极大的优势，在一些领域已经获得成功。将人工神经网络应用在数值计算当中，可以代替传统求解n-s方程的方法，极大的缩短计算时间，同时又保证了计算精度。因此本发明能快速、高效地进行高超声速飞行器气动热预测，对高超声速飞行器研制意义重大。
技术实现要素：
本发明要解决的问题在于：数值计算具有计算时间长、计算资源耗费量大的特点。
为解决上述问题，本发明对不同算例的数值计算数据进行分类，推导对应的n-s方程，使用计算bp神经网络智能求解n-s方程，得到最终的预测结果。具体方案如下：
第一方面，本发明对所有算例数值计算得到的数据集，按照空间维度，对数据集分类，生成不同维度的特定数据集，存储在数据库中。
第二方面，对不同维度的特定数据集计算所采用的n-s方程进一步推导，得到方便构建网络的n-s方程。
第三方面，利用推导后的n-s方程构建计算bp神经网络，并利用不同维度的数据集训练网络，本发明设计的计算bp神经网络结构如图2，包括：
两个网络层和一个模块，两个网络层分别是输入层和输出层；一个模块是物理特征学习模块，根据数据处理的先后次序，其网络架构由底向上分别是输入层、物理特征学习模块、输出层。
第四方面，结合第三方面构建的计算bp神经网络得到智能计算结果，完成气动热预测。从以上技术方案以及实验结果发现，本发明实施例具有以下优点：
对气动热数值计算数据进行预先的分类处理可以方便后续流程进行，降低对最终预测结果的影响，然后对计算时所用的n-s方程进行推导后，构建计算bp神经网络智能求解n-s方程，最终输出气动热预测结果。训练好的计算bp神经网络具有良好的自适应性，而且预测精度高，在计算速度和计算资源消耗方面取得了突破性的成果，节约了大量的时间和成本。
附图说明
为了更清楚地说明本发明实施例中的技术发明方案，下面将对实施例或现有技术描述中所需要使用的附图做简单说明，显而易见，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为气动热预测流程图。
图2为计算bp神经网络结构图。
图3为实施方案气动热预测流程图。
图4为球头飞行器计算网格。
具体实施方式
发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
如图3所示，本实施例的主要步骤包括：第一步，获取数值计算数据；第二步，推导n-s方程；第三步，构建计算bp神经网络；第四步，智能求解n-s方程；第五步，输出预测结果具体实施步骤如下：
步骤1，获取球头飞行器模型的数值计算数据；
步骤1.1：确定球头模型的来流条件和边界条件，来流条件如下：
模型所处流场为层流，飞行马赫数为2ma，飞行器攻角为0度，雷诺数为1.657e+05,等壁温且壁温为300k。
步骤1.2：n-s方程可以描述流体流动情况，通过求解可压缩守恒型n-s方程来计算模拟球头模型表面的气动热环境，可压缩守恒型n-s方程具体形式如下：
步骤1.3：利用数值计算方法对可压缩守恒型n-s方程进行变换，以方便编写计算程序，变换后的方程形式如下：
其中：
步骤1.4：划分计算网格，网格结构如图4，编写计算程序求解n-s方程，得到球头表面的气动热环境的数值计算数据。
步骤2，推导n-s方程；
步骤2.1:步骤1中的n-s方程形式对于数值计算非常适合，但是不适合bp神经网络的搭建，因此，针对球头飞行器模型，推导n-s方程，具体推导过程如下：
连续介质运动方程、流体的本构方程、应变率张量具体形式如下：
联立方程消去应力张量、应变率张量，并设粘度为常数，得到如下方程：
其中，δ为拉普拉斯算子，具体形式为：
将式(4)代入欧拉运动方程，可得到：
将上述方程在x、y方向分离可得
步骤3，根据步骤2推导的n-s方程，构造计算bp神经网络；
步骤3.1：构造计算bp神经网络，其组成结构如图2所示，bp神经网络结构包括：输入层、物理特征学习模块、输出层，其中，物理特征学习模块由多个隐藏层组成，每个隐藏层由50个神经元组成，通过多个隐藏层提取输入的抽象特征，学习输入中所包含的物理特征，方便后续智能计算。
步骤3.2：利用球头模型数值计算数据训练计算bp神经网络，过程如下：
将数据集中的坐标、压力、温度、温度变量提取出来，统一转换成矩阵，作为训练集，坐标作为测试集。对网络进行训练和测试，测试误差5.2e-3。
步骤4，智能求解n-s方程；
步骤4.1：将需要求解的坐标输入到计算bp神经网络中，智能求解n-s方程。
步骤5，将智能求解n-s方程的结果输出，得到最终的气动热智能计算结果。