Swept Wing Research Articles

Sensitivity of flows over three-dimensional swept wings at low Reynolds number

High angle of attack flows over swept three-dimensional wings based on the NACA 0015 profile are studied numerically at low Reynolds numbers. Linear stability analysis is used to compute instability and receptivity of the flow via the respective three-dimensional (triglobal) direct and adjoint eigenmodes. The magnitude of the adjoint eigenvectors is used to identify regions of maximum flow receptivity to momentum forcing. It is found that such regions are located above the primary three-dimensional separation line, their spanwise position varying with wing sweep. The wavemaker region corresponding to the leading global eigenmode is computed and found to lie inside the laminar separation bubble (LSB) at the spanwise location of peak recirculation. Increasing the Reynolds number leads to the wavemaker becoming more compact in the spanwise direction, and concentrated in the top and bottom shear layers of the LSB. As sweep is introduced, the wavemaker moves towards the wing tip, following the spanwise displacement of maximum recirculation. Flow modifications resulting from application of different types of forcing are studied by direct numerical simulation initialised with insights gained from stability analysis. Periodic forcing at the regions of maximum receptivity to momentum forcing results in greater departure from the baseline case compared to same (low, linear) amplitude forcing applied elsewhere, underlining the potential of linear stability analysis to identify optimal regions for actuator positioning.

Plasma-based base flow modification on swept-wing boundary layers: dependence on flow parameters

This work examines the control of cross-flow instabilities (CFIs) and laminar–turbulent transition on a swept wing, through the plasma-based base flow modification (BFM) technique. The effect of experimentally derived plasma body forces on the steady boundary layer base flow is explored through numerical simulations. Linear stability theory is subsequently used to predict the net BFM effect on CFIs. Based on these preliminary predictions, experiments are conducted in a low-turbulence wind tunnel where a spanwise-invariant plasma actuator is installed near the wing leading edge and operated at constant input voltage and frequency. Various flow parameters governing the plasma-based BFM technique are investigated, namely the Reynolds number, angle of attack and wavelength of excited stationary CFI modes. Stationary and travelling CFIs are quantified by planar particle image velocimetry while the transition topology and location are recorded by infrared thermography. The results confirm the stabilising effect of BFM on the swept-wing boundary layer. However, the plasma-based BFM is found to render the boundary layer more susceptible to travelling CFIs. In the presence of both net BFM effect and intrinsic plasma unsteady perturbations, the plasma-based BFM technique achieves transition delay with specific combinations of Reynolds number, angle of attack and wavelength of excited stationary CFI modes. The present findings provide insights into the fundamental principles of operating plasma actuators within the context of BFM control.

АНАЛІЗ ВПЛИВУ КУТА СТРІЛОПОДІБНОСТІ НА ХАРАКТЕРИСТИКИ СТАТИЧНОЇ АЕРОПРУЖНОСТІ КРИЛА НАВЧАЛЬНО-ТРЕНУВАЛЬНОГО ЛІТАКА

The article is dedicated to the analysis of static aeroelastic phenomena in the wing of a trainer aircraft and the general influence of wing sweep on critical speeds. Aeroelasticity studies the interaction between aerodynamic and elastic forces and the impact of this interaction on the aircraft's structure. There are two types of aeroelasticity: static and dynamic. The article reviews and analyzes publications related to this issue. The analysis of the research shows that when designing an aircraft, it is necessary to determine the aeroelastic characteristics of its aerodynamic surfaces. In addition to structural rigidity, the sweep angle also influences these phenomena. The scientific novelty of the study lies in the fact that, unlike the analyzed research, the designed wing structure of a trainer aircraft is considered.For the designed trainer aircraft, the critical speeds of static aeroelastic phenomena were determined. The calculation of the influence of the wing sweep angle on the divergence speeds of the trainer aircraft confirms previously obtained results regarding this influence. When the sweep is increased or decreased for wings with the same span, the torsional rigidity of the wing decreases due to the increase in the length of the line of centers of rigidity. However, the most significant impact on the aeroelastic characteristics is made by the shift of the bending center when the sweep angle changes.With forward sweep, the arm of the aerodynamic forces relative to the bending center increases, thereby reducing the critical divergence speeds. For straight sweep, as the sweep increases, the bending center moves forward, and at a certain point, the bending center is located in front of the aerodynamic center. This means that the wing twist angle decreases as the aerodynamic force increases, which, in the considered planar case, eliminates the occurrence of divergence.The preliminarily determined maximum flight speed does not exceed the critical speeds of divergence and aileron reversal. This result represents the practical value of the study. To achieve higher flight speeds, wings with straight sweep should be used. More precise computational or experimental methods for determining aeroelastic characteristics should be used in the subsequent design stages to prevent aeroelastic phenomena across the entire operational range of speeds and altitudes of the trainer aircraft.

