Search Results

Numerically predicted roughness distributions obtained in in-flight icing simulations with a beading model are used in a quasi-steady manner to compute ice shapes. This approach, called "Multishot," uses a number of steady flow and droplet solutions for computing short intervals (shots) of the total ice accretion time. The iced geometry, the grid, and the surface roughness distribution are updated after each shot, producing a better match with the unsteady ice accretion phenomena. Comparisons to multishot results with uniform roughness show that the evolution of the surface roughness distribution has a strong influence on the final ice shape. The ice horns that form are longer and thinner compared to constant roughness results. The constant roughness approach usually fails to capture the formation of the pressure side horns and under-predicts the thickness of the ice in this region.

The irregular shapes that glaze ice may grow into while accreting over the surface of an aircraft represent a major difficulty in the numerical simulation of long periods of in-flight icing. In the framework of Arbitrary Lagrangian-Eulerian (ALE) formulations, a mesh movement scheme is presented, in which frame and elasticity analogies are loosely coupled. The resulting deformed mesh preserves the quality of elements, especially in the near-wall region, where accurate prediction of heat flux and shear stresses are required. The proposed scheme handles mesh movement in a computationally efficient manner by localizing the mesh deformation. Numerical results of ice shapes and the corresponding aerodynamic coefficients are compared with the experimental results.

The high computational cost of 3-D viscous turbulent aero-icing simulations is one of the main limitations to address in order to more extensively use computational fluid dynamics to explore the wide variety of icing conditions to be tested before achieving aircraft airworthiness. In an attempt to overcome the computational burden of these simulations, a Reduced Order Modeling (ROM) approach, based on Proper Orthogonal Decomposition (POD) and Kriging interpolation techniques, is applied to the computation of the impingement pattern of supercooled large droplets (SLD) on aircraft. Relying on a suitable database of high fidelity full-order simulations, the ROM approach provides a lower-order approximation of the system in terms of a linear combination of appropriate functions. The accuracy of the resulting surrogate solution is successfully compared to experimental and CFD results for sample 2-D problems and then extended to a typical 3-D case.