Selection and Impact of an Aerofoil Leading Edge on Boundary Layer Transition

The choice of leading-edge aspect ratio (AR) plays a crucial role when plan-ning boundary layer wind tunnel tests on a flat plate. Poor selection of the leading-edge profile hampers effectiveness of the experiment and increases testing costs associated with interchanging of leading edges to attain accurate results. Thus, the appropriate selection of the leading edge is a very crucial part of the wind tunnel experiment process. It is argued that the curvature of the leading edge and thus the AR is of paramount importance to achieve accurate results from the wind tunnel testing. In this project, seven different elliptical leading edges were tested, and their performance was compared with an ideal leading edge with zero thickness. Experiments and computation have been done for leading edges ranging from AR6 to AR20. Results were evaluated for boundary layer transition onset location, and it was found that AR20 has the least influence on the flow structure when compared to the ideal leading edge. A study of the flow structure at the stagnation point indicates an increase in adverse pressure gradient with an increase in the AR but also shows a decrease in the size of the stagnation region. The presence of a higher AR leading edge reduces the turbulent spot production rate, which is one of the primary causes of boundary layer transition. This paper presents a correlation that enables aerodynamicists to quantify the impact of the leading-edge AR on transition. A typical case is also presented to compare the relative performance of a wedge and the higher AR leading edge, which provides a choice between an elliptical or a wedge-shaped leading edge.