Design perspectives of a system identification approach of the NATO AVT-316 multi-swept wing aerodynamics

Air supremacy requires enhanced maneuvering capabilities at high angles of attack and even beyond the stall angle. High angle-of-attack (high-α) maneuverability can be achieved using swept wings and tailored leading-edge shapes to replace local and disorganized flow separation with controlled separation-induced leading-edge vortices (LEV). Strong leading-edge vortices delay the stall to higher angles and generate a non-linear lift increment up to the angle of attack where the jet-type flow structure of LEV changes to a reversed-flow bubble (vortex breakdown phenomenon). A multi-swept wing configuration has the potential to delay the vortex breakdown to even higher angles as the vortices formed over the front wing can energize vortices formed over the main wing and lead to vortex merging. This article is focused on understanding the unsteady interactions between multiple leading edge vortices formed over multi-swept wing configurations in the subsonic speed regime. Specifically, this article investigates the aerodynamic characteristics of wings being evaluated under the initiative of the NATO STO AVT-316 Task Group. Highly refined meshes, and the use of hybrid turbulence models such as Detached Eddy Simulations (DES) and Delayed DES are required to accurately resolve the vortical flows over these wings. Simulating all flow conditions of interest is a computationally expensive approach. A system identification method was therefore proposed to rapidly and accurately generate aerodynamic models of these wings. The method uses a novel piece-wise chirp (constant amplitude and increasing frequency) motion as an input signal (training maneuver). The motion computational cost is equivalent to the cost of six static CFD simulations, however it can predict aerodynamic responses of the wings over a wide range of angles of attack. The method has been tested for double and triple delta wings. Some design considerations are provided based on predicted flow features and aerodynamic data. Prediction results with different turbulence models, sting geometries, grid resolutions and an adaptive mesh refinement approach are provided.

Impact of Blowing Location on the Aerodynamic Characteristics Over the Delta Wing

The performance of an aircraft can be enhanced by altering the flow field favourably by adopting flow control techniques. The present study deals with the application of the active flow control methods on a sharp-edged delta wing with a wing sweep of 65°. The concept of blowing was employed as an active flow control technique. The blowing technique is applied on the suction surface of the delta wing by varying its location. The various identified locations of the blowing holes are 1.62 %, 3.24 % and 4.86 % of root chord from the leading edge to the centre of the blowing holes. The computation is performed using the commercial software ANSYS Fluent. An unsteady, incompressible Reynolds-averaged Navier–Stokes equation and the shear-stress transport k-ω turbulence model are employed. The angles of attack varied in the range of 0°<α<35° and Reynolds number is 2.64×106 and the jet momentum coefficient is fixed at 0.05. The blowing of air from the injection region enhances the strength of the leading-edge vortices, resulting in a delay in the vortex breakdown. The performance of the delta wing is greatly improved while using the blowing method specifically for the blowing holes located at 3.24 % of root chord from the leading edge compared to without the blowing method.

Modeling and Analysis of Kamikaze UAV Design with 3 Different Wing Configurations

