Search NASASearch

Engineering topics

Tsao, Jen-Ching

Publications and source records attributed to Tsao, Jen-Ching.

At least 19 records

Analysis of Experimental Ice Accretion Data and Assessment of a Thermodynamic Model During Ice Crystal Icing

This paper evaluates a thermodynamic ice crystal icing model that has been previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. The model functions between two distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. Freeze-dominated icing occurs when liquid water (from melted ice crystals) freezes and accretes on a surface along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized as having strong adhesion to the surface. The amount of ice accretion is partially dictated by a freeze fraction, which is the fraction of impinging liquid water that freezes. Melt-dominated icing occurs as unmelted ice on a surface accumulates. This melt-dominated icing is characterized by weakly bonded surface adhesion. The amount of ice accumulation is partially dictated by a melt fraction, which is the fraction of impinging ice crystals that melts. Experimentally observed ice growth rates suggest that only a small fraction of the impinging ice remains on the surface, implying a mass loss mechanism such as splash, runback, bounce, or erosion. The fraction of mass loss must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface for a given ice growth rate. This mass loss parameter, however, along with the freeze fraction and melt fraction, are the only experimental parameters that are currently not measured directly. Using icing growth rates from ice crystal icing experiments, a methodology that has been previously proposed is used to determine these unknown parameters. This work takes ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) at the Glenn Research Center in 2018 that examined the fundamental physics of ice crystal icing. This paper continues evaluation of the thermodynamic model from a previous effort, with additions to the model that account for sub-freezing temperatures that have been observed at the leading edge of the airfoil during icing. The predicted temperatures were generally in good agreement with measured temperatures. Other key findings include the total wet-bulb temperature being a good first order indicator of whether icing is freeze-dominated (sub-freezing values) or melt-dominated (above freezing). Maximum sticking efficiency values, the fraction of impinging mass that adheres to a surface, was calculated to be about 0.2, and retained this maximum value for a range of melt ratios (0.3 to 0.65 and possibly higher), which is defined as the ratio of liquid water content to total water content. Higher air velocities reduced the maximum sticking efficiency and shifted the icing regime to higher melt ratio values. Finally, the leading edge ice accretion angle was found to be related to ice growth (lower growth rates for smaller angles) and melt ratio (smaller melt ratios resulted in smaller angles, likely due to erosion effects).

Ice-Crystal Icing

Scaling Evaluation of Ice-Crystal Icing on a Modern Turbofan Engine in PSL Using the COMDES-MELT Code

This paper presents preliminary ice-crystal icing (ICI) altitude scaling evaluation results of a Honeywell Uncertified Research Engine (HURE) that was tested in the NASA Glenn Research Center Propulsion Systems Laboratory (PSL) during January of 2018. This engine geometry features a hidden core design to keep the core less exposed. The engine was fitted with internal video cameras to observe various ice buildup processes at multiple selected locations within the engine core flow path covering the fan stator, the splitter-lip/shroud/strut, and the high pressure compressor (HPC) variable inlet guide vane (IGV) regions. The potential ice accretion risk was pre-determined to occur by using NASA’s in-house 1D Engine Icing Risk assessment code, COMDES-MELT. The code was successful in predicting the risk of ice accretion in adiabatic regions like the fan-stator of the HURE at specific engine operating points. However at several operating points during the test, liquid water was observed running along the shroud toward the variable IGV of the HPC regions with an air temperature well below freezing, thus no particle melting could have occurred due to heating from the air alone. It was reasoned that other sources of heat were present in that region. To account for these heat sources the inlet total temperature was adjusted to give a wet bulb temperature of 24 F below the standard minimum wet bulb temperature of 492 R to allow ice to accrete in the splitter-lip/shroud/strut region, which was determined from a reference case where hard ice was observed in that region. With that adjustment the COMDES-MELT code was successful in providing operating points where there was a risk of ice accretion during the test campaign. In addition to calculating possible conditions at different selected lower altitudes, simulations were run to determine potential inlet conditions that could lead to ice-crystal accretion along the prescribed stations where the cameras were available. From there, scaled test conditions were determined by best matching the following three icing related parameters of the reference condition: (1) the local air total wet bulb temperature, (2) the local ice crystal cloud melt ratio and (3) the engine fan face ice/water to air mass flux ratio of the ice crystal cloud. Instantaneous images taken from the time-lapsed movies of ice buildup were used along with the relevant thermodynamic data of air, water vapor and local icing condition to help evaluate how closely the proposed altitude scaling method could be used in ground based test facility to duplicate selected reference ICI features observed at specific location inside this engine at different scale altitudes. Discussions on observed limitation for engine icing scaling application from this test campaign and needed improvement are provided. A scaling test procedure to help identify potential ICI risk conditions and possible ice accretion locations of a new turbofan engine is evaluated in PSL.

