Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Modeling of Commercial Turbofan Engine with Ice Crystal Ingestion; Follow-On

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that have been attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was degraded engine performance, and one or more of the following: loss of thrust control (roll back), compressor surge or stall, and flameout of the combustor. As ice crystals are ingested into the fan and low pressure compression system, the increase in air temperature causes a portion of the ice crystals to melt. It is hypothesized that this allows the ice-water mixture to cover the metal surfaces of the compressor stationary components which leads to ice accretion through evaporative cooling. Ice accretion causes a blockage which subsequently results in the deterioration in performance of the compressor and engine. The focus of this research is to apply an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion. The tool is comprised of an engine system thermodynamic cycle code, a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor flow path, without modeling the actual ice accretion. A commercial turbofan engine which has previously experienced icing events during operation in a high altitude ice crystal environment has been tested in the Propulsion Systems Laboratory (PSL) altitude test facility at NASA Glenn Research Center. The PSL has the capability to produce a continuous ice cloud which is ingested by the engine during operation over a range of altitude conditions. The PSL test results confirmed that there was ice accretion in the engine due to ice crystal ingestion, at the same simulated altitude operating conditions as experienced previously in flight. The computational tool was utilized to help guide a portion of the PSL testing, and was used to predict ice accretion could also occur at significantly lower altitudes. The predictions were qualitatively verified by subsequent testing of the engine in the PSL. In a previous study, analysis of select PSL test data points helped to calibrate the engine icing computational tool to assess the risk of ice accretion. This current study is a continuation of that data analysis effort. The study focused on tracking the variations in wet bulb temperature and ice particle melt ratio through the engine core flow path. The results from this study have identified trends, while also identifying gaps in understanding as to how the local wet bulb temperature and melt ratio affects the risk of ice accretion and subsequent engine behavior.

Jorgenson, Philip C. E.↗
NASA NTRSDate not supplied

Impact Ice Adhesion at NASA Glenn: Current Experimental Methods and Supporting Measurements

When examining the literature on the adhesion strength of impact ice, there have been a wide range of methodologies tried to measure the required stresses to induce interfacial delamination. Utilizing the Icing Research Tunnel at the NASA Glenn Research Center to generate the impact ice required for this work, several different mechanical tests have been and are being developed to determine the stresses along the interface between ice and coupon. This set of tests includes the technical mature modified lap joint test which has been used to conduct ice adhesion studies through a wide sweep of icing conditions. To conduct in situ ice adhesion measurements inside of the Icing Research Tunnel, several new experiments are currently being developed to make ice adhesion measurements during and immediately after ice accretion. In addition to these experimental methods, several supporting measurement techniques have been developed to allow for a better understanding on the influence of icing cloud conditions on the mechanical behavior of impact ice. Digital image correlation has been successfully implemented to augment the data generated by the modified lap joint test with full field surface displacement and strain measurements which allow for insight into the deformation processes present during a test. Both optical microscopy of impact ice samples along with ice replication techniques have been used to study the grain structure of the impact ice. This has led to a deeper understanding of the results from the modified lap joint method and how the structure of impact ice changes as it is accreted during an icing spray. The freezing process of impact ice generated by supercooled liquid water is not a volume conserving process, which leads to the presence of residual strains along the interface between ice and substrate. These strains have been observed using both a simplified flat geometry and a representative airfoil. The data gathered by these experimental adhesion methods and supporting measurements allows for a comprehensive understanding on the behavior of impact ice which will be critical to the development of future ice shedding models.

Ice Adhesion↗
NASA NTRSDate not supplied

Impact Ice Adhesion at NASA Glenn: Current Experimental Methods and Supporting Measurements

When examining the literature on the adhesion strength of impact ice, there have been a wide range of methodologies tried to measure the required stresses to induce interfacial delamination. Utilizing the Icing Research Tunnel at the NASA Glenn Research Center to generate the impact ice required for this work, several different mechanical tests have been and are being developed to determine the stresses along the interface between ice and coupon. This set of tests includes the technical mature modified lap joint test which has been used to conduct ice adhesion studies through a wide sweep of icing conditions. To conduct in situ ice adhesion measurements inside of the Icing Research Tunnel, several new experiments are currently being developed to make ice adhesion measurements during and immediately after ice accretion. In addition to these experimental methods, several supporting measurement techniques have been developed to allow for a better understanding on the influence of icing cloud conditions on the mechanical behavior of impact ice. Digital image correlation has been successfully implemented to augment the data generated by the modified lap joint test with full field surface displacement and strain measurements which allow for insight into the deformation processes present during a test. Both optical microscopy of impact ice samples along with ice replication techniques have been used to study the grain structure of the impact ice. This has led to a deeper understanding of the results from the modified lap join method and how the structure of impact ice changes as it is accreted during an icing spray. The freezing process of impact ice generated by supercooled liquid water is not a volume conserving process, which leads to the presence of residual strains along the interface between ice and substrate. These strains have been observed using both a simplified flat geometry and a representative airfoil. The data gathered by these experimental adhesion methods and supporting measurements allows for a compressive understanding on the behavior of impact ice which will be critical to the development of future ice shedding models.

