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

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 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)

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 collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Collection efficiencies were generated for three spherical particle sizes, 100, 20, and 5 micron in diameter, using the codes LEWICE3D and LEWICE2D. The free stream Mach number was 0.27, representing a velocity of approximately 86 ms. It was observed that a reduction in velocity of about 15-20 occurred as the flow entered the shroud of the probe.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 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 effect of the protective shroud

A Model to Assess the Risk of Ice Accretion Due to Ice Crystal Ingestion in a Turbofan Engine and its Effects on Performance

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that were attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was one or more of the following anomalies: degraded engine performance, engine roll back, compressor surge and stall, and flameout of the combustor. The main focus of this research is the development of a computational tool that can estimate whether there is a risk of ice accretion by tracking key parameters through the compression system blade rows at all engine operating points within the flight trajectory. The tool has an engine system thermodynamic cycle code, coupled with a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Assumptions are made to predict the complex physics involved in engine icing. Specifically, the code does not directly estimate ice accretion and does not have models for particle breakup or erosion. Two key parameters have been suggested as conditions that must be met at the same location for ice accretion to occur: the local wet-bulb temperature to be near freezing or below and the local melt ratio must be above 10%. These parameters were deduced from analyzing laboratory icing test data and are the criteria used to predict the possibility of ice accretion within an engine including the specific blade row where it could occur. Once the possibility of accretion is determined from these parameters, the degree of blockage due to ice accretion on the local stator vane can be estimated from an empirical model of ice growth rate and time spent at that operating point in the flight trajectory. The computational tool can be used to assess specific turbine engines to their susceptibility to ice accretion in an ice crystal environment

Transient Catalytic Combustor Model With Detailed Gas and Surface Chemistry

In this work, we numerically investigate the transient combustion of a premixed gas mixture in a narrow, perfectly-insulated, catalytic channel which can represent an interior channel of a catalytic monolith. The model assumes a quasi-steady gas-phase and a transient, thermally thin solid phase. The gas phase is one-dimensional, but it does account for heat and mass transfer in a direction perpendicular to the flow via appropriate heat and mass transfer coefficients. The model neglects axial conduction in both the gas and in the solid. The model includes both detailed gas-phase reactions and catalytic surface reactions. The reactants modeled so far include lean mixtures of dry CO and CO/H2 mixtures, with pure oxygen as the oxidizer. The results include transient computations of light-off and system response to inlet condition variations. In some cases, the model predicts two different steady-state solutions depending on whether the channel is initially hot or cold. Additionally, the model suggests that the catalytic ignition of CO/O2 mixtures is extremely sensitive to small variations of inlet equivalence ratios and parts per million levels of H2

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

Engine Icing Performance loss: rollback, surge, flameout, and even internal engine damagePartial melting and refreeze of ice inside engine core (Mason et al., 2006). Ingestion of ice crystals and aggregates, mixed-phase droplets, or supercooled liquid dropletsNeed to better understand the conditions and properties that lead to engine icing.Simulation and analysis (physical and computational, and modeling)Test facilities (PSL, NRC, ...). Thermal and computational models and analysisProbesMultiple probes (aerothermal probes and ice cloud characterization probes and techniques). Total temperatureTraditional total temperature probes (vented forward facing)Heated total temperature probes (Goodyear). Rearward facing (developmental). Total temperature relevance. Thermal interaction between the icing cloud and air flow impinging particles contribute to kinetic heating effect (Gent et al., 2000). Measurement considerations Temperature sensor accuracy. Incomplete recovery of total temperature. Thermal surfaces (sources and sinks). Flow effects (viscous losses). Debris contamination, including icing and ice ingestion

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 NASAs predictive ice-accretion codes to include conditions occurring in engine ice-crystal icing