Numerical study of compression corner flowfield using Gao-Yong turbulence model

A numerical simulation of shock wave turbulent boundary layer interaction induced by a 24° compression corner based on Gao-Yong compressible turbulence model was presented.The convection terms and the diffusion terms were calculated using the second-order AUSM (advection upstream splitting method) scheme and the second-order central difference scheme,respectively.The Runge-Kutta time marching method was employed to solve the governing equations for steady state solutions.Significant flow separation-region which indicates highly non-isotropic turbulence structure has been found in the present work due to intensity interaction under the 24° compression corner.Comparisons between the calculated results and experimental data have been carried out,including surface pressure distribution,boundary-layer static pressure profiles and mean velocity profiles.The numerical results agree well with the experimental values,which indicate Gao-Yong compressible turbulence model is suitable for the prediction of shock wave turbulent boundary layer interaction in two-dimensional compression corner flows.      

Experimental Study of Corner Stall in a Linear Compressor Cascade

In order to gain a better knowledge of the mechanisms and to calibrate computational fluid dynamics (CFD) tools including both Reynolds-averaged Navier-Stokes (RANS) and large eddy simulation (LES),a detailed and accurate experimental study of corner stall in a linear compressor cascade has been carried out.Data are taken at a Reynolds number of 382 000 based on blade chord and inlet velocity.At first,inlet flow boundary layer is surveyed using hot-wire anemometry.Then in order to investigate the effects of incidence,measurements are acquired at five incidences,including static pressures on both blade and endwall surfaces measured by pressure taps and the total pressure losses of outlet flow measured by a five-hole pressure probe.The maximum losses as well as the extent of losses of the corner stall are presented as a function of the investigated incidences.     

Experimental study of flow field distribution over a generic cranked double delta wing

The flow fields over a generic cranked double delta wing were investigated. Pressure and velocity distributions were obtained using a Pitot tube and a hot wire anemometer. Two different leading edge shapes, namely ‘‘sharp" and ‘‘round", were applied to the wing. The wing had two sweep angles of 55° and 30°. The experiments were conducted in a closed circuit wind tunnel at velocity 20 m/s and angles of attack of 5°–20° with the step of 5°. The Reynolds number of the model was about 2 ×10~5 according to the root chord. A dual vortex structure was formed above the wing surface. A pressure drop occurred at the vortex core and the root mean square of the measured velocity increased at the core of the vortices, reflecting the instability of the flow in that region. The magnitude of power spectral density increased strongly in spanwise direction and had the maximum value at the vortex core. By increasing the angle of attack, the pressure drop increased and the vortices became wider; the vortices moved inboard along the wing, and away from the surface; the flow separation was initiated from the outer portion of the wing and developed to its inner part. The vortices of the wing of the sharp leading edge were stronger than those of the round one.     

Effect of blade sweep on inlet flow in axial compressor cascades

This paper presents comparative numerical studies to investigate the effects of blade sweep on inlet flow in axial compressor cascades. A series of swept and straight cascades was modeled in order to obtain a general understanding of the inlet flow field that is induced by sweep.A computational fluid dynamics(CFD) package was used to simulate the cascades and obtain the required three-dimensional(3D) flow parameters. A circumferentially averaged method was introduced which provided the circumferential fluctuation(CF) terms in the momentum equation.A program for data reduction was conducted to obtain a circumferentially averaged flow field.The influences of the inlet flow fields of the cascades were studied and spanwise distributions of each term in the momentum equation were analyzed. The results indicate that blade sweep does affect inlet radial equilibrium. The characteristic of radial fluid transfer is changed and thus influencing the axial velocity distributions. The inlet flow field varies mainly due to the combined effect of the radial pressure gradient and the CF component. The axial velocity varies consistently with the incidence variation induced by the sweep, as observed in the previous literature. In addition, factors that might influence the radial equilibrium such as blade camber angles, solidity and the effect of the distance from the leading edge are also taken into consideration and comparatively analyzed.     

Influence of coupled boundary layer suction and bowed blade on flow field and performance of a diffusion cascade

