Modelling of Unsteady Transition in Low-Pressure Turbine Blade Flows with Two Dynamic Intermittency Equations

A transition model for describing wake-induced transition is presented based on the SST turbulence model by Menter and two dynamic equations for intermittency: one for near-wall intermittency and one for free-stream intermittency. In the Navier-Stokes equations, the total intermittency factor, which is the sum of the two, multiplies the turbulent viscosity computed by the turbulence model. The quality of the transition model is illustrated on the T106A test cascade for two Reynolds numbers, using experimental results by Stieger and Hodson for transition mainly due to kinematic wake impact on a separation bubble. The quality of the model is also revealed on the T106D test cascade using experimental results from Hilgenfeld, Stadtmuller and Fottner for wake turbulence induced transition. The test cases differ in pitch to chord ratio, Reynolds number and inlet free-stream turbulence intensity, causing different transition mechanisms. The unsteady results are presented in space-time diagrams of shape factor and wall shear stress on the suction side. The results show the capability of the model to capture the physics of unsteady transition in separated state. Inevitable shortcomings are revealed as well.     

Experimental study on the effect of unsteadiness on boundary layer development on a linear turbine cascade

 The results from an experimental investigation of unsteady boundary layer behavior on a linear turbine cascade are presented in this paper. To perform a detailed study on unsteady cascade aerodynamics and heat transfer, a new large-scale, high-subsonic research facility for simulating the periodic unsteady flow has been developed. It is capable of sequentially generating up to four different unsteady inlet flow conditions that lead to four different passing frequencies, wake structures, and freestream turbulence intensities. For a given Reynolds number, two different unsteady wake formations are utilized. Detailed unsteady boundary layer velocity. turbulence intensity, and pressure measurements are performed along the suction and pressure surfaces of one blade. The results display the transition and development of the boundary layer, ensemble-averaged velocity, and turbulence intensity. Received: 23 September 1996/Accepted: 19 February 1997     

Flow visualization and LDV measurements of laminar flow in a helical square duct with finite pitch

Flow visualization and LDV measurements are performed on laminar flow in a helical square duct with finite pitch. The experimental observations are compared to results of numerical calculations employing the finite-volume method and assuming a fully developed flow. Good agreement is found between measured and computed velocity profiles. This suggests that the physical velocity components used in the numerical calculations are suitable to describe the flow. It is further demonstrated that the contravariant velocity components employed by some authors may lead to results that are difficult to interpret. Two stable solution branches are detected in the numerical calculations. For Reynolds numbers between the stable branches, unsteady and fully developed computations predict an oscillating flow between a two-vortex and a four-vortex structure. In the experiments, the flow normally retained a stable two-vortex structure in the numerically predicted unstable regime. However, by disturbing the flow at the duct inlet, a four-vortex flow that showed similarities to the computed flow could occasionally be obtained. For Reynolds numbers above 600, unsteady flow behavior was observed both experimentally and numerically, which might be an early sign of transition. In the experiments, Gdrtler-like extra vortices emerged spontaneously from the outer wall without disturbing the flow at the inlet. The same phenomenon was observed in the numerical calculations, assuming an unsteady and fully developed flow, but the extra vortices appeared with a lower frequency than in the experiments.     

Large-Eddy Simulation of Transition to Turbulence in Boundary Layers

Large-eddy simulation (LES) results for laminar-to-turbulent transition in a spatially developing boundary layer are presented. The disturbances are ingested into a laminar flow through an unsteady suction-and-blowing strip. The filtered, three-dimensional time-dependent Navier–Stokes equations are integrated numerically using spectral, high-order finite-differences, and a three-stage low-storage Runge–Kutta/Crank–Nicolson time-advancement method. The buffer-domain technique is used for the outflow boundary condition. The localized dynamic model used to parametrize the subgrid-scale (SGS) stresses begins to have a significant impact at the beginning of the nonlinear transition (or intermittency) region. The flow structures commonly found in experiments are also observed in the present simulation; the computed linear instability modes and secondary instability $\Lambda$-vortex structures are in agreement with the experiments, and the streak-like structures and turbulent statistics compare with both the experiments and the theory. The physics captured in the present LES are consistent with the experiments and the full Navier–Stokes simulation (DNS), at a significant fraction of the DNS cost. A comparison of the results obtained with several SGS models shows that the localized model gives accurate results both in a statistical sense and in terms of predicting the dynamics of the energy-carrying eddies, while requiring fewer ad hoc adjustments than the other models. Received: 5 April 1996 and accepted 27 March     

低压涡轮内部流动及其气动设计研究进展

