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 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

NASA NTRS2017

A Fast Hybrid (3-D/1-D) Model for Thermal Radiative Transfer in Cirrus via Successive Orders of Scattering

We investigate the impact of cirrus cloud heterogeneity on the direct emission by cloud or surface and on the scattering by ice particles in the thermal infrared (TIR). Realistic 3-D cirri are modeled with the 3DCLOUD code, and top-of-atmosphere radiances are simulated by the 3-D Monte Carlo radiative transfer (RT) algorithm 3DMCPOL for two (8.65 micrometers and 12.05 micrometers) channels of the Imaging Infrared Radiometer on CALIPSO. At nadir, comparisons of 1-D and 3-D RT show that 3-D radiances are larger than their 1-D counterparts for direct emission but smaller for scattered radiation. For our cirrus cases, 99% of the 3-D total radiance is computed by the third scattering order, which corresponds to 90% of the total computational effort, but larger optical thicknesses need more scattering orders. To radically accelerate the 3-D RT computations (using only few percent of 3-D RT time with a Monte Carlo code), even in the presence of large optical depths, we develop a hybrid model based on exact 3-D direct emission, the first scattering order from 1-D in each homogenized column, and an empirical adjustment linearly dependent on the optical thickness to account for higher scattering orders. Good agreement is found between the hybrid model and the exact 3-D radiances for two very different cirrus models without changing the empirical parameters. We anticipate that a future deterministic implementation of the hybrid model will be fast enough to process multiangle thermal imagery in a practical tomographic reconstruction of 3-D cirrus fields.

radiative transfer↗

DOE OSTI2025

Ice storage model-predictive control in an office building with PV: scenario, error and sensitivity analysis

Thermal energy storage (TES) can enable more building-sited renewable electricity generation and lower utility bill costs for buildings owners and occupants, especially when there are high demand and variable time-of-use (TOU) charges. A model predictive control (MPC) strategy can offer additional savings over a schedule-based control with added complexity and reliance on forecasts. Here, this study examines savings for medium office buildings with chiller plants in three locations with building-installed solar photovoltaics (PV) to understand the impact of MPC. Control setpoints are fixed by a schedule-based control or optimized by nonlinear MPC. These control setpoints are actuated within EnergyPlus building models to simulate the utility cost of the chiller plant. NLP solutions can be unstable or unrealistic, but our results show that by regularizing the NLP, the solutions can be reasonably followed by the building model. MPC models make simplifications that lead to errors once the controller is participating in and changing the operation of the building. These errors average 9 % across the cases, showing that the most important parts of the system are represented. The no-thermal load costs are computed to show that the optimization can in some cases achieve both the minimum TOU and minimum monthly demand costs by demand management while reducing TOU energy costs by energy arbitrage. The MPC saves 35–66 % in the annual chiller plant operating costs, which is an additional savings above the schedule by 1–33 %. PV and TES are complementary and mostly independent, but a load with PV often results in better performance for the schedule. Our case study and sensitivity analysis show the importance of modeling and optimization for complex rates, but also the circumstances wherein a simpler strategy achieves the same performance with less potential for error.

14 SOLAR ENERGY↗

NASA NTRS1982

The interaction of radio frequency electromagnetic fields with atmospheric water droplets and applications to aircraft ice prevention

The feasibility of computerized simulation of the physics of advanced microwave anti-icing systems, which preheat impinging supercooled water droplets prior to impact, was investigated. Theoretical and experimental work performed to create a physically realistic simulation is described. The behavior of the absorption cross section for melting ice particles was measured by a resonant cavity technique and found to agree with theoretical predictions. Values of the dielectric parameters of supercooled water were measured by a similar technique at lambda = 2.82 cm down to -17 C. The hydrodynamic behavior of accelerated water droplets was studied photograhically in a wind tunnel. Droplets were found to initially deform as oblate spheroids and to eventually become unstable and break up in Bessel function modes for large values of acceleration or droplet size. This confirms the theory as to the maximum stable droplet size in the atmosphere. A computer code which predicts droplet trajectories in an arbitrary flow field was written and confirmed experimentally. The results were consolidated into a simulation to study the heating by electromagnetic fields of droplets impinging onto an object such as an airfoil. It was determined that there is sufficient time to heat droplets prior to impact for typical parameter values. Design curves for such a system are presented.

