机翼边界层计算及其在飞机设计中的应用

本文给出一种计算实际机翼可压缩层流-紊流三维边界层的新方法。方法使用了Cebeci提出的一种非正交坐标系,并采用一种二次变换以克服在某处横流发生反号时计算遇到的困难。对Lockheed C-5A机翼和F-8超临界机翼给出了一系列计算结果,并同其它方法的结果以及试验数据作了比较。     

考虑涡轮传热性能的气动设计耦合计算方法研究

为了在涡轮叶片气动型线设计过程中同时提高气动效率并保证传热性能,提出了一种基于一维管道网络法和三维CFD的耦合计算方法,分别采用管网/三维耦合计算方法和全三维耦合计算方法对MARK-II冷却叶片多个工况进行计算,两种数值计算方法计算结果与实验数据交叉对比,以验证本文计算方法可行性。计算结果表明,两种数值计算得到的叶片型面压力、温度、换热系数和实验值都比较吻合,但管网/三维耦合计算得到的壁面温度相比全三维耦合计算结果整体略微偏低,最大偏差不超过3.89%。基于管网/三维耦合计算方法对某航空发动机涡轮第二级动叶叶片型线优化设计,气动效率提高0.34%,壁面平均温度几乎没有变化。     

轴流压气机动态失速过程三维计算方法研究

基于三维非定常欧拉方程和三维激盘模型发展了一种用于研究轴流压气机动态失速过程的三维计算方法。利用该计算方法研究了某型高压压气机第一级的动态失速过程,并就转子总静压升特性的三种不同径向分布对失速过程的影响进行了分析。计算结果表明该方法能够反映出压气机的三维动态失速过程、失速团的三维空间结构,并能有效地表现出转子沿径向的特性变化对压气机气动稳定性的影响,而且还能够用于判断压气机的危险截面。     

三维不可压缩层流、转捩、湍流边界层的数值计算

本文提出了一种计算三维不可压缩边界层的数值解法,其中包括对层流边界层的微分法,对湍流边界层的积分法以及利用Orr-Sommerfeld方程计算转捩的数值解,当扰动振幅达到e~(13.5)时发生转捩。 按照上述方法,本文计算了ONERA M6后掠角为30°的机翼,计算结果分别和Cebeci和Cousteix的结果作了对比,二者吻合较好。     

轴流压气机动态失速过程三维计算

基于三维非定常欧拉方程和三维激盘模型发展了一种用于研究轴流压气机动态失速过程的三维计算方法.利用该计算方法研究了某型高压压气机第一级的动态失速过程,并就转子总静压升特性的三种不同径向分布对失速过程的影响进行了分析.计算结果表明该方法能够反映出压气机的三维动态失速过程、失速团的三维空间结构,并能有效地表现出转子沿径向的特性变化对压气机气动稳定性的影响,而且还能够用于判断压气机的危险截面.      

扫描有限元法求解二元喷管的三维可压粘性流场

发展了一种二元喷管流场计算的三维Ｎ－Ｓ方程扫描有限元方法,给出了Ｎ－Ｓ方程中二阶导数项新的离散方法,以及更完善、简洁的扫描有限元方程的表达式。开发了相应的计算程序,用所发展的方法对二维平面收扩喷管、两种大宽高比的二元喷管及其圆变方转接段的三维粘性流场进行了计算,计算结果表明本方法具有较宽广的适用范围,计算结果良好     

密布多孔构件的低周疲劳寿命予测和试验研究

一、多孔构件应力场有限元分析　　为了计算多孔的涡轮叶片和燃烧室等构件的应力应变 ,我们编制和开发了壳体弹塑性有限元分析程序。用退化方法从三维单元推导出的曲壳单元。这种单元直接离散连续介质力学的三维方程 ,避免了由一般壳理论带来的复杂。曲面等参元具有独立的转角自由度和平移自由度 ,三维应力应变就退化为壳体情形。本文编制的这种退化壳单元程序还可进行各向异性材料的应力应变分析[1 ] 。涡轮叶片因开有很多孔而使结构强度大大降低 ,由于孔间的相互影响 ,应力场相当复杂。我们用下列三种简化模型作计算分析。三种平板尺寸为 …     

三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。     

三维燃烧移动边界计算网格示踪技术

以任意设计的固体火箭发动机装药几何形状为对象,分析研究了燃面燃烧移动的基本规律,提出了一种三维复杂燃烧移动边界的示踪方法。它可以用于流场计算、发动机性能预示计算和侵蚀燃烧计算中,也可以发展成为一种燃面计算的新方法。文中以假想三维局部翼形装药为例进行了数值计算。     

三维装药通用计算方法——体积元素法

本文提出一种新的三维装药通用计算方法——体积元素法.这种计算方法使计算过程的数学处理十分简单、方便.根据该方法编制的计算机程序可以对各种类型的三维装药进行快速、准确的计算.用该程序对一台具有复杂燃面变化规律的实际发动机进行了计算,获得了令人满意的结果.     

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

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

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.     