Engine Testing in PSL

Ice-Crystal Icing Accretion Studies at the NASA Propulsion Systems Laboratory

This paper describes an ice-crystal icing experiment conducted at the NASA Propulsion System Laboratory during June 2018. This test produced ice shape data on an airfoil for different test conditions similar to those inside the compressor region of a turbo-fan jet engine. Mixed-phase icing conditions were generated by partially freezing out a water spray using the relative humidity of flow as the primary parameter to control freeze-out. The paper presents the ice shape data and associated conditions which include pressure, velocity, temperature, humidity, total water content, melt ratio, and particle size distribution. The test featured a new instrument traversing system which allowed surveys of the flow and cloud. The purpose of this work was to provide experimental ice shape data and associated conditions to help develop and validate ice-crystal icing accretion models. The results support previous experimental observations of a minimum melt-ratio threshold for accretion to occur as well as the existence of a plateau region where the icing severity is high for a range of melt ratios. However, a maximum limit for melt ratio, which is suggested in the ice crystal icing literature, was not observed perhaps complicated by the potential for some supercooling of the water at these conditions.

Engine Icing

The Influence of SLD Drop Size Distributions on Ice Accretion in the NASA Icing Research Tunnel

An ice shape database has been created to document ice accretions on a 21-inch chord NACA0012 model and a 72-inch chord NACA 23012 airfoil model resulting from an exposure to a Supercooled Large Drop (SLD) icing cloud with a bimodal drop size distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each accretion. The icing conditions used in the test matrix were based upon previously used conditions on the same models but with an alternate approach to evaluation of drop distribution effects. Ice shapes resulting from the bimodal distribution as well as from equivalent monomodal drop size distributions were obtained and compared. Results indicate that the ice shapes resulting from the monomodal and bimodal drop size distributions had similar shapes, but the bimodal distributions had greater mass and volume measurements and icing limits that extended further back on the chord of the model.

Potapczuk, Mark G.

Evaluation of a Thermodynamic Ice Crystal Icing Model Using Experimental Ice Accretion Data

This paper presents the evaluation of a thermodynamic ice crystal icing model, previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. It has been proposed that there are two types of distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. In the former, ice accretion occurs where a freeze fraction (0 to 1) of melted ice crystals freezes on a surface, along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized by having strong adhesion to the surface. In the latter, icing occurs from accumulated unmelted ice on a surface, where a melt fraction (0 to 1) dictates the amount of unmelted impinged ice. This melt-dominated icing is characterized by weakly bonded surface adhesion. The experimentally observed ice growth rates suggest that only a small fraction of the impinging ice remains on the surface, implying a mass loss mechanism such as splash, runback, bounce, or erosion. This mass loss parameter must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface. This loss parameter, however, along with the freeze and melt fraction, are the only experimental parameters that are currently not measured directly. Using reported icing growth rates from published ice crystal icing experiments, a methodology is proposed to determine these unknown parameters. This work takes reported ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) in 2016 and tests NASA collaborated on with the National Research Council (NRC) of Canada in 2012 that examined the fundamental physics of ice crystal icing. Those research efforts sought to generate icing conditions representative of those that occur inside a jet engine when ingesting ice crystals. This paper presents the fundamental equations of the thermodynamic model, the methodology used to determine the aforementioned unknown icing parameters, and results from model evaluation using experimental data. In addition, this paper builds on the previously proposed model by adding a transient conduction term to explain ice growth behavior at the onset of experimental tests that was observed to be different from steady-state ice growth that occurred later in the test run.With the addition of this energy term, this becomes a quasi-steady model. A key finding from this work suggests that mass loss fractions can exceed 0.90 for steady ice growth periods. In addition, due to conductive heat fluxes when using a warmer-than-freezing airfoil, lower mass loss fraction values were calculated during the initial transient period.

