Enhancing the Performance of an Oscillating Wing Energy Harvester Using a Leading-Edge Flap

: In this study, we investigated the power generation capability of an oscillating wing energy harvester featuring an actively controlled ﬂap positioned at the wing’s leading edge. The ﬁndings revealed that a�aching a leading-edge ﬂap reduces ﬂuid ﬂow separation below the wing’s lower surface at the leading edge, resulting in smoother flow and increased velocity near the hinge region. The leading-edge flap increases the pressure difference across the wing’s surface, thereby enhancing the overall performance. In addition, the introduction of the leading-edge flap effectively elon-gates the wing’s effective projected length in the heaving direction, leading to increased thrust. We examined flap lengths ranging from 10% to 50% of the chord length, with the maximum pitch angles of the wing and flap varying from 75° to 105° and 30° to 55°, respectively. The optimal power generation was achieved using a flap length of 40% of the chord length, combined with maximum wing and flap pitch angles of 95° and 45°, respectively. These conditions yielded a 29.9% overall power output increase and a 20.2% efficiency improvement compared to the case without the leading-edge flap.


Introduction
The global shift toward decreasing reliance on fossil fuels has motivated the investigation of alternative energy sources and methods, with a focus on enhancing their effectiveness.Energy harvesters based on flapping foils, inspired from natural movements, have gained traction as an effective means to extract energy from various natural sources, such as streams, rivers, tidal currents, and wind.This approach has garnered substantial interest in recent years [1,2].The idea of utilizing a flapping motion for energy generation was initially introduced by Wu in 1972 [3], and later, in 1981, McKinney and Delaurier pioneered the extraction of energy from the heaving and pitching motions of a fluid using a flapping foil [4].Subsequently, extensive research has been conducted on optimizing the energy-harvesting capabilities of flapping foils via the investigation of motion parameters [1,5,6], geometric and viscous parameters [2,[7][8][9][10], and flow control methods [11][12][13][14][15].
Xiao et al. [6] introduced an innovative non-sinusoidal motion for oscillating hydrofoils within energy harvesters.They combined sinusoidal plunging movements with a nearly trapezoidal pitching motion to assess the impact on energy extraction performance.They found that a higher angle of a ack results in enhanced power extraction.In comparison to sinusoidal motion, they observed substantial enhancements with the optimal pitching profile, boosting the power coefficient by up to 63% and the total efficiency by up to 50%.Wu et al. [16] investigated the effects of stroke deviation and variables such as the horizontal motion amplitude, phase difference, and frequency on the performance of flapping-foil energy harvesters.Their results emphasized that incorporating horizontal motion enhances the lift force, thereby improving the power generation capacity of the energy harvesters.Their study underscored the enhancement of power output performance owing to flapping-foil energy harvesters.Wang et al. [17] proposed a novel reversed-D trajectory and concluded that this specific trajectory motion can significantly enhance the power output.Swain et al. [18] investigated energy extraction performance by altering the flapping trajectory and the spatial arrangement of foils, identifying favorable wake interaction pa erns that contribute to increased efficiency.
Shanmugam and Sohn [19] investigated various deflector designs aimed at improving power generation in flapping-foil energy harvesters.They examined the impact of upstream deflectors (N), the angle of inclination of these deflectors, and the spacing between the hydrofoil and the deflector on the performance.Their simulations highlighted the considerable influence of these parameters on the power generation efficiency of the energy harvesters.Wang et al. [20] analyzed the impact of density-stratified flow on energy extraction efficiency and its interaction with the pitching amplitude.They [21] proposed a tandem-hydrofoil-based tidal array with improved density and reduced costs.Petikidis and Papadakis [22] studied fully passive flapping foils in free surface flow, determining optimal submergence depths and assessing the influence of monochromatic waves.Dahmani and Sohn [23] introduced a novel approach using oscillating tandem wings inside a convergent duct to enhance power extraction in energy harvesting.The downstream hydrofoil was shown to be more influential, contributing 4-80% to the improved power output owing to the high incoming fluid velocity due to the duct design.In addition, they explored the impact of vertical interfoil spacing on tandem/parallel coupled oscillating wings, achieving a 23% increase in power extraction efficiency compared to traditional arrangements [24].In these systems, maintaining stability is important, as it ensures consistent operation and reduces the risk of performance degradation.Zhao [25][26][27][28] explored the principles of stability as they pertain to the design and analysis of such systems.Building on these insights, recent advancements in piezoelectric harvesters, such as the development of a harvester with a U-shaped geometry [29], have led to significant improvements.This proposed model surpasses traditional models in efficiency, showcasing the evolution of energy-harvesting technology.Further contributing to this field, the creation of a semi-submersible piezoelectric harvester optimized for increased power output [30] underscores the growing potential of these technologies in a range of applications, driven by continuous improvements in design and optimization.
Flow control methods for harnessing energy from flapping-foil systems can be categorized into two types, passive control and active control, based on whether additional energy is required within the flow field [31].Passive control techniques involve incorporating simple structures on or adjacent to the surface of the airfoil to enhance energy harvesting without introducing extra energy.Examples of passive control methods include passive foil deformation [32], Gurney flaps [33,34], and corrugated foils [35].Passive control methods are cost-effective and relatively straightforward to install.However, they operate within predefined states and lack adaptability for real-time adjustments to meet changing requirements.In contrast, active control methods such as circulation control [36] and plasma actuators [37,38] can modify the flow field by adjusting excitation parameters as needed, making them more efficient compared to passive control methods.
Flaps are commonly used in lifting mechanisms in the aerospace industry owing to their simple structure, robustness, and effectiveness [39,40].Flaps can be implemented in both passively and actively controlled forms.Although flaps have been extensively studied for their applications in single blades and vertical-axis wind turbines [41][42][43], the investigation into energy harvesting using flapping wings has mainly focused on Gurney flaps [33,34].For example, Bing [33] revealed the introduction of a Gurney flap-influenced vortex generation at the trailing edge, which increased the pressure difference on the foil surface and improved the lift force, leading to a 21% improvement in energy-harvesting efficiency.On the other hand, any additional energy consumption resulting from the swing of the Gurney flap itself was not mentioned.Totpal et al. [44] studied the impacts of passive leading-edge flaps on flapping-wing energy harvesters at low cutoff frequencies.Although they observed an improvement in the heave force, this improvement was primarily limited to the initial stages of the flapping cycle when the heave speed was low, and the synchronization between force and speed was poor, making overall flappingwing control challenging.Alam and Sohn investigated [45] the impact of an actively controlled trailing-edge flap, accompanied by a comprehensive parametric analysis [46].Their findings indicated that the implementation of the flap yielded advantageous outcomes in terms of enhancing power extraction performance.
In this study, we examined the potential for enhancing energy harvesting in oscillating wing systems via the addition of a controllable leading-edge flap.Our primary goal was to improve the performance of the energy harvester.To achieve this, we conducted a quantitative evaluation of numerous factors, including various flap lengths and maximum pitching angles for the wing and flap, to identify the optimized configurations that maximize power generation and enhance overall efficiency.Flap lengths varying from 10% to 50% of the chord length were systematically investigated.In addition, variations in maximum wing pitch angles spanning from 75° to 105° and maximum flap pitch angles ranging from 30° to 55° were examined.Through this comprehensive investigation, we aimed to reveal the relationships among these diverse parameters and their impacts on energyharvesting performance.

Problem Description
We investigated the dynamics of unsteady fluid flow in an oscillating wing energy harvester using a two-dimensional model under turbulent flow conditions by examining the impact of a leading-edge flap on the performance of energy harvesting.Figure 1a presents a two-dimensional (2D) representation of the energy harvester with and without the flap.In addition, Figure 1b shows three-dimensional (3D) visualizations to view this configuration more clearly.This research primarily focuses on a comparative analysis of the wing's performance with and without the leading-edge flap, shedding light on the role of the leading-edge flap in boosting energy-harvesting efficiency.Table 1 provides an overview of the physical model and the optimal parameter settings needed to maximize the power generation in an oscillating wing energy harvester without a flap, as reported by Usoh et al. [9].The oscillating wing model, when equipped with a leading-edge flap, maintains identical physical characteristics (e.g., thickness and total length).However, factors such as the length of the flap and the peak pitch angles of both the wing and flap are varied.The wing and flap are separated by a gap equivalent to half a percent of the wing's chord length, and the junction between the wing and flap features a radius equal to half the plate thickness.Figure 2 shows the oscillating wing kinematic in two scenarios: a wing without a flap (2a) and another wing with a leading-edge flap (2b).The incoming fluid velocity is U∞, and the wing's chord length is c.The wing's main body is pivotally a ached at one-third of the chord length from the leading edge, enabling both heaving and pitching motions.Meanwhile, the flap, hinged at the wing's leading edge, executes pitching motions relative to the wing's main body.Notably, changes in the flap's pitching profile, dictated by the angle of incidence, significantly affect the energy harvester's performance.To evaluate the impact of a aching a leading-edge flap to an oscillating wing, its power generation capacity is compared with the maximum power output of a wing without a flap based on the optimum conditions obtained for a flapless wing outlined in Table 1 [9].The wing's heaving motion is described by Equation ( 1): In this equation, h denotes the instantaneous amplitude of heaving motion, ω is the angular frequency (expressed as ω = 2πf) associated with the wing's oscillation, t is the time, and φ is the heave and pitch motion phase difference.
The pitch motion of the wing's main body is governed by Equation (2): Here, θ represents the pitch angle of the main body.Similarly, the pitch motion of the flap is expressed by Equation (3): Here, ψ represents the pitch angle of the flap.The force ( ) and the moment ( ) are associated with heaving and pitching relative to the wing's pivot point, respectively.The corresponding average power outputs derived from these forces and moments are represented by and , respectively, and represents the average total power output [33,34].
In this equation, ( ) denotes the instantaneous heaving motion of the wing, and ( ) represents the angular velocity at time .The coefficients of heaving and pitching power are represented by Cpy and Cpm, respectively.The total averaged power extraction coefficient in a cycle is expressed as ̅ and can be calculated by summing up the aver- aged power output coefficients obtained from ( ) and ( ) [47,48]: Considering the combined effects of the main body and flap, the forces, moments, and power generated by both components were summed.The computation of the moment power for the flap incorporated the consideration of relative angular motion between the wing's main body and flap.The efficiency of power generation in the oscillating wing energy harvester is denoted as η, obtained using Equation ( 6) [49]:

Numerical Method
The power extraction efficiency of the oscillating wing was simulated using the incompressible unsteady Reynolds-averaged Navier-Stokes equations.The k-w SST turbulence model was employed, as adopted by other researchers in similar studies [17,48].The computational model was solved using ANSYS commercial software version 20 (ANSYS, Canonsburg, PA, USA) and the overset method by dividing it into stationary and moving grid sections, as illustrated in Figure 3.To ensure domain independence, a computational domain spanning from −20 50 in length to −25 25 in breadth was selected based on the observation that enlarging the domain size had an insignificant impact on the obtained results.At the inlet, outlet, and wing surface, the boundary conditions, velocity inlet, pressure outlet, and no-slip wall conditions were implemented.The pressure outlet was placed at 50 from the oscillating wing pivot point, as depicted in Figure 3a.Separate subgrids (moving grid) were utilized for the main body and the flap, as shown in Figure 3b.The movable grid areas for each part were organized into rectangular grids measuring 6c and 5c.For pressure, momentum, and turbulent viscosity, the spatial discretization was conducted using the second-order upwind approach.The dynamic mesh method facilitated mesh motion, and the airfoil's movement was input through user-defined functions.The presented analysis shows the data collected after the completion of the eighth stable periodic cycle, which included 2000 time steps.
To investigate the precision and dependability of the numerical simulations, studies on grid and time independence were conducted by altering the grid resolution and time step.Three distinct mesh densities were employed: a coarse grid (D1), a medium grid (D2), and a fine grid (D3).These simulations were performed at a Reynolds number of 5 × 10 5 using a chord length of 1 and a pivot point located at c/3.The heave amplitude was set at c = Ho = 1, the maximum pitching amplitude, at 75°; the phase angle, at 90°; and the reduced frequency, at 0.14.Table 2 illustrates the influence of grid resolution and time step on the computed value of ̅ .Based on the outcomes of the grid independence assessment, the D2 grid with a time step of T/2000 was chosen for the simulations to ensure both accuracy and reliability.Further, we validated the results obtained upon comparison with the simulation results reported by Usoh et al. [9] and Kinsey et al. [48].The validation process involved examining a flat plate and a NACA0015 airfoil.For the flat plate scenario, the simulation was conducted at Re = 1100, Ho/c = 1.0, f* = 0.14, and θo = 75°.The calculated Cpt values closely matched the results reported by Usoh et al. [9], as depicted in Figure 4a.For the NACA0015 airfoil, shown in Figure 4b, the simulations at Re = 5 × 10 5 , Ho/c = 1.0, θo = 75°, and f* = 0.14 demonstrated a good match between the pushing force coefficient and the xdirection force coefficient in this study and those documented by Kinsey et al. [50].These comparisons validate the applied numerical model and mesh used in this study.

Results and Discussion
Our investigation presents an active flow control methodology aimed at improving the power generation of an oscillating wing.The effects of a rectangular-plate oscillating wing featuring a flap at the leading edge on the energy-harvesting performance are investigated during coupled heaving and pitching movements.The study further evaluates the wider applicability of this flow control technique (over a range of motion parameters) by varying multiple parameters, such as the maximum pitch angles of the wing and the flap, as well as the flap length.Throughout the study, the following parametric conditions were kept constant: the nondimensional pivot point position (xp) = 0.33c, the heave amplitude (Ho/c) = 1.0, the phase angle (φ) = 90°, and the frequency (f*) = 0.14.

Effect of Varied Maximum Flap Pitch Angles on Power Output
In this section, a comparative analysis of the energy-harvesting performance of a wing with an integrated flap is performed against a baseline foil without a flap at a fixed maximum wing pitch angle (θo) of 75°.The flap length is maintained at 10% of the chord length (c), and the maximum flap pitch angle (ψo) is varied from 20° to 50°.
Table 3 shows the average pushing power coefficient ( ̅ ), average moment power coefficient ( ̅ ), and average total power coefficient ( ̅ ) for the wing with a flap length of 10% of c at θo = 75°.The power coefficients are tabulated against a range of ψo values, along with the percentage change in ̅ relative to the wing without a flap (denoted as ̅ ).In this context, "plate" refers to the baseline wing without a flap, and the varying flap conditions are indicated by their respective maximum pitch angles, such as ψo = 20°, 25°, and so on.The data exhibit a direct correlation between ψo and ̅ , suggesting that an increase in the flap's pitch angle enhances the average heaving power coefficient.Conversely, an inverse relationship is observed between ψo and ̅ , where larger flap pitch angles correlate with more negative values of the average moment power coefficient.̅ displays a nonlinear association with ψo: it shows a modest increase as ψo ascends from 20° to 40°, followed by a decrease in ̅ beyond 40°.Notably, the power output increases by a maximum of 6% with a 10% flap length compared to the plate case.5a,b show the variations in the pushing force coefficient (Cy) and the pushing power coefficient (Cpy) in time (T) for different ψo cases outlined in Table 3. Notably, between the time intervals from 0.1 t/T to 0.3 t/T, both Cy and Cpy increase as ψo increases, indicating that the oscillating foil experiences enhanced pushing force and an increase in pushing power as ψo increases during this interval.Figure 5c,d present the moment power coefficient (Cpm) and the total power coefficient (Cpt) throughout a cycle considering different ψo cases.The variation in Cpm values is similar for different cases until 0.3 t/T, whereas from 0.3 t/T to 0.5 t/T, the Cpm values decrease with ψo.During the interval from 0.15 t/T to 0.35 t/T, the Cpt values also increase with increasing ψo.The total power output of the oscillating foil shows an upward trend with increasing ψo values during this interval.However, the maximum ̅ is achieved at ψo = 40°, as detailed in Table 3.In Figure 6, streamline and vorticity, velocity magnitude, and pressure contour plots around the foil's surface are presented for various ψo cases at 0.15 t/T and θo = 75°.For a flat plate, fluid separation occurs below the leading edge on the lower surface of the wing, as illustrated in Figure 6a.However, this separation is significantly reduced at higher ψo values by a aching a leading-edge flap to the wing, resulting in smoother fluid flow and increased velocity near the hinge region on the lower surface, as depicted in Figure 6b. Figure 6c displays the pressure contour plot, showing that the leading-edge flap reduces the lower surface pressure compared to the wing without a flap.In Figure 7, the pressure coefficient on the wing's surface along the heaving motion is illustrated for various ψo cases at 0.15 t/T and θ0 = 75°.The presence of a flap decreases the surface pressure owing to the expanded high-velocity region in comparison to the plate, resulting in a greater pressure differential between the upper and lower surfaces and generating a stronger pushing force in the wing, with higher ψo values compared to the plate.In addition, the flap notably increases the projected length in the wing's heaving motion, positively influencing the pushing force at this specific moment.

Impact of Flap Length Variation on Power Output
To examine the role of the leading-edge flap length on the power output of the energy harvester, this section maintains constant θο = 75° and ψo = 40°.The study evaluates changes in flap length relative to the chord length considering lengths of 0% (equivalent to a flat plate), 10%, 15%, 20%, 33% (pivot point), and 40%.
Table 4 summarizes the average coefficients of power ( ̅ , ̅ , ̅ ) and the percentage change in ̅ (denoted as ̅ ) for the wing with different flap lengths at fixed θo = 75° and ψo = 40°.The ̅ values are maximum for a flap span of 20% of the chord length, beyond which the values start decreasing.Conversely, the magnitude of the ̅ value steadily increases with the extension of the flap length.For a flap length of 10% of the chord, the ̅ value is higher than the flat plate configuration but decreases progressively as the flap span increases.However, the ̅ value steadily decreases as the flap span extends up to 20% of the chord length, experiencing a significant drop beyond this threshold.This indicates that beyond a certain flap length, an increase in flap length adversely affects the power output generation of the energy harvester system for the same θo and Ψo. Figure 8 illustrates the effects of different flap lengths on the pushing force (Cy), pushing power (Cpy), moment power (Cpm), and total power (Cpt) coefficients within an oscillation cycle at fixed pitch angles (θο = 75° and ψo = 40°).Figure 8a displays the variation in Cy, revealing that adding a flap length of 10% of c increases Cy relative to the baseline "plate" case from 0.15 t/T to 0.4 t/T.This increase becomes more pronounced with a flap span of 15% of c, and it continues to grow up to 20% of c.Beyond this flap length, the Cy values plateau, suggesting a threshold for enhancement through flap extension during this interval.Conversely, the interval from 0.4 t/T to 0.5 t/T is characterized by a marked and consistent reduction in Cy as the flap lengths continue to increase.Given that the wing's vertical velocity (Vy) is higher from 0.1 t/T to 0.4 t/T, even a modest increase in Cy due to changes in flap length has a significant impact on Cpy, as depicted in Figure 8b. Figure 8c indicates that the magnitude of the negative Cpm values intensifies with the increment in the flap length up to 0.25 t/T, suggesting an increase in adverse moment power.Conversely, the positive Cpm values decline from 0.25 t/T to 0.5 t/T, implying a reduction in favorable moment power during this period, particularly for longer flaps.This corresponds to a notable increase in the average negative moment power for larger flap lengths, consistent with the data in Table 4.As for the total power coefficient displayed in Figure 8d, there is an observable increase in Cpt with the lengthening of the flap from 0.1 t/T to 0.3 t/T.However, this trend is reversed from 0.3 t/T to 0.5 t/T, where Cpt begins to decrease.The results suggest that the total power output initially benefits from the increased flap length, which peaks at a certain point before starting to diminish, suggesting that an optimal flap length exists for maximizing the total power output under the specific pitch angles of θο = 75° and ψo = 40°.Figure 10 displays the variation in pressure across the surfaces of the wing for various flap lengths at 0.25t/T.In Figure 10, with a flap span of 10% of c, the wing's lower surface pressure is observed to be lower than that of the base plate configuration, a ributed to fast fluid flow near the hinge position, where the camber becomes the maximum.As the flap length is extended to 15% of c, a further pressure decrease is observed.However, beyond a flap span of 20% of c, the lower surface pressure begins to increase as flow separation becomes more pronounced, suggesting that up to a certain flap length, the incoming flow stays a ached to the wing's surface, reducing pressure on the lower surface.Conversely, when the flap is excessively long, the flow separates earlier, causing an increase in pressure.Therefore, the variations in pressure across the surfaces of the wing decrease beyond a flap span of 20% of c.In addition, with an increase in flap length, there is also an increase in the flap's projected length along the x-axis, which is beneficial to enhancing the pushing force, offse ing the effects of the reduced pressure differential.

Effect of Varied Maximum Wing Pitch Angles
This section examines the influence of altering the maximum wing pitch angle from 75° to 95° while maintaining a constant flap length of 20% of c and a maximum flap pitch angle of 40°.The values in Table 5 reflect the average pushing power ( ̅ py), average moment power ( ̅ pm), and average total power output ( ̅ pt) to assess the performance across varied θo.̅ py increases from θo = 75° to 85°, reaching its peak at θo = 85°, and it decreases thereafter, suggesting that a wing angle of 85° is optimal for generating the pushing power at this operating condition.Conversely, the absolute ̅ pm values exhibit a continuous decrease as θo increases.On the other hand, the variations in ̅ show an upward trend with higher wing angles, indicating that greater wing angles enhance power generation for a flap span of 20% of c.These findings emphasize the importance of selecting an appropriate wing angle to maximize the power output of the oscillating wing with a fixed flap length.A pitch angle θo of 90° strikes a balance between pushing power and moment power requirements, resulting in the optimum power generation at a flap span of 20% of c and ψo = 40°.3.2 10.9 13.9 14.9 9.5 Figure 11 provides the pushing force, pitching moment, and power coefficients in time T for a flap length of 20% of c, ψo = 40°, and varied θo. and Cpy variations over a cycle for different wing pitch angles spanning from 75° to 95° in 5° increments.Notably, from 0.1 t/T to 0.4 t/T, Cy and Cpy exhibit higher values when the wing's maximum pitch angle is lower.However, from 0.4 t/T to 0.5 t/T, Cy and Cpy show an increase with higher maximum wing pitch angles.In addition, the absolute Cpm consistently increases up to 0.25 t/T, indicating an increase in adverse moment power for longer flaps.From 0.25 t/T to 0.5 t/T, the trend of increasing Cpm values continues, suggesting an increase in favorable moment power for the wing, as shown in Figure 11c.Notably, the rise in adverse moment power is less pronounced than that in favorable moment power, resulting in a marked decrease in the average negative moment power with higher θo values, which agrees with the values presented in Table 5.Interestingly, the examina- tion of Cpt values reveals a decreasing trend from 0.1 t/T to 0.4 t/T with a maximum increasing wing pitch angle, followed by an increase from 0.4 t/T to 0.5 t/T as the maximum wing pitch angle increases, as shown in Figure 11d.Figure 12 presents the streamline and velocity magnitude and pressure contour plots for a wing with a leading-edge flap length of 20% of c and ψo = 40° based on varying θo at intervals of 0.25 t/T.In Figure 12a, the streamline and velocity magnitude contours show a notable increase in velocity around the flap's hinge area on the lower surface as θo increases.Figure 12b illustrates that the wing's upper surface pressure significantly increases for θo = 90°.Figure 13 depicts the distribution of pressure coefficients on the foil surface across different θo cases at 0.25 t/T. Figure 13 displays a continuous decrease in the wing's lower surface pressure, particularly near the flap-wing junction, with increasing θo.As θo increases, the wing's projected length along the heaving motion direction decreases, accompanied by an increase in the pressure differential across the wing surface.

Combined Impact of Varying Flap Length and Maximum Pitch Angles of the Wing and Flap
Figure 14 illustrates the combined influence of the span of the flap (ranging from 10% to 50% of c), a maximum wing pitch angle (θo) ranging from 75° to 105°, and a maximum flap pitch angle (ψo) varying from 30° to 55° on the resulting power output ( ̅ ).When the flap's extent is at 10% of c, ̅ is maximized at θo = 75° across various ψo, as shown in Figure 14a.The maximum ̅ for varied θo occurs at ψo = 40°, although the variation in power output across different θo is small.As the span of the flap increases to 20% of c, the optimum θo significantly rises to 90°, and the maximum ̅ values are notably higher compared to the flap span of 10% of c, as depicted in Figure 14b.The optimal θo is found to be the same for flap lengths of 25% and 33% of c, as shown in Figure 14c,d, respectively.However, the maximum power generation is attained for a flap span of 40% of c, beyond which it starts to decrease for various θo and ψo values depicted in Figure 14e,f, respectively.
The optimum power output was achieved for the flap spans ranging from 40% to 50% of c, θo = 95°-100°, and ψo = 45° and 50°.As the length of the flap increases, the optimum θo gradually increases from 75° to 100°, whereas ψo varies between 45° and 50°.The optimal power generation was achieved with a flap span of 40% of c at θo = 95° and ψo = 45°.This configuration resulted in a 29.9% enhancement in power output and a 20.2% boost in efficiency relative to the wing configuration without a leading-edge flap.

Conclusions
We investigated a flat plate oscillating wing with an integrated leading-edge flap by focusing on its impact on energy-harvesting performance.Numerical simulations were conducted using the overset grid method under transient conditions, and the k-ω SST turbulence model was employed.Various parameters such as the flap length (from 10% to 50% of the chord length) and maximum pitch angles for the wing and flap (from 75° to 105° and from 30° to 55°, respectively) were quantitatively assessed to identify the configurations that maximize power generation and efficiency.
The results showed that the addition of the leading-edge flap leads to smoother fluid flow along the wing's lower surface and a high-velocity region compared to the wing without a flap.In addition, an increase in the flap's maximum pitch angle or the flap length enhances the wing's effective projected length in the heaving direction.By analyzing the combined effects of varying flap length, wing angle, and flap angle, we identified the optimal configurations that yielded substantial power output enhancements.In particular, flap lengths within the range of 40%-50% of the chord length, combined with θo values from 95° to 100° and ψo of 45°-50°, deliver peak power generation.The results demonstrate a remarkable 29.9% increase in overall power output and a commendable 23% efficiency enhancement compared to the baseline wing configuration without the flap.

Figure 1 .
Figure 1.Configurations of an oscillating wing with and without a flap: (a) 2D view; (b) 3D view.

Figure 2 .
Figure 2. Kinematics of (a) oscillating wing without a flap; (b) oscillating wing with a leading-edge flap.

Figure 3 .
Figure 3. Illustrations of the computational domain and boundary conditions: (a) computational domain; (b) sub-region grid.

Figure 4 .
Figure 4. Comparative analysis of the numerical results obtained in this study against the validation data: (a) comparison of the instantaneous Cpt values for a rectangular-plate wing at Re = 1100 in comparison to the findings reported by Usoh et al. [9]; (b) comparison of the instantaneous CL and CD values for NACA0015 at Re = 5 × 10 5 , benchmarked against the findings reported by Kinsey et al. [48].

Figure 5 .
Figure 5. Comparative analyses of Cy, Cpy, Cpm, and Cpt values within a cycle for the wing with a flap length of 10% of c at θo = 75° and varied ψo (a) Cy (b) Cpy (c) Cpm (d) Cpt.

Figure 6 .
Figure 6.(a) Streamline and vorticity pa erns, (b) velocity magnitude, and (c) pressure distribution of the wing with varied flap lengths at θo = 75° with varied ψo at 0.15 t/T.

Figure 7 .
Figure 7. Pressure coefficient across the wing surface during heaving motion at θο = 75° with varied ψo at 0.15 t/T.

Figure 9 Figure 9 .
Figure9presents the streamline and vorticity, along with pressure contour plots, around the foil's surface for different flap lengths at specific instances of 0.25t/T.In Figure9a, a flap span of 0% of c (equivalent to a flat plate) leads to fluid flow separation below the lower surface of the wing's leading edge.This flow separation decreases as the flap length increases to 10% of c, and it decreases further when the flap length reaches 15% of c.In Figure9b, the associated pressure contour plots reveal that larger flap lengths result in an expanded region of lower pressures on the lower surface.
Figure 11a,b illustrate the Cy

Figure 11 .
Figure 11.Comparison of Cy, Cpy, Cpm, and Cpt within a cycle for the wing with a leading-edge flap of 20% of c, varied θo and ψo = 40° (a) Cy (b) Cpy (c) Cpm (d) Cpt.

Figure 12 .
Figure 12.(a) Streamline and velocity magnitude and (b) pressure contour plots of the wing with a leading-edge flap of 20% of c, varied θo, and ψo = 40° at 0.25 t/T.

Figure 13 .
Figure 13.Pressure coefficient on the wing surface during heaving motion for a flap length of 20% of c, ψo = 40°, and varied θo at 0.25 t/T.

Figure 14 .
Figure 14.Average total power coefficient ( ̅ ) for an oscillating wing with variations in flap length, θo, and Ψo.

Author Contributions:
Conceptualization, M.A. and C.H.S.; Methodology, M.A. and C.H.S.; Simulations, M.A.; Investigation, M.A. and C.H.S.; Writing and editing, M.A. and C.H.S.All authors have read and agreed to the published version of the manuscript.Funding: This research received funding from the Korean government's Ministry of Science and ICT (MSIT), provided through a grant from the National Research Foundation of Korea (NRF), under grant number 2022R1F1A1061903.Institutional Review Board Statement: Not applicable.Informed Consent Statement: Not applicable.

Table 1 .
Physical model parameters of the oscillating wing energy harvester.

Table 2 .
Analysis of mesh independence in the flat plate case.

Table 3 .
Total power extraction coefficients for the oscillating wing at θo = 75°, associated with different ψo values, and the percentage increase in total power output compared to the wing without a flap.

Table 4 .
Overall power coefficients for a wing with different flap lengths at θο = 75° and ψo = 40°.