Hansman, R. J., Jr.↗

NASA NTRS1985

Measurements in liquid fuel sprays

A ground test facility is being established at NASA Lewis Research Center to simulate the environmental and flight conditions needed to study adverse weather effects. One of the most important components is the water spray system which consists of many nozzles fitted on spray bars. Water is injected through air-assisted atomizers to generate uniform size drops to simulate icing in clouds. The primary objective is to provide experimental data on drop size distribution over a wide range of operating conditions. Correlation equations for mean drop size and initial injection parameters are being determined to assist in the design and modification of the Altitude Wind Tunnel. Special emphasis is being placed on the study of the aerodynamic structure of the air-assisted atomizer sprays. Detailed measurements of the variation of drop size distribution and velocity as a function of time and space are being made. Accurate initial and boundary conditions are being provided for computer model evaluation.

Chigier, N.↗

NASA NTRS1994

The global distribution of ozone destruction rates obtained from 13 years of Nimbus/TOMS data (1979-1991)

Long-term ozone trends (percentage change) have been computed from 13 years of Nimbus/TOMS (Total Ozone Mapping Spectrometer) data as a function of latitude, longitude, and month for the period January 1, 1979 to December 31, 1991. In both hemispheres, the ozone column content has decreased at latitudes above 30 deg by amounts that are larger than predicted by homogeneous chemistry models for the 13-year time period. The largest rates of ozone decrease occur in the Southern Hemisphere during winter and spring, with recovery during the summer and autumn. The large winter ozone loss rates are consistent with observed low stratospheric temperatures, ice-cloud formation, and heterogeneous chemistry at middle and high latitudes. There are similar, but smaller changes observed in the Northern Hemisphere. At midlatitudes, (40 deg N to 50 deg N) there are increased zonal average ozone depletion rates that correspond to 5 geographically localized regions of increased ozone depletion rates. Only the equatorial band between plus or minus 20 deg shows little or no long-term ozone change since January, 1979. The long-term winter ozone depletion rate data for both hemispheres suggests that heterogeneous chemistry processes may operate over a wide range of latitudes during half of the year.

Herman, Jay R.↗

NASA NTRS2003

Radar Reflectivity Simulated by a 2-D Spectra Bin Model: Sensitivity of Cloud-aerosol Interaction

The Goddard Cumulus Ensemble (GCE) model with bin spectra microphysics is used to simulate mesoscale convective systems.The model uses explicit bins to represent size spectra of cloud nuclei, water drops, ice crystals, snow and graupel. Each hydrometeorite category is described by 33 mass bins. The simulations provide a unique data set of simulated raindrop size distribution in a realistic dynamic frame. Calculations of radar parameters using simulated drop size distribution serve as an evaluation of numerical model performance. In addition, the GCE bin spectra modes is a very useful tool to study uncertainties related to radar observations; all the environmental parameters are precisely known. In this presentation, we concentrate on the discussion of Z-R (ZDR-R) relation in the simulated systems. Due to computational limitations, the spectra bin model has been run in two dimensions with 31 stretched vertical layers and 1026 horizontal grid points (1 km resolution). Two different cases, one in midlatitude continent, the other in tropical ocean, have been simulated. The continental case is a strong convection which lasted for two hours. The oceanic case is a persistent system with more than 10 hours' life span. It is shown that the simulated Z-R (ZDR-R) relations generally agree with observations using radar and rain gauge data. The spatial and temporal variations of Z-R relation in different locations are also analyzed. Impact of aerosols on cloud formation and raindrop size distribution was studied. Both clean (low CCN) and dirty (high CCN) cases are simulated. The Z-R relation is shown to vary considerable in the initial CCN concentrations.

