钝后缘风力机翼型的环量控制研究

钝后缘风力机翼型具有结构强度高、对表面污染不敏感等优点,但其较大的阻力系数使得翼型的整体气动特性不够理想. 利用环量控制方法对钝后缘风力机翼型进行了流动控制,以改善钝后缘风力机翼型的气动特性,减弱尾迹区脱体涡强度. 通过对钝后缘风力机翼型环量控制方法进行相关的数值模拟,对比研究了环量控制方法的增升减阻效果, 研究了环量控制下翼型升阻力特性随射流动量系数的变化规律,并对不同射流动量系数下环量控制方法的气动品质因子和控制效率进行了分析. 研究结果表明:环量控制方法能够大幅提升钝后缘风力机翼型的升力系数,同时有效地降低翼型的阻力系数; 翼型的升力系数随射流动量系数的增大而增大,表现出很明显的分离控制阶段和超环量控制阶段的变化规律; 射流能耗的功率系数随射流动量系数的增大而增大,且增长速率逐渐增大;实施环量控制方法后叶片的输出功率同样随射流动量系数增大而增大,但增长速率逐渐降低. 总体来说,环量控制方法可以有效地改善钝后缘风力机翼型的气动特性以及功率输出特性,在大型风力机流动控制中具有很好的应用前景.      

Numerical simulation of an axisymmetric separated and reattached flow over a longitudinal blunt circular cylinder

This article presents a numerical investigation of turbulent flow in an axisymmetric separated and reattached flow over a longitudinal blunt circular cylinder. The governing equations were discretized by the finite-volume method and SIMPLER method was applied to solve the equations on a staggered grid. The turbulent flow was numerically simulated using the standard k–ε, Abe–Kondoh–Nagano (AKN) and Shear Stress Transport (SST) turbulence models. The comparisons made between numerical results and experimental measurements showed that the SST model is superior to other models in the present calculation.Computations were performed for three different Reynolds numbers of 6000, 10 000 and 20 000 based on the cylinder diameter. To our knowledge, this study represents the first numerical investigation of the present flow configuration. The computational results were validated with the available experimental data of reattachment length, mean velocity distribution and wall static pressure coefficient in the turbulent blunt circular cylinder flows. Further, other characteristics of the flow, such as turbulent kinetic energy, pressure, streamlines, and the velocity vectors are discussed.The results show that the main characteristics of the turbulence flow in the separation region, such as reattachment length or velocity profiles, are nearly independent of the Reynolds number. The obtained results showed that a secondary separation bubble may appear in the main separation bubble near the leading edge. Furthermore, it was found that the turbulent kinetic energy has a large effect on the formation of the secondary bubble.     

Direct Numerical Simulation of Turbulent Flow over a Rectangular Trailing Edge

This paper describes a direct numerical simulation (DNS) study of turbulent flow over a rectangular trailing edge at a Reynolds number of 1000, based on the freestream quantities and the trailing edge thickness h; the incoming boundary layer displacement thickness δ^* is approximately equal to h. The time-dependent inflow boundary condition is provided by a separate turbulent boundary layer simulation which is in good agreement with existing computational and experimental data. The turbulent trailing edge flow simulation is carried out using a parallel multi-block code based on finite difference methods and using a multi-grid Poisson solver. The turbulent flow in the near-wake region of the trailing edge has been studied first for the effects of domain size and grid resolution. Then two simulations with a total of 256 × 512 × 64 (∼ 8.4×10⁶) and 512 × 1024 × 128 (∼ 6.7×10⁷) grid points in the computational domain are carried out to investigate the key flow features. Visualization of the instantaneous flow field is used to investigate the complex fluid dynamics taking place in the near-wake region; of particular importance is the interaction between the large-scale spanwise, or Kármán, vortices and the small-scale quasi-streamwise vortices contained within the inflow boundary layer. Comparisons of turbulence statistics including the mean flow quantities are presented, as well as the pressure distributions over the trailing edge. A spectral analysis applied to the force coefficient in the wall normal direction shows that the main shedding frequency is characterized by a Strouhal number based on h of approximately 0.118. Finally, the turbulence kinetic energy budget is analysed. Received 4 March 1999 and accepted 27 October 2000     

Numerical simulation of an impinging jet on a flat plate

Two-dimensional normal impinging jet flowfields, with or without an upper plate, were analysed by employing an implicit bidiagonal numerical method developed by Lavante and Thompkins Jr. The Jones–Launder K–? two-equation turbulent model was employed to study the turbulent effects of the impinging jet flowfield. The upper plate surface pressure, the ground plane pressure and other physical parameters of the momentum flowfield were calculated at various jet exit height and jet inlet Reynolds numbers. These results were compared with those of Beam and Warming's numerical method, Hsiao and Chuang, and others, along with experimental data. The potential core length of the impinging jet without an upper plate is longer than that of the free jet because of the effects of the ground plane, while the potential core length of the impinging jet with an upper plate is shorter than that of the free jet because of the effects of the upper plate. This phenomenon in the present analysis provides a fundamental numerical study of an impinging jet and a basis for further analysis of impinging jet flowfields on a variable angle plate.     

