小型可调Re数弹道靶上的跨声速实验

本文介绍了一种可调 Re 数小型弹道靶设备;讨论了该设备实验中的一些问题.由此开展了两种尖头细长旋成体不同 Re 数下的跨声速实验,并给出了部分流场显示结果.     

Controlling Transonic Flow around Airfoils by Means of Local Pulsed Addition of Energy

The influence of local pulsed-periodic addition of energy into a supersonic region on the flow structure and wave drag of an airfoil in transonic flow regimes is considered by methods of mathematical modeling. The study reveals significant prospects of the considered method of controlling airfoil performance in transonic flow regimes, including wave-drag reduction.     

基于气动力降阶模型的跨音速气动弹性稳定性分析

基于离散型输入输出差分模型,运用非定常CFD方法训练信号,然后运用最小二乘方法进行参数辨识,得到降阶的非定常气动力模型,再将该离散差分模型转换为连续时间域内的状态方程。耦合气动状态方程和结构状态方程,得到耦合系统的气动弹性状态方程。求解不同动压下状态矩阵的特征值,根据根轨迹图分析系统的稳定性特性。分析结果与直接耦合CFD/CSD方法结果相吻合,可以计算跨音速非线性气动弹性问题。其计算效率比直接耦合CFD/CSD方法提高1~2个数量级。针对Isogai wing在跨音速出现的S型颤振边界进行了较为细致的分析,阐述了该现象是由于系统诱发颤振的分支随着速度(来流动压)的提高而发生转移所导致的。     

横向三维化学反应流动数值模拟

对一个三维化学反应流程序作了改进,使之适应于模拟带横向碘蒸汽注入的COIL(化学氧碘激光)喷管流动。其中,主副气流的入流边界条件的改进是基本的,重要的。这里,入流边界附近的气体流动被视为一维均熵流动,用双曲方程的特征理论,分析以流量作变量的EULER方程的特征关系,用特征方法数值求解简化的特征方程,从而确定亚声速(跨声速)流动边界点的入流密度或速度。对相关文献提供的算例的计算,佐证了程序改进的有效性。对氮气作载气的条件也作了模拟,得出了与氦气条件不同的流动特征与增益分布特征。     

A variational approach to transonic potential flow problems

This paper continues considerations of transonic potential flow problems by variational methods. A functional which is associated with a boundary value problem for the (full) potential equation and which possesses a real physical meaning is minimized over a class of admissible functions. These functions have to satisfy a non‐linear local entropy condition and a certain boundness constraint. Though this class is not a compact set of the underlying Hilbert space and though the functional need not be convex, the existence of a solution to the established variational problem can be proved by direct methods of the calculus of variations. Furthermore, some properties of minimizers concerning uniqueness, relation to the potential equation, and behaviour on supersonic regions are derived. Copyright © 2000 John Wiley & Sons, Ltd.     

A matrix‐free preconditioned Newton/GMRES method for unsteady Navier–Stokes solutions

The unsteady compressible Reynolds‐averaged Navier–Stokes equations are discretized using the Osher approximate Riemann solver with fully implicit time stepping. The resulting non‐linear system at each time step is solved iteratively using a Newton/GMRES method. In the solution process, the Jacobian matrix–vector products are replaced by directional derivatives so that the evaluation and storage of the Jacobian matrix is removed from the procedure. An effective matrix‐free preconditioner is proposed to fully avoid matrix storage. Convergence rates, computational costs and computer memory requirements of the present method are compared with those of a matrix Newton/GMRES method, a four stage Runge–Kutta explicit method, and an approximate factorization sub‐iteration method. Effects of convergence tolerances for the GMRES linear solver on the convergence and the efficiency of the Newton iteration for the non‐linear system at each time step are analysed for both matrix‐free and matrix methods. Differences in the performance of the matrix‐free method for laminar and turbulent flows are highlighted and analysed. Unsteady turbulent Navier–Stokes solutions of pitching and combined translation–pitching aerofoil oscillations are presented for unsteady shock‐induced separation problems associated with the rotor blade flows of forward flying helicopters. Copyright © 2000 John Wiley & Sons, Ltd.     

翼型跨音速流动计算和外形设计的变域变分有限元方法

在气体动力学问题研究中经常会碰到诸如激波、翼型设计等未知界面问题。未知界面的存在为该类问题的理论分析和数值求解带来了很大困难。刘高联针对未知界面问题发展了一种变域变分有限元方法,该方法将未知界面看作是一个变化区域的边界,采用变域变分将未知界面结合在变分泛函中,使其与求解流场的控制方程结合起来,从而将未知界面的求解和流场的求解完全耦合进行,因而是一种处理未知界面的独特工具,极适合于气动外形的设计求解。本文运用变域变分有限元方法对翼型跨音速流动正、反命题进行了数值研究。由于在跨音速翼型绕流中存在激波,所以为了得到压缩激波解,采用了“人工密度”办法。几个算例均得到了满意的计算结果和设计结果,证明了本文方法的有效性和优越性。     

对二维分离流涡黏性系数非线性分布的新认识

以弱非线性涡黏性模型为出发点,对Delery分离流动实验结果进行分析并获得了非平衡态分离区涡黏性系数与形状因子J之间的非线性关系. 该非线性关系显示在分离起始阶段,涡黏性系数较平衡态先减小,后增大;再附阶段,涡黏性系数较平衡态数值逐渐增大,并在再附点位置接近最大,而后又逐渐减小,恢复到平衡态水平. 总结涡黏性系数的这种非线性发展数学关系式,并将它应用于BL模型,在不添加微分方程的情况下发展出一种适用于分离流动的改进代数湍流模型. 对低速平板流动,跨声速,超声速以及高超声速分离流动的计算结果表明,该改进湍流模型可以较准确地模拟各类复杂分离流动,计算精度明显优于传统代数模型以及一些两方程模型,而计算工作量仍与BL模型相当. 这表明所提出的涡黏性系数非线性发展规律是正确的,且应用在二维分离流动中具有一定的普适性.      

Aerodynamic drag of bodies in transonic flow. Theory and applications to computational aerodynamics

Formulas for all the components of the aerodynamic drag (total, friction, inductive, wave, pressure, and heat-transfer) are uniformly derived as applied to flows governed by the Navier-Stokes and Reynolds equations. For flows of this type the definition of the aerodynamic drag components is refined and the physical basis of the chosen method of breaking up the total drag into components is discussed. Ways of calculating the aerodynamic drag components using the methods of computational aerodynamics are considered. On the basis of the refined formulas the drag components are calculated for flows around airfoils and a high-aspect-ratio wing in transonic flow.     

翼身组合体气动力计算研究

采用求解Euler方程结合附面层修正的方法在结构网格上对翼身组合体跨音速流场进行了数值模拟.附面层方程的求解应用Whitfield提出的动量积分方程和平均流动能积分方程,为了保持Euler方程求解过程中计算网格的固定性,用加在物面上的溢出速度来模拟附面层效应.针对传统的近场方法计算阻力,计算精度较低、误差较大并且不能给出各阻力分量值的缺点,将基于动量定理的远场方法用于飞机的阻力估算,采用远场法将阻力分解为:粘性阻力,激波阻力,诱导阻力,并对各个分量分别进行了求解,将计算结果与近场法以及风洞实验值做了比较.以DLR-F4翼身组合体为考核算例,对所述方法进行了验证,结果显示远场法的计算结果与风洞实验值吻合的很好.