SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Preparation, Identification, and Low-Temperature Infrared Spectra of Two Elusive Crystalline Nitrile Ices

Infrared (IR) spectra of the alkyl nitrile ices CH3CN and CH3CH2CN are relevant to the study of the low-temperature chemistry of objects within and beyond the Solar System, but the thermodynamically favored low-temperature crystalline phases of these compounds have not been presented and identified in the planetary-sciences literature. Moreover, there is a large variation in the published IR spectra of these two icy solids, including spectra that are used for the analyses of spacecraft data from the Voyager and Cassini missions. Here the IR spectra of the low-temperature crystalline forms of CH3CN and CH3CH2CN, which are the thermodynamically stable phases at Titan temperatures, are presented for the first time with all samples being made by vapor-phase deposition. Conditions are described for producing these ice phases for both compounds, and new measurements are reported of ice density and refractive index, quantities needed to compute IR absorption coefficients, band strengths, optical constants, and, ultimately, nitrile abundances in extraterrestrial environments. Comparisons are drawn between CH3CN and CH3CH2CN, an earlier prediction on their similar crystallization behaviors is verified, and previous work is examined in light of these new results, including a counter-intuitive observation in which the low-temperature phase of an ice is made by heating a high-temperature phase. Applications and extensions are described.

Reggie L Hudson↗

NASA NTRS2019

Preparation, Identification, and Low-Temperature Infrared Spectra of Two Elusive Crystalline Nitrile Ices

Ices, IR spectroscopy, Titan, organic chemistry, i↗

NASA NTRS2017

Modeling of a Turbofan Engine with Ice Crystal Ingestion in the NASA Propulsion System Laboratory

The main focus of this study is to apply a computational tool for the flow analysis of the turbine engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) at NASA Glenn Research Center. The PSL has been used to test a highly instrumented Honeywell ALF502R-5A (LF11) turbofan engine at simulated altitude operating conditions. Test data analysis with an engine cycle code and a compressor flow code was conducted to determine the values of key icing parameters, that can indicate the risk of ice accretion, which can lead to engine rollback (un-commanded loss of engine thrust). The full engine aerothermodynamic performance was modeled with the Honeywell Customer Deck specifically created for the ALF502R-5A engine. The mean-line compressor flow analysis code, which includes a code that models the state of the ice crystal, was used to model the air flow through the fan-core and low pressure compressor. The results of the compressor flow analyses included calculations of the ice-water flow rate to air flow rate ratio (IWAR), the local static wet bulb temperature, and the particle melt ratio throughout the flow field. It was found that the assumed particle size had a large effect on the particle melt ratio, and on the local wet bulb temperature. In this study the particle size was varied parametrically to produce a non-zero calculated melt ratio in the exit guide vane (EGV) region of the low pressure compressor (LPC) for the data points that experienced a growth of blockage there, and a subsequent engine called rollback (CRB). At data points where the engine experienced a CRB having the lowest wet bulb temperature of 492 degrees Rankine at the EGV trailing edge, the smallest particle size that produced a non-zero melt ratio (between 3 percent - 4 percent) was on the order of 1 micron. This value of melt ratio was utilized as the target for all other subsequent data points analyzed, while the particle size was varied from 1 micron - 9.5 microns to achieve the target melt ratio. For data points that did not experience a CRB which had static wet bulb temperatures in the EGV region below 492 degrees Rankine, a non-zero melt ratio could not be achieved even with a 1 micron ice particle size. The highest value of static wet bulb temperature for data points that experienced engine CRB was 498 degrees Rankine with a particle size of 9.5 microns. Based on this study of the LF11 engine test data, the range of static wet bulb temperature at the EGV exit for engine CRB was in the narrow range of 492 degrees Rankine - 498 degrees Rankine , while the minimum value of IWAR was 0.002. The rate of blockage growth due to ice accretion and boundary layer growth was estimated by scaling from a known blockage growth rate that was determined in a previous study. These results obtained from the LF11 engine analysis formed the basis of a unique “icing wedge.”

Turbo engine↗

NASA NTRS1990