Development of a Two-velocities Hybrid RANS-LES Model and its Application to a Trailing Edge Flow

The flow around a trailing edge is computed with a new hybrid method designed to more clearly separate the effects of total and sub-grid turbulent stress-modelling on the time-averaged and instantaneous velocity fields, and in turn, mean momentum and kinetic energy balances. These two velocity fields independently define Reynolds averaged and sub-grid-scale viscosities, and distinct stresses, at the same location. In particular, resolved eddies can emerge, or sweep in and out of the Reynolds averaged near wall layer, without being dampened by higher levels of the viscosity in this RANS dominated layer. The two-field hybrid model, first tested on channel flows, gives accurate predictions of mean velocities and stresses for different Reynolds numbers and coarse meshes. For the trailing edge flow the results of the hybrid model are close to the reference fine LES for mean velocity and turbulent content, whereas the DES-SST on the same coarse mesh gives too early separation.     

网栅后湍流流动沿流程衰减特性的实验研究

应用ＬＤＶ测试技术对垂直方管内三种不同栅距的网所形成的湍流流动做了详细测量,网栅栅距Ｍ分别为１４ｍｍ,１０ｍｍ,６ｍｍ实验给出了管内不同位置湍流流动的各种参数沿流动方向的衰减规律,并就三种网栅进行比较分析,以探讨这种网栅流动所具有共同规律和在不同栅距对湍流结构及脉动大小的不同特点,从而对这种网栅所形成的湍流流结构较为清晰的了解。     

绕振荡水翼流动及其转捩特性的数值计算研究简

通过对比标准k-ω SST 湍流模型和基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型对绕振荡NACA66 水翼流动的数值计算结果与实验结果,对水翼振荡过程的水动力特性和流场结构变化进行了分析研究. 结果表明：与标准k-ω SST 湍流模型的数值计算结果相比,基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型能有效预测绕振荡翼型流场结构和水动力特性,捕捉流场边界层发生的流动分离和转捩现象;绕振荡水翼的流动过程可分为5 个特征阶段,当来流攻角较小时,在水翼前缘发生层流向湍流的转捩现象,水翼动力特征曲线出现变化拐点;随着来流攻角的增大,顺时针尾缘涡逐渐形成并向水翼前缘发展;当攻角较大时,前缘涡分离导致动力失速,水翼的动力特征曲线出现大幅波动;水翼处于顺时针向下旋转阶段,绕水翼的流动状态逐渐由湍流过渡为层流.     

Effects of a vectored trailing edge jet on delta wing vortex breakdown

Flow visualization was used to study the effects of a vectored trailing edge jet on the leading edge vortex breakdown of a 65° delta wing. The experimental results indicated that there is little effect of the jet on the leading edge vortex breakdown when the angle of the vectored jet is less than 10°. With the increase of the vectored angle ß, the effect of the jet on the flow becomes stronger, i.e., the jet delays the leading edge vortex breakdown in the direction of the vectored jet, and accelerates breakdown of the other leading edge vortex. Moreover, the effect of the jet control tends to be weaker with the angle of attack.     

The effect of a single vortex generator jet on the characteristics of a turbulent boundary layer

A flat plate experiment was performed in a water tunnel to determine the effects of a vortex generator jet on the characteristics of a turbulent boundary layer at various wall normal locations. The results show that the characteristic distributions of the turbulent fluctuation quantities are nearly unaffected by the induced vortex structures neither in the steady nor in the dynamic blowing case. The shear layer interaction between the turbulent main flow and the jet flow produces less turbulent fluctuations than it is expected from a turbulent free jet flow. Thus, the mixing process of this flow control strategy is based only on a large-scale momentum transport superimposed by the turbulent fluctuation quantities. This allows a separation of scales for physical interpretation and numerical simulations.     

A discrete-vortex model for the arbitrary motion of a thin airfoil with fluidic control

A low-order model for the arbitrary motion of a thin airfoil with trailing edge fluidic control is derived from basic fluid mechanics principles. The model consist of solving a single ordinary differential equation with a special treatment of a vortex shedding criteria. The model is compared with experimental and high-order numerical simulations and the results give a reasonable means of predicting the lift and moment on a thin airfoil. Furthermore, the model is extended to account for the actuation and control due to the synthetic jet actuation near the trailing edge. The model response is compared to experimental results.     

Asymmetric vortex shedding flow past an inclined flat plate at high incidence

