基于机器学习的结冰风洞温度场预测

2022-12-09张兴焕张平涛彭博易贤

实验流体力学 2022年5期

张兴焕，张平涛，彭博，易贤, *

1. 西南石油大学计算机科学学院，成都 610500

2. 中国空气动力研究与发展中心结冰与防除冰重点实验室，绵阳 621000

0 引言

结冰风洞是开展飞机结冰研究的重要地面试验设备，在飞机防/除冰系统研制等领域中发挥着重要作用[1–3]。其制冷系统通过控制换热器内制冷剂的流量和蒸发压力实现换热器出口气流温度控制，为试验提供精确稳定的低温环境。温度控制性能直接影响风洞热流场品质和试验效率。试验时需要根据不同工况进行压缩机吸气压力的精确预测和在线调整。影响制冷系统换热器换热效率和出口气流温度的因素较多，控制过程复杂，气流温度控制面临多参数耦合问题，研究人员难以建立准确的数学或物理模型来实现换热器的高效换热和气流温度的准确控制。因此，为实现对温度的精确控制，有必要开展压缩机吸气压力的快速准确预测以及风洞降温速率变化规律的研究。

为实现制冷系统气流温度场的精确控制并在此基础上提高制冷系统降温效率，常采用传统基于流场的计算方法、降温试验和工程计算等方法。传统计算方法需要考虑的影响因子较多，计算尺度大，建立的数学或物理方程复杂且换热或控制方程求解困难；降温试验通过在线调整压缩机吸气压力等相关参数的方式来准确控制气流温度；工程计算通过建立仿真模型的方式研究温度场控制方法，从而总结出控制规律。由于结冰风洞结构与换热器换热过程复杂，目前工程计算方法和实时降温试验存在过程复杂、耗时长、成本高等不足。而机器学习方法是以数据为驱动，利用大量试验数据经过不断训练调优总结出风洞内的温度控制规律，在复杂过程控制领域有着应用难度低、计算量小、响应速度快等显著优势，对复杂过程建立数学模型较为容易。目前，机器学习方法在结冰预测[4–9]和流场预测控制[10]领域有少量应用，但在制冷系统研究领域鲜有相关报道。

本文采用自适应粒子群优化（Adaptive Particle Swarm Optimization，APSO）算法、支持向量回归（Support Vector Regression，SVR）算法建立压缩机吸气压力预测模型；利用多层感知机（Multi−Layer Perceptron，MLP）神经网络建立分析模型以研究试验工况参数对结冰风洞降温速率的影响。以结冰风洞中的制冷系统为研究对象，开展了压缩机吸气压力预测和风洞降温速率变化规律研究。

1 大型结冰风洞及其制冷系统

本文基于中国空气动力研究与发展中心（CARDC）3 m×2 m 结冰风洞开展研究。该风洞是一座闭口回流式高亚声速风洞（图1），主要由结冰喷雾系统、制冷系统、高度模拟系统和风机动力系统组成[11-12]，可实现云雾参数模拟、热流场模拟、压力模拟等功能。结冰喷雾系统利用喷雾耙和喷嘴产生结冰云雾，主要模拟飞行器穿越云层飞行时的云雾环境。

图1 3 m×2 m 结冰风洞[11]Fig. 1 The CARDC icing tunnel

制冷系统对结冰风洞热流场的模拟至关重要。其主要由螺杆压缩机组、换热器、冷凝器、氨泵、储液器和低压循环桶组成，管路流程如图2 所示。采用氨液（R717）作为制冷剂，氨液由低压循环桶经氨泵加压后输送至换热器。试验最低气流温度可达–40 ℃，结冰降温速率直接决定了试验效率。风洞校测数据表明：高度模拟系统、风机动力系统会根据试验工况通过改变换热器入口气流场的气流压力和速度，影响换热器的传热过程，从而影响结冰风洞降温速率。

图2 制冷系统管路流程图Fig. 2 Refrigeration system pipeline flow chart

制冷系统的温度控制结构如图3 所示，包括外环和内环两个控制环。外环以压缩机控制器作为主控器，以距蒸发器最近处的低压循环桶上的压力传感器测得的压力值作为反馈，用低压循环桶实际压力与控制压力的差值控制压缩机的能级，从而实现风洞内对应目标温度的冷量控制。内环采用蒸发器回气电动控制阀作为控制器，通过控制蒸发器内的回气压力值实现对温度的精确控制。

图3 温度控制结构图Fig. 3 Temperature control structure diagram

2 压缩机吸气压力预测

为快速获得给定工况下压缩机吸气压力控制值，提高温度控制精度和自动化程度，需对其进行预测。本文采用APSO、SVR 算法进行压力预测，并开展试验对比。

2.1 样本数据

