基于深度学习的非定常周期性流动预测方法

基于深度学习的非定常周期性流动预测方法摘要：为了克服传统CFD计算需要耗费大量的计算时间与成本的缺陷，提出了一种基于深度学习的非定常周期性流场的预测框架，可以实时生成给定状态的高可信度的流场结果。将条件生成对抗网络与卷积神经网络相结合，改进条件生成对抗网络对生成样本的约束方法，建立了基于深度学习策略采用改进的回归生成对抗网络模型,并与常规的条件生成对抗网络模型的预测结果进行对比&研究表明，基于改进的回归生成对抗网络的深度学习策略能准确预测出指定时刻的流场变量，且总时长比CFD数值模拟减少至少1个量级&关键词：深度学习;卷积神经网络；生成对抗网络；回归；非定常流场;预测Amethodofunsteadyperiodicflowfieldprediction

basedonthedeeplearningAbstract:Inordertoovercometheshortagesofthecomputationallyexpensiveandtime!consumingiterativeprocessintraditionalCFDsimulation，aframeworkbasedonthedeeplearningtopredictperiodicunsteadyflowfieldisproposed，whichcanaccuratelypredictreal-timecomplexvortexflowstateatdifferentmoments.TheconditionalgenerativeadversarialnetworkandconvoAutionaAneuraAnetworkarecombinedtoimprovetheconditionaconstraintmethodfromconditionalgenerativeadversarialnetwork.Theimprovedregressiongenerativeadversarialnetworkbasedonthedeeplearnng;sproposed.Thetwoscenarosofcondtonalgeneratveadversaralnetworkandregressongeneratveadversaralnetworkaretestedandcomparedvagvngdfferentperodcmomentstopredctthecorrespondngflowfeldvarables.Thefnalresultsdemonstratethatregressongeneratveadversaralnetworkcanestmatecomplexflowfelds，and;sfasterthantradtonalCFDsmulatonoveroneorderofmagntudes.!Peywords:deeplearnng$convolutonalneuralnetwork$generatveadversaralnetworks$regresson$unsteadyflow$predcton在航空领域中进行与气动相关的设计或优化时，经常使用CFD数值模拟。但CFD往往需要消耗大量的计算资源与时间，在日益追求效率的现今时代,许多学者都对如何节省计算成本进行了大量研究&例如借助降阶模型，提取流场主要特征来降低计算时间。但降阶模型对流场参数变化非常敏感⑵，缺乏足够的鲁棒性，即使是参数微小的变化，也容易导致降阶模型精度降低。近年来，深度学习技术强大的学习隐藏特征的能力，使其在图像识别等领域得到广泛应用与发展$这也为如何进行流场计算提供了新的思路$目前，深度学习在流体力学中的相关应用仍不多见$2016年，美国Sandia国家实验室Ling等，通过将伽利略不变量嵌入深度神经网络结构，构建了对雷诺平均(ReynoldsAveragedNavier-Stokes,RANS)湍流模型的深度学习，对雷诺应力的进行了预测，并以此为结果对管道流和速度场进行了预测$Myanawala等人(2017)首次实现了深度卷积网络通过二维柱体尾迹区流场变量对气动载荷的准确预测$他们仅采用3组不同几何外形柱体的全阶模拟数据来训练卷积神经网络(ConvolutionalNeuralNetworks,CNN)[5]，就成功预测了其他10组不同几何外形柱体的流体载荷"6#$Lee等人(2017)采用生成对抗网络，来学习全阶N!S方程(Navier-StokesEquations)的非定常数值解，通过对两个雷诺数下二维圆柱绕流全阶数值解的学习，就能够准确预测出介于两个训练雷诺数之间的二维圆柱尾流融合的非定常过程"7#$2018年,Lee通过雷诺数100〜400的流场样本集，成功预测了未参与训练的雷诺数为500与3000的流场"8#$上述例子充分说明深度学习具有学习和预测复杂流动现象的强大潜力$然而，现有的研究结果仍然存在一定的局限性，例如Myanawala"6#提出的卷积神经网络框架仅仅预测了圆柱绕流的气动力系数但并未直接预测非定常流场，Lee"7-8」提出的生成对抗网络框架必须将预测时刻的前几个时刻的流场作为输入信息，在事先没有初始流场结果时,很难开展相应的流场预测工作。本文旨在实现仅需给出时间条件，即可基于深度学习预测对应非定常流场。以非定常流场数值仿真结果作为训练样本，将条件生成对抗网络与卷积神经网络相结合，改进条件约束方式。建立了回归生成对抗网络，并比较了本文的约束方式与传统方式的效果$1回归生成对抗网络为了能够生成考虑时间条件约束的非定常流场预测结果，基于条件生成对抗网络(ConditionalGenerativeAdversarialNets,CGAN)"D#,提出了改进的回归生成对抗网络(RegressionGAN)，并将常规的条件生成对抗网络与改进的回归生成对抗网络分别用于二维圆柱的非定常流场预测，进行了结果对比$1.1条件生成对抗网络条件生成对抗网络是一种考虑条件约束的生成对抗网络模型$相比于传统的"无控型"生成对抗网络，条件生成对抗网络可以生成给定条件约束下的预测结果$条件生成对抗网络与传统生成对抗网络"10」的区别是在生成器和判别器中均引入了条件y$将输入噪声z与条件y组合，作为生成器新的输入，组合方式多种多样，可选方式很多"11#；同样的，将原本的输入x和y作为判别器的新输入$这样我们可以通过改变y来指定输出信息，其数学表达式为：minmaxV(D,G)=EG〜0dataG)[lnD(x|y)#+GDEG~"ln(1—D(G(z|y)))# (1)其中,G(z|y)、D(G(z|y))分别为生成器、判别器输出数据,Pdaa表示训练样本分布,6L表示噪声分布,E代表期望值$1.2回归生成对抗网络训练流程!回归生成对抗网络对条件生成对抗网络进行了改进，增加回归器对生成器的优化方向进行约束，其训练流程如图1所示，主要分为三部分：生成器生成流场变量G(z|y)$为了学习流场变量x的分布P,我们定义噪声z的先验分布PZ),与条件生成对抗网络相同，生成器将z|y映射至G(z|y)，即为预测得到的流场变量$判别器判别数据真假与自动编码$本文的判别器采用自动编码器E#的网络结构，该方法在基于能量的生成对抗网络"12#中首次提出，主要作用是判断生成数据的真实性，以及自动编码真实数据与生成数据$以G(z|y)与x分别作为判别器的输入，在训练过程中，D(G(z|y))越小，代表G(z|y)与真实数据x越接近，判别器自动编码数据的能力越强；同理,D(x)越小，判别器自动编码数据的能力越强$图1回归生成对抗网络训练流程图Fig.1TrainingprocessofregressionGAN回归生成对抗网络在条件生成对抗网络的基础上增加回归器，用于回归时间标签R(G(z|y))$对于回归器来说，只有G(z|y)足够准确，才能够正确预测出对应的时间R(G(z|y))$为了将回归器