Li, Kiaowen↗

NASA NTRSDate not supplied

ICME for NASA Aerospace Applications: Batteries for Electric Aviation

NASA’s approach to computational materials modeling is detailed in the NASA Vision 2040 Roadmap for Multiscale Modeling and Simulation of Materials and Systems. This report is in the spirit of national initiatives such as the Material Genome Initiative (MGI), Integrated Computational Materials Engineering (ICME), and others. We utilize a combination of fundamental modeling, computational high-throughput screening, and data science methods, e.g., machine learning, are used to find innovative solutions to NASA or national technology challenges. Applications of interest are wide ranging from advanced alloys to batteries to coatings, among others. In this talk, we present three examples for recent work related to NASA applications. First, doping advanced sulfur battery cathodes with selenium boosts electrical conductivity important for electric aircraft applications. First principles calculations will be discussed that result in compositional design maps for these materials. Second, development of icephobic coatings is important to mitigate safety hazards associated with icing for aircraft. Molecular dynamics simulations are reported for ice-surface interfaces to understand adhesion mechanisms and help screen optimal ice-phobic coatings. Third, shape memory alloys have numerous applications as actuators, superelastic materials, etc. for aerospace. We report machine learning models that predict martensitic transition temperatures across a broad swath of compositional space.

John Lawson↗

NASA NTRS2017

Blowing Snow Sublimation and Transport over Antarctica from 11 Years of CALIPSO Observations

Blowing snow processes commonly occur over the earth's ice sheets when the 10 mile wind speed exceeds a threshold value. These processes play a key role in the sublimation and redistribution of snow thereby influencing the surface mass balance. Prior field studies and modeling results have shown the importance of blowing snow sublimation and transport on the surface mass budget and hydrological cycle of high-latitude regions. For the first time, we present continent-wide estimates of blowing snow sublimation and transport over Antarctica for the period 2006-2016 based on direct observation of blowing snow events. We use an improved version of the blowing snow detection algorithm developed for previous work that uses atmospheric backscatter measurements obtained from the CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) lidar aboard the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) satellite. The blowing snow events identified by CALIPSO and meteorological fields from MERRA-2 are used to compute the blowing snow sublimation and transport rates. Our results show that maximum sublimation occurs along and slightly inland of the coastline. This is contrary to the observed maximum blowing snow frequency which occurs over the interior. The associated temperature and moisture reanalysis fields likely contribute to the spatial distribution of the maximum sublimation values. However, the spatial pattern of the sublimation rate over Antarctica is consistent with modeling studies and precipitation estimates. Overall, our results show that the 2006-2016 Antarctica average integrated blowing snow sublimation is about 393 +/- 196 Gt yr(exp -1), which is considerably larger than previous model-derived estimates. We find maximum blowing snow transport amount of 5 Mt km-1 yr(exp -1) over parts of East Antarctica and estimate that the average snow transport from continent to ocean is about 3.7 Gt yr(exp -1). These continent-wide estimates are the first of their kind and can be used to help model and constrain the surface mass budget over Antarctica.

Antarctica Sublimation Ice-sheet↗

DOE OSTI2026

Time irreversibility as an indicator of approaching tipping points in Earth subsystems

