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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS2000

A Numerical Evaluation of Icing Effects on a Natural Laminar Flow Airfoil

As a part of CFD code validation efforts within the Icing Branch of NASA Glenn Research Center, computations were performed for natural laminar flow (NLF) airfoil, NLF-0414. with 6 and 22.5 minute ice accretions. Both 3-D ice castings and 2-D machine-generated ice shapes were used in wind tunnel tests to study the effects of natural ice is well as simulated ice. They were mounted in the test section of the Low Turbulence Pressure Tunnel (LTPT) at NASA Langley that the 2-dimensionality of the flow can be maintained. Aerodynamic properties predicted by computations were compared to data obtained through the experiment by the authors at the LTPT. Computations were performed only in 2-D and in the case of 3-D ice, the digitized ice shape obtained at one spanwise location was used. The comparisons were mainly concentrated on the lift characteristics over Reynolds numbers ranging from 3 to 10 million and Mach numbers ranging from 0.12 to 0.29. WIND code computations indicated that the predicted stall angles were in agreement with experiment within one or two degrees. The maximum lift values obtained by computations were in good agreement with those of the experiment for the 6 minute ice shapes and the minute 3-D ice, but were somewhat lower in the case of the 22.5 minute 2-D ice. In general, the Reynolds number variation did not cause much change in the lift values while the variation of Mach number showed more change in the lift. The Spalart-Allmaras (S-A) turbulence model was the best performing model for the airfoil with the 22.5 minute ice and the Shear Stress Turbulence (SST) turbulence model was the best for the airfoil with the 6 minute ice and also for the clean airfoil. The pressure distribution on the surface of the iced airfoil showed good agreement for the 6 minute ice. However, relatively poor agreement of the pressure distribution on the upper surface aft of the leading edge horn for the 22.5 minute ice suggests that improvements are needed in the grid or turbulence models.

Chung, James J.↗

NASA NTRS2005

The Influence of Viscous Effects on Ice Accretion Prediction and Airfoil Performance Predictions

A computational study was conducted to evaluate the effectiveness of using a viscous flow solution in an ice accretion code and the resulting accuracy of aerodynamic performance prediction. Ice shapes were obtained for one single-element and one multi-element airfoil using both potential flow and Navier-Stokes flowfields in the LEWICE ice accretion code. Aerodynamics were then calculated using a Navier-Stokes flow solver.

Kreeger, Richard E.↗

NASA NTRS1994

Arctic sea ice concentrations from special sensor microwave imager and advanced very high resolution radiometer satellite data

Nearly coincident data from the special sensor microwave imager (SSM/I) and the advanced very high resolution radiometer (AVHRR) are used to compute and compare Arctic sea ice concentrations for different regions and times of the year. To help determine overall accuracies and to highlight sources of differences between passive microwave, optical wavelength, and thermal wavelength data, ice concentrations are estimated using two operational SSM/I ice concentration algorithms and with visible- and thermal-infrared wavelength AVHRR data. All algorithms capture the seasonal patterns of ice growth and melt. The ranges of differences fall within the general levels of uncertainty expected for each method and are similar to previous accuracy estimates. The estimated ice concentrations are all highly correlated, with uniform biases, although differences between individual pairs of observations can be large. On average, the NASA Team algorithm yielded 5% higher ice concentrations than the Bootstrap algorithm, while during nonmelt periods the two SSM/I algorithms agree to within 0.5%. These seasonal differences are consistent with the ways that the 19-GHz and 37-GHz microwave channels are used in the algorithms. When compared to the AVHRR-derived ice concentrations, the Team-algorithm results are more similar on average in terms of correlation and mean differences. However, the Team algorithm underestimates concentrations relative to the AVHRR output by 6% during cold months and overestimates by 3% during summer. Little seasonal difference exists between the Bootstrap and AVHRR results, with a mean difference of about 5%. Although the mean differences are less between the SSM/I-derived concentrations and concentrations estimated using AVHRR channel 1, the correlations appear substantially better between the SSM/I data and concentrations derived from AVHRR channel 4, particularly for the Team algorithm output.

Emery, W. J.↗

NASA NTRS1973

Numerical studies of unsteady two dimensional subsonic flows using the ICE method

A numerical program was developed to compute transient compressible and incompressible laminar flows in two dimensions with multicomponent mixing and chemical reaction. The algorithm used the Los Alamos Scientific Laboratory ICE (Implicit Continuous-Fluid Eulerian) method as its base. The program can compute both high and low speed compressible flows. The numerical program incorporating the stabilization techniques was quite successful in treating both old and new problems. Detailed calculations of coaxial flow very close to the entry plane were possible. The program treated complex flows such as the formation and downstream growth of a recirculation cell. An implicit solution of the species equation predicted mixing and reaction rates which compared favorably with the literature.

