Engineering Papers⌕ Search

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.

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

At least 433 records · Page 24

Snowfall and snow accumulation during the MOSAiC winter and spring seasons

Data from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition allowed us to investigate the temporal dynamics of snowfall, snow accumulation and erosion in great detail for almost the whole accumulation season (November 2019 to May 2020). We computed cumulative snow water equivalent (SWE) over the sea ice based on snow depth and density retrievals from a SnowMicroPen and approximately weekly measured snow depths along fixed transect paths. We used the derived SWE from the snow cover to compare with precipitation sensors installed during MOSAiC. The data were also compared with ERA5 reanalysis snowfall rates for the drift track. We found an accumulated snow mass of 38 mm SWE between the end of October 2019 and end of April 2020. The initial SWE over first-year ice relative to second-year ice increased from 50 % to 90 % by end of the investigation period. Further, we found that the Vaisala Present Weather Detector 22, an optical precipitation sensor, and installed on a railing on the top deck of research vessel Polarstern, was least affected by blowing snow and showed good agreements with SWE retrievals along the transect. On the contrary, the OTT Pluvio 2 pluviometer and the OTT Parsivel 2 laser disdrometer were largely affected by wind and blowing snow, leading to too high measured precipitation rates. These are largely reduced when eliminating drifting snow periods in the comparison. ERA5 reveals good timing of the snowfall events and good agreement with ground measurements with an overestimation tendency. Retrieved snowfall from the ship-based K a -band ARM zenith radar shows good agreements with SWE of the snow cover and differences comparable to those of ERA5. Based on the results, we suggest the K a -band radar-derived snowfall as an upper limit and the present weather detector on RV Polarstern as a lower limit of a cumulative snowfall range. Based on these findings, we suggest a cumulative snowfall of 72 to 107 mm and a precipitation mass loss of the snow cover due to erosion and sublimation as between 47 % and 68 %, for the time period between 31 October 2019 and 26 April 2020. Extending this period beyond available snow cover measurements, we suggest a cumulative snowfall of 98–114 mm.

54 ENVIRONMENTAL SCIENCES↗

The terminal area simulation system. Volume 1: Theoretical formulation

A three-dimensional numerical cloud model was developed for the general purpose of studying convective phenomena. The model utilizes a time splitting integration procedure in the numerical solution of the compressible nonhydrostatic primitive equations. Turbulence closure is achieved by a conventional first-order diagnostic approximation. Open lateral boundaries are incorporated which minimize wave reflection and which do not induce domain-wide mass trends. Microphysical processes are governed by prognostic equations for potential temperature water vapor, cloud droplets, ice crystals, rain, snow, and hail. Microphysical interactions are computed by numerous Orville-type parameterizations. A diagnostic surface boundary layer is parameterized assuming Monin-Obukhov similarity theory. The governing equation set is approximated on a staggered three-dimensional grid with quadratic-conservative central space differencing. Time differencing is approximated by the second-order Adams-Bashforth method. The vertical grid spacing may be either linear or stretched. The model domain may translate along with a convective cell, even at variable speeds.

Proctor, F. H.↗

Cloud Detection Using Measured and Modeled State Parameters

In this study, hourly RUC analyses were used to examine the differences between RH and temperature values from RUC reanalysis data and from radiosonde atmospheric profiles obtained at the ARM SCF. The results show that the temperature observations from the SONDE and RUC are highly correlated. The RHs are also well-correlated, but the SONDE values generally exceed those from RUC. Inside cloud layers, the RH from RUC is 2-14% lower than the RH from SONDE for all RUC layers. Although the layer mean RH within clouds is much greater than the layer mean RH outside cloud or in the clear-sky, RH thresholds chosen as a function of temperature can more accurately diagnose cloud occurrence for either dataset. For overcast clouds, it was found that the 50% probability RH threshold for diagnosing a cloud, within a given upper tropospheric layer is roughly 90% for the Vaisala RS80-15LH radisonde and 80% for RUC data. While for the partial cloud (cloud amount is less than 90%), the RH thresholds of SONDE are close to RUC for a given probability in upper tropospheric layers. The probabilities of detecting clouds at a given RH and temperature should be useful for a variety of application such as the development of new cloud parameterizations or for estimating the vertical profile of cloudiness underneath a given cloud observed from the satellite to construct a 3-D cloud data set for computing atmospheric radiative heating profiles or determining potential aircraft icing conditions.

Yi, Y.↗

NASA Tech Briefs, March 2007

