高速转动圆盘流动数值模拟研究

重点研究高速离心压气机叶轮与机匣间的间隙流动及其温度分布。研究将离心压气机简化为高速转动圆盘,搭建了相关实验平台,并开展了相应的数值模拟研究。通过改变转动圆盘的转速和轴向进入的冷却流的流量,研究了转速和流量对于间隙内温度和速度分布的影响。结果显示,转速是影响温度变化的最主要因素,转速越大,温度越高;同等幅度的流量变化对温度的影响则较小。研究发现,在实验和模拟对应的大雷诺数条件下,无量纲的速度分布基本不受到圆盘转速、冷却流量和温度场的影响。     

用于高压压气机叶栅设计的影响矩阵法

在表面涡模型的基础上,发展了一种能用于压气机不同叶型,在亚临界工况下叶片反问题设计方法——影响矩阵法。方法利用隐式的残值修正思想构造出一个能够连接叶片表面压力函数变化和叶型坐标变化的影响矩阵,同正问题的二维、无粘、势流的有限元法相结合而成。给出的数值结果表明了影响矩阵法用于叶片设计的有效性和准确性,能够满足工程技术要求。本文发展的这一方法的主要目的是:用于对文献[6]以解析方法得到的相似变换准则所具有的相似精度作出数值修正,从而提高叶片表面附近的流动相似精度。     

Non-intrusive aerodynamic loads analysis of an aircraft propeller blade

The flow field in a cross-sectional plane of a scaled Beaver DHC aircraft propeller has been measured by means of a stereoscopic PIV setup. Phase-locked measurements are obtained in a rotational frequency range from 18,900 to 21,000 rpm, at a relative Mach number of 0.6 at ¾ propeller radius. The use of an adapted formulation of the momentum equation in differential form for rotating frame of references, integrated with isentropic relations as boundary conditions, allowed to compute the pressure field around the blade and the surface pressure distribution directly from the velocity data in the compressible regime. The procedure, extended to the computation of the aerodynamic lift and drag coefficients by a momentum contour integral approach, proved to be able to couple the aerodynamical loads to the flow field on the moving propeller blade, comparing favorably with a numerical simulation of the entire scaled model. Results are presented for two propeller rotation speeds and three different yawing angles.     

Turbulence measurements in a transonic two-passage turbine cascade

This paper presents detailed turbulence measurements in a two-dimensional, transonic, double passage turbine cascade. Particle image velocimetry was used to obtain mean velocity and turbulence measurements all around a single turbine blade within about 2 mm of the blade and wall surfaces. The passage walls were designed using an optimization procedure so that the blade surface pressure distribution matches that of the blade in an infinite cascade. The resulting experimental model captures much of the complexity of a real turbine stage (including high streamline curvature, strong accelerations, and shocks) in a passage with a continuous wall shape, allowing for high measurement resolution and well controlled boundary conditions for comparison to CFD. The measurements show that in the inviscid regions of the passage the absolute level of the turbulent fluctuations does not change significantly as the flow accelerates, while the local turbulence intensity drops rapidly as the flow accelerates. These results provide a benchmark data set that can be used to improve turbulence models.     

喷水推进双级轴流泵流体动力性能CFD分析

采用高质量结构化网格对某喷水推进双级轴流泵整个数值模拟系统的各部件进行区域离散,利用CFD方法预报该泵扬程、功率和效率等外特性,分析其内部流动情况。计算结果与实验数据的较好吻合证明了CFD方法可以作为通过试验求取喷水推进双级轴流泵特性曲线的辅助方法。研究中对有、无轴情况下双级轴流泵外特性的预报和比较结果表明,轴对该泵动力性能预报结果影响不大;通过计算叶片压力分布,验证了该双级泵叶轮叶片设计中采用的线型;通过观察一、二级导叶出口的流动状态,分析了该泵设计中存在的缺陷,并为导叶改进设计提出建议。     

Nonlinear vibrations of blade with varying rotating speed