Ice Adhesion↗
DOE OSTI2020

The Role of Ice Splintering on Microphysics of Deep Convective Clouds Forming Under Different Aerosol Conditions: Simulations Using the Model With Spectral Bin Microphysics

Observations during the Ice in Clouds Experiment-Tropical (ICE-T) field experiment show that the ice particles concentration in a developing deep convective clouds at the level of T = -15 °C reached about 500 L -1 , that is, many orders higher than that of ice-nucleating particle. To simulate microphysics of these clouds, the 2-D Hebrew University Cloud model (HUCM) with spectral bin microphysics was used in which two main types of ice multiplication mechanisms were included in addition to the Hallet-Mossop mechanism. In the first ice multiplication mechanism ice splinters form by drop freezing and drop-ice collisions. Ice multiplication of this type dominates during developing stage of cloud evolution, when liquid water content is significant. At later stage when clouds become nearly glaciated, ice crystals are produced largely by ice splintering during ice-ice collisions (the second ice multiplication mechanism). Simulations show that droplet size distributions, as well as size distributions of ice particles, agree well with the measurements during ICE-T. Simulations with different cloud condensation nuclei concentrations show the existence of the “optimum” cloud condensation nuclei concentration (or droplet concentration), at which concentration of ice splinters reaches maximum. In these simulations ice nucleation caused by the direct formation of ice crystals upon ice-nucleating particles, as well as the Hallett-Mossop process, has a negligible contribution to the ice crystal concentration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Study of ice cluster impacts on amorphous silica using the ReaxFF reactive force field molecular dynamics simulation method

We study the dynamics of the collisions between amorphous silica structures and amorphous and crystal ice clusters with impact velocities of 1 km/s, 4 km/s, and 7 km/s using the ReaxFF reactive molecular dynamics simulation method. The initial ice clusters consist of 150 water molecules for the amorphous ice cluster and 128 water molecules for the crystal ice cluster. The ice clusters are collided on the surface of amorphous fully oxidized and suboxide silica. These simulations show that at 1 km/s impact velocities, all the ice clusters accumulate on the surface and at 4 km/s and 7 km/s impact velocities, some of the ice cluster molecules bounce back from the surface. At 4 km/s and 7 km/s impact velocities, few of the water molecules dissociations are observed. The effect of the second ice cluster impacts on the surfaces which are fully covered with ice, on the mass loss/accumulation is studied. These studies show that at 1 km/s impacts, the entire ice cluster accumulates on the surface at both first and second ice impacts. At higher impact velocities, some ice molecules which after the first ice impacts have been attached to the surface will separate from the surface after the second ice impacts at 7 km/s impact velocity. For the 4 km/s ice cluster impact, ice accumulation is observed for the crystal ice cluster impacts and ice separation is observed for the amorphous ice impacts. Observing the temperatures of the ice clusters during the collisions indicates that the possibility of electron excitement at impact velocities less than 10 km/s is minimal and ReaxFF reactive molecular dynamics simulation can predict the chemistry of these hypervelocity impacts. However, at impact velocities close to 10 km/s the average temperature of the impacting ice clusters increase to about 2000 K, with individual molecules occasionally reaching temperatures of over 8000 K and thus it will be prudent to consider the concept of electron excitation at these higher impact velocities, which goes beyond the current ReaxFF ability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Molecular-Resolution Electron Imaging of Defects and Dynamics at the Ice-Water Interface

