Experimental Investigation of Supercooled Water Droplet Breakup near the Leading Edge of an Airfoil

This paper presents the results of an experimental study on supercooled droplet deformation and breakup near the leading edge of an airfoil. The results are compared to prior room-temperature droplet deformation results to explore the effects of droplet supercooling. The experiments were conducted in the Adverse Environment Rotor Test Stand at The Pennsylvania State University. An airfoil model placed at the end of the rotor blades mounted onto the hub in the Adverse Environment Rotor Test Stand chamber was moved at speeds ranging between 50 and 80 ms. The temperature of the chamber was 20C. A monotonic droplet generator was used to produce droplets that fell perpendicular to the airfoil path. High-speed imaging was employed to observe the interaction between the droplets and the airfoil. Cases with equal slip and initial velocity were selected for the two environmental conditions. The airfoil velocity was 60 ms, and the slip velocity for both sets of data was 40 ms. The deformation of the weakly supercooled and warm droplets did not present different trends. The similar behavior for both conditions indicates that water supercooling has no effect on particle deformation for the range of supercooling of the droplets tested and the selected impact velocity

Total Temperature Measurements in Icing Cloud Flows Using a Rearward Facing Probe

This paper reports on temperature and humidity measurements from a series of ice-crystal icing tunnel experiments conducted in June 2018 at the Propulsion Systems Laboratory at the NASA Glenn Research Center. The tests were fundamental in nature and were aimed at investigating the icing processes on a two-dimensional NACA0012 airfoil subjected to artificially generated icing clouds. Prior to the tests on the airfoil, a suite of instruments, including total temperature and humidity probes, were used to characterize the thermodynamic flow and icing cloud conditions of the facility. Two different total temperature probes were used in these tests which included a custom designed rearward facing probe and a commercial self-heating total temperature probe. The rearward facing probe, the main total temperature probe, was designed to reduce and mitigate the contaminating effects of icing and ingestion of ice crystals and water droplets at the probe's inlet. The probe also serves as an air-sample inlet for a light absorption based humidity measurement. The paper includes a section which discusses total temperature and humidity measurement considerations, and another section which provides an analysis of the main probe's performance characteristics. A computational fluid dynamic model of the flow around the probe was also conducted to gain insight into the trajectory of the flow entering the probe inlet. The experiments included a series of tests in which the relative humidity of the facility flow was swept through with increasingly larger values. The data showed that the rearward facing probe can reasonably capture the flow's total temperature and humidity under mild to moderate icing conditions but produces anomalous results under more intense icing conditions. The experimental data was also compared to an in-house developed thermodynamic model which takes into account the interaction of the main flow with the icing cloud. Comparison to the thermodynamic model showed that the rearward facing probe measured the predicted trends

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)

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

Total Temperature Measurements using a Rearward Facing Probe in Supercooled Liquid Droplet and Ice Crystal Clouds

This paper presents an analysis of local total temperature and humidity experimental measurement taken in atmospheric ice cloud flows. The measurements were obtained in a series of tests in NASA's Propulsion Systems Laboratory. The probe used in the tests is referred to as the Rearward Facing Probe which was designed to mitigate the contamination effects of ice accretion and ingestion into the probe. The data provided important insights in the interaction of the ice cloud and the atmospheric flow. For the majority of the test runs, small temperature drops in the range of 0.6 to 2.8 C and up to 1.5 g/kg of water vapor rise were found as a result of the interaction. Under certain very low temperature or high TWC conditions, the interaction with the cloud produced a warming of the airflow. A thermal model based on evaporative and convective heat transfer mechanisms between the spray droplets and the airflow showed good agreement with the experimental data. Detailed analyses of the response of the probe under various flow, thermodynamic, and cloud conditions, are provided in the paper

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

Pyrolysis of Large Black Liquor Droplets

This paper presents the results of experiments involving the pyrolysis of large black liquor droplets in the NASA KC-135 reduced gravity aircraft. The reduced gravity environment facilitated the study of droplets up to 9 mm in diameter extending the results of previous studies to droplet sizes that are similar to those encountered in recovery boilers. Single black liquor droplets were rapidly inserted into a 923 K oven. The primary independent variables were the initial droplet diameter (0.5 mm to 9 mm), the black liquor solids content (66.12% - 72.9% by mass), and the ambient oxygen mole fraction (0.0 - 0.21). Video records of the experiments provided size and shape of the droplets as a function of time. The results show that the particle diameter at the end of the drying stage (D(sub DRY) ) increases linearly with the initial particle diameter (D(sub O)). The results further show that the ratio of the maximum swollen diameter (D(sub MAX)) to D(sub O) decreases with increasing D(sub O) for droplets with D(sub O) less than 4 mm. This ratio was independent of D(sub O) for droplets with D(sub O) greater than 4 mm. The particle is most spherical at the end of drying, and least spherical at maximum swollen size, regardless of initial sphericity and droplet size

Mechanism of Supercooled Water Droplet Breakup near the Leading Edge of an Airfoil

This work presents the results of an experimental study on supercooled droplet deformation and breakup near the leading edge of an airfoil. The results are compared to prior room temperature droplet deformation results to explore the effects of droplet supercooling. The experiments were conducted in the Adverse Environment Rotor Test Stand (AERTS) at The Pennsylvania State University. An airfoil model placed at the end of the rotor blades mounted onto the hub in the AERTS chamber was moved at speeds ranging between 50 and 80 m/sec. The temperature of the chamber was set at -20C. A monotonic droplet generator was used to produce droplets that fell from above, perpendicular to the path of the airfoil. The supercooled state of the droplets was determined by measurement of the temperature of the drops at various locations below the droplet generator exit. A temperature prediction code was also used to estimate the temperature of the droplets based on vertical velocity and the distance traveled by droplets from the droplet generator to the airfoil stagnation line. High speed imaging was employed to observe the interaction between the droplets and the airfoil. The high speed imaging provided droplet deformation information as the droplet approached the airfoil near the stagnation line. A tracking software program was used to measure the horizontal and vertical displacement of the droplet against time. It was demonstrated that to compare the effects of water supercooling on droplet deformation, the ratio of the slip velocity and the initial droplet velocity must be equal. A case with equal slip velocity to initial velocity ratios was selected for room temperature and supercooled droplet conditions. The airfoil velocity was 60 m/s and the slip velocity for both sets of data was 40 m/s. In these cases, the deformation of the weakly supercooled and warm droplets did not present different trends. The similar behavior for both environmental conditions indicates that water supercooling has no effect on particle deformation for the limited range of the weak supercooling of the droplets tested and the selected impact velocity. The assumption of a constant surface tension value was further supported by the equal trend of the Bond number obtained for supercooled and room temperature droplets

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

Recent Advances in the LEWICE Icing Model

This paper will describe two recent modifications to the Glenn ICE software. First, a capability for modeling ice crystals and mixed phase icing has been modified based on recent experimental data. Modifications have been made to the ice particle bouncing and erosion model. This capability has been added as part of a larger effort to model ice crystal ingestion in aircraft engines. Comparisons have been made to ice crystal ice accretions performed in the NRC Research Altitude Test Facility (RATFac). Second, modifications were made to the run back model based on data and observations from thermal scaling tests performed in the NRC Altitude Icing Tunnel