Ice-Crystal Ice Accretion

Further Evaluation of Swept Wing Icing Scaling with Maximum Combined Cross Section Ice Shape Profiles

This study reports results from recent icing scaling tests conducted in the NASA Glenn Icing Research Tunnel (IRT) using recommended scaling methods that were developed for swept wing icing applications. The aluminum test articles used in this study were NACA 0012 Adjustable Sweep Airfoil Models. The reference and scale airfoils have a 36 inch and 15 inch chord measured normal to the leading edge, and a 60 inch and 24 inch span respectively. Reference tests were run with airspeeds of 100, 140.3 and 150 knot and MVD of 41.4, 44.4, 93 and 150 micron. All tests were set at 0⁰ angle of attack (AOA) and 30⁰ sweep angle. Results will be presented for stagnation freezing fractions ranging from 0.3 to 1.0. For non-dimensional reference & scale ice shape comparison, a new post-scanning ice shape digitization procedure was applied to extract 2-D "hand-tracing like" ice shape profiles, i.e. the Maximum Combined Cross Section or MCCS, at selected span-wise locations from the high fidelity 3-D scanned ice shapes obtained in the IRT. Preliminary assessment of ice shape comparison with MCCS showed that good scaling was achieved for the test conditions by using the recommended scaling methods developed for swept wing icing.

MCCS

Ice Crystal Icing Physics Study using a NACA 0012 Airfoil at the National Research Council of Canada's Research Altitude Test Facility

This paper presents results from a study of the fundamental physics of ice-crystal ice accretion using a NACA 0012 airfoil at the National Research Council of Canada (NRC) Research Altitude Test Facility in August 2017. These tests were a continuation of work which began in 2010 as part of a joint collaboration between NASA and NRC. The research seeks to generate icing conditions representative of those that occur inside a jet engine when ingesting ice crystals. In this test, an airfoil was exposed to mixed-phase icing conditions and the resulting ice accretions were recorded and analyzed. This paper details the specific objectives, procedures, and measurements which included the aero-thermal and cloud measurements. The objectives were built upon observations and hypothesis generated from several previous test campaigns regarding mixed-phase ice-crystal icing. The specific objectives included (A) ice accretions under different wet-bulb temperatures, (B) investigations of steady-state ice shapes previously reported in the literature, (C) total water content variations in search of a threshold for accretion, and (D) probe characterization related to measuring melt fraction which is important to characterize the mixed-phase condition. The resulting ice accretions and conditions leading to such accretions are intended to help extend NASA’s predictive ice-accretion codes to include conditions occurring in engine ice-crystal icing.

Icing

Ice Roughness and Thickness Evolution on a Business Jet Airfoil

Experiments were performed in the Icing Research Tunnel at NASA Glenn Research Center (GRC) to investigate the ice roughness and thickness evolution on a 152.4-cm (60-in.) chord business jet airfoil exposed to both Appendix C and Appendix O (SLD - Super-cooled Large Droplet) icing conditions. The resulting measurements demonstrate that the average non-dimensional roughness and the stagnation point thickness scalings are similar to those demonstrated on symmetric wings. However, the surface variations of roughness and thickness exhibit significant differences from those observed on symmetric airfoils. The source of the roughness and thickness differences is the result of surface pressure, velocity and temperature distribution differences from the suction to the pressure sides of the airfoil. LEWICE (LEWis ICE accretion program - software developed at NASA Lewis Research Center - former name of the GRC) simulations are used to further investigate the influences of local collection efficiency and the local freezing fraction on the resulting ice roughness and thickness spatial variations.