Appropriate design parameters need to be determined for unmanned aerial vehicles that can perform kamikaze missions. In this study, a UAV with 3 different wing configurations and a fuselage and tail wings were designed, and the flow around the wing was examined using computational fluid mechanics. Advanced modeling techniques were employed to simulate and analyze the aerodynamic behavior of these configurations. The effect of angle of attack (AoA), wing positioning on the fuselage, and wing configurations were investigated. Due to the effect of the wing sweep angle, high-pressure values in the arrow-angle wing were lower than in rectangular and trapezoidal wings. In a similar situation, the flow separation on the arrow-angle wing is less advanced towards the wing tip. When the wing type and connection location were examined, the highest ${{C}_{l}}/{{C}_{d}}$ ratio was obtained in the trapezoidal model connected to the fuselage in the middle. The results of numerical wing models compared with the theoretical lift coefficient were consistent. Trapezoidal and rectangular wings had a high lift coefficient, but after ${{15}^{\circ }}$ of AoA, the lift coefficient decreased. At angles of attack beyond ${{15}^{\circ }}$, the arrow-angle wing still has an increasing lift coefficient. As the angle of attack increased, the drag coefficient was also enhanced. The lowest drag coefficient occurred in the arrow-angle wing model. Up to ${{5}^{\circ }}$of AoA, all wing models raised the ${{C}_{l}}/{{C}_{d}}$ ratio. The ${{C}_{l}}/{{C}_{d}}$ ratio decreased at higher angles of attack. As a result of the examination, it would be more correct to choose trapezoidal and arrow-angle wings.

Steady as they hover: kinematics of kestrel wing and tail morphing during hovering flights.

Wind-hovering birds exhibit remarkable steadiness in flight, achieved through the morphing of their wings and tail. We analysed the kinematics of two nankeen kestrels (Falco cenchroides) engaged in steady wind-hovering flights in a smooth flow wind tunnel. Motion-tracking cameras were used to capture the movements of the birds as they maintained their position. The motion of the birds' head and body, and the morphing motions of their wings and tail were tracked and analysed using correlation methods. The results revealed that wing sweep, representing the flexion/extension movement of the wing, played a significant role in wing motion. Additionally, correlations between different independent degrees of freedom (DoF), including wing and tail coupling, were observed. These kinematic couplings indicate balancing of forces and moments necessary for steady wind hovering. Variation in flight behaviour between the two birds highlighted the redundancy of DoF and the versatility of wing morphing in achieving control. This study provides insights into fixed-wing craft flight control from the avian world and may inspire novel flight control strategies for future fixed-wing aircraft.

Unsteady Lifting-Line Theory for Camber Morphing Wings State-Space Aeroelastic Modeling

A novel frequency-domain analytical-numerical model for the aerodynamic solution of camber morphing wings is developed. The model combines an unsteady lifting-line formulation with Küssner–Schwarz aerodynamic theory to provide the pressure distribution and thus the transcendental aerodynamic matrix that relates the generalized aerodynamic forces to the Lagrangian coordinates of a given wing structural dynamics model. The state-space form of the aerodynamic loads is obtained from the rational matrix approximation of the aerodynamic matrix. This, combined with the wing structural dynamics operator, can readily provide the state-space form of the aeroelastic problem. To assess the accuracy of the proposed aerodynamic solution method, numerical investigations are performed. These consist of its application to several conventional wing configurations and wings with camber morphing, followed by the comparisons of the corresponding solutions with the predictions given by a well-validated panel-method solver for potential flows. These results are shown to be in very good agreement, thus demonstrating that the proposed approach can capture the effects due to wing tapering, sweep angle, and camber morphing, while requiring a remarkably lower computational effort. The excellent accuracy of a few-pole finite-state approximation of the aerodynamic matrix is finally proven for a swept wing with camber morphing.

The optimization method of wing plasma ice shape regulation based on quantitative assessment of flight risk

Plasma ice shape regulation is a technology which uses plasma actuator to regulate the continuous ice into safer intermittent ice by its significant thermal effect with limited energy. Whether plasma ice shape regulation could reduce flight risk is a new problem under the wing with continuous ice. The 3D printed ice shapes were arranged on the leading edge of the wing based on NACA0012 airfoil, aiming to simulate the configuration after ice shape regulation. And the aerodynamic parameters were obtained by wind tunnel experiments. The experimental results showed that the ratio of signal regulation ice width d\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$d$$\\end{document} to chord length of the wing bA\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$b_{A}$$\\end{document} determined the aerodynamic characteristics, and the aerodynamic characteristics changed better compared with configuration of the continuous ice. However, the flight risk of the wing under given regulation ratio is unknown. Based on the straight and swept wing after regulating, the flight safety boundaries were simulated by the reachable set method. Further, a method of quantitative assessment of flight risk is proposed. Quantitative values of risk were calculated. The results show that the flight risk all decreases from level 2 to level 4 compared with configuration of the continuous ice when d/bA\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$d/b_{A}$$\\end{document} equals 0.15 under conditions of swept and straight wing.

