Shock–boundary layer interaction and energetics in transonic flutter

2017 ◽  
Vol 832 ◽  
pp. 212-240 ◽  
Author(s):  
Pradeepa T. Karnick ◽  
Kartik Venkatraman

We study the influence of shock and boundary layer interactions in transonic flutter of an aeroelastic system using a Reynolds-averaged Navier–Stokes (RANS) solver together with the Spalart–Allmaras turbulence model. We show that the transonic flutter boundary computed using a viscous flow solver can be divided into three distinct regimes: a low transonic Mach number range wherein viscosity mimics increasing airfoil thickness thereby mildly influencing the flutter boundary; an intermediate region of drastic change in the flutter boundary due to shock-induced separation; and a high transonic Mach number zone of no viscous effects when the shock moves close to the trailing edge. Inviscid transonic flutter simulations are a very good approximation of the aeroelastic system in predicting flutter in the first and third regions: that is when the shock is not strong enough to cause separation, and in regions where the shock-induced separated region is confined to a small region near the trailing edge of the airfoil. However, in the second interval of intermediate transonic Mach numbers, the power distribution on the airfoil surface is significantly influenced by shock-induced flow separation on the upper and lower surfaces leading to oscillations about a new equilibrium position. Though power contribution by viscous forces are three orders of magnitude less than the power due to pressure forces, these viscous effects manipulate the flow by influencing the strength and location of the shock such that the power contribution by pressure forces change significantly. Multiple flutter points that were part of the inviscid solution in this regime are now eliminated by viscous effects. Shock motion on the airfoil, shock reversal due to separation, and separation and reattachment of flow on the airfoil upper surface, also lead to multiple aerodynamic forcing frequencies. These flow features make the flutter boundary quantitatively sensitive to the turbulence model and numerical method adopted, but qualitatively they capture the essence of the physical phenomena.

Author(s):  
P. J. Bryanston-Cross ◽  
J. J. Camus

A simple technique has been developed which samples the dynamic image plane information of a schlieren system using a digital correlator. Measurements have been made in the passages and in the wakes of transonic turbine blades in a linear cascade. The wind tunnel runs continuously and has independently variable Reynolds and Mach number. As expected, strongly correlated vortices were found in the wake and trailing edge region at 50 KHz. Although these are strongly coherent we show that there is only limited cross-correlation from wake to wake over a Mach no. range M = 0.5 to 1.25 and variation of Reynolds number from 3 × 105 to 106. The trailing edge fluctuation cross correlations were extended both upstream and downstream and preliminary measurements indicate that this technique can be used to obtain information on wake velocity. The vortex frequency has also been measured over the same Mach number range for two different cascades. The results have been compared with high speed schlieren photographs.


Author(s):  
Marion Mack ◽  
Roland Brachmanski ◽  
Reinhard Niehuis

The performance of the low pressure turbine (LPT) can vary appreciably, because this component operates under a wide range of Reynolds numbers. At higher Reynolds numbers, mid and aft loaded profiles have the advantage that transition of suction side boundary layer happens further downstream than at front loaded profiles, resulting in lower profile loss. At lower Reynolds numbers, aft loading of the blade can mean that if a suction side separation exists, it may remain open up to the trailing edge. This is especially the case when blade lift is increased via increased pitch to chord ratio. There is a trend in research towards exploring the effect of coupling boundary layer control with highly loaded turbine blades, in order to maximize performance over the full relevant Reynolds number range. In an earlier work, pulsed blowing with fluidic oscillators was shown to be effective in reducing the extent of the separated flow region and to significantly decrease the profile losses caused by separation over a wide range of Reynolds numbers. These experiments were carried out in the High-Speed Cascade Wind Tunnel of the German Federal Armed Forces University Munich, Germany, which allows to capture the effects of pulsed blowing at engine relevant conditions. The assumed control mechanism was the triggering of boundary layer transition by excitation of the Tollmien-Schlichting waves. The current work aims to gain further insight into the effects of pulsed blowing. It investigates the effect of a highly efficient configuration of pulsed blowing at a frequency of 9.5 kHz on the boundary layer at a Reynolds number of 70000 and exit Mach number of 0.6. The boundary layer profiles were measured at five positions between peak Mach number and the trailing edge with hot wire anemometry and pneumatic probes. Experiments were conducted with and without actuation under steady as well as periodically unsteady inflow conditions. The results show the development of the boundary layer and its interaction with incoming wakes. It is shown that pulsed blowing accelerates transition over the separation bubble and drastically reduces the boundary layer thickness.


2012 ◽  
Vol 134 (5) ◽  
Author(s):  
Hans-Ju¨rgen Rehder

