Engineering topics

Bartkus, Tadas P.

Publications and source records attributed to Bartkus, Tadas P..

11 records

NASA NTRS2019

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↗

NASA NTRS2019

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.

Agui, Juan H.↗

NASA NTRS2018

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↗

NASA NTRS2018

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↗

NASA NTRS2018

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.

Agui, Juan H.↗

NASA NTRS2018

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.

ice crystal↗

NASA NTRS2017

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↗

NASA NTRS2017

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 -20˚C. 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.

in-flight icing↗

NASA NTRS2017

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↗

NASA NTRS2007

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.

Bartkus, Tadas P.↗

NASA NTRS2007

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.

Bartkus, Tadas P.↗