Based on the investigation of mid-span local boundary layer suction and positive bowed cascade, a coupled local tailored boundary layer suction and positive bowed blade method is developed to improve the performance of a highly loaded diffusion cascade with less suction slot. The effectiveness of the coupled method under different inlet boundary layers is also investigated.Results show that mid-span local boundary layer suction can effectively remove trailing edge separation, but deteriorate the flow fields near the endwall. The positive bowed cascade is beneficial for reducing open corner separation, but is detrimental to mid-span flow fields. The coupled method can further improve the performance and flow field of the cascade. The mid-span trailing edge separation and open corner separation are eliminated. Compared with linear cascade with suction, the coupled method reduces overall loss of the cascade by 31.4% at most. The mid-span loss of the cascade decreases as the suction coefficient increases, but increases as bow angle increases. The endwall loss increases as the suction coefficient increases. By contrast, the endwall loss decreases significantly as the bow angle increases. The endwall loss of coupled controlled cascade is higher than that of bowed cascade with the same bow angle because of the spanwise inverse ‘‘C" shaped static pressure distribution. Under different inlet boundary layer conditions, the coupled method can also improve the cascade effectively.     

Alleviation of spike stall in axial compressors utilizing grooved casing treatment

This article deals with application of grooved type casing treatment for suppression of spike stall in an isolated axial compressor rotor blade row. The continuous grooved casing treatment covering the whole compressor circumference is of 1.8 mm in depth and located between90% and 108% chord of the blade tip as measured from leading edge. The method of investigation is based on time-accurate three-dimensional full annulus numerical simulations for cases with and without casing treatment. Discretization of the Navier–Stokes equations has been carried out based on an upwind second-order scheme and k-w-SST(Shear Stress Transport) turbulence modeling has been used for estimation of eddy viscosity. Time-dependent flow structure results for the smooth casing reveal that there are two criteria for spike stall inception known as leading edge spillage and trailing edge backflow, which occur at specific mass flow rates in near-stall conditions. In this case, two dominant stall cells of different sizes could be observed. The larger one is caused by the spike stall covering roughly two blade passages in the circumferential direction and about 25% span in the radial direction. Spike stall disturbances are accompanied by lower frequencies and higher amplitudes of the pressure signals. Casing treatment causes flow blockages to reduce due to alleviation of backflow regions, which in turn reduces the total pressure loss and increases the axial velocity in the blade tip gap region, as well as tip leakage flow fluctuation at higher frequencies and lower amplitudes. Eventually, it can be concluded that the casing treatment of the stepped tip gap type could increase the stall margin of the compressor. This fact is basically due to retarding the movement of the interface region between incoming and tip leakage flows towards the rotor leading edge plane and suppressing the reversed flow around the blade trailing edge.     

Characteristics of pulsed plasma synthetic jet and its control effect on supersonic flow

The plasma synthetic jet is a novel flow control approach which is currently being studied. In this paper its characteristic and control effect on supersonic flow is investigated both experimentally and numerically. In the experiment, the formation of plasma synthetic jet and its propagation velocity in quiescent air are recorded and calculated with time resolved schlieren method. The jet velocity is up to 100 m/s and no remarkable difference has been found after changing discharge parameters. When applied in Mach 2 supersonic flow, an obvious shockwave can be observed. In the modeling of electrical heating, the arc domain is not defined as an initial condition with fixed temperature or pressure, but a source term with time-varying input power density, which is expected to better describe the influence of heating process. Velocity variation with different heating efficiencies is presented and discussed and a peak velocity of 850 m/s is achieved in still air with heating power density of 5.0 · 1012W/m3. For more details on the interaction between plasma synthetic jet and supersonic flow, the plasma synthetic jet induced shockwave and the disturbances in the boundary layer are numerically researched. All the results have demonstrated the control authority of plasma synthetic jet onto supersonic flow.     

New constitutive relation for eddy viscosity models

A new turbulent constitutive relation was directly derived from Boussinesq's hypothesis and mixing length theory, and then implemented in the standard k-ε model. The performance of this constitutive relation was validated in zero pressure gradient flat-plate boundary layer flow, fully-developed turbulent channel flow and separated flow in a plane asymmetric diffuser. The investigation demonstrated that, this new constitutive relation gave very accurate results in the former two basic cases and provided significant improvement in prediction of separated and reattachment points in the plane asymmetric diffuser. Separation and reattachment points at x/H=7.5 and 29 were calculated accurately in comparison to experimental results, and the static pressure coefficient of 0.82 was very close to large eddy simulation calculation. These results are very encouraging but further verification and extensive application of the new constitutive relation to other two-equation eddy viscosity model are needed.      

Influence of coupled characteristics on the internally staged combustor under low power operation