This paper investigates the nonlinear dynamic responses of the rotating blade with varying rotating speed under high-temperature supersonic gas flow. The varying rotating speed and centrifugal force are considered during the establishment of the analytical model of the rotating blade. The aerodynamic load is determined using first-order piston theory. The rotating blade is treated as a pretwist, presetting, thin-walled rotating cantilever beam. Using the isotropic constitutive law and Hamilton??s principle, the nonlinear partial differential governing equation of motion is derived for the pretwist, presetting, thin-walled rotating beam. Based on the obtained governing equation of motion, Galerkin??s approach is applied to obtain a two-degree-of-freedom nonlinear system. From the resulting ordinary equation, the method of multiple scales is exploited to derive the four-dimensional averaged equation for the case of 1:1 internal resonance and primary resonance. Numerical simulations are performed to study the nonlinear dynamic response of the rotating blade. In summary, numerical studies suggest that periodic motions and chaotic motions exist in the nonlinear vibrations of the rotating blade with varying speed.     

Investigation of flow separation in a transonic-fan linear cascade using visualization methods

An extensive experimental study into the nature of the separated flows on the blade suction surface of modern transonic fans is described in this paper. The study was a subtask of a larger experimental effort focused on blade flutter excited by flow separation in the blade tip region. The tip sections of airfoils on transonic fan blades are designed for precompression and consequently they differ from sections on the rest of the blade. The blade tip section was modeled by a low aspect ratio blade and therefore most of the blade tested was exposed to the secondary flow effects. The aim of this work was to supply reliable data on flow separation on transonic fan blades for validation of future analytical studies. The experimental study focused on two visualization techniques: surface flow visualization using dye oils and schlieren (and shadowgraph) flow visualization. The following key observations were made during the study. For subsonic inlet flow, the flow on the suction surface of the blade was separated over a large portion of the blade, and the separated area increased with increasing inlet Mach number. For the supersonic inlet flow condition, the flow was attached from the leading edge up to the point where a bow shock from the upper neighboring blade imposed on the blade surface. Downstream, there was a separated flow region in which air flowed in the direction opposite the inlet flow. Finally, past the separated flow region, the flow reattached to the blade surface. For subsonic inlet flow, the low cascade solidity resulted in an increased area of separated flow. For supersonic flow conditions, the low solidity resulted in an improvement in flow over the suction surface.     

考虑法向压力梯度的三维旋转边界层计算

本文以量级分析为基础,建立了一般曲线坐标系上的三维旋转边界层方程。对旋转在边界层中的影响进行分析之后,提出了一个能够处理壁面法向压力梯度不为零问题的压力梯度迭代方法。在传统的Box法的基础上发展了一套完整的求解三维旋转边界层的方法和程序,并对螺旋面、压气机转子叶面以及圆柱面上的旋转边界层进行了计算,与他人的计算和实验的对比分析表明,该方法和程序是正确的,可用于求解任意几何物面上的三维旋转边界层。     

Flow control of hub corner stall in a highly loaded axial compressor cascade

Three-dimensional corner stall is one of the key factors limiting the compressor performance. This paper presents a detailed experimental and computational study of a flow control strategy involving the endwall suction, aiming to eliminate the hub corner stall in a highly loaded axial compressor cascade. Various mass flow suction cases were parametrically tested with the aim of eliminating the corner stall by applying a minimum suction flow ratio. In the experiments, seven-hole pressure probe traverses, different loading distributions and surface oil flow visualizations were applied to address the flow and loss mechanisms in the cascade. Measurements were supplemented with numerical predictions from a commercially available CFD code. It was found that the corner stall, characterized by a large amount of reversed fluid, occupied a large region over the blade suction surface in the highly loaded compressor airfoil, rather than occurring at the junction of a blade suction surface and the endwall as in the conventionally loaded compressor airfoil. By applying flow control, the dominant flow structures, e.g. the flow separations and particularly the corner stall, within the compressor cascade were significantly affected. The maximum spanwise penetration depth of the endwall flow on the suction surface was significantly decreased once the endwall suction flow was applied. Furthermore, the corner stall was completely eliminated by suctioning the mass flow at a specific ratio of the inlet boundary layer flow rate. The midspan flow field was not notably affected, and a further increase in suction mass flow did not benefit the flow field approaching the endwall.     

Numerical simulation of blade coating with short dwell and roll application coaters