Titan and other icy satellites - Dielectric properties of constituent materials and implications for radar sounding

Radar sounding observations of Titan and the other outer solar system icy satellites have been conducted in order to ascertain the materials present, their dielectric properties, and the degree of radar signal penetration. Deep penetration is possible at long wavelengths, while strong reflections may occur due to pockets or layers of H2O-NH3 solution. Attention is given to the case of Titan, using models to calculate radar reflection S/N ratios; the models involve H2O ice covered by a 100-m thick layer of organic sediments that are overlain by an ocean of liquid CH4-C2H6-N2. At 13.6-cm wavelength, a greater-than-10 dB S/N is computed from the ocean-sediment and sediment-ice interfaces, suggesting that the geology of Titan may be investigated by such a sounding mode despite deep ocean coverage of the satellite.

Thompson, W. Reid↗

DOE OSTI2024

Quantifying the compound hazard of freezing rain and wind gusts across CONUS

Abstract The co-occurrence of freezing rain, ice accumulation and wind gusts (FZG) poses a significant hazard to infrastructure and transportation. However, quantification of the frequency and intensity of FZG is challenged by the lack of direct icing measurements. In this work, we evaluate and then apply an energy balance model to high-frequency data collected during 2005–2022 to derive hourly ice accumulation at 883 stations across the contiguous USA. These estimates are combined with wind gust observations to compute time series of hourly FZG hazard magnitude using the Sperry–Piltz Ice Accumulation (SPIA) index. Results are evaluated using Storm Reports of damage and economic disruption. The hourly SPIA estimates are also used to (i) derive a geospatial atlas of the hazard including the 50 yr return period event intensities for each US state derived using superstations, and (ii) describe storylines of significant events in terms of meteorological drivers and socioeconomic impacts. The highest values of SPIA during the 18 yr study period occur in a region extending from the Southern Great Plains over the Midwest into the densely populated Northeast. States in these regions also have high 50 yr return period maximum radial ice accumulation of 3–5 cm and co-occurring wind gusts >30 ms −1 . These values are comparable to past estimates for the 500 yr event which may imply this hazard has been previously underestimated. This atlas can be used to inform optimal FZG hazard mitigation strategies for each state/region.

58 GEOSCIENCES↗

DOE OSTI2025

Strong Scalability Analysis of the Albany Land Ice code on HPC Architectures

Scalability is a critical factor in High-Performance Computing (HPC), where optimizing resource usage has a direct impact on cost-effectiveness and time-efficiency. This report presents a strong scaling performance study of the Albany Land Ice (ALI) code across different HPC architectures, towards determining the best configuration to use when running large-scale simulation ensembles.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗

NASA NTRS1993

Arctic sea ice signatures for passive microwave algorithms

Two new sets of global pure type signatures are defined, and the precision of the concentration estimates made using the Team algorithm and the new signatures is estimated. One set varies temporally but not spatially and is defined by the modes of major peaks in histograms of brightness temperature. The other set is defined as the overall means of the temporally varying signatures; these are called 'new constant signatures'. Both sets of new signatures have essentially the same precision: 0.05 (fractional area) for winter ice concentration, 0.10 for summer ice concentration, and 0.15 to 0.3 for winter multilayer ice concentrations. Both sets result in about 0.6 to 0.7 winter multiyear ice for the Arctic Ocean, substantially more than the 0.3 to 0.4 computed using the original team algorithm signatures, and in the same range as several independent estimates of the multiyear ice, 0.54 to 0.75.

Thomas, D. R.↗

DOE OSTI2024

Disk Failure Dataset from the Campaign Storage System

This dataset consists of 1,389 disk (HDD) failure events collected from the Campaign storage system at LANL. The Campaign system supported various compute platforms throughout its lifespan, including Cielo, Fire, Ice, and notably, the Trinity supercomputer. Each recorded event includes its detection timestamp (in ISO 8601 format) and details such as its location within the storage systemâ€”rack, enclosure, and drive slot number. The data, spanning from May 4, 2021, to July 25, 2023 (2 years, 2 months, and 22 days), represents failure events from the terminal years of Campaign's operational period, accounting for 26% of its total operational time.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1996