Introduction
Boundary layer transition is influenced by numerous factors including changes in the surface curvature, turbulence intensity, surface roughness and the velocity of the flow. In a controlled environment wherein the turbulence intensity, velocity of the flow and the flat plate are unchanged, a change in the aspect ratio (AR) of the leading edge will be the primary cause for changes in the transition onset location. The primary factor in this case is the change in the surface curvature.
[1] through his series of experiments found that the presence of a convex curved surface tends to stabilize the boundary layer while a concave surface tends to destabilize the boundary layer. According to the researcher, a laminar boundary layer on a concave surface becomes unstable because of the centrifugal force and three-dimensional disturbances whereas boundary layers on convex surfaces are free from centrifugal forces. Consequently, the effects of concave surface geometry have an impact on the transition onset location by up to 56% [2]. Thus, the curvature whether convex or concave of the leading edge can influence the location of boundary layer transition. In the case of convex surfaces (as is the case with the elliptical leading edges in the present case), the impact of surface curvature has been well documented in experimental studies, such as the one conducted by [3]. According to the researcher, a change in the surface curvature results in a change in the turbulent spot propagation rate nσ . Larger adverse pressure gradient would lead to greater spot propagation rates. These spot propagation rates have a direct impact on boundary layer transition. Changes in the leading-edge AR also led to changes in the pressure distribution along the leading edge [4] [5]. Likewise, [6] also believes that turbulent spot propagation is directly impacted by changes in the AR due to changes in the surface curvature. Therefore, a higher adverse pressure gradient along the leading edge would result in greater spot propagation which in turn would impact the transition onset location.
Studies of boundary layers over flat plates found that the leading edge was a particularly sensitive region where the disturbance first enters the boundary layer [7] [8]. The turbulent spot propagation rate can thus be thought of as an indicator of the magnitude of disturbance that enters the boundary layer due to the presence of the leading-edge AR. This can further be corroborated using the conclusions proposed by [9] [10], which indicate that the receptivity of the boundary layer to freestream disturbances greatly depends on the leading-edge curvature, discontinuities in the surface curvature and surface inhomogeneities.
Experiments to understand the effect of leading edge bluntness on boundary layer development have shown that an increase in the AR of the leading-edge leads to a decrease in the amplitude of Tollmien-Schlichting (TS) waves, which are the 2) The smoothening of the discontinuity at the ellipse/flat plate results in a reduction in the amplitude of TS waves by about 50%.
3) The disturbance which is responsible for the generation of the TS waves originates from the location of the maximum adverse pressure point and depends on the magnitude of the adverse pressure.
The consequence of an increase in the AR is well documented in the previous paragraphs and indicates that an increase in the AR would lead to a delay in the transition. However, according to [11], a change in the AR does not affect the transition onset location in a predictable manner and the transition onset location in these situations varies greatly with each change in AR.
A lack of database clearly highlights the appropriate choice of leading edge for computational, experimental or numerical research. While all previously documented literature has discussed about the impact of surface curvature on boundary layer transition, very little information about integrated viewpoints is available of the impact of the different aspect ratio leading edge and their influence on the boundary layer and flow transition under uniform flow conditions. This research emphasises on ascertaining and quantifying the impact of the change in curvature on boundary layer transition. Creating a database of different AR leading edges and their impact on boundary layer transition is the other important aim of this paper.
The ERCOFTAC T3A, T3B, and T3C experiments are used as benchmark cases in a few experimental and computational research cases with zero pressure gradients [12]. The T3 test cases make use of a semi-circular leading edge. [13] found out that while using the T3 test cases for simulations, a little or no information was available on the impact of changes in the gradient of elliptical leading edges on boundary layer and transition measurements. Thus, there is a dearth of data available for the optimal leading edge to be used for simulations and experiments. The lack of data is extremely crucial, because each small change in leading edge geometry can influence transition. It is important to note that for an elliptical leading edge, there might be no separation bubble as in the case of a blunt leading edge. Changes of leading-edge aspect ratio will lead the increase of elliptical curvature at the tip and curvature decrease at other places. The impact of changes in the curvature on the stagnation point, on rest part of the elliptic curve and connection of the elliptic and flat plate as well as the correlation between the curvature and the pressure gradient have not been studied in isolation. The changes in the flow structure based on the change in the geometry need to be investigated in detail. The presence of a sharp leading edge theoretically should have the best aerodynamic performance. Therefore, an analysis needs to be done to find out the impact of increasing AR on transition onset and its comparison to a theoretical leading edge with AR∞. Research has also been conducted to identify a "modified super ellipse" that could be used for a few computational and practical test cases. This super ellipse could possibly eliminate the leading edge/flat plate juncture discontinuity, by smoothing out the juncture discontinuity by using a polynomial function [7]. One of the key results from the experiment and computation was that the presence of a modified super ellipse can lead to minimizing the region of non-zero pressure gradients which reduces the receptivity of the boundary layer to disturbances. However, while this super ellipse has been tested in a couple of test-cases its implementation in several different test scenarios could be cumbersome when a number of different aspect ratio leading edges are to be tested. Thus, in such cases, it is necessary to find an optimum leading edge that would serve the purpose for both experimental as well as computational test scenarios.
Many experiments and simulations have been carried out by varying the AR of the leading edge and analysing the effects of these changes on boundary layer transition. All experiments and simulations have been conducted under a uniform set of test conditions and several different test-cases have been simulated keeping all parameters the same. This study addresses the significance of an elliptical leading edge on boundary layer transition and draws results related to appropriate leading edge that could possibly be used by researchers depending on the accuracy sought in individual test cases by means of providing the comparison of different practical aspect ratios with a theoretical leading edge with AR∞. The notion that an elliptical leading edge is the most appropriate geometric configuration has also been challenged with the use of a wedge-shaped leading edge in this study.