Wieber, P. R.↗

NASA NTRS2000

Radiative Transfer and Satellite Remote Sensing of Cirrus Clouds Using FIRE-2-IFO Data

Under the support of the NASA grant, we have developed a new geometric-optics model (GOM2) for the calculation of the single-scattering and polarization properties for arbitrarily oriented hexagonal ice crystals. From comparisons with the results computed by the finite difference time domain (FDTD) method, we show that the novel geometric-optics can be applied to the computation of the extinction cross section and single-scattering albedo for ice crystals with size parameters along the minimum dimension as small as approximately 6. We demonstrate that the present model converges to the conventional ray tracing method for large size parameters and produces single-scattering results close to those computed by the FDTD method for size parameters along the minimum dimension smaller than approximately 20. We demonstrate that neither the conventional geometric optics method nor the Lorenz-Mie theory can be used to approximate the scattering, absorption, and polarization features for hexagonal ice crystals with size parameters from approximately 5 to 20. On the satellite remote sensing algorithm development and validation, we have developed a numerical scheme to identify multilayer cirrus cloud systems using AVHRR data. We have applied this scheme to the satellite data collected over the FIRE-2-IFO area during nine overpasses within seven observation dates. Determination of the threshold values used in the detection scheme are based on statistical analyses of these satellite data.

Source record↗

NASA NTRS2009

NASA Iced Aerodynamics and Controls Current Research

This slide presentation reviews the state of current research in the area of aerodynamics and aircraft control with ice conditions by the Aviation Safety Program, part of the Integrated Resilient Aircraft Controls Project (IRAC). Included in the presentation is a overview of the modeling efforts. The objective of the modeling is to develop experimental and computational methods to model and predict aircraft response during adverse flight conditions, including icing. The Aircraft icing modeling efforts includes the Ice-Contaminated Aerodynamics Modeling, which examines the effects of ice contamination on aircraft aerodynamics, and CFD modeling of ice-contaminated aircraft aerodynamics, and Advanced Ice Accretion Process Modeling which examines the physics of ice accretion, and works on computational modeling of ice accretions. The IRAC testbed, a Generic Transport Model (GTM) and its use in the investigation of the effects of icing on its aerodynamics is also reviewed. This has led to a more thorough understanding and models, both theoretical and empirical of icing physics and ice accretion for airframes, advanced 3D ice accretion prediction codes, CFD methods for iced aerodynamics and better understanding of aircraft iced aerodynamics and its effects on control surface effectiveness.

Addy, Gene↗

NASA NTRS2008

Comparison of LEWICE and GlennICE in the SLD Regime

A research project is underway at the NASA Glenn Research Center (GRC) to produce computer software that can accurately predict ice growth under any meteorological conditions for any aircraft surface. This report will present results from two different computer programs. The first program, LEWICE version 3.2.2, has been reported on previously. The second program is GlennICE version 0.1. An extensive comparison of the results in a quantifiable manner against the database of ice shapes that have been generated in the GRC Icing Research Tunnel (IRT) has also been performed, including additional data taken to extend the database in the Super-cooled Large Drop (SLD) regime. This paper will show the differences in ice shape between LEWICE 3.2.2, GlennICE, and experimental data. This report will also provide a description of both programs. Comparisons are then made to recent additions to the SLD database and selected previous cases. Quantitative comparisons are shown for horn height, horn angle, icing limit, area, and leading edge thickness. The results show that the predicted results for both programs are within the accuracy limits of the experimental data for the majority of cases.

Wright, William B.↗

NASA NTRS2006

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗

DOE OSTI2025

An evaluation of multi-fidelity methods for quantifying uncertainty in projections of ice-sheet mass change

Abstract. This study investigated the computational benefits of using multi-fidelity statistical estimation (MFSE) algorithms to quantify uncertainty in the mass change of Humboldt Glacier, Greenland, between 2007 and 2100 using a single climate change scenario. The goal of this study was to determine whether MFSE can use multiple models of varying cost and accuracy to reduce the computational cost of estimating the mean and variance of the projected mass change of a glacier. The problem size and complexity were chosen to reflect the challenges posed by future continental-scale studies while still facilitating a computationally feasible investigation of MFSE methods. When quantifying uncertainty introduced by a high-dimensional parameterization of the basal friction field, MFSE was able to reduce the mean-squared error in the estimates of the statistics by well over an order of magnitude when compared to a single-fidelity approach that only used the highest-fidelity model. This significant reduction in computational cost was achieved despite the low-fidelity models used being incapable of capturing the local features of the ice-flow fields predicted by the high-fidelity model. The MFSE algorithms were able to effectively leverage the high correlation between each model's predictions of mass change, which all responded similarly to perturbations in the model inputs. Consequently, our results suggest that MFSE could be highly useful for reducing the cost of computing continental-scale probabilistic projections of sea-level rise due to ice-sheet mass change.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1985