Water crystallization into hexagonal ice (type I h ) is one of the most critical processes relevant to the Earth’s environment and human activities. However, despite recent breakthroughs in imaging non-equilibrium condensed ice structures, the ice-water interface has never been imaged at a molecular resolution. This is primarily due to the low stability of the hydrogen bonds in ice under high-resolution microscopy conditions and a lack of methods to prepare compatible samples. Here, this presentation describes the first molecular-resolution imaging of ice crystallized from liquid water and the ice-water interface using high-resolution transmission electron microscopy (HRTEM). By encapsulating deionized (DI) water between two amorphous carbon (a-C) TEM grids and subsequently freezing it with liquid N 2 on a cryo sample stage, we generated two types of ice: non-equilibrium, condensed ice from the atmosphere and encapsulated ice from the DI water (Fig. 1A). Condensed ice usually shows irregular, spherulitic shapes (Fig. 1B). Selected area electron diffraction (SAED) shows that they are a mixture of cubic and hexagonal crystals (Fig. 1C). On the contrary, encapsulated ice forms thin films that contain large-area single-crystalline regions of hexagonal ice oriented along the [0001] zone axis (Fig. 1D). Differential electron energy-loss spectroscopy (EELS) confirmed the high purity of the encapsulated ice free from organic contaminations that are common in other encapsulation methods for HRTEM such as graphene liquid cells. These single-crystalline areas are robust under the electron beam up to ~100 e/Å 2 s. Aberration-corrected HRTEM imaging in these areas achieved a line resolution of ~1.3 Å (Fig. 1E and F). This platform allows us to study near-equilibrium ice structures and dynamics at an unprecedented spatial resolution (Fig. 2). For example, we discovered subdomain-rich regions near the defective crystal edges despite the structure appearing single-crystalline according to diffraction criteria. These subdomains connect via low-angle grain boundaries with flat energy landscapes as a function of tilt angles (according to simulations), showing the high tolerance of ice to defect structures. When we tuned the sample temperature and electron flux rate, we observed radiolysis-controlled bubble generation and dissolution in ice single crystals near a steady state of bubble dynamics. Furthermore, rich beam-induced melting and recrystallization dynamics were observed at the ice-water interface with lattice resolution. These data represent the first observation of the ice-water phase transformation at the sub-nanometer level. In summary, the methods developed in this work enabled molecular-resolution observations of ice and the ice-water interface and shed light on the microstructures and phase transformation pathways. Finer control on the temperature, electron irradiation profile, and imaging detector could eventually lead to real-time observation of ice nucleation in water and address long-standing questions in the nucleation pathways.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

Final scientific report for DOE award title: Improving the Representation of Ice Sedimentation Rates in Global Climate Models

It is well known that cirrus clouds play a major role in regulating the earth’s climate, but the details of how this works are just beginning to be understood. This project targeted the main property of cirrus clouds that influence climate processes; the ice fall speed. That is, this project improves the representation of the mass-weighted ice particle fall velocity, V m , in climate models, used to predict future climate on global and regional scales. Prior to 2007, the dominant sizes of ice particles in cirrus clouds were poorly understood, making it virtually impossible to predict how cirrus clouds interact with sunlight and thermal radiation. Due to several studies investigating the performance of optical probes used to measure the ice particle size distribution (PSD), as well as the remote sensing results from our last ARM project, it is now well established that the anomalously high concentrations of small ice crystals often reported prior to 2007 were measurement artifacts. Advances in the design and data processing of optical probes have greatly reduced these ice artifacts that resulted from the shattering of ice particles on the probe tips and/or inlet tube, and PSD measurements from one of these improved probes (the 2-dimensional Stereo or 2D-S probe) are utilized in this project to parameterize V m for climate models. Our original plan in the proposal was to parameterize the ice PSD (in terms of temperature and ice water content) and ice particle mass and projected area (in terms of mass- and area-dimensional power laws or m-D/A-D expressions) since these are the microphysical properties that determine V m , and then proceed to calculate V m from these parameterized properties. But the 2D-S probe directly measures ice particle projected area and indirectly estimates ice particle mass for each size bin. It soon became apparent that the original plan would introduce more uncertainty in the V m calculations than simply using the 2D-S measurements to directly calculate V m . By calculating V m directly from the measured PSD, ice particle projected area and estimated mass, more accurate estimates of V m are obtained. These V m values were then parameterized for climate models by relating them to (1) sampling temperature and ice water content (IWC) and (2) the effective diameter (D e ) of the ice PSD. Parameterization (1) is appropriate for climate models having single-moment microphysical schemes whereas (2) is appropriate for double-moment microphysical schemes and yields more accurate V m estimates. These parameterizations were developed for tropical cirrus clouds, Arctic cirrus, mid-latitude synoptic cirrus and mid-latitude anvil cirrus clouds based on field campaigns in these regions. An important but unexpected result of this research was the discovery of microphysical evidence indicating the mechanisms by which ice crystals are produced in cirrus clouds. This evidence, derived from PSD measurements, indicates that homogeneous freezing ice nucleation dominates in mid-latitude synoptic cirrus clouds, whereas heterogeneous ice nucleation processes dominate in mid-latitude anvil cirrus. Based on these findings, D e was parameterized in terms of temperature (T) for conditions dominated by (1) homo- and (2) heterogeneous ice nucleation. From this, an experiment was designed for global climate models (GCMs). The net radiative forcing from cirrus clouds may be affected by the means ice is produced (homo- or heterogeneously), and this net forcing contributes to climate sensitivity (i.e. the change in mean global surface temperature resulting from a doubling of CO 2 ). The objective of this GCM experiment was to determine how a change in ice nucleation mode affects the predicted global radiation balance. In the first simulation (Run 1), the D e -T relationship for homogeneous nucleation is used at all latitudes, while in the second simulation (Run 2), the D e -T relationship for heterogeneous nucleation is used at all latitudes. For both runs, V m is calculated from D e . Two GCMs were used; the Community Atmosphere Model version 5 (CAM5) and a European GCM known as ECHAM5 (thanks to our European colleagues who collaborated with us). Similar results were obtained from both GCMs in the Northern Hemisphere mid-latitudes, with a net cooling of ~ 1.0 W m -2 due to heterogeneous nucleation, relative to Run 1. The mean global net cooling was 2.4 W m -2 for the ECHAM5 GCM while CAM5 produced a mean global net cooling of about 0.8 W m -2 . This dependence of the radiation balance on nucleation mode is substantial when one considers the direct radiative forcing from a CO 2 doubling is 4 W m -2 . The differences between GCMs in mean global net cooling estimates may demonstrate a need for improving the representation of cirrus clouds in GCMs, including the coupling between microphysical and radiative properties. Unfortunately, after completing this GCM experiment, we learned from the company that provided the 2D-S microphysical data that the data was corrupted due to a computer program coding problem. Therefore the microphysical data had to be reprocessed and reanalyzed, and the GCM experiments were redone under our current ASR project but using an improved experimental design.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Structural and stratigraphic features and ERS 1 synthetic aperture radar backscatter characteristics of ice growing on shallow lakes in NW Alaska, winter 1991-1992

