基于重叠场源法的非定常气动力计算研究-张辉.pdf

《基于重叠场源法的非定常气动力计算研究-张辉.pdf》由会员分享，可在线阅读，更多相关《基于重叠场源法的非定常气动力计算研究-张辉.pdf（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、2016年6月第34卷第3期西北工业大学学报Joumal of Northwestem Polytechnical UniverSityJune 2016V0134 No3基于重叠场源法的非定常气动力计算研究张辉，李杰，韩杰(西北工业大学航空学院，陕西西安710072)摘要：面向气动弹性工程应用，将基于雷诺平均NS方程的定常流动数值解引入跨声速小扰动方程，在小扰动条件下非定常激波效应由定常激波效应确定，从而发展了一种可用于跨声速流动的非定常气动力计算方法。数值算法利用块三对角近似技术大大提高了算法的计算效率、节省了计算所需内存空间，并采用了重叠场源策略为复杂构型的非定常气动力计算提供了有利保障

2、。文中以F5机翼为计算算例，研究了场源模型参数对非定常气动力计算结果的影响，并验证了块三对角近似算法的可靠性；以cRM wBH翼身尾构型为算例，检验了针对复杂构型重叠场源策略的可行性。关键词：非定常流动；跨声速流动；马赫数；气动力；压力分布；计算效率中图分类号：V215 文献标志码：A 文章编号：10002758(2016)03044306非定常气动力的计算对气动弹性分析、动导数计算及机动、突风载荷计算等都具有特别重要的意义。J。长期以来，广泛应用于气动弹性分析的气动力计算方法是偶极子格网法”。5 o(doublet latticemethod，DLM)，多个商用软件均采用了DLM进行气动弹性

3、分析，如MsCNASTRAN的气动弹性模块。DLM方法能够提供亚声速及超声速范围内较为精确的非定常气动力，其最大的优点还在于能够生成气动力影响系数(AIc)矩阵。AIc矩阵不依赖于飞机结构参数，因此，在结构设计循环中只需计算1次而重复使用。但是DLM是基于线化势流理论的方法，无法解决非线性较强的流场，不适用于跨声速非定常气动力计算。计算流体力学(CFD)方法通过求解Euler或Navie卜Stokes方程能够给出精确的跨声速流场解，不过基于cFD的流场求解不能生成AIc矩阵，因而不能利用目前已经发展较完备的基于AIC矩阵的气动弹性分析方法。而目前用于气动弹性分析的cFDcSD耦合计算方法一1针

4、对复杂飞机构型的适用性还有待发展，另外耦合算法需要巨大的计算资源和时间耗费，这样均不利于工程应用。因此，发展一种能够生成AIc矩阵的高效非定常跨声速气动力计算方法对跨声速范围内飞行器的气动弹性特性研究意义重大。非定常场源方法即是一种能够生成AIC矩阵的方法，并且通过引入CFD定常解考虑了跨声速激波效应，从而适用于非定常跨声速气动力的计算。早在20世纪80年代，国外一些研究者就已经开始了对非定常场源法的研究，Larry LErickson和Shawn MStrande【81给出了利用场源法将面元方法推广以求解跨声速流动的理论基础；MHL Hounjet一1研究了基于场源法和有限差分的混合方法求解

5、非定常跨声速流动问题，该方法综合了有限差分和场源积分方法的优点，从而大幅地减少了计算时间，使得其可用于常规的颤振分析；Lutz Gebhardt和DmitriFokin等刮研究了用于跨声速飞行器气动设计的场源方法，改进了自适应场源网格生成方法使得场源法的实际应用更加现实；chen和Gao等o基于重叠场源策略通过求解关于时间线化的跨声速小扰动方程的积分方程生成非定常跨声速气动力影响系数矩阵，从而实现复杂构型的跨声速气弹分析，并利用块三对角近似技术求解大型的体单元影响系数矩阵，大大地提高了计算效率。然而，国内对基于场源法的跨声速非定常气动力计算的研究目前尚无公开发表的文献。收稿日期：2015-10