icing codes

Ice Accretion Roughness Measurements and Modeling

Roughness on aircraft ice accretions is very important to the overall ice accretion process and to the resulting degradation in aircraft aerodynamic performance. Roughness enhances the local convection leading to more rapid ice accumulation rates, and roughness generates local flow perturbations that lead to higher skin friction. This paper presents 1) a review of the developments in ice shape three-dimensional laser scanning developed at NASA Glenn, 2) a review of the approach of McClain and Kreeger employed to characterize ice roughness evolution on an airfoil surface, and 3) a review of the experimental efforts that have been performed over the last five years to characterize, scale, and model ice roughness evolution physics.

SLD

Numerical Analysis of Mixed-Phase Icing Cloud Simulations in the NASA Propulsion Systems Laboratory

This presentation describes the development of a numerical model that couples the thermal interaction between ice particles, water droplets, and the flowing gas of an icing wind tunnel for simulation of NASA Glenn Research Centers Propulsion Systems Laboratory (PSL). The ultimate goal of the model is to better understand the complex interactions between the test parameters and have greater confidence in the conditions at the test section of the PSL tunnel. The model attempts to explain the observed changes in test conditions by coupling the conservation of mass and energy equations for both the cloud particles and flowing gas mass. Model predictions were compared to measurements taken during May 2015 testing at PSL, where test conditions varied gas temperature, pressure, velocity and humidity levels, as well as the cloud total water content, particle initial temperature, and particle size distribution.

NASA PSL

Numerical Analysis of Mixed-Phase Icing Cloud Simulations in the NASA Propulsion Systems Laboratory

This paper describes the development of a numerical model that couples the thermal interaction between ice particles, water droplets, and the flowing gas of an icing wind tunnel for simulation of NASA Glenn Research Centers Propulsion Systems Laboratory (PSL). The ultimate goal of the model is to better understand the complex interactions between the test parameters and have greater confidence in the conditions at the test section of the PSL tunnel. The model attempts to explain the observed changes in test conditions by coupling the conservation of mass and energy equations for both the cloud particles and flowing gas mass. Model predictions were compared to measurements taken during May 2015 testing at PSL, where test conditions varied gas temperature, pressure, velocity and humidity levels, as well as the cloud total water content, particle initial temperature, and particle size distribution.

NASA PSL

Preliminary Evaluation of Altitude Scaling for Turbofan Engine Ice Crystal Icing

Preliminary evaluation of altitude scaling for turbofan engine ice crystal icing simulation was conducted during the 2015 LF11 engine icing test campaign in PSL.The results showed that a simplified approach for altitude scaling to simulate the key reference engine ice growth feature and associated icing effects to the engine is possible. But special considerations are needed to address the facility operation limitation for lower altitude engine icing simulation.

altitude scaling

Comparisons of Mixed-Phase Icing Cloud Simulations with Experiments Conducted at the NASA Propulsion Systems Laboratory

This paper builds on previous work that compares numerical simulations of mixed-phase icing clouds with experimental data. The model couples the thermal interaction between ice particles and water droplets of the icing cloud with the flowing air of an icing wind tunnel for simulation of NASA Glenn Research Centers (GRC) Propulsion Systems Laboratory (PSL). Measurements were taken during the Fundamentals of Ice Crystal Icing Physics Tests at the PSL tunnel in March 2016. The tests simulated ice-crystal and mixed-phase icing that relate to ice accretions within turbofan engines.

Bartkus, Tadas

Ice Roughness and Thickness Evolution on a Swept NACA 0012 Airfoil