Ice Shape Convergence in Multistep Ice Accretion Simulations

Local errors in the geometrical description of the ice front are amplified in multistep simulations over conformal meshes due to the coupling of aerodynamics, water impingement, ice accretion, and grid deformation. Small perturbations in the initial phase of ice formation possibly result in dramatically different ice shapes, which can hinder the stability of the multistep procedure. This problem is analyzed by investigating the combined effects of space and time discretization on the ice growth over airfoils and three-dimensional wings. We propose an automatic procedure for selecting the time interval for the update of the aerodynamics and particle impingement. A local growth limiter λ is introduced here to bound the local ice thickness growth to be comparable to the local grid spacing on the surface, resulting in an automatic adaptive time step to be used in the multistep simulation. Examples are provided for three-dimensional cases under both rime and glaze conditions over straight and swept wings. These examples highlight the different accretion mechanisms of the two ice regimes and preliminarily indicate that ice-shape convergence can be achieved for low values of λ.

Sensitivity Analysis for Design Parameters of Electric Tilt-Rotor Aircraft

In recent years, there has been rapid development in electric aircraft, particularly electric vertical takeoff and landing (eVTOL) aircraft, as part of efforts to promote green aviation. During the conceptual design stage, it is crucial to select appropriate values for key parameters and conduct sensitivity analysis on these parameters. This study focuses on an electric tilt-rotor aircraft and proposes a performance analysis method for electric aircraft while developing a general design tool specifically for this type of aircraft. Subsequently, the impact of wing incidence angle, sweep angle, span, propeller solidity, battery-specific energy, and battery mass on range, maximum takeoff weight, and hover power are analyzed. The results show that the battery mass, wingspan, and wingtip chord length have great effects on the maximum takeoff weight; among these, battery mass had the greatest influence. In terms of range, the battery energy density has a great positive effect on range, while the increase in wing angle of incidence, wingtip chord length and battery mass have some negative effects on range.

Roll Maneuverability of Transonic High-Aspect-Ratio-Wing Aircraft with Flared Folding Wingtips

This paper investigates the roll maneuvers of a very flexible, transonic high-aspect-ratio-wing aircraft with flared folding wingtips, representing a potential future commercial transport configuration. The computational study leverages a fully coupled nonlinear aeroelastic-flight dynamics framework for simulating roll maneuvers commanded by deflecting control surfaces inboard or outboard of the flared folding wingtip hinge. Releasing the flared folding wingtips during flight makes the aircraft roll faster only when deflecting outboard-of-hinge leading-edge control surfaces. However, the achieved roll angles and rates are much lower than the ones obtained by deflecting inboard-of-hinge trailing-edge control surfaces, whose roll control effectiveness is not impacted or even degrades when releasing the flared folding wingtips. These trends are explained by analyzing the role of sweep-induced washout effects in the very flexible, transonic high-aspect-ratio wing, revealing physical mechanisms not identified in previous studies on straight or swept wings with moderate flexibility. Higher wing out-of-plane bending stiffness increases the effectiveness of flared folding wingtips in enhancing roll maneuvers at low dynamic pressure but plays a slight role at high dynamic pressure. This study highlights the need to consider the impact of washout effects for different flight conditions and levels of wing flexibility when assessing the ability of flared folding wingtips to enhance the roll maneuverability of transonic high-aspect-ratio-wing aircraft.

Impact of Installation on the Performance of an Aero-Engine Exhaust at Wind-Milling Flow Conditions

Abstract This paper presents a numerical investigation of the effect of wing integration on the aerodynamic behavior of a typical large civil aero-engine exhaust system at wind-milling flow conditions. The work is based on the dual stream jet propulsion (DSJP) test rig, as will be tested within the transonic wind tunnel (TWT) located at the aircraft research association (ARA) in the UK. The DSJP rig was designed to measure the impact of the installed pressure field due to the effect of the wing on the aerodynamic performance of separate-jet exhausts. The rig is equipped with the dual separate flow reference nozzle (DSFRN), installed under a swept wing. Computational fluid dynamic simulations were carried out for representative ranges of fan and core nozzle pressure ratios (CNPR) for “engine-out” wind-milling scenarios at end of runway (EOR) takeoff, diversion, and cruise conditions. Analyses were done for both isolated and installed configurations to quantify the impact of the installed pressure field on the fan and core nozzle discharge coefficients. The impact of fan and core nozzle pressure ratios, as well as freestream Mach number and high-lift surfaces on the installed suppression effect, was also evaluated. It is shown that the installed pressure field can reduce the fan nozzle discharge coefficient by up to 16%, relative to the isolated configuration for EOR wind-milling conditions. The results were used to inform the design and setup of the experimental activity which is planned for 2023.

