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

NASA Technical Reports Server (NTRS)

Jorgenson, Philip C. E.; Veres, Joseph P.; Wright, William B.; Struk, Peter M.

2013-01-01

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

Ice Accretion Formations on a NACA 0012 Swept Wing Tip in Natural Icing Conditions

NASA Technical Reports Server (NTRS)

Vargas, Mario; Giriunas, Julius A.; Ratvasky, Thomas P.

2002-01-01

An experiment was conducted in the DeHavilland DHC-6 Twin Otter Icing Research Aircraft at NASA Glenn Research Center to study the formation of ice accretions on swept wings in natural icing conditions. The experiment was designed to obtain ice accretion data to help determine if the mechanisms of ice accretion formation observed in the Icing Research Tunnel are present in natural icing conditions. The experiment in the Twin Otter was conducted using a NACA 0012 swept wing tip. The model enabled data acquisition at 0 deg, 15 deg, 25 deg, 30 deg, and 45 deg sweep angles. Casting data, ice shape tracings, and close-up photographic data were obtained. The results showed that the mechanisms of ice accretion formation observed in-flight agree well with the ones observed in the Icing Research Tunnel. Observations on the end cap of the airfoil showed the same strong effect of the local sweep angle on the formation of scallops as observed in the tunnel.

Ice Accretion Test Results for Three Large-Scale Swept-Wing Models in the NASA Icing Research Tunnel

NASA Technical Reports Server (NTRS)

Broeren, Andy; Potapczuk, Mark; Lee, Sam; Malone, Adam; Paul, Ben; Woodard, Brian

2016-01-01

The design and certification of modern transport airplanes for flight in icing conditions increasing relies on three-dimensional numerical simulation tools for ice accretion prediction. There is currently no publically available, high-quality, ice accretion database upon which to evaluate the performance of icing simulation tools for large-scale swept wings that are representative of modern commercial transport airplanes. The purpose of this presentation is to present the results of a series of icing wind tunnel test campaigns whose aim was to provide an ice accretion database for large-scale, swept wings.

A Database of Supercooled Large Droplet Ice Accretions

NASA Technical Reports Server (NTRS)

VanZante, Judith Foss

2007-01-01

A unique, publicly available database regarding supercooled large droplet ice accretions has been developed in NASA Glenn's Icing Research Tunnel. Identical cloud and flight conditions were generated for five different airfoil models. The models chosen represent a variety of aircraft types from the horizontal stabilizer of a large trans-port aircraft to the wings of regional, business, and general aviation aircraft. In addition to the standard documentation methods of 2D ice shape tracing and imagery, ice mass measurements were also taken. This database will also be used to validate and verify the extension of the ice accretion code, LEWICE, into the SLD realm.

Ice Accretions and Full-Scale Iced Aerodynamic Performance Data for a Two-Dimensional NACA 23012 Airfoil

NASA Technical Reports Server (NTRS)

Addy, Harold E., Jr.; Broeren, Andy P.; Potapczuk, Mark G.; Lee, Sam; Guffond, Didier; Montreuil, Emmanuel; Moens, Frederic

2016-01-01

This report documents the data collected during the large wind tunnel campaigns conducted as part of the SUNSET project (StUdies oN Scaling EffecTs due to ice) also known as the Ice-Accretion Aerodynamics Simulation study: a joint effort by NASA, the Office National d'Etudes et Recherches Aérospatiales (ONERA), and the University of Illinois. These data form a benchmark database of full-scale ice accretions and corresponding ice-contaminated aerodynamic performance data for a two-dimensional (2D) NACA 23012 airfoil. The wider research effort also included an analysis of ice-contaminated aerodynamics that categorized ice accretions by aerodynamic effects and an investigation of subscale, low- Reynolds-number ice-contaminated aerodynamics for the NACA 23012 airfoil. The low-Reynolds-number investigation included an analysis of the geometric fidelity needed to reliably assess aerodynamic effects of airfoil icing using artificial ice shapes. Included herein are records of the ice accreted during campaigns in NASA Glenn Research Center's Icing Research Tunnel (IRT). Two different 2D NACA 23012 airfoil models were used during these campaigns; an 18-in. (45.7-cm) chord (subscale) model and a 72-in. (182.9-cm) chord (full-scale) model. The aircraft icing conditions used during these campaigns were selected from the Federal Aviation Administration's (FAA's) Code of Federal Regulations (CFR) Part 25 Appendix C icing envelopes. The records include the test conditions, photographs of the ice accreted, tracings of the ice, and ice depth measurements. Model coordinates and pressure tap locations are also presented. Also included herein are the data recorded during a wind tunnel campaign conducted in the F1 Subsonic Pressurized Wind Tunnel of ONERA. The F1 tunnel is a pressured, high- Reynolds-number facility that could accommodate the full-scale (72-in. (182.9-cm) chord) 2D NACA 23012 model. Molds were made of the ice accreted during selected test runs of the full-scale model

Helicopter rotor noise investigation during ice accretion

NASA Astrophysics Data System (ADS)

Cheng, Baofeng

An investigation of helicopter rotor noise during ice accretion is conducted using experimental, theoretical, and numerical methods. This research is the acoustic part of a joint helicopter rotor icing physics, modeling, and detection project at The Pennsylvania State University Vertical Lift Research Center of Excellence (VLRCOE). The current research aims to provide acoustic insight and understanding of the rotor icing physics and investigate the feasibility of detecting rotor icing through noise measurements, especially at the early stage of ice accretion. All helicopter main rotor noise source mechanisms and their change during ice accretion are discussed. Changes of the thickness noise, steady loading noise, and especially the turbulent boundary layer - trailing edge (TBL-TE) noise due to ice accretion are identified and studied. The change of the discrete frequency noise (thickness noise and steady loading noise) due to ice accretion is calculated by using PSU-WOPWOP, an advanced rotorcraft acoustic prediction code. The change is noticeable, but too small to be used in icing detection. The small thickness noise change is due to the small volume of the accreted ice compared to that of the entire blade, although a large iced airfoil shape is used. For the loading noise calculation, two simplified methods are used to generate the loading on the rotor blades, which is the input for the loading noise calculation: 1) compact loading from blade element momentum theory, icing effects are considered by increasing the drag coefficient; and 2) pressure loading from the 2-D CFD simulation, icing effects are considered by using the iced airfoil shape. Comprehensive rotor broadband noise measurements are carried out on rotor blades with different roughness sizes and rotation speeds in two facilities: the Adverse Environment Rotor Test Stand (AERTS) facility at The Pennsylvania State University, and The University of Maryland Acoustic Chamber (UMAC). In both facilities the

Development of 3D Ice Accretion Measurement Method

NASA Technical Reports Server (NTRS)

Lee, Sam; Broeren, Andy P.; Addy, Harold E., Jr.; Sills, Robert; Pifer, Ellen M.

2012-01-01

Icing wind tunnels are designed to simulate in-flight icing environments. The chief product of such facilities is the ice accretion that forms on various test articles. Documentation of the resulting ice accretion key piece of data in icing-wind-tunnel tests. Number of currently used options for documenting ice accretion in icing-wind-tunnel testing.

Implementation and Validation of 3-D Ice Accretion Measurement Methodology

NASA Technical Reports Server (NTRS)

Lee, Sam; Broeren, Andy P.; Kreeger, Richard E.; Potapczuk, Mark; Utt, Lloyd

2014-01-01

A research program has been implemented to develop and validate the use of a commercial 3-D laser scanning system to record ice accretion geometry in the NASA Icing Research Tunnel. A main component of the program was the geometric assessment of the 3- D laser scanning system on a 2-D (straight wing) and a 3-D (swept wing) airfoil geometries. This exercise consisted of comparison of scanned ice accretion to castings of the same ice accretion. The scan data were also used to create rapid prototype artificial ice shapes that were scanned and compared to the original ice accretion. The results from geometric comparisons on the straight wing showed that the ice shape models generated through the scan/rapid prototype process compared reasonably well with the cast shapes. Similar results were obtained with the geometric comparisons on the swept wing. It was difficult to precisely compare the scans of the cast shapes to the original ice accretion scans because the cast shapes appear to have shrunk during the mold/casting process by as much as 0.10-inch. However the comparison of the local ice-shape features were possible and produced better results. The rapid prototype manufacturing process was shown to reproduce the original ice accretion scan normally within 0.01-inch.

Investigation of surface water behavior during glaze ice accretion

NASA Technical Reports Server (NTRS)

Hansman, R. John, Jr.; Turnock, Stephen R.

1990-01-01

A series of experimental investigations that focused on isolating the primary factors that control the behavior of unfrozen surface water during glaze ice accretion were conducted. Detailed microvideo observations were made of glaze ice accretions on 2.54 cm diam cylinders in a closed-loop refrigerated wind tunnel. Distinct zones of surface water behavior were observed; a smooth wet zone in the stagnation region with a uniform water film, a rough zone where surface tension effects caused coalescence of surface water into stationary beads, and a zone where surface water ran back as rivulets. The location of the transition from the smooth to the rough zone was found to migrate towards the stagnation point with time. Comparative tests were conducted to study the effect of the substrate thermal and roughness properties on ice accretion. The importance of surface water behavior was evaluated by the addition of a surface tension reducing agent to the icing tunnel water supply, which significantly altered the accreted glaze ice shape. Measurements were made to determine the contact angle behavior of water droplets on ice. A simple multizone modification to current glaze ice accretion models was proposed to include the observed surface roughness behavior.

Icing Simulation Research Supporting the Ice-Accretion Testing of Large-Scale Swept-Wing Models

NASA Technical Reports Server (NTRS)

Yadlin, Yoram; Monnig, Jaime T.; Malone, Adam M.; Paul, Bernard P.

2018-01-01

The work summarized in this report is a continuation of NASA's Large-Scale, Swept-Wing Test Articles Fabrication; Research and Test Support for NASA IRT contract (NNC10BA05 -NNC14TA36T) performed by Boeing under the NASA Research and Technology for Aerospace Propulsion Systems (RTAPS) contract. In the study conducted under RTAPS, a series of icing tests in the Icing Research Tunnel (IRT) have been conducted to characterize ice formations on large-scale swept wings representative of modern commercial transport airplanes. The outcome of that campaign was a large database of ice-accretion geometries that can be used for subsequent aerodynamic evaluation in other experimental facilities and for validation of ice-accretion prediction codes.

Airfoil Ice-Accretion Aerodynamics Simulation

NASA Technical Reports Server (NTRS)

Bragg, Michael B.; Broeren, Andy P.; Addy, Harold E.; Potapczuk, Mark G.; Guffond, Didier; Montreuil, E.

2007-01-01

NASA Glenn Research Center, ONERA, and the University of Illinois are conducting a major research program whose goal is to improve our understanding of the aerodynamic scaling of ice accretions on airfoils. The program when it is completed will result in validated scaled simulation methods that produce the essential aerodynamic features of the full-scale iced-airfoil. This research will provide some of the first, high-fidelity, full-scale, iced-airfoil aerodynamic data. An initial study classified ice accretions based on their aerodynamics into four types: roughness, streamwise ice, horn ice, and spanwise-ridge ice. Subscale testing using a NACA 23012 airfoil was performed in the NASA IRT and University of Illinois wind tunnel to better understand the aerodynamics of these ice types and to test various levels of ice simulation fidelity. These studies are briefly reviewed here and have been presented in more detail in other papers. Based on these results, full-scale testing at the ONERA F1 tunnel using cast ice shapes obtained from molds taken in the IRT will provide full-scale iced airfoil data from full-scale ice accretions. Using these data as a baseline, the final step is to validate the simulation methods in scale in the Illinois wind tunnel. Computational ice accretion methods including LEWICE and ONICE have been used to guide the experiments and are briefly described and results shown. When full-scale and simulation aerodynamic results are available, these data will be used to further develop computational tools. Thus the purpose of the paper is to present an overview of the program and key results to date.

Ice Accretions and Icing Effects for Modern Airfoils

NASA Technical Reports Server (NTRS)

Addy, Harold E., Jr.

2000-01-01

Icing tests were conducted to document ice shapes formed on three different two-dimensional airfoils and to study the effects of the accreted ice on aerodynamic performance. The models tested were representative of airfoil designs in current use for each of the commercial transport, business jet, and general aviation categories of aircraft. The models were subjected to a range of icing conditions in an icing wind tunnel. The conditions were selected primarily from the Federal Aviation Administration's Federal Aviation Regulations 25 Appendix C atmospheric icing conditions. A few large droplet icing conditions were included. To verify the aerodynamic performance measurements, molds were made of selected ice shapes formed in the icing tunnel. Castings of the ice were made from the molds and placed on a model in a dry, low-turbulence wind tunnel where precision aerodynamic performance measurements were made. Documentation of all the ice shapes and the aerodynamic performance measurements made during the icing tunnel tests is included in this report. Results from the dry, low-turbulence wind tunnel tests are also presented.

A review of ice accretion data from a model rotor icing test and comparison with theory

NASA Technical Reports Server (NTRS)

Britton, Randall K.; Bond, Thomas H.

1991-01-01

An experiment was conducted by the Helicopter Icing Consortium (HIC) in the NASA Lewis Icing Research Tunnel (IRT) in which a 1/6 scale fuselage model of a UH-60A Black Hawk helicopter with a generic rotor was subjected to a wide range of icing conditions. The HIC consists of members from NASA, Bell Helicopter, Boeing Helicopter, McDonnell Douglas Helicopters, Sikorsky Aircraft, and Texas A&M University. Data was taken in the form of rotor torque, internal force balance measurements, blade strain gage loading, and two dimensional ice shape tracings. A review of the ice shape data is performed with special attention given to repeatability and correctness of trends in terms of radial variation, rotational speed, icing time, temperature, liquid water content, and volumetric median droplet size. Moreover, an indepth comparison between the experimental data and the analysis of NASA's ice accretion code LEWICE is given. Finally, conclusions are drawn as to the quality of the ice accretion data and the predictability of the data base as a whole. Recommendations are also given for improving data taking technique as well as potential future work.

A Database of Supercooled Large Droplet Ice Accretions [Supplement

NASA Technical Reports Server (NTRS)

VanZante, Judith Foss

2007-01-01

A unique, publicly available database regarding supercooled large droplet (SLD) ice accretions has been developed in NASA Glenn's Icing Research Tunnel. Identical cloud and flight conditions were generated for five different airfoil models. The models chosen represent a variety of aircraft types from the horizontal stabilizer of a large transport aircraft to the wings of regional, business, and general aviation aircraft. In addition to the standard documentation methods of 2D ice shape tracing and imagery, ice mass measurements were also taken. This database will also be used to validate and verify the extension of the ice accretion code, LEWICE, into the SLD realm.

Linking ice accretion and crown structure: towards a model of the effect of freezing rain on tree canopies.

PubMed

Nock, Charles A; Lecigne, Bastien; Taugourdeau, Olivier; Greene, David F; Dauzat, Jean; Delagrange, Sylvain; Messier, Christian

2016-06-01

Despite a longstanding interest in variation in tree species vulnerability to ice storm damage, quantitative analyses of the influence of crown structure on within-crown variation in ice accretion are rare. In particular, the effect of prior interception by higher branches on lower branch accumulation remains unstudied. The aim of this study was to test the hypothesis that intra-crown ice accretion can be predicted by a measure of the degree of sheltering by neighbouring branches. Freezing rain was artificially applied to Acer platanoides L., and in situ branch-ice thickness was measured directly and from LiDAR point clouds. Two models of freezing rain interception were developed: 'IceCube', which uses point clouds to relate ice accretion to a voxel-based index (sheltering factor; SF) of the sheltering effect of branch elements above a measurement point; and 'IceTree', a simulation model for in silico evaluation of the interception pattern of freezing rain in virtual tree crowns. Intra-crown radial ice accretion varied strongly, declining from the tips to the bases of branches and from the top to the base of the crown. SF for branches varied strongly within the crown, and differences among branches were consistent for a range of model parameters. Intra-crown variation in ice accretion on branches was related to SF (R(2) = 0·46), with in silico results from IceTree supporting empirical relationships from IceCube. Empirical results and simulations confirmed a key role for crown architecture in determining intra-crown patterns of ice accretion. As suspected, the concentration of freezing rain droplets is attenuated by passage through the upper crown, and thus higher branches accumulate more ice than lower branches. This is the first step in developing a model that can provide a quantitative basis for investigating intra-crown and inter-specific variation in freezing rain damage. © The Author 2016. Published by Oxford University Press on behalf of the Annals of

Linking ice accretion and crown structure: towards a model of the effect of freezing rain on tree canopies

PubMed Central

Nock, Charles A.; Lecigne, Bastien; Taugourdeau, Olivier; Greene, David F.; Dauzat, Jean; Delagrange, Sylvain; Messier, Christian

2016-01-01

Background and Aims Despite a longstanding interest in variation in tree species vulnerability to ice storm damage, quantitative analyses of the influence of crown structure on within-crown variation in ice accretion are rare. In particular, the effect of prior interception by higher branches on lower branch accumulation remains unstudied. The aim of this study was to test the hypothesis that intra-crown ice accretion can be predicted by a measure of the degree of sheltering by neighbouring branches. Methods Freezing rain was artificially applied to Acer platanoides L., and in situ branch-ice thickness was measured directly and from LiDAR point clouds. Two models of freezing rain interception were developed: ‘IceCube’, which uses point clouds to relate ice accretion to a voxel-based index (sheltering factor; SF) of the sheltering effect of branch elements above a measurement point; and ‘IceTree’, a simulation model for in silico evaluation of the interception pattern of freezing rain in virtual tree crowns. Key Results Intra-crown radial ice accretion varied strongly, declining from the tips to the bases of branches and from the top to the base of the crown. SF for branches varied strongly within the crown, and differences among branches were consistent for a range of model parameters. Intra-crown variation in ice accretion on branches was related to SF (R2 = 0·46), with in silico results from IceTree supporting empirical relationships from IceCube. Conclusions Empirical results and simulations confirmed a key role for crown architecture in determining intra-crown patterns of ice accretion. As suspected, the concentration of freezing rain droplets is attenuated by passage through the upper crown, and thus higher branches accumulate more ice than lower branches. This is the first step in developing a model that can provide a quantitative basis for investigating intra-crown and inter-specific variation in freezing rain damage. PMID:27107412

DRA/NASA/ONERA Collaboration on Icing Research. Part 2; Prediction of Airfoil Ice Accretion

NASA Technical Reports Server (NTRS)

Wright, William B.; Gent, R. W.; Guffond, Didier

1997-01-01

This report presents results from a joint study by DRA, NASA, and ONERA for the purpose of comparing, improving, and validating the aircraft icing computer codes developed by each agency. These codes are of three kinds: (1) water droplet trajectory prediction, (2) ice accretion modeling, and (3) transient electrothermal deicer analysis. In this joint study, the agencies compared their code predictions with each other and with experimental results. These comparison exercises were published in three technical reports, each with joint authorship. DRA published and had first authorship of Part 1 - Droplet Trajectory Calculations, NASA of Part 2 - Ice Accretion Prediction, and ONERA of Part 3 - Electrothermal Deicer Analysis. The results cover work done during the period from August 1986 to late 1991. As a result, all of the information in this report is dated. Where necessary, current information is provided to show the direction of current research. In this present report on ice accretion, each agency predicted ice shapes on two dimensional airfoils under icing conditions for which experimental ice shapes were available. In general, all three codes did a reasonable job of predicting the measured ice shapes. For any given experimental condition, one of the three codes predicted the general ice features (i.e., shape, impingement limits, mass of ice) somewhat better than did the other two. However, no single code consistently did better than the other two over the full range of conditions examined, which included rime, mixed, and glaze ice conditions. In several of the cases, DRA showed that the user's knowledge of icing can significantly improve the accuracy of the code prediction. Rime ice predictions were reasonably accurate and consistent among the codes, because droplets freeze on impact and the freezing model is simple. Glaze ice predictions were less accurate and less consistent among the codes, because the freezing model is more complex and is critically

Aerodynamic Simulation of Runback Ice Accretion

NASA Technical Reports Server (NTRS)

Broeren, Andy P.; Whalen, Edward A.; Busch, Greg T.; Bragg, Michael B.

2010-01-01

This report presents the results of recent investigations into the aerodynamics of simulated runback ice accretion on airfoils. Aerodynamic tests were performed on a full-scale model using a high-fidelity, ice-casting simulation at near-flight Reynolds (Re) number. The ice-casting simulation was attached to the leading edge of a 72-in. (1828.8-mm ) chord NACA 23012 airfoil model. Aerodynamic performance tests were conducted at the ONERA F1 pressurized wind tunnel over a Reynolds number range of 4.7?10(exp 6) to 16.0?10(exp 6) and a Mach (M) number ran ge of 0.10 to 0.28. For Re = 16.0?10(exp 6) and M = 0.20, the simulated runback ice accretion on the airfoil decreased the maximum lift coe fficient from 1.82 to 1.51 and decreased the stalling angle of attack from 18.1deg to 15.0deg. The pitching-moment slope was also increased and the drag coefficient was increased by more than a factor of two. In general, the performance effects were insensitive to Reynolds numb er and Mach number changes over the range tested. Follow-on, subscale aerodynamic tests were conducted on a quarter-scale NACA 23012 model (18-in. (457.2-mm) chord) at Re = 1.8?10(exp 6) and M = 0.18, using low-fidelity, geometrically scaled simulations of the full-scale castin g. It was found that simple, two-dimensional simulations of the upper- and lower-surface runback ridges provided the best representation of the full-scale, high Reynolds number iced-airfoil aerodynamics, whereas higher-fidelity simulations resulted in larger performance degrada tions. The experimental results were used to define a new subclassification of spanwise ridge ice that distinguishes between short and tall ridges. This subclassification is based upon the flow field and resulting aerodynamic characteristics, regardless of the physical size of the ridge and the ice-accretion mechanism.

A Systems-Level Perspective on Engine Ice Accretion

NASA Technical Reports Server (NTRS)

May, Ryan D.; Guo, Ten-Huei; Simon, Donald L.

2013-01-01

The accretion of ice in the compression system of commercial gas turbine engines operating in high ice water content conditions is a safety issue being studied by the aviation sector. While most of the research focuses on the underlying physics of ice accretion and the meteorological conditions in which accretion can occur, a systems-level perspective on the topic lends itself to potential near-term operational improvements. This work focuses on developing an accurate and reliable algorithm for detecting the accretion of ice in the low pressure compressor of a generic 40,000 lbf thrust class engine. The algorithm uses only the two shaft speed sensors and works regardless of engine age, operating condition, and power level. In a 10,000-case Monte Carlo simulation, the detection approach was found to have excellent capability at determining ice accretion from sensor noise with detection occurring when ice blocks an average of 6.8% of the low pressure compressor area. Finally, an initial study highlights a potential mitigation strategy that uses the existing engine actuators to raise the temperature in the low pressure compressor in an effort to reduce the rate at which ice accretes.

Validation of 3-D Ice Accretion Measurement Methodology for Experimental Aerodynamic Simulation

NASA Technical Reports Server (NTRS)

Broeren, Andy P.; Addy, Harold E., Jr.; Lee, Sam; Monastero, Marianne C.

2015-01-01

Determining the adverse aerodynamic effects due to ice accretion often relies on dry-air wind-tunnel testing of artificial, or simulated, ice shapes. Recent developments in ice-accretion documentation methods have yielded a laser-scanning capability that can measure highly three-dimensional (3-D) features of ice accreted in icing wind tunnels. The objective of this paper was to evaluate the aerodynamic accuracy of ice-accretion simulations generated from laser-scan data. Ice-accretion tests were conducted in the NASA Icing Research Tunnel using an 18-in. chord, two-dimensional (2-D) straight wing with NACA 23012 airfoil section. For six ice-accretion cases, a 3-D laser scan was performed to document the ice geometry prior to the molding process. Aerodynamic performance testing was conducted at the University of Illinois low-speed wind tunnel at a Reynolds number of 1.8 × 10(exp 6) and a Mach number of 0.18 with an 18-in. chord NACA 23012 airfoil model that was designed to accommodate the artificial ice shapes. The ice-accretion molds were used to fabricate one set of artificial ice shapes from polyurethane castings. The laser-scan data were used to fabricate another set of artificial ice shapes using rapid prototype manufacturing such as stereolithography. The iced-airfoil results with both sets of artificial ice shapes were compared to evaluate the aerodynamic simulation accuracy of the laser-scan data. For five of the six ice-accretion cases, there was excellent agreement in the iced-airfoil aerodynamic performance between the casting and laser-scan based simulations. For example, typical differences in iced-airfoil maximum lift coefficient were less than 3 percent with corresponding differences in stall angle of approximately 1 deg or less. The aerodynamic simulation accuracy reported in this paper has demonstrated the combined accuracy of the laser-scan and rapid-prototype manufacturing approach to simulating ice accretion for a NACA 23012 airfoil. For several

Validation of 3-D Ice Accretion Measurement Methodology for Experimental Aerodynamic Simulation

NASA Technical Reports Server (NTRS)

Broeren, Andy P.; Addy, Harold E., Jr.; Lee, Sam; Monastero, Marianne C.

2014-01-01

Determining the adverse aerodynamic effects due to ice accretion often relies on dry-air wind-tunnel testing of artificial, or simulated, ice shapes. Recent developments in ice accretion documentation methods have yielded a laser-scanning capability that can measure highly three-dimensional features of ice accreted in icing wind tunnels. The objective of this paper was to evaluate the aerodynamic accuracy of ice-accretion simulations generated from laser-scan data. Ice-accretion tests were conducted in the NASA Icing Research Tunnel using an 18-inch chord, 2-D straight wing with NACA 23012 airfoil section. For six ice accretion cases, a 3-D laser scan was performed to document the ice geometry prior to the molding process. Aerodynamic performance testing was conducted at the University of Illinois low-speed wind tunnel at a Reynolds number of 1.8 x 10(exp 6) and a Mach number of 0.18 with an 18-inch chord NACA 23012 airfoil model that was designed to accommodate the artificial ice shapes. The ice-accretion molds were used to fabricate one set of artificial ice shapes from polyurethane castings. The laser-scan data were used to fabricate another set of artificial ice shapes using rapid prototype manufacturing such as stereolithography. The iced-airfoil results with both sets of artificial ice shapes were compared to evaluate the aerodynamic simulation accuracy of the laser-scan data. For four of the six ice-accretion cases, there was excellent agreement in the iced-airfoil aerodynamic performance between the casting and laser-scan based simulations. For example, typical differences in iced-airfoil maximum lift coefficient were less than 3% with corresponding differences in stall angle of approximately one degree or less. The aerodynamic simulation accuracy reported in this paper has demonstrated the combined accuracy of the laser-scan and rapid-prototype manufacturing approach to simulating ice accretion for a NACA 23012 airfoil. For several of the ice-accretion

In Situ Quantification of Experimental Ice Accretion on Tree Crowns Using Terrestrial Laser Scanning

PubMed Central

Nock, Charles A.; Greene, David; Delagrange, Sylvain; Follett, Matt; Fournier, Richard; Messier, Christian

2013-01-01

In the eastern hardwood forests of North America ice storms are an important disturbance event. Ice storms strongly influence community dynamics as well as urban infrastructure via catastrophic branch failure; further, the severity and frequency of ice storms are likely to increase with climate change. However, despite a long-standing interest into the effects of freezing rain on forests, the process of ice accretion and thus ice loading on branches remains poorly understood. This is because a number of challenges have prevented in situ measurements of ice on branches, including: 1) accessing and measuring branches in tall canopies, 2) limitations to travel during and immediately after events, and 3) the unpredictability of ice storms. Here, utilizing a novel combination of outdoor experimental icing, manual measurements and terrestrial laser scanning (TLS), we perform the first in situ measurements of ice accretion on branches at differing heights in a tree crown and with increasing duration of exposure. We found that TLS can reproduce both branch and iced branch diameters with high fidelity, but some TLS instruments do not detect ice. Contrary to the expectations of ice accretion models, radial accretion varied sharply within tree crowns. Initially, radial ice accretion was similar throughout the crown, but after 6.5 hours of irrigation (second scanning) radial ice accretion was much greater on upper branches than on lower (∼factor of 3). The slope of the change in radial ice accretion along branches increased with duration of exposure and was significantly greater at the second scanning compared to the first. We conclude that outdoor icing experiments coupled with the use of TLS provide a robust basis for evaluation of models of ice accretion and breakage in tree crowns, facilitating estimation of the limiting breaking stress of branches by accurate measurements of ice loads. PMID:23741409

Modeling the Effects of Ice Accretion on the Low Pressure Compressor and the Overall Turbofan Engine System Performance

NASA Technical Reports Server (NTRS)

Veres, Joseph P.; Jorgenson, Philip C. E.; Wright, William B.

2011-01-01

The focus of this study is on utilizing a mean line compressor flow analysis code coupled to an engine system thermodynamic code, to estimate the effects of ice accretion on the low pressure compressor, and quantifying its effects on the engine system throughout a notional flight trajectory. In this paper a temperature range in which engine icing would occur was assumed. This provided a mechanism to locate potential component icing sites and allow the computational tools to add blockages due to ice accretion in a parametric fashion. Ultimately the location and level of blockage due to icing would be provided by an ice accretion code. To proceed, an engine system modeling code and a mean line compressor flow analysis code were utilized to calculate the flow conditions in the fan-core and low pressure compressor and to identify potential locations within the compressor where ice may accrete. In this study, an "additional blockage" due to the accretion of ice on the metal surfaces, has been added to the baseline aerodynamic blockage due to boundary layer, as well as the blade metal blockage. Once the potential locations of ice accretion are identified, the levels of additional blockage due to accretion were parametrically varied to estimate the effects on the low pressure compressor blade row performance operating within the engine system environment. This study includes detailed analysis of compressor and engine performance during cruise and descent operating conditions at several altitudes within the notional flight trajectory. The purpose of this effort is to develop the computer codes to provide a predictive capability to forecast the onset of engine icing events, such that they could ultimately help in the avoidance of these events.

Characteristics of surface roughness associated with leading edge ice accretion

NASA Technical Reports Server (NTRS)

Shin, Jaiwon

1994-01-01

Detailed size measurements of surface roughness associated with leading edge ice accretions are presented to provide information on characteristics of roughness and trends of roughness development with various icing parameters. Data was obtained from icing tests conducted in the Icing Research Tunnel (IRT) at NASA Lewis Research Center (LeRC) using a NACA 0012 airfoil. Measurements include diameters, heights, and spacing of roughness elements along with chordwise icing limits. Results confirm the existence of smooth and rough ice zones and that the boundary between the two zones (surface roughness transition region) moves upstream towards stagnation region with time. The height of roughness grows as the air temperature and the liquid water content increase, however, the airspeed has little effect on the roughness height. Results also show that the roughness in the surface roughness transition region grows during a very early stage of accretion but reaches a critical height and then remains fairly constant. Results also indicate that a uniformly distributed roughness model is only valid at a very initial stage of the ice accretion process.

Numerical simulation of ice accretion phenomena on rotor blade of axial blower

NASA Astrophysics Data System (ADS)

Matsuura, Taiki; Suzuki, Masaya; Yamamoto, Makoto; Shishido, Shinichiro; Murooka, Takeshi; Miyagawa, Hiroshi

2012-08-01

Ice accretion is the phenomenon that super-cooled water droplets impinge and accrete on a body. It is well known that ice accretion on blades and airfoils leads to performance degradation and severe accidents. For this reason, experimental investigations have been carried out using flight tests or icing tunnels. However, it is too expensive, dangerous, and difficult to set actual icing conditions. Hence, computational fluid dynamics is useful to predict ice accretion. A rotor blade is one of jet engine components where ice accretes. Therefore, the authors focus on the ice accretion on a rotor blade in this study. Three-dimensional icing phenomena on the rotor blade of a commercial axial blower are computed here, and ice accretion on the rotor blade is numerically investigated.

Distributed ice accretion sensor for smart aircraft structures

NASA Technical Reports Server (NTRS)

Gerardi, J. J.; Hickman, G. A.

1989-01-01

A distributed ice accretion sensor is presented, based on the concept of smart structures. Ice accretion is determined using spectral techniques to process signals from piezoelectric sensors integral to the airfoil skin. Frequency shifts in the leading edge structural skin modes are correlated to ice thickness. It is suggested that this method may be used to detect ice over large areas with minimal hardware. Results are presented from preliminary tests to measure simulated ice growth.

Fundamental Ice Crystal Accretion Physics Studies

NASA Technical Reports Server (NTRS)

Struk, Peter M.; Broeren, Andy P.; Tsao, Jen-Ching; Vargas, Mario; Wright, William B.; Currie, Tom; Knezevici, Danny; Fuleki, Dan

2012-01-01

Due to numerous engine power-loss events associated with high-altitude convective weather, ice accretion within an engine due to ice crystal ingestion is being investigated. The National Aeronautics and Space Administration (NASA) and the National Research Council (NRC) of Canada are starting to examine the physical mechanisms of ice accretion on surfaces exposed to ice-crystal and mixed-phase conditions. In November 2010, two weeks of testing occurred at the NRC Research Altitude Facility utilizing a single wedge-type airfoil designed to facilitate fundamental studies while retaining critical features of a compressor stator blade or guide vane. The airfoil was placed in the NRC cascade wind tunnel for both aerodynamic and icing tests. Aerodynamic testing showed excellent agreement compared with CFD data on the icing pressure surface and allowed calculation of heat transfer coefficients at various airfoil locations. Icing tests were performed at Mach numbers of 0.2 to 0.3, total pressures from 93 to 45 kPa, and total temperatures from 5 to 15 C. Ice and liquid water contents ranged up to 20 and 3 g/m3, respectively. The ice appeared well adhered to the surface in the lowest pressure tests (45 kPa) and, in a particular case, showed continuous leading-edge ice growth to a thickness greater than 15 mm in 3 min. Such widespread deposits were not observed in the highest pressure tests, where the accretions were limited to a small area around the leading edge. The suction surface was typically ice-free in the tests at high pressure, but not at low pressure. The icing behavior at high and low pressure appeared to be correlated with the wet-bulb temperature, which was estimated to be above 0 C in tests at 93 kPa and below 0 C in tests at lower pressure, the latter enhanced by more evaporative cooling of water. The authors believe that the large ice accretions observed in the low pressure tests would undoubtedly cause the aerodynamic performance of a compressor component

Fundamental Ice Crystal Accretion Physics Studies

NASA Technical Reports Server (NTRS)

Currie, Tom; Knezevici, Danny; Fuleki, Dan; Struk, Peter M.; Broeren, Andy P.; Tsao, Jen-ching; Vargas, Mario; Wright, William

2011-01-01

Due to numerous engine power-loss events associated with high-altitude convective weather, ice accretion within an engine due to ice-crystal ingestion is being investigated. The National Aeronautics and Space Administration (NASA) and the National Research Council (NRC) of Canada are starting to examine the physical mechanisms of ice accretion on surfaces exposed to ice-crystal and mixed-phase conditions. In November 2010, two weeks of testing occurred at the NRC Research Altitude Facility utilizing a single wedge-type airfoil designed to facilitate fundamental studies while retaining critical features of a compressor stator blade or guide vane. The airfoil was placed in the NRC cascade wind tunnel for both aerodynamic and icing tests. Aerodynamic testing showed excellent agreement compared with CFD data on the icing pressure surface and allowed calculation of heat transfer coefficients at various airfoil locations. Icing tests were performed at Mach numbers of 0.2 to 0.3, total pressures from 93 to 45 kPa, and total temperatures from 5 to 15 C. Ice and liquid water contents ranged up to 20 and 3 grams per cubic meter, respectively. The ice appeared well adhered to the surface in the lowest pressure tests (45 kPa) and, in a particular case, showed continuous leading-edge ice growth to a thickness greater than 15 millimeters in 3 minutes. Such widespread deposits were not observed in the highest pressure tests, where the accretions were limited to a small area around the leading edge. The suction surface was typically ice-free in the tests at high pressure, but not at low pressure. The icing behavior at high and low pressure appeared to be correlated with the wet-bulb temperature, which was estimated to be above 0 C in tests at 93 kPa and below 0 C in tests at lower pressure, the latter enhanced by more evaporative cooling of water. The authors believe that the large ice accretions observed in the low pressure tests would undoubtedly cause the aerodynamic

Multiphase Reactive Transport and Platelet Ice Accretion in the Sea Ice of McMurdo Sound, Antarctica

NASA Astrophysics Data System (ADS)

Buffo, J. J.; Schmidt, B. E.; Huber, C.

2018-01-01

Sea ice seasonally to interannually forms a thermal, chemical, and physical boundary between the atmosphere and hydrosphere over tens of millions of square kilometers of ocean. Its presence affects both local and global climate and ocean dynamics, ice shelf processes, and biological communities. Accurate incorporation of sea ice growth and decay, and its associated thermal and physiochemical processes, is underrepresented in large-scale models due to the complex physics that dictate oceanic ice formation and evolution. Two phenomena complicate sea ice simulation, particularly in the Antarctic: the multiphase physics of reactive transport brought about by the inhomogeneous solidification of seawater, and the buoyancy driven accretion of platelet ice formed by supercooled ice shelf water onto the basal surface of the overlying ice. Here a one-dimensional finite difference model capable of simulating both processes is developed and tested against ice core data. Temperature, salinity, liquid fraction, fluid velocity, total salt content, and ice structure are computed during model runs. The model results agree well with empirical observations and simulations highlight the effect platelet ice accretion has on overall ice thickness and characteristics. Results from sensitivity studies emphasize the need to further constrain sea ice microstructure and the associated physics, particularly permeability-porosity relationships, if a complete model of sea ice evolution is to be obtained. Additionally, implications for terrestrial ice shelves and icy moons in the solar system are discussed.

Samarium-Neodymium model age and Geochemical (Sr-Nd) signature of a bedrock inclusion from lake Vostok accretion ice.

NASA Astrophysics Data System (ADS)

Delmonte, B.; Petit, J. R.; Michard, A.; Basile-Doelsch, I.; Lipenkov, V.

2003-04-01

We investigated properties of the basal ice from Vostok ice core as well as the sediment inclusions within the accreted ice. The Vostok ice core preserves climatic information for the last 420 kyrs down to 3310m depth, but below this depth the horizontal layers of the climatic record are disrupted by the glacier dynamics. From 3450 m to 3538 m depth thin bedrock particles, as glacial flour, are entrapped. Glacial flour is released in the northern area lake, where glacier mostly melts and contributes to sediment accumulation. In the southern area, close to Vostok station, the lake water freezes and the upstream glacial flour does not contribute to sedimentation. The accreted ice contains visible sediment inclusions down to 3608 m (accretion ice 1), while below this depth and likely down to the water interface (˜3750 m), the ice is clear (accretion ice 2). The fine inclusions (1-2mm in diameter) from Accretion Ice 1 mostly consist of fine clays and quartz aggregates and we suggest they are entrained into ice as the glacier floats over shallow depth bay then it grounds against a relief rise. Afterward the glacier freely floats over the deep lake before reaching Vostok, and accreted ice 2 is clean. Sm-Nd dating of one of two inclusions at 3570 m depth gives 1.88 (+/-0.13)Ga (DM model age), corresponding to 1.47 Ga (TCHUR), suggesting a Precambrian origin. Also the isotopic signature of such inclusion (87Sr/86Sr= 0.8232 and eNd= -16) and that of a second one (87Sr/86Sr= 0.7999 and eNd= -15) are coherent with the nature of an old continental shield. Sediments that may initially accumulate in the shallow bay prior the Antarctic glaciation, should have been eroded and exported out of the lake by the glacier movement, this assuming processes for ice accretion and for sediment entrapping operate since a long time. As the glacial flour from upstream does not contribute to sedimentation, sediments need to be renewed at the surface of the bedrock rising question about the way

Users Manual for the NASA Lewis Ice Accretion Prediction Code (LEWICE)

NASA Technical Reports Server (NTRS)

Ruff, Gary A.; Berkowitz, Brian M.

1990-01-01

LEWICE is an ice accretion prediction code that applies a time-stepping procedure to calculate the shape of an ice accretion. The potential flow field is calculated in LEWICE using the Douglas Hess-Smith 2-D panel code (S24Y). This potential flow field is then used to calculate the trajectories of particles and the impingement points on the body. These calculations are performed to determine the distribution of liquid water impinging on the body, which then serves as input to the icing thermodynamic code. The icing thermodynamic model is based on the work of Messinger, but contains several major modifications and improvements. This model is used to calculate the ice growth rate at each point on the surface of the geometry. By specifying an icing time increment, the ice growth rate can be interpreted as an ice thickness which is added to the body, resulting in the generation of new coordinates. This procedure is repeated, beginning with the potential flow calculations, until the desired icing time is reached. The operation of LEWICE is illustrated through the use of five examples. These examples are representative of the types of applications expected for LEWICE. All input and output is discussed, along with many of the diagnostic messages contained in the code. Several error conditions that may occur in the code for certain icing conditions are identified, and a course of action is recommended. LEWICE has been used to calculate a variety of ice shapes, but should still be considered a research code. The code should be exercised further to identify any shortcomings and inadequacies. Any modifications identified as a result of these cases, or of additional experimental results, should be incorporated into the model. Using it as a test bed for improvements to the ice accretion model is one important application of LEWICE.

Microbiological and Biogeochemical Investigations of the Accreted Ice Above Subglacial Lake Vostok, Antarctica

NASA Astrophysics Data System (ADS)

Christner, B. C.; Foreman, C. F.; Arnold, B. R.; Welch, K. A.; Lyons, W. B.; Priscu, J. C.

2004-12-01

Subglacial Lake Vostok is located ~4 km beneath the surface of the East Antarctic ice sheet and has been isolated from the atmosphere for at least 15 million years. The lake has a surface area near 14,000 km2 and a depth exceeding 1000 m. While the nature of the environment within Subglacial Lake Vostok remains uncertain, if a sustained microbial ecosystem is present, life in this subsurface environment operates under arguably the most extreme conditions in the biosphere (i.e., high pressure, constant cold, high oxygen concentrations, and no light). The lake represents an analogue for ecosystems that may exist in Europa's ice-covered ocean and also provides an Earthly-based model for the evaluation of technology to search for life in icy extraterrestrial subsurface environments. Concerns for environmental protection have prevented direct sampling of the lake water thus far, as a prudent sampling plan that will not contaminate this pristine environment has yet to be developed and tested. However, an ice core has been retrieved at Vostok Station in which the bottom ~85 meters consists of lake water that has accreted to the bottom of the ice sheet, providing frozen samples of water from the lakes' surface. The ice from 3539 to 3609 mbs (accretion ice I) contains visible inclusions due to accretion in the shallow embayment or western grounding line, whereas ice from 3610-3623 mbs (accretion ice II) is very clean, forming above the deep eastern basin of the main lake. Using a multifaceted protocol to monitor cellular and molecular decontamination of ice cores, we show that the microbiology and geochemistry (i.e., dissolve organic carbon, nutrients, and ions) of accretion ice is very different from the overlying glacial ice. The numbers of cells are 2- to 7-fold higher in accretion ice I than in the overlying glacial ice, and decrease with increasing depth in accretion ice II. Cell viability in accretion ice samples has been confirmed by the measurable respiration of 14C

Structural properties of impact ices accreted on aircraft structures

NASA Technical Reports Server (NTRS)

Scavuzzo, R. J.; Chu, M. L.

1987-01-01

The structural properties of ice accretions formed on aircraft surfaces are studied. The overall objectives are to measure basic structural properties of impact ices and to develop finite element analytical procedures for use in the design of all deicing systems. The Icing Research Tunnel (IRT) was used to produce simulated natural ice accretion over a wide range of icing conditions. Two different test apparatus were used to measure each of the three basic mechanical properties: tensile, shear, and peeling. Data was obtained on both adhesive shear strength of impact ices and peeling forces for various icing conditions. The influences of various icing parameters such as tunnel air temperature and velocity, icing cloud drop size, material substrate, surface temperature at ice/material interface, and ice thickness were studied. A finite element analysis of the shear test apparatus was developed in order to gain more insight in the evaluation of the test data. A comparison with other investigators was made. The result shows that the adhesive shear strength of impact ice typically varies between 40 and 50 psi, with peak strength reaching 120 psi and is not dependent on the kind of substrate used, the thickness of accreted ice, and tunnel temperature below 4 C.

Reynolds-averaged Navier-Stokes based ice accretion for aircraft wings

NASA Astrophysics Data System (ADS)

Lashkajani, Kazem Hasanzadeh

This thesis addresses one of the current issues in flight safety towards increasing icing simulation capabilities for prediction of complex 2D and 3D glaze ice shapes over aircraft surfaces. During the 1980's and 1990's, the field of aero-icing was established to support design and certification of aircraft flying in icing conditions. The multidisciplinary technologies used in such codes were: aerodynamics (panel method), droplet trajectory calculations (Lagrangian framework), thermodynamic module (Messinger model) and geometry module (ice accretion). These are embedded in a quasi-steady module to simulate the time-dependent ice accretion process (multi-step procedure). The objectives of the present research are to upgrade the aerodynamic module from Laplace to Reynolds-Average Navier-Stokes equations solver. The advantages are many. First, the physical model allows accounting for viscous effects in the aerodynamic module. Second, the solution of the aero-icing module directly provides the means for characterizing the aerodynamic effects of icing, such as loss of lift and increased drag. Third, the use of a finite volume approach to solving the Partial Differential Equations allows rigorous mesh and time convergence analysis. Finally, the approaches developed in 2D can be easily transposed to 3D problems. The research was performed in three major steps, each providing insights into the overall numerical approaches. The most important realization comes from the need to develop specific mesh generation algorithms to ensure feasible solutions in very complex multi-step aero-icing calculations. The contributions are presented in chronological order of their realization. First, a new framework for RANS based two-dimensional ice accretion code, CANICE2D-NS, is developed. A multi-block RANS code from U. of Liverpool (named PMB) is providing the aerodynamic field using the Spalart-Allmaras turbulence model. The ICEM-CFD commercial tool is used for the iced airfoil

An experimental investigation of multi-element airfoil ice accretion and resulting performance degradation

NASA Technical Reports Server (NTRS)

Potapczuk, Mark G.; Berkowitz, Brian M.

1989-01-01

An investigation of the ice accretion pattern and performance characteristics of a multi-element airfoil was undertaken in the NASA Lewis 6- by 9-Foot Icing Research Tunnel. Several configurations of main airfoil, slat, and flaps were employed to examine the effects of ice accretion and provide further experimental information for code validation purposes. The text matrix consisted of glaze, rime, and mixed icing conditions. Airflow and icing cloud conditions were set to correspond to those typical of the operating environment anticipated tor a commercial transport vehicle. Results obtained included ice profile tracings, photographs of the ice accretions, and force balance measurements obtained both during the accretion process and in a post-accretion evaluation over a range of angles of attack. The tracings and photographs indicated significant accretions on the slat leading edge, in gaps between slat or flaps and the main wing, on the flap leading-edge surfaces, and on flap lower surfaces. Force measurments indicate the possibility of severe performance degradation, especially near C sub Lmax, for both light and heavy ice accretion and performance analysis codes presently in use. The LEWICE code was used to evaluate the ice accretion shape developed during one of the rime ice tests. The actual ice shape was then evaluated, using a Navier-Strokes code, for changes in performance characteristics. These predicted results were compared to the measured results and indicate very good agreement.

Aerodynamic Effects of Simulated Ice Accretion on a Generic Transport Model

NASA Technical Reports Server (NTRS)

Broeren, Andy P.; Lee, Sam; Shah, Gautam H.; Murphy, Patrick C.

2012-01-01

An experimental research effort was begun to develop a database of airplane aerodynamic characteristics with simulated ice accretion over a large range of incidence and sideslip angles. Wind-tunnel testing was performed at the NASA Langley 12-ft Low-Speed Wind Tunnel using a 3.5 percent scale model of the NASA Langley Generic Transport Model. Aerodynamic data were acquired from a six-component force and moment balance in static-model sweeps from alpha = -5deg to 85deg and beta = -45 deg to 45 deg at a Reynolds number of 0.24 x10(exp 6) and Mach number of 0.06. The 3.5 percent scale GTM was tested in both the clean configuration and with full-span artificial ice shapes attached to the leading edges of the wing, horizontal and vertical tail. Aerodynamic results for the clean airplane configuration compared favorably with similar experiments carried out on a 5.5 percent scale GTM. The addition of the large, glaze-horn type ice shapes did result in an increase in airplane drag coefficient but had little effect on the lift and pitching moment. The lateral-directional characteristics showed mixed results with a small effect of the ice shapes observed in some cases. The flow visualization images revealed the presence and evolution of a spanwise-running vortex on the wing that was the dominant feature of the flowfield for both clean and iced configurations. The lack of ice-induced performance and flowfield effects observed in this effort was likely due to Reynolds number effects for the clean configuration. Estimates of full-scale baseline performance were included in this analysis to illustrate the potential icing effects.

Experimental evidence for modifying the current physical model for ice accretion on aircraft surfaces

NASA Technical Reports Server (NTRS)

Olsen, W.; Walker, E.

1986-01-01

Closeup movies, still photographs, and other experimental data suggest that the current physical model for ice accretion needs significant modification. At aircraft airspeeds there was no flow of liquid over the surface of the ice after a short initial flow, even at barely subfreezing temperatures. Instead, there were very large stationary drops on the ice surface that lose water from their bottoms by freezing and replenish their liquid by catching the microscopic cloud droplets. This observation disagrees with the existing physical model, which assumes there is a thin liquid film continuously flowing over the ice surface. With no such flow, the freezing-fraction concept of the model fails when a mass balance is performed on the surface water. Rime ice does, as the model predicts, form when the air temperature is low enough to cause the cloud droplets to freeze almost immediately on impact. However, the characteristic shapes of horn-glaze ice or rime ice are primarily caused by the ice shape affecting the airflow locally and consequently the droplet catch and the resulting ice shape. Ice roughness greatly increases the heat transfer coefficient, stops the movement of drops along the surface, and may also affect the airflow initially and thereby the droplet catch. At high subreezing temperatures the initial flow and shedding of surface drops have a large effect on the ice shape. At the incipient freezing limit, no ice forms.

Ice Accretion Measurements on an Airfoil and Wedge in Mixed-Phase Conditions

NASA Technical Reports Server (NTRS)

Struk, Peter; Bartkus, Tadas; Tsao, Jen-Ching; Currie, Tom; Fuleki, Dan

2015-01-01

This paper describes ice accretion measurements from experiments conducted at the National Research Council (NRC) of Canada's Research Altitude Test Facility during 2012. Due to numerous engine power loss events associated with high altitude convective weather, potential ice accretion within an engine due to ice crystal ingestion is being investigated collaboratively by NASA and NRC. These investigations examine the physical mechanisms of ice accretion on surfaces exposed to ice crystal and mixed phase conditions, similar to those believed to exist in core compressor regions of jet engines. A further objective of these tests is to examine scaling effects since altitude appears to play a key role in this icing process.

Additional Improvements to the NASA Lewis Ice Accretion Code LEWICE

NASA Technical Reports Server (NTRS)

Wright, William B.; Bidwell, Colin S.

1995-01-01

Due to the feedback of the user community, three major features have been added to the NASA Lewis ice accretion code LEWICE. These features include: first, further improvements to the numerics of the code so that more time steps can be run and so that the code is more stable; second, inclusion and refinement of the roughness prediction model described in an earlier paper; third, inclusion of multi-element trajectory and ice accretion capabilities to LEWICE. This paper will describe each of these advancements in full and make comparisons with the experimental data available. Further refinement of these features and inclusion of additional features will be performed as more feedback is received.

Microwave ice accretion meter

NASA Technical Reports Server (NTRS)

Magenheim, Bertram (Inventor); Rocks, James K. (Inventor)

1984-01-01

A system for indicating ice thickness and rate of ice thickness growth on surfaces is disclosed. The region to be monitored for ice accretion is provided with a resonant surface waveguide which is mounted flush, below the surface being monitored. A controlled oscillator provides microwave energy via a feed point at a controllable frequency. A detector is coupled to the surface waveguide and is responsive to electrical energy. A measuring device indicates the frequency deviation of the controlled oscillator from a quiescent frequency. A control means is provided to control the frequency of oscillation of the controlled oscillator. In a first, open-loop embodiment, the control means is a shaft operated by an operator. In a second, closed-loop embodiment, the control means is a processor which effects automatic control.

Collection Efficiency and Ice Accretion Characteristics of Two Full Scale and One 1/4 Scale Business Jet Horizontal Tails

NASA Technical Reports Server (NTRS)

Bidwell, Colin S.; Papadakis, Michael

2005-01-01

Collection efficiency and ice accretion calculations have been made for a series of business jet horizontal tail configurations using a three-dimensional panel code, an adaptive grid code, and the NASA Glenn LEWICE3D grid based ice accretion code. The horizontal tail models included two full scale wing tips and a 25 percent scale model. Flow solutions for the horizontal tails were generated using the PMARC panel code. Grids used in the ice accretion calculations were generated using the adaptive grid code ICEGRID. The LEWICE3D grid based ice accretion program was used to calculate impingement efficiency and ice shapes. Ice shapes typifying rime and mixed icing conditions were generated for a 30 minute hold condition. All calculations were performed on an SGI Octane computer. The results have been compared to experimental flow and impingement data. In general, the calculated flow and collection efficiencies compared well with experiment, and the ice shapes appeared representative of the rime and mixed icing conditions for which they were calculated.

High-speed imaging of the transient ice accretion process on a NACA 0012 airfoil

NASA Astrophysics Data System (ADS)

Waldman, Rye; Hu, Hui

2014-11-01

Ice accretion on aircraft wings poses a performance and safety threat as aircraft encounter supercooled droplets suspended in the cloud layer. The details of the ice accretion depend on the atmospheric conditions and the fight parameters. We present the measurement results of the experiments conducted in the Iowa State icing wind tunnel on a NACA 0012 airfoil to study the transient ice accretion process under varying icing conditions. The icing process on the wing consists of a complex interaction of water deposition, surface water transport, and freezing. The aerodynamics affects the water deposition, the heat and mass transport, and ice accumulation; meanwhile, the accumulating ice also affects the aerodynamics. High-speed video of the unsteady icing accretion process was acquired under controlled environmental conditions to quantitatively measure the transient water run back, rivulet formation, and accumulated ice growth, and the experiments show how varying the environmental conditions modifies the ice accretion process. Funding support from the Iowa Energy Center with Grant No. 14-008-OG and National Science Foundation (NSF) with Grant No. CBET-1064196 and CBET-1438099 is gratefully acknowledged.

Rime ice accretion and its effect on airfoil performance. Ph.D. Thesis. Final Report

NASA Technical Reports Server (NTRS)

Bragg, M. B.

1982-01-01

A methodology was developed to predict the growth of rime ice, and the resulting aerodynamic penalty on unprotected, subcritical, airfoil surfaces. The system of equations governing the trajectory of a water droplet in the airfoil flowfield is developed and a numerical solution is obtained to predict the mass flux of super cooled water droplets freezing on impact. A rime ice shape is predicted. The effect of time on the ice growth is modeled by a time-stepping procedure where the flowfield and droplet mass flux are updated periodically through the ice accretion process. Two similarity parameters, the trajectory similarity parameter and accumulation parameter, are found to govern the accretion of rime ice. In addition, an analytical solution is presented for Langmuir's classical modified inertia parameter. The aerodynamic evaluation of the effect of the ice accretion on airfoil performance is determined using an existing airfoil analysis code with empirical corrections. The change in maximum lift coefficient is found from an analysis of the new iced airfoil shape. The drag correction needed due to the severe surface roughness is formulated from existing iced airfoil and rough airfoil data. A small scale wind tunnel test was conducted to determine the change in airfoil performance due to a simulated rime ice shape.

Navier-Stokes analysis of airfoils with leading edge ice accretions

NASA Technical Reports Server (NTRS)

Potapczuk, Mark G.

1993-01-01

A numerical analysis of the flowfield characteristics and the performance degradation of an airfoil with leading edge ice accretions was performed. The important fluid dynamic processes were identified and calculated. Among these were the leading edge separation bubble at low angles of attack, complete separation on the low pressure surface resulting in premature shell, drag rise due to the ice shape, and the effects of angle of attack on the separated flow field. Comparisons to experimental results were conducted to confirm these calculations. A computer code which solves the Navier-Stokes equations in two dimensions, ARC2D, was used to perform the calculations. A Modified Mixing Length turbulence model was developed to produce grids for several ice shape and airfoil combinations. Results indicate that the ability to predict overall performance characteristics, such as lift and drag, at low angles of attack is excellent. Transition location is important for accurately determining separation bubble shape. Details of the flowfield in and downstream of the separated regions requires some modifications. Calculations for the stalled airfoil indicate periodic shedding of vorticity that was generated aft of the ice accretion. Time averaged pressure values produce results which compare favorably with experimental information. A turbulence model which accounts for the history effects in the flow may be justified.

Modeling Commercial Turbofan Engine Icing Risk With Ice Crystal Ingestion

NASA Technical Reports Server (NTRS)

Jorgenson, Philip C. E.; Veres, Joseph P.

2013-01-01

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that have been attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was degraded engine performance, and one or more of the following: loss of thrust control (roll back), compressor surge or stall, and flameout of the combustor. As ice crystals are ingested into the fan and low pressure compression system, the increase in air temperature causes a portion of the ice crystals to melt. It is hypothesized that this allows the ice-water mixture to cover the metal surfaces of the compressor stationary components which leads to ice accretion through evaporative cooling. Ice accretion causes a blockage which subsequently results in the deterioration in performance of the compressor and engine. The focus of this research is to apply an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion. The tool is comprised of an engine system thermodynamic cycle code, a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor flow path, without modeling the actual ice accretion. A commercial turbofan engine which has previously experienced icing events during operation in a high altitude ice crystal environment has been tested in the Propulsion Systems Laboratory (PSL) altitude test facility at NASA Glenn Research Center. The PSL has the capability to produce a continuous ice cloud which are ingested by the engine during operation over a range of altitude conditions. The PSL test results confirmed that there was ice accretion in the engine due to ice crystal ingestion, at the same simulated altitude operating conditions as experienced previously in

An Approach to Detect and Mitigate Ice Particle Accretion in Aircraft Engine Compression Systems

NASA Technical Reports Server (NTRS)

May, Ryan D.; Guo, Ten-Huei; Simon, Donald L.

2013-01-01

The accretion of ice in the compression system of commercial gas turbine engines operating in high ice water content conditions is a safety issue being studied by the aviation sector. While most of the research focuses on the underlying physics of ice accretion and the meteorological conditions in which accretion can occur, a systems-level perspective on the topic lends itself to potential near-term operational improvements. This work focuses on developing an accurate and reliable algorithm for detecting the accretion of ice in the low pressure compressor of a generic 40,000 lbf thrust class engine. The algorithm uses only the two shaft speed sensors and works regardless of engine age, operating condition, and power level. In a 10,000-case Monte Carlo simulation, the detection approach was found to have excellent capability at determining ice accretion from sensor noise with detection occurring when ice blocks an average of 6.8 percent of the low pressure compressor area. Finally, an initial study highlights a potential mitigation strategy that uses the existing engine actuators to raise the temperature in the low pressure compressor in an effort to reduce the rate at which ice accretes.

An Approach to Detect and Mitigate Ice Particle Accretion in Aircraft Engine Compression Systems

NASA Technical Reports Server (NTRS)

May, Ryan D.; Guo, Ten-Huei; Simon, Donald L.

2013-01-01

The accretion of ice in the compression system of commercial gas turbine engines operating in high ice water content conditions is a safety issue being studied by the aviation sector. While most of the research focuses on the underlying physics of ice accretion and the meteorological conditions in which accretion can occur, a systems-level perspective on the topic lends itself to potential near-term operational improvements. This work focuses on developing an accurate and reliable algorithm for detecting the accretion of ice in the low pressure compressor of a generic 40,000 lbf thrust class engine. The algorithm uses only the two shaft speed sensors and works regardless of engine age, operating condition, and power level. In a 10,000-case Monte Carlo simulation, the detection approach was found to have excellent capability at determining ice accretion from sensor noise with detection occurring when ice blocks an average of 6.8% of the low pressure compressor area. Finally, an initial study highlights a potential mitigation strategy that uses the existing engine actuators to raise the temperature in the low pressure compressor in an effort to reduce the rate at which ice accretes.

Development and test of a Microwave Ice Accretion Measurement Instrument (MIAMI)

NASA Technical Reports Server (NTRS)

Magenheim, B.; Rocks, J. K.

1982-01-01

The development of an ice accretion measurement instrument that is a highly sensitive, accurate, rugged and reliable microprocessor controlled device using low level microwave energy for non-instrusive real time measurement and recording of ice growth history, including ice thickness and accretion rate is discussed. Data is displayed and recorded digitally. New experimental data is presented, obtained with the instrument, which demonstrates its ability to measure ice growth on a two-dimensional airfoil. The device is suitable for aircraft icing protection. It may be mounted flush, non-intrusively, on any part of an aircraft skin including rotor blades and engine inlets.

Potential flow analysis of glaze ice accretions on an airfoil

NASA Technical Reports Server (NTRS)

Zaguli, R. J.

1984-01-01

The results of an analytical/experimental study of the flow fields about an airfoil with leading edge glaze ice accretion shapes are presented. Tests were conducted in the Icing Research Tunnel to measure surface pressure distributions and boundary layer separation reattachment characteristics on a general aviation wing section to which was affixed wooden ice shapes which approximated typical glaze ice accretions. Comparisons were made with predicted pressure distributions using current airfoil analysis codes as well as the Bristow mixed analysis/design airfoil panel code. The Bristow code was also used to predict the separation reattachment dividing streamline by inputting the appropriate experimental surface pressure distribution.

Prediction of ice accretion on a swept NACA 0012 airfoil and comparisons to flight test results

NASA Technical Reports Server (NTRS)

Reehorst, Andrew L.

1992-01-01

In the winter of 1989-90, an icing research flight project was conducted to obtain swept wing ice accretion data. Utilizing the NASA Lewis Research Center's DHC-6 DeHavilland Twin Otter aircraft, research flights were made into known icing conditions in Northeastern Ohio. The icing cloud environment and aircraft flight data were measured and recorded by an onboard data acquisition system. Upon entry into the icing environment, a 24 inch span, 15 inch chord NACA 0012 airfoil was extended from the aircraft and set to the desired sweep angle. After the growth of a well defined ice shape, the airfoil was retracted into the aircraft cabin for ice shape documentation. The ice accretions were recorded by ice tracings and photographs. Ice accretions were mostly of the glaze type and exhibited scalloping. The ice was accreted at sweep angles of 0, 30, and 45 degrees. A 3-D ice accretion prediction code was used to predict ice profiles for five selected flight test runs, which include sweep angle of zero, 30, and 45 degrees. The code's roughness input parameter was adjusted for best agreement. A simple procedure was added to the code to account for 3-D ice scalloping effects. The predicted ice profiles are compared to their respective flight test counterparts. This is the first attempt to predict ice profiles on swept wings with significant scalloped ice formations.

Experimental Investigation of Ice Accretion Effects on a Swept Wing

NASA Technical Reports Server (NTRS)

Wong, S. C.; Vargas, M.; Papadakis, M.; Yeong, H. W.; Potapczuk, M.

2005-01-01

An experimental investigation was conducted to study the effects of 2-, 5-, 10-, and 22.5-min ice accretions on the aerodynamic performance of a swept finite wing. The ice shapes tested included castings of ice accretions obtained from icing tests at the NASA Glenn Icing Research Tunnel (IRT) and simulated ice shapes obtained with the LEWICE 2.0 ice accretion code. The conditions used for the icing tests were selected to provide five glaze ice shapes with complete and incomplete scallop features and a small rime ice shape. The LEWICE ice shapes were defined for the same conditions as those used in the icing tests. All aerodynamic performance tests were conducted in the 7- x 10-ft Low-Speed Wind Tunnel Facility at Wichita State University. Six component force and moment measurements, aileron hinge moments, and surface pressures were obtained for a Reynolds number of 1.8 million based on mean aerodynamic chord and aileron deflections in the range of -15o to 20o. Tests were performed with the clean wing, six IRT ice shape castings, seven smooth LEWICE ice shapes, and seven rough LEWICE ice shapes. Roughness for the LEWICE ice shapes was simulated with 36-size grit. The experiments conducted showed that the glaze ice castings reduced the maximum lift coefficient of the clean wing by 11.5% to 93.6%, while the 5-min rime ice casting increased maximum lift by 3.4%. Minimum iced wing drag was 133% to 3533% greater with respect to the clean case. The drag of the iced wing near the clean wing stall angle of attack was 17% to 104% higher than that of the clean case. In general, the aileron remained effective in changing the lift of the clean and iced wings for all angles of attack and aileron deflections tested. Aileron hinge moments for the iced wing cases remained within the maximum and minimum limits defined by the clean wing hinge moments. Tests conducted with the LEWICE ice shapes showed that in general the trends in aerodynamic performance degradation of the wing with

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

NASA Technical Reports Server (NTRS)

Kreeger, Richard E.; Wright, William B.

2005-01-01

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.

Ultra-low rare earth element content in accreted ice from sub-glacial Lake Vostok, Antarctica

NASA Astrophysics Data System (ADS)

Gabrielli, Paolo; Planchon, Frederic; Barbante, Carlo; Boutron, Claude F.; Petit, Jean Robert; Bulat, Sergey; Hong, Sungmin; Cozzi, Giulio; Cescon, Paolo

2009-10-01

This paper reports the first rare earth element (REE) concentrations in accreted ice refrozen from sub-glacial Lake Vostok (East Antarctica). REE were determined in various sections of the Vostok ice core in order to geochemically characterize its impurities. Samples were obtained from accreted ice and, for comparison, from the upper glacier ice of atmospheric origin (undisturbed, disturbed and glacial flour ice). REE concentrations ranged between 0.8-56 pg g -1 for Ce and 0.0035-0.24 pg g -1 for Lu in glacier ice, and between <0.1-24 pg g -1 for Ce and <0.0004-0.02 pg g -1 for Lu in accreted ice. Interestingly, the REE concentrations in the upper accreted ice (AC 1; characterized by visible aggregates containing a mixture of very fine terrigenous particles) and in the deeper accreted ice (AC 2; characterized by transparent ice) are lower than those in fresh water and seawater, respectively. We suggest that such ultra-low concentrations are unlikely to be representative of the real REE content in Lake Vostok, but instead may reflect phase exclusion processes occurring at the ice/water interface during refreezing. In particular, the uneven spatial distribution (on the order of a few cm) and the large range of REE concentrations observed in AC 1 are consistent with the occurrence/absence of the aggregates in adjacent ice, and point to the presence of solid-phase concentration/exclusion processes occurring within separate pockets of frazil ice during AC 1 formation. Interestingly, if the LREE enrichment found in AC 1 was not produced by chemical fractionation occurring in Lake Vostok water, this may reflect a contribution of bedrock material, possibly in combination with aeolian dust released into the lake by melting of the glacier ice. Collectively, this valuable information provides new insight into the accreted ice formation processes, the bedrock geology of East Antarctica as well as the water chemistry and circulation of Lake Vostok.

Ultra-low rare earth element content in accreted ice from sub-glacial Lake Vostok, Antarctica

NASA Astrophysics Data System (ADS)

Barbante, C.; Gabrielli, P.; Turetta, C.; Planchon, F.; Boutron, C.; Petit, J. R.; Bulat, S.; Hong, S.; Cozzi, G.; Cescon, P.

2009-12-01

We report the first rare earth element (REE) concentrations in accreted ice refrozen from sub-glacial Lake Vostok (East Antarctica). REE were determined in various sections of the Vostok ice core in order to geochemically characterize its impurities. Samples were obtained from accreted ice and, for comparison, from the upper glacier ice of atmospheric origin (undisturbed, disturbed and glacial flour ice). REE concentrations ranged between 0.8-56 pg g-1 for Ce and 0.0035- 0.24 pg g-1 for Lu in glacier ice, and between <0.1-24 pg g-1 for Ce and <0.0004-0.02 pg g-1 for Lu in accreted ice. Interestingly, the REE concentrations in the upper accreted ice (AC1;characterized by visible aggregates containing a mixture of very fine terrigenous particles) and in the deeper accreted ice (AC2; characterized by transparent ice) are lower than those in fresh water and seawater, respectively. We suggest that such ultra-low concentrations are unlikely to be representative of the real REE content in Lake Vostok, but instead may reflect phase exclusion processes occurring at the ice/water interface during refreezing. In particular, the uneven spatial distribution (on the order of a few cm) and the large range of REE concentrations observed in AC1 are consistent with the occurrence/absence of the aggregates in adjacent ice, and point to the presence of solid-phase concentration/exclusion processes occurring within separate pockets of frazil ice during AC1 formation. Interestingly, if the LREE enrichment found in AC1 was not produced by chemical fractionation occurring in Lake Vostok water, this may reflect a contribution of bedrock material, possibly in combination with aeolian dust released into the lake by melting of the glacier ice. Collectively, this valuable information provides new insight into the accreted ice formation processes, the bedrock geology of East Antarctica as well as the water chemistry and circulation of Lake Vostok.

Detection of the Impact of Ice Crystal Accretion in an Aircraft Engine Compression System During Dynamic Operation

NASA Technical Reports Server (NTRS)

May, Ryan D.; Simon, Donald L.; Guo, Ten-Huei

2014-01-01

The accretion of ice in the compression system of commercial gas turbine engines operating in high ice water content conditions is a safety issue being studied by the aviation community. While most of the research focuses on the underlying physics of ice accretion and the meteorological conditions in which accretion can occur, a systems-level perspective on the topic lends itself to potential near-term operational improvements. Here a detection algorithm is developed which has the capability to detect the impact of ice accretion in the Low Pressure Compressor of an aircraft engine during steady flight as well as during changes in altitude. Unfortunately, the algorithm as implemented was not able to distinguish throttle changes from ice accretion and thus more work remains to be done.

Water Droplet Impingement on Simulated Glaze, Mixed, and Rime Ice Accretions

NASA Technical Reports Server (NTRS)

Papadakis, Michael; Rachman, Arief; Wong, See-Cheuk; Yeong, Hsiung-Wei; Hung, Kuohsing E.; Vu, Giao T.; Bidwell, Colin S.

2007-01-01

Water droplet impingement data were obtained at the NASA Glenn Icing Research Tunnel (IRT) for a 36-in. chord NACA 23012 airfoil with and without simulated ice using a dye-tracer method. The simulated ice shapes were defined with the NASA Glenn LEWICE 2.2 ice accretion program and including one rime, four mixed and five glaze ice shapes. The impingement experiments were performed with spray clouds having median volumetric diameters of 20, 52, 111, 154, and 236 micron. Comparisons to the experimental data were generated which showed good agreement for the rime and mixed shapes at lower drop sizes. For larger drops sizes LEWICE 2.2 over predicted the collection efficiencies due to droplet splashing effects which were not modeled in the program. Also for the more complex glaze ice shapes interpolation errors resulted in the over prediction of collection efficiencies in cove or shadow regions of ice shapes.

Changes in evaporation and potential hazards associated with ice accretion in a "New Arctic"

NASA Astrophysics Data System (ADS)

Boisvert, L.

2016-12-01

The Arctic sea ice acts as a barrier between the ocean and atmosphere inhibiting the exchange of heat, momentum, and moisture. Recently, the Arctic has seen unprecedented declines in the summer sea ice area, changing to a "New Arctic" climate system, one that is dominated by processes affected by large ice-free areas for the majority of the year as the melt season lengthens. Using atmospheric data from the Atmospheric Infrared Sounder (AIRS) instrument, we found that accompanying this loss of sea ice, the Arctic is becoming warmer and wetter. Evaporation, which plays an important role in the Arctic energy budget, water vapor feedback, and Arctic amplification, is also changing. The largest increases seen in evaporation are in the Arctic coastal seas during the spring and fall where there has been a reduction in sea ice cover and an increase in sea surface temperatures. Increases in evaporation also correspond to increases in low-level clouds. In this "New Arctic" transportation and shipping throughout the Arctic Ocean is becoming a more viable option as the areas in which ships can travel and the time period for ship travel continue to increase. There are various hazards associated with Arctic shipping, one being ice accretion. Ice accretion is the build up of ice on the surface of ships as they travel through regions of specific meteorological conditions unique to high-latitude environments. Besides reduced visibility, this build up of ice can cause ships to sink or capsize (by altering the ships center of gravity) depending on the severity and/or the location of ice build-up. With these changing atmospheric conditions in the Arctic, we expect there have been increases in the ice accretion potential over recent years, and an increase in the likelihood of high, and potentially dangerous ice accretion rates. Improved understanding of how this rapid loss of sea ice affects the "New Arctic" climate system, how evaporation is changing and how ice accretion could change

An experimental study of the aerodynamics of a NACA 0012 airfoil with a simulated glaze ice accretion

NASA Technical Reports Server (NTRS)

Bragg, M. B.

1986-01-01

An experimental study was conducted in the Ohio State University subsonic wind tunnel to measure the detailed aerodynamic characteristics of an airfoil with a simulated glaze ice accretion. A NACA 0012 model with interchangeable leading edges and pressure taps every one percent chord was used. Surface pressure and wake data were taken on the airfoil clean, with forced transition and with a simulated glaze ice shape. Lift and drag penalties due to the ice shape were found and the surface pressure clearly showed that large separation bubbles were present. Both total pressure and split-film probes were used to measure velocity profiles, both for the clean model and for the model with a simulated ice accretion. A large region of flow separation was seen in the velocity profiles and was correlated to the pressure measurements. Clean airfoil data were found to compare well to existing airfoil analysis methods.

Ice Accretion and Performance Degradation Calculations with LEWICE/NS

NASA Technical Reports Server (NTRS)

Potapczuk, Mark G.; Al-Khalil, Kamel M.; Velazquez, Matthew T.

1993-01-01

The LEWICE ice accretion computer code has been extended to include the solution of the two-dimensional Navier-Stokes 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, calculate aeroperformance changes, and plot results. The new elements of the code are described. Calculated results are compared to experiment for several cases, including both ice shape and drag rise.

A Study of Large Droplet Ice Accretions in the NASA-Lewis IRT at Near-Freezing Conditions

NASA Technical Reports Server (NTRS)

Miller, Dean R.; Addy, Harold E. , Jr.; Ide, Robert F.

1996-01-01

This report documents the results of an experimental study on large droplet ice accretions which was conducted in the NASA-Lewis Icing Research Tunnel (IRT) with a full-scale 77.25 inch chord Twin-Otter wing section. This study was intended to: (1) document the existing capability of the IRT to produce a large droplet icing cloud, and (2) study the effect of various parameters on large droplet ice accretions. Results are presented from a study of the IRT's capability to produce large droplets with MVD of 99 and 160 microns. The effect of the initial water droplet temperature on the resultant ice accretion was studied for different initial spray bar air and water temperatures. The initial spray bar water temperature was found to have no discernible effect upon the large droplet ice accretions. Also, analytical and experimental results suggest that the water droplet temperature is very nearly the same as the tunnel ambient temperature, thus providing a realistic simulation of the large droplet natural icing condition. The effect of temperature, droplet size, airspeed, angle-of attack, flap setting and de-icer boot cycling time on ice accretion was studied, and will be discussed in this report. It was found that, in almost all of the cases studied, an ice ridge formed immediately aft of the active portion of the de-icer boot. This ridge was irregular in shape, varied in location, and was in some cases discontinuous due to aerodynamic shedding.

Mixed Phase Modeling in GlennICE with Application to Engine Icing

NASA Technical Reports Server (NTRS)

Wright, William B.; Jorgenson, Philip C. E.; Veres, Joseph P.

2011-01-01

A capability for modeling ice crystals and mixed phase icing has been added to GlennICE. Modifications have been made to the particle trajectory algorithm and energy balance to model this behavior. This capability has been added as part of a larger effort to model ice crystal ingestion in aircraft engines. Comparisons have been made to four mixed phase ice accretions performed in the Cox icing tunnel in order to calibrate an ice erosion model. A sample ice ingestion case was performed using the Energy Efficient Engine (E3) model in order to illustrate current capabilities. Engine performance characteristics were supplied using the Numerical Propulsion System Simulation (NPSS) model for this test case.

Aerodynamic Classification of Swept-Wing Ice Accretion

NASA Technical Reports Server (NTRS)

Broeren, Andy; Diebold, Jeff; Bragg, Mike

2013-01-01

The continued design, certification and safe operation of swept-wing airplanes in icing conditions rely on the advancement of computational and experimental simulation methods for higher fidelity results over an increasing range of aircraft configurations and performance, and icing conditions. The current state-of-the-art in icing aerodynamics is mainly built upon a comprehensive understanding of two-dimensional geometries that does not currently exist for fundamentally three-dimensional geometries such as swept wings. The purpose of this report is to describe what is known of iced-swept-wing aerodynamics and to identify the type of research that is required to improve the current understanding. Following the method used in a previous review of iced-airfoil aerodynamics, this report proposes a classification of swept-wing ice accretion into four groups based upon unique flowfield attributes. These four groups are: ice roughness, horn ice, streamwise ice, and spanwise-ridge ice. For all of the proposed ice-shape classifications, relatively little is known about the three-dimensional flowfield and even less about the effect of Reynolds number and Mach number on these flowfields. The classifications and supporting data presented in this report can serve as a starting point as new research explores swept-wing aerodynamics with ice shapes. As further results are available, it is expected that these classifications will need to be updated and revised.

Close-up analysis of inflight ice accretion

NASA Technical Reports Server (NTRS)

Reehorst, Andrew L.; Ratvasky, Thomas P.; Sims, James

1994-01-01

The objective of this effort was to validate in flight, data that has been gathered in the NASA Lewis Research Center's Icing Research Tunnel (IRT) over the past several years. All data was acquired in flight on the NASA Lewis Research Center's Twin Otter Icing Research Aircraft. A faired 3.5 in. diameter metal-clad cylinder exposed to the natural icing environment was observed by a close-up video camera. The grazing angle video footage was recorded to S-VHS video tape and after the icing encounter, the resultant ice shape was documented by 35 mm photography and pencil tracings. The feather growth area was of primary interest; however, all regions of the ice accretion, from the stagnation line to the aft edge of run back were observed and recorded. After analysis of the recorded data several interesting points became evident: (1) the measured flight feather growth rate is consistent with IRT values, (2) the feather growth rate appears to be influenced by droplet size, (3) the feathers were straighter in the lower, spottier LWC of flight in comparison to those observed in the IRT, (4) feather shedding and ice sublimation may be significant to the final ice shape, and (5) the snow encountered on these flights appeared to have little influence on ice growth.

Users manual for the improved NASA Lewis ice accretion code LEWICE 1.6

NASA Technical Reports Server (NTRS)

Wright, William B.

1995-01-01

This report is intended as an update/replacement to NASA CR 185129 'User's Manual for the NASALewis Ice Accretion Prediction Code (LEWICE)' and as an update to NASA CR 195387 'Update to the NASA Lewis Ice Accretion Code LEWICE'. In addition to describing the changes specifically made for this version, information from previous manuals will be duplicated so that the user will not need three manuals to use this code.

An experimental and theoretical study of the ice accretion process during artificial and natural icing conditions

NASA Technical Reports Server (NTRS)

Kirby, Mark S.; Hansman, R. John

1988-01-01

Real-time measurements of ice growth during artificial and natural icing conditions were conducted using an ultrasonic pulse-echo technique. This technique allows ice thickness to be measured with an accuracy of + or - 0.5 mm; in addition, the ultrasonic signal characteristics may be used to detect the presence of liquid on the ice surface and hence discern wet and dry ice growth behavior. Ice growth was measured on the stagnation line of a cylinder exposed to artificial icing conditions in the NASA Lewis Icing Research Tunnel (IRT), and similarly for a cylinder exposed in flight to natural icing conditions. Ice thickness was observed to increase approximately linearly with exposure time during the initial icing period. The ice accretion rate was found to vary with cloud temperature during wet ice growth, and liquid runback from the stagnation region was inferred. A steady-state energy balance model for the icing surface was used to compare heat transfer characteristics for IRT and natural icing conditions. Ultrasonic measurements of wet and dry ice growth observed in the IRT and in flight were compared with icing regimes predicted by a series of heat transfer coefficients. The heat transfer magnitude was generally inferred to be higher for the IRT than for the natural icing conditions encountered in flight. An apparent variation in the heat transfer magnitude was also observed for flights conducted through different natural icing-cloud formations.

Aerodynamic Classification of Swept-Wing Ice Accretion

NASA Technical Reports Server (NTRS)

Diebold, Jeff M.; Broeren, Andy P.; Bragg, Michael B.

2013-01-01

The continued design, certification and safe operation of swept-wing airplanes in icing conditions rely on the advancement of computational and experimental simulation methods for higher fidelity results over an increasing range of aircraft configurations and performance, and icing conditions. The current stateof- the-art in icing aerodynamics is mainly built upon a comprehensive understanding of two-dimensional geometries that does not currently exist for fundamentally three-dimensional geometries such as swept wings. The purpose of this report is to describe what is known of iced-swept-wing aerodynamics and to identify the type of research that is required to improve the current understanding. Following the method used in a previous review of iced-airfoil aerodynamics, this report proposes a classification of swept-wing ice accretion into four groups based upon unique flowfield attributes. These four groups are: ice roughness, horn ice, streamwise ice and spanwise-ridge ice. In the case of horn ice it is shown that a further subclassification of "nominally 3D" or "highly 3D" horn ice may be necessary. For all of the proposed ice-shape classifications, relatively little is known about the three-dimensional flowfield and even less about the effect of Reynolds number and Mach number on these flowfields. The classifications and supporting data presented in this report can serve as a starting point as new research explores swept-wing aerodynamics with ice shapes. As further results are available, it is expected that these classifications will need to be updated and revised.

Aerodynamic Classification of Swept-Wing Ice Accretion

NASA Technical Reports Server (NTRS)

Diebold, Jeff M.; Broeren, Andy P.; Bragg, Michael B.

2013-01-01

The continued design, certification and safe operation of swept-wing airplanes in icing conditions rely on the advancement of computational and experimental simulation methods for higher fidelity results over an increasing range of aircraft configurations and performance, and icing conditions. The current state-of-the-art in icing aerodynamics is mainly built upon a comprehensive understanding of two-dimensional geometries that does not currently exist for fundamentally three-dimensional geometries such as swept wings. The purpose of this report is to describe what is known of iced-swept-wing aerodynamics and to identify the type of research that is required to improve the current understanding. Following the method used in a previous review of iced-airfoil aerodynamics, this report proposes a classification of swept-wing ice accretion into four groups based upon unique flowfield attributes. These four groups are: ice roughness, horn ice, streamwise ice and spanwise-ridge ice. In the case of horn ice it is shown that a further subclassification of nominally 3D or highly 3D horn ice may be necessary. For all of the proposed ice-shape classifications, relatively little is known about the three-dimensional flowfield and even less about the effect of Reynolds number and Mach number on these flowfields. The classifications and supporting data presented in this report can serve as a starting point as new research explores swept-wing aerodynamics with ice shapes. As further results are available, it is expected that these classifications will need to be updated and revised.

Ice Crystal Icing Research at NASA

NASA Technical Reports Server (NTRS)

Flegel, Ashlie B.

2017-01-01

Ice crystals found at high altitude near convective clouds are known to cause jet engine power-loss events. These events occur due to ice crystals entering a propulsion system's core flowpath and accreting ice resulting in events such as uncommanded loss of thrust (rollback), engine stall, surge, and damage due to ice shedding. As part of a community with a growing need to understand the underlying physics of ice crystal icing, NASA has been performing experimental efforts aimed at providing datasets that can be used to generate models to predict the ice accretion inside current and future engine designs. Fundamental icing physics studies on particle impacts, accretion on a single airfoil, and ice accretions observed during a rollback event inside a full-scale engine in the Propulsion Systems Laboratory are summarized. Low fidelity code development using the results from the engine tests which identify key parameters for ice accretion risk and the development of high fidelity codes are described. These activities have been conducted internal to NASA and through collaboration efforts with industry, academia, and other government agencies. The details of the research activities and progress made to date in addressing ice crystal icing research challenges are discussed.

Ice Crystal Icing Research at NASA

NASA Technical Reports Server (NTRS)

Flegel, Ashlie B.

2017-01-01

Ice crystals found at high altitude near convective clouds are known to cause jet engine power-loss events. These events occur due to ice crystals entering a propulsion systems core flowpath and accreting ice resulting in events such as uncommanded loss of thrust (rollback), engine stall, surge, and damage due to ice shedding. As part of a community with a growing need to understand the underlying physics of ice crystal icing, NASA has been performing experimental efforts aimed at providing datasets that can be used to generate models to predict the ice accretion inside current and future engine designs. Fundamental icing physics studies on particle impacts, accretion on a single airfoil, and ice accretions observed during a rollback event inside a full-scale engine in the Propulsion Systems Laboratory are summarized. Low fidelity code development using the results from the engine tests which identify key parameters for ice accretion risk and the development of high fidelity codes are described. These activities have been conducted internal to NASA and through collaboration efforts with industry, academia, and other government agencies. The details of the research activities and progress made to date in addressing ice crystal icing research challenges are discussed.

Development of 3-D Ice Accretion Measurement Method

NASA Technical Reports Server (NTRS)

Lee, Sam; Broeren, Andy P.; Addy, Harold E., Jr.; Sills, Robert; Pifer, Ellen M.

2012-01-01

A research plan is currently being implemented by NASA to develop and validate the use of a commercial laser scanner to record and archive fully three-dimensional (3-D) ice shapes from an icing wind tunnel. The plan focused specifically upon measuring ice accreted in the NASA Icing Research Tunnel (IRT). The plan was divided into two phases. The first phase was the identification and selection of the laser scanning system and the post-processing software to purchase and develop further. The second phase was the implementation and validation of the selected system through a series of icing and aerodynamic tests. Phase I of the research plan has been completed. It consisted of evaluating several scanning hardware and software systems against an established selection criteria through demonstrations in the IRT. The results of Phase I showed that all of the scanning systems that were evaluated were equally capable of scanning ice shapes. The factors that differentiated the scanners were ease of use and the ability to operate in a wide range of IRT environmental conditions.

Implementation of Combined Feather and Surface-Normal Ice Growth Models in LEWICE/X

NASA Technical Reports Server (NTRS)

Velazquez, M. T.; Hansman, R. J., Jr.

1995-01-01

Experimental observations have shown that discrete rime ice growths called feathers, which grow in approximately the direction of water droplet impingement, play an important role in the growth of ice on accreting surfaces for some thermodynamic conditions. An improved physical model of ice accretion has been implemented in the LEWICE 2D panel-based ice accretion code maintained by the NASA Lewis Research Center. The LEWICE/X model of ice accretion explicitly simulates regions of feather growth within the framework of the LEWICE model. Water droplets impinging on an accreting surface are withheld from the normal LEWICE mass/energy balance and handled in a separate routine; ice growth resulting from these droplets is performed with enhanced convective heat transfer approximately along droplet impingement directions. An independent underlying ice shape is grown along surface normals using the unmodified LEWICE method. The resulting dual-surface ice shape models roughness-induced feather growth observed in icing wind tunnel tests. Experiments indicate that the exact direction of feather growth is dependent on external conditions. Data is presented to support a linear variation of growth direction with temperature and cloud water content. Test runs of LEWICE/X indicate that the sizes of surface regions containing feathers are influenced by initial roughness element height. This suggests that a previous argument that feather region size is determined by boundary layer transition may be incorrect. Simulation results for two typical test cases give improved shape agreement over unmodified LEWICE.

The pre- and post-accretion irradiation history of cometary ices

NASA Technical Reports Server (NTRS)

Chyba, Christopher; Sagan, Carl

1989-01-01

Comets Halley and Wilson exhibited similar 3.4 micron emission features at approx. 1 AU from the Sun. A simple model of thermal emission from organic grains fits the feature, provides optical depths in good agreement with spacecraft measurements, and explains the absence of longer-wavelength organic features as due to spectral heliocentric evolution (Chyba and Sagan, 1987). The model utilizes transmission spectra of organics synthesized in the laboratory by irradiation of candidate cometary ices; the authors have long noted that related gas-phase syntheses yield polycyclic aromatic hydrocarbons, among other organic residues (Sagan et al., 1967). The authors previously concluded (Chyba and Sagan, 1987) that Halley's loss of several meters' depth with each perihelion passage, combined with the good fit of the Halley 3.4 micron feature to that of comet Wilson (Allen and Wickramasinghe, 1987), argues for the primordial - but not necessarily interstellar - origin of cometary organics. The authors examine the relative importance to the formation of organics of the variety of radiation environments experienced by comets. They conclude that there is at present no compelling reason to choose any of three contributing mechanisms (pre-accretion UV, pre-accretion cosmic ray, and post-accretion radionuclide processing) as the most important.

Recent Advances in the LEWICE Icing Model

NASA Technical Reports Server (NTRS)

Wright, William B.; Addy, Gene; Struk, Peter; Bartkus, Tadas

2015-01-01

This paper will describe two recent modifications to the Glenn ICE software. First, a capability for modeling ice crystals and mixed phase icing has been modified based on recent experimental data. Modifications have been made to the ice particle bouncing and erosion model. This capability has been added as part of a larger effort to model ice crystal ingestion in aircraft engines. Comparisons have been made to ice crystal ice accretions performed in the NRC Research Altitude Test Facility (RATFac). Second, modifications were made to the run back model based on data and observations from thermal scaling tests performed in the NRC Altitude Icing Tunnel.

Smart skin technology development for measuring ice accretion, stall, and high AOA aircraft performance. Part 1: Capacitive ice detector development

NASA Technical Reports Server (NTRS)

Pruzan, Daniel A.; Khatkhate, Ateen A.; Gerardi, Joseph J.; Hickman, Gail A.

1993-01-01

A reliable way to detect and measure ice accretion during flight is required to reduce the hazards of icing currently threatening present day aircraft. Many of the sensors used for this purpose are invasive (probe) sensors which must be placed in areas of the airframe where ice does not naturally form. Due to the difference in capture efficiency of the exposed surface, difficulties result in correlating the ice accretion on the probe to what is happening on a number of vastly different airfoil sections. Most flush mounted sensors in use must be integrated into the aircraft surface by cutting or drilling the aircraft surface. An alternate type of ice detector which is based on a NASA patent is currently being investigated at Innovative Dynamics, Inc. (IDI). Results of the investigation into the performance of different capacitive type sensor designs, both rigid as well as elastic, are presented.

Ice Accretion Calculations for a Commercial Transport Using the LEWICE3D, ICEGRID3D and CMARC Programs

NASA Technical Reports Server (NTRS)

Bidwell, Colin S.; Pinella, David; Garrison, Peter

1999-01-01

Collection efficiency and ice accretion calculations were made for a commercial transport using the NASA Lewis LEWICE3D ice accretion code, the ICEGRID3D grid code and the CMARC panel code. All of the calculations were made on a Windows 95 based personal computer. The ice accretion calculations were made for the nose, wing, horizontal tail and vertical tail surfaces. Ice shapes typifying those of a 30 minute hold were generated. Collection efficiencies were also generated for the entire aircraft using the newly developed unstructured collection efficiency method. The calculations highlight the flexibility and cost effectiveness of the LEWICE3D, ICEGRID3D, CMARC combination.

NASA Iced Aerodynamics and Controls Current Research

NASA Technical Reports Server (NTRS)

Addy, Gene

2009-01-01

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.

Structural health monitoring approach for detecting ice accretion on bridge cable using the Haar Wavelet Transform

NASA Astrophysics Data System (ADS)

Andre, Julia; Kiremidjian, Anne; Liao, Yizheng; Georgakis, Christos; Rajagopal, Ram

2016-04-01

Ice accretion on cables of bridge structures poses serious risk to the structure as well as to vehicular traffic when the ice falls onto the road. Detection of ice formation, quantification of the amount of ice accumulated, and prediction of icefalls will increase the safety and serviceability of the structure. In this paper, an ice accretion detection algorithm is presented based on the Continuous Wavelet Transform (CWT). In the proposed algorithm, the acceleration signals obtained from bridge cables are transformed using wavelet method. The damage sensitive features (DSFs) are defined as a function of the wavelet energy at specific wavelet scales. It is found that as ice accretes on the cables, the mass of cable increases, thus changing the wavelet energies. Hence, the DSFs can be used to track the change of cables mass. To validate the proposed algorithm, we use the data collected from a laboratory experiment conducted at the Technical University of Denmark (DTU). In this experiment, a cable was placed in a wind tunnel as ice volume grew progressively. Several accelerometers were installed at various locations along the testing cable to collect vibration signals.

Experimental study of performance degradation of a model helicopter main rotor with simulated ice shapes

NASA Technical Reports Server (NTRS)

Korkan, K. D.; Cross, E. J., Jr.; Cornell, C. C.

1984-01-01

An experimental study utilizing a remote controlled model helicopter has been conducted to measure the performance degradation due to simulated ice accretion on the leading edge of the main rotor for hover and forward flight. The 53.375 inch diameter main rotor incorporates a NACA 0012 airfoil with a generic ice shape corresponding to a specified natural ice condition. Thrust coefficients and torque coefficients about the main rotor were measured as a function of velocity, main rotor RPM, angle-of-incidence of the fuselage, collective pitch angle, and extent of spanwise ice accretion. An experimental airfoil data bank has been determined using a two-dimensional twenty-one inch NACA 0012 airfoil with scaled ice accretion shapes identical to that used on the model helicopter main rotor. The corresponding experimental data are discussed with emphasis on Reynolds number effects and ice accretion scale model testing.

Modeling of Commercial Turbofan Engine With Ice Crystal Ingestion: Follow-On

NASA Technical Reports Server (NTRS)

Jorgenson, Philip C. E.; Veres, Joseph P.; Coennen, Ryan

2014-01-01

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that have been attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was degraded engine performance, and one or more of the following: loss of thrust control (roll back), compressor surge or stall, and flameout of the combustor. As ice crystals are ingested into the fan and low pressure compression system, the increase in air temperature causes a portion of the ice crystals to melt. It is hypothesized that this allows the ice-water mixture to cover the metal surfaces of the compressor stationary components which leads to ice accretion through evaporative cooling. Ice accretion causes a blockage which subsequently results in the deterioration in performance of the compressor and engine. The focus of this research is to apply an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion. The tool is comprised of an engine system thermodynamic cycle code, a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor flow path, without modeling the actual ice accretion. A commercial turbofan engine which has previously experienced icing events during operation in a high altitude ice crystal environment has been tested in the Propulsion Systems Laboratory (PSL) altitude test facility at NASA Glenn Research Center. The PSL has the capability to produce a continuous ice cloud which is ingested by the engine during operation over a range of altitude conditions. The PSL test results confirmed that there was ice accretion in the engine due to ice crystal ingestion, at the same simulated altitude operating conditions as experienced previously in

Modeling of Commercial Turbofan Engine with Ice Crystal Ingestion; Follow-On

NASA Technical Reports Server (NTRS)

Jorgenson, Philip C. E.; Veres, Joseph P.; Coennen, Ryan

2014-01-01

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that have been attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was degraded engine performance, and one or more of the following: loss of thrust control (roll back), compressor surge or stall, and flameout of the combustor. As ice crystals are ingested into the fan and low pressure compression system, the increase in air temperature causes a portion of the ice crystals to melt. It is hypothesized that this allows the ice-water mixture to cover the metal surfaces of the compressor stationary components which leads to ice accretion through evaporative cooling. Ice accretion causes a blockage which subsequently results in the deterioration in performance of the compressor and engine. The focus of this research is to apply an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion. The tool is comprised of an engine system thermodynamic cycle code, a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor flow path, without modeling the actual ice accretion. A commercial turbofan engine which has previously experienced icing events during operation in a high altitude ice crystal environment has been tested in the Propulsion Systems Laboratory (PSL) altitude test facility at NASA Glenn Research Center. The PSL has the capability to produce a continuous ice cloud which is ingested by the engine during operation over a range of altitude conditions. The PSL test results confirmed that there was ice accretion in the engine due to ice crystal ingestion, at the same simulated altitude operating conditions as experienced previously in

'Scaling' analysis of the ice accretion process on aircraft surfaces

NASA Technical Reports Server (NTRS)

Keshock, E. G.; Tabrizi, A. H.; Missimer, J. R.

1982-01-01

A comprehensive set of scaling parameters is developed for the ice accretion process by analyzing the energy equations of the dynamic freezing zone and the already frozen ice layer, the continuity equation associated with supercooled liquid droplets entering into and impacting within the dynamic freezing zone, and energy equation of the ice layer. No initial arbitrary judgments are made regarding the relative magnitudes of each of the terms. The method of intrinsic reference variables in employed in order to develop the appropriate scaling parameters and their relative significance in rime icing conditions in an orderly process, rather than utilizing empiricism. The significance of these parameters is examined and the parameters are combined with scaling criteria related to droplet trajectory similitude.

User's manual for the NASA Lewis ice accretion/heat transfer prediction code with electrothermal deicer input

NASA Technical Reports Server (NTRS)

Masiulaniec, Konstanty C.; Wright, William B.

1994-01-01

A version of LEWICE has been developed that incorporates a recently developed electrothermal deicer code, developed at the University of Toledo by William B. Wright. This was accomplished, in essence, by replacing a subroutine in LEWICE, called EBAL, which balanced the energies at the ice surface, with a subroutine called UTICE. UTICE performs this same energy balance, as well as handles all the time-timperature transients below the ice surface, for all of the layers of a composite blade as well as the ice layer itself. This new addition is set up in such a fashion that a user may specify any number of heaters, any heater chordwise length, and any heater gap desired. The heaters may be fired in unison, or they may be cycled with periods independent of each other. The heater intensity may also be varied. In addition, the user may specify any number of layers and thicknesses depthwise into the blade. Thus, the new addition has maximum flexibility in modeling virtually any electrothermal deicer installed into any airfoil. It should be noted that the model simulates both shedding and runback. With the runback capability, it can simulate the anti-icing mode of heater performance, as well as detect icing downstream of the heaters due to runback in unprotected portions of the airfoil. This version of LEWICE can be run in three modes. In mode 1, no conduction heat transfer is modeled (which would be equivalent to the original version of LEWICE). In mode 2, all heat transfer is considered due to conduction but no heaters are firing. In mode 3, conduction heat transfer where the heaters are engaged is modeled, with subsequent ice shedding. When run in the first mode, there is virtually identical agreement with the original version of LEWICE in the prediction of accreted ice shapes. The code may be run in the second mode to determine the effects of conduction on the ice accretion process.

Experimental study of icing accretion on a rotating conical spinner

NASA Astrophysics Data System (ADS)

Chen, Ningli; Ji, Honghu; Hu, Yaping; Wang, Jian; Cao, Guangzhou

2015-12-01

A reduced scale experiment has been conducted to investigate the icing accretion procedure on a rotating spinner of 60° cone angle. The experiment was carried out in a small scale ice wind tunnel with three different rotating speeds of the spinner. The experimental conditions were determined from the actual icing condition of the spinner of a turbofan engine by using the similarity theory, which considers the rotating effects. The ice thickness on the spinner was got from the image taken by the high speed camera, by image processing. The results of this investigation show that under the experimental condition, ice on the spinner's tip of three different rotating speeds are all glaze ice and about the same thick. However, on the downstream surface of the spinner, ice shape on the rotating spinner is different from that on the stationary spinner. It is uneven glaze ice on the stationary spinner while it is `particle ice' when the rotating speed is 8240 rpm and it is `needle ice' when the rotating speed is 15,200 rpm. The experiment also reveals that when the rotating speed is higher, the ice layer is thicker.

A Preliminary Study of Ice-Accretion Scaling for SLD Conditions

NASA Technical Reports Server (NTRS)

Anderson, David N.

2003-01-01

Proposed changes to aircraft icing certification rules are being considered by European, Canadian, and American regulatory agencies to include operation in super-cooled large droplet conditions (SLD). This paper reports results of an experimental study in the NASA Glenn Icing Research Tunnel (IRT) to evaluate how well scaling methods developed for Appendix C conditions might apply to SLD conditions. Until now, scaling studies have been confined to the FAA FAR-25 Appendix C envelope of atmospheric cloud conditions. Tests were made in which it was attempted to scale to a droplet MVD of 50 microns from clouds having droplet MVDs of 175, 120, 100, and 70 microns. Scaling was based on the Ruff method with scale velocities found either by maintaining constant Weber number or by using the average of the velocities obtained by maintaining constant Weber number and constant Reynolds number. Models were unswept NACA 0012 wing sections. The reference model had a chord of 91.4 cm. Scale models had chords of 91.4, 80.0, and 53.3 cm. Tests were conducted with reference airspeeds of 100 and 150 kt (52 and 77 m/s) and with freezing fractions of 1.0, 0.6, and 0.3. It was demonstrated that the scaled 50-micron cloud simulated well the non-dimensional ice shapes accreted in clouds with MVD's of 120 microns or less.

Analytical determination of propeller performance degradation due to ice accretion

NASA Technical Reports Server (NTRS)

Miller, T. L.

1986-01-01

A computer code has been developed which is capable of computing propeller performance for clean, glaze, or rime iced propeller configurations, thereby providing a mechanism for determining the degree of performance degradation which results from a given icing encounter. The inviscid, incompressible flow field at each specified propeller radial location is first computed using the Theodorsen transformation method of conformal mapping. A droplet trajectory computation then calculates droplet impingement points and airfoil collection efficiency for each radial location, at which point several user-selectable empirical correlations are available for determining the aerodynamic penalities which arise due to the ice accretion. Propeller performance is finally computed using strip analysis for either the clean or iced propeller. In the iced mode, the differential 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 good agreement in several cases. The code's capability to properly predict iced thrust coefficient, power coefficient, and propeller efficiency is shown to be dependent on the choice of empirical correlation employed as well as proper specification of radial icing extent.

Users manual for the NASA Lewis three-dimensional ice accretion code (LEWICE 3D)

NASA Technical Reports Server (NTRS)

Bidwell, Colin S.; Potapczuk, Mark G.

1993-01-01

A description of the methodology, the algorithms, and the input and output data along with an example case for the NASA Lewis 3D ice accretion code (LEWICE3D) has been produced. The manual has been designed to help the user understand the capabilities, the methodologies, and the use of the code. The LEWICE3D code is a conglomeration of several codes for the purpose of calculating ice shapes on three-dimensional external surfaces. A three-dimensional external flow panel code is incorporated which has the capability of calculating flow about arbitrary 3D lifting and nonlifting bodies with external flow. A fourth order Runge-Kutta integration scheme is used to calculate arbitrary streamlines. An Adams type predictor-corrector trajectory integration scheme has been included to calculate arbitrary trajectories. Schemes for calculating tangent trajectories, collection efficiencies, and concentration factors for arbitrary regions of interest for single droplets or droplet distributions have been incorporated. A LEWICE 2D based heat transfer algorithm can be used to calculate ice accretions along surface streamlines. A geometry modification scheme is incorporated which calculates the new geometry based on the ice accretions generated at each section of interest. The three-dimensional ice accretion calculation is based on the LEWICE 2D calculation. Both codes calculate the flow, pressure distribution, and collection efficiency distribution along surface streamlines. For both codes the heat transfer calculation is divided into two regions, one above the stagnation point and one below the stagnation point, and solved for each region assuming a flat plate with pressure distribution. Water is assumed to follow the surface streamlines, hence starting at the stagnation zone any water that is not frozen out at a control volume is assumed to run back into the next control volume. After the amount of frozen water at each control volume has been calculated the geometry is modified by

LWC and Temperature Effects on Ice Accretion Formation on Swept Wings at Glaze Ice Conditions

NASA Technical Reports Server (NTRS)

Vargas, Mario; Reshotko, Eli

2000-01-01

An experiment was conducted to study the effect of liquid water content and temperature on the critical distance in ice accretion formation on swept wings at glaze ice conditions. The critical distance is defined as the distance from the attachment line to tile beginning of the zone where roughness elements develop into glaze ice feathers. A baseline case of 150 mph, 25 F, 0.75 g/cu m. Cloud Liquid Water Content (LWC) and 20 micrometers in Water Droplet Median Volume Diameter (MVD) was chosen. Icing runs were performed on a NACA 0012 swept wing tip at 150 mph and MVD of 20 micrometers for liquid water contents of 0.5 g/cu m, 0.75 g/cu m, and 1.0 g/cu m, and for total temperatures of 20 F, 25 F and 30 F. At each tunnel condition, the sweep angle was changed from 0 deg to 45 deg in 5 deg increments. Casting data, ice shape tracings, and close-up photographic data were obtained. The results showed that decreasing the LWC to 0.5 g/cu m decreases the value of the critical distance at a given sweep angle compared to the baseline case, and starts the formation of complete scallops at 30 sweep angle. Increasing the LWC to 1.0 g/cu m increases the value of the critical distance compared to the baseline case, the critical distance remains always above 0 millimeters and complete scallops are not formed. Decreasing the total temperature to 20 F decreases the critical distance with respect to the baseline case and formation of complete scallops begins at 25 deg sweep angle. When the total temperature is increased to 30 F, bumps covered with roughness elements appear on the ice accretion at 25 deg and 30 deg sweep angles, large ice structures appear at 35 deg and 40 deg sweep angles, and complete scallops are formed at 45 deg sweep angle.

Update to the NASA Lewis Ice Accretion Code LEWICE

NASA Technical Reports Server (NTRS)

Wright, William B.

1994-01-01

This report is intended as an update to NASA CR-185129 'User's Manual for the NASA Lewis Ice Accretion Prediction Code (LEWICE).' It describes modifications and improvements made to this code as well as changes to the input and output files, interactive input, and graphics output. The comparison of this code to experimental data is shown to have improved as a result of these modifications.

The Effects of the Critical Ice Accretion on Airfoil and Wing Performance

NASA Technical Reports Server (NTRS)

Selig, Michael S.; Bragg, Michael B.; Saeed, Farooq

1998-01-01

In support of the NASA Lewis Modern Airfoils Ice Accretion Test Program, the University of Illinois at Urbana-Champaign provided expertise in airfoil design and aerodynamic analysis to determine the aerodynamic effect of ice accretion on modern airfoil sections. The effort has concentrated on establishing a design/testing methodology for "hybrid airfoils" or "sub-scale airfoils," that is, airfoils having a full-scale leading edge together with a specially designed and foreshortened aft section. The basic approach of using a full-scale leading edge with a foreshortened aft section was considered to a limited extent over 40 years ago. However, it was believed that the range of application of the method had not been fully exploited. Thus a systematic study was being undertaken to investigate and explore the range of application of the method so as to determine its overall potential.

Ice stream motion facilitated by a shallow-deforming and accreting bed

PubMed Central

Spagnolo, Matteo; Phillips, Emrys; Piotrowski, Jan A.; Rea, Brice R.; Clark, Chris D.; Stokes, Chris R.; Carr, Simon J.; Ely, Jeremy C.; Ribolini, Adriano; Wysota, Wojciech; Szuman, Izabela

2016-01-01

Ice streams drain large portions of ice sheets and play a fundamental role in governing their response to atmospheric and oceanic forcing, with implications for sea-level change. The mechanisms that generate ice stream flow remain elusive. Basal sliding and/or bed deformation have been hypothesized, but ice stream beds are largely inaccessible. Here we present a comprehensive, multi-scale study of the internal structure of mega-scale glacial lineations (MSGLs) formed at the bed of a palaeo ice stream. Analyses were undertaken at macro- and microscales, using multiple techniques including X-ray tomography, thin sections and ground penetrating radar (GPR) acquisitions. Results reveal homogeneity in stratigraphy, kinematics, granulometry and petrography. The consistency of the physical and geological properties demonstrates a continuously accreting, shallow-deforming, bed and invariant basal conditions. This implies that ice stream basal motion on soft sediment beds during MSGL formation is accommodated by plastic deformation, facilitated by continuous sediment supply and an inefficient drainage system. PMID:26898399

Accretion growth of water-ice grains in astrophysically-relevant dusty plasma experiment

NASA Astrophysics Data System (ADS)

Chai, Kil-Byoung; Marshall, Ryan; Bellan, Paul

2016-10-01

The grain growth process in the Caltech water-ice dusty plasma experiment has been studied using a high-speed camera equipped with a long-distance microscope lens. It is found that (i) the ice grain number density decreases four-fold as the average grain length increases from 20 to 80 um, (ii) the ice grain length has a log-normal distribution rather than a power-law dependence, and (iii) no collisions between ice grains are apparent. The grains have a large negative charge so the agglomeration growth is prevented by their strong mutual repulsion. It is concluded that direct accretion of water molecules is in good agreement with the observed ice grain growth. The volumetric packing factor of the ice grains must be less than 0.25 in order for the grain kinetic energy to be sufficiently small to prevent collisions between ice grains; this conclusion is consistent with ice grain images showing a fractal character.

Heat Transfer Measurements on Surfaces with Natural Ice Castings and Modeled Roughness

NASA Technical Reports Server (NTRS)

Breuer, Kenneth S.; Torres, Benjamin E.; Orr, D. J.; Hansman, R. John

1997-01-01

An experimental method is described to measure and compare the convective heat transfer coefficient of natural and simulated ice accretion roughness and to provide a rational means for determining accretion-related enhanced heat transfer coefficients. The natural ice accretion roughness was a sample casting made from accretions at the NASA Lewis Icing Research Tunnel (IRT). One of these castings was modeled using a Spectral Estimation Technique (SET) to produce three roughness elements patterns that simulate the actual accretion. All four samples were tested in a flat-plate boundary layer at angle of attack in a "dry" wind tunnel test. The convective heat transfer coefficient was measured using infrared thermography. It is shown that, dispite some problems in the current data set, the method does show considerable promise in determining roughness-induced heat transfer coefficients, and that, in addition to the roughness height and spacing in the flow direction, the concentration and spacing of elements in the spanwise direction are important parameters.

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

NASA Technical Reports Server (NTRS)

Veres, Joseph P.; Jorgenson, Philip, C. E.; Jones, Scott M.

2014-01-01

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

An Initial Study of the Fundamentals of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory

NASA Technical Reports Server (NTRS)

Struk, Peter; Bartkus, Tadas; Tsao, Jen-Ching; Bencic, Timothy; King, Michael; Ratvasky, Thomas; Van Zante, Judith

2017-01-01

This presentation shows results from an initial study of the fundamental physics of ice-crystal ice accretion using the NASA Propulsion Systems Lab (PSL). Ice accretion due to the ingestion of ice-crystals is being attributed to numerous jet-engine power-loss events. The NASA PSL is an altitude jet-engine test facility which has recently added a capability to inject ice particles into the flow. NASA is evaluating whether this facility, in addition to full-engine and motor-driven-rig tests, can be used for more fundamental ice-accretion studies that simulate the different mixed-phase icing conditions along the core flow passage of a turbo-fan engine compressor. The data from such fundamental accretion tests will be used to help develop and validate models of the accretion process. The present study utilized a NACA0012 airfoil. The mixed-phase conditions were generated by partially freezing the liquid-water droplets ejected from the spray bars. This presentation shows data regarding (1) the freeze out characteristics of the cloud, (2) changes in aerothermal conditions due to the presence of the cloud, and (3) the ice accretion characteristics observed on the airfoil model. The primary variable in this test was the PSL plenum humidity which was systematically varied for two duct-exit-plane velocities (85 and 135 ms) as well as two particle size clouds (15 and 50 m MVDi). The observed clouds ranged from fully glaciated to fully liquid, where the liquid clouds were at least partially supercooled. The air total temperature decreased at the test section when the cloud was activated due to evaporation. The ice accretions observed ranged from sharp arrow-like accretions, characteristic of ice-crystal erosion, to cases with double-horn shapes, characteristic of supercooled water accretions.

An Initial Study of the Fundamentals of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory

NASA Technical Reports Server (NTRS)

Struk, Peter M.; Ratvasky, Thomas P.; Bencic, Timothy J.; Van Zante, Judith F.; King, Michael C.; Tsao, Jen-Ching; Bartkus, Tadas P.

2017-01-01

This paper presents results from an initial study of the fundamental physics of ice-crystal ice accretion using the NASA Propulsion Systems Lab (PSL). Ice accretion due to the ingestion of ice-crystals is being attributed to numerous jet-engine power-loss events. The NASA PSL is an altitude jet-engine test facility which has recently added a capability to inject ice particles into the flow. NASA is evaluating whether this facility, in addition to full-engine and motor-driven-rig tests, can be used for more fundamental ice-accretion studies that simulate the different mixed-phase icing conditions along the core flow passage of a turbo-fan engine compressor. The data from such fundamental accretion tests will be used to help develop and validate models of the accretion process. The present study utilized a NACA0012 airfoil. The mixed-phase conditions were generated by partially freezing the liquid-water droplets ejected from the spray bars. This paper presents data regarding (1) the freeze out characteristics of the cloud, (2) changes in aerothermal conditions due to the presence of the cloud, and (3) the ice accretion characteristics observed on the airfoil model. The primary variable in this test was the PSL plenum humidity which was systematically varied for two duct-exit-plane velocities (85 and 135 ms) as well as two particle size clouds (15 and 50 m MVDi). The observed clouds ranged from fully glaciated to fully liquid, where the liquid clouds were at least partially supercooled. The air total temperature decreased at the test section when the cloud was activated due to evaporation. The ice accretions observed ranged from sharp arrow-like accretions, characteristic of ice-crystal erosion, to cases with double-horn shapes, characteristic of supercooled water accretions.

Effect of High-Fidelity Ice Accretion Simulations on the Performance of a Full-Scale Airfoil Model

NASA Technical Reports Server (NTRS)

Broeren, Andy P.; Bragg, Michael B.; Addy, Harold E., Jr.; Lee, Sam; Moens, Frederic; Guffond, Didier

2010-01-01

The simulation of ice accretion on a wing or other surface is often required for aerodynamic evaluation, particularly at small scale or low-Reynolds number. While there are commonly accepted practices for ice simulation, there are no established and validated guidelines. The purpose of this article is to report the results of an experimental study establishing a high-fidelity, full-scale, iced-airfoil aerodynamic performance database. This research was conducted as a part of a larger program with the goal of developing subscale aerodynamic simulation methods for iced airfoils. Airfoil performance testing was carried out at the ONERA F1 pressurized wind tunnel using a 72-in. (1828.8-mm) chord NACA 23012 airfoil over a Reynolds number range of 4.5x10(exp 6) to 16.0 10(exp 6) and a Mach number range of 0.10 to 0.28. The high-fidelity, ice-casting simulations had a significant impact on the aerodynamic performance. A spanwise-ridge ice shape resulted in a maximum lift coefficient of 0.56 compared to the clean value of 1.85 at Re = 15.9x10(exp 6) and M = 0.20. Two roughness and streamwise shapes yielded maximum lift values in the range of 1.09 to 1.28, which was a relatively small variation compared to the differences in the ice geometry. The stalling characteristics of the two roughness and one streamwise ice simulation maintained the abrupt leading-edge stall type of the clean NACA 23012 airfoil, despite the significant decrease in maximum lift. Changes in Reynolds and Mach number over the large range tested had little effect on the iced-airfoil performance.

Modeling of Antarctic Sea Ice in a General Circulation Model.

NASA Astrophysics Data System (ADS)

Wu, Xingren; Simmonds, Ian; Budd, W. F.

1997-04-01

A dynamic-thermodynamic sea ice model is developed and coupled with the Melbourne University general circulation model to simulate the seasonal cycle of the Antarctic sea ice distribution. The model is efficient, rapid to compute, and useful for a range of climate studies. The thermodynamic part of the sea ice model is similar to that developed by Parkinson and Washington, the dynamics contain a simplified ice rheology that resists compression. The thermodynamics is based on energy conservation at the top surface of the ice/snow, the ice/water interface, and the open water area to determine the ice formation, accretion, and ablation. A lead parameterization is introduced with an effective partitioning scheme for freezing between and under the ice floes. The dynamic calculation determines the motion of ice, which is forced with the atmospheric wind, taking account of ice resistance and rafting. The simulated sea ice distribution compares reasonably well with observations. The seasonal cycle of ice extent is well simulated in phase as well as in magnitude. Simulated sea ice thickness and concentration are also in good agreement with observations over most regions and serve to indicate the importance of advection and ocean drift in the determination of the sea ice distribution.

Modeling of Antarctic sea ice in a general circulation model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wu, Xingren; Budd, W.F.; Simmonds, I.

1997-04-01

A dynamic-thermodynamic sea ice model is developed and coupled with the Melbourne University general circulation model to simulate the seasonal cycle of the Antarctic sea ice distributions The model is efficient, rapid to compute, and useful for a range of climate studies. The thermodynamic part of the sea ice model is similar to that developed by Parkinson and Washington, the dynamics contain a simplified ice rheology that resists compression. The thermodynamics is based on energy conservation at the top surface of the ice/snow, the ice/water interface, and the open water area to determine the ice formation, accretion, and ablation. Amore » lead parameterization is introduced with an effective partitioning scheme for freezing between and under the ice floes. The dynamic calculation determines the motion of ice, which is forced with the atmospheric wind, taking account of ice resistance and rafting. The simulated sea ice distribution compares reasonably well with observations. The seasonal cycle of ice extent is well simulated in phase as well as in magnitude. Simulated sea ice thickness and concentration are also in good agreement with observations over most regions and serve to indicate the importance of advection and ocean drift in the determination of the sea ice distribution. 64 refs., 15 figs., 2 tabs.« less

Method to Generate Full-Span Ice Shape on Swept Wing Using Icing Tunnel Data

NASA Technical Reports Server (NTRS)

Lee, Sam; Camello, Stephanie

2015-01-01

There is a collaborative research program by NASA, FAA, ONERA, and university partners to improve the fidelity of experimental and computational simulation methods for swept-wing ice accretion formulations and resultant aerodynamic effects on large transport aircraft. This research utilizes a 65 scale Common Research Model as the baseline configuration. In order to generate the ice shapes for the aerodynamic testing, ice-accretion testing will be conducted in the NASA Icing Research Tunnel utilizing hybrid model from the 20, 64, and 83 spanwise locations. The models will have full-scale leading edges with truncated chord in order to fit the IRT test section. The ice shapes from the IRT tests will be digitized using a commercially available articulated-arm 3D laser scanning system. The methodology to acquire 3D ice shapes using a laser scanner was developed and validated in a previous research effort. Each of these models will yield a 1.5ft span of ice than can be used. However, a full-span ice accretion will require 75 ft span of ice. This means there will be large gaps between these spanwise ice sections that must be filled, while maintaining all of the important aerodynamic features. A method was developed to generate a full-span ice shape from the three 1.5 ft span ice shapes from the three models.

Icing simulation: A survey of computer models and experimental facilities

NASA Technical Reports Server (NTRS)

Potapczuk, M. G.; Reinmann, J. J.

1991-01-01

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.

Icing simulation: A survey of computer models and experimental facilities

NASA Technical Reports Server (NTRS)

Potapczuk, M. G.; Reinmann, J. J.

1991-01-01

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.

Convection from Hemispherical and Conical Model Ice Roughness Elements in Stagnation Region Flows

NASA Technical Reports Server (NTRS)

Hughes, Michael T.; Shannon, Timothy A.; McClain, Stephen T.; Vargas, Mario; Broeren, Andy

2016-01-01

To improve ice accretion prediction codes, more data regarding ice roughness and its effects on convective heat transfer are required. The Vertical Icing Studies Tunnel (VIST) at NASA Glenn Research was used to model realistic ice roughness in the stagnation region of a NACA 0012 airfoil. In the VIST, a test plate representing the leading 2% chord of the airfoil was subjected to flows of 7.62 m/s (25 ft/s), 12.19 m/s (40 ft/s), and 16.76 m/s (55 ft/s). The test plate was fitted with multiple surfaces or sets of roughness panels, each with a different representation of ice roughness. The sets of roughness panels were constructed using two element distribution patterns that were created based on a laser scan of an iced airfoil acquired in the Icing Research Tunnel at NASA Glenn. For both roughness patterns, surfaces were constructed using plastic hemispherical elements, plastic conical elements, and aluminum conical elements. Infrared surface thermometry data from tests run in the VIST were used to calculate area averaged heat transfer coefficient values. The values from the roughness surfaces were compared to the smooth control surface, showing convective enhancement as high as 400% in some cases. The data gathered during this study will ultimately be used to improve the physical modeling in LEWICE or other ice accretion codes and produce predictions of in-flight ice accretion on aircraft surfaces with greater confidence.

Estimate of temperature change due to ice and snow accretion in the boreal forest regions

NASA Astrophysics Data System (ADS)

Sugiura, K.; Nagai, S.; Suzuki, R.; Eicken, H.; Maximov, T. C.

2016-12-01

Previous research has demonstrated that there is a wide difference between the surface albedo in winter/spring in snow-covered forest regions in various global climate models. If the forest is covered with snow, the surface albedo would increase. In this study, we carried out field observations to monitor the frequency of ice and snow accretion in the boreal forest regions. The time-lapse digital camera was set up on each side of the observation towers at the site located to the north of Fairbanks (USA) and at the site located to the north of Yakutsk (Russia). It was confirmed that both forests were not necessarily covered with snow without a break from the start of continuous snow cover until the end. In addition, the boreal forest at the Yakutsk site is covered with snow in comparison with the boreal forest at the Fairbanks site for a long term such as for about five month. Using a one-dimensional mathematics model about the energy flow including atmospheric multiple scattering, we estimated temperature change due to ice and snow accretion in the boreal forest regions. The result show that the mean surface temperature rises approximately 0.5 [oC] when the boreal forest is not covered with snow. In this presentation, we discuss the snow albedo parameterization in the boreal forest regions and the one-dimensional mathematics model to provide a basis for a better understanding of the role of snow in the climate system.

Subglacial Lake Vostok (Antarctica) Accretion Ice Contains a Diverse Set of Sequences from Aquatic, Marine and Sediment-Inhabiting Bacteria and Eukarya

PubMed Central

Edgar, Robyn; Veerapaneni, Ram S.; D’Elia, Tom; Morris, Paul F.; Rogers, Scott O.

2013-01-01

Lake Vostok, the 7th largest (by volume) and 4th deepest lake on Earth, is covered by more than 3,700 m of ice, making it the largest subglacial lake known. The combination of cold, heat (from possible hydrothermal activity), pressure (from the overriding glacier), limited nutrients and complete darkness presents extreme challenges to life. Here, we report metagenomic/metatranscriptomic sequence analyses from four accretion ice sections from the Vostok 5G ice core. Two sections accreted in the vicinity of an embayment on the southwestern end of the lake, and the other two represented part of the southern main basin. We obtained 3,507 unique gene sequences from concentrates of 500 ml of 0.22 µm-filtered accretion ice meltwater. Taxonomic classifications (to genus and/or species) were possible for 1,623 of the sequences. Species determinations in combination with mRNA gene sequence results allowed deduction of the metabolic pathways represented in the accretion ice and, by extension, in the lake. Approximately 94% of the sequences were from Bacteria and 6% were from Eukarya. Only two sequences were from Archaea. In general, the taxa were similar to organisms previously described from lakes, brackish water, marine environments, soil, glaciers, ice, lake sediments, deep-sea sediments, deep-sea thermal vents, animals and plants. Sequences from aerobic, anaerobic, psychrophilic, thermophilic, halophilic, alkaliphilic, acidophilic, desiccation-resistant, autotrophic and heterotrophic organisms were present, including a number from multicellular eukaryotes. PMID:23843994

Subglacial Lake Vostok (Antarctica) accretion ice contains a diverse set of sequences from aquatic, marine and sediment-inhabiting bacteria and eukarya.

PubMed

Shtarkman, Yury M; Koçer, Zeynep A; Edgar, Robyn; Veerapaneni, Ram S; D'Elia, Tom; Morris, Paul F; Rogers, Scott O

2013-01-01

Lake Vostok, the 7(th) largest (by volume) and 4(th) deepest lake on Earth, is covered by more than 3,700 m of ice, making it the largest subglacial lake known. The combination of cold, heat (from possible hydrothermal activity), pressure (from the overriding glacier), limited nutrients and complete darkness presents extreme challenges to life. Here, we report metagenomic/metatranscriptomic sequence analyses from four accretion ice sections from the Vostok 5G ice core. Two sections accreted in the vicinity of an embayment on the southwestern end of the lake, and the other two represented part of the southern main basin. We obtained 3,507 unique gene sequences from concentrates of 500 ml of 0.22 µm-filtered accretion ice meltwater. Taxonomic classifications (to genus and/or species) were possible for 1,623 of the sequences. Species determinations in combination with mRNA gene sequence results allowed deduction of the metabolic pathways represented in the accretion ice and, by extension, in the lake. Approximately 94% of the sequences were from Bacteria and 6% were from Eukarya. Only two sequences were from Archaea. In general, the taxa were similar to organisms previously described from lakes, brackish water, marine environments, soil, glaciers, ice, lake sediments, deep-sea sediments, deep-sea thermal vents, animals and plants. Sequences from aerobic, anaerobic, psychrophilic, thermophilic, halophilic, alkaliphilic, acidophilic, desiccation-resistant, autotrophic and heterotrophic organisms were present, including a number from multicellular eukaryotes.

Characterization of Ice Roughness From Simulated Icing Encounters

NASA Technical Reports Server (NTRS)

Anderson, David N.; Shin, Jaiwon

1997-01-01

Detailed measurements of the size of roughness elements on ice accreted on models in the NASA Lewis Icing Research Tunnel (IRT) were made in a previous study. Only limited data from that study have been published, but included were the roughness element height, diameter and spacing. In the present study, the height and spacing data were found to correlate with the element diameter, and the diameter was found to be a function primarily of the non-dimensional parameters freezing fraction and accumulation parameter. The width of the smooth zone which forms at the leading edge of the model was found to decrease with increasing accumulation parameter. Although preliminary, the success of these correlations suggests that it may be possible to develop simple relationships between ice roughness and icing conditions for use in ice-accretion-prediction codes. These codes now require an ice-roughness estimate to determine convective heat transfer. Studies using a 7.6-cm-diameter cylinder and a 53.3-cm-chord NACA 0012 airfoil were also performed in which a 1/2-min icing spray at an initial set of conditions was followed by a 9-1/2-min spray at a second set of conditions. The resulting ice shape was compared with that from a full 10-min spray at the second set of conditions. The initial ice accumulation appeared to have no effect on the final ice shape. From this result, it would appear the accreting ice is affected very little by the initial roughness or shape features.

Academic Airframe Icing Perspective

NASA Technical Reports Server (NTRS)

Bragg, Mike; Rothmayer, Alric; Thompson, David

2009-01-01

2-D ice accretion and aerodynamics reasonably well understood for engineering applications To significantly improve our current capabilities we need to understand 3-D: a) Important ice accretion physics and modeling not well understood in 3-D; and b) Aerodynamics unsteady and 3-D especially near stall. Larger systems issues important and require multidisciplinary team approach

Flight Testing an Iced Business Jet for Flight Simulation Model Validation

NASA Technical Reports Server (NTRS)

Ratvasky, Thomas P.; Barnhart, Billy P.; Lee, Sam; Cooper, Jon

2007-01-01

A flight test of a business jet aircraft with various ice accretions was performed to obtain data to validate flight simulation models developed through wind tunnel tests. Three types of ice accretions were tested: pre-activation roughness, runback shapes that form downstream of the thermal wing ice protection system, and a wing ice protection system failure shape. The high fidelity flight simulation models of this business jet aircraft were validated using a software tool called "Overdrive." Through comparisons of flight-extracted aerodynamic forces and moments to simulation-predicted forces and moments, the simulation models were successfully validated. Only minor adjustments in the simulation database were required to obtain adequate match, signifying the process used to develop the simulation models was successful. The simulation models were implemented in the NASA Ice Contamination Effects Flight Training Device (ICEFTD) to enable company pilots to evaluate flight characteristics of the simulation models. By and large, the pilots confirmed good similarities in the flight characteristics when compared to the real airplane. However, pilots noted pitch up tendencies at stall with the flaps extended that were not representative of the airplane and identified some differences in pilot forces. The elevator hinge moment model and implementation of the control forces on the ICEFTD were identified as a driver in the pitch ups and control force issues, and will be an area for future work.

Determination of Local Densities in Accreted Ice Samples Using X-Rays and Digital Imaging

NASA Technical Reports Server (NTRS)

Broughton, Howard; Sims, James; Vargas, Mario

1996-01-01

At the NASA Lewis Research Center's Icing Research Tunnel ice shapes, similar to those which develop in-flight icing conditions, were formed on an airfoil. Under cold room conditions these experimental samples were carefully removed from the airfoil, sliced into thin sections, and x-rayed. The resulting microradiographs were developed and the film digitized using a high resolution scanner to extract fine detail in the radiographs. A procedure was devised to calibrate the scanner and to maintain repeatability during the experiment. The techniques of image acquisition and analysis provide accurate local density measurements and reveal the internal characteristics of the accreted ice with greater detail. This paper will discuss the methodology by which these samples were prepared with emphasis on the digital imaging techniques.

Numerical investigation on super-cooled large droplet icing of fan rotor blade in jet engine

NASA Astrophysics Data System (ADS)

Isobe, Keisuke; Suzuki, Masaya; Yamamoto, Makoto

2014-10-01

Icing (or ice accretion) is a phenomenon in which super-cooled water droplets impinge and accrete on a body. It is well known that ice accretion on blades and vanes leads to performance degradation and has caused severe accidents. Although various anti-icing and deicing systems have been developed, such accidents still occur. Therefore, it is important to clarify the phenomenon of ice accretion on an aircraft and in a jet engine. However, flight tests for ice accretion are very expensive, and in the wind tunnel it is difficult to reproduce all climate conditions where ice accretion can occur. Therefore, it is expected that computational fluid dynamics (CFD), which can estimate ice accretion in various climate conditions, will be a useful way to predict and understand the ice accretion phenomenon. On the other hand, although the icing caused by super-cooled large droplets (SLD) is very dangerous, the numerical method has not been established yet. This is why SLD icing is characterized by splash and bounce phenomena of droplets and they are very complex in nature. In the present study, we develop an ice accretion code considering the splash and bounce phenomena to predict SLD icing, and the code is applied to a fan rotor blade. The numerical results with and without the SLD icing model are compared. Through this study, the influence of the SLD icing model is numerically clarified.

Proceedings of the Airframe Icing Workshop

NASA Technical Reports Server (NTRS)

Colantonio, Ron O. (Editor)

2009-01-01

The NASA Glenn Research Center (GRC) has a long history of working with its partners towards the understanding of ice accretion formation and its associated degradation of aerodynamic performance. The June 9, 2009, Airframe Icing Workshop held at GRC provided an opportunity to examine the current NASA airframe icing research program and to dialogue on remaining and emerging airframe icing issues and research with the external community. Some of the airframe icing gaps identified included, but are not limited to, ice accretion simulation enhancements, three-dimensional benchmark icing database development, three-dimensional iced aerodynamics modeling, and technology development for a smart icing system.

Atmospheric icing of structures: Observations and simulations

NASA Astrophysics Data System (ADS)

Ágústsson, H.; Elíasson, Á. J.; Thorsteins, E.; Rögnvaldsson, Ó.; Ólafsson, H.

2012-04-01

This study compares observed icing in a test span in complex orography at Hallormsstaðaháls (575 m) in East-Iceland with parameterized icing based on an icing model and dynamically downscaled weather at high horizontal resolution. Four icing events have been selected from an extensive dataset of observed atmospheric icing in Iceland. A total of 86 test-spans have been erected since 1972 at 56 locations in complex terrain with more than 1000 icing events documented. The events used here have peak observed ice load between 4 and 36 kg/m. Most of the ice accretion is in-cloud icing but it may partly be mixed with freezing drizzle and wet snow icing. The calculation of atmospheric icing is made in two steps. First the atmospheric data is created by dynamically downscaling the ECMWF-analysis to high resolution using the non-hydrostatic mesoscale Advanced Research WRF-model. The horizontal resolution of 9, 3, 1 and 0.33 km is necessary to allow the atmospheric model to reproduce correctly local weather in the complex terrain of Iceland. Secondly, the Makkonen-model is used to calculate the ice accretion rate on the conductors based on the simulated temperature, wind, cloud and precipitation variables from the atmospheric data. In general, the atmospheric model correctly simulates the atmospheric variables and icing calculations based on the atmospheric variables correctly identify the observed icing events, but underestimate the load due to too slow ice accretion. This is most obvious when the temperature is slightly below 0°C and the observed icing is most intense. The model results improve significantly when additional observations of weather from an upstream weather station are used to nudge the atmospheric model. However, the large variability in the simulated atmospheric variables results in high temporal and spatial variability in the calculated ice accretion. Furthermore, there is high sensitivity of the icing model to the droplet size and the possibility that

Modern Airfoil Ice Accretions

NASA Technical Reports Server (NTRS)

Addy, Harold E., Jr.; Potapczuk, Mark G.; Sheldon, David W.

1997-01-01

This report presents results from the first icing tests performed in the Modem Airfoils program. Two airfoils have been subjected to icing tests in the NASA Lewis Icing Research Tunnel (IRT). Both airfoils were two dimensional airfoils; one was representative of a commercial transport airfoil while the other was representative of a business jet airfoil. The icing test conditions were selected from the FAR Appendix C envelopes. Effects on aerodynamic performance are presented including the effects of varying amounts of glaze ice as well as the effects of approximately the same amounts of glaze, mixed, and rime ice. Actual ice shapes obtained in these tests are also presented for these cases. In addition, comparisons are shown between ice shapes from the tests and ice shapes predicted by the computer code, LEWICE for similar conditions. Significant results from the tests are that relatively small amounts of ice can have nearly as much effect on airfoil lift coefficient as much greater amounts of ice and that glaze ice usually has a more detrimental effect than either rime or mixed ice. LEWICE predictions of ice shapes, in general, compared reasonably well with ice shapes obtained in the IRT, although differences in details of the ice shapes were observed.

Plans and Preliminary Results of Fundamental Studies of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory

NASA Technical Reports Server (NTRS)

Struk, Peter; Tsao, Jen-Ching; Bartkus, Tadas

2017-01-01

This paper describes plans and preliminary results for using the NASA Propulsion Systems Lab (PSL) to experimentally study the fundamental physics of ice-crystal ice accretion. NASA is evaluating whether this facility, in addition to full-engine and motor-driven-rig tests, can be used for more fundamental ice-accretion studies that simulate the different mixed-phase icing conditions along the core flow passage of a turbo-fan engine compressor. The data from such fundamental accretion tests will be used to help develop and validate models of the accretion process. This paper presents data from some preliminary testing performed in May 2015 which examined how a mixed-phase cloud could be generated at PSL using evaporative cooling in a warmer-than-freezing environment.

Model helicopter performance degradation with simulated ice shapes

NASA Technical Reports Server (NTRS)

Tinetti, Ana F.; Korkan, Kenneth D.

1987-01-01

An experimental program using a commercially available model helicopter has been conducted in the Texas A&M University Subsonic Wind Tunnel to investigate main rotor performance degradation due to generic ice. The simulated ice, including both primary and secondary formations, was scaled by chord from previously documented artificial ice accretions. Base and iced performance data were gathered as functions of fuselage incidence, blade collective pitch, main rotor rotational velocity, and freestream velocity. It was observed that the presence of simulated ice tends to decrease the lift to equivalent drag ratio, as well as thrust coefficient for the range of velocity ratios tested. Also, increases in torque coefficient due to the generic ice formations were observed. Evaluation of the data has indicated that the addition of roughness due to secondary ice formations is crucial for proper evaluation of the degradation in main rotor performance.

Characterization of subglacial Lake Vostok as seen from physical and isotope properties of accreted ice.

PubMed

Lipenkov, Vladimir Ya; Ekaykin, Alexey A; Polyakova, Ekaterina V; Raynaud, Dominique

2016-01-28

Deep drilling at the Vostok Station has reached the surface of subglacial Lake Vostok (LV) twice-in February 2012 and January 2015. As a result, three replicate cores from boreholes 5G-1, 5G-2 and 5G-3 became available for detailed and revalidation analyses of the 230 m thickness of the accreted ice, down to its contact with water at 3769 m below the surface. The study reveals that the concentration of gases in the lake water beneath Vostok is unexpectedly low. A clear signature of the melt water in the surface layer of the lake, which is subject to refreezing on the icy ceiling of LV, has been discerned in the three different properties of the accreted ice: the ice texture, the isotopic and the gas content of the ice. These sets of data indicate in concert that poor mixing of the melt (and hydrothermal) water with the resident lake water and pronounced spatial and/or temporal variability of local hydrological conditions are likely to be the characteristics of the southern end of the lake. The latter implies that the surface water may be not representative enough to study LV's behaviour, and that direct sampling of the lake at different depths is needed in order to move ahead with our understanding of the lake's hydrological regime. © 2015 The Author(s).

Icing: Accretion, Detection, Protection

NASA Technical Reports Server (NTRS)

Reinmann, John J.

1994-01-01

The global aircraft industry and its regulatory agencies are currently involved in three major icing efforts: ground icing; advanced technologies for in-flight icing; and tailplane icing. These three major icing topics correspondingly support the three major segments of any aircraft flight profile: takeoff; cruise and hold; and approach and land. This lecture addressess these three topics in the same sequence as they appear in flight, starting with ground deicing, followed by advanced technologies for in-flight ice protection, and ending with tailplane icing.

Ice Accretion with Varying Surface Tension

NASA Technical Reports Server (NTRS)

Bilanin, Alan J.; Anderson, David N.

1995-01-01

During an icing encounter of an aircraft in flight, super-cooled water droplets impinging on an airfoil may splash before freezing. This paper reports tests performed to determine if this effect is significant and uses the results to develop an improved scaling method for use in icing test facilities. Simple laboratory tests showed that drops splash on impact at the Reynolds and Weber numbers typical of icing encounters. Further confirmation of droplet splash came from icing tests performed in the NaSA Lewis Icing Research Tunnel (IRT) with a surfactant added to the spray water to reduce the surface tension. The resulting ice shapes were significantly different from those formed when no surfactant was added to the water. These results suggested that the droplet Weber number must be kept constant to properly scale icing test conditions. Finally, the paper presents a Weber-number-based scaling method and reports results from scaling tests in the IRT in which model size was reduced up to a factor of 3. Scale and reference ice shapes are shown which confirm the effectiveness of this new scaling method.

Helium and Neon in the Accreted Ice of the Subglacial Antarctic Lake Vostok

NASA Astrophysics Data System (ADS)

Jean-Baptiste, P.; Fourré, E.; Petit, J. R.; Lipenkov, V.; Bulat, S.; Chetverikov, Y.; Raynaud, D.

2018-05-01

We analyzed helium and neon in 24 samples from between 3,607 and 3,767 m (i.e., down to 2 m above the lake-ice interface) of the accreted ice frozen to the ceiling of Lake Vostok. Within uncertainties, the neon budget of the lake is balanced, the neon supplied to the lake by the melting of glacier ice being compensated by the neon exported by lake ice. The helium concentration in the lake is about 12 times more than in the glacier ice, with a measured 3He/4He ratio of 0.12 ± 0.01 Ra. This shows that Lake Vostok's waters are enriched by a terrigenic helium source. The 3He/4He isotope ratio of this helium source was determined. Its radiogenic value (0.057 × Ra) is typical of an old continental province, ruling out any magmatic activity associated with the tectonic structure of the lake. It corresponds to a low geothermal heat flow estimated at 51 mW/m2.

Development and applications of a radar-attenuation model for polar ice sheets

NASA Astrophysics Data System (ADS)

MacGregor, Joseph A.

Modern ice sheets are currently responding to significant climatic forcings and undergoing ice-dynamics changes that are not yet well understood. Ice-penetrating radar surveys are often used to infer their basal condition (e.g., is the bed wet or dry?) and internal properties. However, such inferences typically require a model of the electromagnetic attenuation through the ice sheet. Here I first develop and test a radar-attenuation model that is based on a synthesis of existing laboratory measurements of the dielectric properties of ice. This synthesis shows that radar attenuation in polar ice has a strong non-linear temperature dependence and a weaker linear dependence on the concentrations of acid and sea-salt chloride. This model was tested at Siple Dome, West Antarctica, using ice-core-chemistry and borehole-temperature data, and the model agreed well with an existing radar-attenuation measurement. I then use this model to investigate the nature of radar detection of accreted ice over Lake Vostok, East Antarctica. My analysis of ice-core and radar data found that the observed reflection is likely due to a fabric contrast near the boundary between the dirty and clean accreted ices. This reflection mechanism is also consistent with the spatial pattern of detection of the reflection. In anticipation of the requirements of a thermomechanical ice-sheet model to predict the spatial variation of attenuation over Lake Vostok, I develop an accumulation-rate map for the Lake Vostok region using radar data, a steady-state flow-band model, and inverse methods. I found that accumulation rates there are not inversely correlated with surface elevation, that there is a broad maximum above the lake's northwestern corner, and a minimum above most of its eastern shoreline. Finally, I investigate the spatial variability of attenuation in an ice sheet, using the flowline that crosses through the Vostok ice core as an example. I use radar layers and ice-velocity and temperature

Steamworlds: Atmospheric Structure and Critical Mass of Planets Accreting Icy Pebbles

NASA Astrophysics Data System (ADS)

Chambers, John

2017-11-01

In the core accretion model, gas-giant planets first form a solid core, which then accretes gas from a protoplanetary disk when the core exceeds a critical mass. Here, we model the atmosphere of a core that grows by accreting ice-rich pebbles. The ice fraction of pebbles evaporates in warm regions of the atmosphere, saturating it with water vapor. Excess water precipitates to lower altitudes. Beneath an outer radiative region, the atmosphere is convective, following a moist adiabat in saturated regions due to water condensation and precipitation. Atmospheric mass, density, and temperature increase with core mass. For nominal model parameters, planets with core masses (ice + rock) between 0.08 and 0.16 Earth masses have surface temperatures between 273 and 647 K and form an ocean. In more massive planets, water exists as a supercritical convecting fluid mixed with gas from the disk. Typically, the core mass reaches a maximum (the critical mass) as a function of the total mass when the core is 2-5 Earth masses. The critical mass depends in a complicated way on pebble size, mass flux, and dust opacity due to the occasional appearance of multiple core-mass maxima. The core mass for an atmosphere of 50% hydrogen and helium may be a more robust indicator of the onset of gas accretion. This mass is typically 1-3 Earth masses for pebbles that are 50% ice by mass, increasing with opacity and pebble flux and decreasing with pebble ice/rock ratio.

Modeling of Highly Instrumented Honeywell Turbofan Engine Tested with Ice Crystal Ingestion in the NASA Propulsion System Laboratory

NASA Technical Reports Server (NTRS)

Veres, Joseph P.; Jorgenson, Philip C. E.; Jones, Scott M.

2016-01-01

The Propulsion Systems Laboratory (PSL), an altitude test facility at NASA Glenn Research Center, has been used to test a highly instrumented turbine engine at simulated altitude operating conditions. This is a continuation of the PSL testing that successfully duplicated the icing events that were experienced in a previous engine (serial LF01) during flight through ice crystal clouds, which was the first turbofan engine tested in PSL. This second model of the ALF502R-5A serial number LF11 is a highly instrumented version of the previous engine. The PSL facility provides a continuous cloud of ice crystals with controlled characteristics of size and concentration, which are ingested by the engine during operation at simulated altitudes. Several of the previous operating points tested in the LF01 engine were duplicated to confirm repeatability in LF11. The instrumentation included video cameras to visually illustrate the accretion of ice in the low pressure compressor (LPC) exit guide vane region in order to confirm the ice accretion, which was suspected during the testing of the LF01. Traditional instrumentation included static pressure taps in the low pressure compressor inner and outer flow path walls, as well as total pressure and temperature rakes in the low pressure compressor region. The test data was utilized to determine the losses and blockages due to accretion in the exit guide vane region of the LPC. Multiple data points were analyzed with the Honeywell Customer Deck. A full engine roll back point was modeled with the Numerical Propulsion System Simulation (NPSS) code. The mean line compressor flow analysis code with ice crystal modeling was utilized to estimate the parameters that indicate the risk of accretion, as well as to estimate the degree of blockage and losses caused by accretion during a full engine roll back point. The analysis provided additional validation of the icing risk parameters within the LPC, as well as the creation of models for

Identification of Accretion as Grain Growth Mechanism in Astrophysically Relevant Water&ice Dusty Plasma Experiment

NASA Astrophysics Data System (ADS)

Marshall, Ryan S.; Chai, Kil-Byoung; Bellan, Paul M.

2017-03-01

The grain growth process in the Caltech water-ice dusty plasma experiment has been studied using a high-speed camera and a long-distance microscope lens. It is observed that (I) the ice grain number density decreases fourfold as the average grain major axis increases from 20 to 80 μm, (II) the major axis length has a log-normal distribution rather than a power-law dependence, and (III) no collisions between ice grains are apparent. The grains have a large negative charge resulting in strong mutual repulsion and this, combined with the fractal character of the ice grains, prevents them from agglomerating. In order for the grain kinetic energy to be sufficiently small to prevent collisions between ice grains, the volumetric packing factor (I.e., ratio of the actual volume to the volume of a circumscribing ellipsoid) of the ice grains must be less than ˜0.1 depending on the exact relative velocity of the grains in question. Thus, it is concluded that direct accretion of water molecules is very likely to dominate the observed ice grain growth.

Steamworlds: Atmospheric Structure and Critical Mass of Planets Accreting Icy Pebbles

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chambers, John, E-mail: jchambers@carnegiescience.edu

In the core accretion model, gas-giant planets first form a solid core, which then accretes gas from a protoplanetary disk when the core exceeds a critical mass. Here, we model the atmosphere of a core that grows by accreting ice-rich pebbles. The ice fraction of pebbles evaporates in warm regions of the atmosphere, saturating it with water vapor. Excess water precipitates to lower altitudes. Beneath an outer radiative region, the atmosphere is convective, following a moist adiabat in saturated regions due to water condensation and precipitation. Atmospheric mass, density, and temperature increase with core mass. For nominal model parameters, planetsmore » with core masses (ice + rock) between 0.08 and 0.16 Earth masses have surface temperatures between 273 and 647 K and form an ocean. In more massive planets, water exists as a supercritical convecting fluid mixed with gas from the disk. Typically, the core mass reaches a maximum (the critical mass) as a function of the total mass when the core is 2–5 Earth masses. The critical mass depends in a complicated way on pebble size, mass flux, and dust opacity due to the occasional appearance of multiple core-mass maxima. The core mass for an atmosphere of 50% hydrogen and helium may be a more robust indicator of the onset of gas accretion. This mass is typically 1–3 Earth masses for pebbles that are 50% ice by mass, increasing with opacity and pebble flux and decreasing with pebble ice/rock ratio.« less

Ice Roughness in Short Duration SLD Icing Events

NASA Technical Reports Server (NTRS)

McClain, Stephen T.; Reed, Dana; Vargas, Mario; Kreeger, Richard E.; Tsao, Jen-Ching

2014-01-01

Ice accretion codes depend on models of roughness parameters to account for the enhanced heat transfer during the ice accretion process. While mitigating supercooled large droplet (SLD or Appendix O) icing is a significant concern for manufacturers seeking future vehicle certification due to the pending regulation, historical ice roughness studies have been performed using Appendix C icing clouds which exhibit mean volumetric diameters (MVD) much smaller than SLD clouds. Further, the historical studies of roughness focused on extracting parametric representations of ice roughness using multiple images of roughness elements. In this study, the ice roughness developed on a 21-in. NACA 0012 at 0deg angle of attack exposed to short duration SLD icing events was measured in the Icing Research Tunnel at the NASA Glenn Research Center. The MVD's used in the study ranged from 100 micrometer to 200 micrometers, in a 67 m/s flow, with liquid water contents of either 0.6 gm/cubic meters or 0.75 gm/cubic meters. The ice surfaces were measured using a Romer Absolute Arm laser scanning system. The roughness associated with each surface point cloud was measured using the two-dimensional self-organizing map approach developed by McClain and Kreeger (2013) resulting in statistical descriptions of the ice roughness.

Plans and Preliminary Results of Fundamental Studies of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory

NASA Technical Reports Server (NTRS)

Struk, Peter; Tsao, Jen-Ching; Bartkus, Tadas

2016-01-01

This presentation accompanies the paper titled Plans and Preliminary Results of Fundamental Studies of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory. NASA is evaluating whether PSL, in addition to full-engine and motor-driven-rig tests, can be used for more fundamental ice-accretion studies that simulate the different mixed-phase icing conditions along the core flow passage of a turbo-fan engine compressor. The data from such fundamental accretion tests will be used to help develop and validate models of the accretion process. This presentation (and accompanying paper) presents data from some preliminary testing performed in May 2015 which examined how a mixed-phase cloud could be generated at PSL using evaporative cooling in a warmer-than-freezing environment.

Scaling Methods for Simulating Aircraft In-Flight Icing Encounters

NASA Technical Reports Server (NTRS)

Anderson, David N.; Ruff, Gary A.

1997-01-01

This paper discusses scaling methods which permit the use of subscale models in icing wind tunnels to simulate natural flight in icing. Natural icing conditions exist when air temperatures are below freezing but cloud water droplets are super-cooled liquid. Aircraft flying through such clouds are susceptible to the accretion of ice on the leading edges of unprotected components such as wings, tailplane and engine inlets. To establish the aerodynamic penalties of such ice accretion and to determine what parts need to be protected from ice accretion (by heating, for example), extensive flight and wind-tunnel testing is necessary for new aircraft and components. Testing in icing tunnels is less expensive than flight testing, is safer, and permits better control of the test conditions. However, because of limitations on both model size and operating conditions in wind tunnels, it is often necessary to perform tests with either size or test conditions scaled. This paper describes the theoretical background to the development of icing scaling methods, discusses four methods, and presents results of tests to validate them.

Numerical simulations of icing in turbomachinery

NASA Astrophysics Data System (ADS)

Das, Kaushik

Safety concerns over aircraft icing and the high experimental cost of testing have spurred global interest in numerical simulations of the ice accretion process. Extensive experimental and computational studies have been carried out to understand the icing on external surfaces. No parallel initiatives were reported for icing on engine components. However, the supercooled water droplets in moist atmosphere that are ingested into the engine can impinge on the component surfaces and freeze to form ice deposits. Ice accretion could block the engine passage causing reduced airflow. It raises safety and performance concerns such as mechanical damage from ice shedding as well as slow acceleration leading to compressor stall. The current research aims at developing a computational methodology for prediction of icing phenomena on turbofan compression system. Numerical simulation of ice accretion in aircraft engines is highly challenging because of the complex 3-D unsteady turbomachinery flow and the effects of rotation on droplet trajectories. The aim of the present research focuses on (i) Developing a computational methodology for ice accretion in rotating turbomachinery components; (ii) Investigate the effect of inter-phase heat exchange; (iii) Characterize droplet impingement pattern and ice accretion at different operating conditions. The simulations of droplet trajectories are based on a Eulerian-Lagrangian approach for the continuous and discrete phases. The governing equations are solved in the rotating blade frame of reference. The flow field is computed by solving the 3-D solution of the compressible Reynolds Averaged Navier Stokes (RANS) equations. One-way interaction models simulate the effects of aerodynamic forces and the energy exchange between the flow and the droplets. The methodology is implemented in the cool, TURBODROP and applied to the flow field and droplet trajectories in NASA Roto-67r and NASA-GE E3 booster rotor. The results highlight the variation

Application of a High-Fidelity Icing Analysis Method to a Model-Scale Rotor in Forward Flight

NASA Technical Reports Server (NTRS)

Narducci, Robert; Orr, Stanley; Kreeger, Richard E.

2012-01-01

An icing analysis process involving the loose coupling of OVERFLOW-RCAS for rotor performance prediction and with LEWICE3D for thermal analysis and ice accretion is applied to a model-scale rotor for validation. The process offers high-fidelity rotor analysis for the noniced and iced rotor performance evaluation that accounts for the interaction of nonlinear aerodynamics with blade elastic deformations. Ice accumulation prediction also involves loosely coupled data exchanges between OVERFLOW and LEWICE3D to produce accurate ice shapes. Validation of the process uses data collected in the 1993 icing test involving Sikorsky's Powered Force Model. Non-iced and iced rotor performance predictions are compared to experimental measurements as are predicted ice shapes.

A comparison of selected models for estimating cable icing

NASA Astrophysics Data System (ADS)

McComber, Pierre; Druez, Jacques; Laflamme, Jean

In many cold climate countries, it is becoming increasingly important to monitor transmission line icing. Indeed, by knowing in advance of localized danger for icing overloads, electric utilities can take measures in time to prevent generalized failure of the power transmission network. Recently in Canada, a study was made to compare the estimation of a few icing models working from meteorological data in estimating ice loads for freezing rain events. The models tested were using only standard meteorological parameters, i.e. wind speed and direction, temperature and precipitation rate. This study has shown that standard meteorological parameters can only achieve very limited accuracy, especially for longer icing events. However, with the help of an additional instrument monitoring the icing rate intensity, a significant improvement in model prediction might be achieved. The icing rate meter (IRM) which counts icing and de-icing cycles per unit time on a standard probe can be used to estimate the icing intensity. A cable icing estimation is then made by taking into consideration the accretion size, temperature, wind speed and direction, and precipitation rate. In this paper, a comparison is made between the predictions of two previously tested models (one obtained and the other reconstructed from their description in the public literature) and of a model based on the icing rate meter readings. The models are tested against nineteen events recorded on an icing test line at Mt. Valin, Canada, during the winter season 1991-1992. These events are mostly rime resulting from in-cloud icing. However, freezing rain and wet snow events were also recorded. Results indicate that a significant improvement in the estimation is attained by using the icing rate meter data together with the other standard meteorological parameters.

A turbulence model for iced airfoils and its validation

NASA Technical Reports Server (NTRS)

Shin, Jaiwon; Chen, Hsun H.; Cebeci, Tuncer

1992-01-01

A turbulence model based on the extension of the algebraic eddy viscosity formulation of Cebeci and Smith developed for two dimensional flows over smooth and rough surfaces is described for iced airfoils and validated for computed ice shapes obtained for a range of total temperatures varying from 28 to -15 F. The validation is made with an interactive boundary layer method which uses a panel method to compute the inviscid flow and an inverse finite difference boundary layer method to compute the viscous flow. The interaction between inviscid and viscous flows is established by the use of the Hilbert integral. The calculated drag coefficients compare well with recent experimental data taken at the NASA-Lewis Icing Research Tunnel (IRT) and show that, in general, the drag increase due to ice accretion can be predicted well and efficiently.

Analytical ice shape predictions for flight in natural icing conditions

NASA Technical Reports Server (NTRS)

Berkowitz, Brian M.; Riley, James T.

1988-01-01

LEWICE is an analytical ice prediction code that has been evaluated against icing tunnel data, but on a more limited basis against flight data. Ice shapes predicted by LEWICE is compared with experimental ice shapes accreted on the NASA Lewis Icing Research Aircraft. The flight data selected for comparison includes liquid water content recorded using a hot wire device and droplet distribution data from a laser spectrometer; the ice shape is recorded using stereo photography. The main findings are as follows: (1) An equivalent sand grain roughness correlation different from that used for LEWICE tunnel comparisons must be employed to obtain satisfactory results for flight; (2) Using this correlation and making no other changes in the code, the comparisons to ice shapes accreted in flight are in general as good as the comparisons to ice shapes accreted in the tunnel (as in the case of tunnel ice shapes, agreement is least reliable for large glaze ice shapes at high angles of attack); (3) In some cases comparisons can be somewhat improved by utilizing the code so as to take account of the variation of parameters such as liquid water content, which may vary significantly in flight.

Passive infrared ice detection for helicopter applications

NASA Technical Reports Server (NTRS)

Dershowitz, Adam L.; Hansman, R. John, Jr.

1990-01-01

A technique is proposed to remotely detect rotor icing on helicopters by using passive IR thermometry to detect the warming caused by latent heat release as supercooled water freezes. During icing, the ice accretion region will be warmer than the uniced trailing edge, resulting in a characteristic chordwise temperature profile. Preliminary tests were conducted on a static model in the NASA Icing Research Tunnel for a variety of wet (glaze) and dry (rime) ice conditions. The chordwise temperature profiles were confirmed by observation with an IR thermal video system and thermocouple observations. The IR observations were consistent with predictions of the LEWICE ice accretion code, which was used to extrapolate the observations to rotor icing conditions. Based on the static observations, the passive IR ice detection technique appears promising; however, further testing or rotating blades is required.

Preliminary Results From a Heavily Instrumented Engine Ice Crystal Icing Test in a Ground Based Altitude Test Facility

NASA Technical Reports Server (NTRS)

Flegel, Ashlie B.; Oliver, Michael J.

2016-01-01

Preliminary results from the Heavily Instrumented ALF503R-5 Engine test conducted in the NASA Glenn Research Center Propulsion Systems Laboratory will be discussed. The effects of ice crystal icing on a full scale engine is examined and documented. This model engine, serial number LF01, was used during the inaugural icing test in the PSL facility. The reduction of thrust (rollback) events experienced by this engine in flight were replicated in the facility. Limited instrumentation was used to detect icing. Metal temperature on the exit guide vanes and outer shroud and the load measurement were the only indicators of ice formation. The current study features a similar engine, serial number LF11, which is instrumented to characterize the cloud entering the engine, detect characterize ice accretion, and visualize the ice accretion in the region of interest.

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

NASA Technical Reports Server (NTRS)

Thompson, David S.; Soni, Bharat K.

2000-01-01

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.

Astronomical Ice: The Effects of Treating Ice as a Porous Media on the Dynamics and Evolution of Extraterrestrial Ice-Ocean Environments

NASA Astrophysics Data System (ADS)

Buffo, J.; Schmidt, B. E.

2015-12-01

With the prevalence of water and ice rich environments in the solar system, and likely the universe, becoming more apparent, understanding the evolutionary dynamics and physical processes of such locales is of great importance. Piqued interest arises from the understanding that the persistence of all known life depends on the presence of liquid water. As in situ investigation is currently infeasible, accurate numerical modeling is the best technique to demystify these environments. We will discuss an evolving model of ice-ocean interaction aimed at realistically describing the behavior of the ice-ocean interface by treating basal ice as a porous media, and its possible implications on the formation of astrobiological niches. Treating ice as a porous media drastically affects the thermodynamic properties it exhibits. Thus inclusion of this phenomenon is critical in accurately representing the dynamics and evolution of all ice-ocean environments. This model utilizes equations that describe the dynamics of sea ice when it is treated as a porous media (Hunke et. al. 2011), coupled with a basal melt and accretion model (Holland and Jenkins 1999). Combined, these two models produce the most accurate description of the processes occurring at the base of terrestrial sea ice and ice shelves, capable of resolving variations within the ice due to environmental pressures. While these models were designed for application to terrestrial environments, the physics occurring at any ice-water interface is identical, and these models can be used to represent the evolution of a variety of icy astronomical bodies. As terrestrial ice shelves provide a close analog to planetary ice-ocean environments, we truth test the models validity against observations of ice shelves. We apply this model to the ice-ocean interface of the icy Galilean moon Europa. We include profiles of temperature, salinity, solid fraction, and Darcy velocity, as well as temporally and spatially varying melt and

Slush Fund: Modeling the Multiphase Physics of Oceanic Ices

NASA Astrophysics Data System (ADS)

Buffo, J.; Schmidt, B. E.

2016-12-01

The prevalence of ice interacting with an ocean, both on Earth and throughout the solar system, and its crucial role as the mediator of exchange between the hydrosphere below and atmosphere above, have made quantifying the thermodynamic, chemical, and physical properties of the ice highly desirable. While direct observations of these quantities exist, their scarcity increases with the difficulty of obtainment; the basal surfaces of terrestrial ice shelves remain largely unexplored and the icy interiors of moons like Europa and Enceladus have never been directly observed. Our understanding of these entities thus relies on numerical simulation, and the efficacy of their incorporation into larger systems models is dependent on the accuracy of these initial simulations. One characteristic of seawater, likely shared by the oceans of icy moons, is that it is a solution. As such, when it is frozen a majority of the solute is rejected from the forming ice, concentrating in interstitial pockets and channels, producing a two-component reactive porous media known as a mushy layer. The multiphase nature of this layer affects the evolution and dynamics of the overlying ice mass. Additionally ice can form in the water column and accrete onto the basal surface of these ice masses via buoyancy driven sedimentation as frazil or platelet ice. Numerical models hoping to accurately represent ice-ocean interactions should include the multiphase behavior of these two phenomena. While models of sea ice have begun to incorporate multiphase physics into their capabilities, no models of ice shelves/shells explicitly account for the two-phase behavior of the ice-ocean interface. Here we present a 1D multiphase model of floating oceanic ice that includes parameterizations of both density driven advection within the `mushy layer' and buoyancy driven sedimentation. The model is validated against contemporary sea ice models and observational data. Environmental stresses such as supercooling and

An Experimental Investigation on Bio-inspired Icephobic Coatings for Aircraft Icing Mitigation

NASA Astrophysics Data System (ADS)

Hu, Hui; Li, Haixing; Waldman, Rye

2016-11-01

By leveraging the Icing Research Tunnel available at Iowa State University (ISU-IRT), a series of experimental investigations were conducted to elucidate the underlying physics pertinent to aircraft icing phenomena. A suite of advanced flow diagnostic techniques, which include high-speed photographic imaging, digital image projection (DIP), and infrared (IR) imaging thermometry, were developed and applied to quantify the transient behavior of water droplet impingement, wind-driven surface water runback, unsteady heat transfer and dynamic ice accreting process over the surfaces of airfoil/wing models. The icephobic performance of various bio-inspired superhydrophobic coatings were evaluated quantitatively at different icing conditions. The findings derived from the icing physics studies can be used to improve current icing accretion models for more accurate prediction of ice formation and accretion on aircraft wings and to develop effective anti-/deicing strategies for safer and more efficient operation of aircraft in cold weather. The research work is partially supported by NASA with Grant Number NNX12AC21A and National Science Foundation under Award Numbers of CBET-1064196 and CBET-1435590.

Geometry and Reynolds-Number Scaling on an Iced Business-Jet Wing

NASA Technical Reports Server (NTRS)

Lee, Sam; Ratvasky, Thomas P.; Thacker, Michael; Barnhart, Billy P.

2005-01-01

A study was conducted to develop a method to scale the effect of ice accretion on a full-scale business jet wing model to a 1/12-scale model at greatly reduced Reynolds number. Full-scale, 5/12-scale, and 1/12-scale models of identical airfoil section were used in this study. Three types of ice accretion were studied: 22.5-minute ice protection system failure shape, 2-minute initial ice roughness, and a runback shape that forms downstream of a thermal anti-ice system. The results showed that the 22.5-minute failure shape could be scaled from full-scale to 1/12-scale through simple geometric scaling. The 2-minute roughness shape could be scaled by choosing an appropriate grit size. The runback ice shape exhibited greater Reynolds number effects and could not be scaled by simple geometric scaling of the ice shape.

A Revised Validation Process for Ice Accretion Codes

NASA Technical Reports Server (NTRS)

Wright, William B.; Porter, Christopher E.

2017-01-01

A research project is underway at NASA Glenn to produce computer software that can accurately predict ice growth under any meteorological conditions for any aircraft surface. This report will present results from the latest LEWICE release, version 3.5. This program differs from previous releases in its ability to model mixed phase and ice crystal conditions such as those encountered inside an engine. It also has expanded capability to use structured grids and a new capability to use results from unstructured grid flow solvers. A quantitative comparison of the results against a database of ice shapes that have been generated in the NASA Glenn Icing Research Tunnel (IRT) has also been performed. This paper will extend the comparison of ice shapes between LEWICE 3.5 and experimental data from a previous paper. Comparisons of lift and drag are made between experimentally collected data from experimentally obtained ice shapes and simulated (CFD) data on simulated (LEWICE) ice shapes. Comparisons are also made between experimentally collected and simulated performance data on select experimental ice shapes to ensure the CFD solver, FUN3D, is valid within the flight regime. The results show that the predicted results are within the accuracy limits of the experimental data for the majority of cases.

In-flight icing on unmanned aerial vehicle and its aerodynamic penalties

NASA Astrophysics Data System (ADS)

Szilder, K.; Yuan, W.

2017-06-01

A numerical prediction of ice accretion on HQ309, SD7032, and SD7037 airfoils and its aerodynamic penalties is described. Ice accretion prediction on a three-dimensional (3D) swept wing is also presented. In addition to airflow and drop trajectory solvers, NRC's (National Research Council) original, 3D, morphogenetic icing modeling approach has been used. The analysis was performed for a wide range of icing conditions identi¦ed in the FAA (Federal Aviation Administration) Appendix C icing envelope. They cover a range of drop sizes, air temperatures, and liquid water contents. For selected icing conditions, the resulting decrease in lift and increase in drag have been calculated.

Simulation numerique de l'accretion de glace sur une pale d'eolienne

NASA Astrophysics Data System (ADS)

Fernando, Villalpando

specifications, which are different from those that must be met in the wind energy industry. Among these differences are the following: wind turbines operate at subsonic speeds; the cords and angles of attack of wind turbine blades are smaller than those of aircraft wings; and a wind turbine can operate with a larger ice mass on its blades than an aircraft can. So, it is important to provide wind energy researchers with tools specifically validated with the operations parameters of a wind turbine. The main goal of this work is to develop a methodology to simulate ice accretion in 2D using Fluent and Matlab, commercial software programs that are available at nearly all research institutions. In this study, we used Gambit, previously the companion tool of Fluent, for mesh generation, and which has now been replaced by ICEM. We decided to stay with Gambit, because we were already deeply involved with the meshing procedure for our simulation of ice accretion at the time Gambit was removed from the market. We validate the methodology with experimental data consisting of iced airfoil contours obtained in a refrigerated wind tunnel using the parameters of actual ice conditions recorded in northern Quebec. This methodology consists of four steps: airfoil meshing, droplet trajectory calculation, thermodynamic model application, and airfoil contour updating. The total simulation time is divided into several time steps, for each of which the four steps are performed until the total time has elapsed. The time step length depends on the icing conditions. (Abstract shortened by UMI.).

Generation of Fullspan Leading-Edge 3D Ice Shapes for Swept-Wing Aerodynamic Testing

NASA Technical Reports Server (NTRS)

Camello, Stephanie C.; Lee, Sam; Lum, Christopher; Bragg, Michael B.

2016-01-01

The deleterious effect of ice accretion on aircraft is often assessed through dry-air flight and wind tunnel testing with artificial ice shapes. This paper describes a method to create fullspan swept-wing artificial ice shapes from partial span ice segments acquired in the NASA Glenn Icing Reserch Tunnel for aerodynamic wind-tunnel testing. Full-scale ice accretion segments were laser scanned from the Inboard, Midspan, and Outboard wing station models of the 65% scale Common Research Model (CRM65) aircraft configuration. These were interpolated and extrapolated using a weighted averaging method to generate fullspan ice shapes from the root to the tip of the CRM65 wing. The results showed that this interpolation method was able to preserve many of the highly three dimensional features typically found on swept-wing ice accretions. The interpolated fullspan ice shapes were then scaled to fit the leading edge of a 8.9% scale version of the CRM65 wing for aerodynamic wind-tunnel testing. Reduced fidelity versions of the fullspan ice shapes were also created where most of the local three-dimensional features were removed. The fullspan artificial ice shapes and the reduced fidelity versions were manufactured using stereolithography.

Dynamic Wind-Tunnel Testing of a Sub-Scale Iced Business Jet

NASA Technical Reports Server (NTRS)

Lee, Sam; Barnhart, Billy P.; Ratvasky, Thomas P.; Dickes, Edward; Thacker, Michael

2006-01-01

The effect of ice accretion on a 1/12-scale complete aircraft model of a business jet was studied in a rotary-balance wind tunnel. Three types of ice accretions were considered: ice protection system failure shape, pre-activation roughness, and runback shapes that form downstream of the thermal ice protection system. The results were compared with those from a 1/12-scale semi-span wing of the same aircraft at similar Reynolds number. The data showed that the full aircraft and the semi-span wing models showed similar characteristics, especially post stall behavior under iced configuration. However, there were also some discrepancies, such as the magnitude in the reductions in the maximum lift coefficient. Most of the ice-induced effects were limited to longitudinal forces. Rotational and forced oscillation studies showed that the effects of ice on lateral forces were relatively minor.

Design, fabrication, and evaluation of a partially melted ice particle cloud facility

NASA Astrophysics Data System (ADS)

Soltis, Jared T.

% melt to 3°C at 90%. There were two types of ice accretion with a transition zone in between. The first type of ice was opaque in color and had a rough surface. This ice occurred roughly from 6.0°C to 12.0°C duct temperatures (8% to 50% melt). The qualitative characteristics of the ice were produced from the low water content in the cloud. The water that was available froze instantly and trapped ice particle. Duct temperatures greater than 17.5°C (80% melt) produced ice that was clear and smooth. The water in the surface did not freeze instantly due to the high water content creating a water film that froze. A mixed-phase cloud dynamics model from NASA Glenn was used to estimate the percent melt of the cloud exiting the duct. There was no way to validate the model by directly measuring the percent melt of the cloud, so single particle melt experiments were conducted and compared to the model. A 0.05g/L solution of rhodamine b was sprayed into a levitator and droplets formed at the nodes of the wave. A 532nm green laser was used to illuminate the dye, and the water emitted orange 593nm light given the luminescent properties of the ink. The emitted light intensity was recorded, and a linear relationship between the light intensity of ice to the light intensity of water was used to determine the percent melt of a droplet. The droplets were frozen with a cold flow of nitrogen gas via a liquid nitrogen heat exchanger. The droplets melted under natural convection when the cold nitrogen was shut off. Fifteen cases were compared with droplet diameters ranging from 324mum to 1112mum, air temperatures from 16°C to 31°C, and relative humidities from 41% to 100%. The average discrepancy between predictions and results for the cases that melted slower than ten seconds was 13% while the cases that melted faster than 10 second had 64% discrepancy between the model and experiment. To explain the discrepancy between the experiment and model, sensitivity studies of the model were

Discrete Surface Evolution and Mesh Deformation for Aircraft Icing Applications

NASA Technical Reports Server (NTRS)

Thompson, David; Tong, Xiaoling; Arnoldus, Qiuhan; Collins, Eric; McLaurin, David; Luke, Edward; Bidwell, Colin S.

2013-01-01

Robust, automated mesh generation for problems with deforming geometries, such as ice accreting on aerodynamic surfaces, remains a challenging problem. Here we describe a technique to deform a discrete surface as it evolves due to the accretion of ice. The surface evolution algorithm is based on a smoothed, face-offsetting approach. We also describe a fast algebraic technique to propagate the computed surface deformations into the surrounding volume mesh while maintaining geometric mesh quality. Preliminary results presented here demonstrate the ecacy of the approach for a sphere with a prescribed accretion rate, a rime ice accretion, and a more complex glaze ice accretion.

Evaluation of constant-Weber-number scaling for icing tests

NASA Technical Reports Server (NTRS)

Anderson, David N.

1996-01-01

Previous studies showed that for conditions simulating an aircraft encountering super-cooled water droplets the droplets may splash before freezing. Other surface effects dependent on the water surface tension may also influence the ice accretion process. Consequently, the Weber number appears to be important in accurately scaling ice accretion. A scaling method which uses a constant-Weber-number approach has been described previously; this study provides an evaluation of this scaling method. Tests are reported on cylinders of 2.5 to 15-cm diameter and NACA 0012 airfoils with chords of 18 to 53 cm in the NASA Lewis Icing Research Tunnel (IRT). The larger models were used to establish reference ice shapes, the scaling method was applied to determine appropriate scaled test conditions using the smaller models, and the ice shapes were compared. Icing conditions included warm glaze, horn glaze and mixed. The smallest size scaling attempted was 1/3, and scale and reference ice shapes for both cylinders and airfoils indicated that the constant-Weber-number scaling method was effective for the conditions tested.

Wind tunnel tests of rotor blade sections with replications of ice formations accreted in hover

NASA Technical Reports Server (NTRS)

Lee, J. D.; Berger, J. H.; Mcdonald, T. J.

1986-01-01

Full scale reproductions of ice accretions molded during the documentation of a hover test program were fabricated by means of epoxy castings and used for a wind tunnel test program. Surface static pressure distributions were recorded and used to evaluate lift and pitching moment increments while drag was determined by wake surveys. Through the range of the tests, corresponding to those conditions encountered in hover and in flat pitch, integration of the pressure distributions showed negligible changes in lift and in pitching moment, but the drag was significantly increased.

An Explorative Study to Use DBD Plasma Generation for Aircraft Icing Mitigation

NASA Astrophysics Data System (ADS)

Hu, Hui; Zhou, Wenwu; Liu, Yang; Kolbakir, Cem

2017-11-01

An explorative investigation was performed to demonstrate the feasibility of utilizing thermal effect induced by Dielectric-Barrier-Discharge (DBD) plasma generation for aircraft icing mitigation. The experimental study was performed in an Icing Research Tunnel available at Iowa State University (i.e., ISU-IRT). A NACA0012 airfoil/wing model embedded with DBD plasma actuators was installed in ISU-IRT under typical glaze icing conditions pertinent to aircraft inflight icing phenomena. While a high-speed imaging system was used to record the dynamic ice accretion process over the airfoil surface for the test cases with and without switching on the DBD plasma actuators, an infrared (IR) thermal imaging system was utilized to map the corresponding temperature distributions to quantify the unsteady heat transfer and phase changing process over the airfoil surface. The thermal effect induced by DBD plasma generation was demonstrated to be able to keep the airfoil surface staying free of ice during the entire ice accretion experiment. The measured quantitative surface temperature distributions were correlated with the acquired images of the dynamic ice accretion and water runback processes to elucidate the underlying physics. National Science Foundation CBET-1064196 and CBET-1435590.

Aerodynamics of a finite wing with simulated ice

NASA Technical Reports Server (NTRS)

Bragg, M. B.; Khodadoust, A.; Kerho, M.

1992-01-01

The effect of a simulated glaze ice accretion on the aerodynamic performance of a three-dimensional wing is studied experimentally. Results are reviewed from earlier two-dimensional tests which show the character of the large leading-edge separation bubbles caused by the simulated ice accretion. The 2-D bubbles are found to closely resemble well known airfoil laminar separation bubbles. For the 3-D experiments a semispan wing of effective aspect ratio five was mounted from the sidewall of the UIUC subsonic wind tunnel. The model uses a NACA 0012 airfoil section on a rectangular planform with interchangeable tip and root sections to allow for 0- and 30-degree sweep. A three-component sidewall balance was used to measure lift, drag and pitching moment on the clean and iced model. Fluorescent oil flow visualization has been performed on the iced model and reveals extensive spanwise and vortical flow in the separation bubble aft of the upper surface horn. Sidewall interaction and spanwise nonuniformity are also seen on the unswept model. Comparisons to the computed flow fields are shown. Results are also shown for roughness effects on the straight wing. Sand grain roughness on the ice shape is seen to have a different effect than isolated 3-D roughness elements.

Dynamic Wind-Tunnel Testing of a Sub-Scale Iced S-3B Viking

NASA Technical Reports Server (NTRS)

Lee, Sam; Barnhart, Billy; Ratvasky, Thomas P.

2012-01-01

The effect of ice accretion on a 1/12-scale complete aircraft model of S-3B Viking was studied in a rotary-balance wind tunnel. Two types of ice accretions were considered: ice protection system failure shape and runback shapes that form downstream of the thermal ice protection system. The results showed that the ice shapes altered the stall characteristics of the aircraft. The ice shapes also reduced the control surface effectiveness, but mostly near the stall angle of attack. There were some discrepancies with the data with the flaps deflected that were attributed to the low Reynolds number of the test. Rotational and forced-oscillation studies showed that the effects of ice were mostly in the longitudinal forces, and the effects on the lateral forces were relatively minor.

Preliminary Results From a Heavily Instrumented Engine Ice Crystal Icing Test in a Ground Based Altitude Test Facility

NASA Technical Reports Server (NTRS)

Flegel, Ashlie B.; Oliver, Michael J.

2016-01-01

Preliminary results from the heavily instrumented ALF502R-5 engine test conducted in the NASA Glenn Research Center Propulsion Systems Laboratory are discussed. The effects of ice crystal icing on a full scale engine is examined and documented. This same model engine, serial number LF01, was used during the inaugural icing test in the Propulsion Systems Laboratory facility. The uncommanded reduction of thrust (rollback) events experienced by this engine in flight were simulated in the facility. Limited instrumentation was used to detect icing on the LF01 engine. Metal temperatures on the exit guide vanes and outer shroud and the load measurement were the only indicators of ice formation. The current study features a similar engine, serial number LF11, which is instrumented to characterize the cloud entering the engine, detect/ characterize ice accretion, and visualize the ice accretion in the region of interest. Data were acquired at key LF01 test points and additional points that explored: icing threshold regions, low altitude, high altitude, spinner heat effects, and the influence of varying the facility and engine parameters. For each condition of interest, data were obtained from some selected variations of ice particle median volumetric diameter, total water content, fan speed, and ambient temperature. For several cases the NASA in-house engine icing risk assessment code was used to find conditions that would lead to a rollback event. This study further helped NASA develop necessary icing diagnostic instrumentation, expand the capabilities of the Propulsion Systems Laboratory, and generate a dataset that will be used to develop and validate in-house icing prediction and risk mitigation computational tools. The ice accretion on the outer shroud region was acquired by internal video cameras. The heavily instrumented engine showed good repeatability of icing responses when compared to the key LF01 test points and during day-to-day operation. Other noticeable

Preliminary Results From a Heavily Instrumented Engine Ice Crystal Icing Test in a Ground Based Altitude Test Facility

NASA Technical Reports Server (NTRS)

Flegel, Ashlie B.; Oliver, Michael J.

2016-01-01

Preliminary results from the heavily instrumented ALF502R-5 engine test conducted in the NASA Glenn Research Center Propulsion Systems Laboratory are discussed. The effects of ice crystal icing on a full scale engine is examined and documented. This same model engine, serial number LF01, was used during the inaugural icing test in the Propulsion Systems Laboratory facility. The uncommanded reduction of thrust (rollback) events experienced by this engine in flight were simulated in the facility. Limited instrumentation was used to detect icing on the LF01 engine. Metal temperatures on the exit guide vanes and outer shroud and the load measurement were the only indicators of ice formation. The current study features a similar engine, serial number LF11, which is instrumented to characterize the cloud entering the engine, detect/characterize ice accretion, and visualize the ice accretion in the region of interest. Data were acquired at key LF01 test points and additional points that explored: icing threshold regions, low altitude, high altitude, spinner heat effects, and the influence of varying the facility and engine parameters. For each condition of interest, data were obtained from some selected variations of ice particle median volumetric diameter, total water content, fan speed, and ambient temperature. For several cases the NASA in-house engine icing risk assessment code was used to find conditions that would lead to a rollback event. This study further helped NASA develop necessary icing diagnostic instrumentation, expand the capabilities of the Propulsion Systems Laboratory, and generate a dataset that will be used to develop and validate in-house icing prediction and risk mitigation computational tools. The ice accretion on the outer shroud region was acquired by internal video cameras. The heavily instrumented engine showed good repeatability of icing responses when compared to the key LF01 test points and during day-to-day operation. Other noticeable

SmaggIce 2D Version 1.8: Software Toolkit Developed for Aerodynamic Simulation Over Iced Airfoils

NASA Technical Reports Server (NTRS)

Choo, Yung K.; Vickerman, Mary B.

2005-01-01

SmaggIce 2D version 1.8 is a software toolkit developed at the NASA Glenn Research Center that consists of tools for modeling the geometry of and generating the grids for clean and iced airfoils. Plans call for the completed SmaggIce 2D version 2.0 to streamline the entire aerodynamic simulation process--the characterization and modeling of ice shapes, grid generation, and flow simulation--and to be closely coupled with the public-domain application flow solver, WIND. Grid generated using version 1.8, however, can be used by other flow solvers. SmaggIce 2D will help researchers and engineers study the effects of ice accretion on airfoil performance, which is difficult to do with existing software tools because of complex ice shapes. Using SmaggIce 2D, when fully developed, to simulate flow over an iced airfoil will help to reduce the cost of performing flight and wind-tunnel tests for certifying aircraft in natural and simulated icing conditions.

A laser-based ice shape profilometer for use in icing wind tunnels

NASA Technical Reports Server (NTRS)

Hovenac, Edward A.; Vargas, Mario

1995-01-01

A laser-based profilometer was developed to measure the thickness and shape of ice accretions on the leading edge of airfoils and other models in icing wind tunnels. The instrument is a hand held device that is connected to a desk top computer with a 10 meter cable. It projects a laser line onto an ice shape and used solid state cameras to detect the light scattered by the ice. The instrument corrects the image for camera angle distortions, displays an outline of the ice shape on the computer screen, saves the data on a disk, and can print a full scale drawing of the ice shape. The profilometer has undergone extensive testing in the laboratory and in the NASA Lewis Icing Research Tunnel. Results of the tests show very good agreement between profilometer measurements and known simulated ice shapes and fair agreement between profilometer measurements and hand tracing techniques.

Titan's Cold Accretion and its Internal Structure

NASA Astrophysics Data System (ADS)

Estrada, Paul R.; Mosqueira, I.

2010-10-01

Recent Cassini radio tracking data has provided a normalized moment of inertia for Titan of 0.34 (Iess et al. 2010). Given that the quadrupole field is consistent with hydrostatic equilibrium, a two-layer interior model implies incomplete differentiation with a 700 km water-ice shell and an undifferentiated ice and rock-metal interior. We investigate the accretional history of Titan in connection with its internal structure. Our formation model allows for a size distribution of impactors with upper size cut-off constrained by Hyperion's size and a variable power-law exponent (Mosqueira et al. 2010). The burial of impact energy takes place in a lengthscale of order of the impactor radius, as indicated by numerical simulations (e.g., Pierazzo et al. 1997) applied to our energy regime of interest. Our thermal model includes radiogenic heating due to short and long-lived radionuclides, latent heat of melting, gravitational energy release due to sinking rock, heat of accretion and radiative cooling. We find that melting in the interior takes place well before the satellite reaches its final size. As a result, we expect the formation of an ocean overlying a silicate carapace, which may spend a considerable amount of time in contact with the liquid layer. Such a framework not only facilitates the transport of heat from the interior, but also can help both in leaching Ar40 into the ocean and then releasing into the atmosphere. We consider a range of parameters such as the degree of hydration of the rock component, the fraction of the impact energy that is deposited at the surface of the satellite, and accretion times. But we do not yet consider the effects of small admixtures of contaminants. We argue that models that form Titan in a cold environment may have allowed for the interior to remain cold enough as to preclude complete differentiation.

Can a Wind Model Mimic a Convection-Dominated Accretion Flow Model?

NASA Astrophysics Data System (ADS)

Chang, Heon-Young

2001-06-01

In this paper we investigate the properties of advection-dominated accretion flows(ADAFs) in case that outflows carry away infalling matter with its angular momentum and energy. Positive Bernoulli numbers in ADAFs allow a fraction of the gas to be ex-pelled in a form of outflows. The ADAFs are also unstable to convection. We present self-similar solutions for advection-dominated accretion flows in the presence of out-flows from the accretion flows (ADIOS). The axisymmetric flow is treated in variables integrated over polar sections and the effects of outflows on the accretion rlow are parameterized for possible configurations compatible with the one dimensional self-similar ADAF solution. We explicitly derive self-similar solutions of ADAFs in the presence of outflows and show that the strong outflows in the accretion flows result in a flatter density profile, which is similar to that of the convection-dominated accretion flows (CDAFs) in which convection transports the a! ngular momentum inward and the energy outward. These two different versions of the ADAF model should show similar behaviors in X-ray spectrum to some extent. Even though the two models may show similar behaviors, they should be distinguishable due to different physical properties. We suggest that for a central object of which mass is known these two different accretion flows should have different X-ray flux value due to deficient matter in the wind model.

Mass Fluxes of Ice and Oxygen Across the Entire Lid of Lake Vostok from Observations of Englacial Radiowave Attenuation

NASA Astrophysics Data System (ADS)

Winebrenner, D. P.; Kintner, P. M. S.; MacGregor, J. A.

2017-12-01

Over deep Antarctic subglacial lakes, spatially varying ice thickness and the pressure-dependent melting point of ice result in areas of melting and accretion at the ice-water interface, i.e., the lake lid. These ice mass fluxes drive lake circulation and, because basal Antarctic ice contains air-clathrate, affect the input of oxygen to the lake, with implications for subglacial life. Inferences of melting and accretion from radar-layer tracking and geodesy are limited in spatial coverage and resolution. Here we develop a new method to estimate rates of accretion, melting, and the resulting oxygen input at a lake lid, using airborne radar data over Lake Vostok together with ice-temperature and chemistry data from the Vostok ice core. Because the lake lid is a coherent reflector of known reflectivity (at our radar frequency), we can infer depth-averaged radiowave attenuation in the ice, with spatial resolution 1 km along flight lines. Spatial variation in attenuation depends mostly on variation in ice temperature near the lid, which in turn varies strongly with ice mass flux at the lid. We model ice temperature versus depth with ice mass flux as a parameter, thus linking that flux to (observed) depth-averaged attenuation. The resulting map of melt- and accretion-rates independently reproduces features known from earlier studies, but now covers the entire lid. We find that accretion is dominant when integrated over the lid, with an ice imbalance of 0.05 to 0.07 km3 a-1, which is robust against uncertainties.

Airframe Icing Research Gaps: NASA Perspective

NASA Technical Reports Server (NTRS)

Potapczuk, Mark

2009-01-01

qCurrent Airframe Icing Technology Gaps: Development of a full 3D ice accretion simulation model. Development of an improved simulation model for SLD conditions. CFD modeling of stall behavior for ice-contaminated wings/tails. Computational methods for simulation of stability and control parameters. Analysis of thermal ice protection system performance. Quantification of 3D ice shape geometric characteristics Development of accurate ground-based simulation of SLD conditions. Development of scaling methods for SLD conditions. Development of advanced diagnostic techniques for assessment of tunnel cloud conditions. Identification of critical ice shapes for aerodynamic performance degradation. Aerodynamic scaling issues associated with testing scale model ice shape geometries. Development of altitude scaling methods for thermal ice protections systems. Development of accurate parameter identification methods. Measurement of stability and control parameters for an ice-contaminated swept wing aircraft. Creation of control law modifications to prevent loss of control during icing encounters. 3D ice shape geometries. Collection efficiency data for ice shape geometries. SLD ice shape data, in-flight and ground-based, for simulation verification. Aerodynamic performance data for 3D geometries and various icing conditions. Stability and control parameter data for iced aircraft configurations. Thermal ice protection system data for simulation validation.

Overview of Icing Physics Relevant to Scaling

NASA Technical Reports Server (NTRS)

Anderson, David N.; Tsao, Jen-Ching

2005-01-01

An understanding of icing physics is required for the development of both scaling methods and ice-accretion prediction codes. This paper gives an overview of our present understanding of the important physical processes and the associated similarity parameters that determine the shape of Appendix C ice accretions. For many years it has been recognized that ice accretion processes depend on flow effects over the model, on droplet trajectories, on the rate of water collection and time of exposure, and, for glaze ice, on a heat balance. For scaling applications, equations describing these events have been based on analyses at the stagnation line of the model and have resulted in the identification of several non-dimensional similarity parameters. The parameters include the modified inertia parameter of the water drop, the accumulation parameter and the freezing fraction. Other parameters dealing with the leading edge heat balance have also been used for convenience. By equating scale expressions for these parameters to the values to be simulated a set of equations is produced which can be solved for the scale test conditions. Studies in the past few years have shown that at least one parameter in addition to those mentioned above is needed to describe surface-water effects, and some of the traditional parameters may not be as significant as once thought. Insight into the importance of each parameter, and the physical processes it represents, can be made by viewing whether ice shapes change, and the extent of the change, when each parameter is varied. Experimental evidence is presented to establish the importance of each of the traditionally used parameters and to identify the possible form of a new similarity parameter to be used for scaling.

Ice-Accretion Scaling Using Water-Film Thickness Parameters

NASA Technical Reports Server (NTRS)

Anderson, David N.; Feo, Alejandro

2003-01-01

Studies were performed at INTA in Spain to determine water-film thickness on a stagnation-point probe inserted in a simulated cloud. The measurements were correlated with non-dimensional parameters describing the flow and the cloud conditions. Icing scaling tests in the NASA Glenn Icing Research Tunnel were then conducted using the Ruff scaling method with the scale velocity found by matching scale and reference values of either the INTA non-dimensional water-film thickness or a Weber number based on that film thickness. For comparison, tests were also performed using the constant drop-size Weber number and the average-velocity methods. The reference and scale models were both aluminum, 61-cm-span, NACA 0012 airfoil sections at 0 deg. AOA. The reference had a 53-cm-chord and the scale, 27 cm (1/2 size). Both models were mounted vertically in the center of the IRT test section. Tests covered a freezing fraction range of 0.28 to 1.0. Rime ice (n = 1.0) tests showed the consistency of the IRT calibration over a range of velocities. At a freezing fraction of 0.76, there was no significant difference in the scale ice shapes produced by the different methods. For freezing fractions of 0.40, 0.52 and 0.61, somewhat better agreement with the reference horn angles was typically achieved with the average-velocity and constant-film thickness methods than when either of the two Weber numbers was matched to the reference value. At a freezing fraction of 0.28, the four methods were judged equal in providing simulations of the reference shape.

ICEG2D (v2.0) - An Integrated Software Package for Automated Prediction of Flow Fields for Single-Element Airfoils With Ice Accretion

NASA Technical Reports Server (NTRS)

Thompson David S.; Soni, Bharat K.

2001-01-01

An integrated geometry/grid/simulation software package, ICEG2D, is being developed to automate computational fluid dynamics (CFD) simulations for single- and multi-element airfoils with ice accretions. The current version, ICEG213 (v2.0), was designed to automatically perform four primary functions: (1) generate a grid-ready surface definition based on the geometrical characteristics of the iced airfoil surface, (2) generate high-quality structured and generalized grids starting from a defined surface definition, (3) generate the input and restart files needed to run the structured grid CFD solver NPARC or the generalized grid CFD solver HYBFL2D, and (4) using the flow solutions, generate solution-adaptive grids. ICEG2D (v2.0) can be operated in either a batch mode using a script file or in an interactive mode by entering directives from a command line within a Unix shell. This report summarizes activities completed in the first two years of a three-year research and development program to address automation issues related to CFD simulations for airfoils with ice accretions. As well as describing the technology employed in the software, this document serves as a users manual providing installation and operating instructions. An evaluation of the software is also presented.

Possible Mechanisms for Turbofan Engine Ice Crystal Icing at High Altitude

NASA Technical Reports Server (NTRS)

Tsao, Jen-Ching; Struk, Peter M.; Oliver, Michael

2014-01-01

A thermodynamic model is presented to describe possible mechanisms of ice formation on unheated surfaces inside a turbofan engine compression system from fully glaciated ice crystal clouds often formed at high altitude near deep convective weather systems. It is shown from the analysis that generally there could be two distinct types of ice formation: (1) when the "surface freezing fraction" is in the range of 0 to 1, dominated by the freezing of water melt from fully or partially melted ice crystals, the ice structure is formed from accretion with strong adhesion to the surface, and (2) when the "surface melting fraction" is the range of 0 to 1, dominated by the further melting of ice crystals, the ice structure is formed from accumulation of un-melted ice crystals with relatively weak bonding to the surface. The model captures important qualitative trends of the fundamental ice-crystal icing phenomenon reported earlier1,2 from the research collaboration work by NASA and the National Research Council (NRC) of Canada. Further, preliminary analysis of test data from the 2013 full scale turbofan engine ice crystal icing test3 conducted in the NASA Glenn Propulsion Systems Laboratory (PSL) has also suggested that (1) both types of ice formation occurred during the test, and (2) the model has captured some important qualitative trend of turning on (or off) the ice crystal ice formation process in the tested engine low pressure compressor (LPC) targeted area under different icing conditions that ultimately would lead to (or suppress) an engine core roll back (RB) event.

Possible Mechanisms for Turbofan Engine Ice Crystal Icing at High Altitude

NASA Technical Reports Server (NTRS)

Tsao, Jen-Ching; Struk, Peter M.; Oliver, Michael J.

2016-01-01

A thermodynamic model is presented to describe possible mechanisms of ice formation on unheated surfaces inside a turbofan engine compression system from fully glaciated ice crystal clouds often formed at high altitude near deep convective weather systems. It is shown from the analysis that generally there could be two distinct types of ice formation: (1) when the "surface freezing fraction" is in the range of 0 to 1, dominated by the freezing of water melt from fully or partially melted ice crystals, the ice structure is formed from accretion with strong adhesion to the surface, and (2) when the "surface melting fraction" is the range of 0 to 1, dominated by the further melting of ice crystals, the ice structure is formed from accumulation of un-melted ice crystals with relatively weak bonding to the surface. The model captures important qualitative trends of the fundamental ice-crystal icing phenomenon reported earlier (Refs. 1 and 2) from the research collaboration work by NASA and the National Research Council (NRC) of Canada. Further, preliminary analysis of test data from the 2013 full scale turbofan engine ice crystal icing test (Ref. 3) conducted in the NASA Glenn Propulsion Systems Laboratory (PSL) has also suggested that (1) both types of ice formation occurred during the test, and (2) the model has captured some important qualitative trend of turning on (or off) the ice crystal ice formation process in the tested engine low pressure compressor (LPC) targeted area under different icing conditions that ultimately would lead to (or suppress) an engine core roll back (RB) event.

The Dripping Handrail Model: Transient Chaos in Accretion Systems

NASA Technical Reports Server (NTRS)

Young, Karl; Scargle, Jeffrey D.; Cuzzi, Jeffrey (Technical Monitor)

1995-01-01

We define and study a simple dynamical model for accretion systems, the "dripping handrail" (DHR). The time evolution of this spatially extended system is a mixture of periodic and apparently random (but actually deterministic) behavior. The nature of this mixture depends on the values of its physical parameters - the accretion rate, diffusion coefficient, and density threshold. The aperiodic component is a special kind of deterministic chaos called transient chaos. The model can simultaneously exhibit both the quasiperiodic oscillations (QPO) and very low frequency noise (VLFN) that characterize the power spectra of fluctuations of several classes of accretion systems in astronomy. For this reason, our model may be relevant to many such astrophysical systems, including binary stars with accretion onto a compact object - white dwarf, neutron star, or black hole - as well as active galactic nuclei. We describe the systematics of the DHR's temporal behavior, by exploring its physical parameter space using several diagnostics: power spectra, wavelet "scalegrams," and Lyapunov exponents. In addition, we note that for large accretion rates the DHR has periodic modes; the effective pulse shapes for these modes - evaluated by folding the time series at the known period - bear a resemblance to the similarly- determined shapes for some x-ray pulsars. The pulsing observed in some of these systems may be such periodic-mode accretion, and not due to pure rotation as in the standard pulsar model.

Evaluation of Alternative Altitude Scaling Methods for Thermal Ice Protection System in NASA Icing Research Tunnel

NASA Technical Reports Server (NTRS)

Lee, Sam; Addy, Harold; Broeren, Andy P.; Orchard, David M.

2017-01-01

A test was conducted at NASA Icing Research Tunnel to evaluate altitude scaling methods for thermal ice protection system. Two scaling methods based on Weber number were compared against a method based on the Reynolds number. The results generally agreed with the previous set of tests conducted in NRCC Altitude Icing Wind Tunnel. The Weber number based scaling methods resulted in smaller runback ice mass than the Reynolds number based scaling method. The ice accretions from the Weber number based scaling method also formed farther upstream. However there were large differences in the accreted ice mass between the two Weber number based scaling methods. The difference became greater when the speed was increased. This indicated that there may be some Reynolds number effects that isnt fully accounted for and warrants further study.

Self-Consistent Models of Accretion Disks

NASA Technical Reports Server (NTRS)

Narayan, Ramesh

2000-01-01

Research was carried out on several topics in the theory of astrophysical accretion flows around black holes, neutron stars and white dwarfs. The focus of our effort was the advection-dominated accretion flow (ADAF) model which the PI and his collaborators proposed and developed over the last several years. Our group completed a total of 46 papers, of which 36 are in refereed journals and 12 are in conference proceedings. All the papers have either already appeared in print or are in press.

Reduced gas accretion on super-Earths and ice giants

NASA Astrophysics Data System (ADS)

Lambrechts, M.; Lega, E.

2017-10-01

A large fraction of giant planets have gaseous envelopes that are limited to about 10% of their total mass budget. Such planets are present in the solar system (Uranus, Neptune) and are frequently observed in short periods around other stars (the so-called super-Earths). In contrast to these observations, theoretical calculations based on the evolution of hydrostatic envelopes argue that such low-mass envelopes cannot be maintained around cores exceeding five Earth masses. Instead, under nominal disk conditions, these planets would acquire massive envelopes through runaway gas accretion within the lifetime of the protoplanetary disk. In this work we show that planetary envelopes are not in hydrostatic balance, which slows down envelope growth. A series of 3D global, radiative hydrodynamical simulations reveal a steady-state gas flow, which enters through the poles and exits in the disk midplane. Gas is pushed through the outer envelope in about ten orbital timescales. In regions of the disk that are not significantly dust-depleted, envelope accretion onto cores of about five Earth masses can get stalled as the gas flow enters the deep interior. Accreted solids sublimate deep in the convective interior, but small opacity-providing grains are trapped in the flow and do not settle, which further prevents rapid envelope accretion. The transition to runaway gas accretion can however be reached when cores grow larger than typical super-Earths, beyond 15 Earth masses, and preferably when disk opacities are below κ = 1 cm2/g. These findings offer an explanation for the typical low-mass envelopes around the cores of super-Earths.

Study on the glaze ice accretion of wind turbine with various chord lengths

NASA Astrophysics Data System (ADS)

Liang, Jian; Liu, Maolian; Wang, Ruiqi; Wang, Yuhang

2018-02-01

Wind turbine icing often occurs in winter, which changes the aerodynamic characteristics of the blades and reduces the work efficiency of the wind turbine. In this paper, the glaze ice model is established for horizontal-axis wind turbine in 3-D. The model contains the grid generation, two-phase simulation, heat and mass transfer. Results show that smaller wind turbine suffers from more serious icing problem, which reflects on a larger ice thickness. Both the collision efficiency and heat transfer coefficient increase under smaller size condition.

Growth of asteroids, planetary embryos, and Kuiper belt objects by chondrule accretion

PubMed Central

Johansen, Anders; Low, Mordecai-Mark Mac; Lacerda, Pedro; Bizzarro, Martin

2015-01-01

Chondrules are millimeter-sized spherules that dominate primitive meteorites (chondrites) originating from the asteroid belt. The incorporation of chondrules into asteroidal bodies must be an important step in planet formation, but the mechanism is not understood. We show that the main growth of asteroids can result from gas drag–assisted accretion of chondrules. The largest planetesimals of a population with a characteristic radius of 100 km undergo runaway accretion of chondrules within ~3 My, forming planetary embryos up to Mars’s size along with smaller asteroids whose size distribution matches that of main belt asteroids. The aerodynamical accretion leads to size sorting of chondrules consistent with chondrites. Accretion of millimeter-sized chondrules and ice particles drives the growth of planetesimals beyond the ice line as well, but the growth time increases above the disc lifetime outside of 25 AU. The contribution of direct planetesimal accretion to the growth of both asteroids and Kuiper belt objects is minor. In contrast, planetesimal accretion and chondrule accretion play more equal roles in the formation of Moon-sized embryos in the terrestrial planet formation region. These embryos are isolated from each other and accrete planetesimals only at a low rate. However, the continued accretion of chondrules destabilizes the oligarchic configuration and leads to the formation of Mars-sized embryos and terrestrial planets by a combination of direct chondrule accretion and giant impacts. PMID:26601169

Growth of asteroids, planetary embryos, and Kuiper belt objects by chondrule accretion

NASA Astrophysics Data System (ADS)

Johansen, Anders; Mac Low, Mordecai-Mark; Lacerda, Pedro; Bizzarro, Martin

2015-04-01

Chondrules are millimeter-sized spherules that dominate primitive meteorites (chondrites) originating from the asteroid belt. The incorporation of chondrules into asteroidal bodies must be an important step in planet formation, but the mechanism is not understood. We show that the main growth of asteroids can result from gas drag-assisted accretion of chondrules. The largest planetesimals of a population with a characteristic radius of 100 km undergo runaway accretion of chondrules within ~3 My, forming planetary embryos up to Mars's size along with smaller asteroids whose size distribution matches that of main belt asteroids. The aerodynamical accretion leads to size sorting of chondrules consistent with chondrites. Accretion of millimeter-sized chondrules and ice particles drives the growth of planetesimals beyond the ice line as well, but the growth time increases above the disc lifetime outside of 25 AU. The contribution of direct planetesimal accretion to the growth of both asteroids and Kuiper belt objects is minor. In contrast, planetesimal accretion and chondrule accretion play more equal roles in the formation of Moon-sized embryos in the terrestrial planet formation region. These embryos are isolated from each other and accrete planetesimals only at a low rate. However, the continued accretion of chondrules destabilizes the oligarchic configuration and leads to the formation of Mars-sized embryos and terrestrial planets by a combination of direct chondrule accretion and giant impacts.

Growth of asteroids, planetary embryos, and Kuiper belt objects by chondrule accretion.

PubMed

Johansen, Anders; Low, Mordecai-Mark Mac; Lacerda, Pedro; Bizzarro, Martin

2015-04-01

Chondrules are millimeter-sized spherules that dominate primitive meteorites (chondrites) originating from the asteroid belt. The incorporation of chondrules into asteroidal bodies must be an important step in planet formation, but the mechanism is not understood. We show that the main growth of asteroids can result from gas drag-assisted accretion of chondrules. The largest planetesimals of a population with a characteristic radius of 100 km undergo runaway accretion of chondrules within ~3 My, forming planetary embryos up to Mars's size along with smaller asteroids whose size distribution matches that of main belt asteroids. The aerodynamical accretion leads to size sorting of chondrules consistent with chondrites. Accretion of millimeter-sized chondrules and ice particles drives the growth of planetesimals beyond the ice line as well, but the growth time increases above the disc lifetime outside of 25 AU. The contribution of direct planetesimal accretion to the growth of both asteroids and Kuiper belt objects is minor. In contrast, planetesimal accretion and chondrule accretion play more equal roles in the formation of Moon-sized embryos in the terrestrial planet formation region. These embryos are isolated from each other and accrete planetesimals only at a low rate. However, the continued accretion of chondrules destabilizes the oligarchic configuration and leads to the formation of Mars-sized embryos and terrestrial planets by a combination of direct chondrule accretion and giant impacts.

An Overview of NASA Engine Ice-Crystal Icing Research

NASA Technical Reports Server (NTRS)

Addy, Harold E., Jr.; Veres, Joseph P.

2011-01-01

Ice accretions that have formed inside gas turbine engines as a result of flight in clouds of high concentrations of ice crystals in the atmosphere have recently been identified as an aviation safety hazard. NASA s Aviation Safety Program (AvSP) has made plans to conduct research in this area to address the hazard. This paper gives an overview of NASA s engine ice-crystal icing research project plans. Included are the rationale, approach, and details of various aspects of NASA s research.

Massive star formation by accretion. I. Disc accretion

NASA Astrophysics Data System (ADS)

Haemmerlé, L.; Eggenberger, P.; Meynet, G.; Maeder, A.; Charbonnel, C.

2016-01-01

Context. Massive stars likely form by accretion and the evolutionary track of an accreting forming star corresponds to what is called the birthline in the Hertzsprung-Russell (HR) diagram. The shape of this birthline is quite sensitive to the evolution of the entropy in the accreting star. Aims: We first study the reasons why some birthlines published in past years present different behaviours for a given accretion rate. We then revisit the question of the accretion rate, which allows us to understand the distribution of the observed pre-main-sequence (pre-MS) stars in the HR diagram. Finally, we identify the conditions needed to obtain a large inflation of the star along its pre-MS evolution that may push the birthline towards the Hayashi line in the upper part of the HR diagram. Methods: We present new pre-MS models including accretion at various rates and for different initial structures of the accreting core. We compare them with previously published equivalent models. From the observed upper envelope of pre-MS stars in the HR diagram, we deduce the accretion law that best matches the accretion history of most of the intermediate-mass stars. Results: In the numerical computation of the time derivative of the entropy, some treatment leads to an artificial loss of entropy and thus reduces the inflation that the accreting star undergoes along the birthline. In the case of cold disc accretion, the existence of a significant swelling during the accretion phase, which leads to radii ≳ 100 R⊙ and brings the star back to the red part of the HR diagram, depends sensitively on the initial conditions. For an accretion rate of 10-3M⊙ yr-1, only models starting from a core with a significant radiative region evolve back to the red part of the HR diagram. We also obtain that, in order to reproduce the observed upper envelope of pre-MS stars in the HR diagram with an accretion law deduced from the observed mass outflows in ultra-compact HII regions, the fraction of the

The Antartic Ice Borehole Probe

NASA Technical Reports Server (NTRS)

Behar, A.; Carsey, F.; Lane, A.; Engelhardt, H.

2000-01-01

The Antartic Ice Borehole Probe mission is a glaciological investigation, scheduled for November 2000-2001, that will place a probe in a hot-water drilled hole in the West Antartic ice sheet. The objectives of the probe are to observe ice-bed interactions with a downward looking camera, and ice inclusions and structure, including hypothesized ice accretion, with a side-looking camera.

Investigating and Modeling Ecosystem Response to an Experimental and a Natural Ice Storm

NASA Astrophysics Data System (ADS)

Fakhraei, H.; Driscoll, C. T.; Rustad, L.; Campbell, J. L.; Groffman, P.; Fahey, T.; Likens, G.; Swaminathan, R.

2017-12-01

Our understanding of ecosystem response to the extreme events is generally limited to rare observations from the natural historical events. However, investigating extreme events under controlled conditions can improve our understanding of these natural phenomena. A novel field experiment was conducted in a northern hardwood forest at the Hubbard Brook Experimental Forest in New Hampshire in the northeastern United States to quantify the influence of ice storms on the ecological processes. During subfreezing conditions in the winters of 2016 and 2017, water from a nearby stream was pumped and sprayed on the canopy of eight experimental plots to accrete ice to a targeted thickness on the canopy. The experiment was conducted at three levels of icing thickness (0.25, 0.5, 0.75 in.) in 2016 comparable to the naturally occurring 1998 ice storm and a second 0.5 in. treatment 2017 which were compared with reference plots. The most notable response of the icing treatments was a marked increase in fine and course litter fall which increased exponentially with increases in the icing thickness. Post-treatment openings in the canopy caused short-term increases in soil temperature in the ice-treatment plots compared to the reference plots. No response from the ice storm treatments were detected for soil moisture, net N mineralization, net nitrification, or denitrification after both natural and experimental ice storm. In contrast to the marked increase in the stream water nitrate after the natural occurring 1998 ice storm, we have not observed any significant change in soil solution N concentrations in the experimental ice storm treatments. Inconsistency in the response between the natural and experimental ice storm is likely due to differences in geophysical characteristics of the study sites including slope and lateral uptake of nutrient by the trees outside the experimental plots. In order to evaluate the long-term impacts of ice storms on northern hardwood forests, we used

Self-Consistent Models of Accretion Disks

NASA Technical Reports Server (NTRS)

Narayan, Ramesh

1997-01-01

The investigations of advection-dominated accretion flows (ADAFs), with emphasis on applications to X-ray binaries containing black holes and neutron stars is presented. This work is now being recognized as the standard paradigm for understanding the various spectral states of black hole X-ray Binaries (BHXBs). Topics discussed include: (1) Problem in BHXBS, namely that several of these binaries have unusually large concentrations of lithium in their companion stars; (2) A novel test to show that black holes have event horizons; (3) Application of the ADAF model to the puzzling X-ray delay in the recent outburst of the BHXB, GRO J1655-40; (4) Description of the various spectral states in BHXBS; (5) Application of the ADAF model to the famous supermassive black hole at the center of our Galaxy, Sgr A(*); (6) Writing down and solving equations describing steady-state, optically thin, advection-dominated accretion onto a Kerr black hole; (7) The effect of "photon bubble" instability on radiation dominated accretion disks; and (8) Dwarf nova disks in quiescence that have rather low magnetic Reynolds number, of order 10(exp 3).

Close-up analysis of aircraft ice accretion

NASA Technical Reports Server (NTRS)

Hansman, R. John; Breuer, Kenneth S.; Hazan, Didier; Reehorst, Andrew; Vargas, Mario

1993-01-01

Various types of ice formation have been studied by analysis of high magnification video observations. All testing was conducted in the NASA Lewis Icing Research Tunnel (IRT). A faired 8.9 cm (3.5 in.) diameter metal-clad cylinder and a 5.1 (2 in.) aluminum cylinder were observed by close-up and overview video cameras for several wind tunnel conditions. These included close-up grazing angle, close-up side view, as well as overhead and side overview cameras. Still photographs were taken at the end of each spray along with tracings of the subsequent ice shape. While in earlier tests only the stagnation region was observed, the entire area from the stagnation line to the horn region of glaze ice shapes was observed in this test. The modes or horn formation have been identified within the range of conditions observed. In the horn region, Horn Type A ice is formed by 'dry' feather growth into the flow direction and Horn Type B is formed by a 'wet' growth normal to the surface. The feather growth occurs when the freezing fraction is near unity and roughness elements exist to provide an initial growth site.

An Accretion Model for Anomalous X-Ray Pulsars

NASA Astrophysics Data System (ADS)

Chatterjee, Pinaki; Hernquist, Lars; Narayan, Ramesh

2000-05-01

We present a model for the anomalous X-ray pulsars (AXPs) in which the emission is powered by accretion from a fossil disk, established from matter falling back onto the neutron star following its birth. The time-dependent accretion drives the neutron star toward a ``tracking'' solution in which the rotation period of the star increases slowly, in tandem with the declining accretion rate. For appropriate choices of disk mass, neutron star magnetic field strength, and initial spin period, we demonstrate that a rapidly rotating neutron star can be spun down to periods characteristic of AXPs on timescales comparable to the estimated ages of these sources. In other cases, accretion onto the neutron star switches off after a short time and the star becomes an ordinary radio pulsar. Thus, in our picture, radio pulsars and AXPs are drawn from the same underlying population, in contrast to the situation in models involving neutron stars with ultrastrong magnetic fields, which require a new population of stars with very different properties.

Turbulent dispersion of the icing cloud from spray nozzles used in icing tunnels

NASA Technical Reports Server (NTRS)

Marek, C. J.; Olsen, W. A., Jr.

1986-01-01

To correctly simulate flight in natural icing conditions, the turbulence in an icing simulator must be as low as possible. But some turbulence is required to mix the droplets from the spray nozzles and achieve an icing cloud of uniform liquid water content. The goal for any spray system is to obtain the widest possible spray cloud with the lowest possible turbulence in the test section of a icing tunnel. This investigation reports the measurement of turbulence and the three-dimensional spread of the cloud from a single spray nozzle. The task was to determine how the air turbulence and cloud width are affected by spray bars of quite different drag coefficients, by changes in the turbulence upstream of the spray, the droplet size, and the atomizing air. An ice accretion grid, located 6.3 m downstream of the single spray nozzle, was used to measure cloud spread. Both the spray bar and the grid were located in the constant velocity test section. Three spray bar shapes were tested: the short blunt spray bar used in the NASA Lewis Icing Research Tunnel, a thin 14.6 cm chord airfoil, and a 53 cm chord NACA 0012 airfoil. At the low airspeed (56 km/hr) the ice accretion pattern was axisymmetric and was not affected by the shape of the spray bar. At the high airspeed (169 km/hr) the spread was 30 percent smaller than at the low airspeed. For the widest cloud the spray bars should be located as far upstream in the low velocity plenum of the icing tunnel. Good comparison is obtained between the cloud spread data and predicitons from a two-dimensional cloud mixing computer code using the two equation turbulence (k epsilon g) model.

NASA's program on icing research and technology

NASA Technical Reports Server (NTRS)

Reinmann, John J.; Shaw, Robert J.; Ranaudo, Richard J.

1989-01-01

NASA's program in aircraft icing research and technology is reviewed. The program relies heavily on computer codes and modern applied physics technology in seeking icing solutions on a finer scale than those offered in earlier programs. Three major goals of this program are to offer new approaches to ice protection, to improve our ability to model the response of an aircraft to an icing encounter, and to provide improved techniques and facilities for ground and flight testing. This paper reviews the following program elements: (1) new approaches to ice protection; (2) numerical codes for deicer analysis; (3) measurement and prediction of ice accretion and its effect on aircraft and aircraft components; (4) special wind tunnel test techniques for rotorcraft icing; (5) improvements of icing wind tunnels and research aircraft; (6) ground de-icing fluids used in winter operation; (7) fundamental studies in icing; and (8) droplet sizing instruments for icing clouds.

Aircraft Icing Handbook. Volume 2

DTIC Science & Technology

1991-03-01

an airfoil surface. icenhobig - A surface property exhibiting a reduced adhesion to ice; literally, "ice-hating." light icing - The rate of...power, and are a light weight system of reasonable cost. K. ill I-I1 1.I.2 Pneumatic Impulse Ice Protection A Pneumatic Impulse Ice Protection System...should be about 5 to 6 seconds. During moderate icing a 60 second cycle is suggested, while for light icing, longer accretion times of 3 to 4 minutes

Terrane accretion: Insights from numerical modelling

NASA Astrophysics Data System (ADS)

Vogt, Katharina; Gerya, Taras

2016-04-01

The oceanic crust is not homogenous, but contains significantly thicker crust than norm, i.e. extinct arcs, spreading ridges, detached continental fragments, volcanic piles or oceanic swells. These (crustal) fragments may collide with continental crust and form accretionary complexes, contributing to its growth. We analyse this process using a thermo-mechanical computer model (i2vis) of an ocean-continent subduction zone. In this model the oceanic plate can bend spontaneously under the control of visco-plastic rheologies. It moreover incorporates effects such as mineralogical phase changes, fluid release and consumption, partial melting and melt extraction. Based on our 2-D experiments we suggest that the lithospheric buoyancy of the downgoing slab and the rheological strength of crustal material may result in a variety of accretionary processes. In addition to terrane subduction, we are able to identify three distinct modes of terrane accretion: frontal accretion, basal accretion and underplating plateaus. We show that crustal fragments may dock onto continental crust and cease subduction, be scrapped off the downgoing plate, or subduct to greater depth prior to slab break off and subsequent exhumation. Direct consequences of these processes include slab break off, subduction zone transference, structural reworking, formation of high-pressure terranes, partial melting and crustal growth.

Simulation Tools Model Icing for Aircraft Design

NASA Technical Reports Server (NTRS)

2012-01-01

Here s a simple science experiment to try: Place an unopened bottle of distilled water in your freezer. After 2-3 hours, if the water is pure enough, you will notice that it has not frozen. Carefully pour the water into a bowl with a piece of ice in it. When it strikes the ice, the water will instantly freeze. One of the most basic and commonly known scientific facts is that water freezes at around 32 F. But this is not always the case. Water lacking any impurities for ice crystals to form around can be supercooled to even lower temperatures without freezing. High in the atmosphere, water droplets can achieve this delicate, supercooled state. When a plane flies through clouds containing these droplets, the water can strike the airframe and, like the supercooled water hitting the ice in the experiment above, freeze instantly. The ice buildup alters the aerodynamics of the plane - reducing lift and increasing drag - affecting its performance and presenting a safety issue if the plane can no longer fly effectively. In certain circumstances, ice can form inside aircraft engines, another potential hazard. NASA has long studied ways of detecting and countering atmospheric icing conditions as part of the Agency s efforts to enhance aviation safety. To do this, the Icing Branch at Glenn Research Center utilizes a number of world-class tools, including the Center s Icing Research Tunnel and the NASA 607 icing research aircraft, a "flying laboratory" for studying icing conditions. The branch has also developed a suite of software programs to help aircraft and icing protection system designers understand the behavior of ice accumulation on various surfaces and in various conditions. One of these innovations is the LEWICE ice accretion simulation software. Initially developed in the 1980s (when Glenn was known as Lewis Research Center), LEWICE has become one of the most widely used tools in icing research and aircraft design and certification. LEWICE has been transformed over

The evolution of a coupled ice shelf-ocean system under different climate states

NASA Astrophysics Data System (ADS)

Grosfeld, Klaus; Sandhäger, Henner

2004-07-01

Based on a new approach for coupled applications of an ice shelf model and an ocean general circulation model, we investigate the evolution of an ice shelf-ocean system and its sensitivity to changed climatic boundary conditions. Combining established 3D models into a coupled model system enabled us to study the reaction and feedbacks of each component to changes at their interface, the ice shelf base. After calculating the dynamics for prescribed initial ice shelf and bathymetric geometries, the basal mass balance determines the system evolution. In order to explore possible developments for given boundary conditions, an idealized geometry has been chosen, reflecting basic features of the Filchner-Ronne Ice Shelf, Antarctica. The model system is found to be especially sensitive in regions where high ablation or accretion rates occur. Ice Shelf Water formation as well as the build up of a marine ice body, resulting from accretion of marine ice, is simulated, indicating strong interaction processes. To improve consistency between modeled and observed ice shelf behavior, we incorporate the typical cycle of steady ice front advance and sudden retreat due to tabular iceberg calving in our time-dependent simulations. Our basic hypothesis is that iceberg break off is associated with abrupt crack propagation along elongated anomalies of the inherent stress field of the ice body. This new concept yields glaciologically plausible results and represents an auspicious basis for the development of a thorough calving criterion. Experiments under different climatic conditions (ocean warming of 0.2 and 0.5 °C and doubled surface accumulation rates) show the coupled model system to be sensitive especially to ocean warming. Increased basal melt rates of 100% for the 0.5 °C ocean warming scenario and an asymmetric development of ice shelf thicknesses suggest a high vulnerability of ice shelf regions, which represent pivotal areas between the Antarctic Ice Sheet and the Southern

CJ2 Icing Effects Simulator. Delivery Order 0019: Development of an Icing Effects Simulation for a Typical Business Jet Configuration

DTIC Science & Technology

2007-08-01

considered were: - Icing protection system failure ice - Inter-cycle (roughness) ice - Run-back ice. The study entailed wind tunnel tests of different...jet that incorporates the effects of various forms of ice. The ice conditions considered were:  Icing protection system failure ice  Inter-cycle...accretions. These were pre-activation roughness, runback shapes that form downstream of the thermal wing ice protection system , and a wing ice

Frazil-ice growth rate and dynamics in mixed layers and sub-ice-shelf plumes

NASA Astrophysics Data System (ADS)

Rees Jones, David W.; Wells, Andrew J.

2018-01-01

The growth of frazil or granular ice is an important mode of ice formation in the cryosphere. Recent advances have improved our understanding of the microphysical processes that control the rate of ice-crystal growth when water is cooled beneath its freezing temperature. These advances suggest that crystals grow much faster than previously thought. In this paper, we consider models of a population of ice crystals with different sizes to provide insight into the treatment of frazil ice in large-scale models. We consider the role of crystal growth alongside the other physical processes that determine the dynamics of frazil ice. We apply our model to a simple mixed layer (such as at the surface of the ocean) and to a buoyant plume under a floating ice shelf. We provide numerical calculations and scaling arguments to predict the occurrence of frazil-ice explosions, which we show are controlled by crystal growth, nucleation, and gravitational removal. Faster crystal growth, higher secondary nucleation, and slower gravitational removal make frazil-ice explosions more likely. We identify steady-state crystal size distributions, which are largely insensitive to crystal growth rate but are affected by the relative importance of secondary nucleation to gravitational removal. Finally, we show that the fate of plumes underneath ice shelves is dramatically affected by frazil-ice dynamics. Differences in the parameterization of crystal growth and nucleation give rise to radically different predictions of basal accretion and plume dynamics, and can even impact whether a plume reaches the end of the ice shelf or intrudes at depth.

A Wind Tunnel Study of Icing Effects on a Business Jet Airfoil

NASA Technical Reports Server (NTRS)

Addy, Harold E., Jr.; Broeren, Andy P.; Zoeckler, Joesph G.; Lee, Sam

2003-01-01

Aerodynamic wind tunnel tests were conducted to study the effects of various ice accretions on the aerodynamic performance of a 36-inch chord, two-dimensional business jet airfoil. Eight different ice shape configurations were tested. Four were castings made from molds of ice shapes accreted in an icing wind tunnel. Two were made using computationally smoothed tracings of two of the ice shapes accreted in the icing tunnel. These smoothed profiles were then extended in the spanwise direction to form a two-dimensional ice shape. The final two configurations were formed by applying grit to the smoothed ice shapes. The ice shapes resulted in as much as 48% reduction in maximum lift coefficient from that of the clean airfoil. Large increases in drag and changes in pitching moment were also observed. The castings and their corresponding smoothed counterparts yielded similar results. Little change in performance was observed with the addition of grit to the smoothed ice shapes. Changes in the Reynolds number (from 3 x 10(exp 6) to 10.5 x 10(exp 6) and Mach number (from 0.12 to 0.28) did not significantly affect the iced-airfoil performance coefficients.

Comparisons of Mixed-Phase Icing Cloud Simulations with Experiments Conducted at the NASA Propulsion Systems Laboratory

NASA Technical Reports Server (NTRS)

Bartkus, Tadas; Tsao, Jen-Ching; Struk, Peter

2017-01-01

This paper builds on previous work that compares numerical simulations of mixed-phase icing clouds with experimental data. The model couples the thermal interaction between ice particles and water droplets of the icing cloud with the flowing air of an icing wind tunnel for simulation of NASA Glenn Research Centers (GRC) Propulsion Systems Laboratory (PSL). Measurements were taken during the Fundamentals of Ice Crystal Icing Physics Tests at the PSL tunnel in March 2016. The tests simulated ice-crystal and mixed-phase icing that relate to ice accretions within turbofan engines.

Ice shelf fracture parameterization in an ice sheet model

NASA Astrophysics Data System (ADS)

Sun, Sainan; Cornford, Stephen L.; Moore, John C.; Gladstone, Rupert; Zhao, Liyun

2017-11-01

Floating ice shelves exert a stabilizing force onto the inland ice sheet. However, this buttressing effect is diminished by the fracture process, which on large scales effectively softens the ice, accelerating its flow, increasing calving, and potentially leading to ice shelf breakup. We add a continuum damage model (CDM) to the BISICLES ice sheet model, which is intended to model the localized opening of crevasses under stress, the transport of those crevasses through the ice sheet, and the coupling between crevasse depth and the ice flow field and to carry out idealized numerical experiments examining the broad impact on large-scale ice sheet and shelf dynamics. In each case we see a complex pattern of damage evolve over time, with an eventual loss of buttressing approximately equivalent to halving the thickness of the ice shelf. We find that it is possible to achieve a similar ice flow pattern using a simple rule of thumb: introducing an enhancement factor ˜ 10 everywhere in the model domain. However, spatially varying damage (or equivalently, enhancement factor) fields set at the start of prognostic calculations to match velocity observations, as is widely done in ice sheet simulations, ought to evolve in time, or grounding line retreat can be slowed by an order of magnitude.

Surface and basal ice shelf mass balance processes of the Southern McMurdo Ice Shelf determined through radar statistical reconnaissance

NASA Astrophysics Data System (ADS)

Grima, C.; Koch, I.; Greenbaum, J. S.; Soderlund, K. M.; Blankenship, D. D.; Young, D. A.; Fitzsimons, S.

2017-12-01

The McMurdo ice shelves (northern and southern MIS), adjacent to the eponymous station and the Ross Ice Shelf, Antarctica, are known for large gradients in surface snow accumulation and snow/ice impurities. Marine ice accretion and melting are important contributors to MIS's mass balance. Due to erosive winds, the southern MIS (SMIS) shows a locally negative surface mass balance. Thus, marine ice once accreted at the ice shelf base crops out at the surface. However, the exact processes that exert primary control on SMIS mass balance have remained elusive. Radar statistical reconnaissance (RSR) is a recent technique that has been used to characterize the surface properties of the Earth's cryosphere, Mars, and Titan from the stochastic character of energy scattered by the surface. Here, we apply RSR to map the surface density and roughness of the SMIS and extend the technique to derive the basal reflectance and scattering coefficients of the ice-ocean interface. We use an airborne radar survey grid acquired over the SMIS in the 2014-2015 austral summer by the University of Texas Institute for Geophysics with the High Capability Radar Sounder (HiCARS2; 60-MHz center frequency and 15-MHz bandwidth). The RSR-derived snow density values and patterns agree with directly -measured ice shelf surface accumulation rates. We also compare the composition of SMIS ice surface samples to test the ability of RSR to discriminate ices with varying dielectric properties (e.g., marine versus meteoric ice) and hypothesize relationships between the RSR-derived basal reflectance/scattered coefficients and accretion or melting at the ice-ocean interface. This improved knowledge of air-ice and ice-ocean boundaries provides a new perspective on the processes governing SMIS surface and basal mass balance.

A Protosolar Nebula Origin for the Ices Agglomerated by Comet 67P/Churyumov-Gerasimenko

NASA Astrophysics Data System (ADS)

Mousis, O.; Lunine, J. I.; Luspay-Kuti, A.; Guillot, T.; Marty, B.; Ali-Dib, M.; Wurz, P.; Altwegg, K.; Bieler, A.; Hässig, M.; Rubin, M.; Vernazza, P.; Waite, J. H.

2016-03-01

The nature of the icy material accreted by comets during their formation in the outer regions of the protosolar nebula (PSN) is a major open question in planetary science. Some scenarios of comet formation predict that these bodies agglomerated from crystalline ices condensed in the PSN. Concurrently, alternative scenarios suggest that comets accreted amorphous ice originating from the interstellar cloud or from the very distant regions of the PSN. On the basis of existing laboratory and modeling data, we find that the N2/CO and Ar/CO ratios measured in the coma of the Jupiter-family comet 67P/Churyumov-Gerasimenko by the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis instrument on board the European Space Agency’s Rosetta spacecraft match those predicted for gases trapped in clathrates. If these measurements are representative of the bulk N2/CO and Ar/CO ratios in 67P/Churyumov-Gerasimenko, it implies that the ices accreted by the comet formed in the nebula and do not originate from the interstellar medium, supporting the idea that the building blocks of outer solar system bodies have been formed from clathrates and possibly from pure crystalline ices. Moreover, because 67P/Churyumov-Gerasimenko is impoverished in Ar and N2, the volatile enrichments observed in Jupiter’s atmosphere cannot be explained solely via the accretion of building blocks with similar compositions and require an additional delivery source. A potential source may be the accretion of gas from the nebula that has been progressively enriched in heavy elements due to photoevaporation.

Experimental investigation of passive infrared ice detection for helicopter applications

NASA Technical Reports Server (NTRS)

Dershowitz, Adam; Hansman, R. John, Jr.

1991-01-01

A technique is proposed to remotely detect rotor icing on helicopters. Using passive infrared (IR) thermometry it is possible to detect the warming caused by latent heat released as supercooled water freezes. During icing, the ice accretion region on the blade leading edge will be warmer than the uniced trailing edge resulting in a chordwise temperature profile characteristic of icing. Preliminary tests were conducted on a static model in the NASA Icing Research Tunnel for a variety of wet (glaze) and dry (rime) ice conditions. The characteristic chordwise temperature profiles were observed with an IR thermal video system and confirmed with thermocouple measurements. A prototype detector system was built consisting of a single point IR pyrometer, and experiments were run on a small scale rotor model. Again the characteristic chordwise temperature profiles were observed during icing, and the IR system was able to remotely detect icing. Based on the static and subscale rotor tests the passive IR technique is promising for rotor ice detection.

Experimental investigation of passive infrared ice detection for helicopter applications

NASA Technical Reports Server (NTRS)

Dershowitz, Adam; Hansman, R. John, Jr.

1991-01-01

A technique is proposed to remotely detect rotor icing on helicopters. Using passive infrared (IR) thermometry, it is possible to detect the warming caused by latent heat released as supercooled water freezes. During icing, the ice accretion region on the blade leading edge will be warmer than the uniced trailing edge, resulting in a chordwise temperature profile characteristic of icing. Preliminary tests were conducted on a static model in the NASA Icing Research Tunnel for a variety of wet (glaze) and dry (rime) ice conditions. The characteristic chordwise temperature profiles were observed with an IR thermal video system and confirmed with thermocouple measurements. A prototype detector system was built consisting of a single point IR pyrometer. Experiments were run on a small scale rotor model. Again, the characteristic chordwise temperature profiles were observed during icing, and the IR system was able to remotely detect icing. Based on the static and subscale rotor tests, the passive IR technique is promising for rotor ice detection.

Accretion Models for Young Neutron Stars

NASA Astrophysics Data System (ADS)

Alpar, M. A.

2003-07-01

Interaction with possible fallback material, along with the magnetic fields and rotation rates at birth should determine the fates and categories of young neutron stars. This paper addresses some issues related to pure or hybrid accretion models for explaining the properties of young neutron stars.

Towards combined modeling of planetary accretion and differentiation

NASA Astrophysics Data System (ADS)

Golabek, G. J.; Gerya, T. V.; Morishima, R.; Tackley, P. J.; Labrosse, S.

2012-09-01

accretion yield an onion-like thermal structure with very high internal temperatures due to powerful short-lived radiogenic heating in the planetesimals. These lead to extensive silicate melting in the parent bodies. Yet, magma ocean and impact processes are not considered in these models and core formation is, if taken into account, assumed to be instantaneous with no feedback on the mantle evolution. It was pointed out that impacts can not only deposit heat deep into the target body, which is later buried by ejecta of further impacts [1], but also that impacts expose in the crater region originally deep-seated layers, thus cooling the interior [2]. This combination of impact effects becomes even more important when we consider that planetesimals of all masses contribute to planetary accretion. This leads occasionally to collisions between bodies with large ratios between impactor and target mass. Thus, all these processes can be expected to have a profound effect on the thermal evolution during the epoch of planetary accretion and may have implications for the onset of mantle convection and cannot be described properly in 1D geometry. Here we present a new methodology, which can be used to simulate the internal evolution of a planetary body during accretion and differentiation: Using the N-body code PKDGRAV[3] we simulate the accretion of planetary embryos from an initial annulus of several thousand planetesimals. The growth history of the largest resulting planetary embryo is used as an input for the thermomechanical 2D code I2ELVIS [4]. The thermomechanical model takes recent parametrizations of impact processes like impact heating and crater excavation [5] into account. The model also includes both long- and short-lived radiogenic isotopes and a more realistic treatment of largely molten silicates [6]. Results show that late-formed planetesimals do not experience silicate melting and avoid thermal alteration, whereas in early-formed bodies accretion and iron

Determination of Shed Ice Particle Size Using High Speed Digital Imaging

NASA Technical Reports Server (NTRS)

Broughton, Howard; Owens, Jay; Sims, James J.; Bond, Thomas H.

1996-01-01

A full scale model of an aircraft engine inlet was tested at NASA Lewis Research Center's Icing Research Tunnel. Simulated natural ice sheds from the engine inlet lip were studied using high speed digital image acquisition and image analysis. Strategic camera placement integrated at the model design phase allowed the study of ice accretion on the inlet lip and the resulting shed ice particles at the aerodynamic interface plane at the rear of the inlet prior to engine ingestion. The resulting digital images were analyzed using commercial and proprietary software to determine the size of the ice particles that could potentially be ingested by the engine during a natural shedding event. A methodology was developed to calibrate the imaging system and insure consistent and accurate measurements of the ice particles for a wide range of icing conditions.

A coupled ice-ocean model of upwelling in the marginal ice zone

NASA Technical Reports Server (NTRS)

Roed, L. P.; Obrien, J. J.

1983-01-01

A dynamical coupled ice-ocean numerical model for the marginal ice zone (MIZ) is suggested and used to study upwelling dynamics in the MIZ. The nonlinear sea ice model has a variable ice concentration and includes internal ice stress. The model is forced by stresses on the air/ocean and air/ice surfaces. The main coupling between the ice and the ocean is in the form of an interfacial stress on the ice/ocean interface. The ocean model is a linear reduced gravity model. The wind stress exerted by the atmosphere on the ocean is proportional to the fraction of open water, while the interfacial stress ice/ocean is proportional to the concentration of ice. A new mechanism for ice edge upwelling is suggested based on a geostrophic equilibrium solution for the sea ice medium. The upwelling reported in previous models invoking a stationary ice cover is shown to be replaced by a weak downwelling due to the ice motion. Most of the upwelling dynamics can be understood by analysis of the divergence of the across ice edge upper ocean transport. On the basis of numerical model, an analytical model is suggested that reproduces most of the upwelling dynamics of the more complex numerical model.

A coupled ice-ocean model of ice breakup and banding in the marginal ice zone

NASA Technical Reports Server (NTRS)

Smedstad, O. M.; Roed, L. P.

1985-01-01

A coupled ice-ocean numerical model for the marginal ice zone is considered. The model consists of a nonlinear sea ice model and a two-layer (reduced gravity) ocean model. The dependence of the upwelling response on wind stress direction is discussed. The results confirm earlier analytical work. It is shown that there exist directions for which there is no upwelling, while other directions give maximum upwelling in terms of the volume of uplifted water. The ice and ocean is coupled directly through the stress at the ice-ocean interface. An interesting consequence of the coupling is found in cases when the ice edge is almost stationary. In these cases the ice tends to break up a few tenths of kilometers inside of the ice edge.

A PROTOSOLAR NEBULA ORIGIN FOR THE ICES AGGLOMERATED BY COMET 67P/CHURYUMOV–GERASIMENKO

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mousis, O.; Vernazza, P.; Lunine, J. I.

The nature of the icy material accreted by comets during their formation in the outer regions of the protosolar nebula (PSN) is a major open question in planetary science. Some scenarios of comet formation predict that these bodies agglomerated from crystalline ices condensed in the PSN. Concurrently, alternative scenarios suggest that comets accreted amorphous ice originating from the interstellar cloud or from the very distant regions of the PSN. On the basis of existing laboratory and modeling data, we find that the N{sub 2}/CO and Ar/CO ratios measured in the coma of the Jupiter-family comet 67P/Churyumov–Gerasimenko by the Rosetta Orbitermore » Spectrometer for Ion and Neutral Analysis instrument on board the European Space Agency’s Rosetta spacecraft match those predicted for gases trapped in clathrates. If these measurements are representative of the bulk N{sub 2}/CO and Ar/CO ratios in 67P/Churyumov–Gerasimenko, it implies that the ices accreted by the comet formed in the nebula and do not originate from the interstellar medium, supporting the idea that the building blocks of outer solar system bodies have been formed from clathrates and possibly from pure crystalline ices. Moreover, because 67P/Churyumov–Gerasimenko is impoverished in Ar and N{sub 2}, the volatile enrichments observed in Jupiter’s atmosphere cannot be explained solely via the accretion of building blocks with similar compositions and require an additional delivery source. A potential source may be the accretion of gas from the nebula that has been progressively enriched in heavy elements due to photoevaporation.« less

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

NASA Technical Reports Server (NTRS)

Jun, GaRam; Oliden, Daniel; Potapczuk, Mark G.; Tsao, Jen-Ching

2014-01-01

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.

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

NASA Technical Reports Server (NTRS)

Jun, Garam; Oliden, Daniel; Potapczuk, Mark G.; Tsao, Jen-Ching

2014-01-01

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.

Finite wing aerodynamics with simulated glaze ice

NASA Technical Reports Server (NTRS)

Khodadoust, A.; Bragg, M. B.; Kerho, M.; Wells, S.; Soltani, M. R.

1992-01-01

The effect of a simulated glaze ice accretion on the aerodynamic performance of a three-dimensional wing is studied experimentally. The model used for these tests was a semi-span wing of effective aspect ratio five, mounted from the sidewall of the UIUC subsonic wind tunnel. The model has an NACA 0012 airfoil section on a rectangular, untwisted planform with interchangeable leading edges to allow for testing both the baseline and the iced wing geometry. A three-component sidewall balance was used to measure lift, drag and pitching moment on the clean and iced model. A four-beam two-color fiberoptic laser Doppler velocimeter (LDV) was used to map the flowfield along several spanwise cuts on the model. Preliminary results from LDV scans, which will be the bulk of this paper, are presented following the force balance measurement results. Initial comparison of LDV surveys compare favorably with inviscid theory results and 2D split hot-film measurements near the model surface.

Ice Roughness and Thickness Evolution on a Swept NACA 0012 Airfoil

NASA Technical Reports Server (NTRS)

McClain, Stephen T.; Vargas, Mario; Tsao, Jen-Ching

2017-01-01

Several recent studies have been performed in the Icing Research Tunnel (IRT) at NASA Glenn Research Center focusing on the evolution, spatial variations, and proper scaling of ice roughness on airfoils without sweep exposed to icing conditions employed in classical roughness studies. For this study, experiments were performed in the IRT to investigate the ice roughness and thickness evolution on a 91.44-cm (36-in.) chord NACA 0012 airfoil, swept at 30-deg with 0deg angle of attack, and exposed to both Appendix C and Appendix O (SLD) icing conditions. The ice accretion event times used in the study were less than the time required to form substantially three-dimensional structures, such as scallops, on the airfoil surface. Following each ice accretion event, the iced airfoils were scanned using a ROMER Absolute Arm laser-scanning system. The resulting point clouds were then analyzed using the self-organizing map approach of McClain and Kreeger to determine the spatial roughness variations along the surfaces of the iced airfoils. The resulting measurements demonstrate linearly increasing roughness and thickness parameters with ice accretion time. Further, when compared to dimensionless or scaled results from unswept airfoil investigations, the results of this investigation indicate that the mechanisms for early stage roughness and thickness formation on swept wings are similar to those for unswept wings.

Ice Shapes on a Tail Rotor

NASA Technical Reports Server (NTRS)

Kreeger, Richard E.; Tsao, Jen-Ching

2014-01-01

Testing of a thermally-protected helicopter rotor in the Icing Research Tunnel (IRT) was completed. Data included inter-cycle and cold blade ice shapes. Accreted ice shapes were thoroughly documented, including tracing, scanning and photographing. This was the first time this scanning capability was used outside of NASA. This type of data has never been obtained for a rotorcraft before. This data will now be used to validate the latest generation of icing analysis tools.

Observations on the Growth of Roughness Elements Into Icing Feathers

NASA Technical Reports Server (NTRS)

Vargas, Mario; Tsao, Jen, Ching

2007-01-01

This work presents the results of an experiment conducted in the Icing Research Tunnel at NASA Glenn Research Center to understand the process by which icing feathers are formed in the initial stages of ice accretion formation on swept wings. Close-up photographic data were taken on an aluminum NACA 0012 swept wing tip airfoil. Two types of photographic data were obtained: time sequence close-up photographic data during the run and close-up photographic data of the ice accretion at the end of each run. Icing runs were conducted for short ice accretion times from 10 to 180 sec. The time sequence close-up photographic data was used to study the process frame by frame and to create movies of how the process developed. The movies confirmed that at glaze icing conditions in the attachment line area icing feathers develop from roughness elements. The close-up photographic data at the end of each run showed that roughness elements change into a pointed shape with an upstream facet and join on the side with other elements having the same change to form ridges with pointed shape and upstream facet. The ridges develop into feathers when the upstream facet grows away to form the stem of the feather. The ridges and their growth into feathers were observed to form the initial scallop tips present in complete scallops.

Development of Experimental Icing Simulation Capability for Full-Scale Swept Wings: Hybrid Design Process, Years 1 and 2

NASA Technical Reports Server (NTRS)

Fujiwara, Gustavo; Bragg, Mike; Triphahn, Chris; Wiberg, Brock; Woodard, Brian; Loth, Eric; Malone, Adam; Paul, Bernard; Pitera, David; Wilcox, Pete;

A passive infrared ice detection technique for helicopter applications

NASA Technical Reports Server (NTRS)

Dershowitz, Adam L.; Hansman, R. John, Jr.

1991-01-01

A technique has been developed, and successfully tested, to detect icing remotely on helicopter rotor blades. Using passive infrared (IR) thermometry it is possible to detect the warming caused by latent heat released as supercooled water freezes. During icing, the ice accretion region on the leading edge of the blade is found to be warmer than the uniced trailing edge resulting in a chordwise temperature profile characteristic of icing. Preliminary tests, using an IR Thermal video system, were conducted on a static model in the NASA Icing Research Tunnel (IRT) for a variety of wet (glaze) and dry (rime) ice conditions. A prototype detector system was built consisting of a single point IR pyrometer, and experiments were run on a small scale rotor model. Using this prototype detector, the characteristic chordwise temperature profiles were again observed for a range of icing conditions. Several signal processing methods were investigated, to allow automatic recognition of the icing signature. Additionally, several implementation issues were considered. Based on both the static and subscale rotor tests, where ice was successfully detected, the passive IR technique appears to be promising for rotor ice detection.

An Integrated Approach to Swept Wing Icing Simulation

NASA Technical Reports Server (NTRS)

Potapczuk, Mark G.; Broeren, Andy P.

2017-01-01

This paper describes the various elements of a simulation approach used to develop a database of ice shape geometries and the resulting aerodynamic performance data for a representative commercial transport wing model exposed to a variety of icing conditions. This effort included testing in the NASA Icing Research Tunnel, the Wichita State University Walter H. Beech Wind Tunnel, and the ONERA F1 Subsonic Wind Tunnel as well as the use of ice accretion codes, an inviscid design code, and computational fluid dynamics codes. Additionally, methods for capturing full three-dimensional ice shape geometries, geometry interpolation along the span of the wing, and creation of artificial ice shapes based upon that geometric data were developed for this effort. The icing conditions used for this effort were representative of actual ice shape encounter scenarios and run the gamut from ice roughness to full three-dimensional scalloped ice shapes. The effort is still underway so this paper is a status report of work accomplished to date and a description of the remaining elements of the effort.

A model for accretion of the terrestrial planets

NASA Technical Reports Server (NTRS)

Weidenschilling, S. J.

1974-01-01

One possible origin of the terrestrial planets involves their formation by gravitational accretion of particles originally in Keplerian orbits about the sun. Some implications of this theory are considered. A formal expression for the rate of mass accretion by a planet is developed. The formal singularity of the gravitational collision cross section for low relative velocities is shown to be without physical significance when the accreting bodies are in heliocentric orbits. The distribution of particle velocities relative to an accreting planet is considered; the mean velocity increases with time. The internal temperature of an accreting planet is shown to depend simply on the accretion rate. A simple and physically reasonable approximate expression for a planetary accretion rate is proposed.

Capabilities and performance of Elmer/Ice, a new generation ice-sheet model

NASA Astrophysics Data System (ADS)

Gagliardini, O.; Zwinger, T.; Gillet-Chaulet, F.; Durand, G.; Favier, L.; de Fleurian, B.; Greve, R.; Malinen, M.; Martín, C.; Råback, P.; Ruokolainen, J.; Sacchettini, M.; Schäfer, M.; Seddik, H.; Thies, J.

2013-03-01

The Fourth IPCC Assessment Report concluded that ice-sheet flow models are unable to forecast the current increase of polar ice sheet discharge and the associated contribution to sea-level rise. Since then, the glaciological community has undertaken a huge effort to develop and improve a new generation of ice-flow models, and as a result, a significant number of new ice-sheet models have emerged. Among them is the parallel finite-element model Elmer/Ice, based on the open-source multi-physics code Elmer. It was one of the first full-Stokes models used to make projections for the evolution of the whole Greenland ice sheet for the coming two centuries. Originally developed to solve local ice flow problems of high mechanical and physical complexity, Elmer/Ice has today reached the maturity to solve larger scale problems, earning the status of an ice-sheet model. Here, we summarise almost 10 yr of development performed by different groups. We present the components already included in Elmer/Ice, its numerical performance, selected applications, as well as developments planned for the future.

An Experimental Study of Airfoil Icing Characteristics

NASA Technical Reports Server (NTRS)

Shaw, R. J.; Sotos, R. G.; Solano, F. R.

1982-01-01

A full scale general aviation wing with a NACA 63 sub 2 A415 airfoil section was tested to determine icing characteristics for representative rime and glaze icing conditions. Measurements were made of ice accretion shapes and resultant wing section drag coefficient levels. It was found that the NACA 63 sub 2 A415 wing section was less sensitive to rime and glaze icing encounters for climb conditions.

Icing Test Results on an Advanced Two-Dimensional High-Lift Multi-Element Airfoil

NASA Technical Reports Server (NTRS)

Shin, Jaiwon; Wilcox, Peter; Chin, Vincent; Sheldon, David

1994-01-01

An experimental study has been conducted to investigate ice accretions on a high-lift, multi-element airfoil in the Icing Research Tunnel at the NASA Lewis Research Center. The airfoil is representative of an advanced transport wing design. The experimental work was conducted as part of a cooperative program between McDonnell Douglas Aerospace and the NASA Lewis Research Center to improve current understanding of ice accretion characteristics on the multi-element airfoil. The experimental effort also provided ice shapes for future aerodynamic tests at flight Reynolds numbers to ascertain high-lift performance effects. Ice shapes documented for a landing configuration over a variety of icing conditions are presented along with analyses.

Modelling sea ice dynamics

NASA Astrophysics Data System (ADS)

Murawski, Jens; Kleine, Eckhard

2017-04-01

Sea ice remains one of the frontiers of ocean modelling and is of vital importance for the correct forecasts of the northern oceans. At large scale, it is commonly considered a continuous medium whose dynamics is modelled in terms of continuum mechanics. Its specifics are a matter of constitutive behaviour which may be characterised as rigid-plastic. The new developed sea ice dynamic module bases on general principles and follows a systematic approach to the problem. Both drift field and stress field are modelled by a variational property. Rigidity is treated by Lagrangian relaxation. Thus one is led to a sensible numerical method. Modelling fast ice remains to be a challenge. It is understood that ridging and the formation of grounded ice keels plays a role in the process. The ice dynamic model includes a parameterisation of the stress associated with grounded ice keels. Shear against the grounded bottom contact might lead to plastic deformation and the loss of integrity. The numerical scheme involves a potentially large system of linear equations which is solved by pre-conditioned iteration. The entire algorithm consists of several components which result from decomposing the problem. The algorithm has been implemented and tested in practice.

The Role of the Outer Boundary Condition in Accretion Disk Models: Theory and Application

NASA Astrophysics Data System (ADS)

Yuan, Feng; Peng, Qiuhe; Lu, Ju-fu; Wang, Jianmin

2000-07-01

In a previous paper, we find that the outer boundary conditions (OBCs) of an optically thin accretion flow play an important role in determining the structure of the flow. Here in this paper, we further investigate the influence of OBCs on the dynamics and radiation of the accretion flow on a more detailed level. Bremsstrahlung and synchrotron radiations amplified by Comptonization are taken into account, and two-temperature plasma assumption is adopted. The three OBCs we adopted are the temperatures of the electrons and ions and the specific angular momentum of the accretion flow at a certain outer boundary. We investigate the individual role of each of the three OBCs on the dynamical structure and the emergent spectrum. We find that when the general parameters such as the mass accretion rate M and the viscous parameter α are fixed the peak flux at various bands such as radio, IR, and X-ray can differ by as much as several orders of magnitude under different OBCs in our example. Our results indicate that the OBC is both dynamically and radiatively important and therefore should be regarded as a new ``parameter'' in accretion disk models. As an illustrative example, we further apply the above results to the compact radio source Sgr A* located at the center of our Galaxy. The advection-dominated accretion flow (ADAF) model has turned out to be a great success in explaining its luminosity and spectrum. However, there exists a discrepancy between the mass accretion rate favored by ADAF models in the literature and that favored by the three-dimensional hydrodynamical simulation, with the former being 10-20 times smaller than the latter. By seriously considering the outer boundary condition of the accretion flow, we find that because of the low specific angular momentum of the accretion gas the accretion in Sgr A* should belong to a new accretion pattern, which is characterized by the possession of a very large sonic radius. This accretion pattern can significantly

Basal Freeze-on: An Active Component of Hydrology from the Ice Divide to the Margin

NASA Astrophysics Data System (ADS)

Bell, R. E.; Tinto, K. J.; Abdi, A.; Creyts, T. T.; Wolovick, M.; Das, I.; Ferraccioli, F.; Csatho, B. M.

2012-12-01

Subglacial hydrology is considered a key control of ice sheet dynamics. Here we show that basal freeze-on is a process that can terminate basal hydrologic networks both in the interior of East Antarctica and at the margins of the Greenland Ice Sheet. Basal freeze-on modifies the ice thickness, ice structure, and ice rheology and therefore must be considered in developing accurate understanding of how hydrology interacts with ice dynamics. In East Antarctica, the freeze-on process follows well-defined hydrologic networks within Gamburtsev Mountain valleys. The steep mountain topography strongly controls the routing of the subglacial water. Ice surface slope drives the water up the mountain valleys and freeze-on occurs at the valley heads. Freeze-on ice is characterized by distinct basal radar reflectors that emerge from the hydrologic network. Evidence that these spatially coherent reflectors demark accreted ice is the upward deflection of the overlying internal layers accompanied by thickening of base of the ice sheet. Individual accretion bodies can be 25 km wide across flow, 100 km along flow with average thicknesses of ~500m although the maximum thickness is 1100m. Regional accumulation rates near the accretion sites average 4cm/yr with low ice velocity (1.5 m/yr). The volume of the ice enclosed by the accretion ice reflectors is 45-1064 km3. The accretion occurs beneath 2200-3000m thick ice and has been persistent for at least 50,000yr. Other basal reflectors in northern Greenland appear in radar from NASA's Icebridge mission and CRESIS. To identify freeze-on ice, we use specific criteria: reflectors must originate from the bed, must be spatially continuous from line to line and the meteoric stratigraphy is deflected upward. The absence of coincident gravity anomalies indicates these reflectors define distinct packages of ice rather than frozen sediment or off-nadir subglacial topography. In the Petermann Glacier Catchment, one of the largest in northern

Capabilities and performance of the new generation ice-sheet model Elmer/Ice

NASA Astrophysics Data System (ADS)

Gagliardini, O.; Zwinger, T.; Durand, G.; Favier, L.; de Fleurian, B.; Gillet-chaulet, F.; Seddik, H.; Greve, R.; Mallinen, M.; Martin, C.; Raback, P.; Ruokolainen, J.; Schäfer, M.; Thies, J.

2012-12-01

Since the Fourth IPCC Assessment Report, and its conclusion about the inability of ice-sheet flow models to forecast the current increase of polar ice sheet discharge and associated contribution to sea-level rise, a huge development effort has been undertaken by the glaciological community. All around the world, models have been improved and, interestingly, a significant number of new ice-sheet models have emerged. Among them, the parallel finite-element model Elmer/Ice (based on the open-source multi-physics code Elmer) was one of the first full-Stokes models used to make projections of the future of the whole Greenland ice sheet for the coming two centuries. Originally developed to solve dedicated local ice flow problems of high mechanical and physical complexity, Elmer/Ice has today reached the maturity to solve larger scale problems, earning the status of an ice-sheet model. In this presentation, we summarise the almost 10 years of development performed by different groups. We present the components already included in Elmer/Ice, its numerical performance, selected applications, as well as developments planed for the future.

SPH non-Newtonian Model for Ice Sheet and Ice Shelf Dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tartakovsky, Alexandre M.; Pan, Wenxiao; Monaghan, Joseph J.

2012-07-07

We propose a new three-dimensional smoothed particle hydrodynamics (SPH) non-Newtonian model to study coupled ice sheet and ice shelf dynamics. Most existing ice sheet numerical models use a grid-based Eulerian approach, and are usually restricted to shallow ice sheet and ice shelf approximations of the momentum conservation equation. SPH, a fully Lagrangian particle method, solves the full momentum conservation equation. SPH method also allows modeling of free-surface flows, large material deformation, and material fragmentation without employing complex front-tracking schemes, and does not require re-meshing. As a result, SPH codes are highly scalable. Numerical accuracy of the proposed SPH model ismore » first verified by simulating a plane shear flow with a free surface and the propagation of a blob of ice along a horizontal surface. Next, the SPH model is used to investigate the grounding line dynamics of ice sheet/shelf. The steady position of the grounding line, obtained from our SPH simulations, is in good agreement with laboratory observations for a wide range of bedrock slopes, ice-to-fluid density ratios, and flux. We examine the effect of non-Newtonian behavior of ice on the grounding line dynamics. The non-Newtonian constitutive model is based on Glen's law for a creeping flow of a polycrystalline ice. Finally, we investigate the effect of a bedrock geometry on a steady-state position of the grounding line.« less

Incompressible Wind Accretion

NASA Astrophysics Data System (ADS)

Tejeda, E.

2018-04-01

We present a simple, analytic model of an incompressible fluid accreting onto a moving gravitating object. This solution allows us to probe the highly subsonic regime of wind accretion. Moreover, it corresponds to the Newtonian limit of a previously known relativistic model of a stiff fluid accreting onto a black hole. Besides filling this blank in the literature, the new solution should be useful as a benchmark test for numerical hydrodynamics codes. Given its simplicity, it can also be used as an illustrative example in a gas dynamics course.

An Experimental Study of the Flowfield on a Semispan Rectangular Wing with a Simulated Glaze Ice Accretion. Ph.D. Thesis, 1993 Final Report

NASA Technical Reports Server (NTRS)

Khodadoust, Abdollah

1994-01-01

Wind tunnel experiments were conducted in order to study the effect of a simulated glaze ice accretion on the flowfield of a semispan, reflection-plane, rectangular wing at Re = 1.5 million and M = 0.12. A laser Doppler velocimeter was used to map the flowfield on the upper surface of the model in both the clean and iced configurations at alpha = 0, 4, and 8 degrees angle of attack. At low angles of attack, the massive separation bubble aft of the leading edge ice horn was found to behave in a manner similar to laminar separation bubbles. At alpha = 0 and 4 degrees, the locations of transition and reattachment, as deduced from momentum thickness distributions, were found to be in good agreement with transition and reattachment locations in laminar separation bubbles. These values at y/b = 0.470, the centerline measurement location, matched well with data obtained on a similar but two dimensional model. The measured velocity profiles on the iced wing compared reasonably with the predicted profiles from Navier-Stokes computations. The iced-induced separation bubble was also found to have features similar to the recirculating region aft of rearward-facing steps. At alpha = 0 degrees and 4 degrees, reverse flow magnitudes and turbulence intensity levels were typical of those found in the recirculating region aft of rearward-facing steps. The calculated separation streamline aft of the ice horn at alpha = 4 degrees, y/b = 0.470 coincided with the locus of the maximum Reynolds normal stress. The maximum Reynolds normal stress peaked at two locations along the separation streamline. The location of the first peak-value coincided with the transition location, as deduced from the momentum thickness distributions. The location of the second peak was just upstream of reattachment, in good agreement with measurements of flows over similar obstacles. The intermittency factor in the vicinity of reattachment at alpha = 4 degrees, y/b = 0.470, revealed the time-dependent nature of

Ocean-Forced Ice-Shelf Thinning in a Synchronously Coupled Ice-Ocean Model

NASA Astrophysics Data System (ADS)

Jordan, James R.; Holland, Paul R.; Goldberg, Dan; Snow, Kate; Arthern, Robert; Campin, Jean-Michel; Heimbach, Patrick; Jenkins, Adrian

2018-02-01

The first fully synchronous, coupled ice shelf-ocean model with a fixed grounding line and imposed upstream ice velocity has been developed using the MITgcm (Massachusetts Institute of Technology general circulation model). Unlike previous, asynchronous, approaches to coupled modeling our approach is fully conservative of heat, salt, and mass. Synchronous coupling is achieved by continuously updating the ice-shelf thickness on the ocean time step. By simulating an idealized, warm-water ice shelf we show how raising the pycnocline leads to a reduction in both ice-shelf mass and back stress, and hence buttressing. Coupled runs show the formation of a western boundary channel in the ice-shelf base due to increased melting on the western boundary due to Coriolis enhanced flow. Eastern boundary ice thickening is also observed. This is not the case when using a simple depth-dependent parameterized melt, as the ice shelf has relatively thinner sides and a thicker central "bulge" for a given ice-shelf mass. Ice-shelf geometry arising from the parameterized melt rate tends to underestimate backstress (and therefore buttressing) for a given ice-shelf mass due to a thinner ice shelf at the boundaries when compared to coupled model simulations.

CFD Analysis of the Aerodynamics of a Business-Jet Airfoil with Leading-Edge Ice Accretion

NASA Technical Reports Server (NTRS)

Chi, X.; Zhu, B.; Shih, T. I.-P.; Addy, H. E.; Choo, Y. K.

2004-01-01

For rime ice - where the ice buildup has only rough and jagged surfaces but no protruding horns - this study shows two dimensional CFD analysis based on the one-equation Spalart-Almaras (S-A) turbulence model to predict accurately the lift, drag, and pressure coefficients up to near the stall angle. For glaze ice - where the ice buildup has two or more protruding horns near the airfoil's leading edge - CFD predictions were much less satisfactory because of the large separated region produced by the horns even at zero angle of attack. This CFD study, based on the WIND and the Fluent codes, assesses the following turbulence models by comparing predictions with available experimental data: S-A, standard k-epsilon, shear-stress transport, v(exp 2)-f, and differential Reynolds stress.

Icing flight research: Aerodynamic effects of ice and ice shape documentation with stereo photography

NASA Technical Reports Server (NTRS)

Mikkelsen, K. L.; Mcknight, R. C.; Ranaudo, R. J.; Perkins, P. J., Jr.

1985-01-01

Aircraft icing flight research was performed in natural icing conditions. A data base consisting of icing cloud measurements, ice shapes, and aerodynamic measurements is being developed. During research icing encounters the icing cloud was continuously measured. After the encounter, the ice accretion shapes on the wing were documented with a stereo camera system. The increase in wing section drag was measured with a wake survey probe. The overall aircraft performance loss in terms of lift and drag coefficient changes was obtained by steady level speed/power measurements. Selective deicing of the airframe components was performed to determine their contributions to the total drag increase. Engine out capability in terms of power available was analyzed for the iced aircraft. It was shown that the stereo photography system can be used to document ice shapes in flight and that the wake survey probe can measure increases in wing section drag caused by ice. On one flight, the wing section drag coefficient (c sub d) increased approximately 120 percent over the uniced baseline at an aircraft angle of attack of 6 deg. On another flight, the aircraft darg coefficient (c sub d) increased by 75 percent over the uniced baseline at an aircraft lift coefficient (C sub d) of 0.5.

Icing flight research - Aerodynamic effects of ice and ice shape documentation with stereo photography

NASA Technical Reports Server (NTRS)

Mikkelsen, K. L.; Mcknight, R. C.; Ranaudo, R. J.; Perkins, P. J., Jr.

1985-01-01

Aircraft icing flight research was performed in natural icing conditions. A data base consisting of icing cloud measurements, ice shapes, and aerodynamic measurements is being developed. During research icing encounters the icing cloud was continuously measured. After the encounter, the ice accretion shapes on the wing were documented with a stereo camera system. The increase in wing section drag was measured with a wake survey probe. The overall aircraft performance loss in terms of lift and drag coefficient changes were obtained by steady level speed/power measurements. Selective deicing of the airframe components was performed to determine their contributions to the total drag increase. Engine out capability in terms of power available was analyzed for the iced aircraft. It was shown that the stereo photography system can be used to document ice shapes in flight and that the wake survey probe can measure increases in wing section drag caused by ice. On one flight, the wing section drag coefficient (c sub d) increased approximately 120 percent over the uniced baseline at an aircraft angle of attack of 6 deg. On another flight, the aircraft drag coefficient (c sub d) increased by 75 percent over the uniced baseline at an aircraft lift coefficient (c sub d) of 0.5.

Modeling Wave-Ice Interactions in the Marginal Ice Zone

NASA Astrophysics Data System (ADS)

Orzech, Mark; Shi, Fengyan; Bateman, Sam; Veeramony, Jay; Calantoni, Joe

2015-04-01

The small-scale (O(m)) interactions between waves and ice floes in the marginal ice zone (MIZ) are investigated with a coupled model system. Waves are simulated with the non-hydrostatic finite-volume model NHWAVE (Ma et al., 2012) and ice floes are represented as bonded collections of smaller particles with the discrete element system LIGGGHTS (Kloss et al., 2012). The physics of fluid and ice are recreated as authentically as possible, to allow the coupled system to supplement and/or substitute for more costly and demanding field experiments. The presentation will first describe the development and validation of the coupled system, then discuss the results of a series of virtual experiments in which ice floe and wave characteristics are varied to examine their effects on energy dissipation, MIZ floe size distribution, and ice pack retreat rates. Although Wadhams et al. (1986) suggest that only a small portion (roughly 10%) of wave energy entering the MIZ is reflected, dissipation mechanisms for the remaining energy have yet to be delineated or measured. The virtual experiments are designed to focus on specific properties and processes - such as floe size and shape, collision and fracturing events, and variations in wave climate - and measure their relative roles the transfer of energy and momentum from waves to ice. Questions to be examined include: How is energy dissipated by ice floe collisions, fracturing, and drag, and how significant is the wave attenuation associated with each process? Do specific wave/floe length scale ratios cause greater wave attenuation? How does ice material strength affect the rate of wave energy loss? The coupled system will ultimately be used to test and improve upon wave-ice parameterizations for large-scale climate models. References: >Kloss, C., C. Goniva, A. Hager, S. Amberger, and S. Pirker (2012). Models, algorithms and validation for opensource DEM and CFD-DEM. Progress in Computational Fluid Dynamics 12(2/3), 140-152. >Ma, G

Validation of Modelled Ice Dynamics of the Greenland Ice Sheet using Historical Forcing

NASA Astrophysics Data System (ADS)

Hoffman, M. J.; Price, S. F.; Howat, I. M.; Bonin, J. A.; Chambers, D. P.; Tezaur, I.; Kennedy, J. H.; Lenaerts, J.; Lipscomb, W. H.; Neumann, T.; Nowicki, S.; Perego, M.; Saba, J. L.; Salinger, A.; Guerber, J. R.

2015-12-01

Although ice sheet models are used for sea level rise projections, the degree to which these models have been validated by observations is fairly limited, due in part to the limited duration of the satellite observation era and the long adjustment time scales of ice sheets. Here we describe a validation framework for the Greenland Ice Sheet applied to the Community Ice Sheet Model by forcing the model annually with flux anomalies at the major outlet glaciers (Enderlin et al., 2014, observed from Landsat/ASTER/Operation IceBridge) and surface mass balance (van Angelen et al., 2013, calculated from RACMO2) for the period 1991-2012. The ice sheet model output is compared to ice surface elevation observations from ICESat and ice sheet mass change observations from GRACE. Early results show promise for assessing the performance of different model configurations. Additionally, we explore the effect of ice sheet model resolution on validation skill.

The CMIP6 Sea-Ice Model Intercomparison Project (SIMIP): Understanding sea ice through climate-model simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Notz, Dirk; Jahn, Alexandra; Holland, Marika

A better understanding of the role of sea ice for the changing climate of our planet is the central aim of the diagnostic Coupled Model Intercomparison Project 6 (CMIP6)-endorsed Sea-Ice Model Intercomparison Project (SIMIP). To reach this aim, SIMIP requests sea-ice-related variables from climate-model simulations that allow for a better understanding and, ultimately, improvement of biases and errors in sea-ice simulations with large-scale climate models. This then allows us to better understand to what degree CMIP6 model simulations relate to reality, thus improving our confidence in answering sea-ice-related questions based on these simulations. Furthermore, the SIMIP protocol provides a standardmore » for sea-ice model output that will streamline and hence simplify the analysis of the simulated sea-ice evolution in research projects independent of CMIP. To reach its aims, SIMIP provides a structured list of model output that allows for an examination of the three main budgets that govern the evolution of sea ice, namely the heat budget, the momentum budget, and the mass budget. Furthermore, we explain the aims of SIMIP in more detail and outline how its design allows us to answer some of the most pressing questions that sea ice still poses to the international climate-research community.« less

The CMIP6 Sea-Ice Model Intercomparison Project (SIMIP): Understanding sea ice through climate-model simulations

DOE PAGES

Notz, Dirk; Jahn, Alexandra; Holland, Marika; ...

2016-09-23

A better understanding of the role of sea ice for the changing climate of our planet is the central aim of the diagnostic Coupled Model Intercomparison Project 6 (CMIP6)-endorsed Sea-Ice Model Intercomparison Project (SIMIP). To reach this aim, SIMIP requests sea-ice-related variables from climate-model simulations that allow for a better understanding and, ultimately, improvement of biases and errors in sea-ice simulations with large-scale climate models. This then allows us to better understand to what degree CMIP6 model simulations relate to reality, thus improving our confidence in answering sea-ice-related questions based on these simulations. Furthermore, the SIMIP protocol provides a standardmore » for sea-ice model output that will streamline and hence simplify the analysis of the simulated sea-ice evolution in research projects independent of CMIP. To reach its aims, SIMIP provides a structured list of model output that allows for an examination of the three main budgets that govern the evolution of sea ice, namely the heat budget, the momentum budget, and the mass budget. Furthermore, we explain the aims of SIMIP in more detail and outline how its design allows us to answer some of the most pressing questions that sea ice still poses to the international climate-research community.« less

Ice Growth Measurements from Image Data to Support Ice Crystal and Mixed-Phase Accretion Testing

NASA Technical Reports Server (NTRS)

Struk, Peter M.; Lynch, Christopher J.

2012-01-01

This paper describes the imaging techniques as well as the analysis methods used to measure the ice thickness and growth rate in support of ice-crystal icing tests performed at the National Research Council of Canada (NRC) Research Altitude Test Facility (RATFac). A detailed description of the camera setup, which involves both still and video cameras, as well as the analysis methods using the NASA Spotlight software, are presented. Two cases, one from two different test entries, showing significant ice growth are analyzed in detail describing the ice thickness and growth rate which is generally linear. Estimates of the bias uncertainty are presented for all measurements. Finally some of the challenges related to the imaging and analysis methods are discussed as well as methods used to overcome them.

Ice Growth Measurements from Image Data to Support Ice-Crystal and Mixed-Phase Accretion Testing

NASA Technical Reports Server (NTRS)

Struk, Peter, M; Lynch, Christopher, J.

2012-01-01

This paper describes the imaging techniques as well as the analysis methods used to measure the ice thickness and growth rate in support of ice-crystal icing tests performed at the National Research Council of Canada (NRC) Research Altitude Test Facility (RATFac). A detailed description of the camera setup, which involves both still and video cameras, as well as the analysis methods using the NASA Spotlight software, are presented. Two cases, one from two different test entries, showing significant ice growth are analyzed in detail describing the ice thickness and growth rate which is generally linear. Estimates of the bias uncertainty are presented for all measurements. Finally some of the challenges related to the imaging and analysis methods are discussed as well as methods used to overcome them.

CICE, The Los Alamos Sea Ice Model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hunke, Elizabeth; Lipscomb, William; Jones, Philip

The Los Alamos sea ice model (CICE) is the result of an effort to develop a computationally efficient sea ice component for a fully coupled atmosphere–land–ocean–ice global climate model. It was originally designed to be compatible with the Parallel Ocean Program (POP), an ocean circulation model developed at Los Alamos National Laboratory for use on massively parallel computers. CICE has several interacting components: a vertical thermodynamic model that computes local growth rates of snow and ice due to vertical conductive, radiative and turbulent fluxes, along with snowfall; an elastic-viscous-plastic model of ice dynamics, which predicts the velocity field of themore » ice pack based on a model of the material strength of the ice; an incremental remapping transport model that describes horizontal advection of the areal concentration, ice and snow volume and other state variables; and a ridging parameterization that transfers ice among thickness categories based on energetic balances and rates of strain. It also includes a biogeochemical model that describes evolution of the ice ecosystem. The CICE sea ice model is used for climate research as one component of complex global earth system models that include atmosphere, land, ocean and biogeochemistry components. It is also used for operational sea ice forecasting in the polar regions and in numerical weather prediction models.« less

Analysis of Antarctic Ice-Sheet Mass Balance from ICESat Measurements

NASA Technical Reports Server (NTRS)

Zwally, H. Jay; Li, Jun; Robbins, John; Saba, Jack L.; Yi, Donghui

2011-01-01

If protoplanets formed from 10 to 20 kilometer diameter planetesimals in a runaway accretion process prior to their oligarchic growth into the terrestrial planets, it is only logical to ask where these planetesimals may have formed in order to assess the initial composition of the Earth. We have used Weidenschilling's model for the formation of comets (1997) to calculate an efficiency factor for the formation of planetesimals from the solar nebula, then used this factor to calculate the feeding zones that contribute to material contained within 10, 15 and 20 kilometer diameter planetesimals at 1 A.V. as a function of nebular mass. We find that for all reasonable nebular masses, these planetesimals contain a minimum of 3% water as ice by mass. The fraction of ice increases as the planetesimals increase in size and as the nebular mass decreases, since both factors increase the feeding zones from which solids in the final planetesimals are drawn. Is there really a problem with the current accretion scenario that makes the Earth too dry, or is it possible that the nascent Earth lost significant quantities of water in the final stages of accretion?

Kinetic Monte Carlo simulations of water ice porosity: extrapolations of deposition parameters from the laboratory to interstellar space.

PubMed

Clements, Aspen R; Berk, Brandon; Cooke, Ilsa R; Garrod, Robin T

2018-02-21

Dust grains in cold, dense interstellar clouds build up appreciable ice mantles through the accretion and subsequent surface chemistry of atoms and molecules from the gas. These mantles, of thicknesses on the order of 100 monolayers, are primarily composed of H 2 O, CO, and CO 2 . Laboratory experiments using interstellar ice analogues have shown that porosity could be present and can facilitate diffusion of molecules along the inner pore surfaces. However, the movement of molecules within and upon the ice is poorly described by current chemical kinetics models, making it difficult either to reproduce the formation of experimental porous ice structures or to extrapolate generalized laboratory results to interstellar conditions. Here we use the off-lattice Monte Carlo kinetics model MIMICK to investigate the effects that various deposition parameters have on laboratory ice structures. The model treats molecules as isotropic spheres of a uniform size, using a Lennard-Jones potential. We reproduce experimental trends in the density of amorphous solid water (ASW) for varied deposition angle, rate and surface temperature; ice density decreases when the incident angle or deposition rate is increased, while increasing temperature results in a more-compact water ice. The models indicate that the density behaviour at higher temperatures (≥80 K) is dependent on molecular rearrangement resulting from thermal diffusion. To reproduce trends at lower temperatures, it is necessary to take account of non-thermal diffusion by newly-adsorbed molecules, which bring kinetic energy both from the gas phase and from their acceleration into a surface binding site. Extrapolation of the model to conditions appropriate to protoplanetary disks, in which direct accretion of water from the gas-phase may be the dominant ice formation mechanism, indicate that these ices may be less porous than laboratory ices.

Sea Ice Summer Camp: Bringing Together Arctic Sea Ice Modelers and Observers

NASA Astrophysics Data System (ADS)

Perovich, D. K.; Holland, M. M.

2016-12-01

The Arctic sea ice has undergone dramatic change and numerical models project this to continue for the foreseeable future. Understanding the mechanisms behind sea ice loss and its consequences for the larger Arctic and global systems is of critical importance if we are to anticipate and plan for the future. One impediment to progress is a disconnect between the observational and modeling communities. A sea ice summer camp was held in Barrow Alaska from 26 May to 1 June 2016 to overcome this impediment and better integrate the sea ice community. The 25 participants were a mix of modelers and observers from 13 different institutions at career stages from graduate student to senior scientist. The summer camp provided an accelerated program on sea ice observations and models and also fostered future collaborative interdisciplinary activities. Each morning was spent in the classroom with a daily lecture on an aspect of modeling or remote sensing followed by practical exercises. Topics included using models to assess sensitivity, to test hypotheses and to explore sources of uncertainty in future Arctic sea ice loss. The afternoons were spent on the ice making observations. There were four observational activities; albedo observations, ice thickness measurements, ice coring and physical properties, and ice morphology surveys. The last field day consisted of a grand challenge where the group formulated a hypothesis, developed an observational and modeling strategy to test the hypothesis, and then integrated the observations and model results. The impacts of changing sea ice are being felt today in Barrow Alaska. We opened a dialog with Barrow community members to further understand these changes. This included an evening discussion with two Barrow sea ice experts and a community presentation of our work in a public lecture at the Inupiat Heritage Center.

Capabilities and performance of Elmer/Ice, a new-generation ice sheet model

NASA Astrophysics Data System (ADS)

Gagliardini, O.; Zwinger, T.; Gillet-Chaulet, F.; Durand, G.; Favier, L.; de Fleurian, B.; Greve, R.; Malinen, M.; Martín, C.; Råback, P.; Ruokolainen, J.; Sacchettini, M.; Schäfer, M.; Seddik, H.; Thies, J.

2013-08-01

The Fourth IPCC Assessment Report concluded that ice sheet flow models, in their current state, were unable to provide accurate forecast for the increase of polar ice sheet discharge and the associated contribution to sea level rise. Since then, the glaciological community has undertaken a huge effort to develop and improve a new generation of ice flow models, and as a result a significant number of new ice sheet models have emerged. Among them is the parallel finite-element model Elmer/Ice, based on the open-source multi-physics code Elmer. It was one of the first full-Stokes models used to make projections for the evolution of the whole Greenland ice sheet for the coming two centuries. Originally developed to solve local ice flow problems of high mechanical and physical complexity, Elmer/Ice has today reached the maturity to solve larger-scale problems, earning the status of an ice sheet model. Here, we summarise almost 10 yr of development performed by different groups. Elmer/Ice solves the full-Stokes equations, for isotropic but also anisotropic ice rheology, resolves the grounding line dynamics as a contact problem, and contains various basal friction laws. Derived fields, like the age of the ice, the strain rate or stress, can also be computed. Elmer/Ice includes two recently proposed inverse methods to infer badly known parameters. Elmer is a highly parallelised code thanks to recent developments and the implementation of a block preconditioned solver for the Stokes system. In this paper, all these components are presented in detail, as well as the numerical performance of the Stokes solver and developments planned for the future.

Modeling ocean wave propagation under sea ice covers

NASA Astrophysics Data System (ADS)

Zhao, Xin; Shen, Hayley H.; Cheng, Sukun

2015-02-01

Operational ocean wave models need to work globally, yet current ocean wave models can only treat ice-covered regions crudely. The purpose of this paper is to provide a brief overview of ice effects on wave propagation and different research methodology used in studying these effects. Based on its proximity to land or sea, sea ice can be classified as: landfast ice zone, shear zone, and the marginal ice zone. All ice covers attenuate wave energy. Only long swells can penetrate deep into an ice cover. Being closest to open water, wave propagation in the marginal ice zone is the most complex to model. The physical appearance of sea ice in the marginal ice zone varies. Grease ice, pancake ice, brash ice, floe aggregates, and continuous ice sheet may be found in this zone at different times and locations. These types of ice are formed under different thermal-mechanical forcing. There are three classic models that describe wave propagation through an idealized ice cover: mass loading, thin elastic plate, and viscous layer models. From physical arguments we may conjecture that mass loading model is suitable for disjoint aggregates of ice floes much smaller than the wavelength, thin elastic plate model is suitable for a continuous ice sheet, and the viscous layer model is suitable for grease ice. For different sea ice types we may need different wave ice interaction models. A recently proposed viscoelastic model is able to synthesize all three classic models into one. Under suitable limiting conditions it converges to the three previous models. The complete theoretical framework for evaluating wave propagation through various ice covers need to be implemented in the operational ocean wave models. In this review, we introduce the sea ice types, previous wave ice interaction models, wave attenuation mechanisms, the methods to calculate wave reflection and transmission between different ice covers, and the effect of ice floe breaking on shaping the sea ice morphology

Modelling the Climate - Greenland Ice Sheet Interaction in the Coupled Ice-sheet/Climate Model EC-EARTH - PISM

NASA Astrophysics Data System (ADS)

Yang, S.; Madsen, M. S.; Rodehacke, C. B.; Svendsen, S. H.; Adalgeirsdottir, G.

2014-12-01

Recent observations show that the Greenland ice sheet (GrIS) has been losing mass with an increasing speed during the past decades. Predicting the GrIS changes and their climate consequences relies on the understanding of the interaction of the GrIS with the climate system on both global and local scales, and requires climate model systems with an explicit and physically consistent ice sheet module. A fully coupled global climate model with a dynamical ice sheet model for the GrIS has recently been developed. The model system, EC-EARTH - PISM, consists of the EC-EARTH, an atmosphere, ocean and sea ice model system, and the Parallel Ice Sheet Model (PISM). The coupling of PISM includes a modified surface physical parameterization in EC-EARTH adapted to the land ice surface over glaciated regions in Greenland. The PISM ice sheet model is forced with the surface mass balance (SMB) directly computed inside the EC-EARTH atmospheric module and accounting for the precipitation, the surface evaporation, and the melting of snow and ice over land ice. PISM returns the simulated basal melt, ice discharge and ice cover (extent and thickness) as boundary conditions to EC-EARTH. This coupled system is mass and energy conserving without being constrained by any anomaly correction or flux adjustment, and hence is suitable for investigation of ice sheet - climate feedbacks. Three multi-century experiments for warm climate scenarios under (1) the RCP85 climate forcing, (2) an abrupt 4xCO2 and (3) an idealized 1% per year CO2 increase are performed using the coupled model system. The experiments are compared with their counterparts of the standard CMIP5 simulations (without the interactive ice sheet) to evaluate the performance of the coupled system and to quantify the GrIS feedbacks. In particular, the evolution of the Greenland ice sheet under the warm climate and its impacts on the climate system are investigated. Freshwater fluxes from the Greenland ice sheet melt to the Arctic

Evolution of Interstellar Ices

NASA Astrophysics Data System (ADS)

Allamandola, Louis J.; Bernstein, Max P.; Sandford, Scott A.; Walker, Robert L.

1999-10-01

Infrared observations, combined with realistic laboratory simulations, have revolutionized our understanding of interstellar ice and dust, the building blocks of comets. Ices in molecular clouds are dominated by the very simple molecules H2O, CH3OH, NH3, CO, CO2, and probably H2CO and H2. More complex species including nitriles, ketones, and esters are also present, but at lower concentrations. The evidence for these, as well as the abundant, carbon-rich, interstellar, polycyclic aromatic hydrocarbons (PAHs) is reviewed. Other possible contributors to the interstellar/pre-cometary ice composition include accretion of gas-phase molecules and in situ photochemical processing. By virtue of their low abundance, accretion of simple gas-phase species is shown to be the least important of the processes considered in determining ice composition. On the other hand, photochemical processing does play an important role in driving dust evolution and the composition of minor species. Ultraviolet photolysis of realistic laboratory analogs readily produces H2, H2CO, CO2, CO, CH4, HCO, and the moderately complex organic molecules: CH3CH2OH (ethanol), HC(=O)NH2 (formamide), CH3C(=O)NH2 (acetamide), R-CN (nitriles), and hexamethylenetetramine (HMT, C6H12N4), as well as more complex species including amides, ketones, and polyoxymethylenes (POMs). Inclusion of PAHs in the ices produces many species similar to those found in meteorites including aromatic alcohols, quinones and ethers. Photon assisted PAH-ice deuterium exchange also occurs. All of these species are readily formed and are therefore likely cometary constituents.

An experimental study of icing control using DBD plasma actuator

NASA Astrophysics Data System (ADS)

Cai, Jinsheng; Tian, Yongqiang; Meng, Xuanshi; Han, Xuzhao; Zhang, Duo; Hu, Haiyang

2017-08-01

Ice accretion on aircraft or wind turbine has been widely recognized as a big safety threat in the past decades. This study aims to develop a new approach for icing control using an AC-DBD plasma actuator. The experiments of icing control (i.e., anti-/de-icing) on a cylinder model were conducted in an icing wind tunnel with controlled wind speed (i.e., 15 m/s) and temperature (i.e., -10°C). A digital camera was used to record the dynamic processes of plasma anti-icing and de-icing whilst an infrared imaging system was utilized to map the surface temperature variations during the anti-/de-icing processes. It was found that the AC-DBD plasma actuator is very effective in both anti-icing and de-icing operations. While no ice formation was observed when the plasma actuator served as an anti-icing device, a complete removal of the ice layer with a thickness of 5 mm was achieved by activating the plasma actuator for ˜150 s. Such information demonstrated the feasibility of plasma anti-/de-icing, which could potentially provide more effective and safer icing mitigation strategies.

Spectra of black hole accretion models of ultraluminous X-ray sources

NASA Astrophysics Data System (ADS)

Narayan, Ramesh; Sa̧dowski, Aleksander; Soria, Roberto

2017-08-01

We present general relativistic radiation magnetohydrodynamics simulations of super-Eddington accretion on a 10 M⊙ black hole. We consider a range of mass accretion rates, black hole spins and magnetic field configurations. We compute the spectra and images of the models as a function of viewing angle and compare them with the observed properties of ultraluminous X-ray sources (ULXs). The models easily produce apparent luminosities in excess of 1040 erg s-1 for pole-on observers. However, the angle-integrated radiative luminosities rarely exceed 2.5 × 1039 erg s-1 even for mass accretion rates of tens of Eddington. The systems are thus radiatively inefficient, though they are energetically efficient when the energy output in winds and jets is also counted. The simulated models reproduce the main empirical types of spectra - disc-like, supersoft, soft, hard - observed in ultraluminous X-ray sources (ULXs). The magnetic field configuration, whether 'standard and normal evolution' (SANE) or 'magnetically arrested disc' (MAD), has a strong effect on the results. In SANE models, the X-ray spectral hardness is almost independent of accretion rate, but decreases steeply with increasing inclination. MAD models with non-spinning black holes produce significantly softer spectra at higher values of \\dot{M}, even at low inclinations. MAD models with rapidly spinning black holes are unique. They are radiatively efficient (efficiency factor ˜10-20 per cent), superefficient when the mechanical energy output is also included (70 per cent) and produce hard blazar-like spectra. In all models, the emission shows strong geometrical beaming, which disagrees with the more isotropic illumination favoured by observations of ULX bubbles.

Experimental Technique and Assessment for Measuring the Convective Heat Transfer Coefficient from Natural Ice Accretions

NASA Technical Reports Server (NTRS)

Masiulaniec, K. Cyril; Vanfossen, G. James, Jr.; Dewitt, Kenneth J.; Dukhan, Nihad

1995-01-01

A technique was developed to cast frozen ice shapes that had been grown on a metal surface. This technique was applied to a series of ice shapes that were grown in the NASA Lewis Icing Research Tunnel on flat plates. Nine flat plates, 18 inches square, were obtained from which aluminum castings were made that gave good ice shape characterizations. Test strips taken from these plates were outfitted with heat flux gages, such that when placed in a dry wind tunnel, can be used to experimentally map out the convective heat transfer coefficient in the direction of flow from the roughened surfaces. The effects on the heat transfer coefficient for both parallel and accelerating flow will be studied. The smooth plate model verification baseline data as well as one ice roughened test case are presented.

Convective Enhancement of Icing Roughness Elements in Stagnation Region Flows

NASA Technical Reports Server (NTRS)

Hughes, Michael T.; McClain, Stephen T.; Vargas, Mario; Broeren, Andy

2015-01-01

To improve existing ice accretion simulation codes, more data regarding ice roughness and its effects on convective heat transfer are required. To build on existing research on this topic, this study used the Vertical Icing Studies Tunnel (VIST) at NASA Glenn Research to model realistic ice roughness in the stagnation region of a NACA 0012 airfoil. Using the VIST, a test plate representing the leading 2% chord of the airfoil was subjected to flows of 7.62 m/s (25 ft/s), 12.19 m/s (40 ft/s), and 16.76 m/s (55 ft/s). The test plate was fitted with 3 surfaces, each with a different representation of ice roughness: 1) a control surface with no ice roughness, 2) a surface with ice roughness with element height scaled by 10x and streamwise rough zone width from the stagnation point scaled by 10x, and 3) a surface with ice roughness with element height scaled by 10x and streamwise rough zone width from the stagnation point scaled by 25x. Temperature data from the tests were recorded using an infrared camera and thermocouples imbedded in the test plate. From the temperature data, a convective heat transfer coefficient map was created for each case. Additional testing was also performed to validate the VIST's flow quality. These tests included five-hole probe and hot-wire probe velocity traces to provide flow visualization and to study boundary layer formation on the various test surfaces. The knowledge gained during the experiments will help improve ice accretion codes by providing heat transfer coefficient validation data and by providing flow visualization data helping understand current and future experiments performed in the VIST.

An ice sheet model validation framework for the Greenland ice sheet

NASA Astrophysics Data System (ADS)

Price, Stephen F.; Hoffman, Matthew J.; Bonin, Jennifer A.; Howat, Ian M.; Neumann, Thomas; Saba, Jack; Tezaur, Irina; Guerber, Jeffrey; Chambers, Don P.; Evans, Katherine J.; Kennedy, Joseph H.; Lenaerts, Jan; Lipscomb, William H.; Perego, Mauro; Salinger, Andrew G.; Tuminaro, Raymond S.; van den Broeke, Michiel R.; Nowicki, Sophie M. J.

2017-01-01

We propose a new ice sheet model validation framework - the Cryospheric Model Comparison Tool (CmCt) - that takes advantage of ice sheet altimetry and gravimetry observations collected over the past several decades and is applied here to modeling of the Greenland ice sheet. We use realistic simulations performed with the Community Ice Sheet Model (CISM) along with two idealized, non-dynamic models to demonstrate the framework and its use. Dynamic simulations with CISM are forced from 1991 to 2013, using combinations of reanalysis-based surface mass balance and observations of outlet glacier flux change. We propose and demonstrate qualitative and quantitative metrics for use in evaluating the different model simulations against the observations. We find that the altimetry observations used here are largely ambiguous in terms of their ability to distinguish one simulation from another. Based on basin-scale and whole-ice-sheet-scale metrics, we find that simulations using both idealized conceptual models and dynamic, numerical models provide an equally reasonable representation of the ice sheet surface (mean elevation differences of < 1 m). This is likely due to their short period of record, biases inherent to digital elevation models used for model initial conditions, and biases resulting from firn dynamics, which are not explicitly accounted for in the models or observations. On the other hand, we find that the gravimetry observations used here are able to unambiguously distinguish between simulations of varying complexity, and along with the CmCt, can provide a quantitative score for assessing a particular model and/or simulation. The new framework demonstrates that our proposed metrics can distinguish relatively better from relatively worse simulations and that dynamic ice sheet models, when appropriately initialized and forced with the right boundary conditions, demonstrate a predictive skill with respect to observed dynamic changes that have occurred on

An ice sheet model validation framework for the Greenland ice sheet

DOE Office of Scientific and Technical Information (OSTI.GOV)

Price, Stephen F.; Hoffman, Matthew J.; Bonin, Jennifer A.

We propose a new ice sheet model validation framework the Cryospheric Model Comparison Tool (CMCT) that takes advantage of ice sheet altimetry and gravimetry observations collected over the past several decades and is applied here to modeling of the Greenland ice sheet. We use realistic simulations performed with the Community Ice Sheet Model (CISM) along with two idealized, non-dynamic models to demonstrate the framework and its use. Dynamic simulations with CISM are forced from 1991 to 2013 using combinations of reanalysis-based surface mass balance and observations of outlet glacier flux change. We propose and demonstrate qualitative and quanti- tative metricsmore » for use in evaluating the different model simulations against the observations. We find 10 that the altimetry observations used here are largely ambiguous in terms of their ability to distinguish one simulation from another. Based on basin- and whole-ice-sheet scale metrics, the model initial condition as well as output from idealized and dynamic models all provide an equally reasonable representation of the ice sheet surface (mean elevation differences of <1 m). This is likely due to their short period of record, biases inherent to digital elevation models used for model initial conditions, and biases resulting from firn dynamics, which are not explicitly accounted for in the models or observations. On the other hand, we find that the gravimetry observations used here are able to unambiguously distinguish between simulations of varying complexity, and along with the CMCT, can provide a quantitative score for assessing a particular model and/or simulation. The new framework demonstrates that our proposed metrics can distinguish relatively better from relatively worse simulations and that dynamic ice sheet models, when appropriately initialized and forced with the right boundary conditions, demonstrate predictive skill with respect to observed dynamic changes occurring on Greenland over the past few

An ice sheet model validation framework for the Greenland ice sheet

PubMed Central

Price, Stephen F.; Hoffman, Matthew J.; Bonin, Jennifer A.; Howat, Ian M.; Neumann, Thomas; Saba, Jack; Tezaur, Irina; Guerber, Jeffrey; Chambers, Don P.; Evans, Katherine J.; Kennedy, Joseph H.; Lenaerts, Jan; Lipscomb, William H.; Perego, Mauro; Salinger, Andrew G.; Tuminaro, Raymond S.; van den Broeke, Michiel R.; Nowicki, Sophie M. J.

2018-01-01

We propose a new ice sheet model validation framework – the Cryospheric Model Comparison Tool (CmCt) – that takes advantage of ice sheet altimetry and gravimetry observations collected over the past several decades and is applied here to modeling of the Greenland ice sheet. We use realistic simulations performed with the Community Ice Sheet Model (CISM) along with two idealized, non-dynamic models to demonstrate the framework and its use. Dynamic simulations with CISM are forced from 1991 to 2013 using combinations of reanalysis-based surface mass balance and observations of outlet glacier flux change. We propose and demonstrate qualitative and quantitative metrics for use in evaluating the different model simulations against the observations. We find that the altimetry observations used here are largely ambiguous in terms of their ability to distinguish one simulation from another. Based on basin- and whole-ice-sheet scale metrics, we find that simulations using both idealized conceptual models and dynamic, numerical models provide an equally reasonable representation of the ice sheet surface (mean elevation differences of <1 m). This is likely due to their short period of record, biases inherent to digital elevation models used for model initial conditions, and biases resulting from firn dynamics, which are not explicitly accounted for in the models or observations. On the other hand, we find that the gravimetry observations used here are able to unambiguously distinguish between simulations of varying complexity, and along with the CmCt, can provide a quantitative score for assessing a particular model and/or simulation. The new framework demonstrates that our proposed metrics can distinguish relatively better from relatively worse simulations and that dynamic ice sheet models, when appropriately initialized and forced with the right boundary conditions, demonstrate predictive skill with respect to observed dynamic changes occurring on Greenland over the

An ice sheet model validation framework for the Greenland ice sheet.

PubMed

Price, Stephen F; Hoffman, Matthew J; Bonin, Jennifer A; Howat, Ian M; Neumann, Thomas; Saba, Jack; Tezaur, Irina; Guerber, Jeffrey; Chambers, Don P; Evans, Katherine J; Kennedy, Joseph H; Lenaerts, Jan; Lipscomb, William H; Perego, Mauro; Salinger, Andrew G; Tuminaro, Raymond S; van den Broeke, Michiel R; Nowicki, Sophie M J

2017-01-01

We propose a new ice sheet model validation framework - the Cryospheric Model Comparison Tool (CmCt) - that takes advantage of ice sheet altimetry and gravimetry observations collected over the past several decades and is applied here to modeling of the Greenland ice sheet. We use realistic simulations performed with the Community Ice Sheet Model (CISM) along with two idealized, non-dynamic models to demonstrate the framework and its use. Dynamic simulations with CISM are forced from 1991 to 2013 using combinations of reanalysis-based surface mass balance and observations of outlet glacier flux change. We propose and demonstrate qualitative and quantitative metrics for use in evaluating the different model simulations against the observations. We find that the altimetry observations used here are largely ambiguous in terms of their ability to distinguish one simulation from another. Based on basin- and whole-ice-sheet scale metrics, we find that simulations using both idealized conceptual models and dynamic, numerical models provide an equally reasonable representation of the ice sheet surface (mean elevation differences of <1 m). This is likely due to their short period of record, biases inherent to digital elevation models used for model initial conditions, and biases resulting from firn dynamics, which are not explicitly accounted for in the models or observations. On the other hand, we find that the gravimetry observations used here are able to unambiguously distinguish between simulations of varying complexity, and along with the CmCt, can provide a quantitative score for assessing a particular model and/or simulation. The new framework demonstrates that our proposed metrics can distinguish relatively better from relatively worse simulations and that dynamic ice sheet models, when appropriately initialized and forced with the right boundary conditions, demonstrate predictive skill with respect to observed dynamic changes occurring on Greenland over the past

An ice sheet model validation framework for the Greenland ice sheet

DOE PAGES

Price, Stephen F.; Hoffman, Matthew J.; Bonin, Jennifer A.; ...

2017-01-17

We propose a new ice sheet model validation framework the Cryospheric Model Comparison Tool (CMCT) that takes advantage of ice sheet altimetry and gravimetry observations collected over the past several decades and is applied here to modeling of the Greenland ice sheet. We use realistic simulations performed with the Community Ice Sheet Model (CISM) along with two idealized, non-dynamic models to demonstrate the framework and its use. Dynamic simulations with CISM are forced from 1991 to 2013 using combinations of reanalysis-based surface mass balance and observations of outlet glacier flux change. We propose and demonstrate qualitative and quanti- tative metricsmore » for use in evaluating the different model simulations against the observations. We find 10 that the altimetry observations used here are largely ambiguous in terms of their ability to distinguish one simulation from another. Based on basin- and whole-ice-sheet scale metrics, the model initial condition as well as output from idealized and dynamic models all provide an equally reasonable representation of the ice sheet surface (mean elevation differences of <1 m). This is likely due to their short period of record, biases inherent to digital elevation models used for model initial conditions, and biases resulting from firn dynamics, which are not explicitly accounted for in the models or observations. On the other hand, we find that the gravimetry observations used here are able to unambiguously distinguish between simulations of varying complexity, and along with the CMCT, can provide a quantitative score for assessing a particular model and/or simulation. The new framework demonstrates that our proposed metrics can distinguish relatively better from relatively worse simulations and that dynamic ice sheet models, when appropriately initialized and forced with the right boundary conditions, demonstrate predictive skill with respect to observed dynamic changes occurring on Greenland over the past few

An Ice Sheet Model Validation Framework for the Greenland Ice Sheet

NASA Technical Reports Server (NTRS)

Price, Stephen F.; Hoffman, Matthew J.; Bonin, Jennifer A.; Howat, Ian M.; Neumann, Thomas A.; Saba, Jack; Tezaur, Irina; Guerber, Jeffrey R.; Chambers, Don P.; Evans, Katherine J.;

IceBridge Provides Novel Evidence for Thick Units of Basal Freeze-on Ice Along Petermann Glacier, Greenland

NASA Astrophysics Data System (ADS)

Bell, R. E.; Tinto, K. J.; Wolovick, M.; Block, A. E.; Frearson, N.; Das, I.; Abdi, A.; Creyts, T. T.; Cochran, J. R.; Csatho, B. M.; Babonis, G. S.

2011-12-01

in the gravity data, and no such signature exists. We interpret these basal features in the radar data as ice. The geometry of these packages is very similar to the frozen-on ice packages imaged along the margins of the Gamburtsev Mountains, East Antarctica. We have interpreted these packages as bodies of accreted ice sourced from the upstream melting. Similar packages of basal ice up to 1200m thick are found throughout the margins of the Greenland ice sheet north of 70°N. The accretion process is modifying the base of the ice sheet in the onset regions of Petermann and other outlet glaciers in Greenland. Any change in rheology between basal and meteoric ice is likely to influence ice flow fostering enhanced melting and deformation of the basal ice. These frozen-on basal ice packages may be more susceptible to melt when ice sheet goes afloat and may be the origin of the elongate melt channels in the ice shelf. Accretion may influence both ice flow in outlet glaciers and melt rates at the grounding line and the adjacent ice shelves.

An Accretion Model for the Growth of Black Hole in Quasars

NASA Technical Reports Server (NTRS)

Lu, Ye; Cheng, K. S.; Zhang, S. N.

2003-01-01

A possible accretion model associated with the ionization instability of quasar disks is proposed to address the growth of the central black hole harbored in the host galaxy. The evolution of quasars in cosmic time is assumed to change from a highly active state to a quiescent state triggered by the S-shaped ionization instability of the quasar accretion disk. For a given external mass transfer rate ionization instability can modify accretion rate in the disk and separates the accretion flows of the disk into three different phases like a S-shape. We suggest that the bright quasars observed today are those quasars with disks in the upper branch of S-shaped instability and the dormant quasars are the system in the lower branch. The disk is assumed to evolve as ADIOS configuration in the lower branch. The mass ratio between black hole and its host galactic bulge is a nature consequence of ADIOS. Our model also demonstrates that a seed black hole 2 x 10(exp 6) solar masses similar to those found in spiral galaxies today is needed to produce a black hole with a final mass 2 x 10(exp 8) solar masses.

Inception of the Laurentide Ice Sheet using asynchronous coupling of a regional atmospheric model and an ice model

NASA Astrophysics Data System (ADS)

Birch, L.; Cronin, T.; Tziperman, E.

2017-12-01

The climate over the past 0.8 million years has been dominated by ice ages. Ice sheets have grown about every 100 kyrs, starting from warm interglacials, until they spanned continents. State-of-the-art global climate models (GCMs) have difficulty simulating glacial inception, or the transition of Earth's climate from an interglacial to a glacial state. It has been suggested that this failure may be related to their poorly resolved local mountain topography, due to their coarse spatial resolution. We examine this idea as well as the possible role of ice flow dynamics missing in GCMs. We investigate the growth of the Laurentide Ice Sheet at 115 kya by focusing on the mountain glaciers of Canada's Baffin Island, where geologic evidence indicates the last inception occurred. We use the Weather Research and Forecasting model (WRF) in a regional, cloud-resolving configuration with resolved mountain terrain to explore how quickly Baffin Island could become glaciated with the favorable yet realizable conditions of 115 kya insolation, cool summers, and wet winters. Using the model-derived mountain glacier mass balance, we force an ice sheet model based on the shallow-ice approximation, capturing the ice flow that may be critical to the spread of ice sheets away from mountain ice caps. The ice sheet model calculates the surface area newly covered by ice and the change in the ice surface elevation, which we then use to run WRF again. Through this type of iterated asynchronous coupling, we investigate how the regional climate responds to both larger areas of ice cover and changes in ice surface elevation. In addition, we use the NOAH-MP Land model to characterize the importance of land processes, like refreezing. We find that initial ice growth on the Penny Ice Cap causes regional cooling that increases the accumulation on the Barnes Ice Cap. We investigate how ice and topography changes on Baffin Island may impact both the regional climate and the large-scale circulation.

The effect of accretion environment at large radius on hot accretion flows

NASA Astrophysics Data System (ADS)

Yang, Xiao-Hong; Bu, De-Fu

2018-05-01

We study the effects of accretion environment (gas density, temperature, and angular momentum) at large radii (˜10 pc) on luminosity of hot accretion flows. The radiative feedback effects from the accretion flow on the accretion environment are also self-consistently taken into account. We find that the slowly rotating flows at large radii can significantly deviate from Bondi accretion when radiation heating and cooling are considered. We further find that when the temperature of environment gas is low (e.g. T = 2 × 107 K), the luminosity of hot accretion flows is high. When the temperature of gas is high (e.g. T ≥ 4 × 107 K), the luminosity of hot accretion flow significantly deceases. The environment gas density can also significantly influence the luminosity of accretion flows. When density is higher than ˜4 × 10-22 g cm-3 and temperature is lower than 2 × 107 K, hot accretion flow with luminosity lower than 2 per cent LEdd is not present. Therefore, the parsec-scale environment density and temperature are two important parameters to determine the luminosity. The results are also useful for the subgrid models adopted by the cosmological simulations.

Sensitivities of Greenland ice sheet volume inferred from an ice sheet adjoint model

NASA Astrophysics Data System (ADS)

Heimbach, P.; Bugnion, V.

2009-04-01

We present a new and original approach to understanding the sensitivity of the Greenland ice sheet to key model parameters and environmental conditions. At the heart of this approach is the use of an adjoint ice sheet model. Since its introduction by MacAyeal (1992), the adjoint method has become widespread to fit ice stream models to the increasing number and diversity of satellite observations, and to estimate uncertain model parameters such as basal conditions. However, no attempt has been made to extend this method to comprehensive ice sheet models. As a first step toward the use of adjoints of comprehensive three-dimensional ice sheet models we have generated an adjoint of the ice sheet model SICOPOLIS of Greve (1997). The adjoint was generated by means of the automatic differentiation (AD) tool TAF. The AD tool generates exact source code representing the tangent linear and adjoint model of the nonlinear parent model provided. Model sensitivities are given by the partial derivatives of a scalar-valued model diagnostic with respect to the controls, and can be efficiently calculated via the adjoint. By way of example, we determine the sensitivity of the total Greenland ice volume to various control variables, such as spatial fields of basal flow parameters, surface and basal forcings, and initial conditions. Reliability of the adjoint was tested through finite-difference perturbation calculations for various control variables and perturbation regions. Besides confirming qualitative aspects of ice sheet sensitivities, such as expected regional variations, we detect regions where model sensitivities are seemingly unexpected or counter-intuitive, albeit ``real'' in the sense of actual model behavior. An example is inferred regions where sensitivities of ice sheet volume to basal sliding coefficient are positive, i.e. where a local increase in basal sliding parameter increases the ice sheet volume. Similarly, positive ice temperature sensitivities in certain parts

Devon Ice cap's future: results from climate and ice dynamics modelling via surface mass balance modelling

NASA Astrophysics Data System (ADS)

Rodehacke, C. B.; Mottram, R.; Boberg, F.

2017-12-01

The Devon Ice Cap is an example of a relatively well monitored small ice cap in the Canadian Arctic. Close to Greenland, it shows a similar surface mass balance signal to glaciers in western Greenland. Here we various boundary conditions, ranging from ERA-Interim reanalysis data via global climate model high resolution (5km) output from the regional climate model HIRHAM5, to determine the surface mass balance of the Devon ice cap. These SMB estimates are used to drive the PISM glacier model in order to model the present day and future prospects of this small Arctic ice cap. Observational data from the Devon Ice Cap in Arctic Canada is used to evaluate the surface mass balance (SMB) data output from the HIRHAM5 model for simulations forced with the ERA-Interim climate reanalysis data and the historical emissions scenario run by the EC-Earth global climate model. The RCP8.5 scenario simulated by EC-Earth is also downscaled by HIRHAM5 and this output is used to force the PISM model to simulate the likely future evolution of the Devon Ice Cap under a warming climate. We find that the Devon Ice Cap is likely to continue its present day retreat, though in the future increased precipitation partly offsets the enhanced melt rates caused by climate change.

Modeling the response of a standard accretion disc to stochastic viscous fluctuations

NASA Astrophysics Data System (ADS)

Ahmad, Naveel; Misra, Ranjeev; Iqbal, Naseer; Maqbool, Bari; Hamid, Mubashir

2018-01-01

The observed variability of X-ray binaries over a wide range of time-scales can be understood in the framework of a stochastic propagation model, where viscous fluctuations at different radii induce accretion rate variability that propagate inwards to the X-ray producing region. The scenario successfully explains the power spectra, the linear rms-flux relation as well as the time-lag between different energy photons. The predictions of this model have been obtained using approximate analytical solutions or empirically motivated models which take into account the effect of these propagating variability on the radiative process of complex accretion flows. Here, we study the variation of the accretion rate due to such viscous fluctuations using a hydro-dynamical code for the standard geometrically thin, gas pressure dominated α-disc with a zero torque boundary condition. Our results confirm earlier findings that the time-lag between a perturbation and the resultant inner accretion rate variation depends on the frequency (or time-period) of the perturbation. Here we have quantified that the time-lag tlag ∝f-0.54 , for time-periods less than the viscous time-scale of the perturbation radius and is nearly constant otherwise. This, coupled with radiative process would produce the observed frequency dependent time-lag between different energy bands. We also confirm that if there are random Gaussian fluctuations of the α-parameter at different radii, the resultant inner accretion rate has a power spectrum which is a power-law.

Towards Improving Sea Ice Predictabiity: Evaluating Climate Models Against Satellite Sea Ice Observations

NASA Astrophysics Data System (ADS)

Stroeve, J. C.

2014-12-01

The last four decades have seen a remarkable decline in the spatial extent of the Arctic sea ice cover, presenting both challenges and opportunities to Arctic residents, government agencies and industry. After the record low extent in September 2007 effort has increased to improve seasonal, decadal-scale and longer-term predictions of the sea ice cover. Coupled global climate models (GCMs) consistently project that if greenhouse gas concentrations continue to rise, the eventual outcome will be a complete loss of the multiyear ice cover. However, confidence in these projections depends o HoHoweon the models ability to reproduce features of the present-day climate. Comparison between models participating in the World Climate Research Programme Coupled Model Intercomparison Project Phase 5 (CMIP5) and observations of sea ice extent and thickness show that (1) historical trends from 85% of the model ensemble members remain smaller than observed, and (2) spatial patterns of sea ice thickness are poorly represented in most models. Part of the explanation lies with a failure of models to represent details of the mean atmospheric circulation pattern that governs the transport and spatial distribution of sea ice. These results raise concerns regarding the ability of CMIP5 models to realistically represent the processes driving the decline of Arctic sea ice and to project the timing of when a seasonally ice-free Arctic may be realized. On shorter time-scales, seasonal sea ice prediction has been challenged to predict the sea ice extent from Arctic conditions a few months to a year in advance. Efforts such as the Sea Ice Outlook (SIO) project, originally organized through the Study of Environmental Change (SEARCH) and now managed by the Sea Ice Prediction Network project (SIPN) synthesize predictions of the September sea ice extent based on a variety of approaches, including heuristic, statistical and dynamical modeling. Analysis of SIO contributions reveals that when the

Icing Frequencies Experienced During Climb and Descent by Fighter-Interceptor Aircraft

NASA Technical Reports Server (NTRS)

Perkins, Porter J.

1958-01-01

Data and analyses are presented on the relative frequencies of occurrence and severity of icing cloud layers encountered by jet aircraft in the climb and descent phases of flights to high altitudes. Fighter-interceptor aircraft operated by the Air Defense Command (USAF) at bases in the Duluth and Seattle areas collected the data with icing meters installed for a l-year period. The project was part of an extensive program conducted by the NACA to collect Icing cloud data for evaluating the icing problem relevant to routine operations. The average frequency of occurrence of icing was found to be about 5 percent of the number of climbs and descents during 1 year of operations The icing encounters were predominantly in the low and middle cloud layers, decreasing above 15,000 feet to practically none above 25,000 feet. The greatest thickness of ice that would accumulate on any aircraft component (as indicated by the accretion on a small object) was measured with the icing meters. The ice thicknesses on a small sensing probe averaged less than 1/32 inch and did not exceed 1/2 inch. Such accumulations are relatively small when compared with those that can form during horizontal flight in icing clouds. The light accretions resulted from relatively steep angles of flight through generally thin cloud layers. Because of the limited statistical reliability of the results, an analysis was made using previous statistics on icing clouds below an altitude of 20,000 feet to determine the general icing severity probabilities. The calculations were made using adiabatic lifting as a basis to establish the liquid-water content. Probabilities of over-all ice accretions on a small object as a function of airspeed and rate of climb were computed from the derived water contents. These results were then combined with the probability of occurrence of icing in order to give the icing severity that can be expected for routine aircraft operations.

Ice formation in subglacial Lake Vostok, Central Antarctica

NASA Astrophysics Data System (ADS)

Souchez, R.; Petit, J. R.; Tison, J.-L.; Jouzel, J.; Verbeke, V.

2000-09-01

The investigation of chemical and isotopic properties in the lake ice from the Vostok ice core gives clues to the mechanisms involved in ice formation within the lake. A small lake water salinity can be reasonably deduced from the chemical data. Possible implications for the water circulation of Lake Vostok are developed. The characteristics of the isotopic composition of the lake ice indicate that ice formation in Lake Vostok occurred by frazil ice crystal generation due to supercooling as a consequence of rising waters and a possible contrast in water salinity. Subsequent consolidation of the developed loose ice crystals results in the accretion of ice to the ceiling of the lake.

Modeling of a Turbofan Engine with Ice Crystal Ingestion in the NASA Propulsion System Laboratory

NASA Technical Reports Server (NTRS)

Veres, Joseph P.; Jorgenson, Philip C. E.; Jones, Scott M.; Nili, Samaun

2017-01-01

The main focus of this study is to apply a computational tool for the flow analysis of the turbine engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) at NASA Glenn Research Center. The PSL has been used to test a highly instrumented Honeywell ALF502R-5A (LF11) turbofan engine at simulated altitude operating conditions. Test data analysis with an engine cycle code and a compressor flow code was conducted to determine the values of key icing parameters, that can indicate the risk of ice accretion, which can lead to engine rollback (un-commanded loss of engine thrust). The full engine aerothermodynamic performance was modeled with the Honeywell Customer Deck specifically created for the ALF502R-5A engine. The mean-line compressor flow analysis code, which includes a code that models the state of the ice crystal, was used to model the air flow through the fan-core and low pressure compressor. The results of the compressor flow analyses included calculations of the ice-water flow rate to air flow rate ratio (IWAR), the local static wet bulb temperature, and the particle melt ratio throughout the flow field. It was found that the assumed particle size had a large effect on the particle melt ratio, and on the local wet bulb temperature. In this study the particle size was varied parametrically to produce a non-zero calculated melt ratio in the exit guide vane (EGV) region of the low pressure compressor (LPC) for the data points that experienced a growth of blockage there, and a subsequent engine called rollback (CRB). At data points where the engine experienced a CRB having the lowest wet bulb temperature of 492 degrees Rankine at the EGV trailing edge, the smallest particle size that produced a non-zero melt ratio (between 3 percent - 4 percent) was on the order of 1 micron. This value of melt ratio was utilized as the target for all other subsequent data points analyzed, while the particle size was varied from 1 micron - 9

Numerical model of ice melange expansion during abrupt ice-shelf collapse

NASA Astrophysics Data System (ADS)

Guttenberg, N.; Abbot, D. S.; Amundson, J. M.; Burton, J. C.; Cathles, L. M.; Macayeal, D. R.; Zhang, W.

2010-12-01

Satellite imagery of the February 2008 Wilkins Ice-Shelf Collapse event reveals that a large percentage of the involved ice shelf was converted to capsized icebergs and broken fragments of icebergs over a relatively short period of time, possibly less than 24 hours. The extreme violence and short time scale of the event, and the considerable reduction of gravitational potential energy between upright and capsized icebergs, suggests that iceberg capsize might be an important driving mechanism controlling both the rate and spatial extent of ice shelf collapse. To investigate this suggestion, we have constructed an idealized, 2-dimensional model of a disintegrating ice shelf composed of a large number (N~100 to >1000) of initially well-packed icebergs of rectangular cross section. The model geometry consists of a longitudinal cross section of the idealized ice shelf from grounding line (or the upstream extent of ice-shelf fragmentation) to seaward ice front, and includes the region beyond the initial ice front to cover the open, ice-free water into which the collapsing ice shelf expands. The seawater in which the icebergs float is treated as a hydrostatic fluid in the computation of iceberg orientation (e.g., the evaluation of buoyancy forces and torques), thereby eliminating the complexities of free-surface waves, but net horizontal drift of the icebergs is resisted by a linear drag law designed to energy dissipation by viscous forces and surface-gravity-wave radiation. Icebergs interact via both elastic and inelastic contacts (typically a corner of one iceberg will scrape along the face of its neighbor). Ice-shelf collapse in the model is embodied by the mass capsize of a large proportion of the initially packed icebergs and the consequent advancement of the ice front (leading edge). Model simulations are conducted to examine (a) the threshold of stability (e.g., what density of initially capsizable icebergs is needed to allow a small perturbation to the system

Spin Evolution of Accreting Young Stars. II. Effect of Accretion-powered Stellar Winds

NASA Astrophysics Data System (ADS)

Matt, Sean P.; Pinzón, Giovanni; Greene, Thomas P.; Pudritz, Ralph E.

2012-01-01

We present a model for the rotational evolution of a young, solar-mass star interacting magnetically with an accretion disk. As in a previous paper (Paper I), the model includes changes in the star's mass and radius as it descends the Hayashi track, a decreasing accretion rate, and a prescription for the angular momentum transfer between the star and disk. Paper I concluded that, for the relatively strong magnetic coupling expected in real systems, additional processes are necessary to explain the existence of slowly rotating pre-main-sequence stars. In the present paper, we extend the stellar spin model to include the effect of a spin-down torque that arises from an accretion-powered stellar wind (APSW). For a range of magnetic field strengths, accretion rates, initial spin rates, and mass outflow rates, the modeled stars exhibit rotation periods within the range of 1-10 days in the age range of 1-3 Myr. This range coincides with the bulk of the observed rotation periods, with the slow rotators corresponding to stars with the lowest accretion rates, strongest magnetic fields, and/or highest stellar wind mass outflow rates. We also make a direct, quantitative comparison between the APSW scenario and the two types of disk-locking models (namely, the X-wind and Ghosh & Lamb type models) and identify some remaining theoretical issues for understanding young star spins.

Environmental controls on micro fracture processes in shelf ice

NASA Astrophysics Data System (ADS)

Sammonds, Peter

2013-04-01

The recent retreat and collapse of the ice shelves on the Antarctic Peninsula has been associated with regional atmospheric warming, oceanic warming, increased summer melt and shelf flexure. Although the cause of collapse is a matter of active discussion, the process is that of fracture of a creep-brittle material, close to its melting point. The environmental controls on how fracturing initiates, at a micro-scale, strongly determine the macroscopic disintegration of ice shelves. In particular the shelf temperature profile controls the plasticity of the ice shelf; the densification of shelf ice due to melting and re-freezing affects the crack tip stress intensity; the accretion of marine ice at the bottom of the shelf imposes a thermal/mechanical discontinuity; saline environments control crack tip stress corrosion; cyclic loading promotes sub-critical crack propagation. These strong environmental controls on shelf ice fracture means that assessing shelf stability is a non-deterministic problem. How these factors may be parameterized in ice shelf models, through the use of fracture mechanisms maps, is discussed. The findings are discussed in relation to the stability of Larsen C.

High Energy Neutrinos Produced in the Accretion Disks by Neutrons from Nuclei Disintegrated in the AGN Jets

NASA Astrophysics Data System (ADS)

Bednarek, W.

2016-12-01

We investigate the consequences of acceleration of nuclei in jets of active galaxies not far from the surface of an accretion disk. The nuclei can be accelerated in the re-connection regions in the jet and/or at the jet boundary, between the relativistic jet and its cocoon. It is shown that the relativistic nuclei can efficiently fragment onto specific nucleons in collisions with the disk radiation. Neutrons, directed toward the accretion disk, take a significant part of energy from the relativistic nuclei. These neutrons develop a cascade in the dense accretion disk. We calculate the neutrino spectra produced in such a hadronic cascade within the accretion disk. We propose that the neutrinos produced in such a scenario, from the whole population of super-massive black holes in active galaxies, can explain the extragalactic neutrino background recently measured by the IceCube neutrino detector, provided that a 5% fraction of galaxies have an active galactic nucleus and a few percent of neutrons reach the accretion disk. We predict that the neutrino signals in the present neutrino detectors, produced in terms of such a model, will not be detectable even from the nearby radio galaxies similar to M87.

Sensitivity of open-water ice growth and ice concentration evolution in a coupled atmosphere-ocean-sea ice model

NASA Astrophysics Data System (ADS)

Shi, Xiaoxu; Lohmann, Gerrit

2017-09-01

A coupled atmosphere-ocean-sea ice model is applied to investigate to what degree the area-thickness distribution of new ice formed in open water affects the ice and ocean properties. Two sensitivity experiments are performed which modify the horizontal-to-vertical aspect ratio of open-water ice growth. The resulting changes in the Arctic sea-ice concentration strongly affect the surface albedo, the ocean heat release to the atmosphere, and the sea-ice production. The changes are further amplified through a positive feedback mechanism among the Arctic sea ice, the Atlantic Meridional Overturning Circulation (AMOC), and the surface air temperature in the Arctic, as the Fram Strait sea ice import influences the freshwater budget in the North Atlantic Ocean. Anomalies in sea-ice transport lead to changes in sea surface properties of the North Atlantic and the strength of AMOC. For the Southern Ocean, the most pronounced change is a warming along the Antarctic Circumpolar Current (ACC), owing to the interhemispheric bipolar seasaw linked to AMOC weakening. Another insight of this study lies on the improvement of our climate model. The ocean component FESOM is a newly developed ocean-sea ice model with an unstructured mesh and multi-resolution. We find that the subpolar sea-ice boundary in the Northern Hemisphere can be improved by tuning the process of open-water ice growth, which strongly influences the sea ice concentration in the marginal ice zone, the North Atlantic circulation, salinity and Arctic sea ice volume. Since the distribution of new ice on open water relies on many uncertain parameters and the knowledge of the detailed processes is currently too crude, it is a challenge to implement the processes realistically into models. Based on our sensitivity experiments, we conclude a pronounced uncertainty related to open-water sea ice growth which could significantly affect the climate system sensitivity.

Projecting Antarctic ice discharge using response functions from SeaRISE ice-sheet models

NASA Astrophysics Data System (ADS)

Levermann, A.; Winkelmann, R.; Nowicki, S.; Fastook, J. L.; Frieler, K.; Greve, R.; Hellmer, H. H.; Martin, M. A.; Meinshausen, M.; Mengel, M.; Payne, A. J.; Pollard, D.; Sato, T.; Timmermann, R.; Wang, W. L.; Bindschadler, R. A.

2014-08-01

The largest uncertainty in projections of future sea-level change results from the potentially changing dynamical ice discharge from Antarctica. Basal ice-shelf melting induced by a warming ocean has been identified as a major cause for additional ice flow across the grounding line. Here we attempt to estimate the uncertainty range of future ice discharge from Antarctica by combining uncertainty in the climatic forcing, the oceanic response and the ice-sheet model response. The uncertainty in the global mean temperature increase is obtained from historically constrained emulations with the MAGICC-6.0 (Model for the Assessment of Greenhouse gas Induced Climate Change) model. The oceanic forcing is derived from scaling of the subsurface with the atmospheric warming from 19 comprehensive climate models of the Coupled Model Intercomparison Project (CMIP-5) and two ocean models from the EU-project Ice2Sea. The dynamic ice-sheet response is derived from linear response functions for basal ice-shelf melting for four different Antarctic drainage regions using experiments from the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. The resulting uncertainty range for the historic Antarctic contribution to global sea-level rise from 1992 to 2011 agrees with the observed contribution for this period if we use the three ice-sheet models with an explicit representation of ice-shelf dynamics and account for the time-delayed warming of the oceanic subsurface compared to the surface air temperature. The median of the additional ice loss for the 21st century is computed to 0.07 m (66% range: 0.02-0.14 m; 90% range: 0.0-0.23 m) of global sea-level equivalent for the low-emission RCP-2.6 (Representative Concentration Pathway) scenario and 0.09 m (66% range: 0.04-0.21 m; 90% range: 0.01-0.37 m) for the strongest RCP-8.5. Assuming no time delay between the atmospheric warming and the oceanic subsurface, these

Sensitivity of an Antarctic Ice Sheet Model to Sub-Ice-Shelf Melting

NASA Astrophysics Data System (ADS)

Lipscomb, W. H.; Leguy, G.; Urban, N. M.; Berdahl, M.

2017-12-01

Theory and observations suggest that marine-based sectors of the Antarctic ice sheet could retreat rapidly under ocean warming and increased melting beneath ice shelves. Numerical models of marine ice sheets vary widely in sensitivity, depending on grid resolution and the parameterization of key processes (e.g., calving and hydrofracture). Here we study the sensitivity of the Antarctic ice sheet to ocean warming and sub-shelf melting in standalone simulations of the Community Ice Sheet Model (CISM). Melt rates either are prescribed based on observations and high-resolution ocean model output, or are derived from a plume model forced by idealized ocean temperature profiles. In CISM, we vary the model resolution (between 1 and 8 km), Stokes approximation (shallow-shelf, depth-integrated higher-order, or 3D higher-order) and calving scheme to create an ensemble of plausible responses to sub-shelf melting. This work supports a broader goal of building statistical and reduced models that can translate large-scale Earth-system model projections to changes in Antarctic ocean temperatures and ice sheet discharge, thus better quantifying uncertainty in Antarctic-sourced sea-level rise.

Experimental study of performance degradation of a rotating system in the NASA Lewis RC icing tunnel

NASA Technical Reports Server (NTRS)

Korkan, Kenneth

1992-01-01

The Helicopter Icing Consortium (HIC) conducted one of the first U.S. tests of a heavily instrumented model in the controlled environment of a refrigerated tunnel. In the Icing Research Tunnel (IRT) at NASA LeRC, ice was accreted on the main rotor blade of the BMTR-1 Sikorsky model helicopter under a variety of environmental conditions, such that liquid water content (LWC) and volume mean droplet diameter (VMD) ranges reflected the Federal Aviation Agency and Department of Defence icing condition envelopes. This report gives the correlated results of the data provided by NASA LeRC. The method of statistical analysis is discussed. Lift, thrust, and torque coefficients are presented as a function of icing time, as correlated with changes in ambient temperature, LWC, and VMD. The physical significance of these forces is discussed.

Modeling Europa's Ice-Ocean Interface

NASA Astrophysics Data System (ADS)

Elsenousy, A.; Vance, S.; Bills, B. G.

2014-12-01

This work focuses on modeling the ice-ocean interface on Jupiter's Moon (Europa); mainly from the standpoint of heat and salt transfer relationship with emphasis on the basal ice growth rate and its implications to Europa's tidal response. Modeling the heat and salt flux at Europa's ice/ocean interface is necessary to understand the dynamics of Europa's ocean and its interaction with the upper ice shell as well as the history of active turbulence at this area. To achieve this goal, we used McPhee et al., 2008 parameterizations on Earth's ice/ocean interface that was developed to meet Europa's ocean dynamics. We varied one parameter at a time to test its influence on both; "h" the basal ice growth rate and on "R" the double diffusion tendency strength. The double diffusion tendency "R" was calculated as the ratio between the interface heat exchange coefficient αh to the interface salt exchange coefficient αs. Our preliminary results showed a strong double diffusion tendency R ~200 at Europa's ice-ocean interface for plausible changes in the heat flux due to onset or elimination of a hydrothermal activity, suggesting supercooling and a strong tendency for forming frazil ice.

Design, fabrication, and testing of an ultrasonic de-icing system for helicopter rotor blades

NASA Astrophysics Data System (ADS)

Palacios, Jose Luis

A low-power, non-thermal ultrasonic de-icing system is introduced as a possible substitute for current electro-thermal systems. The system generates delaminating ultrasonic transverse shear stresses at the interface of accreted ice. A PZT-4 disk driven at 28.5 KHz (radial resonance of the disk) instantaneously de-bonds 2 mm thick freezer ice layers. The ice layers are accreted to a 0.7 mm thick, 30.4 cm x 30.4 cm steel plate at an environment temperature of -20°C. A power input of 50 Watts is applied to the actuator (50 V, 19.6 KV/m), which translates to a de-icing power of 0.07 W/cm2. A finite element model of the actuator bonded to the isotropic plate is used to guide the design of the system, and predicts the transverse shear stresses at the ice interface. Wind tunnel icing tests were conducted to demonstrate the potential use of the proposed system under impact icing conditions. Both glaze ice and rime ice were generated on steel and composite plates by changing the cloud conditions of the wind tunnel. Continuous ultrasonic vibration prevented impact ice formation around the actuator location at an input power not exceeding 0.18 W/cm 2 (1.2 W/in2). As ice thickness reached a critical thickness of approximately 1.2 mm, shedding occurred on those locations where ultrasonic transverse shear stresses exceeded the shear adhesion strength of the ice. Finite element transverse shear stress predictions correlate with observed experimental impact ice de-bonding behavior. To increase the traveling distance of propagating ultrasonic waves, ultrasonic shear horizontal wave modes are studied. Wave modes providing large modal interface transverse shear stress concentration coefficients (ISCC) between the host structure (0.7 mm thick steel plate) and accreted ice (2.5 mm thick ice layer) are identified and investigated for a potential increase in the wave propagation distance. Ultrasonic actuators able to trigger these optimum wave modes are designed and fabricated. Despite

Inclusion of TCAF model in XSPEC to study accretion flow dynamics around black hole candidates

NASA Astrophysics Data System (ADS)

Debnath, Dipak; Chakrabarti, Sandip Kumar; Mondal, Santanu

Spectral and Temporal properties of black hole candidates can be well understood with the Chakrabarti-Titarchuk solution of two component advective flow (TCAF). This model requires two accretion rates, namely, the Keplerian disk accretion rate and the sub-Keplerian halo accretion rate, the latter being composed of a low angular momentum flow which may or may not develop a shock. In this solution, the relevant parameter is the relative importance of the halo (which creates the Compton cloud region) rate with respect to the Keplerian disk rate (soft photon source). Though this model has been used earlier to manually fit data of several black hole candidates quite satisfactorily, for the first time we are able to create a user friendly version by implementing additive Table model FITS file into GSFC/NASA's spectral analysis software package XSPEC. This enables any user to extract physical parameters of accretion flows, such as two accretion rates, shock location, shock strength etc. for any black hole candidate. Most importantly, unlike any other theoretical model, we show that TCAF is capable of predicting timing properties from spectral fits, since in TCAF, a shock is responsible for deciding spectral slopes as well as QPO frequencies.

On the wind production from hot accretion flows with different accretion rates

NASA Astrophysics Data System (ADS)

Bu, De-Fu; Gan, Zhao-Ming

2018-02-01

We perform two-dimensional simulations to study how the wind strength changes with accretion rate. We take into account bremsstrahlung, synchrotron radiation and the Comptonization. We find that when the accretion rate is low, radiative cooling is not important, and the accretion flow is hot. For the hot accretion flow, wind is very strong. The mass flux of wind can be ˜ 50 per cent of the mass inflow rate. When the accretion rate increases to a value at which radiative cooling rate is roughly equal to or slightly larger than viscous heating rate, cold clumps can form around the equatorial plane. In this case, the gas pressure gradient force is small and wind is very weak. Our results may be useful for the sub-grid model of active galactic nuclear feedback study.

Variational Ridging in Sea Ice Models

NASA Astrophysics Data System (ADS)

Roberts, A.; Hunke, E. C.; Lipscomb, W. H.; Maslowski, W.; Kamal, S.

2017-12-01

This work presents the results of a new development to make basin-scale sea ice models aware of the shape, porosity and extent of individual ridges within the pack. We have derived an analytic solution for the Euler-Lagrange equation of individual ridges that accounts for non-conservative forces, and therefore the compressive strength of individual ridges. Because a region of the pack is simply a collection of paths of individual ridges, we are able to solve the Euler-Lagrange equation for a large-scale sea ice field also, and therefore the compressive strength of a region of the pack that explicitly accounts for the macro-porosity of ridged debris. We make a number of assumptions that have simplified the problem, such as treating sea ice as a granular material in ridges, and assuming that bending moments associated with ridging are perturbations around an isostatic state. Regardless of these simplifications, the ridge model is remarkably predictive of macro-porosity and ridge shape, and, because our equations are analytic, they do not require costly computations to solve the Euler-Lagrange equation of ridges on the large scale. The new ridge model is therefore applicable to large-scale sea ice models. We present results from this theoretical development, as well as plans to apply it to the Regional Arctic System Model and a community sea ice code. Most importantly, the new ridging model is particularly useful for pinpointing gaps in our observational record of sea ice ridges, and points to the need for improved measurements of the evolution of porosity of deformed ice in the Arctic and Antarctic. Such knowledge is not only useful for improving models, but also for improving estimates of sea ice volume derived from altimetric measurements of sea ice freeboard.

Detached-Eddy Simulations of Separated Flow Around Wings With Ice Accretions: Year One Report

NASA Technical Reports Server (NTRS)

Choo, Yung K. (Technical Monitor); Thompson, David; Mogili, Prasad

2004-01-01

A computational investigation was performed to assess the effectiveness of Detached-Eddy Simulation (DES) as a tool for predicting icing effects. The AVUS code, a public domain flow solver, was employed to compute solutions for an iced wing configuration using DES and steady Reynolds Averaged Navier-Stokes (RANS) equation methodologies. The configuration was an extruded GLC305/944-ice shape section with a rectangular planform. The model was mounted between two walls so no tip effects were considered. The numerical results were validated by comparison with experimental data for the same configuration. The time-averaged DES computations showed some improvement in lift and drag results near stall when compared to steady RANS results. However, comparisons of the flow field details did not show the level of agreement suggested by the integrated quantities. Based on our results, we believe that DES may prove useful in a limited sense to provide analysis of iced wing configurations when there is significant flow separation, e.g., near stall, where steady RANS computations are demonstrably ineffective. However, more validation is needed to determine what role DES can play as part of an overall icing effects prediction strategy. We conclude the report with an assessment of existing computational tools for application to the iced wing problem and a discussion of issues that merit further study.

SPIN EVOLUTION OF ACCRETING YOUNG STARS. II. EFFECT OF ACCRETION-POWERED STELLAR WINDS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Matt, Sean P.; Pinzon, Giovanni; Greene, Thomas P.

2012-01-20

We present a model for the rotational evolution of a young, solar-mass star interacting magnetically with an accretion disk. As in a previous paper (Paper I), the model includes changes in the star's mass and radius as it descends the Hayashi track, a decreasing accretion rate, and a prescription for the angular momentum transfer between the star and disk. Paper I concluded that, for the relatively strong magnetic coupling expected in real systems, additional processes are necessary to explain the existence of slowly rotating pre-main-sequence stars. In the present paper, we extend the stellar spin model to include the effectmore » of a spin-down torque that arises from an accretion-powered stellar wind (APSW). For a range of magnetic field strengths, accretion rates, initial spin rates, and mass outflow rates, the modeled stars exhibit rotation periods within the range of 1-10 days in the age range of 1-3 Myr. This range coincides with the bulk of the observed rotation periods, with the slow rotators corresponding to stars with the lowest accretion rates, strongest magnetic fields, and/or highest stellar wind mass outflow rates. We also make a direct, quantitative comparison between the APSW scenario and the two types of disk-locking models (namely, the X-wind and Ghosh and Lamb type models) and identify some remaining theoretical issues for understanding young star spins.« less

Constraint on the black hole spin of M87 from the accretion-jet model

NASA Astrophysics Data System (ADS)

Feng, Jianchao; Wu, Qingwen

2017-09-01

The millimetre bump, as found in high-resolution multiwaveband observations of M87 by Prieto et al., most possibly comes from the synchrotron emission of thermal electrons in advection-dominated accretion flow (ADAF). It is possible to constrain the accretion rate near the horizon if both the nuclear millimetre emission and its polarization are produced by the hot plasma in the accretion flow. The jet power of M87 has been extensively explored, which is around 8_-3^{+7}× 10^{42} erg s-1 based on the analysis of the X-ray cavity. The black hole (BH) spin can be estimated if the jet power and the accretion rate near the horizon are known. We model the multiwavelength spectral energy distribution (SED) of M87 with a coupled ADAF-jet model surrounding a Kerr BH, where the full set of relativistic hydrodynamical equations of the ADAF are solved. The hybrid jet formation model, as a variant of the Blandford-Znajek model, is used to model the jet power. We find that the SMBH should be fast rotating with a dimensionless spin parameter a_{*}˜eq 0.98_-0.02^{+0.012}.

The response of grounded ice to ocean temperature forcing in a coupled ice sheet-ice shelf-ocean cavity model

NASA Astrophysics Data System (ADS)

Goldberg, D. N.; Little, C. M.; Sergienko, O. V.; Gnanadesikan, A.

2010-12-01

Ice shelves provide a pathway for the heat content of the ocean to influence continental ice sheets. Changes in the rate or location of basal melting can alter their geometry and effect changes in stress conditions at the grounding line, leading to a grounded ice response. Recent observations of ice streams and ice shelves in the Amundsen Sea sector of West Antarctica have been consistent with this story. On the other hand, ice dynamics in the grounding zone control flux into the shelf and thus ice shelf geometry, which has a strong influence on the circulation in the cavity beneath the shelf. Thus the coupling between the two systems, ocean and ice sheet-ice shelf, can be quite strong. We examine the response of the ice sheet-ice shelf-ocean cavity system to changes in ocean temperature using a recently developed coupled model. The coupled model consists a 3-D ocean model (GFDL's Generalized Ocean Layered Dynamics model, or GOLD) to a two-dimensional ice sheet-ice shelf model (Goldberg et al, 2009), and allows for changing cavity geometry and a migrating grounding line. Steady states of the coupled system are found even under considerable forcing. The ice shelf morphology and basal melt rate patterns of the steady states exhibit detailed structure, and furthermore seem to be unique and robust. The relationship between temperature forcing and area-averaged melt rate is influenced by the response of ice shelf morphology to thermal forcing, and is found to be sublinear in the range of forcing considered. However, results suggest that area-averaged melt rate is not the best predictor of overall system response, as grounding line stability depends on local aspects of the basal melt field. Goldberg, D N, D M Holland and C G Schoof, 2009. Grounding line movement and ice shelf buttressing in marine ice sheets, Journal of Geophysical Research-Earth Surfaces, 114, F04026.

Wave-Ice interaction in the Marginal Ice Zone: Toward a Wave-Ocean-Ice Coupled Modeling System

DTIC Science & Technology

2015-09-30

MIZ using WW3 (3 frequency bins, ice retreat in August and ice advance in October); Blue (solid): Based on observations near Antarctica by Meylan...1 DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Wave- Ice interaction in the Marginal Ice Zone: Toward a...Wave-Ocean- Ice Coupled Modeling System W. E. Rogers Naval Research Laboratory, Code 7322 Stennis Space Center, MS 39529 phone: (228) 688-4727

Ice Particle Impacts on a Flat Plate

NASA Technical Reports Server (NTRS)

Vargas, Mario; Ruggeri, Charles; Struk, Peter M.; Pereira, Mike; Revilock, Duane; Kreeger, Richard E.

2015-01-01

An experimental study was conducted at the Ballistic Laboratory of NASA Glenn Research Center to study the impact of ice particles on a stationary flat surface target set at 45 degrees with respect to the direction of motion of the impinging particle (Figure 1). The experiment is part of NASA efforts to study the physics involved in engine power-loss events due to ice-crystal ingestion and ice accretion formation inside engines. These events can occur when aircraft encounter high-altitude convective weather.

Further Evaluation of Scaling Methods for Rotorcraft Icing

NASA Technical Reports Server (NTRS)

Tsao, Jen-Ching; Kreeger, Richard E.

2012-01-01

The paper will present experimental results from two recent icing tests in the NASA Glenn Icing Research Tunnel (IRT). The first test, conducted in February 2009, was to evaluate the current recommended scaling methods for fixed wing on representative rotor airfoils at fixed angle of attack. For this test, scaling was based on the modified Ruff method with scale velocity determined by constant Weber number and water film Weber number. Models were un-swept NACA 0012 wing sections. The reference model had a chord of 91.4 cm and scale model had a chord of 35.6 cm. Reference tests were conducted with velocity of 100 kt (52 m/s), droplet medium volume diameter (MVD) 195 m, and stagnation-point freezing fractions of 0.3 and 0.5 at angle of attack of 5deg and 7deg . It was shown that good ice shape scaling was achieved with constant Weber number for NACA 0012 airfoils with angle of attack up to 7deg . The second test, completed in May 2010, was primarily focused on obtaining transient and steady-state iced aerodynamics, ice accretion and shedding, and thermal icing validation data from an oscillating airfoil section over some selected ranges of icing conditions and blade assembly operational configurations. The model used was a 38.1-cm chord Sikorsky SC2110 airfoil section installed on an airfoil test apparatus with oscillating capability in the IRT. For two test conditions, size and condition scaling were performed. It was shown that good ice shape scaling was achieved for SC2110 airfoil at dynamic pitching motion. The data obtained will be applicable for future main rotor blade and tail rotor blade applications.

Accretion onto CO White Dwarfs using MESA

NASA Astrophysics Data System (ADS)

Feng, Wanda; Starrfield, Sumner

2018-06-01

The nature of type Ia Supernovae (SNe Ia) progenitor systems and their underlying mechanism are not well understood. There are two competing progenitor scenarios: the single-degenerate scenario wherein a white dwarf (WD) star accretes material from a companion star, reaching the Chandrasekhar mass limit; and, the double-degenerate scenario wherein two WDs merge. In this study, we investigate the single-degenerate scenario by accretion onto carbon-oxygen (CO) WDs using the Modules for Experiments in Stellar Astrophysics (MESA). We vary the WD mass, composition of the accreting material, and accretion rate in our models. Mixing between the accreted material and the WD core is informed by multidimensional studies that suggest occurance after thermonuclear runaway (TNR) ensues. We compare the accretion of solar composition material onto CO WDs with the accretion of mixed solar and core material after TNR. As many of our models eject less material than accreted, our study supports that accretion onto CO WDs is a feasible channel for SNe I progenitors.

Arctic Sea Ice Classification and Mapping for Surface Albedo Parameterization in Sea Ice Modeling

NASA Astrophysics Data System (ADS)

Nghiem, S. V.; Clemente-Colón, P.; Perovich, D. K.; Polashenski, C.; Simpson, W. R.; Rigor, I. G.; Woods, J. E.; Nguyen, D. T.; Neumann, G.

2016-12-01

A regime shift of Arctic sea ice from predominantly perennial sea ice (multi-year ice or MYI) to seasonal sea ice (first-year ice or FYI) has occurred in recent decades. This shift has profoundly altered the proportional composition of different sea ice classes and the surface albedo distribution pertaining to each sea ice class. Such changes impacts physical, chemical, and biological processes in the Arctic atmosphere-ice-ocean system. The drastic changes upset the traditional geophysical representation of surface albedo of the Arctic sea ice cover in current models. A critical science issue is that these profound changes must be rigorously and systematically observed and characterized to enable a transformative re-parameterization of key model inputs, such as ice surface albedo, to ice-ocean-atmosphere climate modeling in order to obtain re-analyses that accurately reproduce Arctic changes and also to improve sea ice and weather forecast models. Addressing this challenge is a strategy identified by the National Research Council study on "Seasonal to Decadal Predictions of Arctic Sea Ice - Challenges and Strategies" to replicate the new Arctic reality. We review results of albedo characteristics associated with different sea ice classes such as FYI and MYI. Then we demonstrate the capability for sea ice classification and mapping using algorithms developed by the Jet Propulsion Laboratory and by the U.S. National Ice Center for use with multi-sourced satellite radar data at L, C, and Ku bands. Results obtained with independent algorithms for different radar frequencies consistently identify sea ice classes and thereby cross-verify the sea ice classification methods. Moreover, field observations obtained from buoy webcams and along an extensive trek across Elson Lagoon and a sector of the Beaufort Sea during the BRomine, Ozone, and Mercury EXperiment (BROMEX) in March 2012 are used to validate satellite products of sea ice classes. This research enables the mapping

Sea Ice Biogeochemistry: A Guide for Modellers

PubMed Central

Tedesco, Letizia; Vichi, Marcello

2014-01-01

Sea ice is a fundamental component of the climate system and plays a key role in polar trophic food webs. Nonetheless sea ice biogeochemical dynamics at large temporal and spatial scales are still rarely described. Numerical models may potentially contribute integrating among sparse observations, but available models of sea ice biogeochemistry are still scarce, whether their relevance for properly describing the current and future state of the polar oceans has been recently addressed. A general methodology to develop a sea ice biogeochemical model is presented, deriving it from an existing validated model application by extension of generic pelagic biogeochemistry model parameterizations. The described methodology is flexible and considers different levels of ecosystem complexity and vertical representation, while adopting a strategy of coupling that ensures mass conservation. We show how to apply this methodology step by step by building an intermediate complexity model from a published realistic application and applying it to analyze theoretically a typical season of first-year sea ice in the Arctic, the one currently needing the most urgent understanding. The aim is to (1) introduce sea ice biogeochemistry and address its relevance to ocean modelers of polar regions, supporting them in adding a new sea ice component to their modelling framework for a more adequate representation of the sea ice-covered ocean ecosystem as a whole, and (2) extend our knowledge on the relevant controlling factors of sea ice algal production, showing that beyond the light and nutrient availability, the duration of the sea ice season may play a key-role shaping the algal production during the on going and upcoming projected changes. PMID:24586604

Multiscale Models of Melting Arctic Sea Ice

DTIC Science & Technology

2013-09-30

September 29, 2013 LONG-TERM GOALS Sea ice reflectance or albedo , a key parameter in climate modeling, is primarily determined by melt pond...and ice floe configurations. Ice - albedo feedback has played a major role in the recent declines of the summer Arctic sea ice pack. However...understanding the evolution of melt ponds and sea ice albedo remains a significant challenge to improving climate models. Our research is focused on

Changes in the morphology of interstellar ice analogues after hydrogen atom exposure.

PubMed

Accolla, Mario; Congiu, Emanuele; Dulieu, François; Manicò, Giulio; Chaabouni, Henda; Matar, Elie; Mokrane, Hakima; Lemaire, Jean Louis; Pirronello, Valerio

2011-05-07

The morphology of water ice in the interstellar medium is still an open question. Although accretion of gaseous water could not be the only possible origin of the observed icy mantles covering dust grains in cold molecular clouds, it is well known that water accreted from the gas phase on surfaces kept at 10 K forms ice films that exhibit a very high porosity. It is also known that in the dark clouds H(2) formation occurs on the icy surface of dust grains and that part of the energy (4.48 eV) released when adsorbed atoms react to form H(2) is deposited in the ice. The experimental study described in the present work focuses on how relevant changes of the ice morphology result from atomic hydrogen exposure and subsequent recombination. Using the temperature-programmed desorption (TPD) technique and a method of inversion analysis of TPD spectra, we show that there is an exponential decrease in the porosity of the amorphous water ice sample following D-atom irradiation. This decrease is inversely proportional to the thickness of the ice and has a value of ϕ(0) = 2 × 10(16) D-atoms cm(-2) per layer of H(2)O. We also use a model which confirms that the binding sites on the porous ice are destroyed regardless of their energy depth, and that the reduction of the porosity corresponds in fact to a reduction of the effective area. This reduction appears to be compatible with the fraction of D(2) formation energy transferred to the porous ice network. Under interstellar conditions, this effect is likely to be efficient and, together with other compaction processes, provides a good argument to believe that interstellar ice is amorphous and non-porous. This journal is © the Owner Societies 2011