Several recent studies have been performed in the Icing Research Tunnel (IRT) at NASA Glenn Research Center focusing on the evolution, spatial variations, and proper scaling of ice roughness on airfoils without sweep exposed to icing conditions employed in classical roughness studies. For this study, experiments were performed in the IRT to investigate the ice roughness and thickness evolution on a 91.44-cm (36-in.) chord NACA 0012 airfoil, swept at 30-deg with 0deg angle of attack, and exposed to both Appendix C and Appendix O (SLD) icing conditions. The ice accretion event times used in the study were less than the time required to form substantially three-dimensional structures, such as scallops, on the airfoil surface. Following each ice accretion event, the iced airfoils were scanned using a ROMER Absolute Arm laser-scanning system. The resulting point clouds were then analyzed using the self-organizing map approach of McClain and Kreeger to determine the spatial roughness variations along the surfaces of the iced airfoils. The resulting measurements demonstrate linearly increasing roughness and thickness parameters with ice accretion time. Further, when compared to dimensionless or scaled results from unswept airfoil investigations, the results of this investigation indicate that the mechanisms for early stage roughness and thickness formation on swept wings are similar to those for unswept wings.

glaze icing

Preliminary Evaluation of Altitude Scaling for Turbofan Engine Ice Crystal Icing

Preliminary evaluation of altitude scaling for turbofan engine ice crystal icing simulation was conducted during the 2015 LF11 engine icing test campaign in PSL.The results showed that a simplified approach for altitude scaling to simulate the key reference engine ice growth feature and associated icing effects to the engine is possible. But special considerations are needed to address the facility operation limitation for lower altitude engine icing simulation.

engine ice crystal icing

Comparisons of Mixed-Phase Icing Cloud Simulations with Experiments Conducted at the NASA Propulsion Systems Laboratory

This paper builds on previous work that compares numerical simulations of mixed-phase icing clouds with experimental data. The model couples the thermal interaction between ice particles and water droplets of the icing cloud with the flowing air of an icing wind tunnel for simulation of NASA Glenn Research Centers (GRC) Propulsion Systems Laboratory (PSL). Measurements were taken during the Fundamentals of Ice Crystal Icing Physics Tests at the PSL tunnel in March 2016. The tests simulated ice-crystal and mixed-phase icing that relate to ice accretions within turbofan engines. Experimentally measured air temperature, humidity, total water content, liquid and ice water content, as well as cloud particle size, are compared with model predictions. The model showed good trend agreement with experimentally measured values, but often over-predicted aero-thermodynamic changes. This discrepancy is likely attributed to radial variations that this one-dimensional model does not address. One of the key findings of this work is that greater aero-thermodynamic changes occur when humidity conditions are low. In addition a range of mixed-phase clouds can be achieved by varying only the tunnel humidity conditions, but the range of humidities to generate a mixed-phase cloud becomes smaller when clouds are composed of smaller particles. In general, the model predicted melt fraction well, in particular with clouds composed of larger particle sizes.

Thermal model

Bimodal SLD Ice Accretion on a NACA 0012 Airfoil Model

An initial ice shape database has been created to document ice accretions on a 21-inch chord NACA 0012 airfoil model resulting from an exposure to a Supercooled Large Droplet (SLD) icing cloud with a bimodal droplet distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each accretion. The icing conditions used in the test matrix were based upon previously measured ice shapes on the same model to connect the current database to previously measured information. Ice shapes resulting from the bimodal distribution as well as from equivalent standard droplet distributions were obtained and compared. Results indicate that the ice shapes resulting from the bimodal droplet distributions had higher mass and volume values than their standard distribution equivalents as well as having icing limits that extended further back on the chord of the model.

aircraft icing

Plans and Preliminary Results of Fundamental Studies of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory

This paper describes plans and preliminary results for using the NASA Propulsion Systems Lab (PSL) to experimentally study the fundamental physics of ice-crystal ice accretion. NASA is evaluating whether this facility, in addition to full-engine and motor-driven-rig tests, can be used for more fundamental ice-accretion studies that simulate the different mixed-phase icing conditions along the core flow passage of a turbo-fan engine compressor. The data from such fundamental accretion tests will be used to help develop and validate models of the accretion process. This paper presents data from some preliminary testing performed in May 2015 which examined how a mixed-phase cloud could be generated at PSL using evaporative cooling in a warmer-than-freezing environment.

Ice-Crystal Ice Accretion