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 19 records

Three Dimensional Surface Redefinition Method for Computational Ice Accretion Solvers

Computational tools have been increasing in maturity and are thus commonly used in the engineering design process. Advancements in NASA’s current state of the art computational ice accretion tool are required to tackle the icing challenges of tomorrow. GlennICE, a next generation ice accretion solver, is under development at NASA to tackle these challenges. One of the elements of this migration is transitioning from a quasi-three dimensional strip theory based ice accretion methodology to a fully three dimensional methodology. This requires construction of a method to redefine or extrude a discretized or tessellated surface geometry based on the predicted volumetric ice growth for each tessellated surface triangle. This paper describes the methodology that is employed in the GlennICE software, and assesses the performance of the method in replicating a well definied analytical test case.

Computational↗

Three Dimensional Surface Redefinition Method for Computational Ice Accretion Solvers

Computational tools have been increasing in maturity and are thus commonly used in the engineering design process. Advancements in NASA’s current state of the art computational ice accretion tool are required to tackle the icing challenges of tomorrow. GlennICE, a next generation ice accretion solver, is under development at NASA to tackle these challenges. One of the elements of this migration is transitioning from a quasi-three dimensional strip theory based ice accretion methodology to a fully three dimensional methodology. This requires construction of a method to redefine or extrude a discretized or tessellated surface geometry based on the predicted volumetric ice growth for each tessellated surface triangle. This paper describes the methodology that is employed in the GlennICE software, and assesses the performance of the method in replicating a well defined analytical test case.

Computational↗

Computational Icing Analysis on NASA’s SIDRM Geometry to Investigate Collection Efficiency

Computational icing analysis results were compared to experimental icing tunnel data including aerothermal (e.g. dry air) and supercooled water droplet rime-ice conditions from tests conducted in early 2022 at the NASA Icing Research Tunnel (IRT). The Simulated Intercompressor Duct Research Model (SIDRM) test article was used in this study, and its geometry represents the inter-compressor duct region of a turbofan engine. The test article’s purpose is to study the physics of supercooled water icing and ice crystal icing. This study compared three different icing codes: FENSAP-ICE (Eulerian approach), LEWICE3D (Lagrangian approach), and GlennICE (Lagrangian approach). All three icing codes were conducted on SIDRM’s complex body flow-field and compared to different experimental supercooled water rime runs. The test article instrumentation (pressure taps, thermocouples, etc.) and 3D laser scans of final ice shapes were used to compare against the different icing code simulations. The overall objectives are to understand how the icing codes handle capturing collection efficiency on the complex test article’s unheated surfaces. In the aerothermal cases, pressure tap readings matched the CFD results, but dry air CFD underpredicted thermocouple readings. Collection efficiency results from all three icing codes matched well together on the main body leading edge, main body slope, and the strut leading edges. All three icing codes underpredicted the final ice shape using a single-shot constant ice density approach, with more difficulty coming from the strut leading edge ice shape due to the swept wing like flow field. NASA’s overall goal for this effort is to develop computational icing tools to assist in the design and certification of engines for flight in icing conditions.

Supercooled icing, SIDRM, aerodynamics, computatio↗

Software Development Processes Applied to Computational Icing Simulation

The development of computational icing simulation methods is making the transition form the research to common place use in design and certification efforts. As such, standards of code management, design validation, and documentation must be adjusted to accommodate the increased expectations of the user community with respect to accuracy, reliability, capability, and usability. This paper discusses these concepts with regard to current and future icing simulation code development efforts as implemented by the Icing Branch of the NASA Lewis Research Center in collaboration with the NASA Lewis Engineering Design and Analysis Division. With the application of the techniques outlined in this paper, the LEWICE ice accretion code has become a more stable and reliable software product.

Levinson, Laurie H.↗

Introduction to the Glenn Icing Computational Environment (GlennICE)

The NASA John H. Glenn Research Center at Lewis Field is developing the Glenn Icing Computational Environment (GlennICE) tool to aid those evaluating, designing and certifying aircraft, engines, and aircraft components for flight in icing conditions. This short course will walk through some of the underlying physics involved with GlennICE and how we achieve efficient 3D ice accretion predictions. After an introduction of GlennICE, an analysis of the Common Research Model High-Lift will be demonstrated to showcase the typical workflow for an aircraft icing analysis. Within this walkthrough, capabilities will be highlighted with future planned capabilities being discussed. Finally, we will showcase the impact GlennICE is having on NASA’s icing portfolio and how it is advancing aircraft icing research and safety.

Icing↗

Computational Icing Risk Analysis of the D8 "Double Bubble" Aircraft

