SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Overall stage and stator blade element performance with straight stator and tilted stator in transonic axial flow compressor stage
The paper presents a study of low aspect blading for the inlet stages of a high pressure ratio, high-speed core compressor. The basic overall design variables were stage pressure ratio and blade aspect ratio; these four stages represent two levels of total pressure ratio, two levels of rotor blade aspect ratio, and two levels of stator vane aspect ratios. Comparisons of the overall performance, radial distributions of performance parameters, diffusion factors at the near-stall conditions, blade element data, and the axial distribution of rotor tip static pressures yielded the following results: (1) higher peak pressure ratio, high stage and rotor efficiencies, and greater stall margin were obtained with the lower aspect ratio blading, (2) the lower aspect ratio blading showed improved performance over the entire blade span, and (3) the lower aspect ratio rotors operated at higher diffusion factors and higher incidence angles over the entire blade span.
Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.
An investigation has been performed to assess methods used for rotor inflow modeling. A key element of this assessment has been the recent acquisition of high quality experimental measurements of inflow velocities taken in the proximity of a lifting rotor in forward flight. Widely used rotor performance predictive methods are based on blade element strip theory coupled with an inflow model. The inflow prediction models assessed in this paper include the uniform inflow based on momentum, a skewed disk model, and two methods based on a vortex wake structure.
It is well known that airfoils under unsteady flow conditions with a periodically varying angle of attack exhibit aerodynamic characteristics different from those under steady flow conditions, a phenomenon commonly known as dynamic stall. It is also well known that the steady aerodynamic characteristics of airfoils in the inboard region of a rotating blade differ from those under steady two-dimensional (2D) flow conditions, a phenomenon commonly known as rotational augmentation. This paper presents an investigation of these two phenomena together in the inboard parts of wind turbine blades. This analysis is carried out using data from three sources: (1) the National Renewable Energy Laboratory’s Unsteady Aerodynamics Experiment Phase VI experimental data, including constant as well as continuously pitching blade conditions during axial operation, (2) data from unsteady Delayed Detached Eddy Simulations (DDES) carried out using the Technical University of Denmark’s in-house flow solver Ellipsys3D, and (3) data from a simplified model based on the blade element momentum method with a dynamic stall subroutine that uses rotationally augmented steady-state polars obtained from steady Phase VI experimental sequences, instead of the traditional 2D nonrotating data. The aim of this work is twofold. First, the blade loads estimated by the DDES simulations are compared to three select cases of the N sequence experimental data, which serves as a validation of the DDES method. Results show reasonable agreement between the two data in two out of three cases studied. Second, the dynamic time series of the lift and the moment polars obtained from the experiments are compared to those from the dynamic stall subroutine that uses the rotationally augmented steady polars. This allowed the differences between the stall phenomenon on the inboard parts of harmonically pitching blades on a rotating wind turbine and the classic dynamic stall representation in 2D flow to be investigated. Results from the dynamic stall subroutine indicated a good qualitative agreement between the model and the experimental data in many cases, which suggests that the current 2D dynamic stall model as used in BEM-based aeroelastic codes may provide a reasonably accurate representation of three-dimensional rotor aerodynamics when used in combination with a robust rotational augmentation model.
The topic of laminar flow effects on hover performance is introduced with respect to some historical efforts where laminar flow was either measured or attempted. An analysis method is outlined using combined blade element, momentum method coupled to an airfoil analysis method, which includes the full e(sup N) transition model. The analysis results compared well with the measured hover performance including the measured location of transition on both the upper and lower blade surfaces. The analysis method is then used to understand the upper limits of hover efficiency as a function of disk loading. The impact of laminar flow is higher at low disk loading, but significant improvement in terms of power loading appears possible even up to high disk loading approaching 20 ps f. A optimum planform design equation is derived for cases of zero profile drag and finite drag levels. These results are intended to be a guide for design studies and as a benchmark to compare higher fidelity analysis results. The details of the analysis method are given to enable other researchers to use the same approach for comparison to other approaches.
The aerodynamic design and the overall blade element performance of a 51 centimeter diameter fan stage is presented. The stage was designed to minimize the noise generated by rotor stator interactions. The design pressure ratio was 1.20 at a flow of 30.6 kilograms per second and a rotor blade tip speed of 228.6 meters per second. At design speed the rotor peak efficiency was 0.935. The peak efficiency of the stage, however, was 0.824. The radial distribution of rotor performance parameters at peak efficiency and design speed indicated excellent agreement with design values.
Interstage data recorded on a J85-13 engine were used to analyze the internal flow of the compressor. Measured pressures and temperatures were used as input to a streamline analysis program to calculate the velocity diagrams at the inlet and outlet of each blade row. From the velocity diagrams and blade geometry, selected blade-element performance parameters were calculated. From the detailed analysis it is concluded that the compressor is probably hub critical (stall initiates at the hub) in the latter stages for the design speed conditions. As a result, the casing treatment over the blade tips has little or no effect on stall margin at design speed. Radial inlet distortion did not appear to change the flow in the stages that control stall because of the rapid attenuation of the distortion within the compressor.
The coupled flap-lag-torsion equations of motion of an isolated horizontal axis wind turbine blade are formulated. Quasi-steady blade-element strip theory was applied to derive the aerodynamic operator which includes boundary layer type gradient winds. The final equations which have periodic coefficients were solved in order to obtain the aeroelastic response and stability of large horizontal axis wind turbine blade. A new method of generating an appropriate time-dependent equilibrium position (required for the stability analysis) has been implemented. Representative steady-state responses and stability boundaries, applicable mainly to an existing blade design (NASA/-ERDA MOD-0), are presented. The results indicate that the MOD-0 configuration is a basically stable design and that blade stability is not sensitive to offsets between blade elastic axis and aerodynamic center. Blade stability appears to be sensitive to precone. The tower shadow (or wake) has a considerable effect on the flap response but leaves blade stability unchanged. Finally, it was found that non linear terms in the equations of motion can significantly affect the linearized stability boundaries, however, these terms have a negligible effect on blade response at operating conditions.
The Thermal/Structural Tailoring of Engine Blades (T/STAEBL) system is a family of computer programs executed by a control program. The T/STAEBL system performs design optimizations of cooled, hollow turbine blades and vanes. This manual contains an overview of the system, fundamentals of the data block structure, and detailed descriptions of the inputs required by the optimizer. Additionally, the thermal analysis input requirements are described as well as the inputs required to perform a finite element blade vibrations analysis.
A new time-domain method for simulating cyclostationary turbulence as seen by a translating and rotating blade element has recently been developed for the case of one-dimensional spectral distribution. This paper extends the simulation method to the cases of two- and three-dimensional spectral distributions and presents validation results for the two-dimensional case. The statistics of an isolated rigid blade flapping response to turbulence are computed using a two-dimensional spectral representation of the von Karman turbulence model, and the results are compared with those obtained using the conventional space-fixed turbulence analysis. The new turbulence simulation method is used for predicting the Black Hawk helicopter response to atmospheric turbulence.
In the present research, tilt rotor aeroacoustics have been studied experimentally and computationally. Experimental measurements were made on a 1/12.5 scale model. A dimensional analysis showed that the model was a good aeroacoustic approximation to the full-scale aircraft, and scale factors were derived to extrapolate the model measurements to the full-scale XV-15. The experimental measurements included helium bubble flow visualization, silk tuft flow visualization, 2-component hot wire anemometry, 7-hole pressure probe measurements, vorticity measurements, and outdoor far field acoustic measurements. The hot wire measurements were used to estimate the turbulence statistics of the flow field into the rotors, such as length scales, velocity scales, dissipation, and turbulence intermittency. To date, these flow measurements are the only ones in existence for a hovering tilt rotor. Several different configurations of the model were tested: (1) standard configurations (single isolated rotor, two rotors without the aircraft, standard tilt rotor configuration); (2) flow control devices (the 'plate', the 'diagonal fences'); (3) basic configuration changes (increasing the rotor/rotor spacing, reducing the rotor plane/wing clearance, operating the rotors out of phase). Also, an approximation to Sikorsky's Variable Diameter Tilt Rotor (VDTR) configuration was tested, and some flow measurements were made on a semi-span configuration of the model. Acoustic predictions were made using LOWSON.M, a Mathematica code. This hover prediction code, from HOVER.FOR, used blade element theory for the aerodynamics, and Prandtl's Vortex theory to model the wake, along with empirical formulas for the effects of Reynolds number, Mach number, and stall. Aerodynamic models were developed from 7-hole pressure probe measurements of the mean velocity into the model rotors. LOWSON.M modeled a rotor blade as a single force and source/sink combination separated in the chordwise direction, at an effective blade radius. Spanwise, Mach-weighted integrals were used to find the equivalent forces and equivalent source strengths.
Accurate, reliable and robust numerical predictions of wind turbine rotor power remain a challenge to the wind energy industry. The literature reports various methods that compare predictions to experiments. The methods vary from Blade Element Momentum Theory (BEM), Vortex Lattice (VL), to variants of Reynolds-averaged Navier-Stokes (RaNS). The BEM and VL methods consistently show discrepancies in predicting rotor power at higher wind speeds mainly due to inadequacies with inboard stall and stall delay models. The RaNS methodologies show promise in predicting blade stall. However, inaccurate rotor vortex wake convection, boundary layer turbulence modeling and grid resolution has limited their accuracy. In addition, the inherently unsteady stalled flow conditions become computationally expensive for even the best endowed research labs. Although numerical power predictions have been compared to experiment. The availability of good wind turbine data sufficient for code validation experimental data that has been extracted from the IEA Annex XIV download site for the NREL Combined Experiment phase II and phase IV rotor. In addition, the comparisons will show data that has been further reduced into steady wind and zero yaw conditions suitable for comparisons to "steady wind" rotor power predictions. In summary, the paper will present and discuss the capabilities and limitations of the three numerical methods and make available a database of experimental data suitable to help other numerical methods practitioners validate their own work.
Radial-dimension computations in the RSRA's blade-element model are modified for both the acquisition of extensive baseline data and for real-time simulation use. The baseline data, which are for the evaluation of model changes, use very small increments and are of high quality. The modifications to the real-time simulation model are for accuracy improvement, especially when a minimal number of blade segments is required for real-time synchronization. An accurate technique for handling tip loss in discrete blade models is developed. The mathematical consistency and convergence properties of summation algorithms for blade forces and moments are examined and generalized integration coefficients are applied to equal-annuli midpoint spacing. Rotor conditions identified as 'constrained' and 'balanced' are used and the propagation of error is analyzed.
A simple analytical model to account for fuselage-induced velocities at rotor blade elements and at rotor wake nodes is described. The method is applied to three different fuselage configurations. Results obtained with a comprehensive rotor code show the fuselage effect on rotor trim controls, comparing the isolated rotor with inclusion of the fuselage for the same trim. This is compared to a simple analytical estimate of the fuselage effect using blade element/momentum theory. It is found that in forward flight the lateral control is mainly affected by fuselage effects. Rotor thrust can be varied by the presence of the fuselage, depending on its angle of attack, and the fuselage influence generally increases with flight speed.
A turbulence model has been developed for realtime aircraft simulation that accommodates stochastic turbulence and distributed discrete gusts as a function of the terrain. This model is applicable to conventional aircraft, V/STOL aircraft, and disc rotor model helicopter simulations. Vehicle angular activity in response to turbulence is computed from geometrical and temporal relationships rather than by using the conventional continuum approximations that assume uniform gust immersion and low frequency responses. By using techniques similar to those recently developed for blade-element rotor models, the angular-rate filters of conventional turbulence models are not required. The model produces rotational rates as well as air mass translational velocities in response to both stochastic and deterministic disturbances, where the discrete gusts and turbulence magnitudes may be correlated with significant terrain features or ship models. Assuming isotropy, a two-dimensional vertical turbulence field is created. A novel Gaussian interpolation technique is used to distribute vertical turbulence on the wing span or lateral rotor disc, and this distribution is used to compute roll responses. Air mass velocities are applied at significant centers of pressure in the computation of the aircraft's pitch and roll responses.
The aerodynamic test of a small, single stage, highly loaded, axial flow transonic compressor is covered. The stage was modified by fabricating a 24 blade rotor with mis-set blades in a repeating pattern - two degrees closed from nominal, two degrees open from nominal and nominal. The unit was instrumented to determine overall performance and average blade element data. High-response, dynamic pressure probes were installed to record pressure patterns at selected points in the flowpath. Testing was conducted at speeds from 70 to 94% of design equivalent speed with a conventional casing and also with circumferential grooves over the rotor tip. Testing indicated severe performance penalties were incurred as a result of the mis-set blading. Lower flow, pressure ratio, and efficiency were observed for the stage with or without casing treatment. Periodic pressure variations were detected at every location where high response pressure sensors were located and were directly related to blading geometry.
Current fast aeroelastic wind turbine codes suitable for certification lack an induction model for standstill conditions. A trailed vorticity model previously used as addition to a blade element momentum theory based aerodynamic model in normal operation has been extended to allow computing the induced velocities in standstill. The model is validated against analytical results for an elliptical wing in constant inflow and against stand still measurements from the NREL/NASA Phase VI unsteady experiment. The extended model obtains good results in case of the elliptical wing, but underpredicts the steady loading for the Phase VI blade in attached flow. The prediction of the dynamic force coefficient loops from the Phase VI experiment is improved by the trailed vorticity modeling in both attached flow and stall in most cases. The exception is the tangential force coefficient in stall, where the codes and measurements deviate and no clear improvement is visible.