Topics include: Advanced Systems for Monitoring Underwater Sounds; Wireless Data-Acquisition System for Testing Rocket Engines; Processing Raw HST Data With Up-to-Date Calibration Data; Mobile Collection and Automated Interpretation of EEG Data; System for Secure Integration of Aviation Data; Servomotor and Controller Having Large Dynamic Range; Digital Multicasting of Multiple Audio Streams; Translator for Optimizing Fluid-Handling Components; AIRSAR Web-Based Data Processing; Pattern Matcher for Trees Constructed From Lists; Reducing a Knowledge-Base Search Space When Data Are Missing; Ground-Based Correction of Remote-Sensing Spectral Imagery; State-Chart Autocoder; Pointing History Engine for the Spitzer Space Telescope; Low-Friction, High-Stiffness Joint for Uniaxial Load Cell; Magnet-Based System for Docking of Miniature Spacecraft; Electromechanically Actuated Valve for Controlling Flow Rate; Plumbing Fixture for a Microfluidic Cartridge; Camera Mount for a Head-Up Display; Core-Cutoff Tool; Recirculation of Laser Power in an Atomic Fountain; Simplified Generation of High-Angular-Momentum Light Beams; Imaging Spectrometer on a Chip; Interferometric Quantum-Nondemolition Single-Photon Detectors; Ring-Down Spectroscopy for Characterizing a CW Raman Laser; Complex Type-II Interband Cascade MQW Photodetectors; Single-Point Access to Data Distributed on Many Processors; Estimating Dust and Water Ice Content of the Martian Atmosphere From THEMIS Data; Computing a Stability Spectrum by Use of the HHT; Theoretical Studies of Routes to Synthesis of Tetrahedral N4; Estimation Filter for Alignment of the Spitzer Space Telescope; Antenna for Measuring Electric Fields Within the Inner Heliosphere; Improved High-Voltage Gas Isolator for Ion Thruster; and Hybrid Mobile Communication Networks for Planetary Exploration.

Source record↗

Global Sensitivity Analysis Using the Ultra‐Low Resolution Energy Exascale Earth System Model

Abstract For decades, Arctic temperatures have increased twice as fast as average global temperatures. As a first step toward quantifying parametric uncertainty in Arctic climate, we performed a variance‐based global sensitivity analysis (GSA) using a fully coupled, ultra‐low resolution (ULR) configuration of version 1 of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SMv1). Specifically, we quantified the sensitivity of six quantities of interests (QOIs), which characterize changes in Arctic climate over a 75 year period, to uncertainties in nine model parameters spanning the sea ice, atmosphere, and ocean components of E3SMv1. Sensitivity indices for each QOI were computed with a Gaussian process emulator using 139 random realizations of the random parameters and fixed preindustrial forcing. Uncertainties in the atmospheric parameters in the Cloud Layers Unified by Binormals (CLUBB) scheme were found to have the most impact on sea ice status and the larger Arctic climate. Our results demonstrate the importance of conducting sensitivity analyses with fully coupled climate models. The ULR configuration makes such studies computationally feasible today due to its low computational cost. When advances in computational power and modeling algorithms enable the tractable use of higher‐resolution models, our results will provide a baseline that can quantify the impact of model resolution on the accuracy of sensitivity indices. Moreover, the confidence intervals provided by our study, which we used to quantify the impact of the number of model evaluations on the accuracy of sensitivity estimates, have the potential to inform the computational resources needed for future sensitivity studies.

54 ENVIRONMENTAL SCIENCES↗

Computation of Low Speed Cavity Noise

Over the last five years, the Aircraft Icing Project of the NASA Aviation Safety Program has developed a number of in-flight icing education and training aids to support increased awareness for pilots of the hazards associated with atmospheric icing conditions. Through the development of this work, a number of new instructional design approaches and media delivery methods have been introduced to enhance the learning experience, expand user interactivity and participation, and, hopefully, increase the learner retention rates. The goal of using these multimedia techniques is to increase the effectiveness of the training materials. This paper will describe the multimedia technology that has been introduced and give examples of how it was used.

Blech, Richard A.↗

Microphysical Characterization of Ice in Hurricane Erin for Wideband Passive Microwave Comparisons