As part of a European research project, the aerodynamic and thermodynamic performance of a high pressure turbine cascade with different trailing edge cooling configurations was investigated in the wind tunnel for linear cascades at DLR in Go¨ttingen. A transonic rotor profile with a relative thick trailing edge was chosen for the experiments. Three trailing edge cooling configurations were applied, first central trailing edge ejection, second a trailing edge shape with a pressure side cut-back and slot equipped with a diffuser rib array, and third pressure side film cooling through a row of cylindrical holes. For comparison, aerodynamic investigations on a reference cascade with solid blades (no cooling holes or slots) were performed. The experiments covered the subsonic, transonic and supersonic exit Mach number range of the cascade while varying cooling mass flow ratios up to 2 %. This paper analyzes the effect of coolant ejection on the airfoil losses. Emphasis was given on separating the different loss contributions due to shocks, pressure, and suction side boundary layer, trailing edge, and mixing of the coolant flow. Employed measurement techniques are schlieren visualization, blade surface pressure measurements, and traverses by pneumatic probes in the cascade exit flow field and around the trailing edge. The results show that central trailing edge ejection significantly reduces the mixing losses and therefore decreases the overall loss. Higher loss levels are obtained when applying the configurations with pressure side blowing. In particular, the cut-back geometry reveals strong mixing losses due to the low momentum coolant fluid, which is decelerated by the diffuser rib array inside the slot. The influence of coolant flow rate on the trailing edge loss is tremendous, too. Shock and boundary layer losses are major contributions to the overall loss but are less affected by the coolant. Finally a parameter variation changing the temperature ratio of coolant to main flow was performed, resulting in increasing losses with decreasing coolant temperature.


2000 ◽  
Vol 123 (2) ◽  
pp. 409-417 ◽  
Author(s):  
Bjo¨rn Gru¨ber ◽  
Volker Carstens

A parametric study which investigates the influence of viscous effects on the damping behavior of vibrating compressor cascades is presented here. To demonstrate the dependence of unsteady aerodynamic forces on the flow viscosity, a computational study was performed for a transonic compressor cascade of which the blades underwent tuned pitching oscillations while the flow conditions extended from fully subsonic to highly transonic flow. Additionally, the reduced frequency and Reynolds number were varied. In order to check the linear behavior of the aerodynamic forces, all calculations were carried out for three different oscillation amplitudes. Comparisons with inviscid Euler results helped identify the influence of viscous effects. The computations were performed with a Navier-Stokes code, the basic features of which are the use of an AUSM upwind scheme, an implicit time integration, and the implementation of the Baldwin-Lomax turbulence model. In order to demonstrate the possibility of this code to correctly predict the unsteady behavior of strong shock-boundary layer interactions, the experiment of Yamamoto and Tanida on a self-induced shock oscillation due to shock-boundary layer interaction was calculated. A significant improvement in the prediction of the shock amplitude was achieved by a slight modification of the Baldwin Lomax turbulence model. An important result of the presented compressor cascade investigations is that viscous effects may cause a significant change in the aerodynamic damping. This behavior is demonstrated by two cases in which an Euler calculation predicts a damped oscillation whereas a Navier-Stokes computation leads to an excited vibration. It was found that the reason for these contrary results are shock-boundary-layer interactions which dramatically change the aerodynamic damping.


1996 ◽  
Vol 118 (2) ◽  
pp. 240-248 ◽  
Author(s):  
Y. T. Lee ◽  
C. Hah ◽  
J. Loellbach

Steady-state analyses of the incompressible flow past a single-stage stator/rotor propulsion pump are presented and compared to experimental data. The purpose of the current study is to validate a numerical method for the design application of a typical propulsion pump and for the acoustic analysis based on predicted flowfields. A steady multiple-blade-row approach is used to calculate the flowfields of the stator and the rotor. The numerical method is based on a fully conservative control-volume technique. The Reynolds-averaged Navier–Stokes equations are solved along with the standard two-equation k–ε turbulence model. Numerical results for both mean flow and acoustic properties compare well with measurements in the wake of each blade row. The rotor blade has a thick boundary layer in the last quarter of the chord and the flow separates near the trailing edge. These features invalidate many Euler prediction results. Due to the dramatic reduction of the turbulent eddy viscosity in the thick boundary layer, the standard k–ε model cannot predict the correct local flow characteristics near the rotor trailing edge and in its near wake. Thus, a modification of the turbulence length scale in the turbulence model is applied in the thick boundary layer in response to the reduction of the turbulent eddy viscosity.


1973 ◽  
Vol 60 (2) ◽  
pp. 257-271 ◽  
Author(s):  
G. T. Coleman ◽  
C. Osborne ◽  
J. L. Stollery

A hypersonic gun tunnel has been used to measure the heat transfer to a sharpedged flat plate inclined at various incidences to generate local Mach numbers from 3 to 9. The measurements have been compared with a number of theoretical estimates by plotting the Stanton number against the energy-thickness Reynolds number. The prediction giving the most reasonable agreement throughout the above Mach number range is that due to Fernholz (1971).The values of the skin-friction coefficient derived from velocity profiles and Preston tube data are also given.


2018 ◽  
Vol 8 (11) ◽  
pp. 2320 ◽  
Author(s):  
Hidemi Takahashi ◽  
Mitsuru Kurita ◽  
Hidetoshi Iijima ◽  
Monami Sasamori

