Search NASA⌕ Search

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.

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

At least 91 records · Page 5

Icing simulation: A survey of computer models and experimental facilities

A survey of the current methods for simulation of the response of an aircraft or aircraft subsystem to an icing encounter is presented. The topics discussed include a computer code modeling of aircraft icing and performance degradation, an evaluation of experimental facility simulation capabilities, and ice protection system evaluation tests in simulated icing conditions. Current research focussed on upgrading simulation fidelity of both experimental and computational methods is discussed. The need for increased understanding of the physical processes governing ice accretion, ice shedding, and iced airfoil aerodynamics is examined.

Potapczuk, M. G.↗

Icing simulation: A survey of computer models and experimental facilities

A survey of the current methods for simulation of the response of an aircraft or aircraft subsystem to an icing encounter is presented. The topics discussed include a computer code modeling of aircraft icing and performance degradation, an evaluation of experimental facility simulation capabilities, and ice protection system evaluation tests in simulated icing conditions. Current research focused on upgrading simulation fidelity of both experimental and computational methods is discussed. The need for the increased understanding of the physical processes governing ice accretion, ice shedding, and iced aerodynamics is examined.

Potapczuk, M. G.↗

Computational Heat Flux Analysis for HEATheR Electric Aircraft Wing Ice Protection Systems

Increased electrification of aircraft systems, a growing trend in the aircraft industry, affects subsystems that have historically been pneumatically driven by the engine. Optimized design of propulsion and power systems must account for the requirements of these subsystems. One such subsystem, the wing ice protection system, is the focus of this paper. A preliminary computational icing analysis has been conducted using the LEWICE and ANSYS FENSAP-ICE software to characterize the requirements of a thermal ice protection system for three More Electric Aircraft concepts. This paper summarizes the methodology, simulation execution, and icing conditions used for the analysis as well as discusses the results, observed trends, and impacts to system design. This initial analysis indicates that the computed heat requirements for the wing ice protection systems are not within the capabilities of known traditional systems. A discussion of potential methods for reducing the heat flow requirement by the ice protection system is included. Potential refinements to the analysis to improve the accuracy are also discussed.

Nic Heersema↗

LEWICE/E: An Euler based ice accretion code

A new version of the LEWICE ice accretion computer code was developed which calculates the ice growth on two dimensional surfaces, incorporating the effects of compressibility through the solution of the Euler equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, and plot results. This code increases the applicability of ice accretion predictions by allowing calculations at higher Mach numbers. The new elements of the code are described. Calculated results are compared to experiment for several cases, including a LEWICE example case and a thin airfoil section at a Mach number of 0.58.

Potapczuk, Mark G.↗

LEWICE/E - An Euler based ice accretion code

A new version of the LEWICE ice accretion computer code was developed which calculates the ice growth on two dimensional surfaces, incorporating the effects of compressibility through the solution of the Euler equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, and plot results. This code increases the applicability of ice accretion predictions by allowing calculations at higher Mach numbers. The new elements of the code are described. Calculated results are compared to experiment for several cases, including a LEWICE example case and a thin airfoil section at a Mach number of 0.58.

Potapczuk, Mark G.↗

Analytical determination of propeller performance degradation due to ice accretion

A computer code capable of computing the propeller performance for clean, glaze, or rime ice propeller configurations to determine the performance degradation resulting from a given icing encounter has been developed. The inviscid, incompressible flowfield at each specified propeller radial location is first computed using the Theodorsen method. A droplet trajectory computation then calculates the droplet impingement points and airfoil collection efficiency for each radial location. User-selectable empirical correlations are available for determining the aerodynamic penalties due to ice accretion. Propeller performance is finally computed using strip analysis for either the clean or iced propeller. In the iced mode, the thrust and torque coefficient equations are modified by the drag and lift coefficient increments due to ice to obtain the appropriate iced values. Comparison with available experimental propeller icing data shows generally good agreement. The code's capability of properly predicting the thrust coefficient, power coefficient, and propeller efficiency of an iced propeller is shown to be dependent on the choice of empirical correlation employed as well as on the proper specification of the radial icing extent and propeller blade angle.

Miller, T. L.↗

Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation

The ice accrued on airfoils can have enormously complicated shapes with multiple protruded horns and feathers. In this paper, several blocking topologies are proposed and evaluated on their ability to produce high-quality structured multi-block grid systems. A transition layer grid is introduced to ensure that jaggedness on the ice-surface geometry do not to propagate into the domain. This is important for grid-generation methods based on hyperbolic PDEs (Partial Differential Equations) and algebraic transfinite interpolation. A 'thick' wrap-around grid is introduced to ensure that grid lines clustered next to solid walls do not propagate as streaks of tightly packed grid lines into the interior of the domain along block boundaries. For ice shapes that are not too complicated, a method is presented for generating high-quality single-block grids. To demonstrate the usefulness of the methods developed, grids and CFD solutions were generated for two iced airfoils: the NLF0414 airfoil with and without the 623-ice shape and the B575/767 airfoil with and without the 145m-ice shape. To validate the computations, the computed lift coefficients as a function of angle of attack were compared with available experimental data. The ice shapes and the blocking topologies were prepared by NASA Glenn's SmaggIce software. The grid systems were generated by using a four-boundary method based on Hermite interpolation with controls on clustering, orthogonality next to walls, and C continuity across block boundaries. The flow was modeled by the ensemble-averaged compressible Navier-Stokes equations, closed by the shear-stress transport turbulence model in which the integration is to the wall. All solutions were generated by using the NPARC WIND code.