随着高空无人飞行器研究的不断升温, 高空低雷诺数条件下动力装置的研究越来越受到人们的重视.结合近年来国内外相关领域的研究工作, 对低雷诺数低压涡轮内部复杂流动机理的研究进展进行了介绍, 包括低雷诺数情况下低压涡轮内部非定常流动的特点, 叶片边界层分离及转捩现象机理, 上游周期性尾迹与下游叶片边界层相互作用机理等. 在此基础上给出了适合低雷诺数条件的低压涡轮气动设计方法:尾迹通过与边界层的相互作用, 能够抑制分离, 进而减小叶型损失, 在气动设计中有效引入非定常效应可以大幅度提高低压涡轮的气动负荷或降低气动损失, 最终达到提高性能的目的;数值及实验结果验证了这种设计方法的有效性.      

A general approach for computing unsteady 3D thin lifting and/or propulsive systems derived from a complete theory

This paper presents the basis of a numerical method for unsteady aerodynamic computation around thin lifting and/or propulsive systems with arbitrary variable geometries, involving the velocity field, the velocity potential, the pressure field and the wake characteristics (geometry and vortex strength). Most of the corresponding theory actually stems from the unsteady wake model established by Mudry, in which the wake is considered to be a median layer, characterized by a pair of functions on which Mudry founded the concept of continuous vortex particle. The governing relations of the continuous problem are then the flow tangency condition, the wake integro‐differential evolution equation, and a flow regularity condition at the trailing edge. This constitutes a rigorous and complete theoretical formulation of this problem, from which a discretization scheme and a numerical method of solution are derived. The view of the vortex wake is similar to the one in the classical vortex lattice approaches, but uses a discrete vortex particle concept, particularly well suited for the prediction of the unsteady wake deformation. This, together with the continuous theory, ensures the computing method compares favorably with the classical methods in terms of flexibility and computing costs. In order to demonstrate the capabilities of the present method, the calculation of flapping wings of variable geometry is also presented. Copyright © 1999 John Wiley & Sons, Ltd.     

Experimental investigation about the effect of non-axisymmetric wake impact on a low speed axial compressor

Non-axisymmetric wake impact experiments were carried out after the best exciting frequency for a low speed axial compressor had been found by axisymmetric wake impact experiments. When the number and circumferential distribution of inlet guide vanes （IGV） are logical, the wakes of non-axisymmetric IGVs can exert beneficial unsteady exciting effect on their downstream rotor flow fields and improve the compressor＇s performance. In the present paper, four non-axisymmetric wake impact plans were found working better than the axisymmetric wake impact plan. Compared with the base plan, the best non-axisymmetric plan increased the compressor＇s peak efficiency, and the total pressure rise by 1.1 and 2%, and enhanced the stall margin by 4.4%. The main reason why non-axisymmetric plans worked better than the axisymmetric plan was explained as the change of the unsteady exciting signal arising from IGV wakes. Besides the high-frequency components, the nonaxisymmetric plan generated a beneficial low-frequency square-wave exciting signal and other secondary frequency components. Compared with the axisymmetric plan, multifrequency exciting wakes arising from the non-axisymmetric plans are easier to get coupling relation with complex vortices such as clearance vortices, passage vortices and shedding vortices.     

DES evaluation of near-wake characteristics in a shallow flow

Three-dimensional numerical modeling using Detached Eddy Simulation (DES) based on unsteady Reynolds-Averaged Navier–Stokes (RANS) with the k–ω SST (Shear-Stress Transport) turbulence model has been carried out to evaluate the characteristics of a shallow wake flow. The shallow wake is generated by inserting a sharp-edged bluff body in the open channel flow. A horseshoe vortex is captured in front of the body, which stretches downstream and envelops the vortices that form part of the shear layers. The mean and instantaneous flow field characteristics in the wake are examined and compared at different downstream locations to evaluate the three-dimensional features in the flow. Streamwise positive directed velocity is observed in the wake centerline at horizontal planes close to the bed. Flow features hitherto not captured in experimental studies can be identified in sections parallel to the bed and body. A typical signature of three-dimensionality, upward ejection of fluid elements from the bed towards the free surface, is also observed in the wake.     

串列双圆柱绕流问题的数值模拟

本文运用有限体积方法,对绕串列放置的双圆柱的二维不可压缩流动进行了数值计算。为研究两圆柱不同间距对圆柱相互作用和尾流特征的影响,选取间距比Ｌ／Ｄ（Ｌ为两圆柱中心间的距离,Ｄ为圆柱直径）在１．５～５．０之间每隔０．５共八个有代表性的间距进行了计算模拟。计算均在Ｒｅ＝２００条件下进行。计算结果表明：对该绕流问题,流动特征在很大程度上取决于间距的大小。且间距存在一临界值,间距比从小于临界值变化到大于临界     

Transonic Rotor/Stator Interaction Using a 3-D Linearised Unsteady Viscous Flow Model

