基于Euler方程结合附面层修正方法的应用研究

采用欧拉加附面层修正的方法在结构网格上对民用飞机翼身组合体流场进行了数值模拟研究.采用一种适合翼身组合体模型的的网格生成技术来生成计算用的结构网格,该技术是先采用一种基于Higenstock源项修正的椭圆型偏微分方程生成二维网格,然后沿机身轴向合理的光滑性调整成三维网格.在求解Euler方程时采用中心离散格式和显式时间推进,并利用附面层修正的方法来考虑粘性的影响.为了验证此方法的正确性,本文对DLR-F4翼身组合体进行了数值模拟,结果表明此方法和试验结果能够很好的吻合.     

增升状态下的全机位流数值模拟

研究了势流理论下的全机增升装置复杂外形气动力数值计算方法,该方法以低阶面元法为基础,具有计算量小、速度快的优点.复杂流动的粘性效应主要通过附面层粘性修正技术加以体现.采用该方法对某增升装置翼身组合体外形压力分布特性进行计算,并和风洞实验数据相比较,验证了方法的正确性和实用性.在此基础上对某型飞机真实复杂外形的压力分布特性进行了计算,为复杂外形飞行器增升装置的初步设计、方案选型与评估提供了非常理想的工程技术手段.     

桨叶变形对旋翼气动性能的影响分析

针对某型机桨叶出现大面积变形情况,本文通过计算评估,分析了其在旋翼气动性能方面的影响.首先采用三坐标测量仪测量了变形桨叶的外形数据,并与理论外形进行了对比,然后用Euler方程加附面层修正方法,计算了二维翼型气动性能的变化情况,进而采用二阶升力线方法分析桨叶变形对旋翼性能的影响,得到了一些有价值的结论.     

数值模拟二维喷管激波/湍流附面层干扰流动

采用可压缩性修正两方程湍流模型,数值模拟了3种不同波前马赫数的跨声速二维喷管内激波/湍流附面层干扰流动,对流场中时均参数和脉动参数的计算结果与实验值进行了比较。结果表明可压缩性修正的两方程湍流模型准确地模拟了正激波/湍流附面层干扰流动的时均参数和脉动参数,无分离和有分离的激波/湍流附面层干扰流动的基本规律。      

NLR-7301高升力翼型绕流评估计算

探讨了在二维高升力多段翼型绕流计算中附面层网格采用各向异性三角形网格的技术。通过对经典的NLR-7301两段翼型绕流的计算,从数值计算角度评估了自研混合网格雷诺平均N-S方程求解软件对该类型附面层网格的解算能力。通过与确认试验数据对比,在附面层采用各向异性三角形网格进行二维高升力多段翼型绕流计算,其结果具有较高的可信度。在此基础上,采用同样策略,评估了多段翼型绕流计算中翼型后缘的几何处理、湍流模型选取等关键CFD要素对计算结果的影响。     

亚超音速全机纵向气动系数数值计算方法的改善与发展

本报告是在亚超音速位流解基础上改进发展的数值计算方法。用面分布有限基本解计算作用在物面上的压力,用数值积分方法计算作用在物体上的力与力矩。本文采用联合流场、有效剖面吸力概念以及顺流向片条作可压缩紊流附面层的积分解,来改进阻力特性的数值计算。 数值方法已编制成计算程序。 本文提供了用这个计算程序计算的若干算例,并与实验结果作了对比。比较表明,计算结果与实验值有良好的一致性。     

涡轮叶栅叶型损失的数值模拟方法

给出了一个计算亚、跨音涡轮叶栅叶型损失的数值计算方法。主流采用时间推进有限体积法求解积分型欧拉方程,并采用了局部网格修正技术;附面层采用全隐格式求解有限差分方程;在叶栅出口与远后方均匀流之间进行了叶片尾迹与主流的掺混损失计算。算例表明本文的数值方法可准确地预测涡轮叶栅的叶型损失。     

基于CFD结果的轴流压气机附面层堵塞计算方法

