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

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

钝头翼型和机翼跨音速绕流的渐近展开分析

本文利用渐近展开匹配法分析钝头翼型的跨音速绕流,导出了描述前缘附近流动的一级近似、二级近似下的速位方程、边界条件及相应的近似解析解,并构成关于翼表面速度的一致有效合成解,消除了跨音速小扰动近似的前缘奇性,对于大展弦比后掠翼绕流,可利用翼型绕流分析结果,消除机翼前缘奇性。     

Euler flow solutions for transonic shock wave–boundary layer interaction

Steady 2D Euler flow computations have been performed for a wind tunnel section, designed for research on transonic shock wave–boundary layer interaction. For the discretization of the steady Euler equations, an upwind finite volume technique has been applied. The solution method used is collective, symmetric point Gauss–Seidel relaxation, accelerated by non-linear multigrid. Initial finest grid solutions have been obtained by nested iteration. Automatic grid adaptation has been applied for obtaining sharp shocks. An indication is given of the mathematical quality of four different boundary conditions for the outlet flow. Two transonic flow solutions with shock are presented: a choked and a non-choked flow. Both flow solutions show good shock capturing. A comparison is made with experimental results.     

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.     

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.     

Multiple solutions and stability of the steady transonic small-disturbance equation

Numerical solutions of the steady transonic small-disturbance(TSD) potential equation are computed using the conservative Murman-Cole scheme. Multiple solutions are discovered and mapped out for the Mach number range at zero angle of attack and the angle of attack range at Mach number 0.85 for the NACA 0012 airfoil. We present a linear stability analysis method by directly assembling and evaluating the Jacobian matrix of the nonlinear finite-difference equation of the TSD equation. The stability of all the discovered multiple solutions are then determined by the proposed eigen analysis. The relation of stability to convergence of the iterative method for solving the TSD equation is discussed. Computations and the stability analysis demonstrate the possibility of eliminating the multiple solutions and stabilizing the remaining unique solution by adding a sufficiently long splitter plate downstream the airfoil trailing edge. Finally, instability of the solution of the TSD equation is shown to be closely connected to the onset of transonic buffet by comparing with experimental data.     

跨音速大迎角绕流的Euler方程数值模拟

采用保角变换与代数方法相结合，生成全场统一的贴体、正交Ｏ－Ｈ型网格．采用有限体积法求解Ｅｕｌｅｒ方程，模拟具有歼击机外形的全机及翼身组合体大迎角跨音速绕流．计算表明，法向力系数、气动中心位置及压力分布的计算结果与实验结果吻合良好     

Asymptotic analysis of transonic flows

In this paper we derive the equations of the second and third approximations for the stream function of two-dimensional and axisymmetric potential transonic flow of an inviscid gas and find their particular solutions corresponding to certain transonic flows.A similar study concerning the second approximation of subsonic and supersonic flow was made by Van Dyke [1] and Hayes [2]. The second approximation for the velocity potential of transonic flow has been examined in detail by Hayes [3]. Euvrard [4, 5] has investigated the asymptotic behavior of transonic flow far from a body, while Fal'kovich, Chernov, and Gorskii [6] have studied the flow in a nozzle throat.The transonic asymptotic analysis for the stream function is presented in this paper.     

Symmetric Gauss-Seidel multigrid solution of the Euler equations on structured and unstructured grids

The efficient symmetric Gauss-Seidel (SGS) algorithm for solving the Euler equations of inviscid, compressible flow on structured grids, developed in collaboration with Jameson of Stanford University, is extended to unstructured grids. The algorithm uses a nonlinear formulation of an SGS solver, implemented within the framework of multigrid. The earlier form of the algorithm used the natural (lexicographic) ordering of the mesh cells available on structured grids for the SGS sweeps, but a number of features of the method that are believed to contribute to its success can also be implemented for computations on unstructured grids. The present paper reviews, the features of the SGS multigrid solver for structured gr0ids, including its nonlinear implementation, its use of “absolute” Jacobian matrix preconditioning, and its incorporation of multigrid, and then describes the incorporation of these features into an algorithm suitable for computations on unstructured grids. The implementation on unstructured grids is based on the agglomerated multigrid method developed by Sørensen, which uses an explicit Runge-Kutta smoothing algorithm. Results of computations for steady, transonic flows past two-dimensional airfoils are presented, and the efficiency of the method is evaluated for computations on both structured and unstructured meshes.     

Steady shock wave reflections in thermochemical nonequilibrium flows

