SEARCH · Search NASA

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Cirrus optical depth retrieval from satellite-measured bidirectional reflectance

In the framework of investigating the relationship between the visible optical depth and microphysics of cirrus clouds, the doubling scheme for radiative transfer in cirrus cloudy atmospheres, as developed by Takano and Liou (1988), has been used to compute the bidirectional reflectances at the top of the atmosphere. Reflectance data taken at about the same time from GOES and AVHRR over Waussau, as well as particle distributions observed during the FIRE IFO carried out in Wisconsin (1986) were the basis of the experiment. For the cirrus microphysical model, hexagonal ice crystals and area-equivalent ice spheres are considered. The single scattering phase functions for randomly oriented hexagonal columns and area-equivalent ice spheres are computed from geometric optics and Mie theory, respectively.

Takano, Y.↗

NASA NTRS1992

Comparison of two-dimensional and three-dimensional droplet trajectory calculations in the vicinity of finite wings

Computational predictions of ice accretion on flying aircraft most commonly rely on modeling in 2D. These 2D methods treat an aircraft geometry either as wing-like with infinite span, or as an axisymmetric body. Recently, fully 3D methods have been introduced that model an aircraft's true 3D shape. Because 3D methods are more computationally expensive than 2D methods, 2D methods continue to be widely used. However, a 3D method allows investigation of whether it is valid to continue applying 2D methods to a finite wing. The extent of disagreement between LEWICE, a 2D method, and LEWICE3D, a 3D method, in calculating local collection efficiencies at the leading edge of finite wings is investigated.

Mohler, Stanley R., Jr.↗

NASA NTRS1992

Comparison of two-dimensional and three-dimensional droplet trajectory calculations in the vicinity of finite wings

Computational predictions of ice accretion on flying aircraft most commonly rely on modeling in two dimensions (2D). These 2D methods treat an aircraft geometry either as wing-like with infinite span, or as an axisymmetric body. Recently, fully three dimensional (3D) methods have been introduced that model an aircrafts true 3D shape. Because 3D methods are more computationally expensive than 2D methods, 2D methods continue to be widely used. However, a 3D method allows us to investigate whether it is valid to continue applying 2D methods to a finite wing. The extent of disagreement between LEWICE, a 2D method, and LEWICE3D, a 3D method, in calculating local collection efficiencies at the leading edge of finite wings is investigated in this paper.

Mohler, Stanley R., Jr.↗

NASA NTRS2003

Space/Time Statistics of Polar Ice Motion

Ice motions have been computed from passive microwave imagery (SMMR and SSM/I) on a daily basis for both Polar Regions. In the Arctic these daily motions have been merged with daily motions from AVHRR imagery and the Arctic buoy program. In the Antarctic motion only from the AVHRR were available for merging with the passive microwave vectors. Long-term means, monthly means and weekly means have all been computed from the resulting 22-year time series of polar ice motion. Papers are in preparation that present the long term (22 year) means, their variability and show animations of the monthly means over this time period for both Polar Regions. These papers will have links to "enhanced objects" that allow the reader to view the animations as part of the paper. The first paper presents the ice motion results from each of the Polar Regions. The second paper looks only at ice motion in the Arctic in order to develop a time series of ice age in the Arctic. Starting with the first full SMMR year in 1979 we keep track of each individual "ice element" (resolution of the sensor) and track it in the subsequent monthly time series. After a year we "age" each "particle" and we thus can keep track of the age of the ice starting in 1979. We keep track of ice age classes between one and five years and thus we can see the evolution of the ice as it ages after the initial 5-year period. This calculation shows how we are losing the older ice through Fram Strait at a rather alarming rate particularly in the past 15 years. This loss of older ice has resulted in an overall decrease in the thickest, oldest ice, which is now limited to a region just north of the Canadian Archipelago with tongues extending out across the pole towards the Siberian Shelf. This loss of old ice is consistent with the effects of global warming which provides the heat needed to melt, move and disperse this oldest ice through Fram Strait. This is the first step in a progression that may eventually open the Arctic ice pack and lead to an ice-free Arctic Ocean.

Emery, William J.↗

