[Low-Frequency Flow Oscillation

The results of the research conducted under this grant are presented in detail in three Master theses, by Heinrich, Balow, and Broeren. Additional analysis of the experimental data can be found in two AIAA Journal articles and two conference papers. Citations for all of the studies' publications can be found in the bibliography which is attached. The objective of Heinrich's study was to document the low-frequency flow oscillation on the LRN-1007 airfoil, which had been previously observed at low Reynolds number, to determine its origin, and explore the phenomenon at higher Reynolds number. Heinrich performed detailed flow visualization on the airfoil using surface fluorescent oil and laser-sheet off-body visualization. A large leading-edge separation bubble and trailing-edge separation was identified on the airfoil just prior to the onset of the unsteady stall flow oscillation. From the laser-sheet data, the unsteady flow appeared as a massive boundary-layer separation followed by flow reattachment. Hot-wire data were taken in the wake to identify the presence of the flow oscillation and the dominant frequency. The oscillation was found in the flow from a Reynolds number of 0.3 to 1.3 x 10 exp 6. The Strouhal number based on airfoil projected height was nominally 0.02 and increased slightly with increasing Reynolds number and significantly with increasing airfoil angle of attack. Balow focused his research on the leading-edge separation bubble which was hypothesized to be the origin of the low-frequency oscillation. Initially, experimental measurements in the bubble at the onset of the low-frequency oscillation were attempted to study the characteristics of the bubble and explain possible relationships to the shear-layer-flapping phenomena. Unfortunately, the bubble proved to be extremely sensitive to the probe interference and it drastically reduced the size of the bubble. These detailed measurements were then abandoned by Balow. However, this led to a series of tests where the leading-edge bubble and trailing-edge separation were altered and the affect on the flow-oscillation studied. Balow found that by tripping the airfoil boundary-layer with "zigzag" tape ahead of bubble separation, the bubble was effectively eliminated mid the oscillation suppressed. Wake survey drag measurements showed a drastic reduction in airfoil drag when the bubble and oscillation were eliminated. Using the "zigzag" tape, the trailing-edge separation was moved downstream approximately 5 percent chord. This was found to reduce the amplitude of the oscillation, particularly in the onset stage at low angle of attack (around 14 degrees). Through detailed analysis of the wake behind the airfoil during the unsteady flow oscillation, Balow provided a better understanding of the wake flowfield. Broeren studied the oscillating flowfield in detail at Reynolds number equal 3 x 10 exp 5 and an angle of attack of 15 degrees using laser Doppler velocimetry (LDV). Two-dimensional LDV data were acquired at 687 grid points above the model upper surface while hot-wire data were taken simultaneously in the wake. Using the hot-wire signal, the LDV data were phase averaged into 24 bins to represent a single ensemble average of one oscillation cycle. The velocity data showed a flowfield oscillation that could be divided into three flow regimes. In the first regime, the flow over the airfoil was completely separated initially, the flowfield reattached from the leading edge and the reattachment point moved downstream with increasing time or phase. Broeren referred to this as the reattachment regime. The bubble development regime followed, where a leading-edge separation bubble formed at the leading edge and grew with increasing time. During the initial part of this regime the trailing-edge separation continued to move downstream. However, during the last 30 degrees of phase the trailing-edge separation moved rapidly forward and appeared to merge with the leading-edge bubble. During the third regime, the separation regime, the flow was segmented from the airfoil leading edge and did not reattach to the airfoil surface. The reverse flow was seen to grow in vertical extent up from the model surface as the phase increased. Next reattachment began again at the leading edge signaling the start of the reattachment regime, and so the cycle continued. From Broeren's work, the details of the unsteady flowfield over the airfoil were seen for the first time. From this research a great deal has been learned about the low-frequency flow oscillation which naturally occurs on the LRN-1007 airfoil near stall. The oscillation was seen to persist at higher Reynolds number, the dependence of the Strouhal number on angle of attack and Reynolds number were discovered, the critical role played by the laminar bubble was shown and the entire upper surface flowfield during a flow oscillation cycle was measured and analyzed. What still eludes understanding is the scaling of the flow oscillation and why certain airfoils, such as the LRN, have a very strong low-frequency mode and other airfoils exhibit no organized low-frequency oscillation at all

Effects of Ice Accretion on Aircraft Aerodynamics

The primary objective of this research was to support the development of a new ice accretion model by improving our physical understanding of the ice accretion process through experimental measurements. The focus was on the effect of the initial ice roughness (smooth/rough boundary) on the accretion process. This includes understanding the boundary-layer development over the roughness and especially its effect on the heat transfer which is fundamental to the ice accretion process. The research focused on acquiring the experimental data needed to formulate a new ice accretion physical model. Research was conducted to analyze boundary-layer data taken on a NACA 0012 airfoil with roughness to simulate the smooth/rough boundary. The effect of isolated roughness on boundary-layer transition was studied experimentally to determine if the classical critical roughness Reynolds number criteria could be applied to transition in the airfoil leading-edge area. The effect of simulated smooth/rough boundary roughness on convective heat transfer was studied to complete the study. During the course of this research the effect of free-stream wind tunnel turbulence on the boundary layer was measured. Since this quantity was not well known, research to accurately measure the wind tunnel turbulence in an icing cloud was undertaken. Preliminary results were attained and the final data were acquired, reduced and presented under a subsequent grant

An Experimental Study of the Aerodynamics of a Swept and Unswept Semispan Wing with a Simulated Glaze Ice Accretion

Two semispan wings, one with a rectangular planform and one with 30 degrees of leading edge sweep were tested. Both had a NACA 0012 airfoil section, and both were tested clean and with simulated glaze ice shapes on their leading edges. Several surface roughness were tested. Each model geometry is documented and each surface roughness is explained. Aerodynamic performance of the wing in the form of sectional lift and integrated three-dimensional lift is documented through pressure measurements obtained from rows of surface pressure taps placed at five span locations on the wing. For the rectangular wing, sectional drag near the midspan is obtained from wake total pressure profiles. The data is presented in tabular and graphical form and is also available on computer disk

The 3-D LDV Measurements on a 30-Degree Swept Wing with a Simulated Ice Accretion

Three dimensional flowfield measurements have been obtained for a semispan 30-degree swept wing with a simulated glaze ice accretion. The model tested has a NACA 0012 section perpendicular to the leading edge. Measurements were made using a two-component laser Doppler velocimeter (LDV) system. Mean velocity measurements were obtained for all three velocity components. Streamwise turbulence intensities were also obtained. All measurements were taken in the University of Illinois 3 by 4 foot subsonic wind tunnel at a Reynolds number of 1 million and 8 degrees angle of attack. The data is presented in tabular form

The Effects of the Critical Ice Accretion on Airfoil and Wing Performance

In support of the NASA Lewis Modern Airfoils Ice Accretion Test Program, the University of Illinois at Urbana-Champaign provided expertise in airfoil design and aerodynamic analysis to determine the aerodynamic effect of ice accretion on modern airfoil sections. The effort has concentrated on establishing a design/testing methodology for "hybrid airfoils" or "sub-scale airfoils," that is, airfoils having a full-scale leading edge together with a specially designed and foreshortened aft section. The basic approach of using a full-scale leading edge with a foreshortened aft section was considered to a limited extent over 40 years ago. However, it was believed that the range of application of the method had not been fully exploited. Thus a systematic study was being undertaken to investigate and explore the range of application of the method so as to determine its overall potential

Airfoil Boundary-Layer Development and Transition with Large Leading-Edge Roughness

An experimental study of the effects of large distributed roughness located near the leading edge of an airfoil has been performed to determine the effect on boundary-layer development and transition. Boundary-layer measurements were carried out on a two-dimensional NACA 0012 airfoil with a 53.34-cm chord through the use of hot-wire anemometry at Reynolds numbers of 0.75 £ 10 6 , 1.25 £ 10 6 , and 2.25 £ 10 6 . These measurements included mean and¯uctuating velocity, turbulence intensity,¯ow® eld intermittency, and associated integral parameters. The roughness used was of the type and density observed to occur during the initial glaze ice accretion process. Results have shown that the transitional boundary layer induced by large distributed roughness is markedly different from the smooth model Tollmein± Schlicting induced transition process. No fully developed turbulent boundary layers were observed to occur near the roughness location. Instead, the large distributed roughness was observed to trigger a transitional boundary layer at or very near the roughness location. This transitional boundary layer required a substantial chordwise extent to obtain a fully developed turbulent state. Streamwise turbulence intensity levels in the roughness induced transitional region were observed to be relatively low as compared with the smooth model transitional region. Nomenclatur

Aerodynamic Classification of Swept-Wing Ice Accretion

The continued design, certification and safe operation of swept-wing airplanes in icing conditions rely on the advancement of computational and experimental simulation methods for higher fidelity results over an increasing range of aircraft configurations and performance, and icing conditions. The current state-of-the-art in icing aerodynamics is mainly built upon a comprehensive understanding of two-dimensional geometries that does not currently exist for fundamentally three-dimensional geometries such as swept wings. The purpose of this report is to describe what is known of iced-swept-wing aerodynamics and to identify the type of research that is required to improve the current understanding. Following the method used in a previous review of iced-airfoil aerodynamics, this report proposes a classification of swept-wing ice accretion into four groups based upon unique flowfield attributes. These four groups are: ice roughness, horn ice, streamwise ice and spanwise-ridge ice. In the case of horn ice it is shown that a further subclassification of nominally 3D or highly 3D horn ice may be necessary. For all of the proposed ice-shape classifications, relatively little is known about the three-dimensional flowfield and even less about the effect of Reynolds number and Mach number on these flowfields. The classifications and supporting data presented in this report can serve as a starting point as new research explores swept-wing aerodynamics with ice shapes. As further results are available, it is expected that these classifications will need to be updated and revised

Aerodynamic Simulation of Runback Ice Accretion

This report presents the results of recent investigations into the aerodynamics of simulated runback ice accretion on airfoils. Aerodynamic tests were performed on a full-scale model using a high-fidelity, ice-casting simulation at near-flight Reynolds (Re) number. The ice-casting simulation was attached to the leading edge of a 72-in. (1828.8-mm ) chord NACA 23012 airfoil model. Aerodynamic performance tests were conducted at the ONERA F1 pressurized wind tunnel over a Reynolds number range of 4.7?10(exp 6) to 16.0?10(exp 6) and a Mach (M) number ran ge of 0.10 to 0.28. For Re = 16.0?10(exp 6) and M = 0.20, the simulated runback ice accretion on the airfoil decreased the maximum lift coe fficient from 1.82 to 1.51 and decreased the stalling angle of attack from 18.1deg to 15.0deg. The pitching-moment slope was also increased and the drag coefficient was increased by more than a factor of two. In general, the performance effects were insensitive to Reynolds numb er and Mach number changes over the range tested. Follow-on, subscale aerodynamic tests were conducted on a quarter-scale NACA 23012 model (18-in. (457.2-mm) chord) at Re = 1.8?10(exp 6) and M = 0.18, using low-fidelity, geometrically scaled simulations of the full-scale castin g. It was found that simple, two-dimensional simulations of the upper- and lower-surface runback ridges provided the best representation of the full-scale, high Reynolds number iced-airfoil aerodynamics, whereas higher-fidelity simulations resulted in larger performance degrada tions. The experimental results were used to define a new subclassification of spanwise ridge ice that distinguishes between short and tall ridges. This subclassification is based upon the flow field and resulting aerodynamic characteristics, regardless of the physical size of the ridge and the ice-accretion mechanism

Generation of Fullspan Leading-Edge 3D Ice Shapes for Swept-Wing Aerodynamic Testing

The deleterious effect of ice accretion on aircraft is often assessed through dry-air flight and wind tunnel testing with artificial ice shapes. This paper describes a method to create fullspan swept-wing artificial ice shapes from partial span ice segments acquired in the NASA Glenn Icing Reserch Tunnel for aerodynamic wind-tunnel testing. Full-scale ice accretion segments were laser scanned from the Inboard, Midspan, and Outboard wing station models of the 65% scale Common Research Model (CRM65) aircraft configuration. These were interpolated and extrapolated using a weighted averaging method to generate fullspan ice shapes from the root to the tip of the CRM65 wing. The results showed that this interpolation method was able to preserve many of the highly three dimensional features typically found on swept-wing ice accretions. The interpolated fullspan ice shapes were then scaled to fit the leading edge of a 8.9% scale version of the CRM65 wing for aerodynamic wind-tunnel testing. Reduced fidelity versions of the fullspan ice shapes were also created where most of the local three-dimensional features were removed. The fullspan artificial ice shapes and the reduced fidelity versions were manufactured using stereolithography

Effect of Airfoil Geometry on Performance with Simulated Intercycle Ice Accretions

This paper presents the results of an experimental study designed to evaluate the performance effects of intercycle ice accretions on airfoils with different geometries. The intercycle ice accretions were simulated using combinations of various size grit roughness. These simulations were tested on three airfoils: NACA 23012, NACA 3415, and NLF 0414 at a Reynolds number of 1.8 × × 10 6 and a Mach number of 0.18. Results from the NACA 23012 airfoil tests closely matched those from a previous study, validating the ice-shape simulation method. This also showed that a simple geometric (chord-based) scaling of the ice was appropriate. The simulated ice effect, in terms of maximum lift performance, was most severe for the NACA 23012 airfoil. The maximum lift coefficients were in the range of 0.65 to 0.80 for the iced configuration compared to a clean value of 1.47 for the NACA 23012 airfoil at this Reynolds number. In contrast, the maximum lift coefficients for the NLF 0414 airfoil with the same ice simulations were in the range of 0.90 to 1.05, compared to a clean value of 1.34. The results for the NACA 3415 with the simulated intercycle ice shapes were between the other two airfoils