跨音速自然层流翼型反设计研究

!第"#卷第#期郑州大学学报!理学版"$%&’"#(%’# !)*#+年,月-./012340%56278.!(9:.;<7.=>."?9@.)*#+
跨音速自然层流翼型反设计研究
李!静#!!刘艳萍#!!高正红)!!梁天水#
!#.郑州大学力学与工程科学学院!河南郑州E"***#$
).西北工业大学翼型叶栅空气动力国防科技重点实验室!陕西西安D#**D)"
摘要!基于V9LLM U e d降阶模型的数据重构方法提出反设计方法#由已知的翼型表面压力分布对缺失的目标翼型
数据进行重构.在初始快照采样过程中#采用!Y;P方法对初始翼型进行扰动取样#并用反设计得到翼型替换基础
扰动翼型#将新的快照加入原始快照库#反复迭代求解.针对层流翼型反设计#利用边界层转捩数值模拟程序获得
快照向量的翼型压力分布.接着#提出两种不同的跨音速自然层流翼型的压力分布形态#通过构建的反设计方法得
到两组不同的跨音速自然层流翼型#并详细研究了不同雷诺数及不同马赫数下#不同压力分布形态对转捩发生抑
制能力的影响#同时研究其与阻力发散特性的协调性.
关键词!Y;P$自然层流翼型$转捩模型$反设计
中图分类号!$)D#文献标志码!G文章编号!#CD#B CAE#!)*#+"*#B**+"B*C
!"#!#*’#,D*"Q R.7S S2.#CD#B CAE#’)*#A*A*
$%引言
由于现代民机的摩阻占总阻力的E*W i""W#增大机翼表面的层流区域以期降低摩阻成为设计空气动力学近些年的研究热点(#X,).目前#自然层流机翼已经在高空长航时无人机(E)’中小型公务机上得到了广泛的应用("X C)#但针对跨音速的大中型民机尚处于理论与试验阶段#尤其是后掠效应与实现层流对机翼表面压力分布的要求相互矛盾#更限制了自然层流技术在中大型民机上的实际应用(D).尽管如此#自然层流技术依然是提高未来民机性能的发展方向#而对于大多数在跨声速状态巡航的民机来说#对跨音速自然层流翼型的深入研究#必然为真正制造出高性能的层流机翼提供一定的理论依据及技术支持.
维持层流的顺压区长度与激波位置是跨音速自然层流翼型设计的关键#顺压有利于维持层流#但过长的顺压区为翼型后缘的压力恢复带来较强的激波#因此在气动设计中#需要兼顾这两个关键因素.本文旨在通过压力分布的反设计来协调激波的位置’强度以及顺压区#并通过考察翼型在不同雷诺数下及不同马赫数下的气动性能#为跨音速层流翼型的设计提供新方法和新思路.
反设计方法的优点是&在设计人员丰富气动设计经验的基础上#针对设计要求提出合理压力分布或速度分布#然后通过求解空气动力学反问题来确定对应的几何外形.近些年发展的特征正交分解!L@%L1@%@:0%3%B 29&>1S<%I L%S7:7%2#U e d"方法是一种高效的数据处理方法#已经广泛地应用到了诸多领域(A X+)#它通过模态分析的思想提取到数据的主要特征#用在翼型反设计的优势是在给定目标压力分布和对已知基准翼型形状修改而形成的翼型设计空间的条件下#通过基于U e d的数据重构方法!V9LLM U e d"就可以反设计得到与给定的目标压力分布相对应的最优翼型形状.
本文采用了一种基于差量采样解的V9LLM U e d翼型反设计方法#将该方法成功应用在跨音速层流翼型的反设计中#针对两种经典的压力分布形态#进行了精确高效的反设计#评估了压力分布形态对转捩发生抑制能力的影响.
&%基于差量采样解的U4T T G M"!方法
文献(#*)详细介绍了U e d方法#V9LLM U e d方法是U e d方法用于已知系统中的填补某一缺失数据的收稿日期!)*#AB*,B))
基金项目!国家自然科学基金项目!##C*)))C"$河南省高等学校重点科研项目!#CG#,***#".
作者简介!李静!#+A"%"#女#河南平顶山人#讲师#主要从事飞行器气动设计及Y F d研究#=B I97&&&7R723S1&\N445.1>5.<2.
郑州大学学报!理学版"第"#卷元素时的变形#通过在系统中获取合适的足够数量采样解!快照样本"#在采样解的基础上提取此系统的基
模态#采用V
9LLM U e d 方法(##)可以对该系统中任一缺失数据的元素进行精确填补.首先定义向量矩阵及其对应的.标签/向量矩阵G 1.标签/向量矩阵G 用来识别向量矩阵?中的元素
是否缺失1若不缺失7X ;-##若缺失7X ;-*#其中7X ;对应于向量矩阵?中第X 个不完全数据向量中的第;个
元素的.标签/1假设0,;19;-#是通过U e d 方法求得的已知采样集0@;1T ;-#的基向量
(#))1H 表示该系统中待填补的向量#由已经得到的9个U
e d 基来填补近似向量9H 1那么9H -#9
;-#";
,;#!#"
并通过求解(I j J 线性方程组来确定";#其中&(;G -!,;#,;"7$!;-!H #,;"71近似向量9H 在求出";后由公
式!#"求得1
本文研究了一种基于差量采样解的V
9LLMU e d 方法#对原始采样解的结构重新定义#重新定义采样解集合为!5-(!5#.25"#!5).25"#-#!5T .25")#其中&25表示采样解的平均值$5;T ;-#表示第;
个采样解1通过上述方法求解!5的基向量及系数#那么最终修正向量可表示为9H -25##9
;-#
";
,;1+%翼型反设计方法
基于差量采样解V
9LLM U e d 方法构建翼型反设计方法分E 步&#"根据反设计需要定义合理的设计空间#利用计算流体力学!<
%I L5:9:7%29&\&57>>M 29I 7<S #Y F d "方法获得设计空间中样本的解#以此构造快照矩阵$)"对生成的快照矩阵按文献(#*)中的方法求解基向量+#并计算每个基向量的能量$,"根据所求解问题的需要对低能量模态进行适当截断#形成对应设计空间的降阶子空间$E "利用降阶子空间对目标缺失数据进行数据重构.
+,&%8’!数值模拟技术
采用有限体积法求解雷诺平均(9871@B ;:%J1S 方程#其积分形式为--"<’V >’#=
)!,7>)-*#在预测转捩时#采用文献(#)X #,)提出的基于;;P X b #湍流模型中的完全基于流场当地变量的%b _1&"转捩模型#%b _1&
"两方程转捩模型通过%控制湍动能生成来实现当地转捩#无量纲守恒形式的%输运方程为
--%-"#--%@;-I ;---I ;(!1#1"+%
"-%-I ;)#2:#-)F &123:0F %2S 1:*’"%*’"!#.2<#%".2:)-*F :5@Z %!2<)%.#"1!!引进变量_1&"#边界层以外的_1&"通过经验关联函数_1&"计算#边界层之内的_1&
"通过输运方程从边界层以外联立计算#
--_1&"-"#--_1&"@;-I ;---I ;(+&"!1#1""-_1&"-I ;)#2&"-"S <9&1
!_1&"._1&""!#.F &""1+,+%快照获取方法
本文的快照矩阵是由翼型数据和压力分布构成的#每一个翼型数据和对应压力分布组成快照矩阵的列向量#通过在基准翼型的基础上加入扰动量形成多个快照#依次构成翼型反设计的设计空间.
对于基准翼型的扰动是通过!
Y ;P 翼型参数化方法实现的#对于在采用类函数变化!<&9S S \52<:7%2Q S 09L1\52<:7%2:@92S \%@I 9:7%2#Y ;P "参数化(#E )翼型时#定义9!5"N-5*’"
!#.5"#7
;-*
!H ;#!H ;"D ;5;!#.5"7.;#59Z #即在初始设计变量基础上引入一个扰动量!H ;
#将上式展开得9!5"N-9!5"#5*’"
!#.5"#7
;-*!H ;D ;5;!#.5"7.;#!)"
引入定义
!9!5"-5*’"
!#.5"#7
;-*!H ;D ;5;!#.5"7.;#!,"
因此#由式!)"和!,"得到9!5"N-9!5"#!9
!5"1C +
!第#期李!静#等$跨音速自然层流翼型反设计研究
+,-%反设计方法
本文的反设计过程为&#"对快照向量进行U e d分解#得到基向量#用基向量逼近目标压力分布得到反设计翼型$)"以反设计得到的翼型为初始翼型#再进行扰动得到扰动翼型#加上反设计得到翼型#进行Y F d分析获得快照向量#将这些快照向量加入原始快照向量库#如不满足条件#重复第#’)步$,"满足收敛条件#结束.
-%跨音速自然层流翼型反设计
本节主要通过两种不同的压力分布形态反设计得到气动外形#进而在不同雷诺数及马赫数下分析压力分布形态对转捩发生的抑制能力#同时研究其与阻力发散特性的协调性.
!Y;P参数化上下表面各#"个设计变量#共,*个设计变量#这,*个设计变量依次取指定的变量上下限#共产生C*个新翼型#加上基准翼型共C#个.采用文中构建的%b_1
&"
转捩模型进行边界层转捩预测和气动特性分析评估#计算网格为A万.
对生成的C#个翼型形状在如下状态进行Y F d计算#获得采样解&,:
q -*’DE$$
Q
-*’",$_1-C’*/
#*C$飞行环境湍流度Z@-*’)0$湍流黏性比为#*.
图#为目标压力分布反设计结果#可以看出来#本文的反设计方法对于不同的压力分布都可以得到令人满意的结果#记此次反设计得到的翼型为G7@F%7&K];Y).图)为初始基准翼型与两个反设计翼型的比较.选
取一个设计状态为,:
q -*’DE#$
Q
-*’"D"#_1-)’*/#*D#阻力发散特性较好的传统超临界翼型与上述跨音
速自然层流翼型进行对比#标号分别为&超临界翼型!G7@\%7&B;Y"’G7@\%7&B];Y#’G7@\%7&B];Y)#如图,所示
.
图&!G7@\%7&B];Y)反设计压力分布与
目标压力分布对比
’()*&!P01L@1S S5@1>7S:@7Z5:7%2 <%I L9@7S%2Z1:[112:017281@S1B>1S732 @1S5&:92>:01G7@\%7&B];Y)%
21图+!初始翼型与两个反设计翼型比较
’()*+!P01<%I L9@7S%2Z1:[112:01
727:79&97@\%7&92>:01:[%7281@S1B>1S732
%21
S
图-!超临界翼型与反设计翼型比较
’()*-!P01<%I L9@7S%2Z1:[112:01
S5L1@<@7:7<9&97@\%7&92>:01
7281@S1B>1S732%21S
!!数值模拟,个翼型在雷诺数分别为C**万’A**万’+**万的气动性能.图E i A给出了各个翼型之间压力分布和阻力发散特性以及层流区随马赫数’雷诺数大小的变化对比.可以看出在雷诺数低于A**万时#随着马赫数的增大#激波位置后移#同时上表面转捩位置后移.G7@\%7&B;Y翼型从压力分布上看#翼型上表面没有明显的顺压梯度#但依然可以依靠上表面的弱激波维持E*W左右的层流区域#并且层流区域是随着激波位置不断向翼型后缘移动.对于层流理念设计的G7@\%7&B];Y#和G7@\%7&B];Y)来说#翼型上表面层流区域的变化基本与激波位置变化相一致$从下表面的转捩位置与压力分布的比较来看#翼型下表面的层流区域的维持是需要压力分布顺压区域的#顺压区域越长#层流区越长.从不同雷诺数下的阻力发散特性比较可以看出#雷诺数在A**万以下#层流理念设计的跨音速层流翼型在层流区大小’阻力系数大小方面有着无可替代的优势#但阻力发散特性要低于传统超临界翼型#这是因为自然层流翼型由于顺压梯度的存在#随着马赫数的增大#激波不断增强后移#为激波后翼型压力恢复带来过大的逆压梯度#因此带来较大的激波阻力!波阻"#过大的波阻削弱了层流摩阻降低的效益#总阻力的优化在高马赫数时不是特别明显.
从A**万雷诺数及+**万雷诺数翼型气动特性的比较可以看出#G7@\%7&K];Y#’G7@\%7&K];Y)翼型上表面压力分布在)*W i E*W弦长区域出现了微小的波动!逆压梯度"#这个微小的波动导致了转捩的发生$下表
D+
郑州大学学报!理学版"第"#
卷
图.!C**万雷诺数气动特性及层流区比较
’()*.!P 01<%I L9@7S %2%\91@%>M 29I 7<L1@\%@I 92<192>&9I 729@@137%29:C I 7&&7%2_1M 2%&>S 25I Z1@
面的层流区域的变化基本与顺压区吻合#层流理念设计的翼型#可以维持层流区域到C*W 弦长的位置#微小的压力分布波动并未影响层流区域的维持.从A**万雷诺数之下的翼型压力分布也可以观察到这个小波动#但在A**万雷诺数之下#这个微小的逆压梯度并未造成转捩的发生#转捩的发生与激波有很大的关系#但超过A**万雷诺数后#翼型上表面的转捩对顺压梯度要求变得更加苛刻#需要对压力分布进行更细致的设计.从A**万及+**万雷诺数的阻力发散特性可以看出#在低马赫数时#层流设计的翼型依然能带来气动效率的
提高#但在马赫数大于*’
D"时#由于翼型前部顺压梯度带来过大压力恢复梯度#陡增了波阻#降低了翼型的阻力发散特性
.
图/!A**万雷诺数压力分布比较
’()*/!P 01<%I L9@7S %2%\L@1S S 5@1>7S :@7Z5:7%29:A I 7&&7%2_1M 2%&>S 25I Z1
@
图0!A**万雷诺数气动特性及层流区比较
’()*0!P 01<%I L9@7S %2%\91@%>M 29I 7<L1@\%@I 92<192>&9I 729@@137%29:A I 7&&7%2_1M 2%&>S 25I Z1@
综上#跨音速自然层流设计的翼型典型的问题是在中高雷诺数下阻力发散特性明显下降#对顺压区域要求更加严苛$与传统的超临界翼型对比表明#在中低等雷诺数下顺压设计的层流翼型阻力特性明显#阻力发散特性不是很差#转捩区域与激波位置有很大关系#微小的逆压梯度并未造成转捩的发生$在中高等雷诺数
A +
!第#期李!静#等$跨音速自然层流翼型反设计
研究图1!+**万雷诺数压力分布比较
’()*1!P 01<%I L9@7S %2%\L@1S S 5@1>7S :@7Z5:7%29:+I 7&&7%2_1M 2%&>S 25I Z1
@
图2!+**万雷诺数气动特性及层流区比较
’()*2!P 01<%I L9@7S %2%\91@%>M 29I 7<L1@\%@I 92<192>&9I 729@@137%29:+I 7&&7%2_1M 2%&>S 25I Z1@
下#抑制转捩的发生需要更严苛的顺压梯度#阻力发散特性变得更差#在这种情况下#可以考虑采用严苛的弱顺压梯度’短顺压长度和及早的弱激波恢复等设计理念进行设计.
.%结论
本文基于差量采样解V
9LLM U e d 的翼型反设计方法#对跨声速自然层流翼型进行了反设计#通过分析设计结果可得到以下结论&
#"层流理念设计的跨音速翼型在A**万雷诺数之下#层流区域的维持与激波位置有很大的关系#对顺
压区域的要求不太明显#小的逆压梯度并未造成翼型上下表面转捩的发生$马赫数小于*’
D"时#气动效率明显优于传统的超临界翼型$马赫数大于等于*’
D"时#由于波阻的陡增#降低了层流带来的气动效率的增益#因此若考虑阻力发散特性#需要控制翼型上表面的层流区域#牺牲设计点的气动效率#以带来相对较好的阻力发散特性.
)"层流理念设计的跨音速翼型在雷诺数大于A**万时#层流区域的维持对顺压梯度要求变得严苛#微小的逆压梯度会造成翼型上表面转捩的发生#需要对逆压梯度进行更细致的设计#阻力特性随着马赫数的增大变得更差#在这种情况下#需要考虑采用严苛的弱顺压梯度’短顺压长度和及早的弱激波恢复等设计理念进行设计.
参考文献!
(#)!;P 6_d /GU .G 291@%>M 29I 7<>1S 732I 1:0%>\%@S 5L1@S %27<29:5@9&&9I 729@\&%[97@<@9\:(d ).;:92\%@>&;:92\%@>62781@S 7:M #
)**E.
())!乔志德.自然层流超临界翼型的设计研究(-).流体力学实验与测量##++A ##)!E "&),X ,*.
(,)!/T G (Vf #F G (Vg #Y T =(T #1:9&.;5L1@<@7:7<9&29:5@9&&9I 729@\&%[97@\%7&%L:7I 749:7%2\%@@137%29&97@<@9\:[723>1S 732
+
+
**#
郑州大学学报!理学版"第"#卷(-).G@1%S L9<1S<712<192>:1<02%&%3M#)*#"#E,&#")X#CE.
(E)!黄江涛.飞行器气动外形优化设计方法及应用研究(d).西安&西北工业大学#)*##.
(")!?^;;=]a^P/F#H=_b e a^P/;d#U=_/H e_(=.d1S73292>>181&%LI12:%\:01T%2>9-1:(-).-%5@29&%\97@<@9\:#)**"#
E)!,"&D""X DCE.
(C)!李权#段卓毅#张彦军.民用飞机自然层流机翼研究进展(-).航空工程进展#)*#,#E!E"&,++X E*C.
(D)!朱自强#陈迎春#王晓璐#等.现代飞机的空气动力设计(?).北京&国防工业出版社#)*##.
(A)!]=V=;]=fU#G]e(;e-.^281S:739:7%2%\2%2B&7219@L@%R1<:7%2\%@U e dZ9S1>@1>5<1>%@>1@I%>1&S\%@91@%>M29I7<S(Y)"
,+:0G1@%S L9<1;<712<1S?11:72392>=‘07Z7:._12%#($#)**#&+)C.
(+)!H6^P T G(T P#d G?e d G_G(?#a^]]Y e gb.G1@%>M29I7<>9:9@1<%2S:@5<:7%292>7281@S1>1S7325S723L@%L1@%@:0%3%29& >1<%I L%S7:7%2(-).G^G G-%5@29&#)**E#E)!A"&#"*"X#"#C.
(#*)刘浩#徐敏#叶茂.基于特征正交分解的跨声速流场重构和翼型反设计方法研究(-).空气动力学学报#)*#)#,*!E"& ",+X"E".
(##)段焰辉#蔡晋生#刘秋洪.基于代理模型方法的翼型优化设计(-).航空学报#)*###,)!E"&C#D X C)D.
(#))]G(V P_f_H.G<%@@1&9:7%2Z9S1>:@92S7:7%2I%>1&5S723&%<9&89@79Z&1S%L9@:^^&:1S:<9S1S92>72>5S:@79&9LL&7<9:7%2S(-).
-%5@29&%\:5@Z%I9<0721@M#)**C##)A!,"&E),X E,E.
(#,)?=(P=_F_#]G(V P_f_H#]^b b^;_.P@92S7:7%2I%>1&&723\%@3121@9&L5@L%S1Y F d<%>1S(-).F&%[:5@Z5&12<192> <%I Z5S:7%2#)**C#DD!#Q E"&)DD X,*,.
(#E)b6]F G(H?.62781@S9&L9@9I1:@7<31%I1:@M@1L@1S12:9:7%2I1:0%>(-).-%5@29&%\97@<@9\:#)**A#E"!#"&#E)X#"A.
#:Q<B><!<>():9H;N<3B4:>9:(E R4;D B4?
V4@(:4B’?9K7(B H9(?
]^-723##]^6f92L723##V G e/01230%23)#]^G(VP792S057#
!#’)234454!,<23:7;2B:7=>7?;7<<8;7?)2;<72<#O3<7?J34@&7;A<8B;"C#O3<7?J34@E"***##$3;7:$ )’*:";47:5D<C Q:P48:"48C4!H<84=C7:9;2R<B;?7:7=+<B<:823#*48"3E<B"<87
(45C"<237;2:5&7;A<8B;"C#M;N:7D#**D)#$3;7:"
75>;B4E;&P017281@S1>1S732I1:0%>[9SS:5>71>Z9S1>%239LLML@%L1@%@:0%3%29&>1<%I L%S7:7%2
!U e d".G<<%@>723:%S5@\9<1L@1S S5@1>7S:@7Z5:7%2#:0139LLM U e d@1>5<1>%@>1@I%>1&[9S5S1>:%@1B <%2S:@5<:I7S S723>9:9%\9297@\%7&.^2:01727:79&S29LS0%:S9I L&723L@%<1S S#!Y;P[9S5S1>:%31::01 S29LS0%:#:012:01Z9S7<97@\%7&[9S@1L&9<1>ZM:017281@S1>1S7321>#92>:0121[S29LS0%:>1S7321>ZM U e d[9S9>>1>:%S29LS0%:.F%@&9I729@\&%[97@\%7&7281@S1>1S732#L@1S S5@1>7S:@7Z5:7%2S29LS0%:[1@1 3%:ZM%b_1
:@92S7:7%2I%>1&.P[%>7\\1@12:L@1S S5@1>7S:@7Z5:7%2S%\:01:@92S%27<(]F97@\%7&[1@1L@%B
&
L%S1>.P012:[%:@92S%27<(]F97@\%7&S<%5&>Z1%Z:9721>ZM:017281@S1>1S732I1:0%>.P019Z7&7:M\%@ L%S:L%2723:01:@92S7:7%2%\>7\\1@12:L@1S S5@1>7S:@7Z5:7%2S52>1@>7\\1@12:_1M2%&>S25I Z1@92>?9<0 25I Z1@S[1@1S:5>71>72>1:97&.
J<G K9B C>&Y;P$(]F97@\%7&$:@92S7:7%2I%>1&$7281@S1>1S732
!责任编辑&方惠敏"。