User Interface Developed for Controls/CFD Interdisciplinary Research

The NASA Lewis Research Center, in conjunction with the University of Akron, is developing analytical methods and software tools to create a cross-discipline "bridge" between controls and computational fluid dynamics (CFD) technologies. Traditionally, the controls analyst has used simulations based on large lumping techniques to generate low-order linear models convenient for designing propulsion system controls. For complex, high-speed vehicles such as the High Speed Civil Transport (HSCT), simulations based on CFD methods are required to capture the relevant flow physics. The use of CFD should also help reduce the development time and costs associated with experimentally tuning the control system. The initial application for this research is the High Speed Civil Transport inlet control problem. A major aspect of this research is the development of a controls/CFD interface for non-CFD experts, to facilitate the interactive operation of CFD simulations and the extraction of reduced-order, time-accurate models from CFD results. A distributed computing approach for implementing the interface is being explored. Software being developed as part of the Integrated CFD and Experiments (ICE) project provides the basis for the operating environment, including run-time displays and information (data base) management. Message-passing software is used to communicate between the ICE system and the CFD simulation, which can reside on distributed, parallel computing systems. Initially, the one-dimensional Large-Perturbation Inlet (LAPIN) code is being used to simulate a High Speed Civil Transport type inlet. LAPIN can model real supersonic inlet features, including bleeds, bypasses, and variable geometry, such as translating or variable-ramp-angle centerbodies. Work is in progress to use parallel versions of the multidimensional NPARC code.

Source record↗

DOE OSTI2020

Arctic Tipping Points Triggering Global Change (LDRD Final Report)

The Arctic is warming and feedbacks in the coupled Earth system may be driving the Arctic to tipping events that could have critical downstream impacts for the rest of the globe. In this project we have focused on analyzing sea ice variability and loss in the coupled Earth system Summer sea ice loss is happening rapidly and although the loss may be smooth and reversible, it has significant consequences for other Arctic systems as well as geopolitical and economic implications. Accurate seasonal predictions of sea ice minimum extent and long-term estimates of timing for a seasonally ice-free Arctic depend on a better understanding of the factors influencing sea ice dynamics and variation in this strongly coupled system. Under this project we have investigated the most influential factors in accurate predictions of September Arctic sea ice extent using machine learning models trained separately on observational data and on simulation data from five E3SM historical ensembles. Monthly averaged data from June, July, and August for a selection of ice, ocean, and atmosphere variables were used to train a random forest regression model. Gini importance measures were computed for each input feature with the testing data. We found that sea ice volume is most important earlier in the season (June) and sea ice extent became a more important predictor closer to September. Results from this study provide insight into how feature importance changes with forecast length and illustrates differences between observational data and simulated Earth system data. We have additionally performed a global sensitivity analysis (GSA) using a fully coupled ultra- low resolution configuration E3SM. To our knowledge, this is the first global sensitivity analysis involving the fully-coupled E3SM Earth system model. We have found that parameter variations show significant impact on the Arctic climate state and atmospheric parameters related to cloud parameterizations are the most significant. We also find significant interactions between parameters from different components of E3SM. The results of this study provide invaluable insight into the relative importance of various parameters from the sea ice, atmosphere and ocean components of the E3SM (including cross-component parameter interactions) on various Arctic-focused quantities of interest (QOIs).

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1972

ERTS computer compatible tape data processing and analysis. Appendix 1: The utility of imaging radars for the study of lake ice

There are no author-identified significant results in this report. Remotely sensed multispectral scanner and return beam vidicon imagery from ERTS-1 is being used for: (1) water depth measurements in the Virgin Islands and Upper Lake Michigan areas; (2) mapping of the Yellowstone National Park; (3) assessment of atmospheric effects in Colorado; (4) lake ice surveillance in Canada and Great Lakes areas; (5) recreational land use in Southeast Michigan; (6) International Field Year on the Great Lakes investigations of Lake Ontario; (7) image enhancement of multispectral scanner data using existing techniques; (8) water quality monitoring of the New York Bight, Tampa Bay, Lake Michigan, Santa Barbara Channel, and Lake Erie; (9) oil pollution detection in the Chesapeake Bay, Gulf of Mexico southwest of New Orleans, and Santa Barbara Channel; and (10) mapping iron compounds in the Wind River Mountains.

