A finite volume method for two-dimensional transonic potential flow through turbomachinery blade rows

To predict inviscid transonic flow through turbomachinery blade rows, the exact transonic potential flow equation is solved on a mesh constructed from small area elements. A transformation is introduced through which distorted squares of the physical plane are mapped into computational squares. Two sets of overlapping elements are used; while the thermodynamic properties are calculated at the primary element centres, the flux balance is established on the secondary elements. For transonic flows an artificial compressibility term (upwind density gradient) is added to density in order to produce the desired directional bias in the hyperbolic region. while the entropy does not increase across mass conservative shock jump regions. Comparisons withexperiments and with other numerical and analytical solutions for various turbomachinery configurations show that this approach is comparatively accurate, reliable, and fast.     

An Euler solver for three-dimensional turbomachinery flows

A time-marching finite volume numerical procedure is presented for three-dimensional Euler analysis of turbomachinery flows. The proposed scheme is applied to the conservative form of the Euler equations written in general curvilinear co-ordinates. A simple but computationally efficient grid is constructed. Numerical solution results for three 3D turbine cascade flows have been presented and compared with available measurements as well as with another state-of-the-art 3D Euler analysis numerical solution in order to demonstrate the accuracy and computational efficiency of the analysis method. Also, the predicted results are compared with a 3D potential flow solver and comparison is made with the analytical solution. The proposed method is an accurate and reliable technique for solving the compressible flow equations in turbomachinery geometries.     

A PRESSURE-BASED ALGORITHM FOR HIGH-SPEED TURBOMACHINERY FLOWS

The steady state Navier–Stokes equations are solved in transonic flows using an elliptic formulation. A segregated solution algorithm is established in which the pressure correction equation is utilized to enforce the divergence-free mass flux constraint. The momentum equations are solved in terms of the primitive variables, while the pressure correction field is used to update both the convecting mass flux components and the pressure itself. The velocity components are deduced from the corrected mass fluxes on the basis of an upwind-biased density, which is a mechanism capable of overcoming the ellipticity of the system of equations, in the transonic flow regime. An incomplete LU decomposition is used for the solution of the transport-type equations and a globally minimized residual method resolves the pressure correction equation. Turbulence is resolved through the k–ε model. Dealing with turbomachinery applications, results are presented in two-dimensional compressor and turbine cascades under design and off-design conditions. © 1997 John Wiley & Sons, Ltd.     

2-D/3-D finite-element solution of the steady Euler equations for transonic lifting flow by stream vector correction

In this paper we study the validation of the new formulation (potential-stream vector) of the steady Euler equations in 2-D/3-D transonic lifting regime flow. This approach, which is based on the Helmholtz decomposition of a velocity vector field, is designed to extend the potential approximation of Euler equations for severe situations such as high transonic or rotational subsonic flows. Different results computed by a fixed point algorithm on the stream vector correction are shown and discussed by comparing them with those obtained by the full potential approach.     

Boundary linear integral method for compressible potential flows

A boundary linear integral method based on Green function theory has been developed to solve the full potential equation for subsonic and transonic flows. In this integral method, potential values in the flow region are determined by potential values represented by boundary integrals and a volume integral. The boundary potential values are obtained by implementing the boundary integrals along boundary segments where a linear potential relation is assumed. The volume integral is evaluated in a grid generated by finite element discretization. The volume integral is evaluated only outside the body. Therefore there is no extra boundary treatment required for evaluation of the volume integral. The source term is assumed to be constant in an element integral volume. The volume integral needs to be evaluated only once and can be stored in computer memory for further usage.     

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.     

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.     

A numerical modelling of stator–rotor interaction in a turbine stage with oscillating blades

In real flows unsteady phenomena connected with the circumferential non-uniformity of the main flow and those caused by oscillations of blades are observed only jointly. An understanding of the physics of the mutual interaction between gas flow and oscillating blades and the development of predictive capabilities are essential for improved overall efficiency, durability and reliability. In the study presented, the algorithm proposed involves the coupled solution of 3D unsteady flow through a turbine stage and the dynamics problem for rotor-blade motion by the action of aerodynamic forces, without separating the outer and inner flow fluctuations. The partially integrated method involves the solution of the fluid and structural equations separately, but information is exchanged at each time step, so that solution from one domain is used as a boundary condition for the other domain. 3-D transonic gas flow through the stator and rotor blades in relative motion with periodicity on the whole annulus is described by the unsteady Euler conservation equations, which are integrated using the explicit monotonous finite volume difference scheme of Godunov–Kolgan. The structural analysis uses the modal approach and a 3-D finite element model of a blade. The blade motion is assumed to be constituted as a linear combination of the first natural modes of blade oscillations, with the modal coefficients depending on time. A calculation has been done for the last stage of the steam turbine, under design and off-design regimes. The numerical results for unsteady aerodynamic forces due to stator–rotor interaction are compared with results obtained while taking into account blade oscillations. The mutual influence of both outer flow non-uniformity and blade oscillations has been investigated. It is shown that the amplitude-frequency spectrum of blade oscillations contains the high-frequency harmonics, corresponding to the rotor moving past one stator blade pitch, and low-frequency harmonics caused by blade oscillations and flow non-uniformity downstream from the blade row; moreover, the spectrum involves the harmonics which are not multiples of the rotation frequency.     

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

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