Chi, X.↗

Normality of I-V Measurements Using ML

There is an increased interest in instrument-computing ecosystems (ICEs) that support science workflows empowered by AI-automated experiments and computations in diverse areas. In particular, electrochemistry ICEs are promising for accelerating the design and discovery of electrochemical systems for energy storage and conversion, by automating significant parts of workflows that combine synthesis and characterization experiments with computations. They require the integration of flow controllers, solvent containers, pumps, fraction collectors, and potentiostats, all connected to an electrochemical cell, as illustrated in Fig. 1. These are specialized instruments with custom software that is not originally designed for network integration. We developed network and software solutions for electrochemical workflows that adapt system and instrument settings in real-time for multiple rounds of experiments. In particular, we developed Python wrappers for Application Programming Interfaces (APIs) of instrument commands and Pyro client-server modules that enable them to be executed from remote computers. The entire workflow is orchestrated by a Jupyter notebook running on a remote computer.

Al Najjar, Anees↗

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Center's Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Centers Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗

Shear zone ice deformation using supervised analysis of Seasat data

Spacecraft and aircraft Synthetic-Aperture Radar (SAR) images are commonly used to compute sea ice conditions including especially the drift velocities of ice features. Recent developments, described in this paper, in computer-aided analysis and in digital image processing have extended the analysis speed and accuracy of registration to the point that ice deformation can be examined on the geophysically important scales of about 100 m over tens of kilometers between two 'snapshots.' Examples are shown for the ice motion and rotation in the shear zone near the Alaska coast from Seasat data in 1978.

Carsey, F.↗

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Icing Physics Studies Using the 3D SIDRM Test Article: Ice Crystal Icing Analysis

In-flight icing is an important safety issue and is a factor that affects aircraft design and performance. Newer regulations are driving a need for improvements in airframe and engine icing simulation capability. Experimental data is required for the development of icing physics models and simulation validation. This paper presents the analysis of the ice crystal icing data subset from tests conducted in 2022 at the NASA Icing Research Tunnel that studied both supercooled liquid and ice-crystal icing. The test article that was utilized replicated 3D geometrical features of an inter-compressor duct and strut region of a turbofan engine. The surfaces of the Simulated Inter-compressor Duct Research Model (SIDRM) can be heated to simulate the warm surfaces of the turbofan inter-compressor duct. The test article is instrumented with pressure taps, heaters, heat flux gauges, and thermocouples, while a 3D laser scanner, cameras, and a weighing scale to measure ice mass were utilized to characterize the icing behavior. The aim of these tests was to generate ice accretions on the SIDRM test article under well-characterized icing conditions. To that end, ice crystal cloud characterization tests were conducted that provided cloud property values at the test section. The ice crystal accretion tests investigated how different test condition parameters impacted ice accretion size, location, and characteristics (such as physical attributes). A key finding was that a sufficiently high surface heat flux at the leading edge can provide a continuous supply of liquid melt to critical accretion zones downstream, resulting in elevated ice mass accretions. Another key finding was that warmer heater settings in critical accretion zones can suppress ice accretion. These findings are consistent with icing behavior observed in full-scale engine icing tests conducted previously at the NASA Glenn Research Center. Test runs investigating the contribution of liquid melt generated at the leading edge found that the melt can be transported downstream to critical accretion regions in two forms, via runback and splash. The form of melt supplied downstream affected accretion behavior and features. Sharkteeth-shaped ice accretions were generated at certain test conditions, a physical attribute that has been observed in full-scale engine icing tests. The ice accretion data will be used to develop and validate 3D computational engine icing tools, such as GlennICE, that predictively assesses the onset and growth of ice. One of the goals of the sponsoring NASA project is to develop simulation models and tools that can assist in the design and certification of engines for flight in icing conditions in a cost‑effective way.

ice crystal icing↗

Icing Physics Studies Using the 3D SIDRM Test Article: Ice Crystal Icing Analysis

