NACA 2412 Drag Reduction Using V-Shaped Riblets

: This research focuses on addressing a significant concern in the aviation industry, which is drag. The primary objective of this project is to achieve drag reduction through the implementation of riblets on a wing featuring the NACA 2412 aerofoil, operating at subsonic speeds. Riblets, with the flow direction on wing surfaces, have demonstrated the potential to effectively decrease drag in diverse applications. This investigation includes computational analysis within the ANSYS Workbench framework, employing a polyhedral mesh model. The scope of this research encompasses the analysis of both a conventional wing and a modified wing with riblets. A comparative analysis is conducted to assess variations in drag values between the two configurations. Parameters, including geometry, dimensions, and riblet placement at varying angles of attack, are explored to comprehend their impact on drag reduction. Notably, 15.6% and 23% reductions in drag were identified at a 16-degree angle of attack with midspan and three-riblet models, separately. The computational mesh and method were validated using appropriate techniques.


Introduction
The energy crisis in the current environment can be managed by developing renewable energy resources and reducing aircraft emissions [1].Kuhn et al. mentioned in their work that the aviation industry has some limitations on using renewable energies, such as hydrogen fuel and solar energy, and hence, it is essential to protect the available fuels through their proper and planned usage [2].In aircrafts, the fuel consumption can be reduced by implementing proper lift increment techniques by reducing the drag components and thereby increasing the lift-to-drag ratio, which is a measure of an aircraft's aerodynamic efficiency.Reducing aircraft drag is a crucial area of research in aerospace engineering due to its direct impact on aircraft performance, efficiency, and sustainability.Since data show that one count of drag increase (∆C d = 0.0001) is equal to two passenger weights in concord aircrafts [3], drag reduction methods are more relevant all the time.Among the various sources of drag on an aircraft, the wing has been a primary focus for developing drag reduction methods.Through extensive research, innovative techniques have been explored to reduce drag on aircraft wings, resulting in improved fuel efficiency, increased range, and reduced environmental impact.
There are many ways to control the drag, for example, increasing the Reynolds number can reduce the drag by minimizing the size of the boundary layer [4].However, there is a possibility of instabilities in the flow due to the implementation of boundary-layer decrement.The other way of reducing the drag is by changing the design of the aircraft wing surfaces.The study performed by Zeno et al. illustrated that the design change performed with anti-fairing reduced the drag [5].With the use of surface roughness elements, such as dimples, drag can be reduced by controlling the flow separation [6].There was a significant decrease in skin friction drag when a NACA 2412 aerofoil with Eng 2024, 5 945 variable-aspect-ratio dimples was examined.The three-dimensional inwardly dimpled wing at various angles of attack demonstrated that the air in the dimple surface accelerates, changing the boundary layer from laminar to turbulent.This delay in flow separation lowers drag.Even though it affects lift and drag, there are certain problems with the viability of creating dimples in an aircraft wing, and it also depends on where and how big the dimples are.It is not possible to apply this strategy to every aerofoil profile [7].
The main drag component, which contributes half of the total drag, is skin friction drag.It can be limited by employing laminar flow control that can be performed by using the hybrid laminar flow concept, which was conducted by applying suction in the leading-edge region to control the development of cross flow, and it helped to improve the aircraft performance.While it failed, the maintenance of this setup became a challenge since it was affected by insect debris and roughness, which affects its performance [8].Surface modifiers, such as vortex generators, were employed to mitigate drag by inducing turbulence without enhancing the coefficient of lift, which helped in postponing flow separation.However, the placement of these vortex generators is crucial, and they might not function effectively at lower angles of attack.Moreover, there is a potential for induced drag caused by trailing-edge vortices stemming from the use of vortex generators [9].
Another passive-separation-control technique using bumps employed in an SD7003 rectangular wing with a low Reynolds number provided a good improvement in the lift-todrag ratio.The bumps were positioned in close proximity to the leading edge of the wing to enhance aerodynamic performance by managing the bubbles resulting from laminar separation to avoid flow separation.Although it was possible to minimise the pressure drag forces, the degree of that reduction was heavily influenced by other factors, including the angle of attack and the number, placement, and size of the bumps [10].
The experiment conducted in a wind tunnel on a blunt wing with surface roughness gave a good picture of how flow behaviour influences drag.However, the amount of drag decrement and flow behaviour have a strong correlation with the leading-edge contour and roughness height [11].Large tripping in turbulent flow was seen in another study on flow behaviour using a spanwise row of cylindrical roughness components in a laminar, compressible three-dimensional boundary layer on a wing profile.It changed the location of the transition, causing an advance transition.However, it deliberately increased the flow's instabilities and turbulent effect [12].The study was carried out to examine how an aerofoil's flow separates when its upper surface is partially uneven.Tests were conducted in a subsonic wind tunnel at a 14-degree angle of attack.This demonstrated how the bumpy surface can be used to postpone the aerofoil's flow separation while increasing lift force.This experimental study has shown that modifications involving frequent disturbances, or "bumps," can postpone the flow separation on the aerofoil's surface.However, the actual implementation in aeroplanes is a little more difficult [13].Based on the review, it can be concluded that there is a significant lack of efficient drag reduction methods.
Utilizing riblets to close this gap is promising and has drawn a lot of interest.Due to their capacity in flow control dynamics and minimizing resistance, riblets have become recognized as a possible tool for lowering drag.This idea is extracted from sharkskin, which is made of tiny keels, as shown in Figure 1.Riblets come in many shapes and sizes based on the shark species.For example, Isurus oxyrinchus has a height of 0.012 mm and a spacing of 0.041 mm, and Carcharhinus limbatus has a height of 0.029 mm and a spacing of 0.065 mm [14].The flow separation on a staggered array of scales can lead to the scales becoming disrupted, causing a disturbance in the flow separation process and resulting in a delay in transition from laminar to turbulent.This helps to reduce the skin friction drag that emerges from the no-slip boundary effect.When the shark swims at speeds exceeding 10 m/s, the riblet induces turbulent flow over its skin due to the high Reynolds number.Even this turbulent boundary layer is less resistant to separation, leading to increased skin friction drag compared to a purely laminar flow scenario.The riblet mechanism can reduce 8 to 10% of the total skin friction drag [14][15][16][17][18][19].The effectiveness of drag reduction using riblet surfaces is impacted by a range of elements, encompassing factors like riblet geometry, Reynolds number, surface roughness, riblet alignment, interaction with the boundary layer, and many more.The main concept behind the implementation of riblets is making the surface rough, which reduces the amount of drag in turbulent regions, especially in stalling [21].Another study noted that the riblets lessen momentum transfer by reducing the vortices and cross flow velocity in the turbulent region.This influenced the turbulence to become less intense and postponed the transition to turbulence [22,23].Hence, in this research, an attempt was made to reduce the drag on a subsonic NACA 2412 wing using riblets.This aerofoil was selected because it is an asymmetrical aerofoil that offers lift, even at zero angle of attack, and still requires improvement.