Changes in Earth Remote-Sensing Satellite (ERS) 1 C band synthetic aperture radar (SAR) backscatter intensity (sigma(exp 0)) from ice growing on shallow tundra lakes at three locations in NW Alaska are described. Ice core analysis shows that all lakes on the coast at Barrow the ice, whether floating or frozen to the bottom, includes an inclusion-free layer overlying a layer of ice with tubular bubbles oriented parallel to the direction of growth. The clear ice may also be overlain by a discontinuous layer of bubbly snow ice. Backscatter is low (-16 to -22 dB) at the time of initial ice formation, probably due to the specular nature of the upper and lower ice surfaces causing the radar pulse to be reflected away from the radar. As the ice thickens during the autumn, backscatter rises steadily. Once the ice freezes to the lake bottom, regardless of the presence of foward scattering tubular bubbles, low backscatter values of -17 to -18 dB are caused by absorption of the radar signal in the lake bed. For ice that remains afloat all winter the ice-water interface and the tubular bubbles combine, presumably via an incoherent double-bounce mechanism, to cause maximum backscatter values of the order of -6 to -7 dB. The sigma(exp 0) saturates at -6 to -7 dB before maximum ice thickness and tubular bubble content are attained. A simple ice growth model suggests that the layer of ice with tubular bubbles need be only a few centimeters thick midway through the growth season to cause maximum backscatter from floating ice. During the spring thaw a previously unreported backscatter reversal is observed on the floating and grounded portions of the coastal lakes but not on the lakes farther inland. This reversal may be related to the ice surface topography and wetness plus the effects of a longer, cooler melt period by the coast. Time series of backscatter variations from shallow tundra lakes are a record of (1) the development of tubular bubbles in the ice and, by association, changes in the gas content of the underlying water and (2) the freezing of ice to the bottoms of the lakes and therefore lake bathymetry and water availability. SAR is also able to detect the onset of lake ice growth in autumn and the initiation of the spring thaw and thus has potential for monitoring high-altitude lake ice growth and decay processes in relation to climate variability.

Jeffries, M. O.↗
NASA NTRS1996

Seasonal Characteristics of the Perennial Ice Cover of the Beaufort Sea

By definition, ice which survives the summer is classified as multiyear ice. Thus the area covered by multiyear ice during the winter should be nearly equivalent to the ice area during the previous summer's minima. This condition provides a reasonable criterion for the evaluation of ice concentration and ice type retrieval algorithms using remote-sensing data sets. From special sensor microwave imager (SSM/I) data the NASA Team algorithm estimates the multiyear, first-year, and total ice concentrations during the winter using combinations of the polarization and spectral gradient ratios. The Team algorithm provides only estimates of ice concentration in the summer. From ERS 1 synthetic aperture radar (SAR) data the remarkably stable contrast between multiyear ice and first-year ice in winter provides consistent estimates of multiyear ice concentrations. In the summer, multiyear ice concentration cannot be estimated from SAR or SSM/I data because free water on the surface effectively masks the backscatter and emissivity signature of this ice type. From SAR data a technique which takes advantage of the high backscatter of wind-roughened open water as a discrimination feature is used to estimate the total ice concentration in the summer. With a year-long (January 1992 to January 1993) data set from the Beaufort Sea we found that the multiyear ice concentration estimates from the SAR data are stable and are nearly equivalent to the ice concentration estimated at the end of the previous summer. We contrast this with the variability of the multiyear ice concentration and ice friction estimates obtained using SSM/I data. The Team algorithm produces ice concentration and multiyear ice estimates which are consistently lower than those from the SAR data. We discuss reasons for these discrepancies and the implications of the higher than previously noted multiyear ice concentrations.