In-flight icing is an important safety issue and is a factor that affects aircraft design and performance. Newer regulations are driving a need for improvements in airframe and engine icing simulation capability. Experimental data is required for the development of icing physics models and simulation validation. This paper presents the analysis of the ice crystal icing data subset from tests conducted in 2022 at the NASA Icing Research Tunnel that studied both supercooled liquid and ice-crystal icing. The test article that was utilized replicated 3D geometrical features of an inter-compressor duct and strut region of a turbofan engine. The surfaces of the Simulated Inter-compressor Duct Research Model (SIDRM) can be heated to simulate the warm surfaces of the turbofan inter-compressor duct. The test article is instrumented with pressure taps, heaters, heat flux gauges, and thermocouples, while a 3D laser scanner, cameras, and a weighing scale to measure ice mass were utilized to characterize the icing behavior. The aim of these tests was to generate ice accretions on the SIDRM test article under well-characterized icing conditions. To that end, ice crystal cloud characterization tests were conducted that provided cloud property values at the test section. The ice crystal accretion tests investigated how different test condition parameters impacted ice accretion size, location, and characteristics (such as physical attributes). A key finding was that a sufficiently high surface heat flux at the leading edge can provide a continuous supply of liquid melt to critical accretion zones downstream, resulting in elevated ice mass accretions. Another key finding was that warmer heater settings in critical accretion zones can suppress ice accretion. These findings are consistent with icing behavior observed in full-scale engine icing tests conducted previously at the NASA Glenn Research Center. Test runs investigating the contribution of liquid melt generated at the leading edge found that the melt can be transported downstream to critical accretion regions in two forms, via runback and splash. The form of melt supplied downstream affected accretion behavior and features. Sharkteeth-shaped ice accretions were generated at certain test conditions, a physical attribute that has been observed in full-scale engine icing tests. The ice accretion data will be used to develop and validate 3D computational engine icing tools, such as GlennICE, that predictively assesses the onset and growth of ice. One of the goals of the sponsoring NASA project is to develop simulation models and tools that can assist in the design and certification of engines for flight in icing conditions in a cost‑effective way.

ice crystal icing↗

Solar radiative transfer in cirrus clouds. I - Single-scattering and optical properties of hexagonal ice crystals. II - Theory and computations of multiple scattering in an anisotropic medium

The light scattering and absorption programs of Cai and Liou (1982) and Takano and Jaweera (1985) are extended to include hexagonal ice crystals randomly and horizontally oriented in space. The scattering and polarization results for the ice crystals are calculated. The results are compared with measurement data. The single-scattering properties for horizontally oriented columns and plates are presented and used to explain halos and arcs observed in the atmopshere. In the second section, the theory and computations for multiple scattering in cirrus clouds containing oriented ice crystals are presented. The radiative transfer in clouds composed of horizontally oriented ice crystals is formulated. Also, reflected and transmitted intensities, planetary albedo, and polarization in multiple scattering by ice crystals are discussed.

Takano, Yoshihide↗

Aerodynamic characteristics of airfoils with ice accretions

Results of a wind tunnel test to evaluate the performance of an airfoil with simulated rime ice are presented with theoretical comparisons. A NACA 65A413 airfoil was tested in the OSU 6 x 22 inch Transonic Airfoil Wind Tunnel at a Reynolds number near three million and Mach numbers from 0.20 to 0.80. The model was tested in four configurations to determine the aero-dynamic effects of the roughness and shape of a rime ice accretion. The simulated rime ice shape was obtained analytically using a time-stepping dry ice accretion computer code. Lift, drag, moment coefficients, and pressure distributions for the clean and simulated rime ice cases are reported. The measured degradation in airfoil performance is compared to an analytical method which uses existing airfoil analysis computer codes with empirical corrections for the surface roughness. A discussion of the empirical surface roughness correction and uses of other airfoil computer methods is included.

Bragg, M. B.↗

A hybrid deep neural operator/finite element method for ice-sheet modeling

One of the most challenging and consequential problems in climate modeling is to provide probabilistic projections of sea level rise. A large part of the uncertainty of sea level projections is due to uncertainty in ice sheet dynamics. At the moment, accurate quantification of the uncertainty is hindered by the cost of ice sheet computational models. In this work we develop a hybrid approach to approximate existing ice sheet models at a fraction of their cost. Our approach consists of replacing the finite element model for the momentum equations for the ice velocity, the most expensive part of an ice sheet model, with a Deep Operator Network, while we retain a classic finite element discretization for the evolution of the ice thickness. We show that the resulting hybrid model is very accurate and it is an order of magnitude faster than the traditional finite element model. Further, a distinctive feature of the proposed model, compared to other neural network approaches, is that it can handle high-dimensional parameter spaces (parameter fields) such as the basal friction at the bed of the glacier and can therefore be used for generating samples for uncertainty quantification. Further, we study the impact of hyper-parameters, number of unknowns and correlation length of the parameter distribution on the training and accuracy of the Deep Operator Network on a synthetic ice sheet model. We then target the evolution of the Humboldt glacier in Greenland and show that our hybrid model can provide accurate statistics of the glacier mass loss and can be effectively used to accelerate the quantification of uncertainty.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