Colin S. Bidwell

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) multi-element probe are presented. Simulation of the flow field was produced using the Glenn-HT Navier-Stokes solver. Three-dimensional unsteady results were produced and then time averaged for the heat transfer and collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Simulations were completed for free stream velocities ranging from 85-135 meters per second, and free stream total pressure of 44.8 and 93.1 kilopascals (6.5 and 13.5 pounds per square inch absolute). In addition, the effect of angle of attack and yaw were investigated by including 5 degree deviations from straight for one of the flow conditions. All but one of the cases simulated a probe in isolation (i.e. in a very large domain without any support strut). One case is included which represents a probe mounted on a support strut within a finite sized wind tunnel. Collection efficiencies were generated, using the LEWICE3D code, for four spherical particle sizes, 100, 50, 20, and 5 micron in diameter. It was observed that a reduction in velocity of about 20% occurred, for all cases, as the flow entered the shroud of the probe. The reduction in velocity within the shroud is not indicative of any error in the probe measurement accuracy. Heat transfer results are presented which agree quite well with a correlation for the circular cross section heated elements. Collection efficiency results indicate a reduction in collection efficiency as particle size is reduced. The reduction with particle size is expected, however, the results tended to be lower than the previous results generated for isolated two-dimensional elements. The deviation from the two-dimensional results is more pronounced for the smaller particles and is likely due to the reduced flow within the protective shroud. As particle size increases differences between the two-dimensional and three dimensional results become negligible. Taken as a group, the total collection efficiency of the elements including the effects of the shroud has been shown to be in the range of 0.93 to 0.99 for particles above 20 microns. The 3D model has improved the estimated collection efficiency for smaller particles where errors in previous estimates were more significant.

Collection Efficiency and Ice Accretion Calculations for a Boeing 737-300 Inlet

Collection efficiency and ice accretion calculations have been made for a Boeing 737-300 inlet using a three-dimensional panel code, an adaptive grid code, the NASA Lewis LEWICE3D grid based ice accretion code. Flow solutions for the inlet were generated using the VSAERO panel code. Grids used in the ice accretion calculations were generated using the newly developed adaptive grid code ICEGRID3D. The LEWICE3D grid based ice accretion program was used to calculate impingement efficiency and ice shapes. Ice shapes typifying rime and mixed icing conditions were generated for a 30 minute hold condition. All calculations were performed on an SGI Power Challenge computer. The results have been compared to experimental flow and impingement data. In general, the calculated flow and collection efficiencies compared well with experiment, and the ice shapes looked reasonable and appeared representative of the rime and mixed icing conditions for which they were calculated.

Three Dimensional Simulation of Flow in an Axial Low Pressure Compressor at Engine Icing Operating Points

Viscous three-dimensional simulations of the Honeywell ALF502R-5 low pressure compressor (sometimes called a booster) using the NASA Glenn code GlennHT have been carried out. A total of ten simulations were produced. Five operating points are investigated, with each point run with two different wall thermal conditions. These operating points are at, or near, points where engine icing has been determined to be likely. In the future, the results of this study will be used for further analysis such as predicting collection efficiency of ice particles and ice growth rates at various locations in the compressor. A mixing plane boundary condition is used between each blade row, resulting in convergence to steady state within each blade row. The k-omega turbulence model of Wilcox, combined with viscous grid spacing near the wall on the order of one, is used to resolve the turbulent boundary layers. For each of the operating points, heat transfer coefficients are generated on the blades and walls. The heat transfer coefficients are produced by running the operating point with two different wall thermal conditions and then solving simultaneously for the heat transfer coefficient and adiabatic wall temperature at each point. Average Nusselt numbers are calculated for the most relevant surfaces. The values are seen to scale with Reynolds number to approximately a power of 0.7. Additionally, images of surface distribution of Nusselt number are presented. Qualitative comparison between the five operating points show that there is relatively little change in the character of the distribution. The dominant observed effect is that of an overall scaling, which is expected due to Reynolds number differences. One interesting aspect about the Nusselt number distribution is observed on the casing (outer diameter) downstream of the exit guide vanes (EGVs). The Nusselt number is relatively high between the pairs of EGVs, with two lower troughs downstream of each EGV trailing edge. This is of particular interest since rather complex ice shapes have been observed in that region. 1 Senior Research Engineer, Vantage Partners, LLC. 2 Research Scientist, Department of Mechanical and Aerospace Engineering, The Ohio State University. 3 Aerospace Engineer, Turbomachinery and Heat Transfer Branch, NASA Glenn Research Center. 4 Aerospace Engineer, Turbomachinery and Heat Transfer Branch, NASA Glenn Research Center. https://ntrs.nasa.gov/search.jsp?R=20170009137 2019-11-11T12:03:45+00:00Z

Ice Particle Analysis of the Honeywell AL502 Engine Booster