Progress in development of a Navier-Stokes solver for evaluation of iced airfoil performance

A method is being developed for evaluation of the flow field behavior about an airfoil with significant ice accretion on the leading edge. The computer code, being evaluated for this purpose, solves the Navier-Stokes equations in a body-fitted curvilinear coordinate system. This requires the use of a grid generation code to transform the x-y coordinates of the physical space into xi-eta coordinates of the computational space. Evaluation of the suitability of these two codes for predicting iced airfoil performance is presently being carried out in anticipation of use in an overall icing analysis effort. Results of this evaluation to date indicate good correlation with known information on clean airfoils. Preliminary results for rime and glaze, iced airfoil shapes are also presented.

Potapczuk, M. G.↗

DOE OSTI2022

Simulating ice-shelf extent using damage mechanics

Abstract Inaccurate representations of iceberg calving from ice shelves are a large source of uncertainty in mass-loss projections from the Antarctic ice sheet. Here, we address this limitation by implementing and testing a continuum damage-mechanics model in a continental scale ice-sheet model. The damage-mechanics formulation, based on a linear stability analysis and subsequent long-wavelength approximation of crevasses that evolve in a viscous medium, links damage evolution to climate forcing and the large-scale stresses within an ice shelf. We incorporate this model into the BISICLES ice-sheet model and test it by applying it to idealized (1) ice tongues, for which we present analytical solutions and (2) buttressed ice-shelf geometries. Our simulations show that the model reproduces the large disparity in lengths of ice shelves with geometries and melt rates broadly similar to those of four Antarctic ice shelves: Erebus Glacier Tongue (length ~ 13 km), the unembayed portion of Drygalski Ice Tongue (~ 65 km), the Amery Ice Shelf (~ 350 km) and the Ross Ice Shelf (~ 500 km). These results demonstrate that our simple continuum model holds promise for constraining realistic ice-shelf extents in large-scale ice-sheet models in a computationally tractable manner.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2000

ICEG2D: An Integrated Software Package for Automated Prediction of Flow Fields for Single-Element Airfoils with Ice Accretion

An integrated software package, ICEG2D, was developed to automate computational fluid dynamics (CFD) simulations for single-element airfoils with ice accretion. ICEG2D is designed to automatically perform three primary functions: (1) generating a grid-ready, surface definition based on the geometrical characteristics of the iced airfoil surface, (2) generating a high-quality grid using the generated surface point distribution, and (3) generating the input and restart files needed to run the general purpose CFD solver NPARC. ICEG2D can be executed in batch mode using a script file or in an interactive mode by entering directives from a command line. This report summarizes activities completed in the first year of a three-year research and development program to address issues related to CFD simulations for aircraft components with ice accretion. Specifically, this document describes the technology employed in the software, the installation procedure, and a description of the operation of the software package. Validation of the geometry and grid generation modules of ICEG2D is also discussed.

Thompson, David S.↗

NASA NTRS1989

Late Pleistocene and Holocene Glacial Evolution and Isotasy in the Antarctic Peninsula

Employing a numerical model of Payne et al. that simulates the late-Pleistocene evolution of the former Antarctic Peninsula Ice Sheet (APIS) as a basis, we compute the present-day postglacial vertical isostasy of this region. The region may also experience significant mid-to late-Holocene glacial mass changes. Climate and oceanographic studies indicate that the ice mass imbalance of this region may be of larger magnitude that elsewhere in Antarctica. We compute the crustal response to these more recent ice mass changes and Holocene fluctuations with a simple gravitating Earth model consisting of an elastic lithosphere and a viscoelastic mantle (half-space). The calculations demonstrate that the present-day response could be significant, possibly at the level of about 4 - 11 mm/yr. Such significant crustal motion could be driven by glacial mass changes integrated over the last 1000 years if the regional mantle viscosity is below about 2 x 10(exp 20) Pa sec. In this lower viscosity range, present-day crustal motion has a significant phase-lagged character and the composite lithosphere/mantle viscoelastic response to late-Holocene events dominates over purely elastic (instantaneous) responses to present-day ice mass changes. For a higher mantle viscosity, greater than about 5 x 10(exp 20) Pa sec, the predicted present-day vertical isostasy is dominated by gravitational response to glacial unloading during the 18 - 6 kyr BP collapse of the APIS, and is analogous to that known to be occurring in the Gulf of Bothnia and Hudson Bay.

Ivins, E. R.↗

NASA NTRS1994

A hybrid Navier-Stokes/full-potential method for the prediction of iced wing aerodynamics