The shock wave reflection phenomena in hypersonic steady air flows, including thermochemical nonequilibrium effects, are investigated. The main objectives are to study the influence of these effects on the two shock wave reflections (regular and Mach reflections), on the Mach stem height and on the hysteresis behavior. The air computations are performed using a multi-block MUSCL-TVD finite-volume scheme. The computational results with and without thermochemical effects in the air mixture flow at an upstream Mach number equal to 7 are compared. The comparison reveals a strong dependence of the transition angles, of the height and location of the Mach stem on the physical modeling of the gas flow. Received 17 February 2000 / Accepted 30 August 2000     

Computation of cascade flutter by uncoupled and coupled methods

A computational method for flutter prediction of turbomachinery cascades is presented. The flow through multiple blade passages is calculated using a time-domain approach with coupled aerodynamic and structural models. The unsteady Euler/Navier-Stokes equations are solved in quasi-three-dimensions using a second-order implicit scheme with dual time-stepping and a multigrid method. A structural model for the blades with bending and torsion degrees of freedom is integrated in time together with the flow field. Information between structural and aerodynamic models is exchanged until convergence in each real-time step. Computational results for a cascade are presented and compared with those obtained by the conventional energy method and with experimental and numerical data by other authors. Significant differences are found between the coupled and uncoupled methods at low mass ratios. A transonic test case with strong nonlinear phenomena is investigated with the fluid-structure coupled method. Results for inviscid flow are compared with results of Navier-Stokes computations.     

二维定常湍流计算中的GMRES算法

在以前工作的基础上将广义极小残差（Generalized Minimum RESidual (GMRES)算法发展到用于求解二维可压Favier平均Navier-Stokes方程组。控制方程经Newton线化处理后构成近似的线性系统,然后采用分别耦合了LUSGS和ILU两种预处理矩阵的GMRES算法求解。Spalart-Allmaras湍流模型被用来封闭流体控制方程组,采用与流体控制方程非耦合的方式,使用LUSGS方法求解。对GMRES算法中矩阵－向量的乘积采用了有限插分方法,从而避免了精确的左端系数矩阵的计算和存储。对预处理矩阵的两种使用方法（左预处理和右预处理）进行了分析和讨论。用两个算例对LUSGS和ILU两各预处理矩阵进行了比较,同时探讨了左预处理和右预处理各自的优缺点。通过对Sajben扩压器和NACA0012有攻角流动的计算,表明带有预处理的GMRES算法在二维定常跨音黏性流动计算中相比于得到广泛应用的DDADI方法具有很大优势,左预处理要优于右预处理。     

An efficient implicit mesh‐less method for compressible flow calculations

A dual‐time implicit mesh‐less scheme is presented for calculation of compressible inviscid flow equations. The Taylor series least‐square method is used for approximation of spatial derivatives at each node which leads to a central difference discretization. Several convergence acceleration techniques such as local time stepping, enthalpy damping and residual smoothing are adopted in this approach. The capabilities of the method are demonstrated by flow computations around single and multi‐element airfoils at subsonic, transonic and supersonic flow conditions. Results are presented which indicate good agreements with other reliable finite‐volume results. The computational time is considerably reduced when using the proposed mesh‐less method compared with the explicit mesh‐less and finite‐volume schemes using the same point distributions. Copyright © 2010 John Wiley & Sons, Ltd.     

A SOLVER USING NEWTON‘S METHOD FOR UNSTEADY VISCOUS FLOWS

A solver is developed for time-accurate computations of viscous flows based on the conception of Newton‘s method. A set of pseudo-time derivatives are added into governing equations and the discretized system is solved using GMRES algorithm. Due to some special properties of GMRES algorithm, the solution procedure for unsteady flows could be regarded as a kind of Newton iteration. The physical-time derivatives of governing equations are discretized using two different approaches, I.e., 3-point Euler backward, and Crank-Nicolson formulas, both with 2nd-order accuracy in time but with different truncation errors. The turbulent eddy viscosity is calculated by using a version of Spalart~Allmaras one-equation model modified by authors for turbulent flows. Two cases of unsteady viscous flow are investigated to validate and assess the solver, I.e., low Reynolds number flow around a row of cylinders and transonic bi-circular-arc airfoil flow featuring the vortex shedding and shock buffeting problems, respectively. Meanwhile, comparisons between the two schemes of timederivative discretizations are carefully made. It is illustrated that the developed unsteady flow solver shows a considerable efficiency and the Crank-Nicolson scheme gives better results compared with Euler method.     

A meshless finite point method for three‐dimensional analysis of compressible flow problems involving moving boundaries and adaptivity

A finite point method for solving compressible flow problems involving moving boundaries and adaptivity is presented. The numerical methodology is based on an upwind‐biased discretization of the Euler equations, written in arbitrary Lagrangian–Eulerian form and integrated in time by means of a dual‐time steeping technique. In order to exploit the meshless potential of the method, a domain deformation approach based on the spring network analogy is implemented, and h‐adaptivity is also employed in the computations. Typical movable boundary problems in transonic flow regime are solved to assess the performance of the proposed technique. In addition, an application to a fluid–structure interaction problem involving static aeroelasticity illustrates the capability of the method to deal with practical engineering analyses. The computational cost and multi‐core performance of the proposed technique is also discussed through the examples provided. Copyright © 2013 John Wiley & Sons, Ltd.     