The paper first summarizes the forced response problem in turbomachinery and reviews various numerical methods for the simulation of unsteady flows. A particular technique, based on the linearisation of the unsteady Favre-averaged Navier-Stokes equations on three-dimensional mixed-element grids of tetrahedra, hexahedra and wedges, is described in some detail. The methodology was applied to a NGV/rotor interaction benchmark case for which detailed steady and unsteady flow measurements are available. The steady-state flow, calculated using a non-linear viscous representation, was described in detail with emphasis on features such as separation, horseshoe and passage vortices, tip leakage and shock structure since these are likely to influence the unsteady flow. The sources of unsteadiness on the rotor passage were evaluated from the steady-state solution at the NGV outlet. The disturbances were split into vortical, entropic and potential waves, the Fourier components of which were considered separately. The summation of the vortical and entropic waves was used as a rotor inlet boundary condition in order to assess the wake/rotor unsteady interaction. Similarly, potential waves were used to study the potential/rotor interaction. The results obtained from these two types of unsteady interactions were superimposed and compared with experimental data. Good qualitative and, in most cases, quantitative agreement was obtained, a finding which suggests that the unsteady flowfield generated by the relative blade motion can be considered to be a quasi-linear phenomenon for the particular HP turbine studied. Finally, the mechanisms of wake/rotor and potential/rotor interactions were studied in some detail and it was concluded that the former was strong in the crown of the blade while the latter was dominant in the leading edge region.     

非结构动网格在三维可动边界问题中的应用

研究用于非结构动网格的弹簧近似方法，采用顶点弹簧描述，导出并讨论了弹簧倔强系数的取值。通过引入边界修正和扭转效应修正，对标准弹簧近似方法进行了改进，转动翼型算例的结果表明，改进后的方法大大提高了网格变形能力和网格质量，应用该动网格方法耦合求解基于（Arbitrary Lagrangian-Eulerian,ALE）描述的三维Euler方程，模拟了作俯抑振动的矩形机翼绕流，计算结果与实验数据及文献计算结果十分一致，作为多个自由刚与流体耦合运动问题的简单例证，耦合刚体动力学方程，模拟了激波与双立方体的相互作用，得到了非定常流场结构，研究表明，基于弹簧近似的非结构动网格与有限体积流式流场解算器相结合，是模拟包含运动边界的非定常流动问题的有效方法。     

Computational analysis of impinging shock-wave boundary layer interaction under conditions of incipient separation

The interaction of an oblique shock wave with a turbulent boundary layer under conditions of incipient separation is analyzed by means of large-eddy simulation (LES) and Reynolds-averaged Navier–Stokes (RANS) turbulence models, with the objective to explore their predictive capabilities, in particular with respect to the unsteady features of the interaction. Consistent with earlier direct numerical simulations, we have found that the flow dynamics in the interaction zone is characterized by strong intermittency associated with the formation of scattered spots of flow reversal near the nominal position of the reflected shock. Comparison with experimental results (at much larger Reynolds number) show that the qualitative features of the interaction are predicted reasonably well by both LES and RANS models. RANS models supplemented with a semi-empirical closure are also found to provide reasonable estimate of the fluctuating pressure loads at the wall.     

脉冲爆震发动机进气道气动性能的数值研究

采用有限体积法计算了脉冲爆震发动机某轴对称超音速进气道在3种 不同出口条件(单个正弦扰动压力、某脉冲爆震发动机爆震室头部表压和进气道出口堵塞) 下的进气道内结尾正激波的运动情况，得出了进气道内结尾正激波运动特性和不同出口条件 的关系. 在计算中，采用了多块结构化网格，控制体积的界面无黏通量采用三阶迎风格 式插值获得，同时采用了minmod通量限制器以确保在激波处的解的物理特性；扩散通量采 用二阶中心差分格式插值获得. 定常计算采用当地时间步法，非定常计算采用双时间步法. 离散的代数方程采用交替方向迭代法求解。     

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

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

Simulation and Modelling of Turbulent Trailing-Edge Flow

