Ice accretion simulation on multi-element airfoils using extended Messinger model

In the current article, the problem of in-flight ice accumulation on multi-element airfoils is studied numerically. The analysis starts with flow field computation using the Hess-Smith panel method. The second step is the calculation of droplet trajectories and droplet collection efficiencies. In the next step, convective heat transfer coefficient distributions around the airfoil elements are calculated using the Integral Boundary-Layer Method. The formulation accounts for the surface roughness due to ice accretion. The fourth step consists of establishing the thermodynamic balance and computing ice accretion rates using the Extended Messinger Model. At low temperatures and low liquid water contents, rime ice occurs for which the ice shape is determined by a simple mass balance. At warmer temperatures and high liquid water contents, glaze ice forms for which the energy and mass conservation equations are combined to yield a single first order ordinary differential equation, solved numerically. Predicted ice shapes are compared with experimental shapes reported in the literature and good agreement is observed both for rime and glaze ice. Ice shapes and masses are also computed for realistic flight scenarios. The results indicate that the smaller elements in multielement configurations accumulate comparable and often greater amount of ice compared to larger elements. The results also indicate that the multi-layer approach yields more accurate results compared to the one-layer approach, especially for glaze ice conditions.     

Rime- and glaze-ice accretion due to freezing rain falling vertically on a horizontal thermally insulated overhead line conductor

A theory of atmospheric icing due to freezing rain on an overhead line conductor (OHLC) is developed. The rain falls vertically on a horizontal OHLC that is thermally insulated. It is assumed that the collection efficiency of the accretion surface is unity and that this surface is in thermodynamic equilibrium with the environment.

For air temperature T_A 0°C and raindrop temperature T_D 0°C, the freezing rain accretes as rime ice, provided that the temperature of the ice surface T_l < 0°C. The evolution equation governing the mass transfer at the accretion surface is solved analytically, yielding the shape of the rime-ice surface. Equations governing the thermal state of the rime-ice deposit are also given. These determine the onset of wet growth or glaze accretion at the upper stagnation line during suitable environmental conditions.

For environmental conditions producing an ice surface at temperature T_l = 0°gC, the freezing accretes as glaze. Equations governing the heat and mass transfer at the surface determine the shape of the glaze surface and the downward viscous motion of the unfrozen water. For T_D < 0°C, glaze evolution equations are developed for T_A 0°C and T_A 0°C. Analytical solutions of these equations are obtained. In particular, when T_D < −T_A < 0°C, the evolution equation predicts a novel limiting growth that is triangular in shape. Further study of the mass and heat transfer conditions, in the neighborhood of this final stage of glaze accretion, shows that it is maintained in thermodynamic equilibrium with its warm air environment.     

Modeling of Ice Accretion on the Airfoil Surface in an Air Flow Containing Ice Particles

The process of airfoil icing caused by incidence of ice crystals is considered. A physicomathematical model of motion of spheroidal crystals in the gas flow and their interaction with the body is formulated. The model takes into account the non-spherical particle orientation with respect to the velocity vector of the gas flow. It is assumed that the particle impact onto the body surface leads to partial destruction of the particle under the action of normal and tangential stresses, and some part of the particle mass remains in the vicinity of the impact point. An inviscid flow around the airfoil with a time-dependent shape is calculated by the method of approximate conformal mapping.     

Buoyancy effects on mixed convection heat and mass transfer in a duct with sudden expansions

The present paper deals with the study of heat and mass transfer by mixed convection in an inclined duct preceded with a double step expansion. The control volume based finite element method was used to solve the set of non-dimensional equations for the vorticity, stream function, energy and species conservation. Numerical simulations are carried out for different combinations of the Lewis number, thermal and mass diffusion Grashof numbers for different inclinations. Streamlines, temperature and concentration distributions are presented and discussed. The results show the effect of the secondary flow induced by buoyancy forces and the presence of the double step expansion on the heat and mass transfer mechanism. It is found that the recirculation vortices induced by the expansion can be present along the channel and the flow structure can be wavy. For the vertical orientation, asymmetric fields are observed for the different simulated cases.     

An icing physics study by using lifetime-based molecular tagging thermometry technique