An ultra low emissions combustor, namely low emission stirred swirl (LESS) combustor was studied, based on a scheme of internally staged/lean premixed and prevaporized (LPP) combustion. The LESS combustor consists of central pilot stage and outer surrounded coaxially main stage, between which there exists a physical isolation, namely the step height. The existence of step height delayed the pilot and main jets mixing. Experimental and numerical studies were carried out to investigate the influence of the step height on the combustion performance. A single dome rectangular combustor was utilized to conduct the lean lightoff and blowout experiments, and pollutant emission experiments. The experimental results showed that with the increase of step height by 38%, the lean lightoff and blowout fuel air ratio decreased by 574% and 375%, the NOx emission increased by 35.1%, and the combustion efficiency increased by 1.78%; while the CO,unburned hydrocarbons (UHC) emissions decreased. Furthermore, the total pressure loss was kept nearly constant. Non reacting and reacting flow fields were numerically investigated to analyze the coupled characteristics of pilot and main jets with different step heights. A comparison of flow characteristics, spray structure, and combustion component as well as temperature field with different step heights was conducted. The numerical results showed that the increase of the step height shifted the peak velocity outwards. The enlargement of the primary recirculation zone (PRZ) resulted in the increase of the combustion efficiency and NOx emission, while the CO, UHC emissions decreeased.      

Simulation of underexpanded supersonic jet flows with chemical reactions

To achieve a detailed understanding of underexpanded supersonic jet structures influenced by afterburning and other flow conditions, the underexpanded turbulent supersonic jet with and without combustions are investigated by computational fluid dynamics(CFD) method.A program based on a total variation diminishing(TVD) methodology capable of predicting complex shocks is created to solve the axisymmetric expanded Navier–Stokes equations containing transport equations of species. The finite-rate ratio model is employed to handle species sources in chemical reactions. CFD solutions indicate that the structure of underexpanded jet is typically influenced by the pressure ratio and afterburning. The shock reflection distance and maximum value of Mach number in the first shock cell increase with pressure ratio. Chemical reactions for the rocket exhaust mostly exist in the mixing layer of supersonic jet flows. This tends to reduce the intensity of shocks existing in the jet, responding to the variation of thermal parameters.     

Radical recombination in a hydrocarbon-fueled scramjet nozzle

To reveal the radical recombination process in the scramjet nozzle flow and study the effects of various factors of the recombination, weighted essentially non-oscillatory(WENO)schemes are applied to solve the decoupled two-dimensional Euler equations with chemical reactions to simulate the hydrocarbon-fueled scramjet nozzle flow. The accuracy of the numerical method is verified with the measurements obtained by a shock tunnel experiment. The overall model length is nearly 0.5 m, with inlet static temperatures ranging from 2000 K to 3000 K, inlet static pressures ranging from 75 k Pa to 175 k Pa, and inlet Mach numbers of 2.0 ± 0.4 are involved.The fraction Damkohler number is defined as functions of static temperature and pressure to analyze the radical recombination progresses. Preliminary results indicate that the energy releasing process depends on different chemical reaction processes and species group contributions. In hydrocarbon-fueled scramjet nozzle flow, reactions with H have the greatest contribution during the chemical equilibrium shift. The contrast and analysis of the simulation results show that the radical recombination processes influenced by inflow conditions and nozzle scales are consistent with Damkohler numbers and potential dissociation energy release. The increase of inlet static temperature improves both of them, thus making the chemical non-equilibrium effects on the nozzle performance more significant. While the increase of inlet static pressure improves the former one and reduces the latter, it exerts little influence on the chemical non-equilibrium effects.     

Low-frequency unsteadiness of vortex wakes over slender bodies at high angle of attack

A type of flow unsteadiness with low frequencies and large amplitude was investigated experimentally for vortex wakes around an ogive-tangent cylinder. The experiments were carried out at angles of attack of 60–80 and subcritical Reynolds numbers of 0.6–1.8×105. The reduced frequencies of the unsteadiness are between 0.038 and 0.072, much less than the frequency of Karman vortex shedding. The unsteady flow induces large fluctuations of sectional side forces. The results of pressure measurements and particle image velocimetry indicate that the flow unsteadiness comes from periodic oscillation of the vortex wakes over the slender body. The time-averaged vortex patterns over the slender body are asymmetric, whose orientation is dependent on azimuthal locations of tip perturbations. Therefore, the vortex oscillation is a type of unsteady oscillation around a time-averaged asymmetric vortex structure.     

等离子体气动激励抑制压气机叶栅角区流动分离的仿真与实验

进行了等离子体气动激励抑制低速压气机叶栅角区流动分离的数值仿真研究,并进行了实验验证.小攻角情况下,叶片吸力面角区流动分离导致显著的尾迹总压损失.来流速度为50 m/s(雷诺数为223 000)时,等离子体气动激励可以有效的抑制角区流动分离,降低总压损失.激励电压、频率分别为10 kV和22 kHz时,50%叶高处的尾迹压力分布基本不变,60%和70%叶高处的最大总压损失分别减小了13.83%和10.74%.增加激励电极组数或激励电压,可以增强抑制效果.      