为研究轴流压气机附面层堵塞计算方法,提出一种基于计算流体力学结果自动提取附面层厚度的多参数组合方法,该方法整合了多种附面层判断准则,可有效避免单参数判断准则的局限性。通过多个算例,包括二维平板、多级轴流压气机端壁、进口导流叶片以及带有分离区的静子叶片等,对该多参数判断准则进行了检验。结果表明,此附面层提取方法具有很好的有效性和普适性。该研究建立了多级轴流压气机中端壁堵塞从高维向低维模型的反馈,为压气机设计体系改进提供了更多可能。     

RANS方程和附面层方程耦合求解转捩位置的方法

研究了基于线性稳定性分析的RANS（Reynolds-Averaged Navier-Stokes）方程和附面层方程耦合求解转捩位置的方法。RANS方程计算得附面层外边界的速度分布作为附面层方程计算的输入,求解附面层方程得到稳定性分析需要的附面层参数。由于附面层的解在分离点处的奇异性,用交互式附面层方程求解方法。用基于线性稳定性分析的e^N方法计算转捩位置,并考虑了Tollmien-Schliching波和层流附面层分离造成的转捩。迭代求解RANS方程和附面层方程,直到转捩位置收敛。RANS求解中使用转捩过渡区模型,避免了点转捩模型引入的数值扰动。通过对NA-CA0012翼型的计算并和实验结果及Xfoil计算的结果进行比较,吻合较好,有较好的工程实用价值。     

复杂外形网格生成技术

复杂外形网格生成是进行流场数值模拟的一个关键技术及前提。本文论述了采用保角变换方法，求解偏微分方程方法以及它们在代数方法结合生成二维及三维具有多块网格结构的贴体且与边界正交网格的技术。     

A Structured Mesh Euler and Interactive Boundary Layer Method for Wing/Body Configurations

To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.     

超临界层流机翼边界层及气动特性分析

高空长航时无人机设计巡航状态的雷诺数较小,黏性边界层对气动特性的影响较大。详细分析了雷诺数对机翼边界层和气动力的影响,用数值方法对超临界层流机翼三维层流-转捩-湍流混合边界层特性进行了研究,分析比较了高空小雷诺数和中空大雷诺数情况下机翼三维边界层的特性,尤其是边界层转捩点位置、表面摩阻和气动特性的雷诺数效应。研究表明雷诺数对于高空无人机机翼边界层厚度、摩擦阻力和升阻比影响较大;对层流机翼的转捩点位置和升力系数影响较小;自然层流机翼技术可以应用于高空无人机设计。     

旋成面叶栅气动反命题的新解法(Ⅱ)

本文利用映象平面ξ-η的特点, 将旋成面叶栅气动正命题的微分-积分解法与叶型修正的“喷气模型”结合, 建立了一种适用于旋成面叶栅气动反命题的新解法。叶片表面给定的气动参数分布作为反命题的边界条件, 在求解过程中直接得到满足。通过对一例实际透平叶型的计算验证了本文解法的可靠性和实用性。      

高超声速进气道壁面开缝对边界层分离影响研究

以二维高超声速进气道GK01[1]为计算模型,对进气道进行无粘、粘性以及边界层分离控制方法数值研究.结果发现,利用文中设计的壁面开缝措施可以有效减小甚至消除激波/边界层干扰带来的边界层分离现象,降低粘性效应造成的负面影响;相对开缝前,开缝措施在设计点和非设计点均能明显提高进气道的总压恢复系数和动能效率,而对流量捕获系数的影响很小.同时文中给出了开缝措施在非设计点下对进气道性能指标的影响规律.     

基于混合笛卡儿网格方法的可压流动数值模拟