样本数据来源于3 m×2 m 结冰风洞2021 年5 月26 日至2021 年8 月14 日共71 天的试验数据。经过筛选和处理，一共获得了143 个有效样本数据。每个样本数据包括5 个参数，分别为换热器入口温度、换热器出口温度、气流压力、试验风速和压缩机吸气压力。

2.2 预测算法

2.2.1 APSO 算法

粒子群优化（Particle Swarm Optimization， PSO）算法是由Kennedy[13]和Eberhart[14]等以鸟群捕食行为为参照提出的一种基于群体智能的全局随机搜索算法。在搜索过程中，粒子不断更新自己的速度和位置。为了改善算法的收敛性能，在速度进化方程中引入惯性权重。粒子的速度进化和位置更新表示方法分别如式（1）、（2）所示：

为提高该算法的全局搜索和局部探索性能，避免陷入局部最优解，采用在PSO 算法基础上引入变异思想的APSO 算法，通过线性递减惯性权重[15]的方法对权重因子进行变异更新，优化公式如下：

式中：wmax和wmin分别为权重因子的最大值和最小值，分别设为0.9 和0.4；kmax为最大迭代次数。APSO算法流程如图4 所示。

图4 APSO 算法流程图Fig. 4 Flow chart of APSO algorithm

2.2.2 SVR 算法

SVR 算法是由支持向量机（Support Vector Machine，SVM）算法发展而来，是SVM 算法在回归问题中的应用。SVM 是由Cortes 和Vapnik[16]提出的基于统计学理论的原理性算法，可用于线性和非线性回归问题，是一种比较好地实现了结构风险最小化思想的算法。SVM 算法的基本思想是通过非线性变换将实际问题转换到高维的特征空间, 在高维空间中构造线性决策函数来解决非线性问题，使得输入空间中的决策超曲面模型与特征空间中的决策超平面模型对应[17]，如图5 所示，图中ε为模型预测输出与真实输出的偏差。超平面的表达式如下：

图5 SVR 算法示意图Fig. 5 The schematic diagram of SVR

式中：w为权重向量，决定了超平面的方向；x为非线性函数；b为阈值常量，决定了超平面与原点之间的距离。

与传统的预测方法相比，SVR 算法克服了局部最小化、过度学习以及过度依赖经验等问题，在预测结果的精确性和一致性上有明显优势。除此之外，SVR 算法还是一种有坚实理论基础的小样本学习方法，在小样本学习上具有较强的鲁棒性。SVR 问题的目标函数可以形式化为：

2.3 建立预测模型

为克服SVR 算法存在的参数选择敏感问题，本文采用APSO 算法对SVR 进行优化（APSO–SVR），建立压缩机吸气压力的预测模型。APSO 算法简单易行，设置参数少，具有较好的全局搜索能力。模型的输入为气流压力p1、试验风速v风、换热器入口温度Tin、换热器出口温度Tout，输出为压缩机吸气压力p吸。其中，p1值由模拟高度（H）转换而来，转换公式如下：

本文使用径向基函数（Radial Basis Function，RBF）作为SVR 算法的核函数。另外，针对正则化常数C、核函数系数γ以及ε的选取，采用APSO 算法搜寻最优组合来建立最终的APSO–SVR 预测模型以提高模型泛化能力。

采用APSO–SVR 预测模型对压缩机吸气压力预测的基本步骤如下：

步骤1 划分训练数据并验证数据的比例为130∶13，对训练集添加均值为0、标准差为1 的随机高斯噪声，再分别对输入、输出数据做标准化处理，表达式如下：

式中：xij，nor为第i个样本的第j个特征标准化后的值，xij为第i个样本第j个特征的真实值，xj，mean为样本中第j个特征的均值，xj，std为样本中第j个特征的标准差。

步骤2 初始化APSO 算法的参数，随机生成一组粒子的初始位置和速度，设置SVR 算法各个参数的范围。

步骤3 建立最初的SVR 预测模型，采用R2分数作为适应度函数值。R2分数既考虑了预测值与真实值之间的差异，也考虑了问题本身真实值之间的差异，是一个归一化的度量标准，其计算函数如下：

步骤6 将步骤5 得到的最优C、γ和ε值带入SVR 算法进行训练，得到最终的模型。

步骤7 对步骤6 得到的模型在测试集上进行测试得到预测值，并对预测值进行反归一化。

2.4 误差与预测结果分析

采用绝对百分比误差（Absolute Percentage Error， EAP）、平均绝对百分比误差（Mean Absolute Percentage Error， EMAP）和均方误差（Mean Squared Error， EMS）作为评价指标，表达式为：

随机选取143 个样本数据之外的13 组工况数据作为预测模型的验证集。每个验证样本数据和训练样本数据相同，都包含有4 个相同的工况条件并且都进行了标准化处理。结果表明：APSO–SVR 预测模型在验证集上的EMAP低于4%，EMS低于0.003。