Paper and board are often coated at high speeds with a mineral-based aqueous suspension in order to improve their printing properties. This suspension is usually called coating colour. The flow behaviour of the coating colour in the cavity of the short dwell coater (SDC) and in the vicinity of the blade tip when paper is coated with a stiff blade has been analysed using the finite element method. The models used to simulate the flow incorporated free surfaces and shear-thinning colours. The Newtonian case was in some cases also included in the modelling. The viscosity level and the shear-thinning character of the coating colour had a significant influence on the flow in the SDC cavity, although the overall behaviour was to a large extent governed by the speed of the coater. The pressure distribution along the paper surface in the SDC cavity was also analysed.In agreement with earlier reported results, increasing the machine speed raised the pressure level in the colour close to the blade tip. The rheological properties of the coating colour also affected the flow field and the pressure distribution in this region; e.g. it was found here that a colour with a high viscosity level at low shear rates developed a high pressure level close to the blade tip. The most interesting result revealed by the analysis was that changes in the configuration close to the blade tip (converging flow between the blade tip and the paper and compression of the substrate under the blade tip) and boundary conditions at the blade tip surface had a very significant effect on the pressure distribution. The predictions of the numerical simulation were to some extent compared with experience from practical coating trials.     

Progresses in application of computational fluid dynamic methods to large scale wind turbine aerodynamics

The computational fluid dynamics (CFD) methods are applied to aerodynamic problems for large scale wind turbines. The progresses including the aerodynamic analyses of wind turbine profiles, numerical flow simulation of wind turbine blades, evaluation of aerodynamic performance, and multi-objective blade optimization are discussed. Based on the CFD methods, significant improvements are obtained to predict two/three dimensional aerodynamic characteristics of wind turbine airfoils and blades, and the vortical structure in their wake flows is accurately captured. Combining with a multi-objective genetic algorithm, a 1.5 MW NH-1500 optimized blade is designed with high efficiency in wind energy conversion.     

Investigation of horizontal-axis wind turbine (HAWT) blade threedimensional rotational effect based on field experiments

Field experiments are performed on a two-bladed 33 kW horizontal-axis wind turbine (HAWT). The pressures are measured with 191 pressure sensors positioned around the surfaces of seven spanwise section airfoils on one of the two blades. Three-dimensional (3D) and two-dimensional (2D) numerical simulations are performed, respectively, on the rotor and the seven airfoils of the blade. The results are compared with the experimental results of the pressure distribution on the seven airfoils and the lift coefficients. The 3D rotational effect on the blade aerodynamic characteristics is then studied with a numerical approach. Finally, some conclusions are drawn as follows. From the tip to the root of the blade, the experimental differential pressure of the blade section airfoil increases at first and then decreases gradually. The calculated 3D result of the pressure distribution on the blade surface is closer to that of the experiment than the 2D result. The 3D rotational effect has a significant impact on the blade surface flow and the aerodynamic load, leading to an increase of the differential pressure on the airfoils and their lift coefficient than that with the 2D one because of the stall delay. The influence of the 3D rotational effect on the wind turbine blade especially takes place on the sections with flow separation.     

Flow field studies on a micro-air-vehicle-scale cycloidal rotor in forward flight

This paper examines the flow physics and principles of force production on a cycloidal rotor (cyclorotor) in forward flight. The cyclorotor considered here consists of two blades rotating about a horizontal axis, with cyclic pitch angle variation about the blade quarter-chord. The flow field at the rotor mid-span is analyzed using smoke flow visualization and particle image velocimeV are compared with flow fields predicted using 2D CFD and time-averaged force measurements acquired in an open-jet wind tunnel at three advance ratios. It is shown that the experimental flow field is nearly two dimensional at μ = 0.73 allowing for qualitative comparisons to be made with CFD. The incoming flow velocity decreases in magnitude as the flow passes through the retreating (upper) half of the rotor and is attributed to power extraction by the blades. A significant increase in flow velocity is observed across the advancing (lower) half of the rotor. The aerodynamic analysis demonstrates that the blades accelerate the flow through the lower aft region of the rotor, where they operate in a high dynamic pressure environment. This is consistent with CFD-predicted values of instantaneous aerodynamic forces which reveal that the aft section of the rotor is the primary region of force production. Phase-averaged flow field measurements showed two blade wakes in the flow, formed by each of the two blades. Analysis of the blades at several azimuthal positions revealed two significant blade-wake interactions. The locations of these blade-wake interactions are correlated with force peaks in the CFD-predicted instantaneous blade forces and highlight their importance to the generation of lift and propulsive force of the cyclorotor.     

Implementation of all‐Mach Roe‐type schemes in fully implicit CFD solvers – demonstration for wind turbine flows