Computations of turbulent trailing-edge flow have been carried out at a Reynolds number of 1000 (based on the free-stream quantities and the trailing-edge thickness) using an unsteady 3D Reynolds-Averaged Navier–Stokes (URANS) code, in which two-equation (k–ε) turbulence models with various low-Re near wall treatments were implemented. Results from a direct numerical simulation (DNS) of the same flow are available for comparison and assessment of the turbulence models used in the URANS code. Two-dimensional URANS calculations are carried out with turbulence mean properties from the DNS used at the inlet; the inflow boundary-layer thickness is 6.42 times the trailing-edge thickness, close to typical turbine blade flow applications. Many of the key flow features observed in DNS are also predicted by the modelling; the flow oscillates in a similar way to that found in bluff-body flow with a von Kármán vortex street produced downstream. The recirculation bubble predicted by unsteady RANS has a similar shape to DNS, but with a length only half that of the DNS. It is found that the unsteadiness plays an important role in the near wake, comparable to the modelled turbulence, but that far downstream the modelled turbulence dominates. A spectral analysis applied to the force coefficient in the wall normal direction shows that a Strouhal number based on the trailing-edge thickness is 0.23, approximately twice that observed in DNS. To assess the modelling approximations, an a priori analysis has been applied using DNS data for the key individual terms in the turbulence model equations. A possible refinement to account for pressure transport is discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Bubble size distribution prediction for large-scale ship flows: Model evaluation and numerical issues

Prediction of the bubble size distribution in the wake of a ship is important to analyze its acoustic signature. To achieve CFD simulation of dynamic ships with moving control surfaces and rotating propellers in waves, a robust implementation is paramount. In this work a mass conserving multigroup discretization strategy of the Boltzmann transport equation for polydispersed bubbly flows is presented, as well as an analysis of available breakup and coalescence models. Modifications of the discrete equations for the fixed pivot method at the boundaries are introduced that guarantee exact bubble mass conservation. The role of the time stepping scheme in the conservation of mass and number of bubbles is discussed. Though the conservation properties of the discrete system of equations are satisfied provided they are solved exactly, in practice an iterative procedure must be used since the ODE’s are non-linear. Three iterative schemes are proposed and they are analyzed in terms of robustness and efficiency. Breakup, coalescence and dissolution models are analyzed from the numerical point of view. Available models of breakup and coalescence are studied finding appropriate choices for ship applications. Other models are appropriate as well, but are more costly numerically. As appropriate for ship applications, an extension to the model of Prince and Blanch for salt water is proposed and analyzed. The final model is tested against experimental data and computations by other researchers, and convergence properties in bubble size discretization is studied. It is found that for salt water the final steady state is dependent on the initial condition since there is a range of sizes for which coalescence and breakup are both negligible.     

A three-dimensional simulation of a steady approach flow past a circular cylinder at low Reynolds number

The three-dimensional (3D) unsteady viscous wake of a circular cylinder exposed to a steady approach flow is calculated using a fractional-step finite-difference/spectral-element method. The calculated flow fields at Reynolds numbers of 100 (2D) and 200 (3D) are examined in detail. The flow field at Re = 100 is 2D as expected, while the flow field at Re = 200 has distinct 3D features, with spanwise wavelengths of about 3.75 cylinder diameters. The calculated results produce drag and lift coefficients and Strouhal numbers that agree extremely well with the experimental values. These 3D values at Re = 200 are in better agreement with experimental values than the results of a 2D calculation at Re = 200, which is expected. © 1998 John Wiley & Sons, Ltd.     

Unsteady response of airfoils due to small-scale pitching motion with considerations for foil thickness and wake motion

Unsteady pressures, forces, and pitching moments generated by foils experiencing vibratory motion in an incompressible, attached flow configuration are studied within this work. Specifically, two-dimensional, unsteady potential flow and unsteady Reynolds-Averaged Navier–Stokes calculations are performed on various Joukowski foils undergoing sinusoidal, variable amplitude, small-scale pitching motion at a chord-based Reynolds number of 10${}^6$ over a range of reduced frequencies between 0.01–100. These calculated results from both approaches are compared directly to predictions from implementing the Theodorsen model, which treats foils as infinitely thin, flat plates that shed a planar sheet of vorticity. The effects of relaxing these seemingly strict conditions are explored, and the particular terms which control the unsteady responses are identified and discussed. For increasing pitch amplitudes and reduced frequencies the shed wake is seen to become quite non-planar and to form coherent vortex structures. Despite this wake behavior, the normalized airfoil responses at the disturbance reduced frequency are seen to be largely unaffected. However, non-negligible responses are generated across a wide range of other frequencies. Potential flow calculations for symmetric Joukowski foils show that there is marginal effect of foil thickness at reduced frequencies less than one. For higher reduced frequency conditions however, the unsteady lift response is seen to experience both an amplification of level and a phase shift relative to the Theodorsen model. A specific augmenting expression is developed for this behavior through analysis within the potential flow framework.     

Comparison of numerical solvers for cavitating flows

Different computational fluid dynamics (CFD) strategies have been developed to simulate, analyse and better understand cavitating flows. Based on homogeneous models, two numerical approaches using compressible and incompressible codes are applied to capture large density variations and unsteady behaviours of cavitating flows. Simulations are performed on two-dimensional Venturi geometries and compared with experimental data. Local and global analyses are proposed and the necessity to account for compressibility phenomena is discussed.