This paper reports an experimental investigation of the vortex shedding wake behind a long flat plate inclined at a small angle of attack to a main flow stream. Detailed velocity fields are obtained with particle-image velocimetry (PIV) at successive phases in a vortex shedding cycle at three angles of attack, α=20°, 25° and 30°, at a Reynolds number Re≈5,300. Coherent patterns and dynamics of the vortices in the wake are revealed by the phase-averaged PIV vectors and derived turbulent properties. A vortex street pattern comprising a train of leading edge vortices alternating with a train of trailing edge vortices is found in the wake. The trailing edge vortex is shed directly from the sharp trailing edge while there are evidences that the formation and shedding of the leading edge vortex involve a more complicated mechanism. The leading edge vortex seems to be shed into the wake from an axial location near the trailing edge. After shedding, the vortices are convected downstream in the wake with a convection speed roughly equal to 0.8 the free-stream velocity. On reaching the same axial location, the trailing edge vortex, as compared to the leading edge vortex, is found to possess a higher peak vorticity level at its centre and induce more intense fluid circulation and Reynolds stresses production around it. It is found that the results at the three angles of attack can be collapsed into similar trends by using the projected plate width as the characteristic length of the flow.     

Experimental study of transient forced turbulent separation and reattachment on a bevelled trailing edge

The purpose of this experimental study is to analyze the transient processes of the separation and reattachment of a turbulent boundary layer at the trailing edge of a splitter plate. The separation is driven by a steady pneumatic injection, considered here as an actuator controlling the mixing layer. Particle image velocimetry and hot-film measurements are performed at various stages of the processes. The results highlight the behavior of each transient. Analysis of the time scale of the processes is realized by means of bi-orthogonal decomposition. This gives essential information (e.g., time scale, qualitative features) for future applications using duty cycle mode to excite instabilities of the mixing layer.     

Flow and mixing characteristics of swirling double-concentric jets subject to acoustic excitation

Characteristic flow modes, flow evolution processes, jet spread width, turbulence properties, and dispersion characteristics of swirling double-concentric jets were studied experimentally. Jet pulsations were induced by means of acoustic excitation. Streak pictures of smoke flow patterns, illuminated by a laser-light sheet, were recorded by a high-speed digital camera. A hot-wire anemometer was used to digitize instantaneous velocity instabilities in the flow. Jet spread width was obtained through a binary edge identification technique. Tracer-gas concentrations were measured for information on jet dispersions. Two characteristic flow patterns were observed: (1) synchronized vortex rings appeared in the low excitation intensity regime (the excitation intensity less than one) and (2) synchronized puffing turbulent jets appeared in the high excitation intensity regime (the excitation intensity greater than one). In the high excitation intensity regime, the “suction back” phenomenon occurred and therefore induced in-tube mixing. The jet spread width and turbulent fluctuation intensity exhibited particularly large values in the high excitation intensity regime at the excitation Strouhal numbers smaller than 0.85. At the excitation Strouhal numbers >0.85, the high-frequency effect caused significant decay of jet breakup and dispersion—the jet spread width and fluctuation intensity decreased sharply and may, at very high Strouhal numbers, asymptotically approach values almost the same as the values associated with unexcited jets. Exciting the jets at the high excitation intensity regime, the effects of puffing motion and in-tube mixing caused breakup of the jet in the near field and therefore resulted in a small Lagrangian integral time and small length scales of fluctuating eddies. This effect, in turn, caused drastic dispersion of the central jet fluids. It is possible that the excited jets can attain 90 % more improvements than the unexcited jets. We provide a domain regarding excitation intensity and Strouhal number to facilitate identification of characteristic flow modes.     

Vortical structures behind a transverse jet in a supersonic flow at high jet to crossflow pressure ratios

A three-dimensional supersonic turbulent flow with symmetric normal injection of circular jets from the channel walls is numerically simulated. The initial Favre-averaged Navier–Stokes equations closed by the k–ω turbulence model are solved by an algorithm based on an ENO scheme. The mechanism of the formation of vortical structures due to the interaction of the jet with the free stream is studied for jet to crossflow total pressure ratios ranging from 3 to 50. It is known from experiments reported in the literature that, for n ? 10, mixing of the jet with the high-velocity flow leads to the formation of a pair of vortices and of an additional separation zone near the wall behind the jet. It is demonstrated that the present numerical results are consistent with such findings and that the pressure distribution on the wall ahead of the jet in the plane of symmetry is also in reasonable agreement with available experimental data.     

Influence of the Trailing Edge of a Cavity on the Flow Fluctuation Intensity

Unsteady supersonic turbulent gas flow over a plane cavity is simulated numerically. The solution is determined by means of the joint integration of the Reynolds equations and the two-parameter turbulence model using the law of the wall. The calculation results obtained for various cavity trailing edge contours are considered. The strong effect of the geometry of the trailing edge on the flow parameter fluctuation intensity is demonstrated. A surface shape which makes it possible considerably to reduce these oscillations is determined.     