Computational Design
The NACA 2412 aerofoil with a chord length of 127 mm [24,25] was used in both the wing root and the wing tip to develop a rectangular-shaped wing of wingspan 381 mm (762/2 = 381 mm), as it is one of the wing forms that produces more lift and less drag in the subsonic flow regime, and a constant taper ratio of 1:1 produces more lift than other wings [26].The wing model was designed using aerofoil coordinates with the help of AN-SYS Design Modeler 2022 R2 and Autodesk Inventor 2023 .A riblet's dimension selection holds significant importance, as heightened surface roughness corresponds to elevated skin friction.This outcome contributed to an overall rise in drag [27].The largest drag reduction occurred when the h/s ratio (h is the riblet's height, and s is the separation) for sawtooth riblets is 0.5, which was evidenced by Heidarian and Ghassemi, in which they looked at the spacing, height, and microchannels of the riblets in various flow velocities [28].The scalloped riblets with a h/s of 0.7 showed a maximum drag reduction of 5.1%, and the sawtooth riblets with a h/s of 0.5 showed a maximum drag reduction of 11% [28].Since this study shows the h/s = 0.5 ratio showed a good drag reduction, the initial design was performed with h/s = 0.5 (h = 0.002C = 0.254 mm, S = 0.00397C = 0.50419 mm).However, while performing the analysis, there was no proper improvement in drag reduction.So, a riblet with a h/s of 0.8 (h = 0.002c = 0.254 mm and s = 0.0025c = 0.3175 mm) was developed and used in this research.Yufei et al. noted in their research that the drag was reduced more effectively when the riblets were inclined at 15 degrees than when they were inclined at 30.They also observed that the drag increased when the angle increased from 20 to 30 degrees.Hence, in this project, the position and geometry of the riblet were focused on to improve the drag reduction.The riblet design was placed between 0.3c and 0.95c based on the previous study on the best location for drag reduction [29] and aligned parallel to the flow direction [30].Figure 2 shows the wing surface with four riblets positioned at the centre of the wingspan, each adhering to the dimensions previously described.The wing surface with riblets placed at three distinct positions, as depicted in Figure 3, was used to observe the variations in drag and lift values.Figure 4 provides a The effectiveness of drag reduction using riblet surfaces is impacted by a range of elements, encompassing factors like riblet geometry, Reynolds number, surface roughness, riblet alignment, interaction with the boundary layer, and many more.The main concept behind the implementation of riblets is making the surface rough, which reduces the amount of drag in turbulent regions, especially in stalling [21].Another study noted that the riblets lessen momentum transfer by reducing the vortices and cross flow velocity in the turbulent region.This influenced the turbulence to become less intense and postponed the transition to turbulence [22,23].Hence, in this research, an attempt was made to reduce the drag on a subsonic NACA 2412 wing using riblets.This aerofoil was selected because it is an asymmetrical aerofoil that offers lift, even at zero angle of attack, and still requires improvement.

