An automatic multi-stepping approach to aircraft ice prediction

Flying an aircraft in icing conditions may seriously degrade its aerodynamical performance and threaten the flight safety. Over the years, new technologies and improved procedures have limited the potential risks caused by aircraft icing. Experimental studies being very expensive, numerous computer codes have been developed to simulate ice shapes and tackle the problem. Typically in these codes, a flow solution and key icing parameters are evaluated around a clean un-iced geometry and their values remain constant during the entire simulation. This approach may be acceptable for short exposure times or when the ice shape only slightly deforms the initial geometry. However, in other cases, the values of the icing parameters may vary and the simulation will loose its accuracy: for large shapes, the presence of the ice influences the surrounding airflow significantly, altering the value of icing parameters and ultimately the ice accretion. Calculating more accurate ice shapes therefore requires to periodically recompute the flow field around the body during the simulation and determine updated values for icing parameters. This procedure, known as multi-stepping, is investigated in this thesis and adapted to the new threedimensional icing code ICECREMO2. Several multi-step algorithms are presented and tested on cylinders and airfoils. When possible, the ice shapes simulated are compared with experimental results. The first multi-step calculations were generally performed manually. The user had to perform a rather tedious work and inappropriate instructions could lead to severe inaccuracies in the simulations. To avoid these difficulties, a fully automated procedure will be developed including all stages of a multi-step computation. This significantly reduces user interaction and the overall computing time. The present research work forms part of the ICECREMO2 project. ICECREMO2 is a three-dimensional ice accretion and water flow code developed collaboratively by Airbus UK, BAe Systems, Dunlop Aerospace, Rolls-Royce, GKN Westland Helicopters, QinetiQ and Cranfield University under the auspices of the UK Department of Trade and Industry.

[1]  Tim G. Myers,et al.  Thin Films with High Surface Tension , 1998, SIAM Rev..

[2]  Michael B. Bragg,et al.  Unsteady aerodynamic measurements on an iced airfoil , 2002 .

[3]  W. G. Cleveland,et al.  On the evolution of a turbulent boundary layer induced by a three-dimensional roughness element , 1992, Journal of Fluid Mechanics.

[4]  Wagdi G. Habashi,et al.  3D thermal analysis of an anti-icing device using FENSAP-ICE , 1998 .

[5]  R. J. Hansman,et al.  Deterministic multi-zone ice accretion modeling , 1991 .

[6]  Steven C. Chapra,et al.  Numerical Methods for Engineers , 1986 .

[7]  T. Myers.  Extension to the Messinger Model for Aircraft Icing , 2001 .

[8]  Charles Pereira.  Status of NTSB aircraft icing certification-related safety recommendations issued as a result of the 1994 ATR-72 accident at Roselawn, IN , 1997 .

[9]  W. Olsen,et al.  Experimental evidence for modifying the current physical model for ice accretion on aircraft surfaces , 1986 .

[10]  P. Spalart.  A One-Equation Turbulence Model for Aerodynamic Flows , 1992 .

[11]  Droplet impact on water layers: post-impact analysis and computations , 2005, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[12]  Thomas P. Ratvasky,et al.  Freezing Rain as an In-Flight Icing Hazard , 2000 .

[13]  S. A. Sherif,et al.  A semi-empirical model for heat transfer and ice accretion on aircraft wings in supercooled clouds , 1997 .

[14]  Mark G. Potapczuk,et al.  LEWICE/E: An Euler based ice accretion code , 1992 .

[15]  T. Hedde,et al.  Improvement of the ONERA 3D icing code, comparison with 3D experimental shapes , 1993 .

[16]  M. Velazquez,et al.  Ice accretion and performance degradation calculations with LEWICE/NS , 1993 .

[17]  A. Rothmayer,et al.  Transient effects in predicting bulk surface properties for aircraft ice accretion , 2001 .

[18]  Steven Caruso.  NEARICE - An unstructured-mesh Navier-Stokes-based ice accretion prediction method , 1994 .

[19]  C. Maliska,et al.  NUMERICAL SIMULATION OF ICE ACCRETION ON THE LEADING EDGE OF AERODYNAMIC PROFILES , 2001 .

[20]  Tim G. Myers,et al.  The flow and solidification of a thin fluid film on an arbitrary three-dimensional surface , 2002 .

