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 361 records · Page 20

Numerical simulation of an electrothermal deicer pad

A numerical simulation is developed to investigate the removal of ice from composite aircraft blades by means of electrothermal deicing. The model considers one dimensional, unsteady state heat transfer in the composite blade-ice body. The heat conduction equations are approximated by using the Crank-Nicolson finite difference scheme, and the phase change in the ice layer is handled using the Enthalpy method. To solve the system of equations which result, Gauss-Seidel iteration is used. The simulation computes the temperature profile in the composite blade-ice body, as well as the movement of the ice-water interface, as a function of time. This information can be used to evaluate deicer performance. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.

Marano, J. J.↗

Evaluation of capillary reinforced composites for anti-icing

This paper discusses the evaluation of glass capillary reinforced advanced composite structures for anti-icing purposes. The concept involves embedding glass capillary tubes on the surface of a composite structure and ducting heated air through the tubes. A computer program was developed to predict the anti-icing performance of such tubes and a test program was conducted to demonstrate the actual performance of this system. Test data and analytical code results were in excellent agreement. Both indicate the feasibility of using capillary tubes for surface heating in order to combat ice accumulation on advanced composite structures.

Ciardullo, Samuel W.↗

Should Sea-Ice Modeling Tools Designed for Climate Research Be Used for Short-Term Forecasting?

In theory, the same sea-ice models could be used for both research and operations, but in practice, differences in scientific and software requirements and computational and human resources complicate the matter. Although sea-ice modeling tools developed for climate studies and other research applications produce output of interest to operational forecast users, such as ice motion, convergence, and internal ice pressure, the relevant spatial and temporal scales may not be sufficiently resolved. For instance, sea-ice research codes are typically run with horizontal resolution of more than 3 km, while mariners need information on scales less than 300 m. Certain sea-ice processes and coupled feedbacks that are critical to simulating the Earth system may not be relevant on these scales; and therefore, the most important model upgrades for improving sea-ice predictions might be made in the atmosphere and ocean components of coupled models or in their coupling mechanisms, rather than in the sea-ice model itself. This paper discusses some of the challenges in applying sea-ice modeling tools developed for research purposes for operational forecasting on short time scales, and highlights promising new directions in sea-ice modeling.

54 ENVIRONMENTAL SCIENCES↗

Synergistic Retrievals of Ice in High Clouds From Elastic Backscatter Lidar, Ku-band Radar and Submillimeter Wave Radiometer Observations

In this study, we investigate the synergy of elastic backscatter lidar, Ku-band radar, and sub-millimeter-wave radiometer measurements in the retrieval of ice from satellite observations. The synergy is analyzed through the generation of a large dataset of IceWater Content (IWC) profiles and simulated lidar, radar and radiometer observations. The characteristics of the instruments e.g. frequencies, sensitivities, etc. are set based on the expected characteristics of instruments of the Atmosphere Observing System (AOS) mission. A hold-out validation methodology is used to assess the accuracy of the IWC profiles retrieved from various combinations of observations from the three instruments. Specifically, the IWC and associated observations are randomly divided into two datasets, one for training and the other for evaluation. The training dataset is used to train the retrieval algorithm, while the evaluation dataset is used to assess the retrieval performance. The dataset of IWC profiles is derived from CloudSat reflectivity and CALIOP lidar observations. The retrieval of the ice water content IWC profiles from the computed observations is achieved in two steps. In the first step, a class, out of 18 potential classes characterized by different vertical distribution of IWC, is estimated from the observations. The 18 classes are predetermined based on the k-Means clustering algorithm. In the second step, the IWC profile is estimated using an Ensemble Kalman Smoother (EKS) algorithm that uses the estimated class as a priori information. The results of the study show that the synergy of lidar, radar, and radiometer observations is significant in the retrieval of the IWC profiles. Nevertheless, it should be mentioned that this synergy was found under idealized conditions, and additional work might be required to materialize it in practice. The inclusion of the lidar backscatter observations in the retrieval process has a larger impact on the retrieval performance than the inclusion of the radar observations. As ice clouds have a significant impact on atmospheric radiative processes, this work is relevant to ongoing efforts to reduce uncertainties in climate analyses and projections.

Mircea Grecu↗

A Model to Assess the Risk of Ice Accretion Due to Ice Crystal Ingestion in a Turbofan Engine and its Effects on Performance

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 were attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was one or more of the following anomalies: degraded engine performance, engine roll back, compressor surge and stall, and flameout of the combustor. The main focus of this research is the development of a computational tool that can estimate whether there is a risk of ice accretion by tracking key parameters through the compression system blade rows at all engine operating points within the flight trajectory. The tool has an engine system thermodynamic cycle code, coupled with 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 blade rows. Assumptions are made to predict the complex physics involved in engine icing. Specifically, the code does not directly estimate ice accretion and does not have models for particle breakup or erosion. Two key parameters have been suggested as conditions that must be met at the same location for ice accretion to occur: the local wet-bulb temperature to be near freezing or below and the local melt ratio must be above 10%. These parameters were deduced from analyzing laboratory icing test data and are the criteria used to predict the possibility of ice accretion within an engine including the specific blade row where it could occur. Once the possibility of accretion is determined from these parameters, the degree of blockage due to ice accretion on the local stator vane can be estimated from an empirical model of ice growth rate and time spent at that operating point in the flight trajectory. The computational tool can be used to assess specific turbine engines to their susceptibility to ice accretion in an ice crystal environment.