This paper presents the implementation of all‐Mach Roe‐type schemes in a fully implicit CFD solver. Simple 2D cases, such as the flow around inviscid and viscous aerofoils, were used for an initial validation of these methods, along with more challenging computations consisting of the 3D flow around the Model Experiments in Controlled Conditions wind turbine, in parked and rotating conditions. This work is motivated by the increased interest of the wind turbine industry in larger diameter wind turbines where compressibility effects near the blade tips may be important. Instead of using an incompressible flow solver, this paper explores the option of modifying an existing, efficient, compressible flow solver for use at lower Mach numbers. The good performance of the Roe solver and its popularity influenced the selection of schemes for this work. The results suggest that effective all‐Mach solutions are possible with implicit solvers, and the paper defines the implementation of the new fluxes and Jacobian, including an investigation of some numerical parameters, using as platform the Helicopter Multi‐Block solver of Liverpool University. Copyright © 2013 John Wiley & Sons, Ltd.     

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

Experimental evaluation of numerical simulation of cavitating flow around hydrofoil

Cavitation in hydraulic machines causes different problems that can be related to its unsteady nature. An experimental and numerical study of developed cavitating flow was performed. Until now simulations of cavitating flow were limited to the self developed “in house” CFD codes. The goal of the work was to experimentally evaluate the capabilities of a commercial CFD code (Fluent) for simulation of a developed cavitating flow. Two simple hydrofoils that feature some 3D effects of cavitation were used for the experiments. A relatively new technique where PIV method combined with LIF technique was used to experimentally determine the instantaneous and average velocity and void ratio fields (cavity shapes) around the hydrofoils. Distribution of static pressure on the hydrofoil surface was determined. For the numerical simulation of cavitating flow a bubble dynamics cavitation model was used to describe the generation and evaporation of vapour phase. An unsteady RANS 3D simulation was performed. Comparison between numerical and experimental results shows good correlation. The distribution and size of vapour structures and the velocity fields agree well. The distribution of pressure on the hydrofoil surface is correctly predicted. The numerically predicted shedding frequencies are in fair agreement with the experimental data.     

过冷水滴撞击螺旋桨桨叶表面的数值模拟

本文介绍了某型螺旋桨在不同飞行状态下.求解桨叶水滴撞击特性的数值方法.该方法对桨叶运动模型进行简化,并在对绕桨叶运动的气流场计算的基础上,采用拉格朗日方法求解气流场中水滴运动方程,得到水滴运动轨迹.进而,确定了水滴对桨叶的撞击特性参数,为桨叶防冰系统设计提供条件.主要结论如下:(1)在巡航状态下,桨叶沿展向方向上总收集系数Em和局部收集系数β不断增大;(2)在爬升状态下,随着爬升高度H不断增大,飞行速度V0不断增大,水滴撞击在桨叶表面的范围有所增加,而且β随之增大;(3)随着水滴平均有效直径(MVD)的增大,水滴撞击在桨叶表面的范围明显增加同时,β在桨叶表面同一位置的值也随之增大.     

Numerical Simulation of the 3D Laminar Viscous Flow on a Horizontal-axis Wind Turbine Blade

The aim of this paper is to describe the methodology followed in order to determine the viscous effects of a uniform wind on the blades of small horizontal-axis wind turbines that rotate at a constant angular speed. The numerical calculation of the development of the three-dimensional boundary layer on the surface of the blades is carried out under laminar conditions and considering flow rotation, airfoil curvature and blade twist effects. The adopted geometry for the twisted blades is given by cambered thin blade sections conformed by circular are airfoils with constant chords. The blade is working under stationary conditions at a given tip speed ratio, so that an extensive laminar boundary layer without flow separation is expected. The boundary layer growth is determined on a non-orthogonal curvilinear coordinate system related to the geometry of the blade surface. Since the thickness of the boundary layer grows from the leading edge of the blade and also from the tip to the blade root, a domain transformation is proposed in order to solve the discretized equations in a regular computational 3D domain. The non-linear system of partial differential coupled equations that governs the boundary layer development is numerically solved applying a finite difference technique using the Krause zig-zag scheme. The resulting coupled equations of motion are linearized, leading to a tridiagonal system of equations that is iteratively solved for the velocity components inside the viscous layer applying the Thomas algorithm, procedure that allows the subsequent numerical determination of the shear stress distribution on the blade surface.