Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Modeling and Grid Generation of Iced Airfoils

SmaggIce Version 2.0 is a software toolkit for geometric modeling and grid generation for two-dimensional, singleand multi-element, clean and iced airfoils. A previous version of SmaggIce was described in Preparing and Analyzing Iced Airfoils, NASA Tech Briefs, Vol. 28, No. 8 (August 2004), page 32. To recapitulate: Ice shapes make it difficult to generate quality grids around airfoils, yet these grids are essential for predicting ice-induced complex flow. This software efficiently creates high-quality structured grids with tools that are uniquely tailored for various ice shapes. SmaggIce Version 2.0 significantly enhances the previous version primarily by adding the capability to generate grids for multi-element airfoils. This version of the software is an important step in streamlining the aeronautical analysis of ice airfoils using computational fluid dynamics (CFD) tools. The user may prepare the ice shape, define the flow domain, decompose it into blocks, generate grids, modify/divide/merge blocks, and control grid density and smoothness. All these steps may be performed efficiently even for the difficult glaze and rime ice shapes. Providing the means to generate highly controlled grids near rough ice, the software includes the creation of a wrap-around block (called the "viscous sublayer block"), which is a thin, C-type block around the wake line and iced airfoil. For multi-element airfoils, the software makes use of grids that wrap around and fill in the areas between the viscous sub-layer blocks for all elements that make up the airfoil. A scripting feature records the history of interactive steps, which can be edited and replayed later to produce other grids. Using this version of SmaggIce, ice shape handling and grid generation can become a practical engineering process, rather than a laborious research effort.

Vickerman, Mary B.↗
NASA NTRS2003

Predicting Accumulations of Ice on Aerodynamic Surfaces

LEWICE is a computer program that predicts the accumulation of ice on two-dimensional aerodynamic surfaces under conditions representative of the flight of an aircraft through an icing cloud. The software first calculates the airflow surrounding the body of interest, then uses the airflow to compute the trajectories of water droplets that impinge on the surface of the body. The droplet trajectories are also used to compute impingement limits and local collection efficiencies, which are used in subsequent ice-growth calculations and are also useful for designing systems to protect against icing. Next, the software predicts the shape of accumulating ice by modeling transfers of mass and energy in small control volumes. The foregoing computations are repeated over several computational time steps until the total icing exposure time is reached. Results of computations by LEWICE have been compared with an extensive database of measured ice shapes obtained from experiments, and have been shown to closely approximate those shapes under most conditions of interest to the aviation community.

Bidwell, Colin↗
NASA NTRS1993

Roughness of Weddell Sea ice and estimates of the air-ice drag coefficient

Snow surface, ice surface, and ice underside roughness computed from 47 surface elevation profiles collected during a transect of the Weddell Sea are reported. Roughness spectra for all three surfaces and for both deformed and undeformed ice roll off roughly as k exp -1 when the wavenumber k is between 0.1 and 3 rad/m. The snow surface and underside spectra roll off somewhat faster than k exp -1, and the ice surface spectra roll off somewhat slower than k exp -1. It is argued that the excess spectral intensity at high wavenumbers in the Antarctic ice surface spectra results from the small-scale roughness that the ice sheet had on consolidation. The snow surface spectra are used to compute a roughness scale that parameterizes the air-ice momentum coupling and makes it possible to estimate the neutral stability drag coefficient referenced to a height of 10 m, C sub DN10. Typical C sub DN10 values are 1.1-1.4 x 10 exp -3 over undeformed ice and 1.3-1.8 x 10 exp -3 over deformed ice.

Andreas, Edgar L.↗
DOE OSTI2025

A Machine Learning Framework for Predicting Microphysical Properties of Ice Crystals From Cloud Particle Imagery

The microphysical properties of ice crystals are important because they significantly alter the radiative properties and spatiotemporal distributions of clouds, which in turn strongly affect Earth's climate. However, it is challenging to measure key properties of ice crystals, such as mass or morphological features. Here, we present a proof-of-concept framework for predicting three-dimensional (3D) microphysical properties of ice crystals from in situ two-dimensional (2D) imagery. First, we computationally generated synthetic ice crystals using 3D modeling software along with geometric parameters estimated from the 2021 Ice Cryo-Encapsulation Balloon (ICEBall) field campaign. Then, we used synthetic crystals to train machine learning (ML) models to predict effective density ($ρ_e$), effective surface area ($A_e$), and number of bullets ($N_b$) from synthetic rosette imagery. On unseen synthetic images, our ML models accurately predicted ice crystal properties. ResNet-18 performed best, achieving $R^2$ values of 0.99 and 0.98 for $ρ_e$ and $A_e$, respectively, and MAE of 0.10 for mathematical equation in single view tasks. Stereo view ResNet-18 further reduced RMSE by 40% for $ρ_e$ and $A_e$ and reduced MAE by 0.08 for $N_b$. This work provides a novel ML-driven framework for estimating ice microphysical properties from in situ imagery, which will allow for downstream constraints on microphysical parameterizations, such as the mass-size relationship.