以可压缩黏性流动的数值模拟为研究背景,发展了一套自适应混合笛卡儿网格(AHCG)方法以及基于有限体积方法的雷诺数平均Navier-Stokes(RANS)的数值求解方法.为更好地模拟边界层的黏性流动在近壁面处采用贴体结构网格,剩余计算区域自动生成与之相重叠的笛卡儿网格,并同时发展了基于流场特征的笛卡儿网格自适应技术.结合ADT(alternating digital tree)算法显著减少了网格生成中“挖洞”和“贡献单元”搜索的消耗机时,50万左右的网格数目下,搜索耗时为0.062s,仅为普通遍历方法的1/1847.通过二维圆柱与两段翼型绕流的数值算例显示,定常AHCG方法能够准确地预测物面压力分布与升阻力系数并且具备处理复杂外形的能力;同时通过二维非定常圆柱绕流问题与NACA0015矩形机翼翼尖尾涡的捕捉算例显示,结合了动态自适应网格加密的非定常AHCG方法尤其适用于旋涡主导流动.      

Shock control on a swept wing

Shock and boundary layer control by contour bumps and local boundary layer suction have been investigated experimentally and numerically on a transonic swept wing. Additional 2-D numerical investigations were performed for the airfoil, corresponding to the wing. The investigations were primarily stimulated by the question concerned with the influence of sweep on the bump effectiveness. This influence has been found to be rather small; the drag reduction by the bump is slightly lower for the swept wing than for the airfoil. A location of the bump in the shock region has shown its effectiveness for reducing shock strength and hence wave drag. A position of the bump downstream of the shock wave has been shown to reduce viscous drag and to postpone buffet-onset to higher lift coefficients. Furthermore, the results indicate that boundary layer suction is a powerful device for drag reduction, but the effectiveness decreases with increasing Reynolds number. Higher effectiveness of suction can be attained, when it is coupled with a contour bump. The parameters height and position (relative to the shock) of the bump, optimized in terms of drag, depend on the shock strength; an influence of the boundary layer thickness upstream of the shock on the optimal bump parameters has not been found. A possibility to control an adaptive bump, mounted on an aircraft wing, is to employ the trailing edge pressure.     

粘流与无粘流的相互作用计算

本文总结了粘流/无粘流的各种计算方法和结果。重点在于介绍定常流动中的弱相互作用。首先叙述了弱相互作用的数学模型。给出了不可压流动和跨音速流动中粘流/无粘流相互作用的某些正耦合的计算结果。讨论了在分离区附近边界层正方法失效的原因。然后介绍了边界层反方法和适用于带分离的流动中半反方法耦合的粘流/无粘流的相互作用方法。文中也简单地总结了三维情况的应用和强相互作用。     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

超燃冲压发动机进气道低马赫数起动性能研究

为探讨高超声速进气道在低马赫数普遍存在的起动问题,采用等激波强度法设计了超燃冲压发动机二维混压式前体／进气道,给出了前体／进气道的几何尺寸,对所设计的进气道在设计状态、非设计状态的性能与流场进行了数值模拟,对低马赫数下进气道的起动问题进行了研究。研究表明：设计的进气道附加阻力较小,总压恢复系数较高,在低马赫数下通过附面...     

电磁力控制翼型失速的实验研究

将条状电极和磁极交错布置的电磁场激活板置于弱电解质溶液中,在流体边界层上产生具有明显分布特征的电磁力（Lorentz力）,利用电磁力可以改变流体边界层的结构,控制边界层的流动脱落与分离,可以增加翼型升力,减少其阻力,实现对翼型失速的控制。利用电磁力控制流体边界层的方法属流体主动控制方法之一。本文首先基于电磁场和流体的基本方程,对置于弱电解质中的不同极板宽度的电磁激活板周围的电磁场及产生的Lorentz力进行了数值模拟;其次,通过实验来验证此方法的有效性。将包覆有电磁激活板的翼型置于弱电解质溶液中,利用基于TMS320F2812（DSP芯片）组建的翼型失速实验控制系统来灵活改变翼型的迎角和转速,测量升力和阻力的变化;实验结果表明,正向电磁力能够有效地抑制和延缓翼型失速现象的发生。