表1 是从13 组验证工况中随机选取的6 组工况的压缩机吸气压力试验值（p吸，试）和预测值（p吸，预）。从表1 可以看出，该模型在随机选取的6 组验证工况下的EAP较小，预测值接近试验值。

表1 压缩机吸气压力试验值和预测值Table 1 Compressor suction pressure test value and prediction result

图6 为该预测模型在13 组验证工况下的预测值与试验值对比。从图中可以看出：预测值紧密分布在试验值周围，说明该预测模型是有效的。

图6 APSO-SVR 模型预测值与试验值对比Fig. 6 The comparison between predicted result of APSO-SVR model and experimental result

3 试验工况参数对降温速率的影响

上述吸气压力预测工作为温度稳态控制提供了基础，在风洞试验中得到了较好的应用。制冷系统结构复杂、换热过程影响因素多元，为进一步提高制冷系统降温效率，本文利用MLP 神经网络模型，开展试验工况参数对制冷系统降温效率影响的研究。

3.1 样本数据

构建MLP 神经网络模型的数据与压缩机吸气压力预测模型数据来源相同，经过筛选得到了143 个样本数据，但数据参数选择标准不同，每个样本数据包含7 个参数，分别是压缩机吸气压力、试验风速、气流压力、换热器入口温度、换热器出口初始温度、特定工况下的降温持续时间以及目标温度。

训练数据、测试数据和验证数据的比例为104∶26∶13。训练前，对训练数据添加均值为0、标准差为1 的随机高斯噪声来扩大训练集。为加快模型的收敛速度，对输入、输出数据分别进行标准化。

3.2 MLP 神经网络

MLP 神经网络是由多层M–P 单元组成的[17]一类广泛应用反向传播（Back Propagation， BP）算法的前馈神经网络，包含输出层、隐藏层和输出层，不同层之间为全连接[20]。

MLP 神经网络由多层BP 网络组合而成。与单层BP 网络相比，MLP 神经网络具有较强的自适应性和容错性，可解决复杂非线性问题，可用范围更广，训练出的模型鲁棒性更强，模型在已有数据样本基础上学习到的信息也更加丰富。

3.3 网络模型

采用定性分析法研究结冰风洞降温速率变化规律，利用MLP 神经网络搭建分析模型，网络输入为换热器出口初始温度T0、试验风速v风、气流压力p1、换热器入口温度Tin、特定工况下的降温时间t，网络输出为降温的目标温度Tout。

MLP 神经网络结构如图7 所示，图中H1、H2 分别表示第一个隐藏层和第二个隐藏层。可以看出，该网络模型有2 个隐藏层，隐藏层神经元个数分别为14 和7。激活函数使用softplus，相比于早期的激活函数，其更接近脑神经元的激活模型。Softplus的表达式如（13）所示。设置学习率为0.0001，使用EMS作为损失函数。训练过程中使用早停（early stopping）策略防止网络训练次数过多导致模型出现过拟合。

图7 MLP 网络结构图Fig. 7 The network structure of MLP

3.4 结果分析

采用EMAP作为MLP 神经网络模型评价指标。随机选取13 组工况作为降温速率变化规律分析验证模型的验证集。结果表明：模型的EMAP在15%左右。图8 为预测模型在13 组验证工况下的预测值与试验值对比。从图中可以看出：预测值紧密分布在试验值周围。说明该模型有较好的预测能力，为后续对降温速率随试验工况参数变化规律的准确分析提供了良好基础。

图8 MLP 模型预测值与试验值对比Fig. 8 The comparison between predicted result of MLP model and experimental result

3.4.1 消融试验

为研究输入变量对降温速率的影响程度，设置了消融试验，试验结果如表2 所示。其中，MLP–v风、MLP–p1、MLP–p吸和MLP–T0分别表示在训练数据中去除试验风速、气流压力、压缩机吸气压力和换热器出口初始温度训练得到的模型。

表2 消融试验结果Table 2 The result of ablation experiment

由表2 可以看出，各个输入变量对风洞降温速率的影响程度为：压缩机吸气压力>试验风速>初始温度>气流压力。

3.4.2 降温速率变化规律分析

研究降温速率[vcooling=（Tout−T0）/t]随试验工况参数的变化规律时，模拟高度分别设置为0.7、1、3 和6 km，即气流压力分别为92.65、83.88、69.64 和46.37 kPa，试验风速设置为20～140 m/s。图9、10 是在压缩机吸气压力分别为191 和202 kpa 2 个工况下降温速率随气流压力和试验风速的变化曲线。从图9、10 可以看出：试验风速为40～60 m/s 时，降温速率相对最大；当风速大于60 m/s 时，降温速率呈快速下降趋势。这是因为随着试验风速增大，风机功率增大，给风洞系统带来的热负荷增大，导致在低风速条件下降温速率较大，而在高风速条件下降温速率较小。