In order to better understand the characteristics of frozen hydrometeors in hurricane systems, computed brightness temperatures were compared with observations of Hurricane Erin from CAMEX-4 radiometers on the ER-2 aircraft. These wideband observations ranged from 10.7 to 183 plus or minus 10 GHz. In order to initialize the cloud profiles used in the radiative transfer calculations, radar reflectivities from an ER-2 instrument were converted to mass contents at specified height levels. Dropsonde data from the ER-2 and modeled Hurricane Andrew data were used to generate profiles of temperature, pressure, and relative humidity and to provide surface wind speed and ocean temperature. Four different ice (size and density) parameterizations were used to compute brightness temperature values. The comparisons between the calculated and observed brightness temperatures show that in general the Grenfell and Warren particle size distribution does the best job of matching the observations for most of the high frequencies. These high frequencies (greater than or equal to 85 GHz) are particularly sensitive to frozen particles. The Sehkon and Srivastava particle size distribution produced brightness temperatures that were too cold in most cases. On the other hand, the fluffy snow particle parameterizations with and without a melting layer produced brightness temperatures that were much warmer than the observations at these high frequencies. All low frequency calculations track the observations well except in the eye of Hurricane Erin where the surface wind speed and/or near surface temperatures were not properly adjusted to warm the brightness temperature calculations to values closer to the observations. Overall, this work shows that selecting a proper ice particle parameterization can be difficult and additional in situ measurements are needed to define appropriate parameterizations.

Skofronick-Jackson, G.↗

The SPECTRAL Ice Chamber: Application to Titan's Stratospheric Ice Clouds

An ice chamber for obtaining thin-film transmission spectra of Titan-relevant organic ices between 50 and 11,700 cm−1 (200–0.85 μm) is described in detail in this work. The ice chamber, called the SPECtroscopy of Titan-Related ice AnaLogs chamber, is located in the Spectroscopy for Planetary ICes Environments laboratory at NASA Goddard Space Flight Center. Organic vapors are typically deposited between 30 and 150 K via a variable leak valve onto a chemical vapor deposition diamond substrate, where condensation directly into the solid state occurs. Quantitative thin-film infrared transmission spectra of the given ice or ice mixture are then measured, from which optical constants (complex indices of refraction) are computed. These optical constants in turn provide the necessary input for interpreting radiative transfer analyses of Cassini’s Composite InfraRed Spectrometer observations of Titan’s organic stratospheric ice clouds to determine their chemical compositions, abundances, and vertical distributions as functions of Titan latitude and season.

Anderson, C. M.↗

The SPECTRAL Ice Chamber: Application to Titan's Stratospheric Ice Clouds

An ice chamber for obtaining thin-film transmission spectra of Titan-relevant organic ices between 50 and 11,700 cm(exp −1) (200–0.85 μm) is described in detail in this work. The ice chamber, called the SPECtroscopy of Titan-Related ice AnaLogs chamber, is located in the Spectroscopy for Planetary ICes Environments laboratory at NASA Goddard Space Flight Center. Organic vapors are typically deposited between 30 and 150 K via a variable leak valve onto a chemical vapor deposition diamond substrate, where condensation directly into the solid state occurs. Quantitative thin-film infrared transmission spectra of the given ice or ice mixture are then measured, from which optical constants (complex indices of refraction) are computed. These optical constants in turn provide the necessary input for interpreting radiative transfer analyses of Cassini’s Composite InfraRed Spectrometer observations of Titan’s organic stratospheric ice clouds to determine their chemical compositions, abundances, and vertical distributions as functions of Titan latitude and season.

C. M. Anderson↗

Capturing the interactions between ice sheets, sea level and the solid Earth on a range of timescales: a new “time window” algorithm

Retreat and advance of ice sheets perturb the gravitational field, solid surface and rotation of the Earth, leading to spatially variable sea-level changes over a range of timescales $\textit{O}$ (10 0–6 years), which in turn feed back onto ice-sheet dynamics. Coupled ice-sheet–sea-level models have been developed to capture the interactive processes between ice sheets, sea level and the solid Earth, but it is computationally challenging to capture short-term interactions $\textit{O}$ (10 0–2 years) precisely within longer $\textit{O}$ (10 3–6 years) simulations. The standard forward sea-level modelling algorithm assigns a uniform temporal resolution in the sea-level model, causing a quadratic increase in total CPU time with the total number of input ice history steps, which increases with either the length or temporal resolution of the simulation. In this study, we introduce a new “time window” algorithm for 1D pseudo-spectral sea-level models based on the normal mode method that enables users to define the temporal resolution at which the ice loading history is captured during different time intervals before the current simulation time. Utilizing the time window, we assign a fine temporal resolution $\textit{O}$ (10 0–2 years) for the period of ongoing and recent history of surface ice and ocean loading changes and a coarser temporal resolution $\textit{O}$ (10 3–6 years) for earlier periods in the simulation. This reduces the total CPU time and memory required per model time step while maintaining the precision of the model results. We explore the sensitivity of sea-level model results to the model temporal resolution and show how this sensitivity feeds back onto ice-sheet dynamics in coupled modelling. We apply the new algorithm to simulate sea-level changes in response to global ice-sheet evolution over two glacial cycles and the rapid collapse of marine sectors of the West Antarctic Ice Sheet in the coming centuries and provide appropriate time window profiles for each application. The time window algorithm reduces the total CPU time by ~ 50 % in each of these examples and changes the trend of the total CPU time increase from quadratic to linear. This improvement would increase with longer simulations than those considered here. Our algorithm also allows for coupling time intervals of annual temporal scale for coupled ice-sheet–sea-level modelling of regions such as West Antarctica that are characterized by rapid solid Earth response to ice changes due to the thin lithosphere and low mantle viscosities.