An experimental investigation on the aerodynamic characteristics and vortex dynamics of a flying wing

AbstractIn this paper, we present a detailed experimental investigation mainly on the vortical flow fields and the associated vortex breakdown phenomena over a non-slender flying wing (sweep angle, ${\rm{\Lambda }}$ = 53°). In the process, the aerodynamic coefficients were also determined using a six-component force balance. Surface oil flow visualisation, surface pressure measurements and particle image velocimetry (PIV) measurements, in various crossflow planes and in a longitudinal plane passing through the leading-edge vortex core, were carried out at various Reynolds numbers to understand the flow field over the non-slender flying wing. Aerodynamic characteristics of the flying wing show local peaks and valleys in the pitching moment coefficient. The surface flow visualisation reveals that the nonlinearity of the pitching moment curve is due to the complex nature of vortical flow structures. The flow visualisation also demonstrates the presence of a wave-like surface pattern, and its size is found to reduce with increasing Reynolds numbers. The present PIV measurements confirm that this wave-like surface pattern is associated with vortex breakdown phenomena. These measurements also reveal that the vortex breakdown has not reached the apex of the wing, even at post-stall angle-of-attack. For pre-stall ( $\alpha $ = 20°) flow regimes, it is observed that the location of the vortex breakdown moves downstream as the Reynolds number increases, but this influence is minimised at near-stall ( $\alpha $ = 25°) and post-stall ( $\alpha $ = 30°) flow regimes. Reconstructed velocity field using the first 10 dominant proper orthogonal decomposition (POD) modes reveals that the nature of the vortex breakdown over the flying wing is a spiral-type vortex breakdown.

Wind tunnel experiments of bending-torsion and body-freedom flutter on flying wing unmanned aerial vehicles

Wind tunnel experiments are performed on a flying wing model to investigate the effects of key design parameters and their cross-interactions on the flutter speed and frequency for bending-torsion flutter (BTF) and body-freedom flutter (BFF). These design parameters include wing sweep angle, mass of tip weights, and location of weights along the wingspan. A slider and rail setup for BFF is designed to enable rigid-body pitching and plunging degrees-of-freedom, while a custom fixed clamp is used to perform the BTF experiments. BTF and BFF experiences different dominant modes while undergoing flutter, with the former's dominant mode being torsion and the latter being the rigid-body short-period mode. For BTF, flutter speeds increase as weights are moved towards the wingtip or increasing tip weights, which increases the inertia of the system to enhance stability. The effect is opposite for BFF as the movement of weights outboard or increase of tip weights result in a rearwards shift of the centre of gravity that destabilises the system. In general, a higher sweep angle leads to increased sensitivity of flutter speeds to changes in mass of tip weights or changes to spanwise location of the weights. The findings from this study provide insights to the design of flying wing unmanned aerial vehicles for increased stability margins.

Laminar post-stall wakes of tapered swept wings

While tapered swept wings are widely used, the influence of taper on their post-stall wake characteristics remains largely unexplored. To address this issue, we conduct an extensive study using direct numerical simulations to characterize the wing taper and sweep effects on laminar separated wakes. We analyse flows behind NACA 0015 cross-sectional profile wings at post-stall angles of attack $\alpha =14^\circ$ – $22^\circ$ with taper ratios $\lambda =0.27$ – $1$ , leading-edge sweep angles $0^\circ$ – $50^\circ$ and semi aspect ratios $sAR =1$ and $2$ at a mean-chord-based Reynolds number of $600$ . Tapered wings have smaller tip chord length, which generates a weaker tip vortex, and attenuates inboard downwash. This results in the development of unsteadiness over a large portion of the wingspan at high angles of attack. For tapered wings with backward-swept leading edges, unsteadiness emerges near the wing tip. On the other hand, wings with forward-swept trailing edges are shown to concentrate wake-shedding structures near the wing root. For highly swept untapered wings, the wake is steady, while unsteady shedding vortices appear near the tip for tapered wings with high leading-edge sweep angles. For such wings, larger wake oscillations emerge near the root as the taper ratio decreases. While the combination of taper and sweep increases flow unsteadiness, we find that tapered swept wings have more enhanced aerodynamic performance than untapered and unswept wings, exhibiting higher time-averaged lift and lift-to-drag ratio. The current findings shed light on the fundamental aspects of flow separation over tapered wings in the absence of turbulent flow effects.