6、22基金项目：国家自然科学基金(11172240)、国家重点基础研究发展计划(2015cB755800)与航空科学基金(2014zA53002)资助作者简介：张辉(1986一)，西北工业大学博士研究生，主要从事理论与计算流体力学研究。万方数据西北工业大学学报 第34卷本文深入研究了非定常场源法的算法实现，将基于雷诺平均NS方程的定常流动数值解引入跨声速小扰动方程，在小扰动条件下非定常激波效应由定常激波效应确定，从而发展了一种可用于跨声速流动的非定常气动力计算方法。非定常场源法的求解方程旷“品劫卜卜n未(G岫出(6)得到妒。=li品去盯盯。G嶝：dJs+去够盯dy(7)式中，戈，表示激波位置，占

7、表示激波面无限小厚度。由(5)式和(7)式可得关于结构振荡幅值线化非线性跨声速小扰动方程得到时间线化的跨声速小扰动速势方程如下：一2等妒=扣29“+妒+纩2丁妒2孟盯”2式中，I|为减缩频率，卢=1一M：，盯。=脚。币。，K=(y+1)M：佃2，7为比热比。方程右侧项咖。为沿来流方向的定常扰动速度分量，由定常流动数值解确定，包含了跨声速流动中的定常非线性跨声速激波效应，假定方程(1)右侧为体源，则方程(1)在积分方程的积分域内某点(z。，)，。，z。)上的积分解包括3部分：妒(戈o，yo，zo)=9，+9r+妒。hock (2)9。表示表面面元强度对速度势的影响，可写为旷捌争埘卜 式中，C。为

8、升力面上非定常压力差分布，盯为体单元上非定常源分布，K、G分别为加速度势核函数和非定常源核函数，如下露：一广ei愉当丁G=枷渺)(网R式中，尺=(zoz)2+()，oy)2+(zoz)2。妒P为体源对速度势的影响，妒。为激波面对速度势的影响，可分别写为旷一劫熹(Gdy (4)。一寺峨GdS (5)k一石掣舸”洲3 )式中，以表示穿过激波面体源强度的跳跃。当不存在跨声速激波时，盯，=O，9。自动消失。另外，当激波出现时，也可通过对吼完成如下分部积分运算以消掉妒。弧=劫邓，dy (8)k+妒”2石删矿rG，洲 (8)定常流动数值解的引入本文用于求解定常流动数值解的控制方程如下去巧Qdy+甄F。露d

9、s=。 c9，式中，Q=(p，pM，p秽，p训，pe)1为守恒向量；p、e、(“，”，似)分别为密度、单位质量气体的总能量和直角坐标系下的速度分量；为控制体体积；s为控制体表面积；F为通过表面s的黏性通量和无黏通量的和；以为边界外法向单位矢量。黏性项采用中心差分格式离散，无黏项采用Roe三阶迎风偏置通量差分裂方法离散，采用隐式近似因子分解(AF)方法进行时间推进，选用SA湍流模型，通过多重网格技术来加速收敛。采用场源法求解非定常气动力时，需要建立场源模型，场源模型包括面单元和体单元。首先建立面元气动模型，通常将飞机部件分为翼面类和机身类部件处理，翼面类部件由位于翼面类部件均平面的非定常涡面来模

10、拟，机身类部件表面离散为体表面单元，每个体表面单元上布置非定常源，以模拟由于体的体积效应产生的气动力分布；其次围绕翼面类部件(升力面)或机身类部件的面元气动模型定义体块，然后将体块分割为若干体单元。场源模型建立后，即可将上述定常流场数值结果(包括当地马赫数及扰动速度分量)插值到每个体单元上。建立合适的场源模型直接决定着插值精度，进而影响最终非定常气动力计算的精度。在建立场源模型时，确保场源模型中的物面边界尽可能与CFD计算模型的物面边界重合是建立合适的场源模型的必要条件。此外，体块高度和体块分割层数是场源模型的2个重要参数，这2个参数的选取对计算结果有万方数据第3期 张辉，等：基于重叠场源法的