A computational icing risk analysis utilizing LEWICE3D was performed for the D8 Double Bubble aircraft. A variety of discrete drop sizes spanning the Appendix C and O regimes were simulated. For computational efficiency a 50-bin global discretization was produced and projected onto the distributions of interest, eliminating redundant simulations. The trajectory and impingement characteristics for discrete drop diameters were analyzed to help understand the behavior of the water drops in the presence of a complex flow field. The collection efficiency results for the discrete drop diameters were then weighted by their contributions to the total water content of six different continuous distributions and subsequently superposed to approximate these curves. Results indicate that significant variation in impingement exists as a function of drop diameter for complex wing body geometries, and that current discretization practices may be insufficient to accurately predict water collection on certain regions of the aircraft. Results also indicate that the Appendix O distributions, specifically those with considerable water content at large drops, generates water collection patterns that are markedly different from distributions representative of Appendix C.

icing↗

Computational Icing Risk Analysis of the D8 "Double Bubble" Aircraft

A computational icing risk analysis utilizing LEWICE3D was performed for the D8 Double Bubble aircraft. A variety of discrete drop sizes spanning the Appendix C and O regimes were simulated. For computational efficiency a 50-bin global discretization was produced and projected onto the distributions of interest, eliminating redundant simulations. The trajectory and impingement characteristics for discrete drop diameters were analyzed to help understand the behavior of the water drops in the presence of a complex flow field. The collection efficiency results for the discrete drop diameters were then weighted by their contributions to the total water content of six different continuous distributions and subsequently superposed to approximate these curves. Results indicate that significant variation in impingement exists as a function of drop diameter for complex wing body geometries, and that current discretization practices may be insufficient to accurately predict water collection on certain regions of the aircraft. Results also indicate that the Appendix O distributions, specifically those with considerable water content at large drops, generates water collection patterns that are markedly different from distributions representative of Appendix C.

Porter, Christopher E.↗

Experimental and Computational Icing Simulation for Large Swept Wings

The purpose of this final report is to summarize the key results, findings and contributions of a large, multi-organization, multi-year effort focused on swept wing icing and aerodynamics. This research effort was jointly sponsored by NASA, FAA and ONERA and was supported by Boeing and the Universities of Illinois, Virginia, and Washington. The overall goal of this research was to improve the experimental and computational simulation capability for icing on large swept wings typical of commercial transports. This research included both ice accretion and aerodynamic studies using the NASA Common Research Model (CRM) as the reference geometry. For this work, a 65% scaled version—CRM65—was used as the full-scale baseline airplane geometry. This program marks the first non-proprietary research to generate and document full-scale swept wing ice accretions in an icing tunnel and then conduct semispan wing aerodynamic testing on scaled ice accretions at high Reynolds number (i.e., up to 11.9×106). During this work, several new experimental and computational techniques were developed or documented and new knowledge about swept wing icing aerodynamics was identified. All are documented in this report.

Aircraft icing↗

A Computational Icing Effects Study for a Three-Dimensional Wing: Comparison with Experimental Data and Investigation of Spanwise Variation

Steady-state solutions of the Reynolds-averaged Navier-Stokes (RANS) equations were computed using the Colbalt flow solver for a constant-section, rectangular wing based on an extruded two-dimensional glaze ice shape. The one equation Spalart-Allmaras turbulence model was used. The results were compared with data obtained from a recent wind tunnel test. Computed results indicate that the steady RANS solutions do not accurately capture the recirculating region downstream of the ice accretion, even after a mesh refinement. The resulting predicted reattachment is farther downstream than indicated by the experimental data. Additionally, the solutions computed on a relatively coarse baseline mesh had detailed flow characteristics that were different from those computed on the refined mesh or the experimental data. Steady RANS solutions were also computed to investigate the effects of spanwise variation in the ice shape. The spanwise variation was obtained via a bleeding function that merged the ice shape with the clean wing using a sinusoidal spanwise variation. For these configurations, the results predicted for the extruded shape provided conservative estimates for the performance degradation of the wing. Additionally, the spanwise variation in the ice shape and the resulting differences in the flow fields did not significantly change the location of the primary reattachment.

Thompson, D.↗

Experimental and computational ice shapes and resulting drag increase for a NACA 0012 airfoil

Tests were conducted in the Icing Research Tunnel (IRT) at LeRC to document the repeatability of the ice shape over the range of temperatures varying from -15 to 28 F. Measurements of drag increase due to the ice accretion were also made. The ice shape and drag coefficient data, with varying total temperatures at two different airspeeds, were compared with the computational predictions. The calculations were made with the 2D LEWICE/IBL code which is a combined code of LEWICE and the interactive boundary layer method developed for iced airfoils. Comparisons show good agreement with the experimental data in ice shapes. The calculations show the ability of the code to predict drag increases as the ice shape changes from a rime shape to a glaze shape.

Shin, Jaiwon↗

Experimental and computational ice shapes and resulting drag increase for a NACA 0012 airfoil