The use of turbulent inflow conditions for the modelling of a high aspect ratio jet

The aim of this work is to investigate the capabilities of a turbulent inflow method. The application selected for this study is the high aspect ratio jet. The complexities associated to the numerical modelling of a high aspect ratio jet are embedded in its physical complexity. Consequently, the numerical modelling does not only require a high mesh resolution, but furthermore it requires a careful mesh construction, inflow conditions and subgrid-scale modelling to make an accurate computation of the unsteady flow phenomena. The results indicate that increased grid resolution and enhanced turbulence modelling reduce the effect of the imposed flow fluctuations. It is concluded that for a high aspect ratio free-jet turbulent inflow conditions are effective if the mesh resolution is insufficient to trigger shear-layer instabilities. Applied with sufficient mesh resolution the onset of vortex motions will occur in the shear layer, hence there is limited inflow sensitivity.     

Spatiotemporal representation of the dynamic modes in turbulent cavity flows

Proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) were used to extract the coherent structures in turbulent cavity flows. The spatiotemporal representation of the modes was achieved by performing the circular convolution of a change of basis on the data sequence, wherein the transformation function was extracted from the POD or DMD. The spatiotemporal representation of the modes provided significant insight into the evolutionary behavior of the structures. Self-sustained oscillations arise in turbulent cavity flows due to unsteady separation at the leading edge. The turbulent cavity flow at Re_D = 12,000 and a length to depth ratio L/D = 2 was analyzed. The dynamic modes extracted from the data clarified the presence of self-sustained oscillations. The spatiotemporal representation of the POD and DMD modes that caused self-sustained oscillations revealed the prevalent dynamics and evolutionary behavior of the coherent structures from their formation at the leading edge to their impingement at the trailing edge. A local minimum in the mode amplitude representing the energy contributions to the flow was observed upon the impingement of coherent structure at the trailing edge. The modal energy associated with the periodic formation of organized coherent structures followed by their dissipation upon impingement revealed the oscillatory behavior over time.     

Numerical simulations of non‐equilibrium turbulent boundary layer flowing over a bump

Large‐eddy simulation (LES) and Reynolds‐averaged Navier–Stokes simulation (RANS) with different turbulence models (including the standard k?ε, the standard k?ω, the shear stress transport k?ω (SST k?ω), and Spalart–Allmaras (S–A) turbulence models) have been employed to compute the turbulent flow of a two‐dimensional turbulent boundary layer over an unswept bump. The predictions of the simulations were compared with available experimental measurements in the literature. The comparisons of the LES and the SST k?ω model including the mean flow and turbulence stresses are in satisfied agreements with the available measurements. Although the flow experiences a strong adverse pressure gradient along the rear surface, the boundary layer is unique in that intermittent detachment occurring near the wall. The numerical results indicate that the boundary layer is not followed by mean‐flow separation or incipient separation as shown from the numerical results. The resolved turbulent shear stress is in a reasonable agreement with the experimental data, though the computational result of LES shows that its peak is overpredicted near the trailing edge of the bump, while the other used turbulence models, except the standard k?ε, underpredicts it. Analysis of the numerical results from LES confirms the experimental data, in which the existence of internal layers over the bump surface upstream of the summit and along the downstream flat plate. It also demonstrates that the quasi‐step increase in skin friction is due to perturbations in pressure gradient. The surface curvature enhances the near‐wall shear production of turbulent stresses, and is responsible for the formation of the internal layers. The aim of the present work is to examine the response and prediction capability of LES with the dynamic eddy viscosity model as a sub‐grid scale to the complex turbulence structure with the presence of streamline curvature generated by a bumpy surface. Aiming to reduce the computational costs with focus on the mean behavior of the non‐equilibrium turbulent boundary layer of flow over the bump surface, the present investigation also explains the best capability of one of the used RANS turbulence models to capture the driving mechanism for the surprisingly rapid return to equilibrium over the trailing flat plate found in the measurements. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulation numérique directe de l'influence de la forme aval d'une plaque séparatrice sur une couche de mélangeDirect-numerical simulation of the splitting-plate downstream-shape influence upon a mixing layer

In this Note, we present direct numerical simulation results of a spatial mixing layer generated behind an upstream plate separating two boundary layers. The effect of the shape of the trailing edge of the plate is considered through comparisons between flows obtained from a bevelled or a blunt plate. In the former case, a spatial mixing layer consistent with previous experimental and numerical observations is obtained. In the latter case, the self-excited state that establishes in the near wake region dominates primary and secondary instability mechanisms while understating the importance of inflow perturbations. This behaviour is interpreted in terms of convective or absolute instability. The effects on turbulent statistics are also discussed. To cite this article: S. Laizet, E. Lamballais, C. R. Mecanique 334 (2006).