NASA NTRS1975

The size dependence of sublimation rates for interplanetary ice particles

The sublimation rates for water ice have been computed as a function of particle size for various solar distances. Because of the size dependence of the absorption and emission properties of the particles, a sublimation-rate minimum evolves whose depth and position are sensitive to the spectral-absorption properties of the particle in combination with the spectral distribution of solar radiation. As a consequence, a quasistable size of interplanetary ice particles is predicted which is independent of solar distance.

Patashnick, H.↗

NASA NTRS1996

Development of Three-Dimensional Flow Code Package to Predict Performance and Stability of Aircraft with Leading Edge Ice Contamination

In the work reported herein, a simplified, uncoupled, zonal procedure is utilized to assess the capability of numerically simulating icing effects on a Boeing 727-200 aircraft. The computational approach combines potential flow plus boundary layer simulations by VSAERO for the un-iced aircraft forces and moments with Navier-Stokes simulations by NPARC for the incremental forces and moments due to iced components. These are compared with wind tunnel force and moment data, supplied by the Boeing Company, examining longitudinal flight characteristics. Grid refinement improved the local flow features over previously reported work with no appreciable difference in the incremental ice effect. The computed lift curve slope with and without empennage ice matches the experimental value to within 1%, and the zero lift angle agrees to within 0.2 of a degree. The computed slope of the un-iced and iced aircraft longitudinal stability curve is within about 2% of the test data. This work demonstrates the feasibility of a zonal method for the icing analysis of complete aircraft or isolated components within the linear angle of attack range. In fact, this zonal technique has allowed for the viscous analysis of a complete aircraft with ice which is currently not otherwise considered tractable.

Aircraft icing↗

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 NTRS1982

Mathematical modeling of ice accretion on airfoils

The progress toward development of a computer model suitable for predicting icing behavior on airfoils over a wide range of environmental conditions and airfoils shapes is reported. The LEWICE program was formulated to solve a set of equations which describe the physical processes which occur during accretion of ice on an airfoil, including heat transfer in a time dependent mode, with the restriction that the flow must be describable by a two-dimensional flow code. Input data comprises the cloud liquid water content, mean droplet diameter, ambient air temperature, air velocity, and relative humidity. A potential flowfield around the airfoil is calculated, along with the droplet trajectories within the flowfield, followed by local values of water droplet collection efficiency at the impact points. Both glaze and rime ice conditions are reproduced, and comparisons with test results on icing of circular cylinders showed good agreement with the physical situation.

Macarthur, C. D.↗

NASA NTRS2002

Experimental Investigation of Water Droplet Impingement on Airfoils, Finite Wings, and an S-duct Engine Inlet

Validation of trajectory computer codes, for icing analysis, requires experimental water droplet impingement data for a wide range of aircraft geometries as well as flow and icing conditions. This report presents improved experimental and data reduction methods for obtaining water droplet impingement data and provides a comprehensive water droplet impingement database for a range of test geometries including an MS(1)-0317 airfoil, a GLC-305 airfoil, an NACA 65(sub 2)-415 airfoil, a commercial transport tail section, a 36-inch chord natural laminar flow NLF(1)-0414 airfoil, a 48-inch NLF(1)-0414 section with a 25 percent chord simple flap, a state-of-the-art three-element high lift system, a NACA 64A008 finite span swept business jet tail, a full-scale business jet horizontal tail section, a 25 percent-scale business jet empennage, and an S-duct turboprop engine inlet. The experimental results were obtained at the NASA Glenn Icing Research Tunnel (IRT) for spray clouds with median volumetric diameter (MVD) of 11, 11.5, 21, 92, and 94 microns and for a range of angles of attack. The majority of the impingement experiments were conducted at an air speed of 175 mph corresponding to a Reynolds number of approximately 1.6 million per foot. The maximum difference of repeated tests from the average ranged from 0.24 to 12 percent for most of the experimental results presented. This represents a significant improvement in test repeatability compared to previous experimental studies. The increase in test repeatability was attributed to improvements made to the experimental and data reduction methods. Computations performed with the LEWICE-2D and LEWICE-3D computer codes for all test configurations are presented in this report. For the test cases involving median volumetric diameters of 11 and 21 microns, the correlation between the analytical and experimental impingement efficiency distributions was good. For the median volumetric diameters of 92 and 94-micron cases, however, the analysis produced higher impingement efficiencies and larger impingement limits than the experiment. It is speculated that this discrepancy is due to droplet splashing and breakup experienced by large droplets during impingement.

