湿空气非平衡凝结对压气机叶栅性能影响的数值研究

跨音速流动条件下湿空气中的水蒸气由于快速膨胀而发生非平衡凝结,凝结潜热与跨音速气流耦合会显著改变流场结构。本文建立了湿空气非平衡凝结流动的数学物理模型,对三种工况下压气机叶栅内湿空气凝结流动进行计算分析。计算结果表明:湿空气压气机叶栅流动中,湿空气中水蒸气有可能发生非平衡凝结从而对叶栅通道内的流动产生强烈影响,改变流场内的激波结构,叶栅内压力系数减小最大幅度由1.38减小到0.76。     

喷管内水蒸气非平衡凝结流动数值模拟研究

基于水蒸气非平衡凝结动力学模型和水蒸气真实热力学性质模型,建立了水蒸气跨声速非平衡凝结流动数值模型,对喷管内水蒸气非平衡凝结流动进行了数值模拟,模型成功地捕捉到了凝结激波,并且模拟结果与实验数据吻合较好,证明了所建立数值模型的准确性和可靠性。在此基础上研究了水蒸气在跨声速流动过程中发生的非平衡凝结现象以及对流场产生的影响,发现水蒸气凝结现象对流场的影响非常显著;并对凝结激波的产生过程和原因进行了详细分析。     

跨音速有旋流动正问题的赝势函数变分有限元法

本文运用赝势函数变分有限元方法数值模拟了绕翼型的跨音速有旋流动。在含有激波的跨音速有旋流动中，势函数已不存在，但为了保留势函数模型在求解方面的优越性，上海大学的刘高联引入了一个通用函数一赝势函数，可以看出该赝势函数保持了势函数的所有好处，又突破了流动有势的限制，是势函数对有旋流动的一个自然的、物理上相容、数学上求解简便的推广，进一步地，刘高联还得到了赝势函数的变分原理族，为变分有限元法求解有旋流动打下了基础。另一方面，为了提高数值求解的收敛性和有效地捕捉流场中的激波，本文还采用了“人工密度”办法。绕翼型的跨音速有旋流动的计算实践证明了赝势函数的有效性。     

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

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

流体动力系统激波的数值模拟

本文为了对流体动力系统含有激波的跨音速流动进行数值模拟,给出了曲线坐标系下的LU-AUSMLW算法。该算法综合利用了AUSMPW格式与LU-SGS方法的优点。本文通过采用三阶MUSCL格式获得高阶精度。为了验证该算法的精度,对通道和叶栅中含有激波的跨音速流动进行了数值模拟。本文数值试验的计算结果与文献计算结果相符很好。     

高温激波管化学非平衡流动数值模拟研究

为了预测氢氧定容燃烧驱动的高温激波管性能,需要准确分析激波管非定常化学非平衡流动过程.本文在破膜前的驱动段定容燃烧以及破膜后的化学非平衡流动数值模拟中,引入双时间步长方法,发展高温激波管化学非平衡流动数值模拟方法,该方法在时间上具有二阶精度.计算结果与目前存在的激波管流动解析解以及零维化学反应系统的数值解进行了比较,吻合较好.对于典型高温激波管状态,采用有限体积方法离散准一维流动Euler控制方程,并通过将流动过程和化学反应动力学过程耦合求解,获得了激波管内部的化学非平衡流动特征.     

以压力为基本求解变量数值模拟粘性超、跨音速流动

应用以压力为基本求解变量的SIMPLE方法 ,对一双喉喷管中的层流超音速流动和一扩压器中的紊流跨音速流动进行了数值计算。计算结果显示 ,本文的计算结果与文献数据及实验结果相符很好。表明本文方法对可压缩流动有很高的模拟精度。进而表明经过可压缩推广的SIMPLE方法适用于任何马赫数的流动计算     

TVD格式在超音速喷管三维粘性流动求解中的应用

详细给出了任意三维曲线坐标系中Ｎｏｖｉｅｒ－Ｓｔｏｋｅｓ方程的对流项ＴＶＤ格式的构造过程，建立了数值求解三维粘性流动的计算方法，应用该方法对三维超音速喷管中有激波及无激波情况下的两种工况的层流流场进行了数值求解，并与实验做了对比。结果表明本文建立的计算方法具有较高的精度，同时也证明ＴＶＤ格式具有分辩率高，稳定收敛等优点，为进一步开展叶栅流场及紊流的研究打下了基础。     

跨音速翼型上的激波／边界层干扰自适应控制计算

在激波区使用自适应壁对跨音速翼型的激波／边界层的相互作用（干扰）进行控制，可改变机翼的气动性能，这种被动控制可通过在翼型的激波区开一凹腔，其上覆盖一弹性橡胶膜柔壁来，本文给出用Ｎａｖｉｅｒ－Ｓｔｏｋｅｒ方程数值模拟这一自适应控制翼型的跨音速粘性绕流，提出了一个适应于本特殊情况（物面边界局部地区在求解过程中有变化）的处理办法。并探讨了自适应柔壁对当代跨音速翼绕流的影响。     

翼型、机翼非定常运动模式下动态数值模拟

