Ice accretion on the inlet of an aero-engine could adversely affect the characteristics of the flow field into the engine and degrade the engine's performance. When the accrued ice sheds, the portion sucked into the engine could induce serious damages. One method to reduce hazards caused by in-flight icing is an anti-icing system that utilizes hot lubricating oil from the engine rather than utilizes the hot gas. This method not only decreases the amount of anti-icing hot air required but also cools the lubrication oil. This paper presents a computational study of the temperature distribution of an aero-engine strut under icing condition. The flow field around the strut is simulated using ANSYS FLUENT CFD code. The trajectories of the super-cooled water droplets and the local collection coefficient are calculated through the Lagrangian approach after the computation of the air flow field. The coupling effects of heat transfer and mass transfer are considered in the temperature calculation of the strut studied. The thermal model takes both the mass balance of water and the energy balance on the surface of the strut into account. The convection heat transfer coefficients on the surface of the strut are obtained through two methods: correlation equations and CFD calculations. Finally, comparisons between the computational results and the experimental data are presented.
NomenclatureA = area, m 2 C D = drag coefficient D = leading edge diameter, m D r = drag force, N G v = gravity force, N I = evaporation latent heat, J/kg LWC = liquid water content, g/m 3 MVD = mean volumn diameter, µm Nu = Nusselt number P = pressure, Pa P v = saturated vapor pressure, Pa 2 Q = heat flux, W/m 2 Re = Reynolds number S = stream length T = temperature, K θ = angle degree c p = specific heat, J/(kg·K) d = droplet diameter, µm h = convective heat transfer coefficient, W/(m 2 ·K) k = conduction heat transfer coefficient, W/(m·K) m = mass, kg m & = mass flow rate, kg/s n = water droplet number r = radius, m v r = velocity vector, m/s β = local collection coefficient µ = coefficient of dynamic viscosity, N·s/m 2 ρ = density, kg/m 3 τ = time, s Subscripts S = stream length a = air d = water droplet e = external edge of boundary layer evap = evaporation ht = impingement surface area imp = droplets impingement in = inlet of control volume m = impinging location of droplet o = original location of droplet out = outlet of control volume rel = relative s = solid t = total w = water wall = solid wall ∞ = free stream