Search NASA⌕ Search

SEARCH · Search NASA

Results for “AIRCRAFT ICING”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Aircraft Icing Severity Evaluation

Aircraft icing refers to the ice buildup on the surface of an aircraft flying in icing conditions. The ice accretion on the aircraft alters the original aerodynamic configuration and degrades the aerodynamic performances and may lead to unsafe flight conditions. Evaluating the flow structure, icing mechanism and consequences is of great importance to the development of an anti/deicing technique. Studies have shown computational fluid dynamics (CFD) and machine learning (ML) to be effective in predicting the ice shape and icing severity under different flight conditions. CFD solves a set of partial differential equations to obtain the air flow fields, water droplets trajectories and ice shape. ML is a branch of artificial intelligence and, based on the data, the self-improved computer algorithms can be effective in finding the nonlinear mapping relationship between the input flight conditions and the output aircraft icing severity features.

Li, Sibo (ORCID:000000021705844X)↗
DOE OSTI2020

Fast Evaluation of Aircraft Icing Severity Using Machine Learning Based on XGBoost

Aircraft icing represents a serious hazard in aviation which has caused a number of fatal accidents over the years. In addition, it can lead to substantial increase in drag and weight, thus reducing the aerodynamics performance of the airplane. The process of ice accretion on a solid surface is a complex interaction of aerodynamic and environmental variables. The complex relationship makes machine learning-based methods an attractive alternative to traditional numerical simulation-based approaches. In this study, we introduce a purely data-driven approach to find the complex pattern between different flight conditions and aircraft icing severity prediction. The supervised learning algorithm Extreme Gradient Boosting (XGBoost) is applied to establish the prediction framework which makes prediction based on any set of observations. The input flight conditions for the proposed prediction framework are liquid water content, droplet diameter and exposure time. The proposed approach is demonstrated in three cases: maximum ice thickness prediction, icing area prediction and icing severity level evaluation. Performance comparison studies and error analysis are also conducted to verify the effectiveness and performance of the proposed method. Results show that the proposed method has reasonable capability in evaluating aircraft icing severity.

Li, Sibo (ORCID:000000021705844X)↗
DOE OSTI2025

A rigorous framework for an improved Messinger/Myers model of ice accretion under conditions of variable property and unsteady aircraft icing

We analyse the Messinger/Myers model by critically evaluating simplifying assumptions through a rigorous formulation of the rime ice accretion process. We explore the effects of both constant and variable ice density and thermal conductivity, along with the effects of sublimation from the ice surface. The effects of key factors such as droplet impact rate, ambient temperature relative to the freezing temperature and the temperature difference between the ambient air and the airfoil surface are examined. Under these varying conditions, the present rigorous formulation is used to assess the significance of unsteady effects, variable ice properties and sublimation. We observe that the Myers model performs remarkably well in certain icing situations and analyse the reasons for this strong performance. We also show that partially relaxing the model’s assumptions can lead to poorer performance. The Myers model can lead to overprediction of ice surface temperature and correspondingly underprediction of transition time under conditions of relatively weak sublimation and surface cooling. A modified Myers model is presented, which can be used to recover near-perfect results under widely varying icing conditions of relevance. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗
DOE OSTI2016

Effect of drop size on the impact thermodynamics for supercooled large droplet in aircraft icing

Supercooled large droplet (SLD), which can cause abnormal icing, is a well-known issue in aerospace engineering. Although efforts have been exerted to understand large droplet impact dynamics and the supercooled feature in the film/substrate interface, respectively, the thermodynamic effect during the SLD impact process has not received sufficient attention. This work conducts experimental studies to determine the effects of drop size on the thermodynamics for supercooled large droplet impingement. Through phenomenological reproduction, the rapid-freezing characteristics are observed in diameters of 400, 800, and 1300 μm. The experimental analysis provides information on the maximum spreading rate and the shrinkage rate of the drop, the supercooled diffusive rate, and the freezing time. A physical explanation of this unsteady heat transfer process is proposed theoretically, which indicates that the drop size is a critical factor influencing the supercooled heat exchange and effective heat transfer duration between the film/substrate interface. On the basis of the present experimental data and theoretical analysis, an impinging heating model is developed and applied to typical SLD cases. The model behaves as anticipated, which underlines the wide applicability to SLD icing problems in related fields.

42 ENGINEERING↗
DOE OSTI2022

Numerical Modeling of Atmospheric Rime Ice Accretion on an Airfoil Using an Eulerian Approach