An experimental investigation was conducted to quantify the unsteady heat transfer and phase changing process within small icing water droplets in order to elucidate underlying physics to improve our understanding of the important micro-physical process of icing phenomena. A novel, lifetime-based molecular tagging thermometry (MTT) technique was developed and implemented to achieve temporally-and-spatially resolved temperature distribution measurements to reveal the time evolution of the unsteady heat transfer and dynamic phase changing process within micro-sized water droplets in the course of icing process. It was found that, after a water droplet impinged onto a frozen cold surface, the liquid water at the bottom of the droplet would be frozen and turned to solid ice rapidly, while the upper portion of the droplet was still in liquid state. As the time goes by, the interface between the liquid phase water and solid phase ice was found to move upward continuously with more and more liquid water within the droplet turned to solid ice. Interestingly, the averaged temperature of the remaining liquid water within the small icing droplet was found to increase, rather than decrease, continuously in the course of icing process. The temperature increase of the remaining liquid water is believed to be due to the heat release of the latent heat during solidification process. The volume expansion of the water droplet during the icing process was found to be mainly upward to cause droplet height growth rather than radial to enlarge the contact area of the droplet on the test plate. As a result, the spherical-cap-shaped water droplet was found to turn to a prolate-spheroid-shaped ice crystal with cusp-like top at the end of the icing process. The required freezing time for the water droplets to turn to ice crystals completely was found to depend on the surface temperature of the test plate strongly, which would decrease exponentially as the surface temperature of the frozen cold test plate decreases.     

Mixed convection effects on heat and mass transfer in a non Newtonian fluid with chemical reaction over a vertical plate

This paper studies mixed convection, double dispersion and chemical reaction effects on heat and mass transfer in a non-Darcy non-Newtonian fluid over a vertical surface in a porous medium under the constant temperature and concentration. The governing boundary layer equations, namely, momentum, energy and concentration, are converted to ordinary differential equations by introducing similarity variables and then are solved numerically by means of fourth-order Runge-Kutta method coupled with double-shooting technique. The velocity, temperature concentration, heat and mass transfer profiles are presented graphically for various values of the parameters, and the influence of viscosity index n, thermal and solute dispersion, chemical reaction parameter χ are observed.     

Melting heat transfer in steady laminar flow over a moving surface

The steady laminar boundary layer flow and heat transfer from a warm, laminar liquid flow to a melting surface moving parallel to a constant free stream is studied in this paper. The continuity, momentum and energy equations, which are coupled nonlinear partial differential equations are reduced to a set of two nonlinear ordinary differential equations, before being solved numerically using the Runge–Kutta–Fehlberg method. Results for the skin friction coefficient, local Nusselt number, velocity profiles as well as temperature profiles are presented for different values of the governing parameters. Effects of the melting parameter, moving parameter and Prandtl number on the flow and heat transfer characteristics are thoroughly examined. It is found that the problem admits dual solutions.     

Lie group analysis for thermophoretic and radiative augmentation of heat and mass transfer in a Brinkman-Darcy flow over a flat surface with heat generation

A boundary layer analysis is presented to investigate numerically the effects of radiation,thermophoresis and the dimensionless heat generation or absorption on hydromagnetic flow with heat and mass transfer over a flat surface in a porous medium.The boundary layer equations are transformed to non-linear ordinary differential equations using scaling group of transformations and they are solved numerically by using the fourth order Runge-Kutta method with shooting technique for some values of physical parameters.Comparisons with previously published work are performed and the results are found to be in very good agreement.Many results are obtained and a representative set is displayed graphically to illustrate the influence of the various parameters on the dimensionless velocity,temperature and concentration profiles as well as the local skin-friction coefficient,wall heat transfer,particle deposition rate and wall thermophoretic deposition velocity.The results show that the magnetic field induces acceleration of the flow,rather than deceleration(as in classical magnetohydrodynamics(MHD) boundary layer flow) but to reduce temperature and increase concentration of particles in boundary layer.Also,there is a strong dependency of the concentration in the boundary layer on both the Schmidt number and mass transfer parameter.     

The effect of radiation on mixed convection flow past a stretching surface

A boundary layer analysis is presented for a study of the influence of radiation and buoyancy on heat and mass transfer characteristics of continuous surfaces having a prescribed variable surface temperature and stretched with rapidly decreasing power law velocities. The effects of suction in the presence of a quiescent fluid medium of constant temperature are considered. Rosseland approximation is used to describe the radiative heat flux in the energy equation. The transformed governing equations are solved numerically and the velocity and temperature profiles as well as the local Nusselt number and skin friction coefficient are presented. Results show that the effect of radiation is to keep the molten mass away from the slot warmer, reduces the friction factor and increases the heat transfer rate compared to the case with no radiation.     

Heat transfer coefficient measurements on the pressure surface of a transonic airfoil