11、非定常气动力计算研究重要影响，选择的依据是确保场源模型包含非线性流动区域。算例验证部分详细研究了体块高度和体块分割层数对计算结果的影响，并归纳总结了体块高度和体块分割层数选取的一些原则。3 AIC矩阵的生成及算法策略31 AIC矩阵的生成位于面单元控制点处的法向扰动速度为fCPlnV妒=形=Aa妒+曰盯。【磊J(10)体单元控制点处的速度势为rCPl=ca9+Dh (11)【磊J式中，矩阵A和C分别为面奇点对面单元和体单元的影响系数矩阵，矩阵曰和D分别为体源对面单元和体单元的影响系数矩阵。引入如下有限差分算子r盯。=r妒 (12)将方程(12)代入(10)式和(11)式，得到坼cp) (13)

12、式中，五=A+BrE。1C，E=，一Dr，矩阵五求逆即得到非定常气动力影响系数矩阵。32块三对角近似技术上节中矩阵E为满系数复数矩阵，其阶数与体单元数相同。对于简单的三维问题需要较多体单元数，并且体单元数随构型复杂程度的增加而急剧增加。这种大型矩阵求逆对常规的气弹及载荷分析是不实用的。因此，有必要采用块三对角近似技术解决大型满系数矩阵层求逆的问题。块三对角近似技术的思路11是：首先体单元被分为许多子块，同一个子块内的体单元组合在一起。这样矩阵E可写为E=E。+E。 (14)式中，E。为块三角矩阵，其三角块包含自身块和相邻子块的影响系数；E。包含三角块处的零元素和非相邻块的影响系数。下一步比较矩

13、阵E。和矩阵E。内系数的量级，可以看出E，内所有系数均是小量，原因是方程(4)的积分方程的积分函数包含1尺函数，当点(戈。，y。，)远离非相邻子块时，l尺快速衰减。因此，矩阵E的逆可由下式近似E一E。一一层。1E。E。1(15)因为E。为块三角矩阵，所以E。叫可利用块三角矩阵求解技术有效计算。33针对复杂构型的重叠场源策略针对复杂飞机构型，本文采用重叠场源策略以减小网格生成的难度，减少计算所需内存并节省计算时间。首先针对复杂构型的每一个部件独立生成体单元网格，然后基于以下原则构建这种重叠体单元模型的影响系数矩阵：不同体块内的体单元互不影响；同一个体块内的体单元仅影响与其相关联的面元；所有面元影

14、响所有体单元。4算例验证41 F5机翼本节以F5机翼1引为算例，建立不同参数下的场源模型，取计算马赫数为090，减缩频率为0275，进行计算分析并与实验值对比以研究场源模型不同参数对非定常气动力计算精度的影响。主要研究了场源模型的体块高度日、体块分割层数f、体块前(后)缘向前(后)延伸量ExT及分割数M对非定常气动力计算精度的影响。以机翼根弦长C为参考，分别取体块高度025C、050C、100c和150C，保持其他参数不变，其中体块层数为7。计算所得展向515位置处非定常压力系数结果与实验值对比如图1所示，当体块高度较小，如H=025C和050C时，计算结果明显小于实验值，计算对激波强度的预测

15、明显偏小；随着体块高度的增加，在H=100C时，计算结果与实验值比较接近，计算精度较好；随着体块高度的继续增加，在H=150c时，计算结果则高于实验值，显然这不是一个理想的结果，因为这与之前的理论分析相悖，按照理论分析当日=150C时其场源模型比其他场源模型大，所以合理的结果应该是此时的计算结果应更加接近实验值，但图2结果并非如此。进一步，保持体块高度为H=150c，选取体块层数分别为5、7、9、1l完成计算，同样给出展向515位置处非定常压力系数的计算结果与实验值对比如图2所示，对比结果表明，随着体万方数据446 西北工业大学学报 第34卷块层数的增加，计算结果不断逼近实验值，并且，=11对