This paper presents an approach to numerically simulate the inherently unsteady rime ice accretion problem on a two-dimensional airfoil and elucidate the associated variations under different icing conditions. The airflow field and the water impingement on the airfoil are obtained based on an Eulerian two-phase model. A dynamic mesh strategy is employed to unsteadily account for the changes in the ice profile and its impact on the air and droplet flow by continuously reconstructing the computational grid at each time-step through smoothing and layering mechanisms. All main icing modules including the airflow field, droplet trajectory, icing thickness profile, and mesh management are fully coupled within the same computational framework without resorting to any external tools. Classical icing theory is employed to model the rime ice roughness, and it is assumed that the ice accretes in a direction normal to the airfoil surface. The governing Reynolds-averaged Navier–Stokes (RANS) conservation equations along with the energy and continuity equations are solved to produce the velocity and temperature fields. A convective film heat transfer coefficient is computed based on the surface heat flux and a recovery temperature which takes into account the dissipative heat release in the boundary layer in the vicinity of the airfoil surface. With the implemented strategy and calculating the convective heat transfer coefficient, the water film thickness is also calculated along with the ice shape. Furthermore, the model is validated by comparing the local collection efficiency distribution and ice shape with experimental data, and the results show that the implemented approach provides acceptable predictions of ice accretion profiles and rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2019

Modeling of Ice Accretion over Aircraft Wings Using a Compressible OpenFOAM Solver

A method to simulate ice accretion on an aircraft wing using a three-dimensional compressible Navier-Stokes solver, a Eulerian droplet flow field model, a mesh morphing model, and a thermodynamic model, is presented in this paper. The above models are combined together into one solver and implemented in OpenFOAM. Two-way coupling is achieved between airflow field calculation and ice simulation. The density-based solver rhoEnergyFoam is used to calculate the airflow field. The roughness wall function is proposed to simulate the roughness effect caused by ice accretion. For droplet flow field calculation, the Eulerian model is applied and the permeable wall boundary condition is used on the wing to simulate the droplet impingement. The icing thermodynamic model is built based on the Messinger model. The mesh morphing model adjusts the wing’s shape every time step based on the amount of accreted ice so that the airflow field is updated during the simulation. The effect of the ice accretion on the airflow is studied by comparing the aerodynamic performance—with and without ice. The ice accretion on the ONERA M6 wing model under a specific condition has been simulated to validate the solver’s performance and investigate the effect of the accreted ice on the aerodynamic performance.

Li, Sibo↗
DOE OSTI2017

Modeling of Aircraft Measurements of Ice Crystal Concentration in the Arctic and a Parameterization for Mixed-Phase Cloud

Abstract In this study, two parameterizations of ice nucleation rate on dust particles are used in a parcel model to simulate aircraft measurements of ice crystal number concentration Ni in the Arctic. The parcel model has detailed microphysics for droplet and ice nucleation, growth, and evaporation with prescribed vertical air velocities. Three dynamic regimes are considered, including large-scale ascent, cloud-top generating cells, and their combination. With observed meteorological conditions and aerosol concentrations, the parcel model predicts the number concentrations of size-resolved ice crystals, which may be compared to aircraft measurements. Model results show rapid changes with height/time in relative humidity, Ni, and thermodynamic phase partitioning, which is not resolved in current climate and weather forecasting models. Parameterizations for ice number and nucleation rate in mixed-phase stratus clouds are thus developed based on the parcel model results to represent the time-integrated effect of some microphysical processes in large-scale models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Dynamical conditions of ice supersaturation and ice nucleation in convective systems: A comparative analysis between in situ aircraft observations and WRF simulations: Ice Supersaturation in Convective System