基于RANS方程,通过刚性动网格技术实现对翼型和机翼典型运动模式的描述,采用双时间推进方法和Roe空间离散格式对流场求解,构建了一个非定常气动计算平台;以NACA0012翼型为算倒进行了动态数值模拟可信度验证。数值模拟结果与试验数据吻合较好,升力和俯仰力矩的最大计算误差分别为3％和10％,表明了该平台的可靠性。另外,还数值模拟了M6机翼的动态非定常流场,并分析了两种湍流模型对非定常流场激波的捕捉能力。结果表明非定常流动中S-A湍流模型对激波的捕捉较B-L模型更敏感。文中开发的非定常计算平台对进一步解决三维复杂流场的流动问题有很高的工程应用价值。     

Theoretical and Numerical Studies of Transonic Flow of Moist Air Around a Thin Airfoil

Numerical studies of a two-dimensional and steady transonic flow of moist air around a thin airfoil with condensation are presented. The computations are guided by a recent transonic small-disturbance (TSD) theory of Rusak and Lee (2000) on this topic. The asymptotic model provides a simplified framework to investigate the changes in the flow field caused by the heat addition from a nonequilibrium process of condensation of water vapor in the air by homogeneous nucleation. An iterative method which is based on a type-sensitive difference scheme is applied to solve the governing equations. The results demonstrate the similarity rules for transonic flow of moist air and the effects of energy supply by condensation on the flow behavior. They provide a method to formulate various cases with different flow properties that have a sufficiently close behavior and that can be used in future computations, experiments, and design of flow systems operating with moist air. Also, the computations show that the TSD solutions of moist air flows represent the essence of the flow character computed from the inviscid fluid flow equations. Received 5 October 2000 and accepted 21 March 2002     

Numerical study of the spontaneous nucleation of self-rotational moist gas in a converging–diverging nozzle

Spontaneous nucleation is the primary way of droplet formation in the supersonic gas separation technology, and the converging–diverging nozzle is the condensation and separation unit of supersonic gas separation devices. A three-dimensional geometrical model for the generation of self-rotational transonic gas flow is set up, based on which, the spontaneous nucleation of self-rotational transonic moist gas in the converging–diverging nozzle is carried out using an Eulerian multi-fluid model. The simulated results of the main flow and nucleation parameters indicate that the spontaneous nucleation can occur in the diverging part of the nozzle. However, different from the nucleation flow without self-rotation, the distributions of these parameters are unsymmetrical about the nozzle axis due to the irregular flow form caused by the self-rotation of gas flow. The nucleation region is located on the position where gas flows with intense rotation and the self-rotation impacts much on the nucleation process. Stronger rotation delays the onset of spontaneous nucleation and yields lower nucleation rate and narrow nucleation region. In addition, influences of other factors such as inlet total pressure p ₀, inlet total temperature T ₀, the nozzle-expanding ratio ? and the inlet relative humidity ф ₀ on the nucleation of self-rotational moist gas flow in the nozzle are also discussed.     

Nonlinear dynamics of an aeroelastic airfoil with free-play in transonic flow

Nonlinear dynamic behaviors of an aeroelastic airfoil with free-play in transonic air flow are studied. The aeroelastic response is obtained by using time-marching approach with computational fluid dynamics (CFD) and reduced order model (ROM) techniques. Several standardized tests of transonic flutter are presented to validate numerical approaches. It is found that in time-marching approach with CFD technique, the time-step size has a significant effect on the calculated aeroelastic response, especially for cases considering both structural and aerodynamic nonlinearities. The nonlinear dynamic behavior for the present model in transonic air flow is greatly different from that in subsonic regime where only simple harmonic oscillations are observed. Major features of the responses in transonic air flow at different flow speeds can be summarized as follows. The aeroelastic responses with the amplitude near the free-play are dominated by single degree of freedom flutter mechanism, and snap-though phenomenon can be observed when the air speed is low. The bifurcation diagram can be captured by using ROM technique, and it is observed that the route to chaos for the present model is via period-doubling, which is essentially caused by the free-play nonlinearity. When the flow speed approaches the linear flutter speed, the aeroelastic system vibrates with large amplitude, which is dominated by the aerodynamic nonlinearity. Effects of boundary layer and airfoil profile on the nonlinear responses of the aeroelastic system are also discussed.     

Exergy efficiency of a counterflow air/air heat exchanger with vapour condensation

The exergy efficiency of a counterflow hot moist air/ambient air heat exchanger with vapour condensation as a function of the temperature and humidity of the hot moist air, the ratio of ambient air flow rate to hot moist air flow rate and the minimum temperature differences between the flows is determined. The ranges 20 to 120°C and 0.01 to 1.0 kg·kg^–1 for the temperature and humidity ratio of hot moist air respectively are considered. Constant ambient air conditions of 20°C and 0.01 kg·kg^–1 were assumed.

---

Exergie-Wirkungsgrad eines Luft/Luft-Gegenstrom-Wärmetauschers mit Dampfkondensation