[21]  Thomas J. Vukits.  Overview and Risk Assessment of Icing for Transport Category Aircraft and Components , 2002 .

[22]  H. Ockendon,et al.  Slow viscous flow , 1995 .

[23]  William B. Wright,et al.  User Manual for the NASA Glenn Ice Accretion Code LEWICE: Version 2.0 , 1999 .

[24]  Gary A. Ruff,et al.  Users Manual for the NASA Lewis Ice Accretion Prediction Code (LEWICE) , 1990 .

[25]  G. Naterer.  Dispersed multiphase flow with air-driven runback of a liquid layer at a moving boundary , 2003 .

[26]  B. L. Messinger.  Equilibrium Temperature of an Unheated Icing Surface as a Function of Air Speed , 1953 .

[27]  P. Havugimana,et al.  A Comparison of Skin Friction and Heat Transfer Prediction by Various Roughness Models , 2002 .

[28]  Michael B. Bragg,et al.  Experimental Investigation of Simulated Large-Droplet Ice Shapes on Airfoil Aerodynamics , 1999 .

[29]  E. White,et al.  Experiments on Surface Roughness Effects in Ice Accretion , 2005 .

[30]  Paul C. Ivey,et al.  Analysis and experimental aspects of the impact of supercooled water droplets into thin water films , 2005 .

[31]  A. Rothmayer,et al.  A mechanism for ice roughness formation on an airfoil leading edge, contributing to glaze ice accretion , 1998 .

[32]  H. Fernholz.  Boundary Layer Theory , 2001 .

[33]  R. Hansman,et al.  Modeling of surface roughness effects on glaze ice accretion , 1989 .

[34]  Patrick G. Verdin,et al.  Explicit Finite Volume Modeling of Aircraft Anti-Icing and De-Icing , 2008 .

[35]  Robert P. Cassoni,et al.  Aircraft Anti-Icing and De-Icing Techniques and Modeling , 1996 .

[36]  W. H. Allen.  Dictionary of Technical Terms for Aerospace Use , 1965 .

[37]  E. Lozowski,et al.  The Icing of an Unheated, Nonrotating Cylinder. Part II. Icing Wind Tunnel Experiments , 1983 .

[38]  Yves Soufflet,et al.  Improvement of runback water film calculation and its impact on ice prediction , 2008 .

[39]  W B Wright,et al.  DRA/NASA/ONERA COLLABORATION ON ICING RESEARCH, PART 2 , 1997 .

[40]  R. Hansman,et al.  Heat transfer on accreting ice surfaces , 1990 .

[41]  Adam Rutkowski,et al.  Validation Results for LEWICE 2.0. [Supplement] , 1999 .

[42]  Tim G. Myers,et al.  Slowly accreting ice due to supercooled water impacting on a cold surface , 2002 .

[43]  Michael B. Bragg,et al.  Heat-transfer and freestream turbulence measurements for improvement of the ice accretion physical model , 1997 .

[44]  R. J. Kind,et al.  Experimental and computational simulation of in-flight icing phenomena , 1998 .

[45]  Andy P. Broeren,et al.  Iced-airfoil aerodynamics , 2005 .

[46]  W Sproles Darrell,et al.  Computer Program To Obtain Ordinates for NACA Airfoils , 1996 .

[47]  Wagdi G. Habashi,et al.  Development of a Second Generation In-Flight Icing Simulation Code , 2006 .

[48]  Scott K. Thomas,et al.  Aircraft anti-icing and deicing techniques and modeling , 1996 .

[49]  S. Bankoff,et al.  Long-scale evolution of thin liquid films , 1997 .

[50]  Frank Lynch,et al.  Effects of ice accretions on aircraft aerodynamics , 2001 .

[51]  Jaiwon Shin,et al.  Experimental and computational ice shapes and resulting drag increase for a NACA 0012 airfoil , 1992 .

[52]  E. Lozowski,et al.  Discrete Modelling of Ice Accretion Shape and Structure Forming on a Cylinder Under In-Flight Icing Conditions , 2003 .

[53]  Lasse Makkonen,et al.  Models for the growth of rime, glaze, icicles and wet snow on structures , 2000, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[54]  R. Hansman,et al.  Investigation of surface water behavior during glaze ice accretion , 1989 .

