FENSAP-ICE: Unsteady Conjugate Heat Transfer Simulation of Electrothermal De-Icing

DOI: 10.2514/1.C031607 This paper presents a truly unsteady approach for the numerical simulation of in-flight electrothermal anti-icing orde-icing, using a conjugate heat transfer technique. This numerical approach has been implemented in FENSAPICE to compute the complex heat transfer phenomena occurring during in-flight de-icing with multiple heating elements following independent cycling. At each time step, the energy fluxes through the aircraft’s solid skin, the melting ice layer, the liquid water film, and the external fluid are computed. The ice shape is then modified by taking intoaccounttheopposingmassbalanceeffectsoficeaccretingduetotheimpactofsupercooleddropletsand/orwater runback, and the partial or total melting of the existing ice layer due to heating. The results of the verification of this phase-changeconductioncodearepresented,followedbyastudyofintercyclede-icingonawing,showingintercycle ice growth. Nomenclature ch = convective heat flux coefficient d = droplet diameter E = internal energy e = volumetric internal energy H = enthalpy