Kwok, R.↗
NASA NTRS2009

The Satellite Passive-Microwave Record of Sea Ice in the Ross Sea Since Late 1978

Satellites have provided us with a remarkable ability to monitor many aspects of the globe day-in and day-out and sea ice is one of numerous variables that by now have quite substantial satellite records. Passive-microwave data have been particularly valuable in sea ice monitoring, with a record that extends back to August 1987 on daily basis (for most of the period), to November 1970 on a less complete basis (again for most of the period), and to December 1972 on a less complete basis. For the period since November 1970, Ross Sea sea ice imagery is available at spatial resolution of approximately 25 km. This allows good depictions of the seasonal advance and retreat of the ice cover each year, along with its marked interannual variability. The Ross Sea ice extent typically reaches a minimum of approximately 0.7 x 10(exp 6) square kilometers in February, rising to a maximum of approximately 4.0 x 10(exp 6) square kilometers in September, with much variability among years for both those numbers. The Ross Sea images show clearly the day-by-day activity greatly from year to year. Animations of the data help to highlight the dynamic nature of the Ross Sea ice cover. The satellite data also allow calculation of trends in the ice cover over the period of the satellite record. Using linear least-squares fits, the Ross Sea ice extent increased at an average rate of 12,600 plus or minus 1,800 square kilometers per year between November 1978 and December 2007, with every month exhibiting increased ice extent and the rates of increase ranging from a low of 7,500 plus or minus 5,000 square kilometers per year for the February ice extents to a high of 20,300 plus or minus 6,100 kilometers per year for the October ice extents. On a yearly average basis, for 1979-2007 the Ross Sea ice extent increased at a rate of 4.8 plus or minus 1.6 % per decade. Placing the Ross Sea in the context of the Southern Ocean as a whole, over the November 1978-December 2007 period the Ross Sea had the highest rate of increase in sea ice coverage of any of five standard divisions of the Southern Ocean, although the Weddell Sea, Indian Ocean, and Western Pacific Ocean all also had sea ice increases, while only the Bellingshausen/Smundsen Seas experienced overall sea ice decreases. Overall, the Southern Ocean sea ice cover increased at an average rate of 10,800 plus or minus 2,500 square kilometers per year between November 1978 and December 2007, with every month showing positive values although with some of these values not being statistically significant. The sea ice increase since November 1978 was preceded by a sharp decrease in Southern Ocean ice coverage in the 1970's and is in marked contrast to the decrease in Arctic sea ice coverage that has occurred both in the period since November 1978 and since earlier in the 1970's. On a yearly average bases, for 1979-2007 the Southern Ocean sea ice extent increased at a rate of 1.0 plus or minus 0.4% per decade, whereas the Arctic ice extent decreased at the much greater rate of 4.0 plus or minus 0.4 percent per decade (closer to the % per decade rate of increase in the Ross Sea). Considerable research is ongoing to explain the differences.

Parkinson, Claire L.↗
NASA NTRSDate not supplied

Computational Icing Analysis on NASA’s SIDRM Geometry to Investigate Collection Efficiency

Computational icing analysis results were compared to experimental icing tunnel data including aerothermal (e.g. dry air) and supercooled water droplet rime-ice conditions from tests conducted in early 2022 at the NASA Icing Research Tunnel (IRT). The Simulated Intercompressor Duct Research Model (SIDRM) test article was used in this study, and its geometry represents the inter-compressor duct region of a turbofan engine. The test article’s purpose is to study the physics of supercooled water icing and ice crystal icing. This study compared three different icing codes: FENSAP-ICE (Eulerian approach), LEWICE3D (Lagrangian approach), and GlennICE (Lagrangian approach). All three icing codes were conducted on SIDRM’s complex body flow-field and compared to different experimental supercooled water rime runs. The test article instrumentation (pressure taps, thermocouples, etc.) and 3D laser scans of final ice shapes were used to compare against the different icing code simulations. The overall objectives are to understand how the icing codes handle capturing collection efficiency on the complex test article’s unheated surfaces. In the aerothermal cases, pressure tap readings matched the CFD results, but dry air CFD underpredicted thermocouple readings. Collection efficiency results from all three icing codes matched well together on the main body leading edge, main body slope, and the strut leading edges. All three icing codes underpredicted the final ice shape using a single-shot constant ice density approach, with more difficulty coming from the strut leading edge ice shape due to the swept wing like flow field. NASA’s overall goal for this effort is to develop computational icing tools to assist in the design and certification of engines for flight in icing conditions.