59 BASIC BIOLOGICAL SCIENCES↗

Percentage Contributions from Atmospheric and Surface Features to Computed Brightness Temperatures

Over the past few years, there has become an increasing interest in the use of millimeter-wave (mm-wave) and sub-millimeter-wave (submm-wave) radiometer observations to investigate the properties of ice particles in clouds. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. One methodology used since the 1960's to discern the relationship between the physical state observed and the brightness temperature (TB) is through the temperature weighting function profile. In this research, the temperature weighting function concept is exploited to analyze the sensitivity of various characteristics of the cloud profile, such as relative humidity, ice water path, liquid water path, and surface emissivity. In our numerical analysis, we compute the contribution (in Kelvin) from each of these cloud and surface characteristics, so that the sum of these various parts equals the computed TB. Furthermore, the percentage contribution from each of these characteristics is assessed. There is some intermingling/contamination of the contributions from various components due to the integrated nature of passive observations and the absorption and scattering between the vertical layers, but all in all the knowledge gained is useful. This investigation probes the sensitivity over several cloud classifications, such as cirrus, blizzards, light snow, anvil clouds, and heavy rain. The focus is on mm-wave and submm-wave frequencies, however discussions of the effects of cloud variations to frequencies as low as 10 GHz and up to 874 GHz will also be presented. The results show that nearly 60% of the TB value at 89 GHz comes from the earth's surface for even the heaviest blizzard snow rates. On the other hand, a significant percentage of the TB value comes from the snow in the cloud for 166, and 183 plus or minus 7 GHz for the heavy and medium snow rates. For submm-wave channels, there is no contribution from the surface because these channels cannot probe through clouds, nor normal water vapor amounts in clear air regions. This work is extremely valuable in physically-based retrieval algorithm development research.

Gail Skofronick-Jackson↗

Quantitative and interpretable order parameters for phase transitions from persistent homology

Here, we apply modern methods in computational topology to the task of discovering and characterizing phase transitions. As illustrations, we apply our method to four two-dimensional lattice spin models: the Ising, square ice, XY, and fully frustrated XY models. In particular, we use persistent homology, which computes the births and deaths of individual topological features as a coarse-graining scale or sublevel threshold is increased, to summarize multiscale and high-point correlations in a spin configuration. We employ vector representations of this information called persistence images to formulate and perform the statistical task of distinguishing phases. For the models we consider, a simple logistic regression on these images is sufficient to identify the phase transition. Interpretable order parameters are then read from the weights of the regression. This method suffices to identify magnetization, frustration, and vortex-antivortex structure as relevant features for phase transitions in our models. We also define “persistence” critical exponents and study how they are related to those critical exponents usually considered.

36 MATERIALS SCIENCE↗

The SPECTRAL Ice Chamber: Application to Titan's Stratospheric Ice Clouds

An ice chamber for obtaining thin-film transmission spectra of Titan-relevant organic ices between 50 and 11,700 cm−1 (200 to 0.85 µm) is described in detail in this work. The ice chamber, called the SPECtroscopy of Titan-Related ice AnaLogs (SPECTRAL) chamber, is located in the Spectroscopy for Planetary ICes Environments (SPICE) laboratory at NASA Goddard Space Flight Center (GSFC). Organic vapors are typically deposited between 30 K and 150 K via a variable leak valve onto a chemical vapor deposition (CVD) diamond substrate, where condensation directly into the solid state occurs. Quantitative thin-film infrared transmission spectra of the given ice or ice mixture are then measured, from which optical constants (complex indices of refraction) are computed. These optical constants in turn provide the necessary input for interpreting radiative transfer analyses of Cassini’s Composite InfraRed Spectrometer (CIRS) observations of Titan’s organic stratospheric ice clouds to determine their chemical compositions, abundances, and vertical distributions as functions of Titan latitude and season.

C. M. Anderson↗

Impact Ice Microstructure Segmentation Using Transfer Learned Model