Computational Design
The NACA 2412 aerofoil with a chord length of 127 mm [24,25] was used in both the wing root and the wing tip to develop a rectangular-shaped wing of wingspan 381 mm (762/2 = 381 mm), as it is one of the wing forms that produces more lift and less drag in the subsonic flow regime, and a constant taper ratio of 1:1 produces more lift than other wings [26].The wing model was designed using aerofoil coordinates with the help of ANSYS Design Modeler 2022 R2 and Autodesk Inventor 2023.A riblet's dimension selection holds significant importance, as heightened surface roughness corresponds to elevated skin friction.This outcome contributed to an overall rise in drag [27].The largest drag reduction occurred when the h/s ratio (h is the riblet's height, and s is the separation) for sawtooth riblets is 0.5, which was evidenced by Heidarian and Ghassemi, in which they looked at the spacing, height, and microchannels of the riblets in various flow velocities [28].The scalloped riblets with a h/s of 0.7 showed a maximum drag reduction of 5.1%, and the sawtooth riblets with a h/s of 0.5 showed a maximum drag reduction of 11% [28].Since this study shows the h/s = 0.5 ratio showed a good drag reduction, the initial design was performed with h/s = 0.5 (h = 0.002C = 0.254 mm, S = 0.00397C = 0.50419 mm).However, while performing the analysis, there was no proper improvement in drag reduction.So, a riblet with a h/s of 0.8 (h = 0.002c = 0.254 mm and s = 0.0025c = 0.3175 mm) was developed and used in this research.Yufei et al. noted in their research that the drag was reduced more effectively when the riblets were inclined at 15 degrees than when they were inclined at 30.They also observed that the drag increased when the angle increased from 20 to 30 degrees.Hence, in this project, the position and geometry of the riblet were focused on to improve the drag reduction.The riblet design was placed between 0.3c and 0.95c based on the previous study on the best location for drag reduction [29] and aligned parallel to the flow direction [30].Figure 2 shows the wing surface with four riblets positioned at the centre of the wingspan, each adhering to the dimensions previously described.The wing surface with riblets placed at three distinct positions, as depicted in Figure 3, was used to observe the variations in drag and lift values.Figure 4 provides a close view of the riblets, which looks in the form of isosceles triangles arranged spanwise on the wing surface.
close view of the riblets, which looks in the form of isosceles triangles arranged spanwise on the wing surface.

Computational Mesh
The impact of a computational domain was studied by comparing C, O, and rectangular types, and the C-type computational domain was chosen, as it suited the model's blunt leading edge and sharp trailing edge [31].This C domain was designed around the model based on 20 times the chord length of the wing [24] at the front and back, as shown