Occurrence frequency and dynamical conditions of ice supersaturation (ISS, where relative humidity with respect to ice (RHi) > 100%) are examined in the upper troposphere around convective activity. Comparisons are conducted between in situ airborne observations and the Weather Research and Forecasting model simulations using four double-moment microphysical schemes at temperatures ≤-40°C. All four schemes capture both clear-sky and in-cloud ISS conditions. However, the clear-sky (in-cloud) ISS conditions are completely (significantly) limited to the RHi thresholds of the Cooper parameterization. In all of the simulations, ISS occurrence frequencies are higher by ~3–4 orders of magnitude at higher updraft speeds (>1 m s -1 ) than those at the lower updraft speeds when ice water content (IWC) > 0.01 g m -3 , while observations show smaller differences up to ~1–2 orders of magnitude. The simulated ISS also occurs less frequently at weaker updrafts and downdrafts than observed. These results indicate that the simulations have a greater dependence on stronger updrafts to maintain/generate ISS at higher IWC. At lower IWC (≤0.01 g m -3 ), simulations unexpectedly show lower ISS frequencies at stronger updrafts. Overall, the Thompson aerosol-aware scheme has the closest magnitudes and frequencies of ISS >20% to the observations, and the modified Morrison has the closest correlations between ISS frequencies and vertical velocity at higher IWC and number density. The Cooper parameterization often generates excessive ice crystals and therefore suppresses the frequency and magnitude of ISS, indicating that it should be initiated at higher ISS (e.g., ≥25%).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Physicochemical characterization of free troposphere and marine boundary layer ice-nucleating particles collected by aircraft in the eastern North Atlantic

Abstract. Atmospheric ice nucleation impacts the hydrological cycle and climate by modifying the radiative properties of clouds. To improve our predictive understanding of ice formation, ambient ice-nucleating particles (INPs) need to be collected and characterized. Measurements of INPs at lower latitudes in a remote marine region are scarce. The Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign, in the region of the Azores islands, provided the opportunity to collect particles in the marine boundary layer (MBL) and free troposphere (FT) by aircraft during the campaign's summer and winter intensive operation period. The particle population in samples collected was examined by scanning transmission X-ray microscopy with near-edge X-ray absorption fine structure spectroscopy. The identified INPs were analyzed by scanning electron microscopy with energy-dispersive X-ray analysis. We observed differences in the particle population characteristics in terms of particle diversity, mixing state, and organic volume fraction between seasons, mostly due to dry intrusion events during winter, as well as between the sampling locations of the MBL and FT. These differences are also reflected in the temperature and humidity conditions under which water uptake, immersion freezing (IMF), and deposition ice nucleation (DIN) proceed. Identified INPs reflect typical particle types within the particle population on the samples and include sea salt, sea salt with sulfates, and mineral dust, all associated with organic matter, as well as carbonaceous particles. IMF and DIN kinetics are analyzed with respect to heterogeneous ice nucleation rate coefficients, Jhet, and ice nucleation active site density, ns, as a function of the water criterion Δaw. DIN is also analyzed in terms of contact angles following classical nucleation theory. Derived MBL IMF kinetics agree with previous ACE-ENA ground-site INP measurements. FT particle samples show greater ice nucleation propensity compared to MBL particle samples. This study emphasizes that the types of INPs can vary seasonally and with altitude depending on sampling location, thereby showing different ice nucleation propensities, which is crucial information when representing mixed-phase cloud and cirrus cloud microphysics in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Columnar Vertical Profile (CVP) Methodology for Validating Polarimetric Radar Retrievals in Ice Using In Situ Aircraft Measurements

Abstract A novel way to process polarimetric radar data collected via plan position indicator (PPI) scans and display those data in a time–height format is introduced. The columnar vertical profile (CVP) methodology uses radar data collected via multiple elevation scans, limited to data within a set region in range and azimuth relative to the radar, to create vertical profiles of polarimetric radar data representative of that limited region in space. This technique is compared to others existing in the literature, and various applications are discussed. Polarimetric ice microphysical retrievals are performed on CVPs created within the stratiform rain region of two mesoscale convective systems sampled during two field campaigns, where CVPs follow the track of research aircraft. Aircraft in situ data are collocated to microphysical retrieval data, and the accuracy of these retrievals is tested against other retrieval techniques in the literature.

Engineering↗
DOE OSTI2023

A rotorcraft in-flight ice detection framework using computational aeroacoustics and Bayesian neural networks

Abstract This work develops a novel ice detection framework specifically suitable for rotorcraft using computational aeroacoustics and Bayesian neural networks. In an offline phase of the work, the acoustic signature of glaze and rime ice shapes on an oscillating wing are computed. In addition, the aerodynamic performance indicators corresponding to the ice shapes are also monitored. These performance indicators include the lift, drag, and moment coefficients. A Bayesian neural network is subsequently trained using projected Stein variational gradient descent to create a mapping from the acoustic signature generated by the iced wings to predict their performance indicators along with quantified uncertainty that is highly important for time- and safety-critical decision-making scenarios. While the training is carried out fully offline, usage of the Bayesian neural network to make predictions can be conducted rapidly online allowing for an ice detection system that can be used in real time and in-flight.

42 ENGINEERING↗
DOE OSTI2015