用扰动涡方法模拟叶轮机内动静叶的相互作用

发展了一种研究叶轮机内动、静叶间的相互作用的新方法———扰动涡方法，它利用全三维的定常解为基础解，并由此给出非定常扰动场的初始解．为计算叶片对扰动场的响应过程，采用拉格朗日方法追踪扰动涡团的对流流动过程，用确定性涡方法来描述流体的粘性扩散过程．发展了代数湍流模型（Ｂａｌｄｗｉｎ Ｌｏｍａｘ湍流模型）在尾迹中的应用方法，克服了其它数值方法中无法准确捕捉尾迹中心线的运动轨迹，以及计算出的边界层外的湍流涡粘性系数偏大的缺陷．利用该方法计算轴流叶轮机内由于动、静叶间的相互作用而引起的非定常流动过程，与实验的对照表明，模拟结果与实验数据吻合得相当好，从而说明本文发展的方法是可信的，为更直观地描述尾迹等非定常因素的流动及叶轮机内的掺混问题提供依据．     

REDUCED-ORDER MODELS OF UNSTEADY TRANSONIC VISCOUS FLOWS IN TURBOMACHINERY

The proper orthogonal decomposition (POD) technique is applied in the frequency domain to obtain a reduced-order model of the unsteady flow in a transonic turbomachinery cascade of oscillating blades. The flow is described by a inviscid—viscous model, i.e. a full potential equation outer flow model and an integral equation boundary layer model. The nonlinear transonic steady flow is computed first and then the unsteady flow is determined by a small perturbation linearization about the nonlinear steady solution. Solutions are determined for a full range of frequencies and validated. The full model results and the POD method are used to construct a reduced-order model in the frequency domain. A cascade of airfoils forming the Tenth Standard Configuration is investigated to show that the reduced-order model with only 15–75 degrees of freedom accurately predicts the unsteady response of the full system with approximately 15 000 degrees of freedom.     

A parallel pressure implicit splitting of operators algorithm applied to flows at all speeds

A parallel implementation of the pressure‐based implicit splitting of operators (PISO) method is described and applied to both compressible and incompressible flows. The treatment of variables at the interfaces between adjacent blocks is highlighted, and, for compressible flow, a straightforward method for the implicit handling of density is described. Steady state and oscillatory flow through a sudden expansion are considered at low speeds for both two‐ and three‐dimensional geometries. Extension of the incompressible method to compressible flow is assessed for subsonic, transonic and supersonic flow through a two‐dimensional bump. Although good accuracy is achieved in these high‐speed flows, including the automatic capturing of shock waves, the method is deemed unsuitable for simulating steady state high‐speed flows on fine grids due to the requirement of very small time steps. Copyright © 2001 John Wiley & Sons, Ltd.     

Turbulence microscale variation due to interaction with shock wave

A speckle photographic technique is used for visualizing the planar distribution of the refractive deflection angles of light transmitted through a compressible turbulent flow. Double exposure speckle photography is applied for the statistical analysis of such flows. Quantitative visualization of the planar distribution of the refractive deflection angles of the light transmitted through the compressible turbulent flow is performed. Turbulent structures are visible in the patterns of the deflection angles isolines. The 2-D correlation functions of these deflection angles are constructed and analyzed. A mathematical procedure for solving Erbeck-Merzkirch integral transformation is described and the 3-D density correlation functions are evaluated using 2-D correlation functions of deflection angles of the light passed through the turbulence. Received 14 November 1999 / Accepted 3 June 2000     

Numerical simulation of viscous transonic flows over an airfoil

High Reynolds number viscous transonic flow is described based on an interaction of the potential outer flow with the boundary layer and wake. Following the procedure of Lighthill (1958), the solutions in these domains are matched to each other through boundary conditions. The solution to the complete problem is obtained iteratively through successive computations of the flows in the outer and inner domains. Both old and new algorithms are used for the iteration process and subsequent problem solution. Results are given for all the airfoils from the Experimental Data Base for Computer Program Assessment (AGARD-AR-138, 1979). A comparison of these results with experimental data shows the degree of agreement between these unbounded airfoil flow simulations and real transonic flow over the central part of a straight wing.     