Computational Mesh
The impact of a computational domain was studied by comparing C, O, and rectangular types, and the C-type computational domain was chosen, as it suited the model's blunt leading edge and sharp trailing edge [31].This C domain was designed around the model based on 20 times the chord length of the wing [24] at the front and back, as shown in Figure 5, to capture the changes in the fluid flow.A rectangular-shaped body of influence (BOI) was introduced near the wing, approximately 100 mm away, with the intention of enhancing flow visualization by improving the accuracy of mesh within a particular region in the geometry.This involves refining the mesh in specific, significant areas, particularly those that exhibit intricate flow patterns or require heightened precision.This helps obtain a more precise representation of essential flow characteristics and gradients [32].

Computational Mesh
The impact of a computational domain was studied by comparing C, O, and rectangular types, and the C-type computational domain was chosen, as it suited the model's blunt leading edge and sharp trailing edge [31].This C domain was designed around the model based on 20 times the chord length of the wing [24] at the front and back, as shown in Figure 5, to capture the changes in the fluid flow.A rectangular-shaped body of influence (BOI) was introduced near the wing, approximately 100 mm away, with the intention of enhancing flow visualization by improving the accuracy of mesh within a particular region in the geometry.This involves refining the mesh in specific, significant areas, particularly those that exhibit intricate flow patterns or require heightened precision.This helps obtain a more precise representation of essential flow characteristics and gradients [32].After several attempts, it was found that the tetrahedral mesh could not provide the necessary quality of orthogonality and skewness, particularly in the trailing edge.When the mesh type changed to hexahedron, it was challenging to generate structure mesh in the corners, and the computational time of hexahedron mesh was high.Due to these issues encountered, the polyhedral mesh was chosen for the 3D analysis, as it can provide a high level of quality, can increase the numerical stability of the solution when compared to tetrahedral, and has fewer cells than other mesh types and hence displays fewer elements overall [33].The Y+ distance (between the boundary and the first grid's centre) is an additional parameter that must be taken into consideration when meshing because this helps find whether a mesh is fine or coarse.The wall Y+ was maintained between 30 and 300 based on the previous research [34] to suit the K-ε turbulence model with standard wall function [35] used to mesh the model.The mesh obtained had a skewness (deviation of the cell from its optimal cell size) in the range of 0.25 to 0.85 [36] and an orthogonal quality of above 0.1 [37].These values represented that a good mesh was made of reasonable quality.Figure 6 shows the mesh on the wing surface with riblets.In polyhedral mesh, the boundary conditions were required to be set before converting the surface mesh to solid mesh.
Setting the boundary conditions was necessary in polyhedral mesh before turning the surface mesh into solid mesh.Inlets, outlets, and symmetry boundaries were categorised using the ANSYS Update Boundaries option.Riblet surfaces, such as top, right, left, and aerofoil surfaces (leading edge and mid-surface without riblet), were considered walls.This enables the air to flow over the wall surface.If the rate of flow is known, the inlet can be called a velocity inlet and the pressure, which is unknown, can be called a pressure outlet [38].Thus, for this investigation, the pressure outlets and velocity inlet were used.
rised using the ANSYS Update Boundaries option.Riblet surfaces, such as top, right, left, and aerofoil surfaces (leading edge and mid-surface without riblet), were considered walls.This enables the air to flow over the wall surface.If the rate of flow is known, the inlet can be called a velocity inlet and the pressure, which is unknown, can be called a pressure outlet [38].Thus, for this investigation, the pressure outlets and velocity inlet were used.

Computational Analysis
Computational fluid dynamics (CFD) offers a numerical approximation to the equations that dictate the movement of fluids.The group of partial differential equations, such as Navier-Stokes equations, Euler equations, and Reynolds-averaged Navier-Stokes (RANS) equations, that characterize fluid flow are formulated [39].The simulation of the riblet wing is performed using a pressure-based solver because it is more effective in incompressible flows.The SIMPLE solution method was applied, as it effectively attains a steady-state solution when coupled with an implicit time treatment for flow variables.This method enhances the stability and convergence of the iterative process [40].
In this research, the K-Ꜫ turbulence model is used because it has been shown to be suitable for representing pressure distribution and illustrating the wake area [41], is easy to apply, and is valid for all types of flows [42].The standard K-Ꜫ model is a two-equation model (K = turbulent kinetic energy, Ꜫ = rate of dissipation of turbulent kinetic energy).When examining regular aerodynamic behaviours, steady flow is employed.In such instances, the right-hand components of continuity and diffusivity equations become zero,

