基于深度学习的非定常周期性流动预测方法
发布时间:2021-10-20 00:06
为了克服传统CFD计算需要耗费大量的计算时间与成本的缺陷,提出了一种基于深度学习的非定常周期性流场的预测框架,可以实时生成给定状态的高可信度的流场结果。将条件生成对抗网络与卷积神经网络相结合,改进条件生成对抗网络对生成样本的约束方法,建立了基于深度学习策略采用改进的回归生成对抗网络模型,并与常规的条件生成对抗网络模型的预测结果进行对比。研究表明,基于改进的回归生成对抗网络的深度学习策略能准确预测出指定时刻的流场变量,且总时长比CFD数值模拟减少至少1个量级。
【文章来源】:空气动力学学报. 2019,37(03)北大核心CSCD
【文章页数】:8 页
【部分图文】:
图1回归生成对抗网络训练流程图Fig.1TrainingprocessofregressionGAN
误差作为损失函数,作为生成器损失函数的一部分,用于约束生成器的优化方向。训练过程中不断重复上述三部分工作,最终可以得到有良好预测能力的生成器。测试过程仅需进行第一部分,直接输入对应的时间标签y,即可得到需要预测时刻的流场。1.3样本集采用CFD求解器对Re=1×104、Ma=0.2[13]、圆柱直径D=1m的圆柱绕流进行了数值模拟。通过RANS方程求解,求解器为NASACFL3D[14]。计算网格如图2所示,网格尺寸为128×128(周向×法向),附面层第一层高度设置为5×10-4,y+=0.3,远场高度为50倍特征长度(圆柱直径)。图2计算与训练网格Fig.2Sketchofcomputation/traininggrid圆柱绕流是以时间T为一周期的周期性流动,本文以T/40作为一个时间间隔,去除t=0、t=T时刻的流场数据,以剩余的39个时刻的流场变量u、v、p作为训练数据。标签y是该样本在T中对应的时间t。因为神经网络参数通过梯度下降方法[15]进行优化,为了令参数优化更加准确、加速梯度收敛,训练前对数据进行标准化处理。设共有m个训练样本,以流场变量u为例:u-=u-μuσu(2)μu=1m∑mi=1ui(3)σu=1m∑mi=1(ui-μu)槡2(4)v、p与u的标准化方法相同。即训练样本集由38个尺寸为
始kt=0,权重?=0.97;γ=0.5为平衡超参数,其值越小表示生成样本的多样性越低;t为训练步数。kt决定D(G(z|y))所占权重,防止D或G相较另一方能力太强,出现过拟合或模型坍塌。1.5网络结构参数使用的回归生成对抗网络结构设计参数为:生成器以隐变量(随机噪声z=128×1)与真实标签(y=1×1)叠加作为输入,即输入数据尺寸为129×1。使用的生成器结构为多层反向卷积神经网络,其框架如图3所示。训练过程中,生成器生成尺寸为128×128×3(周向点数×法向点数×流场变量个数)的生成样本G(z|y)。生成器模型结构参数见表1。(a)生成器/解码器(b)编码器图3自动编码器网络结构Fig.3Designedautoencodernetworkarchitecture表1生成器/解码器网络结构与参数设置Table1Generator/decodernetworkusedinexperiments层类型核激活函数隐变量输入(h=128+1)--隐含层1卷积3×3×128ELU隐含层2卷积3×3×128ELU上采样1-2×2-…(隐含层1~上采样1重复3次)………隐含层9卷积3×3×128ELU隐含层10卷积3×3×128ELU隐含层11卷积3×3×3ELU判别器结构为自动编码器,分为编码器和解码器两部分,两部分的结构均可视为卷积神经网络,如图3所示。判别器的输入为生成器预测的流场
【参考文献】:
期刊论文
[1]基于特征正交分解的跨声速流场重构和翼型反设计方法研究[J]. 刘浩,徐敏,叶茂. 空气动力学学报. 2012(04)
[2]模拟退火算法和POD降阶模态计算在翼型反设计中的应用[J]. 赵松原,黄明恪. 空气动力学学报. 2007(02)
本文编号:3445859
【文章来源】:空气动力学学报. 2019,37(03)北大核心CSCD
【文章页数】:8 页
【部分图文】:
图1回归生成对抗网络训练流程图Fig.1TrainingprocessofregressionGAN
误差作为损失函数,作为生成器损失函数的一部分,用于约束生成器的优化方向。训练过程中不断重复上述三部分工作,最终可以得到有良好预测能力的生成器。测试过程仅需进行第一部分,直接输入对应的时间标签y,即可得到需要预测时刻的流场。1.3样本集采用CFD求解器对Re=1×104、Ma=0.2[13]、圆柱直径D=1m的圆柱绕流进行了数值模拟。通过RANS方程求解,求解器为NASACFL3D[14]。计算网格如图2所示,网格尺寸为128×128(周向×法向),附面层第一层高度设置为5×10-4,y+=0.3,远场高度为50倍特征长度(圆柱直径)。图2计算与训练网格Fig.2Sketchofcomputation/traininggrid圆柱绕流是以时间T为一周期的周期性流动,本文以T/40作为一个时间间隔,去除t=0、t=T时刻的流场数据,以剩余的39个时刻的流场变量u、v、p作为训练数据。标签y是该样本在T中对应的时间t。因为神经网络参数通过梯度下降方法[15]进行优化,为了令参数优化更加准确、加速梯度收敛,训练前对数据进行标准化处理。设共有m个训练样本,以流场变量u为例:u-=u-μuσu(2)μu=1m∑mi=1ui(3)σu=1m∑mi=1(ui-μu)槡2(4)v、p与u的标准化方法相同。即训练样本集由38个尺寸为
始kt=0,权重?=0.97;γ=0.5为平衡超参数,其值越小表示生成样本的多样性越低;t为训练步数。kt决定D(G(z|y))所占权重,防止D或G相较另一方能力太强,出现过拟合或模型坍塌。1.5网络结构参数使用的回归生成对抗网络结构设计参数为:生成器以隐变量(随机噪声z=128×1)与真实标签(y=1×1)叠加作为输入,即输入数据尺寸为129×1。使用的生成器结构为多层反向卷积神经网络,其框架如图3所示。训练过程中,生成器生成尺寸为128×128×3(周向点数×法向点数×流场变量个数)的生成样本G(z|y)。生成器模型结构参数见表1。(a)生成器/解码器(b)编码器图3自动编码器网络结构Fig.3Designedautoencodernetworkarchitecture表1生成器/解码器网络结构与参数设置Table1Generator/decodernetworkusedinexperiments层类型核激活函数隐变量输入(h=128+1)--隐含层1卷积3×3×128ELU隐含层2卷积3×3×128ELU上采样1-2×2-…(隐含层1~上采样1重复3次)………隐含层9卷积3×3×128ELU隐含层10卷积3×3×128ELU隐含层11卷积3×3×3ELU判别器结构为自动编码器,分为编码器和解码器两部分,两部分的结构均可视为卷积神经网络,如图3所示。判别器的输入为生成器预测的流场
【参考文献】:
期刊论文
[1]基于特征正交分解的跨声速流场重构和翼型反设计方法研究[J]. 刘浩,徐敏,叶茂. 空气动力学学报. 2012(04)
[2]模拟退火算法和POD降阶模态计算在翼型反设计中的应用[J]. 赵松原,黄明恪. 空气动力学学报. 2007(02)
本文编号:3445859
本文链接:https://www.wllwen.com/kejilunwen/lxlw/3445859.html