Jorgenson, Philip C. E.↗

An Approach for Simulating Transport of Nonspherical, Sublimating Particles in Rarefied Gas Flows via the Direct Simulation Monte Carlo Method

This paper extends the particle Green’s function approach of Gallis et al. [Phys. Fluids 13, 3482–3492 (2001)] to compute the force and heat transfer on nonspherical, sublimating particles in rarefied gas flows. Green’s functions are derived for rods and plates, which are representative of ice particle shapes, and a numerical procedure is developed for computing the Green’s functions for a general discretized convex shape. Additionally, a sublimation model is developed to simulate the evolution of ice particles in rarefied gas flows. These models are implemented in a direct simulation Monte Carlo code and demonstrated for the case of ice-particle-laden hypersonic flow over a ramp.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Emulator-based Bayesian calibration of a subglacial drainage model

Subglacial drainage models, often motivated by the relationship between hydrology and ice flow, sensitively depend on numerous unconstrained parameters. We explore using borehole water-pressure time series to calibrate the uncertain parameters of a popular subglacial drainage model, taking a Bayesian perspective to quantify the uncertainty in parameter estimates and in the calibrated model predictions. To reduce the computation time associated with Markov Chain Monte Carlo sampling, we construct a fast Gaussian process emulator to stand in for the subglacial drainage model. We first carry out a calibration experiment using synthetic observations consisting of model simulations with hidden parameter values as a demonstration of the method. Using real borehole water pressures measured in western Greenland, we find meaningful constraints on four of the eight model parameters and a factor-of-three reduction in uncertainty of the calibrated model predictions. These experiments illustrate Gaussian process-based Bayesian inference as a useful tool for calibration and uncertainty quantification of complex glaciological models using field data. However, significant differences between the calibrated model and the borehole data suggest that structural limitations of the model, rather than poorly constrained parameters or computational cost, remain the most important constraint on subglacial drainage modelling.

58 GEOSCIENCES↗

Evaluation of capillary reinforced composites

Anti-icing of the inlet of jet engines is generally performed with high pressure heated air that is directed forward from the compressor through a series of pipes to various manifolds located near the structures to be anti-iced. From these manifolds, the air is directed to all flowpath surfaces that may be susceptible to ice formation. There the anti-icing function may be performed by either heat conduction or film heating. Unfortunately, the prospect of utilizing lighweight, high strength composites for inlet structures of jet engines has been frustrated by the low transverse thermal conductivity of such materials. It was the objective of this program to develop an advanced materials and design concept for anti-icing composite structures. The concept that was evaluated used capillary glass tubes embedded on the surface of a composite structure with heated air ducted through the tubes. An analytical computer program was developed to predict the anti-icing performance of such tubes and a test program was conducted to demonstrate actual performance of this system. Test data and analytical code results were in excellent agreement. Both indicate feasibility of using capillary tubes for surface heating as a means for composite engine structures to combat ice accumulation.

Cahill, J. E.↗

The Effect of Surface Chemical Functionality Upon Ice Adhesion

In nature, anti-freeze proteins present in fish utilize specific organic functionalities to disrupt ice crystal formation and propagation. Based on these structures, surfaces with controlled chemical functionality and chain length were evaluated both experimentally and computationally to assess the effect of both parameters in mitigating ice formation. Linear aliphatic dimethylethoxysilanes terminated with methyl or hydroxyl groups were prepared, characterized, and used to coat aluminum. The effect upon icing using a microdroplet freezing apparatus and the Adverse Environment Rotor Test Stand found hydroxyl-terminated materials exhibited a greater propensity for ice formation and adhesion. Molecular dynamics simulations of a silica substrate bearing functionalized species of similar composition were brought into contact with a pre-equilibrated ice crystal. Several parameters including chain mobility were monitored to ascertain the size of a quasi-liquid layer. The studies suggested that chain mobility affected the interface between ice and the surface more than terminal group chemical composition.

Smith, Joseph G., Jr.↗

The Laboratory for Oceans Computing Facility