Computational Analysis
Computational fluid dynamics (CFD) offers a numerical approximation to the equations that dictate the movement of fluids.The group of partial differential equations, such as Navier-Stokes equations, Euler equations, and Reynolds-averaged Navier-Stokes (RANS) equations, that characterize fluid flow are formulated [39].The simulation of the riblet wing is performed using a pressure-based solver because it is more effective in incompressible flows.The SIMPLE solution method was applied, as it effectively attains a steady-state solution when coupled with an implicit time treatment for flow variables.This method enhances the stability and convergence of the iterative process [40].
In this research, the K-ε turbulence model is used because it has been shown to be suitable for representing pressure distribution and illustrating the wake area [41], is easy to apply, and is valid for all types of flows [42].The standard K-ε model is a two-equation model (K = turbulent kinetic energy, ε = rate of dissipation of turbulent kinetic energy).When examining regular aerodynamic behaviours, steady flow is employed.In such instances, the right-hand components of continuity and diffusivity equations become zero, as there is no change with respect to time.This means that the flow parameters, such as pressure and velocity, remain constant and do not vary over time [43].By integrating and assessing the boundary condition constants, one can formulate steady-state flow equations.The steady-state flow equations are built upon the principles of uniform thickness and incompressible flow [41].The standard wall approach helps capture the effects of viscosity and turbulence near the wall, which are crucial factors in aerodynamic analysis [44].
The transport equation for turbulent kinetic energy is shown below: where P k = production of turbulent kinetic energy (TKE) due to mean velocity shear, P b = production of TKE due to buoyancy, S k = user-defined source, and σ k = turbulent Prandtl number for k.The transport equation for the turbulent dissipation rate is shown below: where C 1 , C 2 , and C 3 are coefficients specific to K-ϵ turbulence models that exhibit variation, S ϵ = user-defined source, and σ ϵ = turbulent Prandtl number for ϵ [45].Establishing a boundary condition during a CFD analysis influences the outcome since the choice of boundary condition relies on the given parameters.When the flow rate is a known factor, the inlet can be designated as a velocity inlet, while the pressure, which is to be determined, is assigned as a pressure outlet [38].The free stream velocity of this study is decided as 20.73 m/s, with a higher Reynolds number of 1.8 × 10 5 , as the flow is incompressible.The Mach number of this flow is calculated as M = velocity/speed of sound = 20.73/340= 0.06 [46].The flow is incompressible since the Mach number is less than 0.3.Additionally, it was identified that the flow is turbulent based on the calculation mentioned below [47]: The designed NACA 2412 wing and the wing with riblets placed at the midspan have a specific h/s ratio of 0.8 (as discussed in the computational design section) and were analysed using the K-ε viscous model, computationally, in airflow at different angles of attack, such as 0 • , 4 • , 8 • , 12 • , 16 • , 20 • , and 24 • .

