This content is not included in your SAE MOBILUS subscription, or you are not logged in.
FENSAP-ICE: A CFD Monte Carlo Approach to Shed-Ice Trajectory and Impact
ISSN: 0148-7191, e-ISSN: 2688-3627
Published June 13, 2011 by SAE International in United States
Annotation ability available
A fully CFD-based methodology for ice particle tracking based on a Monte Carlo statistical approach and a six-degrees-of-freedom particle-tracking module has been developed within the FENSAP-ICE in-flight icing system. A one-way aerodynamic coupling between the airflow and the ice particle has been adopted, such that the flowfield determines the forces and moments on the particle at each location on its track, but the particle, being much smaller, has no aerodynamic effect on the aircraft's flowfield. A complete envelope of force and moment coefficients has been computed for a representative ice shape, in order to generate a permanent database. At each time step during the integration of the particle track, the angles of the local flow velocity vector with the principal axes of the particle are determined and used to interpolate the corresponding force and moment coefficients from the particle's database. These 6-DOFs are then used to compute the next particle location. The linear displacement is obtained from the integration of Newton's Second Law, whereas angular displacement is obtained from the numerical integration of Euler's equations of motion using quaternions and the 4-stage Runge-Kutta algorithm. A comparison of the aerodynamic force and moment envelope for one ice shape with experimental data will be presented to demonstrate the reliability of this approach. Shed ice particle trajectory computation is demonstrated for a complete aircraft geometry.
CitationDeschenes, J., Baruzzi, G., Lagace, P., and Habashi, W., "FENSAP-ICE: A CFD Monte Carlo Approach to Shed-Ice Trajectory and Impact," SAE Technical Paper 2011-38-0089, 2011, https://doi.org/10.4271/2011-38-0089.
- Murman, S.M., Aftmosmis, M.J. and Berger, M.J., Simulations of 6-DOF Motion with a Cartesian Method. AIAA Paper 2003-1246, 41st AIAA Aerospace Sciences Meeting, January 6-9, 2003, Reno, Nevada.
- Prewitt, N.C., Belk, D.M. and Shyy, W., Parallel Computing of Overset Grids for Aerodynamic Problems with Moving Objects. Progress in Aerospace Sciences, vol. 36, 2000, pp. 117-172.
- Baruzzi, G.S., Lagacé, P., Aubé, M.S. and Habashi, W.G. Development of a Shed-Ice Trajectory Simulation in FENSAP-ICE. SAE Paper 07ICE34,. SAE Aircraft and Engine Icing International Conference, Seville, Spain, 2007.
- Chandrasekharan, R., Hinson, M., “Trajectory Simulation of Ice Shed from a Business Jet,” SAE Technical Paper 2003-01-3032, 2003, doi:10.4271/2003-01-3032.
- Jacob, J. and Papadakis, M. Experimental and Computational Investigation of Ice Shedding from Aircraft Surfaces. Master Thesis, Wichita: Wichita State University, 2006.
- Papadakis, M., Yeong, H.-W. and Suares, I.G. Simulation of Ice Shedding from a Business Jet Aircraft. AIAA Paper 2007-506, 45th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 2007.
- Barber, C.J. Qhull, The geometry Center, University of Minnesota, 2010.
- Kuipers, J.B. Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality, Princeton University Press, 2002.