A process of using machine learning to segment impact ice microstructure is presented and analyzed. The segmentation was conducted with the goal of obtaining average grain size estimations. The model was trained on a set of micrographs of impact ice grown at NASA Glenn’s Icing Research Tunnel. The model leveraged a model pre-trained on a large set of micrographs of various materials as a starting point. Post-processing of the segmented images was done to connect broken boundaries. An automatic method of determining grain size following an ASTM standard was implemented. Segmentation results using different training sets as well as different encoder and decoder pairs are presented. Calculated sizes are compared to manual grain size measurement methods. Results show promise in accuracy as well as a possible improvement in repeatability and consistency. Next steps for improving the model are suggested.

Machine learning↗

Debris/Ice/TPS Assessment and Photographic Analysis for Shuttle Mission STS-40

A debris, ice, Thermal Protection System (TPS) assessment and photographic analysis for Space Shuttle Mission STS-40 was conducted. Debris inspections of the flight elements and launch pad were performed before and after launch. Ice and frost conditions on the External Tank were assessed by the use of computer programs, nomographs, and infrared scanner data during cryogenic loading of the vehicle, followed by on-pad visual inspection. High speed photography was analyzed after launch to identify ice and debris sources and to evaluate potential vehicle damage and/or in-flight anomalies.

Katnik, Gregory N.↗

A double-moment multiple-phase four-class bulk ice scheme, Part 1: Description

A detailed ice-phase bulk microphysical scheme has been developed for simulating the hydrometeor distributions of convective and stratiform precipitation in different large-scale environmental conditions. The proposed scheme involves 90 distinct microphysical processes, which predict the mixing ratios and the number concentrations of small ice crystals, snow, graupel, and frozen drops/hail, as well as the mixing ratios of liquid water on wet precipitation ice (snow, graupel, frozen drops). The number of adjustable coefficients has been significantly reduced in comparison with other bulk schemes. Additional improvements have been made to the parameterization in the following areas: (1) representing small ice crystals with nonzero terminal fall velocities and dispersive size distributions, (2) accurate and computationally efficient calculations of precipitation collection processes, (3) reformulating the collection equation to prevent unrealistically large accretion rates, (4) more realistic conversion by riming between different classes of precipitation ice, (5) preventing unrealistically large rates of raindrop freezing and freezing of liquid water on ice, (6) detailed treatment of various rime-splintering ice multiplication mechanisms, (7) a simple representation of the Hobbs-Rangno ice enhancement process, (8) aggregation of small ice crystals and snow, and (9) allowing explicit competition between cloud water condensation and ice deposition rates rather than using saturation adjustment techniques. For the purposes of conserving the higher moments of the particle distributions, preserving the spectral widths (or slopes) of the particle spectra is shown to be more important than strict conservation of particle number concentration when parameterizing changes in ice-particle number concentrations due to melting, vapor transfer processes (sublimation of dry ice, evaporation from wet ice), and conversion between different hydrometeor species. The microphysical scheme is incorporated into a nonhydrostatic cloud model in Part 2 of this study. The model performed well in simulating the radar and microphysical structures of a midlatitude-continental squall lines and a tropical-maritime squall system with minimal tuning of the parameterization, even though the vertical profiles of radar reflectivity differed substantially between these storms.

Ferrier, Brad Shoenberg↗

General equations for the motions of ice crystals and water drops in gravitational and electric fields

General equations for the Reynolds number of a variety of types of ice crystals and water drops are given in terms of the Davies, Bond, and Knudsen numbers. The equations are in terms of the basic physical parameters of the system and are valid for calculating velocities in gravitational and electric fields over a very wide range of sizes and atmospheric conditions. The equations are asymptotically matched at the bottom and top of the size spectrum, useful when checking large computer codes. A numerical system for specifying the dimensional properties of ice crystals is introduced. Within the limits imposed by such variables as particle density, which have large deviations, the accuracy of velocities appears to be within 10 percent over the entire range of sizes of interest.

Nisbet, John S.↗

General equations for the motions of ice crystals and water drops in gravitational and electric fields

General equations for the Reynolds number of a variety of types of ice crystals and water drops are given in terms of the Davies, Bond, and Knudsen numbers. The equations are in terms of the basic physical parameters of the system and are valid for calculating velocities in gravitational and electric fields over a very wide range of sizes and atmospheric conditions. The equations are asymptotically matched at the bottom and top of the size spectrum, useful when checking large computer codes. A numerical system for specifying the dimensional properties of ice crystals is introduced. Within the limits imposed by such variables as particle density, which have large deviations, the accuracy of velocities appears to be within 10 percent over the entire range of sizes of interest.

Nisbet, John S.↗