Polcyn, F. C.↗

NASA NTRS1983

An analytical evaluation of the icing properties of several low and medium speed airfoils

A droplet trajectory computer code is used to analyze the water droplet impingement characteristics of several low and medium speed airfoils. The results are reported in terms of the airfoil maximum impingement efficiency, and limits of impingement. The airfoils are analyzed for angles of attack of -4 to 12 degrees, and modified inertia parameters of 0.01 to 1.0. These data are then used to evaluate the effect of airfoil geometry: leading edge radius, maximum thickness, maximum camber, and angle of attack on airfoil icing characteristics. Airfoils from the NACA four digit series to the new NASA NLF airfoils are considered. The strongest correlations were seen between airfoil leading edge radius and maximum local collection efficiency, thickness and total collection efficiency, and camber and the limits of impingement.

Bragg, M. B.↗

NASA NTRS2008

Wide-Band Radar for Measuring Thickness of Sea Ice

A wide-band penetrating radar system for measuring the thickness of sea ice is under development. The need for this or a similar system arises as follows: Spatial and temporal variations in the thickness of sea ice are important indicators of heat fluxes between the ocean and atmosphere and, hence, are important indicators of climate change in polar regions. A remote-sensing system that could directly measure the thickness of sea ice over a wide thickness range from aboard an aircraft or satellite would be of great scientific value. Obtaining thickness measurements over a wide region at weekly or monthly time intervals would contribute significantly to understanding of changes in the spatial distribution and of the mass balance of sea ice. A prototype of the system was designed on the basis of computational simulations directed toward understanding what signal frequencies are needed to satisfy partly competing requirements to detect both bottom and top ice surfaces, obtain adequate penetration despite high attenuation in the lossy sea-ice medium, and obtain adequate resolution, all over a wide thickness range. The prototype of the system is of the frequency-modulation, continuous-wave (FM-CW) type. At a given time, the prototype functions in either of two frequency-band/operational-mode combinations that correspond to two thickness ranges: a lower-frequency (50 to 250 MHz) mode for measuring thickness greater than about 1 m, and a higher frequency (300 to 1,300 MHz) mode for measuring thickness less than about 1 m. The bandwidth in the higher-frequency (lesser-thickness) mode is adequate for a thickness resolution of 15 cm; the bandwidth in the lower-frequency (greater-thickness) mode is adequate for a thickness resolution of 75 cm. Although a thickness resolution of no more than 25 cm is desired for scientific purposes, the 75-cm resolution was deemed acceptable for the purpose of demonstrating feasibility. The prototype was constructed as a modified version of a 500-to-2,000-MHz FM-CW radar system developed previously for mapping near-surface internal layers of the Greenland ice sheet. The prototype included two sets of antennas: one for each frequency-band/mode. For Arctic and Antarctic field tests, the prototype was mounted on a sled that was towed across the ice. The Arctic field test was performed in the lower-frequency mode on ice ranging in thickness from 1 to 4 m. In the analysis of the results of the Arctic field test, a comparison of the radar-determined ice thicknesses with actual ice thicknesses yielded an overall mean difference of 14 cm and standard deviation of 30 cm. The Antarctic field test was performed in the higher-frequency mode; analysis of the results led to the conclusion that this mode is useful for measuring thicknesses between 0.5 and 1 m. Several modifications have been conceived for implementation in further development toward an improved practical system: The system would function in a single frequency-band/mode (100 to 1,200 MHz) that would afford a resolution of about 15 cm. There would be a single antenna system that would be optimized for the entire 100-to-1,000-MHz frequency band. To enable ice-thickness surveys over larger areas, the system would be made capable of operating aboard a low-flying aircraft that could be either piloted or robotic. Data-processing techniques to deconvolve the system response have been developed on the basis of impulse-response measurements over a calm ocean. Implementation of these techniques in the system would enable correction for imperfections of the system and would thereby increase the effective sensitivity of the system.