[55]  Roy Rasmussen,et al.  Winter Icing and Storms Project (WISP). , 1992 .

[56]  Mario Vargas,et al.  Close-up analysis of aircraft ice accretion , 1993 .

[57]  M. Papadakis,et al.  Experimental Study of Large Droplet Splashing and Breakup , 2007 .

[58]  William B. Wright,et al.  Users manual for the improved NASA Lewis ice accretion code LEWICE 1.6 , 1995 .

[59]  William B. Wright,et al.  Computational Simulation of Large Droplet Icing , 1996 .

[60]  Michael B. Bragg,et al.  AERODYNAMICS OF SUPERCOOLED-LARGE-DROPLET ICE ACCRETIONS AND THE EFFECT ON AIRCRAFT CONTROL , 1996 .

[61]  D. W. Hammond,et al.  Ice and water film growth from incoming supercooled droplets , 1999 .

[62]  Bharat K. Soni,et al.  ICEG2D: A Software Package for Ice Accretion Prediction , 2003 .

[63]  Giuseppe Mingione,et al.  Ice Accretion Prediction on Multielement Airfoils , 1998 .

[64]  On the efficient use of predictor-corrector methods in the numerical solution of differential equations , 1965 .

[65]  M. B. Bragg,et al.  Aerodynamic measurements of an airfoil with simulated glaze ice , 1986 .

[66]  Richard Purvis,et al.  Analysis of super-cooled water droplet impact on a thin water layer and ice growth , 2006 .

[67]  Roger Gent,et al.  Results from SLD Mass Loss Tests in the ACT Luton Icing Research Wind Tunnel , 2003 .

[68]  R. W. Gent,et al.  Aircraft icing , 2000, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[69]  Thomas H. Bond,et al.  Results of an icing test on a NACA 0012 airfoil in the NASA Lewis Icing Research Tunnel , 1992 .

[70]  Harold E. Addy,et al.  A Numerical Evaluation of Icing Effects on a Natural Laminar Flow Airfoil , 2000 .

[71]  Michael B. Bragg,et al.  Effect of ice accretion on aircraft flight dynamics , 2000 .

[72]  François Morency,et al.  HOT AIR ANTI-ICING SYSTEM MODELIZATION IN THE ICE PREDICTION CODE CANICE , 1998 .

[73]  S. Singh,et al.  Determination of the local heat transfer coefficients of three simulated smooth ice formation characteristics on a cylinder , 1985 .

[74]  Ion Paraschivoiu,et al.  Ice accretion on aircraft wings with thermodynamic effects , 1994 .

[75]  Theo G. Keith,et al.  Overview of numerical codes developed for predicted electrothermal deicing of aircraft blades , 1988 .

[76]  Michael B. Bragg,et al.  Computational Investigation of Simulated Large-Droplet Ice Shapes on Airfoil Aerodynamics , 1999 .

[77]  B. Šarler.  Stefan's work on solid-liquid phase changes , 1995 .

[78]  Brian M. Berkowitz,et al.  An experimental investigation of multi-element airfoil ice accretion and resulting performance degradation , 1989 .

[79]  G. Mingione,et al.  Ice accretion prediction on multi-element airfoils , 1997 .

[80]  C. Thompson,et al.  Automatic Multi-Stepping Approach for Ice Predictions , 2004 .

[81]  Peter M. Render,et al.  Certification and Operation of Helicopters in Icing Environments , 1998 .

[82]  T. Hedde,et al.  Development of a three-dimensional icing code - Comparison with experimental shapes , 1992 .

[83]  A. Rothmayer,et al.  Thin liquid films flowing over external aerodynamic surfaces , 2002 .

[84]  Mark G. Potapczuk,et al.  Swept wing ice accretion modeling , 1990 .

[85]  Farooq Saeed,et al.  Surface Heat Transfer Study for Ice Accretion and Anti-Icing Prediction in Three Dimension , 2004 .

[86]  M. G. Potapczuk,et al.  Experimental investigation of multielement airfoil ice accretion and resulting performance degradation , 1990 .

[87]  David N. Anderson,et al.  Ice Accretion with Varying Surface Tension , 1995 .