Zusammenfassung Der Exergie-Wirkungsgrad eines Gegenstrom-Wärmetauschers wird bestimmt für die Kombination: heiße Feuchtluft zu Umgebungsluft mit Dampfkondensation, und zwar als Funktion von Temperatur und Feuchtegrad der heißen Feuchtluft, des Verhältnisses von Umgebungsluft zu heißen Feuchtluftstrom und der minimalen Temperaturdifferenz zwischen beiden Strömen. Die Temperatur des heißen Feuchtluftstroms variiert im Bereich von 20 bis 120°C, sein Feuchtegrad im Bereich 0,01 bis 1,0. Bezüglich der Umgebungsluft wurden konstante Bedingungen angenommen, und zwar 20°C und ein Feuchtegrad von 0,01.

     

Steam condensing flow modeling in turbine channels

The numerical method for modeling of the transonic steam flows with homogeneous and/or heterogeneous condensation has been presented. The experiments carried out for the Laval nozzles, for 2-D turbine cascades and for a 3-D flow in real turbine were selected to validate an in-house CFD code adjusted to the calculations of the steam condensing flows in complicated geometries. The sensitivity of the condensation model and difficulties in the validation process of the CFD code have been discussed. These difficulties limit the possibilities of verification and improvement of the condensation theory based on the existing experimental data.     

Multi-dimensional simulation of thermal non-equilibrium channel flow

The Homogenous Relaxation Model (HRM) is used to study thermal non-equilibrium, two-phase flows with flash-boiling and condensation. Typically, such non-equilibrium phase-change models have been studied in one-dimensional flow, but the goal of the present work is to create and utilize a multi-dimensional CFD implementation. The simulations are able to handle general polyhedral meshes, an important convenience for irregular channel or nozzle shapes. The model is applied to flash-boiling flow in short channels and validated against experimental measurements. The simulations predict the multi-dimensional features that have been observed in the past in experiments. Nozzle choking is also observed in the calculations.     

Numerical investigation of the effect of inlet condition on self‐excited oscillation of wet steam flow in a supersonic turbine cascade

Self‐excited oscillation can be induced due to the interaction between condensation process and local transonic condition in condensing flow, which is an important problem in wet steam turbine. With an Eulerian/Eulerian numerical model, the self‐excited oscillation of wet steam flow is investigated in a supersonic turbine cascade. Owing to supercritical heat addition to the subsonic flow in the convergent part of the cascade, the oscillation frequency decreases with increased inlet supercooling. Mass flow rate increases in the oscillating flow due to the greater supersaturation in condensation process, while the increase will be suppressed with the flow oscillation. Higher inlet supercooling leads to the fact that the condensation process moves upstream and the loss increases. Moreover, some predictions of oscillation effects on outflow angle and aerodynamic force are also presented. Finally, heterogeneous condensations with inlet wetness and periodic inlet conditions, as a result of the interference between stator and rotor, are discussed. Copyright © 2008 John Wiley & Sons, Ltd.     

k-ω SST湍流模式三维激波分离流动修正

“激波?边界层分离”是航空气动领域的典型湍流非平衡流动问题, 准确模拟激波分离对于跨声速飞行器气动性能评估和优化设计具有重要意义. 然而传统涡黏性湍流模式中涡黏性系数的定义方式并不适用于非平衡流动, k-ω SST湍流模式为此引入的Bradshaw假设在应用于三维强逆压梯度和较大分离流动时反而限制了雷诺应力的生成, 导致包括k-ω SST在内的常用涡黏性湍流模式均无法对此类流动进行准确模拟. 同时, 现有的非线性雷诺应力本构关系也并不能有效提高模拟精度. 为此, 针对k-ω SST模式分别提出了基于Bradshaw假设和基于长度尺度的两种激波分离流动修正方法. 前者通过提高Bradshaw常数的方式放宽了对雷诺应力生成的限制, 后者则从湍流长度尺度概念出发, 利用混合长度理论、湍动能生成/耗散之比和一种新定义的长度尺度之比构造了ω方程耗散项修正函数, 提高了模式在三维激波分离流动中的建模长度尺度. 两种方法对ONERA M6机翼跨声速大攻角流动均能得到较雷诺应力模式更好的模拟结果. 进一步的雷诺应力分析表明, 三维激波分离流动中“主雷诺应力分量”的概念不再成立, 各雷诺应力分量大小接近. 网格收敛性分析、对其他攻角状态的验证以及湍流平板边界层壁面律验证进一步确认了所提出的两种修正方法的合理性、有效性和通用性.      

Specific gasdynamical features of spontaneous condensation in an unsteady rarefaction wave

On the basis of numerical modeling, the formation of an unsteady shock wave induced by a condensation shock in a rarefaction wave moving in the high-pressure channel of a shock tube filled with moist air is demonstrated. It is shown that in a fairly long channel a periodic structure consisting of an alternating sequence of condensation shocks and the shock waves they generate may be formed. This structure is a linear unsteady analog of the self-oscillation regime of type IV in the classification [1] for condensing medium flows in the subsonic section of a Laval nozzle. The specific features detected are important for planning and interpreting experiments aimed at investigating spontaneous condensation using a “condensation shock tube”.