With shifting environmental trends, many Earth system elements may be poised to undergo critical transitions or ‘tipping’. Reliable anticipation of these tipping elements is vital to inform policy decisions. Many of the current methods for tipping point detection are based on loss of resilience or ‘critical slowdown’ of the system as it approaches a tipping point. However, these methods are prone to false alarms; the detected slowdown may be an artifact of nonstationary noise unrelated to tipping behavior. Here, we explore the efficacy of early warning signs based on a nonequilibrium thermodynamics framework. The model-free detection method relies on the increased intrinsic time-irreversibility due to detailed balance breaking, preceding the onset of tipping or instabilities. We demonstrate that these EWSs are effective for tipping point detection and robust against false alarms due to nonstationary noise, using idealized models for two key elements of the Earth system that are prone to tipping: the Atlantic Meridional Overturning Circulation and Arctic sea-ice loss.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Extensively Microtwinned Diamond with Nanolaminates of Lonsdaleite Formed by Flash Laser Heating of Glassy Carbon

Diamond’s unique properties on the nanoscale make it one of the most important materials for use in biosensors and quantum computing and for components that can withstand the harsh environments of space. We synthesize oriented, faceted diamond particles by flash laser heating of glassy carbon at 16 GPa and 2300 K. Detailed transmission electron microscopy shows them to consist of a mosaic of diamond nanocrystals frequently joined at twin boundaries forming microtwins. Striking 3-fold translational periodicity was observed in both imaging and diffraction. This periodicity was shown to originate from nanodimensional wedge-shaped overlapping regions of twinned diamond and not from a possible 9R polytype, which has also been reported in other group IVa elements and water ice. Extended bilayers of hexagonal layer stacking were observed, forming lonsdaleite nanolaminates. In conclusion, the particles exhibited optical fluorescence with a rapid quench time (<1 ns) attributed to their unique twinned microstructure.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

NASA NTRS1992

The Global Ocean Observing System

A Global Ocean Observing System (GOOS) should be established now with international coordination (1) to address issues of global change, (2) to implement operational ENSO forecasts, (3) to provide the data required to apply global ocean circulation models, and (4) to extract the greatest value from the one billion dollar investment over the next ten years in ocean remote sensing by the world's space agencies. The objectives of GOOS will focus on climatic and oceanic predictions, on assessing coastal pollution, and in determining the sustainability of living marine resources and ecosystems. GOOS will be a complete system including satellite observations, in situ observations, numerical modeling of ocean processes, and data exchange and management. A series of practical and economic benefits will be derived from the information generated by GOOS. In addition to the marine science community, these benefits will be realized by the energy industries of the world, and by the world's fisheries. The basic oceanic variables that are required to meet the oceanic and predictability objectives of GOOS include wind velocity over the ocean, sea surface temperature and salinity, oceanic profiles of temperature and salinity, surface current, sea level, the extent and thickness of sea ice, the partial pressure of CO2 in surface waters, and the chlorophyll concentration of surface waters. Ocean circulation models and coupled ocean-atmosphere models can be used to evaluate observing system design, to assimilate diverse data sets from in situ and remotely sensed observations, and ultimately to predict future states of the system. The volume of ocean data will increase enormously over the next decade as new satellite systems are launched and as complementary in situ measuring systems are deployed. These data must be transmitted, quality controlled, exchanged, analyzed, and archived with the best state-of-the-art computational methods.

Kester, Dana↗

NASA NTRSDate not supplied

Evaluating Lunar Water Processing System Model Configurations for Small Scale Oxygen and Hydrogen Production Within JAXA’s ISRU Technology Demonstration Plant

In-Situ Resource Utilization (ISRU) refers to novel methods of extracting and processing local resources for use in life support and propulsion systems, reducing or eliminating the required consumables to be transferred from Earth. Current estimates of water-ice availability embedded in regolith within the Moon’s permanently shadowed regions (PSR’s)range between 1-5% by weight. However, the composition and characteristics of the “wet” regolith is unknown. Alternate ISRU excavation techniques and Concept of Operations (ConOps) must be explored to optimize surface system operations based on these factors. To assess the feasibility of different ISRU subsystem technologies and compare system architecture configurations, an interchangeable system model was generated to incorporate technologies spanning excavation of raw materials to storage of products and determine optimal arrangement of total system processing needs. Total Mass, Volume, and Power (M/V/P) requirements were computed for 168 design iterations of this water processing plant.