This paper presents steady-state recovery temperature and heat transfer coefficient measurements on the pressure surface of a modern, highly cambered transonic airfoil. These measurements were collected with a peak Mach number of 1.5 and a maximum turbulence intensity of 30%. We used a single passage model to simulate the idealized two-dimensional flow path between rotor blades in a modern transonic turbine. This set up offered a simpler construction than a linear cascade, yet produced an equivalent flow condition. We performed validated high accuracy (±0.2°C) surface temperature measurements using wide-band thermochromic liquid crystals allowing separate measurements of the previously listed parameters with the same heat transfer surface. We achieved maximum heat transfer coefficient uncertainties that were equivalent to similar investigations (±10%). Two key observations are the heat transfer coefficient along the aft portion of the airfoil is sensitive to the surface heat flux and is highly insensitive to the level of freestream turbulence. Possible explanations for these observations are discussed.     

Comparison of interfacial heat transfer coefficient estimated by two different techniques during solidification of cylindrical aluminum alloy casting

The interfacial heat transfer coefficient (IHTC) is necessary for accurate simulation of the casting process. In this study, a cylindrical geometry is selected for the determination of the IHTC between aluminum alloy casting and the surrounding sand mold. The mold surface heat flux and temperature are estimated by two inverse heat conduction techniques, namely Beck’s algorithm and control volume technique. The instantaneous cast and mold temperatures are measured experimentally and these values are used in the theoretical investigations. In the control volume technique, partial differential heat conduction equation is reduced to ordinary differential equations in time, which are then solved sequentially. In Beck’s method, solution algorithm is developed under the function specification method to solve the inverse heat conduction equations. The IHTC was determined from the surface heat flux and the mold surface temperature by both the techniques and the results are compared.     

Free convection through conducting non-Newtonian fluids over a thin axisymmetric body

An analysis of free convection heat transfer in electrically conducting power law non-Newtonian fluid over a thin axisymmetric body of constant temperature is carried out. The uniform external magnetic field acts normally to the surface through the induced boundary layer. In view of the fact that most of the non-Newtonian fluids have large Prandtl number, the momentum equation is simplified. The equations of conservation of mass, momentum and energy which govern and describe the flow and heat transfer are solved numerically. The effect of the magnetic field on the velocity, temperature, the coefficient of friction and the Nusselt number are investigated. Numerical results are tabulated, presented graphically and discussed.     

Unsteady flow and heat transfer past an axisymmetric permeable shrinking sheet with radiation effect

In this paper, the problem of unsteady axisymmetric boundary layer flow and heat transfer induced by a permeable shrinking sheet in the presence of radiation effect is studied. The transformed boundary layer equations are solved numerically by an implicit finite‐difference scheme known as the Keller‐box method. The influence of radiation, unsteadiness and mass suction parameters on the reduced skin friction coefficient f′′(0) and the heat transfer coefficient ?θ′(0), as well as the velocity and temperature profiles are presented and discussed in detail. It is found that dual solutions exist and suction parameter delays the separation of boundary layer. Copyright © 2010 John Wiley & Sons, Ltd.     

A numerical investigation of primary surface rounded cross wavy ducts

 A three-dimensional numerical study was conducted to assess the hydraulic and heat transfer performance of primary surface type heat exchanger surfaces, called cross wavy (CW) ducts aimed for recuperators. The governing equations, i.e., the mass conservation equation, Navier–Stokes equations and the energy equation, are solved numerically by a finite volume method for boundary fitted coordinates. Periodic boundary conditions are imposed in the main flow direction. In this particular case laminar convective flow and heat transfer prevail. Details of the recuperator ducts and the numerical method as well as relevant results are presented. Received on 4 January 2001 / Published online: 29 November 2001     

Kinetic model of heat transfer processes in a fluidized bed

In the present paper equations are obtained for determining the temperature field in a fluidized layer. The heat and mass transfer processes in a fluidized bed depend significantly on the motion of the solid particles which form the bed. In any small volume of a fluidized bed with nonuniform thermal conditions there are particles with different average temperatures. Therefore it is natural to resort to the statistical representation of such a system, developed previously in [1, 2], for the study of the heat transfer processes. The expression obtained here for the heat conductivity coefficient of the bed is in good qualitative agreement with the experimental data.The author wishes to thank V. G. Levich for his interest and valuable discussions.     

Dispersed multiphase flow with air-driven runback of a liquid layer at a moving boundary