Papadakis, Michael↗

NASA NTRSDate not supplied

An Automated Refinement Process for Particle Trajectory Methods in GlennICE

Computational methods for ice accretion can simulate the impact of water drops and ice crystals on an aircraft surface in a Lagrangian reference frame or in the Eulerian reference frame. In the Eulerian reference frame, particles are considered a continuous fluid while in the Lagrangian frame individual particle trajectories are calculated. Methods that use the Eulerian reference frame are typically easier to develop as established modules used for continuum mechanics can be leveraged. The Eulerian systems can also be faster since the user does not have to simulate millions of particles in order to achieve good results. It is imperative therefore that a Lagrangian method optimize the release points of trajectories such that accurate solutions can be obtained while minimizing as much as possible the number of trajectories computed. This paper will present a methodology for this refinement process and demonstrate its effectiveness on sample three dimensional test cases.

William B Wright↗

NASA NTRSDate not supplied

An Automated Refinement Process for Particle Trajectory Methods in GlennICE

William Wright↗

NASA NTRS2018

Implementation of a Comprehensive Ice Crystal Formation Parameterization for Cirrus and Mixed-Phase Clouds in the EMAC Model (Based on MESSy 2.53)

A comprehensive ice nucleation parameterization has been implemented in the global chemistry-climate model EMAC to realistically represent ice crystal number concentrations. The parameterization of Barahona and Nenes (2009, hereafter BN09) allows the treatment of ice nucleation, taking into account the competition for water vapour between homogeneous and heterogeneous nucleation and pre-existing ice crystals in cold clouds. Furthermore, the influence of chemically-heterogeneous, polydisperse aerosols is considered via multiple ice nucleating particle spectra, which are included in the parameterization to compute the heterogeneously formed ice crystals. BN09 has been implemented to operate both in the cirrus and in the mixed-phase regimes. Compared to the standard EMAC results, BN09 produces fewer ice crystals in the upper troposphere but higher ice crystal number concentrations in the middle troposphere, especially in the Northern Hemisphere where ice nucleating mineral dust particles are relatively abundant. The comparison with a climatological data set of aircraft measurements shows that BN09 used in the cirrus regime improves the model results and, therefore, is recommended for future EMAC simulations.

EMAC↗

NASA NTRS2014

Assessment of the Accuracy of the Conventional Ray-Tracing Technique: Implications in Remote Sensing and Radiative Transfer Involving Ice Clouds.

A fundamental problem in remote sensing and radiative transfer simulations involving ice clouds is the ability to compute accurate optical properties for individual ice particles. While relatively simple and intuitively appealing, the conventional geometric-optics method (CGOM) is used frequently for the solution of light scattering by ice crystals. Due to the approximations in the ray-tracing technique, the CGOM accuracy is not well quantified. The result is that the uncertainties are introduced that can impact many applications. Improvements in the Invariant Imbedding T-matrix method (II-TM) and the Improved Geometric-Optics Method (IGOM) provide a mechanism to assess the aforementioned uncertainties. The results computed by the II-TMþIGOM are considered as a benchmark because the IITM solves Maxwell's equations from first principles and is applicable to particle size parameters ranging into the domain at which the IGOM has reasonable accuracy. To assess the uncertainties with the CGOM in remote sensing and radiative transfer simulations, two independent optical property datasets of hexagonal columns are developed for sensitivity studies by using the CGOM and the II-TMþIGOM, respectively. Ice cloud bulk optical properties obtained from the two datasets are compared and subsequently applied to retrieve the optical thickness and effective diameter from Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Additionally, the bulk optical properties are tested in broadband radiative transfer (RT) simulations using the general circulation model (GCM) version of the Rapid Radiative Transfer Model (RRTMG) that is adopted in the National Center for Atmospheric Research (NCAR) Community Atmosphere Model (CAM, version 5.1). For MODIS retrievals, the mean bias of uncertainties of applying the CGOM in shortwave bands (0.86 and 2.13 micrometers) can be up to 5% in the optical thickness and as high as 20% in the effective diameter, depending on cloud optical thickness and effective diameter. In the MODIS infrared window bands centered at 8.5, 11, and 12 micrometers biases in the optical thickness and effective diameter are up to 12% and 10%, respectively. The CGOM-based simulation errors in ice cloud radiative forcing calculations are on the order of 10Wm(exp 2).