ISRU↗

DOE OSTI2025

Evaluating the Trustworthiness of Explainable Artificial Intelligence (XAI) Methods Applied to Regression Predictions of Arctic Sea Ice Motion

Abstract Recent advances in explainable artificial intelligence (XAI) methods show promise for understanding predictions made by machine learning (ML) models. XAI explains how the input features are relevant or important for the model predictions. We train linear regression (LR) and convolutional neural network (CNN) models to make 1-day predictions of sea ice velocity in the Arctic from inputs of present-day wind velocity and previous-day ice velocity and concentration. We apply XAI methods to the CNN and compare explanations to variance explained by LR. We confirm the feasibility of using a novel XAI method [i.e., global layerwise relevance propagation (LRP)] to understand ML model predictions of sea ice motion by comparing it to established techniques. We investigate a suite of linear, perturbation-based, and propagation-based XAI methods in both local and global forms. Outputs from different explainability methods are generally consistent in showing that wind speed is the input feature with the highest contribution to ML predictions of ice motion, and we discuss inconsistencies in the spatial variability of the explanations. Additionally, we show that the CNN relies on both linear and nonlinear relationships between the inputs and uses nonlocal information to make predictions. LRP shows that wind speed over land is highly relevant for predicting ice motion offshore. This provides a framework to show how knowledge of environmental variables (i.e., wind) on land could be useful for predicting other properties (i.e., sea ice velocity) elsewhere. Significance Statement Explainable artificial intelligence (XAI) is useful for understanding predictions made by machine learning models. Our research establishes trustability in a novel implementation of an explainable AI method known as layerwise relevance propagation for Earth science applications. To do this, we provide a comparative evaluation of a suite of explainable AI methods applied to machine learning models that make 1-day predictions of Arctic sea ice velocity. We use explainable AI outputs to understand how the input features are used by the machine learning to predict ice motion. Additionally, we show that a convolutional neural network uses nonlinear and nonlocal information in making its predictions. We take advantage of the nonlocality to investigate the extent to which knowledge of wind on land is useful for predicting sea ice velocity elsewhere.