Theoretical Concept
To study the impact of changes in AR of elliptical leading edges, the authors have focused on the impact of changes in geometry of the ellipse and identified four key areas mainly the stagnation region, the change of curvature along the length of the leading edge, the pressure gradient, and its changes along the length of the leading edge, and the turbulent spot propagation rate due to a change in the aspect ratio.
To compare different leading edges, an ideal leading edge with zero thickness and infinite length has been considered i.e. AR∞. Comparison to AR∞ would indicate the receptivity of the boundary layer to disturbances due to the presence of a change in geometry. In addition, comparison to AR∞ will provide the required information to select a practical and optimum leading edge for future wind tunnel tests as well as simulations. curvature discontinuity/sharp changes in curvature exist (e.g. the joint of the elliptical leading edge to the flat plate section) tend to provide regions of high receptivity [7]. The rate of change of curvature along the length of the leading edge can give insights into the flow of the fluid over leading edge. At the start of the leading edge, the curvature for each different AR leading edge is different. A study of the change in the curvature of each individual leading edge is a possible indicator as to the flow development on the leading edge. The pressure gradient will also change with a change in curvature. When the surface is convex as is the case for all elliptical leading edges in the present research, the turbulent spot propagation rates would be directly correlated to the curvature. The spot propagation rate can be used as an indicator of the impact of the leading edge on boundary layer transition. The curvature change can also be linked to a change in the pressure gradient. The aerodynamic performance of test specimen depends on the impact of the leading edge and a change in the curvature would alter the pressure gradient. The curvature changes at different points on each leading edge have been studied numerically and theoretically. A comparison of the curvature change at specific points on each leading edge will highlight the impact of changes in the aspect ratio on the fluid flow. Thus, it can be said that the curvature will have a direct impact on the stagnation region, the pressure gradient as well as the turbulent spot propagation rate.
As all leading edges are elliptical in nature, the fundamental equations for calculating the radius R, curvature K as well as gradient for each individual leading edge are derived from the fundamental equation of the ellipse [14]: where, x, y: X and Y co-ordinates for the ellipse respectively; a: Major Radius for the ellipse; b: Minor radius for the ellipse. The curvature κ and the radius R are calculated using the equations [15]: : First and second order differential equations for Equation (1). And,

Stagnation Region
The stagnation region is the region of low-velocity fluid just fore of the nose of could be the first region where the disturbance might start to enter the boundary layer. Figure 1 shows the individual profile of each leading edge and indicates that the radius of the curve at the start of the leading edge reduces with an increase in AR.

Pressure Gradient
One effect of the change in the curvature is the alteration of the pressure gradient. In addition, destabilization of the boundary layer takes place due to presence of adverse pressure gradients and other factors. The presence of adverse pressure gradients increases the rate of the growth of the disturbance and enhances the receptivity of the boundary layer to disturbances [16]. The presence of a pressure gradient leads to an acceleration or deceleration of flow around the flat plate. This acceleration or deceleration (change in pressure gradient) affects the transition to a certain extent. Negative acceleration takes place at adverse pressure gradients [3]. Research conducted by [17] further exemplify that the increase in the adverse pressure gradient results in a much larger spot propagation rate and hence a quicker transition and vice-versa. The pressure gradients, in this case, are directly linked to changes in the curvature of different elliptical leading edges. The pressure gradients have been studied at several key points along different leading edges and a direct link can be established between the curvature change and the pressure gradient. The critical point at which the favourable pressure gradient on the leading-edge changes to adverse pressure gradient will have a direct impact on the transition onset as all other factors apart from the curvature of the leading edges are unchanged.

Turbulent Spot Propagation Rate (nσ )
The changes in the curvature enhance the receptivity of the boundary layer to disturbances. The spot propagation rate is an indicator of the disturbances in the boundary layer. The turbulent spot propagation rates increase with an increase in adverse pressure gradient as discussed in the preceding section. With all other parameters being kept constant, it can be said that the spot propagation depends solely on the curvature and hence the aspect ratio of the leading edge. The turbulent spot propagation rate can be calculated using the equation [2]:

Methodology
To investigate the proposed issues, wind tunnel experiments were conducted first to see if there are significant effects of curvature change on transition. Then, analytic formulae and methods were applied to study characteristics of the curvature of the ellipse. Then, CFD was used to compute the effects of changes of curvature on flow transition. Finally, efforts were devoted to identify the trends in the data and correlate relevant parameters to provide insight of proposed investigation.