A hybrid method for computing compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. The resulting CPU times are less than 50 percent of the required for a full-blown Navier-Stokes analysis. Sample applications of the method to an unswept iced wing at 4 deg and 8 deg angle of attack are presented. Surface pressures are in good agreement with the measurements obtained by Bragg et al. at the University of Illinois.

Mello, O. A. F.↗

NASA NTRS2014

Modeling the Deterioration of Engine and Low Pressure Compressor Performance During a Roll Back Event Due to Ice Accretion

The main focus of this study is to apply a computational tool for the flow analysis of the engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) of NASA Glenn Research Center. A data point was selected for analysis during which the engine experienced a full roll back event due to the ice accretion on the blades and flow path of the low pressure compressor. The computational tool consists of the Numerical Propulsion System Simulation (NPSS) engine system thermodynamic cycle code, and an Euler-based compressor flow analysis code, that has an ice particle melt estimation code with the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Decreasing the performance characteristics of the low pressure compressor (LPC) within the NPSS cycle analysis resulted in matching the overall engine performance parameters measured during testing at data points in short time intervals through the progression of the roll back event. Detailed analysis of the fan-core and LPC with the compressor flow analysis code simulated the effects of ice accretion by increasing the aerodynamic blockage and pressure losses through the low pressure compressor until achieving a match with the NPSS cycle analysis results, at each scan. With the additional blockages and losses in the LPC, the compressor flow analysis code results were able to numerically reproduce the performance that was determined by the NPSS cycle analysis, which was in agreement with the PSL engine test data. The compressor flow analysis indicated that the blockage due to ice accretion in the LPC exit guide vane stators caused the exit guide vane (EGV) to be nearly choked, significantly reducing the air flow rate into the core. This caused the LPC to eventually be in stall due to increasing levels of diffusion in the rotors and high incidence angles in the inlet guide vane (IGV) and EGV stators. The flow analysis indicating compressor stall is substantiated by the video images of the IGV taken during the PSL test, which showed water on the surface of the IGV flowing upstream out of the engine, indicating flow reversal, which is characteristic of a stalled compressor.

Veres, Joseph P.↗

NASA NTRS2014

Modeling the Deterioration of Engine and Low Pressure Compressor Performance During a Roll Back Event Due to Ice Accretion

Jorgenson, Philip C. E.↗

NASA NTRS2021

Evaluation of Higher-order Quadrature Schemes in Improving Computational Efficiency for Orientation-averaged Single-Scattering Properties of Nonspherical Ice Particles

We evaluate several high-order quadrature schemes for accuracy and efficacy in obtaining orientation-averaged single-scattering properties (SSPs). We use the highly efficient MIDAS to perform electromagnetic scattering calculations to evaluate the gain in efficiency from these schemes. MIDAS is shown to be superior to DDSCAT, a popular discrete dipole approximation (DDA) method. This study is motivated by the fact that quality physical precipitation retrievals rely on using accurate orientation-averaged SSPs derived from realistic hydrometeors as input to radiative transfer Models (RTMs). The DDA has been a popular choice for single-scattering calculations, due to its versatility with respect to target geometry. However, being iterative-solver-based (ISB), the most used DDA codes, e.g. DDSCAT and ADDA, must solve the scattering problem for each orientation of the target separately. As the size parameter and geometric anisotropy of the hydrometeor increase, the number of orientations needed to obtain accurate orientation-averages can increase drastically and so does the computation cost incurred by the ISB-DDA methods. MIDAS is a Direct-Solver-Based (DSB) code, its decomposition of the original large matrix with a high rank into multiple more manageable smaller matrices of lower ranks makes it much more computationally efficient while maintaining excellent accuracy. In addition, direct solvers consider all requested orientations at once, giving MIDAS further advantage over popular ISB-DDA methods. MIDAS, when combined with high-order quadrature for orientation averaging, can be greater than three orders of magnitude more efficient in obtaining RTM-ready SSPs of complex-shaped hydrometeors than existing ISB-DDA methods, with the native quadrature schemes they offer.

Ines Fenni↗

NASA NTRS1987

Numerical analysis of a NACA0012 airfoil with leading edge ice accretions

Analysis of a NACA0012 airfoil with leading edge ice has been performed using a Navier-Stokes code coupled with a grid generation code. The computed results were compared to experimental information obtained for an airfoil with a well defined artificial ice shape. The computations were performed at angles of attack ranging from zero to ten degrees. This range is sufficient to show the development of the separation bubble aft of the ice shape on both the upper and lower surfaces. Velocity profile plots in the separation bubble are examined in order to determine if recirculation patterns are predicted properly and if separation and reattachment points are found within the resolution of the experimental information. Also, the massive separation near the point of stall is examined in order to more accurately evaluate the lift coefficient curve in that region. Lift, drag, and moment coefficients are computed and compared to experiment.

Potapczuk, Mark G.↗