Multiphase flow with impinging droplets on an icing surface with a flowing supercooled surface layer is investigated. The air-assisted flowing layer is modelled with the cross-phase shear stresses imparted at the moving liquid/air interface. Runback and runoff of the surface layer are predicted by mass flow across the boundaries between adjacent elements in the numerical formulation. This liquid runoff is determined by coupled heat and momentum balances for the unfrozen water layer. The numerical analysis is developed with a control-volume-based finite element method (CVFEM). An Eulerian formulation with volume averaging is developed to accommodate the near-wall elements containing both dispersed and continuous phases. The predicted results are successfully validated through comparisons with analytical solutions and measured data.     

Controlling buffet through local mechanical action and heat addition near a lambda-shaped shock

Unsteady transonic flow past a two-dimensional airfoil with heat and momentum addition is numerically investigated. The flow analysis is based on the solution of the unsteady Reynolds equations closed by the k-ω SST turbulence model. The equations are integrated using the finite volume method. Several positions and shapes of the heat and momentum addition zones are considered for the purpose of determining an optimal means for controlling buffet. It is established that the most considerable variations in the buffet parameters are achieved when heat addition and mechanical action are realized on the upper wing surface. The thermal energy supply always increases the buffet frequency, while the mechanical action can both increase and reduce it.     

Solution of boundary layer flow and heat transfer of an electrically conducting micropolar fluid in a non-Darcian porous medium

In this paper, we have studied the effects of radiation on the boundary layer flow and heat transfer of an electrically conducting micropolar fluid over a continuously moving stretching surface embedded in a non-Darcian porous medium with a uniform magnetic field has been analyzed analytically. The governing fundamental equations are approximated by a system of nonlinear locally similar ordinary differential equations which are solved analytically by applying homotopy analysis method (HAM). The effects of Darcy number, heat generation parameter and inertia coefficient parameter are determined on the flow. Convergence of the obtained series solution is discussed. The homotopy analysis method provides us with a new way to obtain series solutions of such problems. This method contains the auxiliary parameter which provides us with a simple way to adjust and control the convergence region of series solution. By suitable choice of the auxiliary parameter, we can obtain reasonable solutions for large modulus.     

A two-dimensional splashing model for investigating impingement characteristics of supercooled large droplets

In this article, a two-dimensional (2D) splashing model is proposed to investigate the dynamics when supercooled large droplets (SLD) impinging on a wall surface in the aircraft-icing field. Energy conservation for droplet motion and impingement is used to capture the properties of the splashed droplets. A new statistical treatment of the droplet impinging energy and angle during the droplet-wall interaction is introduced in order to calculate the average dynamics of the SLD within a micro-control volume on wall surface. Based on the LEWICE predictions of droplet collection efficiencies and the available experimental ones, a new criterion for droplet splashing/deposition as well as a new formulation for the splashed mass is suggested. Lagrangian approach is adopted to describe the movement and impingement of SLD. The proposed model together with the previously developed droplet tracking method (DTM) for calculating droplet collection efficiency with the effect of droplet reimpingement constitute a relatively complete predicting approach of SLD impingement characteristics. Comparisons between the current predictions and the experimental observations, including SLD impingement over clean and contaminated airfoil surfaces as well as shapes of ice accretion in typical icing conditions, are carried out. Further, results obtained with the LEWICE splashing model are also plotted on the same graphs in order to assess the accuracy of the current splashing model in predicting SLD impingement. Results show that good agreement is achieved between the current predictions, including SLD impingement and ice accretion shapes, and the experimental ones. The predictions of the impingement distribution over contaminated surfaces obtained with the current splashing model show a much closer agreement with the experimental results than the ones obtained with LEWICE splashing model. For further investigation of SLD impingement, the properties of the droplet splashing and reimpingement during the ice accretion process are also addressed.     

Effect of thermo‐solutal stratification on recirculation flow patterns in a backward‐facing step channel flow

This paper presents results on the combined effect of thermo‐solutal buoyancy forces on the recirculatory flow behavior in a horizontal channel with backward‐facing step and the ensuing impact on heat and mass transfer phenomena. The governing equations for double diffusive mixed convection are represented in velocity–vorticity form of momentum equations, velocity Poisson equations, energy and concentration equations. Galerkin's finite‐element method has been employed to solve the governing equations. Recirculatory flow fields with heat and mass transfer are simulated for opposing and aiding thermo‐solutal buoyancy forces by assuming suitable boundary conditions for energy and concentration equations. The effect of Richardson number (0.1?Ri?10) and buoyancy ratio (?10?N?10) on the recirculation bubble and Nusselt and Sherwood numbers are studied in detail. For Richardson number greater than unity, distinct variations in the gradients of Nusselt number and Sherwood number with buoyancy ratio are observed for flow regimes with opposing and aiding buoyancy forces. Copyright © 2009 John Wiley & Sons, Ltd.