Search Results

This paper studies the level of confidence and applicability of CFD simulations using steady-state Reynolds-Averaged Navier-Stokes (RANS) in predicting aerodynamic performance losses on swept-wings due to contamination with ice accreted in-flight. The wing geometry selected for the study is the 65%-scale Common Research Model (CRM65) main wing, for which NASA Glenn Research Center’s Icing Research Tunnel has generated experimental ice shapes for the inboard, mid-span, and outboard sections. The reproductions at various levels of fidelity from detailed 3D scans of these ice shapes have been used in recent aerodynamic testing at the Office National d’Etudes et Recherches Aérospatiales (ONERA) and Wichita State University (WSU) wind tunnels. The ONERA tests were at higher Reynolds number range in the order of 10 million, while the WSU tests were in the order of 1 million.

Two main approaches are available when studying droplet dynamics for in-flight icing simulations: the Lagrangian approach, in which each droplet trajectory is integrated until it impacts the vehicle under study or when it leaves it behind without impact, and the Eulerian approach, where the droplet dynamics is solved as a continuum. In both cases, the same momentum equations are solved. Each approach has its advantages. In 2D, the Lagrangian approach is easy to code and it is very efficient, particularly when used in combination with a panel method flow solver. However, it is a far less practical approach for 3D simulations, particularly on complex geometries, as it is not an easy task to accurately determine the droplet seeding region without a great number of droplet trajectories, dramatically increasing the computing cost. Converting the impact locations into a water collection distribution is also a complex task, since droplet trajectories in 3D can follow convoluted paths.

3D ice accretion codes have been available for a few decades but, depending on the specific application, their use may be cumbersome, time consuming and requiring a great deal of expertise in using the code. In particular, simulations of large 3D glaze ice accretions using multiple layers of ice is a very challenging and time consuming task. There are several reasons why 2D icing simulations tools are still widely used in the aircraft industry to produce realistic glaze ice shapes. 2D codes are very fast and robust, with a very short turn-around time. They produce adequate results in areas of the aircraft where 3D effects on airflow or droplets concentration can be neglected. Their use can be extended to other areas of the aircraft if relevant 3D effects can be taken into account. This paper proposes a simulation methodology that includes three levels of corrections to extend the use of 2D icing codes to most of the aircraft surfaces.

There are many different numerical approaches to ice accretion simulation. Little comparison has been made between those approaches to identify the best tool for a given application. This paper presents a comparison exercise between 2D codes (CANICE-BA and LEWICE) and 3D codes (CANICE3D-BA, LEWICE3D and FENSAP-ICE). It also compares the 3D first generation code (panel method with Lagrangian droplet trajectory tracking) CANICE3D-BA to the 3D second generation code (Navier-Stokes with Eulerian droplet tracking) FENSAP-ICE. The paper includes a description of the different methodologies. The first comparison exercise is done using three 2D cases for which experimental ice shapes are available. The second exercise addresses a water collection efficiency over an isolated tail for which experimental data is available. Finally, an ice accretion comparison is presented in a DLR4 wing-body configuration.

This paper presents CHT2D, a 2D hot air anti-icing simulation tool developed by the Advanced Aerodynamics group of Bombardier Aerospace. The tool has been developed from two main modules: the ice prediction code CANICE and the Navier-Stokes solver NSU2D, which is used to solve the hot air internal flow. A “weak” coupling beween the two modules based on function calls and information exchange has been priviledged. Three validation test cases are presented: for dry air conditions. Predictions from CHT2D agree quite well with the experiments. Preliminary results are also presented for a test case in icing conditions for different heat loads from the anti-icing system, to study the effect on the accumulated ice.