Supercooled icing, SIDRM, aerodynamics, computatio↗
DOE OSTI2016

A surface ice module for wind turbine dynamic response simulation using FAST

It is a fact that developing offshore wind energy has become more and more serious worldwide in recent years. Many of the promising offshore wind farm locations are in cold regions that may have ice cover during wintertime. The challenge of possible ice loads on offshore wind turbines raises the demand of modeling capacity of dynamic wind turbine response under the joint action of ice, wind, wave, and current. The simulation software FAST is an open source computer-aided engineering (CAE) package maintained by the National Renewable Energy Laboratory. In this paper, a new module of FAST for assessing the dynamic response of offshore wind turbines subjected to ice forcing is presented. In the ice module, several models are presented which involve both prescribed forcing and coupled response. For conditions in which the ice forcing is essentially decoupled from the structural response, ice forces are established from existing models for brittle and ductile ice failure. For conditions in which the ice failure and the structural response are coupled, such as lock-in conditions, a rate-dependent ice model is described, which is developed in conjunction with a new modularization framework for FAST. In this paper, analytical ice mechanics models are presented that incorporate ice floe forcing, deformation, and failure. For lower speeds, forces slowly build until the ice strength is reached and ice fails resulting in a quasi-static condition. For intermediate speeds, the ice failure can be coupled with the structural response and resulting in coinciding periods of the ice failure and the structural response. A third regime occurs at high speeds of encounter in which brittle fracturing of the ice feature occurs in a random pattern, which results in a random vibration excitation of the structure. An example wind turbine response is simulated under ice loading of each of the presented models. This module adds to FAST the capabilities for analyzing the response of wind turbines subjected to forces resulting from ice impact on the turbine support structure. The conditions considered in this module are specifically addressed in the International Organization for Standardization (ISO) standard 19906:2010 for arctic offshore structures design consideration. Special consideration of lock-in vibrations is required due to the detrimental effects of such response with regard to fatigue and foundation/soil response. Finally, the use of FAST for transient, time domain simulation with the new ice module is well suited for such analyses.

17 WIND ENERGY↗
DOE OSTI2024

Approximation of ice phenology of Maine lakes using Aqua MODIS surface temperature data

Studies of lake ice phenology have historically relied on limited in situ data. Relatively few observations exist for ice out and fewer still for ice in, both of which are necessary to determine the temporal extent of ice cover. Satellite data provide an opportunity to better document patterns of ice phenology across landscapes and relate them to the climatological drivers behind changing ice phenology. We developed a model, the Cumulative Sum Method (CSM), that uses daytime and nighttime surface temperature observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on board the Earth-observing Aqua satellite to approximate ice in (the onset of ice cover) and ice out from training datasets of 13 and 58 Maine lakes, respectively, during the 2002/2003 through 2017/2018 ice seasons. Ice in was signaled by reaching a threshold of cumulative negative degrees following the first day of the season below 0°C. Ice out was signaled by reaching a threshold of cumulative positive degrees following the first day of the year above 0°C. The comparison of observed and remotely sensed ice-in dates showed relative agreement with a correlation coefficient of 0.71 and a mean absolute error (MAE) of 9.8 days. Ice-out approximations had a correlation coefficient of 0.67 and an MAE of 8.8 days. Lakes smaller in surface area and nearer the Atlantic coast had the greatest error in approximation. Application of the CSM to 20 additional lakes in Maine produced a comparable ice-out MAE of 8.9 days. Ice-out model performance was weaker for the warmest years; there was a larger MAE of 12.0 days when the model was applied to the years 2019–2023 for the original 58 lakes. The development of this model, which utilizes daily satellite data, demonstrates the promise of remote sensing for quantifying ice phenology over short, temporal scales, and wider geographic regions than can be observed in situ, and allows exploration of the influence of surface temperature patterns on the process and timing of ice in and ice out.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

Analysis of the Honeywell Uncertified Research Engine with Ice Crystal Cloud Ingestion at Simulated Altitudes