Ko, J. [Columbia Univ., New York, NY (United State↗
NASA NTRS1994

An Experimental Study of the Flowfield on a Semispan Rectangular Wing with a Simulated Glaze Ice Accretion

Wind tunnel experiments were conducted in order to study the effect of a simulated glaze ice accretion on the flowfield of a semispan, reflection-plane, rectangular wing at Re = 1.5 million and M = 0.12. A laser Doppler velocimeter was used to map the flowfield on the upper surface of the model in both the clean and iced configurations at alpha = 0, 4, and 8 degrees angle of attack. At low angles of attack, the massive separation bubble aft of the leading edge ice horn was found to behave in a manner similar to laminar separation bubbles. At alpha = 0 and 4 degrees, the locations of transition and reattachment, as deduced from momentum thickness distributions, were found to be in good agreement with transition and reattachment locations in laminar separation bubbles. These values at y/b = 0.470, the centerline measurement location, matched well with data obtained on a similar but two dimensional model. The measured velocity profiles on the iced wing compared reasonably with the predicted profiles from Navier-Stokes computations. The iced-induced separation bubble was also found to have features similar to the recirculating region aft of rearward-facing steps. At alpha = 0 degrees and 4 degrees, reverse flow magnitudes and turbulence intensity levels were typical of those found in the recirculating region aft of rearward-facing steps. The calculated separation streamline aft of the ice horn at alpha = 4 degrees, y/b = 0.470 coincided with the locus of the maximum Reynolds normal stress. The maximum Reynolds normal stress peaked at two locations along the separation streamline. The location of the first peak-value coincided with the transition location, as deduced from the momentum thickness distributions. The location of the second peak was just upstream of reattachment, in good agreement with measurements of flows over similar obstacles. The intermittency factor in the vicinity of reattachment at alpha = 4 degrees, y/b = 0.470, revealed the time-dependent nature of the reattachment process. The size and extent of the separation bubble were found to be a function of angle of attack and the spanwise location. Three dimensional effects were found to be strongest at alpha = 8 degrees. The calculated separation and stagnation streamlines were found to vary little with spanwise location at alpha = 0 degrees. The calculated separation streamlines at alpha = 4 degrees revealed that the bubble was largest near the centerline measurement plane, whereas the tip-induced vortex flow and the model root-tunnel wall boundary-layer interaction reduced the size of the bubble. These effects were found to be most dramatic at alpha = 8 degrees.

Khodadoust, Abdollah↗
NASA NTRS2009

Airframe Icing Research Gaps: NASA Perspective

qCurrent Airframe Icing Technology Gaps: Development of a full 3D ice accretion simulation model. Development of an improved simulation model for SLD conditions. CFD modeling of stall behavior for ice-contaminated wings/tails. Computational methods for simulation of stability and control parameters. Analysis of thermal ice protection system performance. Quantification of 3D ice shape geometric characteristics Development of accurate ground-based simulation of SLD conditions. Development of scaling methods for SLD conditions. Development of advanced diagnostic techniques for assessment of tunnel cloud conditions. Identification of critical ice shapes for aerodynamic performance degradation. Aerodynamic scaling issues associated with testing scale model ice shape geometries. Development of altitude scaling methods for thermal ice protections systems. Development of accurate parameter identification methods. Measurement of stability and control parameters for an ice-contaminated swept wing aircraft. Creation of control law modifications to prevent loss of control during icing encounters. 3D ice shape geometries. Collection efficiency data for ice shape geometries. SLD ice shape data, in-flight and ground-based, for simulation verification. Aerodynamic performance data for 3D geometries and various icing conditions. Stability and control parameter data for iced aircraft configurations. Thermal ice protection system data for simulation validation.

Potapczuk, Mark↗
NASA NTRS2015

A Fast Visible-Infrared Imaging Radiometer Suite Simulator for Cloudy Atmopheres

A fast instrument simulator is developed to simulate the observations made in cloudy atmospheres by the Visible Infrared Imaging Radiometer Suite (VIIRS). The correlated k-distribution (CKD) technique is used to compute the transmissivity of absorbing atmospheric gases. The bulk scattering properties of ice clouds used in this study are based on the ice model used for the MODIS Collection 6 ice cloud products. Two fast radiative transfer models based on pre-computed ice cloud look-up-tables are used for the VIIRS solar and infrared channels. The accuracy and efficiency of the fast simulator are quantify in comparison with a combination of the rigorous line-by-line (LBLRTM) and discrete ordinate radiative transfer (DISORT) models. Relative errors are less than 2 for simulated TOA reflectances for the solar channels and the brightness temperature differences for the infrared channels are less than 0.2 K. The simulator is over three orders of magnitude faster than the benchmark LBLRTM+DISORT model. Furthermore, the cloudy atmosphere reflectances and brightness temperatures from the fast VIIRS simulator compare favorably with those from VIIRS observations.

VIIRS Instrument Simulator↗
NASA NTRSDate not supplied

Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Model

In-flight icing is an important consideration that affects aircraft design, performance and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of ice accretion testing conducted in the NASA Icing Research Tunnel on a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small MVD = 25 μm and a large MVD = 110 μm. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. The results for the small MVD cloud conditions were consistent with previous work and showed the changing ice morphology from rime ice at colder temperatures to highly 3D scallop ice in the range of -8.7 ˚C to -3.8 ˚C. The results for the large MVD cloud conditions exhibited some differences in the ice accretion morphology from the small MVD conditions. Scallop-like ice features were observed at total temperature of -3.8˚C. The large MVD ice accretion at lower temperatures tended to be smoother than the corresponding ice accretion for small MVD. The thickness of the main ice shape tended to be larger for the small MVD conditions compared to the corresponding large MVD ice shape while the latter had more ice farther downstream. The measured ice mass was approximately equal between the corresponding small and large MVD ice shapes up to a freezing fraction of 0.6. For higher freezing fraction, the large MVD ice shapes weighed more. Cloud MVD variation from 50 μm to 230 μm while holding the scaling parameters constant resulted in very minor differences in the ice shapes, mostly in the upper and lower surface chordwise extents. Ice volume was computed from 3D scan data and used to calculate a ratio of ice mass to volume. The resulting values were in the range of 240 to 455 Kg/m3. While this is consistent with analogous values previously reported in the literature, more data are needed to determine a specific trend in this ratio as a function of freezing fraction and the corresponding governing physical phenomena. This work has resulted in a significant experimental database of ice accretion for small and large MVD conditions applicable to large-scale swept wings.

Icing↗
NASA NTRSDate not supplied

Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section

In-flight icing is an important consideration that affects aircraft design, performance and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of ice accretion testing conducted in the NASA Icing Research Tunnel on a large swept wing typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small MVD = 25 μm and a large MVD = 110 μm. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. The results for the small MVD cloud conditions were consistent with previous work and showed the changing ice morphology from rime ice at colder temperatures to highly 3D scallop ice in the range of -8.7 ˚C to -3.8 ˚C. The results for the large MVD cloud conditions exhibited some differences in the ice accretion morphology. Scallop-like ice features were observed at total temperature of -3.8˚C. The large MVD ice accretion at lower temperatures tended to be smoother than the corresponding ice accretion for small MVD. The thickness of the main ice shape tended to be larger for the small MVD conditions compared to the corresponding large MVD ice shape while the latter had more ice farther downstream. The measured ice mass was approximately equal between the corresponding small and large MVD ice shapes up to a freezing fraction of 0.6. For higher freezing fraction, the large MVD ice shapes weighed more. Cloud MVD variation from 50 μm to 230 μm while holding the scaling parameters constant resulted in very minor differences in the ice shapes, mostly in the upper and lower surface chordwise extents. Ice volume was computed from 3D scan data and used to calculate a ratio of ice mass to volume. The resulting values were in the range of 240 to 455 Kg/m3. While this is consistent with analogous values previously reported in the literature, more data are needed to determine a specific trend in this ratio as a function of freezing fraction and the corresponding governing physical phenomena. This work has resulted in a significant experimental database of ice accretion for small and large MVD conditions applicable to large-scale swept wings.

Icing↗
NASA NTRS1995

Predictions Of Accretion Of Ice On A Commercial Airplane

Report presents computational study of accretion of ice on wings and on horizontal tail surfaces of simplified version of Boeing 737 airplane. Calculations performed by use of LEWICE3D computer code, which incorporates panel-method mathematical model of three-dimensional potential flow and mathematical models of trajectories of drops of water entrained in flow, impingement of drops on surfaces, and transfer of heat.

Bidwell, C. S.↗
NASA NTRS2004

Preparing and Analyzing Iced Airfoils

SmaggIce version 1.2 is a computer program for preparing and analyzing iced airfoils. It includes interactive tools for (1) measuring ice-shape characteristics, (2) controlled smoothing of ice shapes, (3) curve discretization, (4) generation of artificial ice shapes, and (5) detection and correction of input errors. Measurements of ice shapes are essential for establishing relationships between characteristics of ice and effects of ice on airfoil performance. The shape-smoothing tool helps prepare ice shapes for use with already available grid-generation and computational-fluid-dynamics software for studying the aerodynamic effects of smoothed ice on airfoils. The artificial ice-shape generation tool supports parametric studies since ice-shape parameters can easily be controlled with the artificial ice. In such studies, artificial shapes generated by this program can supplement simulated ice obtained from icing research tunnels and real ice obtained from flight test under icing weather condition. SmaggIce also automatically detects geometry errors such as tangles or duplicate points in the boundary which may be introduced by digitization and provides tools to correct these. By use of interactive tools included in SmaggIce version 1.2, one can easily characterize ice shapes and prepare iced airfoils for grid generation and flow simulations.

Vickerman, Mary B.↗
NASA NTRS1972

On radiational cooling computations in clouds.

Expressions are derived for the cooling or heating rates for clouds in which the condensed phase is either ice or water. Values are computed for ice and water clouds over a reasonable temperature range for pressures of 1000, 500, and 200 mb. The importance is shown of adequately allowing for the latent load in computations of radiative cooling or heating rates based on measurements of radiative divergence in clouds. It is shown that, other things being equal, the effect of the latent load is always greatest at low levels and higher temperatures.

Knollenberg, R. G.↗
NASA NTRS2009

Interactive, Automated Management of Icing Data

IceVal DatAssistant is software (see figure) that provides an automated, interactive solution for the management of data from research on aircraft icing. This software consists primarily of (1) a relational database component used to store ice shape and airfoil coordinates and associated data on operational and environmental test conditions and (2) a graphically oriented database access utility, used to upload, download, process, and/or display data selected by the user. The relational database component consists of a Microsoft Access 2003 database file with nine tables containing data of different types. Included in the database are the data for all publicly releasable ice tracings with complete and verifiable test conditions from experiments conducted to date in the Glenn Research Center Icing Research Tunnel. Ice shapes from computational simulations with the correspond ing conditions performed utilizing the latest version of the LEWICE ice shape prediction code are likewise included, and are linked to the equivalent experimental runs. The database access component includes ten Microsoft Visual Basic 6.0 (VB) form modules and three VB support modules. Together, these modules enable uploading, downloading, processing, and display of all data contained in the database. This component also affords the capability to perform various database maintenance functions for example, compacting the database or creating a new, fully initialized but empty database file.

Levinson, Laurie H.↗
NASA NTRS2004

Evaluation of Droplet Splashing Algorithm in LEWICE 3.0

The Icing Branch at NASA Glenn Research has developed a computer program to simulate ice formation on the leading edge of an aircraft wing during flight through cold, moist air. As part of the branch's current research, members have developed software known as LEWICE. This program is capable of predicting the formation of ice under designated weather conditions. The success of LEWICE is an asset to airplane manufacturers, ice protection system manufacturers, and the airline industry. Simulations of ice formation conducted in the tunnel and in flight is costly and time consuming. However, the danger of in-flight icing continues to be a concern for both commercial and military pilots. The LEWICE software is a step towards inexpensive and time efficient prediction of ice collection. In the most recent version of the program, LEWICE contains an algorithm for droplet splashing. Droplet splashing is a natural occurrence that causes accumulation of ice on aircraft surfaces. At impingement water droplets lose a portion of their mass to splashing. With part of each droplet joining the airflow and failing to freeze, early versions of LEWICE without the splashing algorithm over-predicted the collection of ice on the leading edge. The objective of my project was to determine whether the revised version of LEWICE accurately reflected the ice collection data obtained from the Icing Research Tunnel (IRT). The experimental data from the IRT was collected by Mark Potapczuk in January, March and July of 2001 and April and December of 2002. Experimental data points were the result of ice tracings conducted shortly after testing in the tunnel. Run sheets, which included a record of velocity, temperature, liquid water content and droplet diameter, served as the input of the LEWICE computer program. Parameters identical to the tunnel conditions were used to run LEWICE 2.0 and LEWICE 3.0. The results from IRT and versions of LEWICE were compared graphically. After entering the raw experimental data and computer output into a spread sheet, I mapped each ice formation onto a clean airfoil. The LEWICE output provided the data points to graphically depict ice formations developed by the program. weather conditions of runs conducted in January 2001, it was evident that the splashing algorithm of LEWICE 3.0 predicts ice formations more accurately than LEWICE 2.0. Especially at conditions with droplet size between 80 and 160 microns, the splashing algorithm of the new LEWICE version compensated for the loss of droplet mass as a result of splashing. In contrast, LEWICE 2.0 consistently over-predicted the mass of the ice in conditions with droplet size exceeding 80 microns. This evidence confirms that changes made to algorithms of LEWICE 3.0 have increased the accuracy of predicting ice collection.

Homenko, Hilary N.↗
NASA NTRS2004

Detached-Eddy Simulations of Separated Flow Around Wings With Ice Accretions: Year One Report

A computational investigation was performed to assess the effectiveness of Detached-Eddy Simulation (DES) as a tool for predicting icing effects. The AVUS code, a public domain flow solver, was employed to compute solutions for an iced wing configuration using DES and steady Reynolds Averaged Navier-Stokes (RANS) equation methodologies. The configuration was an extruded GLC305/944-ice shape section with a rectangular planform. The model was mounted between two walls so no tip effects were considered. The numerical results were validated by comparison with experimental data for the same configuration. The time-averaged DES computations showed some improvement in lift and drag results near stall when compared to steady RANS results. However, comparisons of the flow field details did not show the level of agreement suggested by the integrated quantities. Based on our results, we believe that DES may prove useful in a limited sense to provide analysis of iced wing configurations when there is significant flow separation, e.g., near stall, where steady RANS computations are demonstrably ineffective. However, more validation is needed to determine what role DES can play as part of an overall icing effects prediction strategy. We conclude the report with an assessment of existing computational tools for application to the iced wing problem and a discussion of issues that merit further study.

Choo, Yung K.↗
NASA NTRS2015

The Influence of Sea Ice on Arctic Low Cloud Properties and Radiative Effects

The Arctic is one of the most climatically sensitive regions of the Earth. Climate models robustly project the Arctic to warm 2-3 times faster than the global mean surface temperature, termed polar warming amplification (PWA), but also display the widest range of surface temperature projections in this region. The response of the Arctic to increased CO2 modulates the response in tropical and extra-tropical regions through teleconnections in the atmospheric circulation. An increased frequency of extreme precipitation events in the northern mid-latitudes, for example, has been linked to the change in the background equator-to-pole temperature gradient implied by PWA. Understanding the Arctic climate system is therefore important for predicting global climate change. The ice albedo feedback is the primary mechanism driving PWA, however cloud and dynamical feedbacks significantly contribute. These feedback mechanisms, however, do not operate independently. How do clouds respond to variations in sea ice? This critical question is addressed by combining sea ice, cloud, and radiation observations from satellites, including CERES, CloudSAT, CALIPSO, MODIS, and microwave radiometers, to investigate sea ice-cloud interactions at the interannual timescale in the Arctic. Cloud characteristics are strongly tied to the atmospheric dynamic and thermodynamic state. Therefore, the sensitivity of Arctic cloud characteristics, vertical distribution and optical properties, to sea ice anomalies is computed within atmospheric dynamic and thermodynamic regimes. Results indicate that the cloud response to changes in sea ice concentration differs significantly between atmospheric state regimes. This suggests that (1) the atmospheric dynamic and thermodynamic characteristics and (2) the characteristics of the marginal ice zone are important for determining the seasonal forcing by cloud on sea ice variability.

Taylor, Patrick C.↗
NASA NTRS1954

Experimental Determination of Thermal Conductivity of Low-Density Ice

The thermal conductivity of low-density ice has been computed from data obtained in an experimental investigation of the heat transfer and mass transfer by sublimation for an iced surface on a flat plate in a high-velocity tangential air stream. The results are compared with data from several sources on the thermal conductivity of packed snow and solid glaze ice. The results show good agreement with the equations for the thermal conductivity of packed snow as a function of snow density. The agreement of the curves for packed snow near the solid ice regime with the values of thermal conductivity, of ice indicates that the curves are applicable over the entire-ice-density range.

Coles, Willard D.↗