Turbulent boundary layer profiles on the aircraft surface were characterized by pitot-rake measurements conducted in flight experiments at high subsonic Mach number ranges. Due to slight variations in atmospheric air conditions or aircraft attitudes, such as angles of attack and absolute flight speeds at different flights even under the same premised flight conditions, the boundary layer profiles measured at different flights can exhibit different shape and velocity values. This concern leads to difficulty in evaluating the efficiency of using some kind of drag-controlling device such as riblets in the flight test, since the evaluation would be conducted by comparing the profiles measured with and without using riblets at different flights. An approach was implemented to interpolate the boundary layer profile for a flight condition of interest based on the response surface method, in order to eliminate the influence of the flight conditional difference. Results showed that the interpolation with the 3rd-degree response surface model with a combination of two independent variables of flight Mach number and total pressure successfully eliminated the influence of the flight conditional difference, and interpolated the boundary layer profiles measured at different flights within an inaccuracy of 4.1% for the flight Mach number range of 0.5 to 0.78.


1995 ◽  
Vol 282 ◽  
pp. 339-371 ◽  
Author(s):  
S. J. Leib ◽  
Sang Soo Lee

We study the nonlinear evolution of a pair of oblique instability waves in a supersonic boundary layer over a flat plate in the nonlinear non-equilibrium viscous critical layer regime. The instability wave amplitude is governed by the same integro-differential equation as that derived by Goldstein & Choi (1989) in the inviscid limit and by Wu, Lee & Cowley (1993) with viscous effects included, but the coefficient appearing in this equation depends on the mean flow and linear neutral stability solution of the supersonic boundary layer. This coefficient is evaluated numerically for the Mach number range over which the (inviscid) first mode is the dominant instability. Numerical solutions to the amplitude equation using these values of the coefficient are obtained. It is found that, for insulated and cooled wall conditions and angles corresponding to the most rapidly growing waves, the amplitude ends in a singularity at a finite downstream position over the entire Mach number range regardless of the size of the viscous parameter. The explosive growth of the instability waves provides a mechanism by which the boundary layer can break down. A new feature of the compressible problem is the nonlinear generation of a spanwise-dependent mean distortion of the temperature along with that of the velocity found in the incompressible case.


1966 ◽  
Vol 24 (1) ◽  
pp. 1-31 ◽  
Author(s):  
H. T. Nagamatsu ◽  
B. C. Graber ◽  
R. E. Sheer

An investigation was conducted in a hypersonic shock tunnel to study the laminar boundary-layer transition on a highly cooled 10° cone of 4 ft. length over the Mach-number range of 8·5 to 10·5 with a stagnation temperature of 1400 °K. The effects on transition of tip surface roughness, tip bluntness, and ± 2° angle of attack were investigated. With fast-response, thin film surface heat-transfer gauges, it was possible to detect the passage of turbulent bursts which appeared at the beginning of transition. Pitot-tube surveys and schlieren photographs of the boundary layer were obtained to verify the interpretation of the heat-transfer data. It was found that the surface roughness greatly promoted transition in the proper Reynolds-number range. The Reynolds numbers for the beginning and end of transition at the 8·5 Mach-number location were 3·8 × 106−9·6 × 106and 2·2 × 106−4·2 × 106for the smooth sharp tip and rough sharp tip respectively. The local skin-friction data, determined from the Pitot-tube survey, agreed with the heat-transfer data obtained through the modified Reynolds analogy. The tip-bluntness data showed a strong delay in the beginning of transition for a cone base-to-tip diameter ratio of 20, approximately a 35% increase in Reynolds number over that of the smooth sharp-tip case. The angle-of-attack data indicated the cross flow to have a strong influence on transition by promoting it on the sheltered side of the cone and delaying it on the windward side.


Author(s):  
Shun He ◽  
Shijun Guo ◽  
Wenhao Li

An investigation into transonic flutter characteristic of an airfoil conceived with the morphing leading and trailing edges has been carried out. Computational fluid dynamics (CFD) is used to calculate the unsteady aerodynamic force in transonic flow. An aerodynamic reduced order model (ROM) based on autoregressive model with exogenous input (ARX) is used in the numerical simulation. The flutter solution is determined by eigenvalue analysis at specific Mach number. The approach is validated by comparing the transonic flutter characteristics of the Isogai wing with relevant literatures before applied to a morphing airfoil. The study reveals that by employing the morphing trailing edge, the shock wave forms and shifts to the trailing edge at a lower Mach number, and aerodynamic force stabilization happens earlier. Meanwhile, the minimum flutter speed increases and transonic dip occurs at a lower Mach number. It is also noted that leading edge morphing has negligible effect on the appearance of the shock wave and transonic flutter. The mechanism of improving the transonic flutter characteristics by morphing technology is discussed by correlating shock wave location on airfoil surface, unsteady aerodynamics with flutter solution.


Sign in / Sign up

Export Citation Format

Share Document