The performance and application of sweep wing aircraft

With the development of science and technology, more and more civil and military aircraft have adopted swept-back wing shapes. this paper explores the relationship between the angle of swept-wing and the lift, drag and lift-drag ratio of the aircraft, and expounds it in combination with the actual swept-wing aircraft, such as Boeing737 and MIG-23 to find the best sweep angle of the aircraft in the appropriate range and analyze the effect of swept wing angle on flight speed. This paper studies and analyzes the data through literature retrieval and data processing and cites data models from many academic journals. The experimental data are mainly derived from computational fluid dynamics and wind tunnel simulation. This paper finds that a wider sweep angle can result in better aerodynamic performance, suited for supersonic flight, while a smaller sweep angle can result in a better lift-drag ratio, it is suitable for takeoff and landing of aircraft and low-speed flight attitude.

Workflow for predictor–corrector simulations of in-flight ice accretion, with applications on swept wings

Workflow for predictor–corrector simulations of in-flight ice accretion, with applications on swept wings

Noise reduction for high-lift devices on a swept wing model by droop nose

The noise of slat in the high-lift device of large aircraft is one of the main noise sources of large aircraft during take-off and landing. This paper investigates a noise reduction measure which consists in droop nose for high-lift device of three-dimensional sweep wing instead of leading-edge slat. The unsteady flow field and the far-field noise are predicted by a hybrid SA-IDDES method coupled with the integration of the Ffowcs Williams–Hawkings equation (FW-H). Based on the optimal droop nose configuration from 2D airfoil, the results of the 3D swept wing model with the droop nose on the leading edge of the inner wing are given. Numerical results show that the droop nose configuration can reduce the overall sound pressure level (OASPL) of the intermediate frequency noise, and the peak noise can be reduced by 15 dB, and the sound pressure level of the broadband noise at low and high frequencies is also partially suppressed. The maximum lift-to-drag ratio of the droop nose configuration increased by 2.7 %. In addition, the noise reduction mechanism of leading-edge droop is analyzed by combining the vortex sound theory and the generation and evolution process of the vortex structure. The leading-edge droop nose configuration fundamentally not only eliminates the vortex structure of the slat cove of the leading-edge slat configuration, but also avoids the interaction between the slat slot flow and the upper surface of the main wing. These directly lead to a thinner boundary layer vortex near the top of the main wing, and the wall pressure pulsation is smaller than that of the leading-edge slat high-lift device, thereby reducing noise. The reduction of the BEF distribution on the upper surface boundary shows again that the droop nose inhibits the vortex on the upper surface of the main wing. The change of the Lamb vector divergence reveals the noise reduction mechanism of the drooping leading edge configuration directly from the sound source.

Roughness implementation and convective heat transfer coefficient computation toward ice accretion simulation

In the simulation of icing, the calculation of the convective heat transfer coefficient for rough surfaces has significant importance. Our study investigates the extension of the γ-Reθt transition model for rough icing surfaces and aims to provide an alternative to the integral boundary layer method for facilitating the calculation of the heat transfer coefficient in three-dimensional icing simulations. Initially, the performance of the γ-Reθt transition model in predicting the transition characteristics and heat transfer for smooth surfaces is examined. Subsequently, to compensate for the deficiency of the γ-Reθt transition model in predicting the premature transition caused by roughness, an additional model equation Ar is introduced to simulate the influence of the surface roughness. Three test cases, including a clean cylinder, an iced cylinder, and a swept wing, are used to validate the applicability of the γ-Reθt-Ar transition model for heat transfer computations on rough surfaces. The results demonstrate that the γ-Reθt-Ar model is superior to other turbulence models in computing convective heat transfer for two-dimensional and three-dimensional icing rough surfaces and shows its potential for ice accretion simulation.