Experiment Setup
The present set of experimental test cases were performed in a LTWT (Low Turbulence Wind Tunnel) which was characterised by an extremely low turbulence intensity of 0.018% in the velocity interval of 5 m•s −1 -75 m•s −1 . The experimental set-up is very similar to the ones conducted by previous researchers such as Westin et al. [18] and Klingmann et al. [6] and uses a flat plate which is undisturbed and has zero pressure gradient. The details of the wind tunnel and its construction are given by Bhatia et al. [19]. The experiments were carried out on a plate that was 2.  An X-axis moving system, a Z axis and a Y-axis moving system along with control system used along with a DC motor as a servo system to drive each axis along with a ball screw pair to transmit the motion. The effective stroke of the X axis is 1500 mm with a minimum moving distance of 1.0 mm and an accuracy of  The test cases (sample models) involved in the experimental work were AR1 (1:1), AR2 (2:1), AR4 (4:1), AR6 (6:1), AR9 (9:1) and AR12 (12:1) although the current paper highlights the results for AR6, AR9, and AR12. The schematics of the test setup including constant temperature anemometry (CTA) and the linear traverse along with the some sample leading edges are shown in Figure 2 and found that the flow in the wind tunnel was affected by the presence of the linear traverse arm. The blockage effect was found to affect the transition onset location by about 20% [21]. The impact of the surface roughness was also quantified, and it was found that the surface roughness of the plate can impact the boundary layer transition by up to 9% [21]. Thus, it can be assumed that the impact of the blockage and the flat plate surface roughness can impact the boundary layer transition by up to 29%. This is mentioned in the experimental results section of this paper (Section 3.1).

Simulation Setup
The flat plate geometry used for the simulation is a 2D 1.8 m long and 0.012 m thick flat plate with AR12 as the default leading edge for all reliability calculations. This plate is identical to the plate used by [13] in their series of simulations. Consequently, the flat plate is modified for changes in the AR. All other parameters are kept constant. The simulation setup and the test cases in this paper are similar to the experiments conducted by [19].  The mesh used for the simulations is a structured multi-block mesh with the mesh being clustered around the leading edge as shown in Figure 4. A maximum y+ value of 0.1 is maintained throughout the grid and the wall normal expansion ratio is set to 1.07. A mesh with 400 k elements is chosen for the simulation based on grid independence studies and as optimum balance between accuracy and computational processing requirements. A no-slip wall boundary condition was used for the flat plate along with a velocity inlet and a pressure outlet.
All other planes of the domain are given a symmetry boundary condition.
This study uses the Transition SST model in FLUENT for a flat plate with different geometric configurations. The Transition SST model has been proven to be suitable for accurate prediction of transition onset for high curvature, non-flat and wavy geometries making it an appropriate choice for the leading-edge simulations [22]. All model constants for the simulations were kept exactly at their default values. It was found from preliminary simulations and previous work conducted by [13] that natural transition was not sensitive to the viscosity ratio between a ranges of 5 -10. It was decided to take the higher value of 10

Results
The results in this section will be divided into two main subsections viz. the experimental results which show a comparison between three leading edges AR6, AR9 and AR12 and computational results which will be used to provide an insight into the flow structure around different AR leading edges.