Investigation of ice cloud microphysical properties of DCSs using aircraft in situ measurements during MC3E over the ARM SGP site

In this study, six deep convective systems (DCSs) with a total of 5589 five-second samples and a range of temperatures from -41°C to 0°C during the Midlatitude Continental Convective Clouds Experiment (MC3E) were selected to investigate the ice cloud microphysical properties of DCSs over the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site. The ice cloud measurements of the DCS cases were made by the University of North Dakota Citation II research aircraft, and the ice cloud properties were derived through the following processes. First, the instances of supercooled liquid water in the ice-dominated cloud layers of DCSs have been eliminated using multisensor detection, including the Rosemount Icing Detector, King and Cloud Droplet Probes, as well as 2DC and Cloud Imaging Probe images. Then the Nevzorov-measured ice water contents (IWCs) at maximum diameter D max < 4000 µm are used as the best estimation to determine a new mass-dimensional relationship. Finally, the newly derived mass-dimensional relationship (a = 0.00365, b = 2.1) has been applied to a full spectrum of particle size distributions (PSDs, 120–30,000 µm) constructed from both 2DC and High-Volume Precipitation Spectrometer measurements to calculate the best-estimated IWCs of DCSs during MC3E. The averages of the total number concentrations (N t ), median mass diameter (D m ), maximum diameter (D max ), and IWC from six selected cases are 0.035 cm -3 , 1666 µm, 8841 µm, and 0.45 g m -3 , respectively. The gamma-type-size distributions are then generated matching the observed PSDs (120–30,000 µm), and the fitted gamma parameters are compared with the observed PSDs through multimoment assessments including first moment (D m ), third moment (IWC), and sixth moment (equivalent radar reflectivity, Z e ). Lastly, for application of observed PSDs to the remote sensing community, a series of empirical relationships between fitted parameters and Z e values has been derived, and the bullet rosette ice crystal backscattering relationship has been suggested for ground-based remote sensing.

58 GEOSCIENCES↗
DOE OSTI2023

Ice Supersaturation Variability in Cirrus Clouds: Role of Vertical Wind Speeds and Deposition Coefficients

Abstract Aircraft measurements reveal ice supersaturation statistics in cirrus (ISSs) with broad maxima around ice saturation and pronounced variance. In this study, processes shaping ISSs in midlatitude and tropical upper tropospheric conditions are systematically investigated. Water vapor deposition and sublimation of size‐resolved ice crystal populations are simulated in an air parcel framework. Mesoscale temperature fluctuations (MTFs) due to gravity waves force the temporal evolution of supersaturation. Various levels of background wave forcing and cirrus thickness are distinguished in stochastic ensemble simulations. Kinetic limitations to ice mass growth are brought about by supersaturation‐dependent deposition coefficients that represent efficient and inefficient growth modes as a function of ice crystal size. The simulations identify a wide range of deposition coefficients in cirrus, but most values stay above 0.01 such that kinetic limitations to water uptake remain moderate. Supersaturation quenching times are long, typically 0.5–2 hr. The wave forcing thus causes a remarkably large variability in ISSs and cirrus microphysical properties except in the thickest cirrus, producing ensemble‐mean ISSs in line with in‐situ measurements. ISS variance is controlled by MTFs and increases with decreasing cirrus integral radii. In comparison, the impact of ice crystal growth rates on ISSs is small. These results contribute to efforts directed at identifying and solving issues associated with ice‐supersaturated areas and non‐equilibrium cirrus physics in global models.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Stokes-dependent droplet collection efficiency on a NACA 0012 airfoil from droplet-informed simulations with statistical overloading

Accurate modelling of ice accretion on aircraft wings requires analysing droplet impingement on the surface to optimize the design of ice-protection systems. We perform Euler–Lagrange simulations of a droplet-laden flow impinging on a NACA 0012 airfoil. Our study includes water droplets with eight discrete sizes ranging from 1 to 160 microns. We vary the free-stream velocity of the incoming airflow in the range 60 ≤ U ≤ 240 m s −1 and the chord length of the airfoil in the range 0.5 ≤ c ≤ 2 m. Due to the dilute nature of supercooled clouds, one-way coupling is used in the simulations. The effects of droplet breakup and collision are also neglected. To reduce the computational cost, we employ statistical overloading of droplets, allowing us to simulate millions of impinging droplets in a time span on the order of milliseconds. Our results show that the droplet collection efficiency, which measures the likelihood of droplet impingement on the airfoil surface, increases with droplet size and free-stream velocity but decreases with airfoil size. We demonstrate that collection efficiency, impingement velocity and impingement angle are primarily dictated by a single non-dimensional parameter, the droplet Stokes number. We also identify a critical stagnation-streamline Stokes number below which impingements do not occur and use it to estimate the minimum droplet size for impingement. In addition, we observe droplet behaviour to become Stokes number independent at large values of the Stokes number. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗
DOE OSTI2024