Computational Result Analysis
The riblet wing model was analysed using both density-and pressure-based solvers, although the latter is more suitable for incompressible flows.The outcomes of these analyses were very close, as shown in Figure 7. Figure 8 shows a significant amount of drag reduction attained in the riblet wing between the 10 • and 20 • angles of attack.In order to analyse the drag value, the default output parameter, called the drag coefficient, was monitored and studied in the case of fluid analysis.The drag coefficient is derived from the drag equation [48], which is where D-drag force acting on the object; ρ-fluid density (here, this is considered to be the air density); v-the velocity between the object and fluid; C d -drag Coefficient S-surface area of the object.Despite being a dimensionless parameter, the drag coefficient scales the drag force according to the object's surface area, flow velocity, and flow density.The shape of the object, the roughness of the surface, and the fluid characteristics are other variables that a drag coefficient may impact.
The streamlined shape of the wing was the main focus of the design since a low drag coefficient reduces air resistance between the item and its surroundings [49].While using a riblet wing, a larger amount of 15.6% drag reduction was observed in the riblet wing at the 16 • angle of attack.
Because of the relatively low velocity between the riblets and the high velocity experienced at the riblet tip due to its specific localised location, the riblet behaviour at the 0 • angle of attack demonstrates a 5.7% reduction in total drag.Figure 9 clearly shows that the level of velocities has been increased on the top surface, especially the areas where the riblet was implemented.This increment in the velocity on the top surface of a riblet-implemented wing favoured the drag reduction, as expected.This reduction in drag can be attributed not only to a decrease in turbulent friction but also to the mitigation of the separation area and otherwise smooth trailing section of the object [50].
rienced at the riblet tip due to its specific localised location, the riblet behaviour at the 0° angle of attack demonstrates a 5.7% reduction in total drag.Figure 9 clearly shows that the level of velocities has been increased on the top surface, especially the areas where the riblet was implemented.This increment in the velocity on the top surface of a riblet-implemented wing favoured the drag reduction, as expected.This reduction in drag can be attributed not only to a decrease in turbulent friction but also to the mitigation of the separation area and otherwise smooth trailing section of the object [50].Several studies demonstrate that symmetrical V-shaped riblets have the capability to enhance lift [51].Therefore, the V-shaped riblets were placed at the midsection of the wingspan, ranging from the 0.3C to 0.95C chord position.However, this did not improve the lift coefficients as expected (shown in Figure 10) due to variations in the height, spacing between the riblet tips, utilization of lower stream velocities, and riblet dimensions when compared to the referenced literature [51].As a result, a change was made by relocating the riblets from midspan to three different positions, as depicted in Figure 3 under the computational design section.Identical meshing and analysis approaches were  Several studies demonstrate that symmetrical V-shaped riblets have the capability to enhance lift [51].Therefore, the V-shaped riblets were placed at the midsection of the wingspan, ranging from the 0.3C to 0.95C chord position.However, this did not improve the lift coefficients as expected (shown in Figure 10) due to variations in the height, spacing between the riblet tips, utilization of lower stream velocities, and riblet dimensions when compared to the referenced literature [51].As a result, a change was made by relocating the riblets from midspan to three different positions, as depicted in Figure 3 under the computational design section.Identical meshing and analysis approaches were Several studies demonstrate that symmetrical V-shaped riblets have the capability to enhance lift [51].Therefore, the V-shaped riblets were placed at the midsection of the wingspan, from the 0.3C to 0.95C chord position.However, this did not improve the lift coefficients as expected (shown in Figure 10) due to variations in the height, spacing between the riblet tips, utilization of lower stream velocities, and riblet dimensions when compared to the referenced literature [51].As a result, a change was made by relocating the riblets from midspan to three different positions, as depicted in Figure 3 under the computational design section.Identical meshing and analysis approaches were adopted for this modified configuration, which resulted in improved lift values, as shown in Figure 11.A substantial lift enhancement of 19%, 9%, and 20% were noted at the 4 • , 8 • , and 12 • angles of attack, respectively.Moreover, a remarkable lift increment of 23% was achieved at the 16 • angle of attack.It is essential to point out that the results may vary at different flow conditions, as th pressure distribution varies at different flow speeds.Additionally, it is not possible to gen eralize the results obtained for all 3D aerofoils, as the chord and camber varies from on aerofoil to another.

Method Validation
It is essential to verify the computational method used in this research project.Hence    It is essential to point out that the results may vary at different flow conditions, as the pressure distribution varies at different flow speeds.Additionally, it is not possible to generalize the results obtained for all 3D aerofoils, as the chord and camber varies from one aerofoil to another.

Method Validation
It is essential to verify the computational method used in this research project.Hence, relevant literature [24] (where the actual lift and drag values were noted) was taken as a reference to validate the same method at the start of this project.A 3D design of the NACA  It is essential to point out that the results may vary at different flow conditions, as the pressure distribution varies at different flow speeds.Additionally, it is not possible to generalize the results obtained for all 3D aerofoils, as the chord and camber varies from one aerofoil to another.