16、应的计算结果略好于图1中体块高度日=100C的计算结果。因此，可以得到如下结论：增大场源模型以确保场源模型包含非线性流动区域可以提高非定常气动力的计算精度；若保持其他参数不变，只是增大体块高度会使非定常气动力的计算精度变差，原因是随着体块高度的不断增加，使得体单元分布越来越稀疏，过于稀疏的体单元分布降低了其计算精度；保持其他参数不变，增加体块层数，计算结果趋于收敛；另外，图3给出了一个状态所需的CPU计算时间随体单元个数的变化曲线，可见CPU计算时间随着体单元数的增加呈指数式增加。因此，为了兼顾精度与效率，对于该算例，选取体块高度为1OC，体块层数为7较合适。图l 体块高度的影响 图2体块分割

17、层数的影响 图3计算时间随体单元数的变化取体块高度为10C，体块层数为7，分别选取体块前(后)缘向前(后)延伸量ExT为000C、010C、025C、050C的场源模型进行计算，计算结果与实验值对比如图4所示，可见随着延伸量ExT的增加，计算结果向实验值靠近，但变化量很小。图5给出了延伸段不同分割数的场源模型计算结果对比，可以看到分割数分别为2、4、6时，计算结果基本重合，表明该处分割数对计算结果几乎没有影响。针对F5机翼，分别建立单块场源模型和多块场源模型，对于多块场源模型，利用块三对角近似技术进行求解，对比2种场源模型计算结果，验证块三对角近似的可靠性。图6给出了机翼展向515位置处压力分

18、布对比结果，可见三对角近似求解结果与单独体块求解结果基本重合，这表明了三对角近似求解是可靠的。此外，与单独体块求解相比，块三对角近似求解节省了大约40的CPu计算时间。图4场源模型前后延伸量的影响 图5延伸段分割数的影响 图6块三对角近似技术的验证，=027542 CIU订WBH翼身尾构型采用重叠场源策略建立CRM wBH飞机构型131场源模型，分别针对机翼、平尾、机身独立生成体块，每个体块被分割为若干子块生成体单元网格以应用块三对角近似技术求解。根据上节研究结论，机翼和平尾对应体块高度取各自的根弦长，体块层数为7，机身对应体块高度取为机身中段半径长，体块层数为8，建立场源模型如图7所示。万方

19、数据第3期 张辉，等：基于重叠场源法的非定常气动力计算研究图7 cRM wBH构型场源模型取计算马赫数085，减缩频率015进行计算。图8给出了CRM wBH构型俯仰振荡模态时机翼展向502处翼面非定常压力分布对比结果。由于没有实验结果，本文以zAERO软件4 o的跨声速方法计算结果作为参考验证本文计算方法的可靠性，对比可见本文计算结果与zAERO计算结果基本一致，仅在对激波强度的预测方面稍有差异，这说明了针对复杂飞机构型本文采用的重叠场源策略是可靠的。同时图中给出了偶极子格网法的计算结果，显然，由于偶极子格网法不能考虑跨声速激波效应的影响，而对于当前计算状态下机翼翼面展向502位置处存在较强

20、激波的状况，偶极子格网法计算结果必然与考虑了跨声速激波效应的计算方法所得结果存在较大差异。换言之，如果翼面不存在激波，也不存在流动分离，那么偶极子格网法计算结果就应该与本文气动力方法计算结果接近。平尾翼面压力分布对比结果有利地证实了这一点。图9给出了俯仰振荡模态时平尾展向50处翼面非定常压力分布对比结果。图9中本文气动力方法、偶极子格网法及zAERO跨声速方法计算结果可以看到，3种计算方法所得结果基本一致。参考文献： zAERO跨声速方法计算结果一偶板子格网法计算结果i本文气动力方法计算结果、0 O4 08xl c图8机翼展向502处非定常压力系数图9平尾展向50处非定常压力系数5 结 论本文