The Honeywell Uncertified Research Engine (HURE), a research version of a turbofan engine that never entered production, was tested in the NASA Propulsion System Laboratory (PSL), an altitude test facility at the NASA Glenn Research Center. The PSL is a facility that is equipped with water spray bars capable of producing an ice cloud consisting of ice particles, having a controlled particle diameter and concentration in the air flow. To develop the test matrix of the HURE, numerical analysis of flow and ice particle thermodynamics was performed on the compression system of the turbofan engine to predict operating conditions that could potentially result in a risk of ice accretion due to ice crystal ingestion. The goal of the test matrix was to provide operating conditions such that ice would accrete in either the fan-stator through the inlet guide vane region of the compression system or within the first stator of the high pressure compressor. The predictive analyses were performed with the mean line compressor flow modeling code (COMDESMELT) which includes an ice particle model. The HURE engine was tested in PSL with the ice cloud over the range of operating conditions of altitude, ambient temperature, simulated flight Mach number, and fan speed with guidance from the analytical predictions. The engine was fitted with video cameras at strategic locations within the engine compression system flow path where ice was predicted to accrete, in order to visually confirm ice accretion when it occurred. In addition, traditional compressor instrumentation such as total pressure and temperature probes, static pressure taps, and metal temperature thermocouples were installed in targeted areas where the risk of ice accretion was expected. The current research focuses on the analysis of the data that was obtained after testing the HURE engine in PSL with ice crystal ingestion. The computational method (COMDES-MELT) was enhanced by computing key parameters through the fan-stator at multiple span wise locations, in order to increase the fidelity with the current mean-line method. The Icing Wedge static wet-bulb temperature thresholds were applied for determining the risk of ice accretion in the fan-stator, which is thought to be an adiabatic region. At some operating conditions near the splitter-lip region, other sources of heat (non-adiabatic walls) were suspected to be the cause of accretion, and the Icing Wedge was not applied to predict accretion at that location. A simple order-of-magnitude heat transfer model was implemented into the COMDES-MELT code to estimate the wall temperature minimum and maximum thresholds that support ice accretion, as observed by video confirmation. The results from this model spanned the range of wall temperatures measured on a previous engine that experienced ice accretion at certain operating conditions. The goal of this study is to show that the computational process developed on earlier engine icing tests can be used to provide an icing risk assessment in adiabatic regions for other engines.

Turbofan engine↗
NASA NTRS1993

Polarimetric backscattering from sea ice

Polarimetric backscattering from sea ice is presented. Theoretical models for backscattering are first presented for various ice types. Then, theoretical results are compared with experimental data for new thin ice, first-year ice, and multi-year ice. Sea ice is modeled as a layer medium containing random scatterers and rough interfaces. For multi-year sea ice with snow cover, the sea ice layer is modeled as an ice background with embedded spheroidal air bubbles and the snow layer as air with ice grains. The hummocky topography on multi-year ice is characterized by a Gaussian distribution which has an averaging effect on backscattering coefficients. First-year sea ice is described by an ice medium hosting ellipsoidal brine inclusions. These inclusions are oriented preferentially in the vertical direction due to the columnar structure of first-year sea ice. Azimuthally, the orientation of the brine pockets are random corresponding to the random c-axes in sea ice, unless the axes are oriented by sea currents. For thin ice in newly opened leads, it has been observed that there exists a thin brine layer with very high salinity on the top surface of the new ice. This brine layer is depicted as a medium with high permittivity which can significantly affect electromagnetic scattering signatures from the lower thin ice layer, with a higher fractional volume of brine inclusions due to higher salinity as compared to thick first-year sea ice. The rough medium interfaces are described as Gaussian rough surfaces characterized by root-mean-square heights and surface correlation lengths. The contribution from rough surface, calculated under the Kirchhoff approximation or small perturbation method, is assumed to be independent from volume scattering. The total loss including absorption and scattering losses in the scattering media is represented by the imaginary part of effective permittivity obtained from the strong fluctuation theory. The polarimetric scattering coefficients for different ice types are then derived under the distorted Born approximation.

Nghiem, S. V.↗
NASA NTRS2012

Simulation Tools Model Icing for Aircraft Design

Here s a simple science experiment to try: Place an unopened bottle of distilled water in your freezer. After 2-3 hours, if the water is pure enough, you will notice that it has not frozen. Carefully pour the water into a bowl with a piece of ice in it. When it strikes the ice, the water will instantly freeze. One of the most basic and commonly known scientific facts is that water freezes at around 32 F. But this is not always the case. Water lacking any impurities for ice crystals to form around can be supercooled to even lower temperatures without freezing. High in the atmosphere, water droplets can achieve this delicate, supercooled state. When a plane flies through clouds containing these droplets, the water can strike the airframe and, like the supercooled water hitting the ice in the experiment above, freeze instantly. The ice buildup alters the aerodynamics of the plane - reducing lift and increasing drag - affecting its performance and presenting a safety issue if the plane can no longer fly effectively. In certain circumstances, ice can form inside aircraft engines, another potential hazard. NASA has long studied ways of detecting and countering atmospheric icing conditions as part of the Agency s efforts to enhance aviation safety. To do this, the Icing Branch at Glenn Research Center utilizes a number of world-class tools, including the Center s Icing Research Tunnel and the NASA 607 icing research aircraft, a "flying laboratory" for studying icing conditions. The branch has also developed a suite of software programs to help aircraft and icing protection system designers understand the behavior of ice accumulation on various surfaces and in various conditions. One of these innovations is the LEWICE ice accretion simulation software. Initially developed in the 1980s (when Glenn was known as Lewis Research Center), LEWICE has become one of the most widely used tools in icing research and aircraft design and certification. LEWICE has been transformed over the years from strictly a research tool to one used routinely by industry and other government agencies. Glenn contractor William Wright has been the architect of this development, supported by a team of researchers investigating icing physics, creating validation data, and ensuring development according to standard software engineering practices. The program provides a virtual simulation environment for determining where water droplets strike an airfoil in flight, what kind of ice would result, and what shape that ice would take. Users can enter geometries for specific, two-dimensional cross sections of an airfoil or other airframe surface and then apply a range of inputs - different droplet sizes, temperatures, airspeeds, and more - to model how ice would build up on the surface in various conditions. The program s versatility, ease of use, and speed - LEWICE can run through complex icing simulations in only a few minutes - have contributed to it becoming a popular resource in the aviation industry.