Experiment Results
Displacement thickness for the experimental case at AR12 and its comparison to the theoretically obtained Blasius displacement thickness is shown in Figure 6.
The experimental displacement thickness closely matches the Blasius value. The displacement thickness at the end of the flat plate is roughly 1.49 mm for the experimentally obtained value while it is roughly 1.5 mm for the Blasius value.
However, while there is a close correlation, transition onset cannot be ascertained.
The shape factor H is used for ascertaining the transition onset in the boundary layer.
The shape factor H gives an exact indicator of transition onset location and length of the transition region. Theoretical calculations for the shape factor show that for transition onset in the boundary layer, the value for H should be approximately 2.6. The point for the turbulent reattachment of the boundary layer is indicated by a shape factor value of H = 1.4 [26]. As the shape factor reduces, the flow tends to indicate a highly turbulent flow structure [27]. For AR12 as shown in (c), the length of the transition region is the largest indicating that the amplification of disturbances is the lowest. This can also be demonstrated by the turbulent spot propagation rates (Table 1).   can be observed that there is a significant difference between AR6 and AR12.
The turbulent spot propagation rates indicate that AR6 would have an extremely receptive boundary layer and amplification of disturbances happen at a much faster rate. Figure 8 shows the percentage difference between the transition onset for AR6, 9 and 12. Transition onset for AR6 is 28% earlier than AR12 while it is 15% earlier for AR9.
Primary trend emerges from the experiment shows that as the aspect ratio increases, the transition onset gets delayed. This corroborates with existing literature presented in this study.

Computational Results
The computational results presented here aim at highlighting the flow structure around different elliptical leading edges as well as highlight the performance differentials between higher AR leading edges that could not be tested experimentally. A plot of C f vs. the local Reynolds number Re x shown in Figure 9 indicates that the transition onset follows a sequential pattern with the transition onset being delayed as the aspect ratio increases. An interesting trend emerges from the plot shows that there are three distinct groups of elliptical leading edges. The first group consists of the standalone AR4. The second consists of AR6, AR9 and AR12 while the third group consists of the remaining two leading edges i.e.
AR16 and AR20. The key outtake from this trend is that the performance deficit between AR6, 9 and 12 is not significant and there might be no justification in choosing one over the other.  However, there is a much larger gap between the transition onset Reynolds number of AR12 and AR16 while AR16 and AR20 have a very marginal performance deficit. Thus, there would be a strong case in choosing AR16 over AR12. However, the choice between AR16 and AR20 is a difficult one due to the relatively small performance difference. Additionally, the transition onset Reynolds number for AR16 and AR20 are the closest to the transition onset Reynolds number for the ideal flat plate AR∞. This indicates that AR16 and AR20 have superior aerodynamic performance as compared to the other elliptical leading edges. The difference between the transition onset location between AR∞ and each individual leading edge as well as the spot propagation rates are given in Table 2.
The turbulent spot propagation rates given in Table 2 further exemplify that AR20 demonstrates the lowest spot propagation rates of all elliptical leading edges. With all other parameters being constant in all test cases, it can be stated that the spot propagation rate is inversely proportional to the aspect ratio of the leading edge. Higher AR leads to lower spot propagation rate and vice versa. Figure 10 shows the percentage difference between the transition onset for AR∞ and other leading edges. It is assumed that the transition onset percentage for AR∞ is 0. Hence it is not reflected on the graph. A logarithmic equation has been used to give a rough estimate of the percentage of early transition for different leading edges (Equation (5)). This will enable researchers to choose an appropriate leading edge based on their requirement. Journal of Transportation Technologies   edge provide regions of high receptivity [28]. The other key region for receptivity is the leading edge-flat plate juncture [7]. Pressure contours at nose of the leading edge indicate that the region foreward of the leading edge also known as the stagnation region is greatly impacted by the presence of the leading edge as shown in Figure 11. However, the influence of the leading edge diminishes with an increase in the AR. The stagnation region is the largest for AR9 and the smallest for AR20. This can be further verified by evaluating the flow structure through the analysis of the pressure gradient dP/dx for different leading edges as shown in Figure 12. From Figure 12(a), it can be observed that the adverse pressure increases with an increase in AR. However, the pressure gradient tends to zero rapidly as the AR increases. This is an indicator of the impact of the leading edge at the stagnation point. AR20 has the lowest peak favourable pressure gradient indicating the largest impact on the stagnation region albeit within a very small region. Conversely, AR9 has the largest favourable pressure gradient, but with a larger stagnation region. AR∞ being a plate with zero thickness and an infinite leading edge, the pressure gradient tends to remain constant with no stagnation point. Similarly, at the end of the leading edge as the AR increases the adverse pressure gradient decreases (Figure 12(b)). AR9 has the highest adverse pressure gradient while AR20 has the lowest adverse pressure gradient amongst all the elliptical leading edge. Higher adverse pressure gradient reduces the aerody-    (Figure 13(a)). Amongst all the leading edges, AR20 has the highest curvature at the start i.e. the nose of the leading edge. The curvature at the start decreases with a decrease in AR. Thus, the detrimental effect on the flow structure at the nose due to an increase in AR can be explained with the higher curvature values. At x/L = 0.05., AR 9 has the largest curvature. The curvature decreases with an increase in AR. By x/L = 0.5, AR 16 and AR20 show a very gradual curvature change and tend to indicate curvature values closer to the ideal leading edge with AR∞ as shown in Figure 13(b). By the end of the leading edge i.e. at x/L = 1, AR20 has a curvature value closest to AR∞.
The change in curvature dk will further highlight the differences between different aspect ratio leading edges and is shown in Figure 14. The curvature change is plotted along the length of the curve ds which is calculated using Equation (6).  From Figure 14, shows that at the start of the leading edge, the curvature change is more drastic for the higher AR leading edge, such as AR16 and AR20, but AR9 and AR12 show a more gradual change in the curvature along length of the leading edge. As the AR increases, the rate of change of curvature reduces along length of the leading edge. That means, as the flow progresses over the leading edge, the influence of curvature changes diminishes with an increase in AR. These curvature changes are high receptivity areas for disturbances to enter the boundary layer. Thus, the impact of a higher AR leading edge is to reduce the receptivity of the boundary layer to disturbances owing to a reduction in curvature discontinuities.
To further analyse the impact of the leading edge on the receptivity of the boundary layer, a plot of Cp* vs. x/L is shown in Figure 15. Cp* representing the ratio of the Cp of each individual leading edge to that of AR∞. The variable Cp* can be used to perform a comparative analysis between the relative performances of each individual leading edge. From Figure 15, it can be observed that as compared to the ideal leading edge AR4 has a very high peak Cp* value. This value reduces with an increase in the AR. The data for AR20 indicate that the impact of AR20 is the least and having a relatively low Cp* value indicating that its performance is closer to that of AR∞. Also, the data shows that the peak value for AR20 is obtained very close to the start of the leading