Experimental investigations on the power extraction of a turbine driven by a pulse detonation combustor

In order to grasp the interaction mechanism between the pulse detonation combustor（PDC）and the turbine,the experimental work in this paper investigates the key factors on the power extraction of a turbocharger turbine driven by a PDC.A PDC consisting of an unvalved tube is integrated with a turbocharger turbine which has a nominal mass flow rate of 0.6 kg/s and50000 r/min.The PDC-turbine hybrid engine is operated on gasoline-air mixtures and runs for6＋min to achieve a thermal steady state,and then the engine performance is evaluated under different operating conditions.Results show that the momentum difference per unit area between the turbine inlet and outlet plays an important role in the power extraction,while the pressure peak of the detonation has little effect.The equivalence ratio of fuel and air mixture and the transition structure between PDC and turbine are also important to the power extraction of the turbine.The present work is promising as it suggests that the performance beneft of a PDC-turbine hybrid engine can be realized by increasing the momentum difference per unit area through the optimal design of transition section between the PDC and turbine.     

扫频式射流对设计工况压气机叶栅流动分离影响的数值研究

高负荷压气机中的大尺度流动分离是导致其性能下降的主要原因,通过数值方法研究了扫频式射流控制角区分离、减小气动损失的效果,并以模型方程代替实际射流器,讨论了扫频式射流的基本控制参数对压气机叶栅流场控制效果及气动性能的影响。结果表明:扫频式射流使得流场呈现出稳定的周期性变化趋势,且存在一个扫频频率,使得超过该频率后的控制效果趋于稳定;合理选择扫频激励参数对实现流动分离的控制至关重要。在本文的方案中,采用较小的扫频射流角和射流流速能取得较好的控制效果,而更大的最大扫频摆角能强化这种控制效果,时均总压损失最大减小6.1%;扫频式射流能够在更大范围内提高吸力面边界层低能流体的动能,以更好地限制角区分离沿叶高方向发展,从而改善对角区分离的控制效果。     

基于PaSR模型的低旋流燃烧大涡模拟

为了深入理解低旋流流场特征和燃烧稳定性,基于OpenFOAM平台,采用动态k方程模型和有限速率PaSR模型对甲烷/空气预混低旋流燃烧进行了大涡模拟,研究了气流入口速度、当量比和压力等流场参数对流场结构和燃烧非稳态特性的影响,分析了流场大尺度结构与火焰相互作用。结果表明,流场结构和火焰抬升高度受入口速度影响较小,流场和火焰形态能够保持自相似性;随着当量比和压力提高,流场扩张性增强并在燃烧区下游产生回流区,火焰稳定不依赖回流区,根部火焰锋面形状由U形转变为W形,火焰抬升高度降低。火焰锋面稳定在剪切层,剪切层产生的周期性有序涡结构引起当地流场速度脉动和火焰表面褶皱,反映了流场非稳态特性;通过剪切层监测点瞬时轴向速度分析,涡结构特征频率随速度增大而提高,由250Hz提高至300Hz,随当量比和压力提高而降低,由250Hz降低至125Hz。     

引入压力梯度的针对角区分离流动的Spalart-Allmaras模型的改进

针对Spalart-Allmaras（S-A）模型在角区分离计算中的问题,将无量纲的压力梯度引入其涡黏性输运方程的生成项,得到了改进的S-A模型.通过对两套含角区分离的低速压气机叶栅进行验证计算发现:与实验结果相比,原始S-A模型所得的分离区偏大,分离区内壁面压力偏低;而改进模型得到了与实验一致的分离区尺寸以及吸力面、压力面压力系数分布等结果.针对S-A模型涡黏性生成项和耗散项的分析表明:引入的无量纲压力梯度有效的识别了角区分离,在分离区内改变了涡黏性的生成、耗散关系,增大了涡黏性,从而缩小了计算所得分离区,同时在主流区保留了原始S-A模型的计算结果,进而带来了良好的改进效果.      

用SPIV测量结果分析叶轮机内流动损失的方法初探

基于大涡模拟和PIV(particle image velocimetry)测量的相似性,利用SPIV(stereoscopic PIV)测量结果,可以使用亚格子应力模型求解流场中的湍流耗散率.对比不同亚格子应力模型求出的结果,混合模型求出的结果较为准确,但需要合理选择模型中的系数C_Sm和CAm.利用该方法分析低速大尺寸压气机试验台转子近叶尖区域的SPIV测量结果,发现在设计状态流场中的损失主要源于叶尖泄漏涡,而在近失速状态则主要源于叶尖泄漏涡和角区旋涡.