A flow and ice particle trajectory analysis was performed for the booster of the Honeywell ALF502 engine. The analysis focused on two closely related conditions one of which produced an icing event and another which did not during testing of the ALF502 engine in the Propulsion Systems Lab (PSL) at NASA Glenn Research Center. The flow analysis was generated using the NASA Glenn GlennHT flow solver and the particle analysis was generated using the NASA Glenn LEWICE3D v3.63 ice accretion software. The inflow conditions for the two conditions were similar with the main differences being that the condition that produced the icing event was 6.8 K colder than the non-icing event case and the inflow ice water content (IWC) for the non-icing event case was 50% less than for the icing event case. The particle analysis, which considered sublimation, evaporation and phase change, was generated for a 5 micron ice particle with a sticky impact model and for a 24 micron median volume diameter (MVD), 7 bin ice particle distribution with a supercooled large droplet (SLD) splash model used to simulate ice particle breakup. The particle analysis did not consider the effect of the runback and re-impingement of water resulting from the heated spinner and anti-icing system. The results from the analysis showed that the amount of impingement for the components were similar for the same particle size and impact model for the icing and non-icing event conditions. This was attributed to the similar aerodynamic conditions in the booster for the two cases. The particle temperature and melt fraction were higher at the same location and particle size for the non-icing event than for the icing event case due to the higher incoming inflow temperature for the non-event case. The 5 micron ice particle case produced higher impact temperatures and higher melt fractions on the components downstream of the fan than the 24 micron MVD case because the average particle size generated by the particle breakup was larger than 5 microns which yielded less warming and melting. The analysis also showed that the melt fraction and wet bulb temperature icing criterion developed during tests in the Research Altitude Test Facility (RATFac) at the National Research Council (NRC) of Canada were useful in predicting icing events in the ALF502 engine. The development of an ice particle impact model which includes the effects of particle breakup, phase change, and surface state is necessary to further improve the prediction of ice particle transport with phase change through turbomachinery.

Collection Efficiency and Ice Accretion Characteristics of Two Full Scale And One 1/4 Scale Business Jet Horizontal Tails

Collection efficiency and ice accretion calculations have been made for a series of business jet horizontal tail configurations using a three-dimensional panel code, an adaptive grid code, and the NASA Glenn LEWICE3D grid based ice accretion code. The horizontal tail models included two full scale wing tips and a 25 percent scale model. Flow solutions for the horizontal tails were generated using the PMARC panel code. Grids used in the ice accretion calculations were generated using the adaptive grid code ICEGRID. The LEWICE3D grid based ice accretion program was used to calculate impingement efficiency and ice shapes. Ice shapes typifying rime and mixed icing conditions were generated for a 30 minute hold condition. All calculations were performed on an SGI Octane computer. The results have been compared to experimental flow and impingement data. In general, the calculated flow and collection efficiencies compared well with experiment, and the ice shapes appeared representative of the rime and mixed icing conditions for which they were calculated.

Experimental water droplet impingement data on modern aircraft surfaces

An experimental method has been developed to determine the water droplet impingement characteristics on two- and three-dimensional aircraft surfaces. The experimental water droplet impingement data are used to validate particle trajectory analysis codes that are used in aircraft icing analyses and engine inlet particle separator analyses. The aircraft surface is covered with thin strips of blotter paper in areas of interest. The surface is then exposed to an airstream that contains a dyed-water spray cloud. The water droplet impingement data are extracted from the dyed blotter paper strips by measuring the optical reflectance of each strip with an automated reflectometer. Preliminary experimental and analytical impingement efficiency data are presented for a NLF(1)-0414F airfoil, s swept MS(1)-0317 airfoil, a swept NACA 0012 wingtip and for a Boeing 737-300 engine inlet model.

Progress in Measuring Water Impingement Characteristics on Aircraft Surfaces

Validation of trajectory computer codes, for icing analysis, requires experimental water droplet impingement data for a wide range of aircraft geometries as well as flow and icing conditions. This paper provides a summary of experimental water droplet impingement research starting with the early efforts of NACA in the 1950s. Industry requirements for additional impingement data are discussed based on results from a recent industry survey. A new research program aimed at expanding and modernizing the existing water droplet impingement data is outlined. Improved experimental and data reduction methods for obtaining water impingement data are presented. A discussion on sources of error in the experimental and data reduction methods is provided. Relative humidity is shown to have a significant effect on the repeatability of the experimental data. Recent wind tunnel test results for an MS(1)-0317 airfoil and a three element McDonnell Douglas high lift system are compared with computational results obtained with the LEWICE code. In most cases, the correlation between the analytical and experimental impingement efficiency distributions is good. The experimental and analysis results include impingement characteristics for median volumetric diameters of 11.5, 21 and 92 microns. Test repeatability for the experimental results presented ranged from 4% to 10% with respect to the average data. This represents a significant improvement compared to previous experimental results.