The mathematical simulation of various condensation regimes in turbulent isobaric jets

Mathematical models are considered and calculations made for flows in turbulent isobaric steam—air jets in the presence of condensation of the water vapor they contain. The models consist of gasdynamic equations for a turbulent jet, equations for a differential two-parameter model of turbulence, thermodynamic relations, and kinetic equations. A study is made of steam—air jets in a regime of condensation in equilibrium, when the flow region is broken down into zones of frozen flow and flow in equilibrium, described by the equations for a turbulent jet with the use of the traditional thermodynamic relations and of the thermodynamic relations for condensation in equilibrium. An analysis is made of the influence of pulsating motion on the kinetic parameters: rate of nucleation, the critical size of the nuclei, and rate of growth of the drops. It is shown that the rate of nucleation, determined from a quasilaminar averaging model, is several orders of magnitude less than the mean value obtained by averaging using the density distribution of the passive admixture concentration probability. A numerical study is made of the heterogeneous condensation in turbulent jets on extraneous particles entering from the nozzle. Kinetic equations are written down for the case when the rate of growth of the drops does not depend on their radius. A study is made of the dynamics of the transition of heterogeneous condensation from disequilibrium to equilibriumTranslated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 59–67, January–February, 1985.     

Finite volume method for three-dimensional transonic potential flow through turbomachinery blade rows

A finite-volume method has been developed for the calculation of transonic, potential flows through 3-D turbomachinery blades with complex geometries. The exact transonic potential flow equation is solved on a mesh constructed from small volume elements. A transformation is introduced through which cuboids of the physical plane are mapped into computational cubes. Two sets of overlapping volumes are used. While the thermodynamic properties are calculated at the primary volume centres, the flux balance is established on the secondary volumes. For transonic flows an artificial compressibility term (upwind density gradient) is added to density to produce the necessary directional bias in the hyperbolic region. The successive point over-relaxation Gauss-Seidel method has been used to solve the non-linear partial differential equations. Comparisons with experiments and/or other numerical solutions for various turbomachinery configurations show that the 3-D finite-volume approach is a relatively accurate, reliable and fast method for inviscid, transonic flow predictions through turbomachinery blade rows     

高超声速三维化学非平衡流动N－S方程数值解

详细地叙述了模拟三维高超声速电离空气化学非平衡粘性流全流场（前体和底部近尾迹）的数值方法，化学反应粘性流求解是基于带有化学源项的Ｎ－Ｓ方程的守恒形式，总的连续方程由组元守恒方程组所代替．对于高温电离空气流动，存在如下７个主要组元，它们依次为Ｎ＿２，Ｏ＿２，ＮＯ，ＮＯ￣＋，Ｎ，Ｏ和ｅ￣－．研究发展了求解完全耦合的，与时间相关的偏微分方程组的数值方法．利用一类新的迎风ＴＶＤ激波捕捉有限差分格式求解守恒形式的控制方程组，对高超声速层流有攻角钝锥体绕流流场，在非催化壁条件下、对包括底部近尾迹在内的全流场各参数，如组元浓度，电子密度和温度，以及物面压力分布和热流率分布得出了计算结果．对化学反应气体和完全气体模型的计算结果进行了比较，计算的流场电子密度与飞行实验结果符合很好．     

Performance Improvement of a Supercritical Airfoil by a Multi-point Optimized Shock Control Channel

A shock control channel (SCC) is a flow control method introduced here to control the shock wave/boundarylayer interaction (SWBLI) in order to reduce the resulting wave drag in transonic flows. An SCC transfers an appropriate amount of mass and momentum from downstream of the shock wave location to its upstream to decrease the pressure gradient across the shock wave and as a result the shock-wave strength is reduced. Here, a multi-point optimization method under a constant-lift-coefficient constraint is used to find the optimum design of the SCC. This flow control method is implemented on a RAE-2822 supercritical airfoil for a wide range of off-design transonic Mach numbers. The RANS flow equations are solved using Roe’s averages scheme and a gradient-based adjoint algorithm is used to find the optimum location and shape of the SCC. The solver is validated against experimental works studying effect of cavities in transonic airfoils. It is shown that the application of an SCC improves the average aerodynamic efficiency in a range of off-design conditions by 13.2% in comparison with the original airfoil. The SCC is shown to be an effective tool also for higher angle of attack at transonic flows. We have also studied the SWBLI and how the optimization algorithm makes the flow wave structure and interactions of the shock wave with the boundary layer favorable.