Gogineni, Prasad↗

NASA NTRS2000

Scattering Properties of Needle-Like and plate-like Ice Spheroids with Moderate Size Parameters

We use the current advanced version of the T-matrix method to compute the optical cross sections, the asymmetry parameter of the phase function, and the scattering matrix elements of ice spheroids with aspect ratios up to 20 and surface-equivalent-sphere size parameters up to 12. We demonstrate that plate-like and needle-like particles with moderate size parameters possess unique scattering properties: their asymmetry parameters and phase functions are similar to those of surface-equivalent spheres, whereas all other elements of the scattering matrix are typical of particles much smaller than the wavelength (Rayleigh scatterers). This result may have important implications for optical particle sizing and remote sensing of the terrestrial and planetary atmospheres.

Zakharova, Nadia T.↗

NASA NTRS1990

Microwave radiative transfer studies using combined multiparameter radar and radiometer measurements during Cohmex

Theoretical calculations of the upwelling microwave radiances from clouds containing layers of rain, ice, and a melting region were performed at frequencies of 18, 37, and 92 GHz. These frequencies coincide with high-resolution microwave radiometer measurements taken aboard the NASA ER-2 high-altitude aircraft during the summer 1986 Cohmex (Cooperative Huntsville Meteorological Experiment) in Alabama. For purposes of brightness temperature computations, the storms were modeled with rain, melting phase, and ice layers. The melting phase region was composed of water-coated ice spheres defined by a 'melt index' in terms of the volume fraction of water. Single scatter albedo, scattering, and extinction coefficients were computed at the above frequencies as a function of the rain rate and melt index. In addition, multiparameter radar observations of the storm were mapped into a Cartesian space and averaged over regions comparable to the radiometer footprint.

Vivekanandan, J.↗

NASA NTRS1989

A laboratory investigation into microwave backscattering from sea ice

The sources of scattering of artificial sea ice were determined, backscatter measurements semi-quantitatively were compared with theoretical predictions, and inexpensive polarimetric radars were developed for sea ice backscatter studies. A brief review of the dielectric properties of sea ice and of commonly used surface and volume scattering theories is presented. A description is provided of the backscatter measurements performed and experimental techniques used. The development of inexpensive short-range polarimetric radars is discussed. The steps taken to add polarimetric capability to a simple FM-W radar are considered as are sample polarimetric phase measurements of the radar. Ice surface characterization data and techniques are discussed, including computation of surface rms height and correlation length and air bubble distribution statistics. A method is also presented of estimating the standard deviation of rms height and correlation length for cases of few data points. Comparisons were made of backscatter measurements and theory. It was determined that backscatter from an extremely smooth saline ice surface at C band cannot be attributed only to surface scatter. It was found that snow cover had a significant influence on backscatter from extremely smooth saline ice at C band.

Bredow, Jonathan W.↗

NASA NTRS1994

Numerical study of the effects of icing on viscous flow over wings

An improved hybrid method for computing unsteady compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved using a diagonal form of an alternating-direction implicit (ADI) approximate factorization procedure. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. Characteristic-based viscous/inviscid interface boundary conditions are employed to avoid spurious reflections at that interface. The resulting CPU times are less than 60 percent of that required for a full-blown Navier-Stokes analysis for steady flow applications and about 60 percent of the Navier-Stokes CPU times for unsteady flows in non-vector processing machines. Applications of the method are presented for a rectangular NACA 0012 wing in low subsonic steady flow at moderate and high angles of attack, and for an F-5 wing in steady and unsteady subsonic and transonic flows. Steady surface pressures are in very good agreement with experimental data and are essentially identical to Navier-Stokes predictions. Density contours show that shocks cross the viscous/inviscid interface smoothly, so that the accuracy of full Navier-Stokes equations can be retained with a significant savings in computational time.

Sankar, L. N.↗