remote sensing↗

NASA NTRS1987

Studies of the scattering/absorption of minerals

Reflectance spectra were computed for water ice and ammonia ice mixtures as functions of weight fraction, grain size, and viewing geometry to simulate possible outer solar system satellite surfaces. Reflectance spectra of planetary surfaces are most affected by the weight fraction and grain sizes of the minerals in the surface. The reflectance can range from 1.0 to about 0.01 by changing the grain size or weight fraction, a factor of 100. Viewing geometry changes the reflectance by about 25 percent or less.

Clark, Roger N.↗

NASA NTRS1991

Statistical structural analysis of rotor impact ice shedding

The statistical characteristics of impact ice shear strength are analyzed, with emphasis placed on the most probable shear strength and statistical distribution of an ice deposit. Several distribution types are considered: the Weibull, two-parameter Weibull, and exponential distributions, as well as the Gumbell distribution of the smallest extreme and the Gumbell distribution of the largest extreme. It is concluded that the Weibull distribution yields the best results; however, the expected life, shape parameter, and scale parameter should be determined separately for each case of varying wind speed and droplet size. The theoretical predictions of shear stresses in a specific rotating ice shape are compared, and it is noted that when the effects of lift are added to the theoretical model and the interference is calculated with a new mean and standard deviation, the probability of ice shed is computed as 36.64 pct.

Kellacky, C. J.↗

NASA NTRS1991

Reports to the Nation on our changing planet. Winter 1991 No. 1: The climate system

This pamphlet is a semi-annual report for the period April 1992 through September 1992 dealing with the study of climatology. This introductory pamphlet treats in the simplest of terms such topics as climate change, the greenhouse effect, the global heat engine, hydrology, oceanography, computers and climatology, ice ages, and global warming. Informative charts are included.

Source record↗

NASA NTRS1993

Blade Assessment for Ice Impact (BLASIM). User's manual, version 1.0

The Blade Assessment Ice Impact (BLASIM) computer code can analyze solid, hollow, composite, and super hybrid blades. The solid blade is made up of a single material where hollow, composite, and super hybrid blades are constructed with prescribed composite layup. The properties of a composite blade can be specified by inputting one of two options: (1) individual ply properties, or (2) fiber/matrix combinations. When the second option is selected, BLASIM utilizes ICAN (Integrated Composite ANalyzer) to generate the temperature/moisture dependent ply properties of the composite blade. Two types of geometry input can be given: airfoil coordinates or NASTRAN type finite element model. These features increase the flexibility of the program. The user's manual provides sample cases to facilitate efficient use of the code while gaining familiarity.

Reddy, E. S.↗

NASA NTRS1985

Full two-dimensional transient solutions of electrothermal aircraft blade deicing

Two finite difference methods are presented for the analysis of transient, two-dimensional responses of an electrothermal de-icer pad of an aircraft wing or blade with attached variable ice layer thickness. Both models employ a Crank-Nicholson iterative scheme, and use an enthalpy formulation to handle the phase change in the ice layer. The first technique makes use of a 'staircase' approach, fitting the irregular ice boundary with square computational cells. The second technique uses a body fitted coordinate transform, and maps the exact shape of the irregular boundary into a rectangular body, with uniformally square computational cells. The numerical solution takes place in the transformed plane. Initial results accounting for variable ice layer thickness are presented. Details of planned de-icing tests at NASA-Lewis, which will provide empirical verification for the above two methods, are also presented.

Masiulaniec, K. C.↗