21、通过引入CFD定常解考虑跨声速激波效应，面向气动弹性分析发展了一种跨声速非定常气动力计算方法。文中以F5机翼为算例研究了场源模型参数对非定常气动力计算结果的影响，得到了选择合适场源模型的一般规律，并验证了块三对角近似算法的可靠性；以CRM wBH翼身尾构型为计算算例，单块和多块场源模型计算结果的对比表明针对复杂构型重叠场源策略是可行性的。算例计算结果与实验及参考文献结果的对比分析表明本文发展的方法适用于跨声速非定常气动力的计算。1 Rafael N，Joao Lcompason of Viscous and InViscid unsteady Aerodynamic Loads for Aer

22、oeLastic An出ysisRAIAA一201507672 Hyungro Lee，Beom-s00 Kim，et al，computational study of the stab订ity DeriVatiVes for the Standard Dynamic ModelRAIAA一20 l 326583Leonard M，Ilhan TTime simulations of the Response of Maneuvering Flexible AircraftJJoumal of Guidance，and Dynamics，2004，27(5)：814828万方数据456789

23、10121314Albano ERodden W PA DoubletLattice Method for Calculating Lift Distributions on 0scillating Surfaces in Subsonic FlowsJAIAA Joumal，1969，7(2)：279285管德非定常空气动力学计算J北京：北京航空航天大学出版社，1991Guan DeCalculation on unsteady Aerodyn帅icsMBeijing，Beijing university of Aeronautics and Astmnautics Press，1991(i

24、n Chinese)徐敏，安效民，陈士橹一种cFDcSD耦合计算方法J航空学报，2006，27(1)：3337xu Min，An xiaomin，chen shilucFDcsD couping Numerical computational MethodologyJActa Aemnautica et Astmnautica Sinica，2006，27(1)：3337(in Chinese)Giulio Romanelli，Michele Castellani，Paolo Mantegazza，et al，Coupled CSDCFD Non-Linear Aemelastic Trim

25、of FreeFlyingF1exible AircrafIRAIAA一20121562Erickson IJany L，stmnde Shawn MA Theoretical Basis for Extending surface Paneling Methods to Transonic FlowJAIAAJoumal，1985，23(12)：18601867Houniet M H L_A Field PaneLFinite Di任色rence Method for Potential Unsteady T眦sonic FlowJAIAA Joumal，1985，23(4)：537545G

26、ebhardt Lutz，Fokin Dmit而，Lutz Thorsten，el a1An ImplicitExplicit DirichletBased Field Panel Method for Transonic Aircraft DesignRAIAA一20023 145Chen P C，Gao X W，Tang L0verset FieldPanel Method for Unsteady Transonic Aemdyn枷ic Irlnuence Coemcient MatrixGeneralionJAIAA Joumal，2004，42(9)：1775一1787Malone

27、J B，sankar N L，et a1Unsteady Aemdynamic Modeling of a Fighter wing in Transonic FlowRAIAA19841566John c，Edward Nsummary of the Fourth AIAA CFD Dmg Prediction workshopRAIAA一2010-4547zAER0 Theoretical ManualMz0NA Technology Inc Scottsdale，2008Calculation of Unsteady Aerodynamics UsingOverset FieldPane

28、l MethodZhang Hui，Li Jie，Han Jie(College of Aeronautics，Northwestem Pol”echnical UniVersity，Xian 710072，China)Abstract：The paper present a method of calculating unsteady aerodynamics in transonic now，in which steady nowsolution based on RANS solver code are imported into transonic small disturbance

29、equations， thus， the unsteadyshock location is dominated by the steady shock location in the smaU amplitude senseBlock-tridiagonal approximation is used to improve the computational emciency and save the memory，and the overset 6eld panel scheme is de-veloped for complex connguration to calculate uns

30、teady aemdynamicsTaking F5 wing as an example，the efI-ect ofparameters of 6eldpanel model on the calculated results of unsteady aerodynamics is inVestigated， and thereliability of block-tridiagonal approximation is validatedIn addition，the validation of oVerset field-panel scheme isperfbImed by comparing the results between the one block modeling and the oVerset modelingKeywords：unsteady flow；transonic now；Mach number；aerodynamics；pressure distribution；computational em-ciency万方数据