Hoffman, Lauren [Scripps Institution of Oceanograp↗

NASA NTRS1975

The effects of localized damping on structural response

The effect of localized structural damping on the excitability of higher order normal modes of the large space telescope was investigated. A preprocessor computer program was developed to incorporate Voigt structural joint damping models in a NASTRAN finite-element dynamic model. A postprocessor computer program was developed to select critical modes for low-frequency attitude control problems and for higher frequency fine-stabilization problems. The mode selection is accomplished by ranking the flexible modes based on coefficients for rate gyro, position gyro, and optical sensors, and on image-plane motions due to sinusoidal or random power spectral density force and torque inputs.

Merchant, D. H.↗

NASA NTRS1990

Stratospheric constituent trends from ER-2 profile data

Trend analysis of the Airborne Arctic Stratospheric Expedition's ER-2 profile data reveals an average decrease in N2O on potential temperature isentropes which can be attributed to diabatic cooling of inner vortex air. This conclusion is independently supported by radiative transfer computations. Trends in ozone and water vapor over the same period are not consistent with the magnitude of the diabatic descent. After accounting for the diabatic motion (estimated from N2O), an additional 0.44 + or - 0.3 percent/day average anomalous O3 decrease above 440 K (about 20 km) is needed to balance the continuity equation. This ozone decrease suggests additional photochemical destruction of ozone in the presence of the high amounts of ClO observed during the mission. A 0.4 + or - 0.3 percent/day average anomalous increase in H2O is also observed near 420 K (about 18 km) which may be due to the evaporation of ice crystals falling from higher, colder stratospheric layers.

Schoeberl, Mark R.↗

DOE OSTI2023

Multirate partitioned Runge–Kutta methods for coupled Navier–Stokes equations

Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. Further, this study explores multirate partitioned Runge-Kutta methods for the fluid-fluid interaction problem and demonstrates its parallel performance by using the PETSc library. We consider compressible Navier-Stokes equations with gravity coupled through a rigid-lid interface. Our large-scale numerical experiments reveal that multirate partitioned Runge-Kutta coupling schemes (1) can conserve total mass; (2) have second-order accuracy in time; and (3) provide favorable strong- and weak-scaling performance on modern computing architectures. We also show that the speedup factors of multirate partitioned Runge-Kutta methods match theoretical expectations over their base (single-rate) method.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Intercomparison of Thermal Regime Algorithms in 1-D Lake Models

Lakes are an important component of the global weather and climate system, but the modeling of their thermal regimes has shown large uncertainties due to the highly diverse lake properties and model configurations. Here we evaluate the algorithms of four key lake thermal processes including turbulent heat fluxes, wind-driven mixing, light extinction, and snow density, using a highly diverse lake dataset provided by the Inter-sectoral Impact Model Intercomparison Project (ISIMIP) 2a lake sector. Algorithm codes are configured and run separately within the same parent model to rule out any interference from factors apart from the algorithms examined. Evaluations are based on both simulation accuracy and recalibration complexity for application to global lakes. For turbulent heat fluxes, the non-Monin-Obukhov similarity (MOS) based, more simplified algorithms perform better in predicting lake epilimnion temperatures and achieve high convergence in the values of the calibrated parameters. For wind-driven mixing, a two-algorithm strategy considering lake shape and season is suggested with the regular mixing algorithm used for spring and earlier summer and the mixing-enhanced algorithm for summer steady stratification and fall overturn periods. There are no evident differences in the simulated thermocline depths using different light extinction algorithms or the observation. Finally, for lake ice phenology, a constant snow density at around 110 kg m-3 is found to be sufficient for most northern lakes while the Arctic lakes require a higher value. Our study provides highly practical guides for improving 1-D lake models and feasible parameterization strategies to better simulate global lake thermal regimes.

58 GEOSCIENCES↗

NASA NTRS2000

A Melting Layer Model for Passive/Active Microwave Remote Sensing Applications: Simulation of TRMM Observations - Part 2

The one-dimensional, steady-state melting layer model developed in Part I of this study is used to calculate both the microphysical and radiative properties of melting precipitation, based upon the computed concentrations of snow and graupel just above the freezing level at applicable horizontal gridpoints of 3-dimensional cloud resolving model simulations. The modified 3-dimensional distributions of precipitation properties serve as input to radiative transfer calculations of upwelling radiances and radar extinction/reflectivities at the TRMM Microwave Imager (TMI) and Precipitation Radar (PR) frequencies, respectively. At the resolution of the cloud resolving model grids (approx. 1 km), upwelling radiances generally increase if mixed-phase precipitation is included in the model atmosphere. The magnitude of the increase depends upon the optical thickness of the cloud and precipitation, as well as the scattering characteristics of ice-phase precipitation aloft. Over the set of cloud resolving model simulations utilized in this study, maximum radiance increases of 43, 28, 18, and 10 K are simulated at 10.65, 19.35 GHz, 37.0, and 85.5 GHz, respectively. The impact of melting on TMI-measured radiances is determined not only by the physics of the melting particles but also by the horizontal extent of the melting precipitation, since the lower-frequency channels have footprints that extend over 10''s of kilometers. At TMI resolution, the maximum radiance increases are 16, 15, 12, and 9 K at the same frequencies. Simulated PR extinction and reflectivities in the melting layer can increase dramatically if mixed-phase precipitation is included, a result consistent with previous studies. Maximum increases of 0.46 (-2 dB) in extinction optical depth and 5 dBZ in reflectivity are simulated based upon the set of cloud resolving model simulations.

Olson, William S.↗