The first VAX computer in the Laboratory for Oceans Computing Facility (LOCF) was installed and the facility was largely expanded. The growth is not only in hardware and software, but also in the number of users and in supporting research and development projects. The LOCF serves as a general purpose computing facility for: ocean color research projects, sea ice research projects, processing of the Nimbus-7 Coastal Zone Color Scanner data set, real time ingest and analysis of TIROS-N satellite data, study of the Synthetic Aperture Radar data, study of LANDSAT data, and many others. The physical space and the electrical power layout of the computing room were modified to accommodate all the equipment. The LOCF has several image processing stations which include two International Imaging Systems (IIS) model 75 processors and one Adage processor. The facility has the capability of ingesting the TIROS-N HRPT satellite data on a real time basis. More than 30 software packages were installed on the systems. System software packages, network software, FORTRAN and C compilers, database management software, image processing software, graphics, mathematics and statistics packages, TAE, Catalog Manager, GEMPAK, LAS and many other software developed on the LOCF computers such as SEAPAK have greatly advanced the capability of the LOCF.

Kao, R.↗

Atmospheric and cloud structures of the Jovian planets

Wet adiabatic atmospheric models are computed for the four Jovian planets. Possible cloud-forming condensates considered are H2O ice, aqueous NH3 solution, NH4SH, and solid NH3, CH4, and Ar. Some techniques for computer calculation of condensation in a multiphase system of variable composition are discussed. Nominal models are computed for all four planets, and the effects of variations in composition and different pressure-temperature regimes are investigated in detail for Jupiter. The meaning of the calculated closed densities and their relation to the actual densities are discussed. The computed atmospheric profiles represent average conditions, subject to considerable local variations.

Weidenschilling, S. J.↗

Properties of aerosols in the Martian atmosphere, as inferred from Viking Lander imaging data

Three types of aerosol were detected from observations of the Martian sky, Phobos, and the sun with the Viking imaging cameras. Atmospheric optical depths were derived from observations of the objects' brightness. Data on the absorption coefficient, mean size, and shape of aerosols were obtained from studies of the sky brightness, using a multiple-scattering computer code. The aerosol types were water ice ground fog, a higher-level ice cloud (polar hood), and soil particles suspended at heights up to 30 km. The properties of each type of aerosol are discussed. The data are subjected to thorough analysis, and the results are summarized.

Pollack, J. B.↗

Investigation of water droplet trajectories within the NASA icing research tunnel

Water droplet trajectories within the NASA Lewis Research Center's Icing Research Tunnel (IRT) were studied through computer analysis. Of interest was the influence of the wind tunnel contraction and wind tunnel model blockage on the water droplet trajectories. The computer analysis was carried out with a program package consisting of a three-dimensional potential panel code and a three-dimensional droplet trajectory code. The wind tunnel contraction was found to influence the droplet size distribution and liquid water content distribution across the test section from that at the inlet. The wind tunnel walls were found to have negligible influence upon the impingement of water droplets upon a wing model.

Reehorst, Andrew↗

Cloud Microphysics and Absorption Validation

Vertical distributions of particle size and habit were developed from in-situ data collected from three midlatitude cirrus field campaigns (FIRE-1, FIRE-2, and ARM IOP). These new midlatitude microphysical models were used to develop new cirrus scattering models at a number of wavelengths appropriate for use with the MODIS imager (Nasiri et al. 2002). This was the first successful collaborative effort between all the investigators on this proposal. Recent efforts have extended the midlatitude cirrus cloud analyses to tropical cirrus, using in-situ data collected during the Tropical Rainfall Measurement Mission (TRMM) Kwajalein field campaign in 1999. We note that there are critical aspects to the work: a) Improvement in computing the scattering and radiative properties of ice crystals; b) Requirement for copious amounts of cirrus in-situ data, presented in terms of both particle size and habit distributions; c) Development of cirrus microphysical and optical models for various satellite, aircraft, and ground-based instruments based on the theoretical calculations and in-situ measurements; d) Application to satellite data.

Ackerman, Steven↗

LEWICE 2.2 Capabilities and Thermal Validation

A computational model of bleed air anti-icing and electrothermal de-icing have been added to the LEWICE 2.0 software by integrating the capabilities of two previous programs, ANTICE and LEWICE/ Thermal. This combined model has been released as LEWICE version 2.2. Several advancements have also been added to the previous capabilities of each module. This report will present the capabilities of the software package and provide results for both bleed air and electrothermal cases. A comprehensive validation effort has also been performed to compare the predictions to an existing electrothermal database. A quantitative comparison shows that for deicing cases, the average difference is 9.4 F (26%) compared to 3 F for the experimental data while for evaporative cases the average difference is 2 F (32%) compared to an experimental error of 4 F.

Wright, William B.↗

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

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

Hudson, Reggie L.↗

A fast particle-based approach for calibrating a 3-D model of the Antarctic ice sheet

We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here, we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high-performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters. For example, based on the same ice sheet model and data set, the 99th percentile of the Antarctic ice sheet contribution to sea level rise in 2300 increases from 6.5 m to 13.1 m when we increase the number of calibrated parameters from three to 11. With previous calibration methods, it would be challenging to go beyond five parameters. Here, this work provides an important next step toward improving the uncertainty quantification of complex, computationally intensive and decision-relevant models.

54 ENVIRONMENTAL SCIENCES↗

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

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

Hudson, Reggie L.↗