的效果反馈加入训练，以R(G(<Iy))与y之间的误差作为损失函数，作为生成器损失函数的一部分，用于约束生成器的优化方向$训练过程中不断重复上述三部分工作，最终可以得到有良好预测能力的生成器$测试过程仅需进行第一部分，直接输入对应的时间标签y，即可得到需要预测时刻的流场$1.3样本集采用CFD求解器对Re(1X10*、Ma(0.2&13'、圆柱直径D=1m的圆柱绕流进行了数值模拟。通过RANS方程求解，求解器为NASACFL3D&1*'$计算网格如图2所示，网格尺寸为128X128(周向X法向)，附面层第一层高度设置为5X10!*,》+(。.3,远场高度为50倍特征长度(圆柱直径)$图2计算与训练网格Fig.2Sketchofcomputation/traininggrid圆柱绕流是以时间T为一周期的周期性流动,本文以T/40作为一个时间间隔，去除t(0、t(T时刻的流场数据，以剩余的39个时刻的流场变量7F、6作为训练数据$标签y是该样本在T中对应的时间£$因为神经网络参数通过梯度下降方法［15］进行优化，为了令参数优化更加准确、加速梯度收敛，训练前对数据进行标准化处理$设共有5个训练样本,以流场变量7为例：171&7F、6与7的标准化方法相同$即训练样本集由38个尺寸为128X128X3(周向点数X法向点数以流场变量7为例：171&7X*7,F,6+y={'}(5)(5)(6)回归生成对抗网络的损失函数为：7Ld=D(x')—$D(G(<Iy))6lg(z)==D(G(zIy))+(1—=)L (7)［kl+1=$'J/k&)D(x)—D(G(<Iy))'L!="MSE(y,R(G(<Iy)))i-+"—R(G(<Iy))1 (8)其中:为一固定比例，取=0.001；初始$,=0,权重==0.97；)=0.5为平衡超参数，其值越小表示生成样本的多样性越低,'为训练步数。$决定D(G<Iy))所占权重，防止D或G相较另一方能力太强，出现过拟合或模型坍塌$1.5网络结构参数使用的回归生成对抗网络结构设计参数为：生成器以隐变量(随机噪声<(128X1)与真实标签(y(1)1)叠加作为输入，即输入数据尺寸为129X1$使用的生成器结构为多层反向卷积神经网络，其框架如图3所示$训练过程中，生成器生成尺寸为128X128X3(周向点数X法向点数X流场变量个数)的生成样本G(<\y)$生成器模型结构参数见表1$(a)生成器/解码器 (b)编码器图3自动编码器网络结构Fig.3Designedautoencodernetworkarchitecture层类型核激活函数隐变量输入0(128+1)--隐含层1卷积3X3X128ELU隐含层2卷积3X3X128ELU上采样1-2X2-(隐含层1〜上采样1重复3次层类型核激活函数隐变量输入0(128+1)--隐含层1卷积3X3X128ELU隐含层2卷积3X3X128ELU上采样1-2X2-(隐含层1〜上采样1重复3次)隐含层9卷积3X3X128ELU隐含层10卷积3X3X128ELU隐含层11卷积3X3X3ELU表1生成器/解码器网络结构与参数设置Table1Generator/decodernetworkusedinexperiments(2)(3)(4)表2编码器网络结构与参数设置Table2Encodernetworkusedinexperiments层类型核激活函数输入层输入(128X128X3)--隐含层0卷积3X3X128ELU隐含层1卷积3X3X128ELU隐含层2卷积3X3X128ELUNN降采样1-2X2-隐含层3卷积3X3X128ELU隐含层4卷积3X3X256ELUNN降采样2-2X2-隐含层5卷积3X3X256ELU隐含层6卷积3X3X384ELUNN降采样3-2X2-隐含层7卷积3X3X384ELU隐含层8卷积3X3X512ELUNN降采样4-2X2-隐含层9卷积3X3X512ELU隐含层10卷积3X3X512ELU隐含层11全连接层128-使用的回归器结构为卷积神经网络，如图4所示$生成器生成的流场变量作为输入，提取与映射输入数据的特征回归标签R!G<Iy))$回归器模型结构参数见表3$图4本文中使用的回归器网络结构Fig.4DesigEedregressornetworkarchitecture表!回归器网络结构与参数设置Table3Regressornetworkusedinexperiments层类型核激活函数输入层输入(128X128X3)--隐含层1卷积3X3X20ReLU隐含层2池化2X2-隐含层3卷积3X3X40ReLU隐含层4池化2X2-隐含层5卷积3X3X60ReLU隐含层6池化2X2-隐含层7卷积3X3X80ReLU隐含层8池化2X2-隐含层9全连接层1024-隐含层10全连接层1tanh一般卷积神经网络主要由卷积层、激活函数、上采样、降采样、池化层和全连接层等组成$卷积层$卷积层是卷积神经网络的核心组成部分，每个卷积层由若干个卷积核组成，使用卷积核的目的是提取数据的不同特征$激活函数$神经网络本质是线性函数的叠加，无法解决非线性问题。激活函数给神经元引入了非线性因素，使得神经网络可以逼近任意非线性函数$本文在自动编码器中的卷积层后都增加了指数线性单元(ExponentialLinearUnits,ELU)&1句激活函数，在回归器中使用了修正线性单元(RectifiedLinearUnt,ReLU)[17]和双切正切函数[18]$上采样。上采样的目的是放大数据尺寸，-般采用内插方法。插值方法有多种，本文采用最邻近插值一一在待求点的四个相邻点中，将距离待求点最近的点上的值直接赋予待求点$降采样。降采样的目的是缩小数据尺寸，获取不同尺度信息$例如对于尺寸为HXW的数据,对其进行；倍降采样，即得到(H/；)X(W/；)尺寸的数据$本文通过卷积层直接对数据进行降采样$池化层$池化层与降采样的作用相似，主要用于特征降维，减少运算数据量，在一定程度上防止过拟合$池化分为最大池化和平均池化。本文使用最大池化$全连接层$全连接层是由许多神经元组成的，起到将学习到的“分布式特征表示”映射到样本标记空间的作用$这样做的目的是高度提纯数据特征,减少特征位置带来的影响$2结果讨论生成的训练样本集总数为39个，搭配对应的周期性时间进行训练$优化器选择Adam&19'最速下降法，—&(0.5，-2(0.999，学习率设置为0.0001,每批输入2个样本进行训练，训练步数为1I0000步$训练完成后，将隐变量与未参与训练的39个标签(随机选取如'(50/4、'(6/4等时刻)进行叠加，输入回归生成对抗网络、条件生成对抗网络中的生成器进行流场预测，并挑选'(0、'(4/3、'(24/3、'(4时刻作为算例展示，如图5〜图8所示$图5〜图8,将回归生成对抗网络与条件生成对抗网络与真实样本进行比较，可见回归生成对抗网络预测效果较为准确，虽然在一些小尺度结构上仍有差别，但已经能够捕捉三个流场变量7F6的绝大部分特征；'(0、'(4/3与'(4时刻比'(24/3时刻的预测效果更好,'(24/3时刻可以明显观察到等值线略微抖动，而不是像真实样本中光滑连续$条件生成对抗网络出现模型坍塌现象，不同时刻的流场预测结果几乎完全一致，仅能够捕捉到单一时刻的流场特征。以上两项实验均表明，深度学习具有模拟复杂流动的潜力$此时N-S方程已经封闭，可通过7F6计算涡

量。涡量图相较速度云图与压力云图可以更加清晰地反映各个时刻不同的流动状态和旋涡的强度及方向，从涡量图中可以直接观察到涡在一个周期内的变化情况。图9展示了计算得到的涡量图。与预测流场变量UFP类似,回归生成对抗网络成功地捕捉到流场特征，与真实样本吻合良好。由于'(24/3时刻的预测效果相对较差，造成涡量虽能够反映流动特征，但涡的边界抖动，出现不符合实际的流动现象。(a)回归生成对抗网络预测结果(b)条件生成对抗网络预测结果(c)真实样本图5非定常周期性流动4=0流场变量预测结果对比(a)回归生成对抗网络预测结果Fig.5Comparisonbetweenthepredictedandgroundtruthflowfieldsat(a)回归生成对抗网络预测结果Fig.5Comparisonbetweenthepredictedandgroundtruthflowfieldsat4=0(b)条件生成对抗网络预测结果(c)真实样本图6非定常周期性流动4=)/3流场变量预测结果对比(a)回归生成对抗网络预测结果Fig.6Comparisonbetweenthepredictedandgroundtruthflowfieldsat4=(a)回归生成对抗网络预测结果Fig.6Comparisonbetweenthepredictedandgroundtruthflowfieldsat4=)/3(b)条件生成对抗网络预测结果(c)真实样本图7非定常周期性流动4=2)/3流场变量预测结果对比Fig.7Comparisonbetweenthepredictedandgroundtruthflowfieldsat4=2)/3(a)回归生成对抗网络预测结果(b)条件生成对抗网络预测结果(c)真实样本图8非定常周期性流动t=TFig.8Comparisonbetweenthepredictedandgroundtruthflowfieldsatt=Tduringtest(a)方(0时刻涡量云图(b)t(T/3时刻涡量云图(c)'(2T/3时刻涡量云图(d)(d)'(T时刻涡量云图图9非定常周期性流动涡量图PredictedvorticityfromregressionGAN,conditionalGANandtruevalueFig.9边界层信息对于计算气动力至关重要$本文收集了位于g(0、g(0.18D、g(0.36：处通过真实样本与回归生成对抗网络预测样本计算的边界层速度型(图10)$将真实样本与回归生成对抗网络预测结果计算的速度型作差，得到相对误差曲线，如图11$从图11中可以明显观察到，'(2T/3时刻的相对误差较其余3个时刻的误差更大，与前文的7、f、6预测效果相符$速度型预测的相对误差维持在4%以内，大部分误差不超过0.5%，反映出本文预测方法效果良好，可以正确捕捉边界层特征$本文统计了CFD计算时长、回归生成对抗网络训练与测试时长如表4(CPU-h表示CPU单核计算1h,GPU・h表示GPU单核计算1h)$可以发现，经过GPU加速，训练完成后，只须运行数秒即可生成所求流场变量，相比于CFD数值模拟，其计算时间至少减少一个量级，且可以非常方便地预测任意时刻的流场变量，大大降低了计算时间与成本，缩短计算周期$

GroundTruth:x=0GroundTruth:x=0.1GroundTruth:x=0GroundTruth:x=0.18DGroundTruth:x=0.36DPredicted:x=0Predicted:x=0.18DPredicted:x=0.36Z)(c)t(2T/3时刻(5)t(T时刻(a)'=%时刻GroundTruth:x=0GroundTruth:x=0.18DGroundTruth:x=0.36DPredicted:x=0(a)'=%时刻GroundTruth:x=0GroundTruth:x=0.18DGroun