Application of a 3-D FEM/FDM Overlapping Scheme to Heat and Mass Transfer and Moving-Boundary Problems

A 3-D FEM/FDM overlapping scheme for viscous, incompressible flow problems is presented that combines the finite element method, which is best suited to analyze flow in any arbitrarily shaped flow geometry, with the finite difference method, which is advantageous in both computing time and computer storage. The combination of both methods enables large-scale viscous flow to be analyzed, which is crucial both for detailed analysis of 3-D flows and for solving flow problems around moving bodies, A modified ABM AC method is used as the basic algorithm, to which a sophisticated time integration scheme, proposed by the present authors, has been applied. In this paper, some numerical results including 3-D heat and mass transfer problem and moving-boundary problems are presented.     

透平叶栅跨音速流动计算中的新型有限元法

将Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法和多级有限元的思想结合起来，构成了在收敛速度和稳定性两方面均较好的新型有限元算法：多级广义有限元。利用这一方法，分别基于Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和Ｅｕｌｅｒ方程，研究了透平跨音速叶栅无粘流动和粘性流动，并将计算结果与实验结果作了比较。计算结果表明，本方法是透平机械内部跨音速流动计算的强有力的手段。     

Adaptive grids in numerical fluid dynamics

A solution adaptive grid (SAG) method which redistributes the nodal points of a function according to its curvature is presented. A single, user-selected step parameter, P, is available for controlling the maximum step size, allowing the application of the technique to a wide variety of problems. Three test cases are cited: (1) the 1-dimensional inviscid Burgers equation, (2) the Falkner-Skan equation and (3) the finite-volume form of the Navier-Stokes equations for transonic aerofoil flows. In all three cases, significant solution improvement in terms of accuracy and convergence acceleration were achieved.     

Transonic viscous-inviscid interaction by a finite element method

A method is outlined for solving two-dimensional transonic viscous flow problems, in which the velocity vector is split into the gradient of a potential and a rotational component. The approach takes advantage of the fact that for high-Reynolds-number flows the viscous terms of the Navier-Stokes equations are important only in a thin shear layer and therefore solution of the full equations may not be needed everywhere. Most of the flow can be considered inviscid and, neglecting the entropy and vorticity effects, a potential model is a good approximation in the flow core. The rotational part of the flow can then be calculated by solution of the potential, streamfunction and vorticity transport equations. Implementation of the no-slip and no-penetration boundary conditions at the walls provides a simple mechanism for the interaction between the viscous and inviscid solutions and no extra coupling procedures are needed. Results are presented for turbulent transonic internal choked flows.     

One-dimensional turbulence modeling for cylindrical and spherical flows: model formulation and application

The one-dimensional turbulence (ODT) model resolves a full range of time and length scales and is computationally efficient. ODT has been applied to a wide range of complex multi-scale flows, such as turbulent combustion. Previous ODT comparisons to experimental data have focused mainly on planar flows. Applications to cylindrical flows, such as round jets, have been based on rough analogies, e.g., by exploiting the fortuitous consistency of the similarity scalings of temporally developing planar jets and spatially developing round jets. To obtain a more systematic treatment, a new formulation of the ODT model in cylindrical and spherical coordinates is presented here. The model is written in terms of a geometric factor so that planar, cylindrical, and spherical configurations are represented in the same way. Temporal and spatial versions of the model are presented. A Lagrangian finite-volume implementation is used with a dynamically adaptive mesh. The adaptive mesh facilitates the implementation of cylindrical and spherical versions of the triplet map, which is used to model turbulent advection (eddy events) in the one-dimensional flow coordinate. In cylindrical and spherical coordinates, geometric stretching of the three triplet map images occurs due to the radial dependence of volume, with the stretching being strongest near the centerline. Two triplet map variants, TMA and TMB, are presented. In TMA, the three map images have the same volume, but different radial segment lengths. In TMB, the three map images have the same radial segment lengths, but different segment volumes. Cylindrical results are presented for temporal pipe flow, a spatial nonreacting jet, and a spatial nonreacting jet flame. These results compare very well to direct numerical simulation for the pipe flow, and to experimental data for the jets. The nonreacting jet treatment overpredicts velocity fluctuations near the centerline, due to the geometric stretching of the triplet maps and its effect on the eddy event rate distribution. TMB performs better than TMA. A hybrid planar-TMB (PTMB) approach is also presented, which further improves the results. TMA, TMB, and PTMB are nearly identical in the pipe flow where the key dynamics occur near the wall away from the centerline. The jet flame illustrates effects of variable density and viscosity, including dilatational effects.     

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.