Aerosol-Cloud Interactions From Aviation Soot Emissions

Current models estimate global aviation contributes approximately 5% to the total anthropogenic climate forcing, with aerosol-cloud interactions having the greatest effect. However, radiative forcing estimates from aviation aerosol-cloud interactions remain undetermined. There is an expected significant increase in aircraft emissions with aviation demand expected to rise by over 4% per year. Soot may play an important role in the ice nucleation of aircraft-induced cirrus formation due to a high emission rate, but the ice nucleating properties are poorly constrained. Understanding the microphysical processes leading to atmospheric ice crystal formation is crucial for the reliable parameterization of aerosol-cloud interactions in climate models due to their impact on precipitation and cloud radiative properties. Ice nucleation of aircraft-emitted soot is potentially affected by particle morphology with condensation of supercooled water occurring in pores followed by ice nucleation. However, soot has heterogeneous properties and undergoes atmospheric aging and oxidation that could change surface properties and contribute to complex ice nucleation processes. Further, this review synthesizes current knowledge of ice nucleation catalyzed by aviation in the cirrus regime and its effects on global radiative forcing. Further research is required to determine the ice nucleation and microphysical processes of cirrus cloud formation from aviation emissions in both controlled laboratory and field investigations to inform models for more accurate climate predictions and to provide efficient mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Relationships between summertime surface albedo and melt pond fraction in the central Arctic Ocean: The aggregate scale of albedo obtained on the MOSAiC floe

As part of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC), the HELiX uncrewed aircraft system (UAS) was deployed over the sea ice in the central Arctic Ocean during summer 2020. Albedo measurements were obtained with stabilized pyranometers, and melt pond fraction was calculated from orthomosaic imagery from a surface-imaging multispectral camera. This study analyzed HELiX flight data to provide insights on the temporal and spatial evolution of albedo and melt pond fraction of the MOSAiC floe during the melt season as it drifted south through Fram Strait. The surface albedo distributions showed peak values changing from high albedo (0.55–0.6) to lower values (0.3) as the season advanced. Inspired by methods developed for satellite data, an algorithm was established to retrieve melt pond fraction from the orthomosaic images. We demonstrate that the near-surface observations of melt pond fraction were highly dependent on sample area, offering insight into the influence of subgrid scale features and spatial heterogeneity in satellite observations. Vertical observations conducted with the HELiX were used to quantify the influence of melt pond scales on observed surface albedo as a function of sensor footprint. These scaling results were used to link surface-based measurements collected during MOSAiC to broader-scale satellite data to investigate the influence of surface features on observed albedo. Albedo values blend underlying features within the sensor footprint, as determined by the melt pond size and concentration. This study framed the downscaling (upscaling) problem related to the airborne (surface) observations of surface albedo across a variety of spatial scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Liquid–Ice Mass Partition in Tropical Maritime Convective Clouds

The liquid–ice mass partitioning in tropical maritime convective clouds is studied in this paper using data collected by the National Center for Atmospheric Research C-130 research aircraft during the Ice in Clouds Experiment–Tropical project. The clouds investigated by the C-130 in this study generally contained weak to moderate updrafts. The liquid water content (LWC) is calculated using a combination of hot-wire and imaging probes. The total condensed water content (CWC) is measured by a counterflow virtual impactor. The ice water content (IWC) is calculated as CWC minus LWC. Taking into account potential significant measurement uncertainties, the liquid fraction [i.e., LWC/(LWC + IWC)] between 0° and -15°C appears to decrease by a factor of about 3 in updrafts near (<500 m) cloud top and a factor of 2 in updrafts far below (>500 m) cloud top. The decrease in liquid fraction as a function of temperature is also correlated with cloud life cycle. In dissipating clouds, ice dominates in all temperature ranges. A comparison between this study and two parameterizations shows that at different geographic locations the liquid fraction in convective clouds differs. Because of the sampling bias and the limitations of instruments, more measurements, especially with more advanced instruments, are needed in the future.

54 ENVIRONMENTAL SCIENCES↗