Source record↗
NASA NTRSDate not supplied

Estimating the Arctic Cloud-Sea Ice Feedback With Observations during the EOS Period

Arctic sea ice responds to and drives Arctic climate change. The interactions between Arctic sea ice and clouds represent a mechanism through which sea ice can drive climate change. We composite active remote sensing satellite cloud properties for ice-free, marginal ice zone (MIZ), and ice-covered surfaces during MIZ crossing events to investigate the influence of the transition from an ice-covered to an ice-free surface on low-level clouds. We demonstrate that the event-based methodology controls for large-scale meteorological factors and isolates the sea ice effect on clouds. We find larger cloud fraction and total water content below ~1.5 km over ice-free relative to ice-covered surfaces during non-summer months, indicating a low-level cloud sensitivity to Arctic sea ice decline. During summer, results show larger cloud fraction and water content over ice-free surfaces, however the differences are statistically indistinguishable. Evidence is provided that atmospheric thermodynamic profile differences cause the cloud property differences, namely that ice-free footprints are warmer, moister, have more positive surface turbulent fluxes and are less stable than their ice-covered counterparts. Ice-free and ice-covered surface cloud property differences scale with surface temperature differences such that cloud property differences are only found in the presence of a surface temperature difference. We conclude that surface temperature differences modulate the cloud response to sea ice loss through influences on surface turbulent fluxes and lower tropospheric stability. The results imply a positive non-summer sea ice-cloud feedback and that up to a 0.02 cloud fraction and 0.05 g m 3 total water content increase in fall are due to the observed sea ice decline.

Patrick C. Taylor↗
NASA NTRSDate not supplied

Simultaneous Ice Water Content Measurements at Multiple Locations on the NASA DC-8 Aircraft during the 2018 HIWC RADAR Flight Campaign

Ice water content measurements were made simultaneously at three locations on the NASA DC-8 in natural, glaciated conditions during the 2018 High Ice Water Content RADAR flight campaign. The purpose of these measurements was to further evaluate efficiency factors of hot-wire total water content probes in glaciated conditions and investigate the enhancement of ice crystal concentrations near fuselage surfaces due to flow field inertial effects, and ice crystals impacting the nose, breaking up, and flowing downstream. The total water content measurements were made using either Science Engineering Associates Ice Crystal Detectors or Robust Probes. Three common sensors were mounted on an underwing canister considered to be in near free-flow conditions, a standoff from a fuselage window, and the nose of the fuselage near the pitot probes. The Ice Crystal Detector total water content sensor and Robust Probe sensor collection and retention efficiencies were evaluated through comparisons with the underwing Ice Crystal Detector and Robust Probe measurements to the IKP2 isokinetic evaporator probe, which provided the reference ice water content measurement. Local ice water content at the nose position was evaluated by comparing ratios of the nose and underwing ice crystal detectors to the IKP2. Local ice water content at the window-standoff location was also evaluated by comparing total water content sensor measurements from the window ice crystal detector to the measurements made with the underwing and nose Ice Crystal Detectors The key findings were: 1) the Ice Crystal Detector total water content sensor ice water content efficiency factor was similar to previous estimates, but reduced with increased ice crystal median mass diameter; 2) the ice water content at the fuselage nose location near the DC-8 pitot probes was approximately 2.5 times the freestream values – although this estimate is affected by a higher probe efficiency factor due to smaller particles in the debris cloud from impacts upstream of the probe; and 3) the ice water content at the 17” standoff from the port window varied from about 50% to 3 times freestream values in a complicated manner. Similar measurement locations are not uncommon on cloud research aircraft, where ice particle measurements may be subject to similar uncertainties.

Aircraft Icing↗