Tests were conducted in the Icing Research Tunnel (IRT) at the NASA Lewis Research Center to document the repeatability of the ice shape over the range of temperatures varying from -15 F to 28 F. Measurements of drag increase due to the ice accretion were also made. The ice shape and drag coefficient data, with varying total temperatures at two different airspeeds, were compared with the computational predictions. The calculations were made with the 2D LEWICE/IBL code which is a combined code of LEWICE and the interactive boundary layer method developed for iced airfoils. Comparisons show good agreement with the experimental data in ice shapes. The calculations show the ability of the code to predict drag increases as the ice shape changes from a rime shape to a glaze shape.

Shin, Jaiwon↗

Computational Simulation of the Formation and Material Behavior of Ice

Computational methods are described for simulating the formation and the material behavior of ice in prevailing transient environments. The methodology developed at the NASA Lewis Research Center was adopted. A three dimensional finite-element heat transfer analyzer was used to predict the thickness of ice formed under prevailing environmental conditions. A multi-factor interaction model for simulating the material behavior of time-variant ice layers is presented. The model, used in conjunction with laminated composite mechanics, updates the material properties of an ice block as its thickness increases with time. A sample case of ice formation in a body of water was used to demonstrate the methodology. The results showed that the formation and the material behavior of ice can be computationally simulated using the available composites technology.

Tong, Michael T.↗

Validation of NASA Thermal Ice Protection Computer Codes: Program Overview - Part 1

The Icing Technology Branch at NASA Lewis has been involved in an effort to validate two thermal ice protection codes developed at the NASA Lewis Research Center. LEWICE/Thermal (electrothermal deicing & anti-icing), and ANTICE (hot-gas & electrothermal anti-icing). The Thermal Code Validation effort was designated as a priority during a 1994 'peer review' of the NASA Lewis Icing program, and was implemented as a cooperative effort with industry. During April 1996, the first of a series of experimental validation tests was conducted in the NASA Lewis Icing Research Tunnel(IRT). The purpose of the April 96 test was to validate the electrothermal predictive capabilities of both LEWICE/Thermal, and ANTICE. A heavily instrumented test article was designed and fabricated for this test, with the capability of simulating electrothermal de-icing and anti-icing modes of operation. Thermal measurements were then obtained over a range of test conditions, for comparison with analytical predictions. This paper will present an overview of the test, including a detailed description of: (1) the validation process; (2) test article design; (3) test matrix development; and (4) test procedures. Selected experimental results will be presented for de-icing and anti-icing modes of operation. Finally, the status of the validation effort at this point will be summarized. Detailed comparisons between analytical predictions and experimental results are contained in the following two papers: 'Validation of NASA Thermal Ice Protection Computer Codes: Part 2- The Validation of LEWICE/Thermal' and 'Validation of NASA Thermal Ice Protection Computer Codes: Part 3-The Validation of ANTICE'

Miller, Dean↗

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗

Probing the structure of cometary ice

Computer simulations of bulk and vapor deposited amorphous ices are presented. The structure of the bulk low density amorphous ice is in good agreement with experiments on pressure disordered amorphous ice. Both the low density bulk ice and the vapor deposited ices exhibit strong ordering. Vapor deposition of hot (300 K) water molecules onto a cold (77 K) substrate yields less porous ices than deposition of cold (77 K) water molecules onto a cold substrate. Both vapor deposited ices are more porous than the bulk amorphous ice. The structure of bulk high density amorphous ice is only in fair agreement with experimental results. Attempts to simulate high density amorphous ice via vapor deposition were not successful. Electron diffraction results on vapor deposited amorphous ice indicate that the temperature of the nucleation of the cubic phase depends upon the amount of time between the deposition and the onset of crystallization, suggesting that freshly deposited ice layers reconstruct on time of the order of hours. The temperature dependence of the microporosity of the vapor deposited amorphous ices might affect laboratory experiments that are aimed at simulating astrophysical ices in the context of the origin of prebiotic organic material and its transport to the Earth.

Wilson, Michael A.↗

Analysis of iced wings

A method for computing ice shapes along the leading edge of a wing and a method for predicting its aerodynamic performance degradation due to icing is described. Ice shapes are computed using an extension of the LEWICE code which was developed for airfoils. The aerodynamic properties of the iced wing are determined with an interactive scheme in which the solutions of the inviscid flow equations are obtained from a panel method and the solutions of the viscous flow equations are obtained from an inverse three-dimensional finite-difference boundary-layer method. A new interaction law is used to couple the inviscid and viscous flow solutions. The application of the LEWICE wing code to the calculation of ice shapes on a MS-317 swept wing shows good agreement with measurements. The interactive boundary-layer method is applied to a tapered ice wing in order to study the effect of icing on the aerodynamic properties of the wing at several angles of attack.

Cebeci, T.↗