Validation 6.1. Method Validation
It is essential to verify the computational method used in this research project.Hence, relevant literature [24] (where the actual lift and drag values were noted) was taken as a reference to validate the same method at the start of this project.A 3D design of the NACA 2412 wing was developed and analysed in airflow using the K-epsilon and the RSM models, similar to the procedure used in the reference journal.Mesh convergence was achieved at 819,854 elements, which is very close to the number of elements (fewer than 800,000 elements) noted in the reference journal/literature.
In terms of drag coefficient, similar values were obtained at the 0 • and 4 • angles of attack, as shown in Figure 12.For example, at 4 • , 0.059 and 0.065 drag coefficients were noted in the recreated design and the literature.Although slightly deviated values were observed at the 8 • , 12 • , and 16 • angles of attack, still those values are within the acceptable range.This verifies that the computational method used in this project is reliable and relatable, as the results are validated against actual values.
Eng 2024, 5, FOR PEER REVIEW 11 observed at the 8 0 , 12 0 , and 16 0 angles of attack, still those values are within the acceptable range.This verifies that the computational method used in this project is reliable and relatable, as the results are validated against actual values.

Mesh Validation
For precise findings, an ideal grid design is needed, and to create the best grid, the grid independence test is typically run [52].Thus, it was performed for the validation of the NACA2412 aerofoil in 3D.The mesh element size was changed repeatedly to obtain a stable value for the lift coefficient.Figure 13 interprets the lift coefficient noted at different mesh element sizes.Exactly six mesh independence studies (MISs) were carried out to obtain the expected results.Two MISs were performed under 50,000 elements, which resulted in different lift coefficients.Another two MISs were performed between 50,000 and 100,000 elements and obtained unstable results.Finally, the solution reached a state of stability, resulting in a lift coefficient of 0.18, obtained with 819,854 elements for a smooth wing, as shown in Figure 14.Similarly, an MIS was performed for the riblet wing, as shown in Figure 14, and noted stable results.The elements that showed stable results (in normal and riblet wings) were used in this research project.

Mesh Validation
For precise findings, an ideal grid design is needed, and to create the best grid, the grid independence test is typically run [52].Thus, it was performed for the validation of the NACA2412 aerofoil in 3D.The mesh element size was changed repeatedly to obtain a stable value for the lift coefficient.Figure 13 interprets the lift coefficient noted at different mesh element sizes.Exactly six mesh independence studies (MISs) were carried out to obtain the expected results.Two MISs were performed under 50,000 elements, which resulted in different lift coefficients.Another two MISs were performed between 50,000 and 100,000 elements and obtained unstable results.Finally, the solution reached a state of stability, resulting in a lift coefficient of 0.18, obtained with 819,854 elements for a smooth wing, as shown in Figure 14.Similarly, an MIS was performed for the riblet wing, as shown in Figure 14, and noted stable results.The elements that showed stable results (in normal and riblet wings) were used in this research project.

Conclusions
The drag reduction efficiency of a V-shaped triangular-type riblet was studied with a pressure-based solver in turbulent flow.Computational analysis was performed at different angles of attack, and the results were compared.The riblet design at the centre of the wingspan from the 0.3 to 0.95 chord position showed effective drag reduction from 5.7 to 15.6% at several angles of attack.This experiment aimed to enhance lift by incorporating

Conclusions
The drag reduction efficiency of a V-shaped triangular-type riblet was studied with a pressure-based solver in turbulent flow.Computational analysis was performed at different angles of attack, and the results were compared.The riblet design at the centre of the wingspan from the 0.3 to 0.95 chord position showed effective drag reduction from 5.7 to 15.6% at several angles of attack.This experiment aimed to enhance lift by incorporating

Figure 2 .
Figure 2. Riblets at the centre of the wingspan.

Figure 3 .
Figure 3. Riblets at three different places of the wingspan.

Figure 2 .
Figure 2. Riblets at the centre of the wingspan.

Figure 2 .
Figure 2. Riblets at the centre of the wingspan.

Figure 3 .
Figure 3. Riblets at three different places of the wingspan.Figure 3. Riblets at three different places of the wingspan.

Figure 3 .Figure 4 .
Figure 3. Riblets at three different places of the wingspan.Figure 3. Riblets at three different places of the wingspan.Eng 2024, 5, FOR PEER REVIEW 5