The Case of the Wedge vs. the Ellipse
One of the key observations from the result is that the aerodynamic performance of the leading-edge increases with an increase in AR. However, in practical experiment situations, higher AR leading edges, such as AR16 and AR20 might be difficult to manufacture. To justify the use of an elliptical leading edge, it is essential to compare it with the simplest leading-edge geometry i.e., the wedge.
Consequently, a wedge having the same length as AR20 is used and the performance of the wedge is compared against AR20 and AR∞ ( Figure 16). This study compares a wedge having the same dimensions as that of an AR20 elliptical leading edge. In case of the wedge, the two areas of discontinuities that introduce disturbances in flow would be at the nose and at the leading-edge flat plate juncture. The leading-edge flat plate juncture could possibly have the highest impact owing to a sharp edge at the leading-edge flat plate juncture. Table 3 provides enough evidence as to the merits of having an elliptical leading edge. The pressure gradient along the length of the leading edge shows that while the wedge has a very small stagnation region, the favourable pressure gradient for the wedge reduces very rapidly and it tends to zero after the stagnation region as shown in Figure 17. The impact of the elliptical leading edge is to maintain the favourable pressure gradient for longer.
A plot of the Cp* vs. x/L for the wedge along with its comparison to different elliptical leading edges is shown in Figure 18. An analysis of the plot indicates that at the start of the leading edge, the wedge has the smallest peak Cp* value    indicating that its impact on the flow at the start is closer to AR∞ as compared to the elliptical leading edges. For much of the length of the leading edge, the wedge maintains a competitive advantage over the elliptical leading edges. However, by x/L = 0.85, Cp* for the wedge is comparable to AR20 and AR16. By x/L = 1, the wedge shows Cp* values that are comparable to AR6. The trend for the Journal of Transportation Technologies wedge indicates that while, for most of the length of the leading edge, it tends to hold a performance advantage over other leading edges, but by end of the plate its advantage diminishes rapidly, and the elliptical leading edges tend to perform better. The advantage that the wedge has over the other AR leading edges over most of the leading edge ultimately enables it to have a performance which is better than AR12 but worse than AR16.

