一先进跨音速离心压气机叶轮内部粘性流动数值模拟

随着计算机和计算技术的发展，三元粘性数值计算被广泛地应用于叶轮机械内部流场分析，本文使用三元Ｎ－Ｓ方程对这一先进跨音速离心压气机叶轮内部流场进行数值模拟，一是对这一成功设计进行分析验证。二是使用数值方法进行内部流场细节的模拟掌握其内部粘性流动特点，为今后设计和改进高压比跨音速离心压气机提供经验。     

跨音速湿空气非平衡凝结流动的数值研究

跨音速流动条件下湿空气中的水蒸气由于快速膨胀而发生非平衡凝结,凝结潜热对跨音速气流进行加热,会显著改变气流的流动特性。通过对商用计算流体动力学软件FLUENT进行二次开发,建立了湿空气非平衡凝结流动的数值求解方法。该方法可用于二维或三维、粘性或无粘、内流或外流的求解中。采用该方法分剐对缩放喷管、透平叶栅以及绕CA-0.1圆弧翼型的湿空气非平衡凝结流动进行了数值分析。计算结果表明：湿空气凝结手l起缩放喷管中的凝结激波、导致叶橱流动中总压降低;对于翼型周围的流动,在相对湿度分别为50％、57．1％、64.1％时,依次计算得到了单激波、五激波、双激波。     

某跨音转子设计中叶型中线分布特征研究

本文通过调整叶型中线分布设计了几种不同的跨音转子,并采用数值模拟的方法研究了其对跨音转子流场性能的影响。研究结果表明,高来流马赫数下跨音转子流场对叶型中线分布比较敏感,中线分布应尽可能保证负荷的平稳加载和预期的气流折转;跨音转子叶根区域仅在前缘附近存在局部超音区,较大的叶型折转对流场影响不大;其叶尖区域主要依靠激波增压...     

Triple-deck analysis of formation of supersonic and local separation regions in transonic steady flow over a roughness element on the surface of a body of revolution

Transonic axisymmetric flow over a body of rotation with a small roughness element located on its surface is considered. The body is manly cylindrical. The roughness height is assumed to be much smaller than the radius of the cylinder and such that a triple-deck flow is induced in its neighborhood. The goal of the work is to study the effect of the cylinder radius and the roughness shape on the triple-deck flow when the cylinder radius is of the same order as the transverse size of the triple-deck interaction region. In this case, the effect of three-dimensionality of the flow is exhibited even in the first approximation. Special attention is given to the structure of supersonic regions and closing shock waves arising in the outer potential region, as well as to local separation regions if they develop in the lower viscous boundary sublayer. Specifically, it is shown that, as the radius of the cylinder increases at a fixed roughness height, the shock intensity grows considerably, whereas the position of the main shock varies little.     

三维机翼跨音速颤振的分布式并行计算

以多台微机组成分布式计算系统,利用无限插值方法(TFI)生成三维多块贴体运动网格,以Navier Stokes方程为控制方程,求解三维机翼的跨音速非定常气动力,并与颤振方程耦合迭代计算,求解机翼广义位移响应的时间历程,根据广义位移的时间历程的衰减、等幅和发散振荡等情况确定机翼跨音速颤振临界条件.通过算例验证,计算结果与实验结果和理论分析相吻合.     

跨声速平面叶栅多工况点反命题变分理论：人工来流振荡模型

本文首先提出一个新的‘人工来流振荡’模型,在此基础上建立了跨声速叶栅多工况点气动反命题的变域变分理论,可以保证叶栅在相当广阔的变工况范围内都能保持优良的运行特性。本理论的突出优点是可较直捷地推广到全三维流动和有旋流动中去,因而具有广阔的发展与应用前景。     

通量限制器对算法性能的影响

在非正交曹线坐标系下,本文给出了求解非线性双曲型Ｅｕｌｅｒ方程的ＬＵ－ＡＵＳＭＬＷ算法。为了改进该算法的性能,将高分辨率ＡＵＳＭＰＷ格式的空间精度由一阶精度扩展到三阶精度。分析了选择通量限制器对算法稳定性、收敛性和精度的影响,并构造了一种新的能量限制器。本文数值结果表明,通量限制器是决定ＬＵ－ＡＵＳＭＰＷ算法性能的关键因素,并且该算法采用本文构造的通量限制器,求解非线性双曲型Ｅｕｌｅｒ方程,具有较     

非定常跨音速流显式与隐式数值分析方法的比较研究

分别运用显式有限体积法与隐式LU-NND方法在同一网格与相同参数条件下对NACA0012翼型跨音速振荡绕流进行了数值分析。并就两种方法的计算结果、计算效率、对计算机内存的要求以及对网格质量的敏感性等做了对比研究。     

Computation of unsteady transonic flow over a fighter wing using a zonal Navier–Stokes/full‐potential method

An improved hybrid method for computing unsteady compressible viscous flows is presented. This method divides the computational domain into two zones. In the inner zone, the Navier–Stokes equations are solved using a diagonal form of an alternating‐direction implicit (ADI) approximate factorisation procedure. In the outer zone, the unsteady full‐potential equation (FPE) is solved. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. Characteristic‐based viscous/inviscid interface boundary conditions are employed to avoid spurious reflections at that interface. The resulting CPU times are about 60% of the full Navier–Stokes CPU times for unsteady flows in non‐vector processing machines. Applications of the method are presented for a F‐5 wing in steady and unsteady transonic flows. Steady surface pressures are in very good agreement with experimental data and are essentially identical to the full Navier–Stokes predictions. Density contours show that shocks cross the viscous/inviscid interface smoothly, so that the accuracy of full Navier–Stokes equations can be retained with significant savings in computational time. Copyright © 1999 John Wiley & Sons, Ltd.     

ROBUST NUMERICAL METHODS FOR TRANSONIC FLOWS

In this paper, numerical methods for solving the transonic full potential equation are developed. The governing equation is discretized by a flux-biasing finite volume method. The resulting non-linear algebraic system is solved by using a continuation method with full Newton iteration. The continuation method is based on solving a highly ‘upstream-weighted’ discretization and then gradually reducing the upstream weighting. A general PCG-like sparse matrix iterative solver is used to solve the Jacobians at each non-linear step. Various types of incomplete LU (ILU) preconditioners and ordering techniques are compared. Numerical results are presented to demonstrate that these methods are efficient and robust for solving the transonic potential equation in the workstation computing environment. © 1997 by John Wiley & Sons, Ltd.