跨音速平面叶栅全位势无粘-粘性迭代计算

本文发展了一种全位势与附面层、尾迹流正逆模式积分法迭代计算程序,给出了计算结果与国内外四个不同的叶栅风洞试验结果的对比分析。结果表明,对比是较为满意的,本文在网格生成,无粘跨音流求解和无粘流与附面层、尾迹流迭代等方面所作的发展和改进适用,使整个程序系统具有工程实用价值。      

跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

火箭发动机塞式喷管流场的数值研究

用时间相关法和张涵信的无振荡、无自由参数ＮＮＤ显式格式,对火箭发动机塞式喷管二维轴对称定常流场进行了数值模拟。计算得到了有粘和无粘中两种情况下,喷管流场的马赫数压强、温度、流线以及速度矢量分布图,计算网络采用代数方法生成,并通过求解Ｌａｐａｃｅ方程对网格进行修正。计算结果表明：ＮＮＤ格式用于计算塞式喷管流场具有较高的效率和良好的分辨率;无粘和有粘情况得到的流场结构有较大的差别;基础喷管的出口马赫数     

Mach reflection in steady supersonic flow considering wedge boundary-layer correction

Mach reflection in steady supersonic flow is an important phenomenon having received extensive studies, among which simplified theoretical models to predict the size of Mach stem and other flow structure are of particular interest. Past efforts for such models were based on inviscid assumption while in real cases the flow is viscous. Here in this paper we consider the influence of wedge boundary layer on the Mach stem height. This is done by including a simplified boundary layer model into a recently published inviscid model. In this viscous model, the wedge angle and the trailing edge height, which control the Mach stem height, are replaced by their equivalent ones accounting for the displacement effect of the wedge boundary layer, with the boundary layer assumed to be laminar or fully turbulent. This viscous model is shown to compare well with numerical results by computational fluid dynamics and gives a Mach stem height as function of the Reynolds number and Mach number. It is shown that due to the viscous effect, the Mach stem height is increased, through increasing the effective wedge angle.     

Viscosity effects on weak irregular reflection of shock waves in steady flow

The viscosity effects on weak shock wave reflection are investigated with the Navier–Stokes and DSMC flow solvers. It is shown that the viscosity plays a crucial role in the vicinity of three-shock intersection. Instead of a singular triple point, in viscous flow there is a smooth three shock transition zone, where one-dimensional shock jump relations cannot be applied. At the flow parameters corresponding to the von Neumann reflection, when no inviscid three-shock solution exists, the computations predict an irregular shock-wave configuration similar to that observed previously in experiments. The existence of a viscous zone in the region of shock-wave interaction allows a continuous transition from the parameters behind the Mach stem to the parameters behind the reflected shock, which is impossible in the inviscid three-shock theory. Results of numerical simulations agree qualitatively with the conclusions of the approximate viscous theory, which describes the flow in the vicinity of the three-shock intersection.     

DCS5在计算气动声学中的应用研究

对流项采用五阶线性耗散紧致格式(DCS5)离散,粘性项采用6阶中心紧致格式离散,时间项采用5级5阶龙格-库塔时间推进格式求解流体动力学方程组,并由此建立近声场的直接模拟方法。为了应用结构网格在复杂几何形状条件下对流场进行计算,发展了高精度的多块特征对接边界技术。完成了低速无粘振荡圆柱和粘性圆柱绕流等标准计算声学算例,获得了令人满意的结果。     

Euler方程与边界层积分方程耦合求解跨音速翼型绕流

 应用Euler方程求解跨音速翼型特性时考虑了粘性影响,粘性影响是通过边界层动量和能量积分方程求解的,即粘流/无粘流迭代方法。其中Euler方程采用LU-ADI方法求解;边界层方程均由正解法过渡到反解法,以解决强激波干扰区出现小分离泡的计算问题。计算中使用了贴体C网格,通过一定变换使其保持基本正交。计算结果表明,压力分布、摩阻系数分布与实验结果符合较好。     

不起动流场对超声速/高超声速进气道自起动性能的影响

对7个典型速域的二维超声速/高超声速进气道加速自起动过程进行了准定常数值仿真,分析了真实临界不起动流场对进气道自起动性能的影响,研究发现:存在一种介于超声速和高超声速临界不起动模式之间的过渡临界不起动模式。当真实不起动流场处于超声速临界不起动模式时,自起动马赫数略大于无黏设计自起动马赫数;处于过渡临界不起动模式时自起动马赫数小于无黏设计自起动马赫数;而该研究中处于高超声速临界不起动模式的进气道,自起动马赫数明显大于无黏设计自起动马赫数。高超声速临界不起动模式下的喉道截面特征气流参数显著偏离无黏临界不起动流场,所以Kantrowitz理论以及基于该理论发展而来的系列方法不适用于预测高超声速进气道自起动性能。