Discussions
The results presented in Section 3 demonstrate that as the AR increases, the The simulations show that AR16 and AR20 have better performance as compared to AR9 and AR12. The difference between AR16 and AR20 is a mere 0.3% while the difference between AR9 and AR20 is roughly 4% ( Table 2). When compared to an ideal case scenario of AR∞, AR12 has an earlier transition onset by roughly 6.2%. In the case of AR16 and AR20, this is around 3%. However, AR20 has the highest adverse pressure gradient at the stagnation region, but its effect is confined to a very small region.
The rate of change of curvature dK/ds and the pressure gradient dP/dx are interlinked. It can be observed that in case of AR20, the curvature is rapid in the beginning, but stabilises later. It is known that the receptivity of the boundary increases with an increase in curvature discontinuities. In the case of AR20, the sharp change in curvature is confined only to start of the leading edge, but the rate of curvature change is gradual later. Conversely, in the case of AR9, the curvature change at start of the leading edge is more gradual and the rate of curvature change does not approach a constant value until very late on the leading edge, thereby having a larger influence on the flow structure. This tends to Journal of Transportation Technologies have a knock-on effect on the transition onset location. The logarithmic function given in Figure 10 can be used to predict the transition offset between AR∞ and a particular leading edge. This would give prospective researchers a choice of an appropriate AR leading edge based on their research objectives.
Research requiring a high degree of accuracy can choose a higher AR leading edge, such as AR16 while experimental work having a higher tolerance could possibly use AR9 and AR12, because the difference to the ideal case scenario varies from 7% -15%.
The wedge discussed in this study (Section 3.3) also demonstrates that its performance is better than AR9 and AR12, but worse than AR16 and AR20. This could potentially provide a choice between choosing an elliptical leading edge or a wedge depending on the application. The wedge has an advantage owing to its simplistic geometric configuration that allows for easy manufacture and maintenance. Further study is needed to find the optimum leading edge for all test cases though. A potentially optimal solution would involve a wedge-ellipse combination however, because the wedge has a small stagnation region, and the elliptical leading edge has better fluid flow properties at the leading edge-flat plate juncture.

Conclusion
The effects of different leading edges on boundary layers subjected to a transitional regime at low free stream turbulence have been studied comprehensively.
The impact of changes in the AR of the elliptical leading edge is quite significant and the presence of a higher AR tends to delay the transition onset. In terms of transition onset location, it can be concluded that AR20 > AR16 > AR12 > AR9 > AR6 > AR4. However, at the start of the leading edge, the lower AR tend to have a larger stagnation region albeit with a lower impact on the fluid. As the flow develops over the leading edge, an increase in AR leads to better aerodynamic performance. From the present set of results, it can be stated that while AR20 is the best elliptical leading edge, AR16 would be the most optimum choice, because the difference in performance is negligible as compared to AR20, but it is better than AR12. Likewise, AR16 would be easy to machine as well. However, the presence of a wedge-shaped leading edge presents a strong case as an alternative to an elliptical leading edge. The authors recommend the use of wedgeshaped leading edge instead of lower AR leading edges in most experimental testing scenarios. Consequently, the choice of the leading edge is a trade-off between machinability and accuracy. Journal of Transportation Technologies simulations.

Availability of Data and Materials
The data and materials for this manuscript have been embedded within the manuscript itself. Datasets can be provided where necessary.