Conjugate Heat Transfer and Aerothermal Wall Modeling for In-Flight Icing Conditions

In-flight ice accretion severely degrades the aerodynamics of lifting surfaces andis closely tied to the wall heat transfer that governs phase change at the fluid–iceinterface. This paper presents a conjugate heat transfer (CHT) methodology thatcouples a 2D compressible RANS solver (k–ω SST) on an evolving unstructured meshwith the solution of heat conduction within the growing ice layer, whose consistencyis checked by a 1D thermal model with adaptive accretion rate. Applied to a NACA-0012 airfoil (M∞ = 0.5, Re = 8.96 × 106) over a 300 s exposure, the simulationpredicts a maximum ice thickness of 9.27 mm, leading to a 32.4 % drop in maximumlift coefficient (CL,max) and a 2.9 % increase in drag at zero angle of attack. This worklays the groundwork for the future integration of electrothermal de-icing systems andfor experimental validation in an icing wind tunnel.

Authors

Publication Details

Journal
Zenodo (CERN European Organization for Nuclear Research)
Published
2026-10-26
DOI
https://doi.org/10.5281/zenodo.23660850
Primary Topic
Icing and De-icing Technologies
Type
preprint
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
preprint

Conjugate Heat Transfer and Aerothermal Wall Modeling for In-Flight Icing Conditions

Valentin Noël Maurice GITTA, Mathieu Daniel Pierre Granger
Zenodo (CERN European Organization for Nuclear Research)
Icing and De-icing Technologies
preprint

Conjugate Heat Transfer and Aerothermal Wall Modeling for In-Flight Icing Conditions

Valentin Noël Maurice GITTA, Mathieu Daniel Pierre Granger
preprint en

Abstract

In-flight ice accretion severely degrades the aerodynamics of lifting surfaces andis closely tied to the wall heat transfer that governs phase change at the fluid–iceinterface. This paper presents a conjugate heat transfer (CHT) methodology thatcouples a 2D compressible RANS solver (k–ω SST) on an evolving unstructured meshwith the solution of heat conduction within the growing ice layer, whose consistencyis checked by a 1D thermal model with adaptive accretion rate. Applied to a NACA-0012 airfoil (M∞ = 0.5, Re = 8.96 × 106) over a 300 s exposure, the simulationpredicts a maximum ice thickness of 9.27 mm, leading to a 32.4 % drop in maximumlift coefficient (CL,max) and a 2.9 % increase in drag at zero angle of attack. This worklays the groundwork for the future integration of electrothermal de-icing systems andfor experimental validation in an icing wind tunnel.

Zenodo (CERN European Organization for Nuclear Research)
Icing and De-icing Technologies
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.

Laika access coming soon