Power Extraction Performance by a Hybrid Non-Sinusoidal Pitching Motion of an Oscillating Energy Harvester

: This study proposes a hybrid pitching motion for oscillating flat plates aimed at augmenting the energy extraction efficiency of an energy harvester. The proposed hybrid pitching motion, within the first half cycle, integrates a non − sinusoidal movement starting at t/T = 0 and progressing to t/T = 0.25, with a sinusoidal movement initiating after t/T > 0.25 and continuing to t/T = 0.5. The second half of the cycle is symmetric to the first half but in the opposite direction. The calculated results show that the proposed hybrid pitching motion outperforms both the sinusoidal and the non − sinusoidal motions. The hybrid pitching motion merges the merits of both the sinusoidal and non − sinusoidal motions to optimize pitch angle variation. This integration is pivotal for enhancing the overall power output performance of an oscillating energy harvester characterized by momentum change that enhances the orientation of the heaving movement, smoother motion transitions, and consistent energy harvesting. The power generation is obtained at wing pitch angles of 55°, 65°, 70°, 75°, and 80° during a hybrid pitching motion. The proposed hybrid pitching motion, set at a pitch angle of 70°, achieves a maximum power output that exceeds the oscillating flat plate using a sinusoidal pitching motion by 16.0% at the same angle.


Introduction
The global emphasis on renewable and pollution−free energy sources, such as solar, wind, and hydropower, is essential for environmental conservation, enhancing human life quality, and promoting the renewable energy sector [1,2].This transition away from traditional fossil fuels, such as natural gas, coal, and oil, toward renewable energy sources has catalyzed the development of numerous innovative technologies, among which flapping foil technology has emerged as a notable solution due to its eco−efficiency in harnessing energy from natural water and air currents.Although the exploration of hydro/aerodynamics in energy extraction through oscillating devices is relatively young, the foundational concept dates back to 1981 with the work by [3].Their research focused on harmonically oscillating wings and combining pitch and plunge motions to extract wind energy, achieving efficiencies comparable to those of conventional windmills.In recent times, there has been a surge in interest regarding the energy−harvesting capabilities of oscillating foils.
Conducted experimental and numerical studies on a NACA0012 foil using an unsteady panel method revealed that optimal plunging amplitude and frequency could yield total energy extraction efficiencies reaching up to 30% [4].Following this, a comprehensive parametric analysis revealed that enhancing pitching amplitude under constant plunging conditions could transition the system from consuming to generating power [5].Research on a NACA0015 airfoil, utilizing FLUENT, demonstrated a peak efficiency of 34% within certain frequency and pitching amplitude ranges by adjusting the pitching motion, proving it to be more effective than plunging alone in improving energy extraction efficiency [6].Recent advancements in piezoelectric energy harvesting have led to more efficient designs, such as the U−shaped harvester [7] and the semi−submersible model [8].Another study has shown that the efficiency of an oscillating energy harvester is influenced by foil kinematics, system geometry, flow physics, and motion trajectory [9,10].
Energy harvesting from oscillating wing arrangements can be managed through two primary flow management strategies, passive and active management, depending on the need for extra energy within the flow field.Passive management enhances energy harvesting by incorporating simple modifications onto the airfoil surface without the need for additional energy input, offering a cost−effective and easily implementable solution.The flexible deformation enhances vortex movement on the suction side of the foil, boosting heaving force and energy efficiency by modifying pressure differentials around the foil [11].Incorporating a Gurney flap increases the creation of vortices at the trailing edge and amplifies the pressure difference across the wing surface, enhancing lift force and leading to a substantial 21% enhancement in energy harvester [12].A fully passive oscillating wing in free surface flow was explored, aiming to identify ideal submersion levels and the effects of single-frequency waves [13].A feedback control system was proposed to regulate the angle of attack, aiming to boost lift force and thereby enhance power output [14].A shroud enhances oscillating foil power extraction efficiency by 35.8% over an unshrouded counterpart [15].A study revealed that inflexible plates enhance performance at low Reynolds numbers while flexible plates maintain high efficiency at higher Reynolds numbers [16].
Conversely, active management techniques, like flow regulation [17] and plasma devices [18], excel in real−time adjustments, surpassing passive techniques.The use of a fully deformable flapping foil showed a notable 2% increase in power output over the conventional foil setup [19].Flaps have been used to increase the energy performance of oscillating energy harvesters.Energy extraction from an oscillating harvester featuring a trailing edge length of 20% of the chord revealed an 11% power increase compared to an oscillating wing without a flap by optimizing flap and wing angles to 35° and 70° [20].A study conducted to enhance energy harvesting with an oscillating flat plate and trailing−flap found that enlarging the trailing flap on an oscillating flat plate enhances both power output and efficiency, with an increase of 26.9% in output power and 21% in efficiency occurring at a flap length of 60% [21].The authors studied integrating a leading−edge flap in oscillating wing harvester and revealed configurations that boost power by 29.9% and efficiency by 23% [22].
The connection between the efficiency of flapping foil energy harvesters and wake stability has been emphasized by [23].The energy performance of oscillating wing harvesters has been investigated in tandem arrangements by many authors.An inter−foil spacing within a tandem hydrofoil configuration has been found to provide stable energy extraction for the upstream foil, regardless of spacing, whereas the downstream foil's energy extraction showed periodic variations with increasing separation, despite a constant motion phase shift [24].A numerical analysis showed that a novel tandem airfoil configuration oscillating within a convergent duct enhances power extraction efficiency due to the effect of wall confinement [25].In another study, adjusting the downstream spacing in tandem/parallel coupled oscillating wings can enhance energy extraction, resulting in a 23% increase in efficiency by adjusting the vertical inter−foil spacing.This highlights the importance of spatial optimization [26].A concept of a tandem hydrofoil−based tidal array has been introduced, characterized by enhanced density and cost−effectiveness [27].The principles of stability have been investigated, particularly in their relevance for designing and analyzing such configurations [28][29][30][31].
Hydraulic energy holds significant appeal due to its vast potential for reliable power delivery.In advancing sustainable energy solutions for unmanned marine devices, significant progress has been made in efficiently harnessing ultra−low−frequency pitch vibrations, achieving an energy extraction efficiency of 33.23% [32].The effects of angle of attack variations on flapping foil propulsion efficiency in aquatic environments have been investigated, offering a detailed analysis of the hydrodynamic mechanisms that enhance our understanding of biomimetic propulsion technologies [33].Research on optimal propulsion profiles for Micro Air Vehicles and Autonomous Underwater Vehicles has underscored the correlation between angle of attack variations and propulsion efficiency [34].Additionally, the analysis of oscillating foil dynamics in fluids has shown how adjustments to the angle of attack affect propulsion efficiency, which is crucial for the design and optimization of propulsion systems for Micro Air Vehicles and Autonomous Underwater Vehicles [35].Various deflector designs have been explored to enhance power output in flapping foil energy harvesters, focusing on the angle of attack and the distance of upstream deflectors to the hydrofoil, highlighting their impact on power generation efficiency [36].The influence of stratified density flow on energy output was analyzed, particularly in relation to pitching movement amplitude [37].
Energy extraction performance can be significantly improved with an appropriate combination of pitching and heaving motions.The energy extraction performance has been investigated by modifying foil trajectories and positions, emphasizing favorable wake interactions for efficiency improvement [38].The performance of power extraction by an oscillating hydrofoil was evaluated across three distinct movements at a Reynolds number of 1.7 × 10 6 : left swing, right swing, and linear.It was determined that the leftswing motion yielded superior power extraction compared to both the right−swing and linear motions [39].The non−sinusoidal motion showed maximum output power over the sinusoidal motion [40].The non-sinusoidal pitching and plunging motions have been examined to show how flapping foil energy extraction markedly improved [41].Research on a sinusoidal motion found optimal efficiency for a flapping foil at a pitching angle θ o = 75° and a reduced frequency  * = 0.16, surpassing 33% energy efficiency [42].The influence of the non-sinusoidal motion on efficiency enhancement has been investigated, finding θ o = 75° and  * = 0.16 to be optimal, yielding 32% efficiency [43].Semi−active flapping foils with cosinusoidal pitching motions enhanced energy harvesting efficiency [44].
In exploring the practical application of the proposed hybrid pitching motion, this study draws parallels with findings from related aerodynamic studies, such as those by [45], who illustrate how variable amplitude and phase angles can enhance aerodynamic efficiency, which is central to the design principles of this study.Similarly, strategic structural modifications have been shown to amplify energy harvesting outputs, as evidenced by the empirical findings of [46].The methodology of this study is complemented by [47], which demonstrates how specific motion parameters affect fluid dynamics and energy extraction efficiency.Insights into the aerodynamic behaviors at low Reynolds numbers, crucial for refining our energy harvesting approaches, are provided by [48].Exploration of the enhancement of vibrational energies through flow modifications is provided by the findings of [49].Collectively, these studies corroborate the viability of our hybrid pitching motion, affirming that our theoretical innovations have a strong foundation in practical, empirically tested aerodynamics.
The aim is to investigate the potential of using a proposed hybrid pitching motion combining non−sinusoidal and sinusoidal pitching motions with heaving motion in oscillating flat plates to enhance energy harvesting performance.A quantitative assessment was undertaken to examine the impact of various pitch angles 55°, 65°, 70°, 75°, and 80°.This comprehensive investigation provided insights into the importance of modifying flat plate trajectories for improving the performance of energy harvesters.

Numerical Methodology
The simple geometry of a rectangular flat plate emphasizes its superior performance over the profiled airfoil [50].It also provides maximum power and efficiency conditions for an oscillating wing and has the following key values: the chord length of the wing, c = 1; the location of the pivot, X p = 0.3333; the amplitude of the heaving motion h o /c = 1; the maximum of the pitch angle θ o = 75°; the difference of the ∅ = 90°; and a reduced frequency of  * = 0.14.
Figure 1 shows a schematic of a 2D flat plate with a foil chord span of c, thickness and   4% of the chord length.The proposed hybrid pitching motion investigates energy extraction efficiency at the same conditions of reduced frequencies  * = 0.14, pitch angle θ o = 75°, X p = c/3, and h o /c = 1.0 plunging amplitude that provides the optimal power and efficiency as mentioned in [50].Figure 2 shows the variation of instantaneous pitching angle θ(t) in one period with sinusoidal pitching profiles, as depicted in Figure 2a. Figure 2b,c show the non−sinusoidal and the proposed hybrid pitching motions as a function of the parameter β, which systematically transition from β = 1.0, 1.5, 2.0, 2.5, 3.0, and 4.0 with a 90° phase lag the heaving motion.The sinusoidal, the non−sinusoidal, and the proposed hybrid pitching motions are expressed as Equations ( 1), (2), and (3), respectively; ω = 2π is angular oscillating frequency, and T = 2π/ω is time.
Figure 2a and Equation ( 1) show the profiles of the instantaneous pitching amplitude θ(t) for the sinusoidal motion [50], where the pitch angle increases and decreases in a regular pattern over time.The flat plate reaches its highest point on the graph at (t/T = 0.25), representing the maximum positive pitch angle (+ θ o ) in the first half of the cycle.Subsequently, the flat plate reaches the horizontal position at (t/T = 0.5), while the lowest point represents the maximum negative pitch angle (− θ o ) at (t/T = 0.75) in the second half of the cycle.
Figure 2b and Equation (2) show the characteristic behavior of the instantaneous pitching amplitude, θ(t), for non−sinusoidal motion that is engineered by integrating a specially designed trapezoidal−like pitching profile modulated by an adaptable parameter β, as described in [51].Notably, when the modulation parameter β = 1.0, the resultant pitching motion closely approximates a sinusoidal profile with a phase lag of 90º behind the heaving movement.As the amplitude θ o remains constant, deviations of β from 1.0 result in an extended duration, where θ(t) sustains the peak values of ± θ o .This effect becomes increasingly pronounced as β ascends toward infinity (β → ∞), with the profile tending toward a square wave.
Figure 2c and Equation ( 3) illustrate the proposed hybrid pitching motion, ingeniously synthesizing the attributes of both the non−sinusoidal and sinusoidal motions delineated previously in Equations ( 2) and (1), respectively.Within this framework, the flat plate exhibits an early swift ascent to the maximum pitch angle via the non−sinusoidal trajectory in the initial quarter of the cycle (0 ≤ t/T ≤ 0.25).Thereafter, from the midpoint of the ascent to the return to the horizontal position (0.25 < t/T ≤ 0.5), the motion adheres to a sinusoidal path.Figure 2c displays the progression of the pitching motion, θ(t), across a range of β values, marking a shift from an almost sinusoidal motion when (β = 1.0) to a more pronounced, half-square waveform as the value of β increases.This transition is particularly notable at the t = 0.25T mark, where the maximum pitch angle is observed.The parameter β acts as a key factor for non-sinusoidal and hybrid motions influencing the time intervals at which the pitch amplitude θ o reaches its extreme values.For instance, in a sinusoidal pitching motion, the maximum pitch angle occurs at t/T = 0.25 of time.However, with the proposed hybrid motion (β = 1.5), the maximum pitch angle reaches slightly earlier, at t/T = 0.23, indicative of an enhanced power generation performance.
Figure 3 shows the time variations of the pitching profile θ(t) for an oscillating flat plate.The pivot point is 0.33c from the leading edge, in line with the centerline of the flat plate.The upstream flow velocity is represented as U ∞ .Due to the rapid pitching motion at the non−sinusoidal and proposed hybrid motions, the flat plate attains higher angles compared to the sinusoidal pitching motion.The flat plate undergoes heaving motions, represented by h(t) in Equation ( 4), with h and h 0 denoting the immediate and maximum vertical oscillations of the wing, respectively.

Re
where  is the fluid density, c is the chord length,  ∞ is the free-stream velocity, and μ is the dynamic viscosity.The two-dimensional unsteady numerical simulation was conducted at Reynolds numbers   500,000 o analyze the performance of the energy harvester.
where * = c/U ∞ is reduced frequency.
The instantaneous power (P) is integrated over a cycle to determine the average power (), which is used to measure the energy production.
here, Y(t) denotes the force in the vertical direction, M(t) represents the pitching moment, and T stands for the period of oscillation.
here, C L (t) represents the instantaneous pushing coefficient and C M (t) indicates the instantaneous momentum coefficient, which are determined as follows: The average power coefficient C pt over a cycle is calculated by the integration of instantaneous C pt over a cycle. or Equations ( 10) and ( 14) reveal that C pl and C pl represent the pushing-based energy extraction, while C pm and C pm signify the moment-based contribution.
The calculation of the power extraction efficiency η can be performed utilizing Equation (15).
where d is the greatest displacement in the y direction of a flapping plate.
The overset mesh method was utilized, which is well−suited for modeling massive body movements, intricate forms, and multi−component structures [52].To simulate flow over the oscillating wing, the unsteady incompressible Reynolds-averaged Navier-Stokes (RANS) flow equations were utilized in Fluent software version 21 R1 [53].The simulation employed the finite volume method and a pressure−based solver.The governing equations are defined by Equations ( 16) and (17).
and   represent the velocities.In Equation ( 17), the Reynolds-stress term    ′   ′ ̅̅̅̅̅̅ is computed through the k−ω SST turbulence model by previous investigators in comparable investigations [52,54].Spatial discretization for pressure, momentum, and turbulent viscosity resolution is accomplished using the second−order upwind method.Flat plate motion was input through defined functions of the user (UDF), and grid movement was executed using the method of dynamic mesh.The y+ value was set to meet the requirements of the turbulence model used, ensuring precise resolution of the boundary layer.In this study, a y+ value of 5.0 yielded a total output power coefficient of 1.04.This closely aligns with the results obtained using a y+ value of 0.98, where the total power coefficient was 1.05.Figure 4 presents the rectangular computational region, meshing structure, subdivided mesh, and external constraints utilized for simulating this case study.An investigation into domain independence was conducted with the length, l, and the width, b, of the numerical domain.The range for  is from −20c to 50c, and for b, from −25c to 25c was chosen, as further increasing the simulation area had no influence on the results.The pressure for the exit point of the boundary was positioned at 50c, while the top and bottom boundaries were subjected to a slip boundary condition, as depicted in Figure 4a.The wing's main body (highlighted by the rectangular region in blue) in Figure 4b was associated with moving grids and the detailed view of mesh near wing wall corners.The moving grid region took the form of a rectangular grid with a length of 6c and a width of 5c.The model utilized 2000 time steps for each oscillation period, with the outcomes being charted following the completion of seven cycles.Mesh independences were carried out with three different groups of structured meshes: coarse, medium, and fine mesh.The effect of grid resolution and the time step on the calculated C ̅ Pt value is presented in Tables 1 and 2. For grid independence tests, three sets of structured grids were utilized within the overset mesh framework.Firstly, for the moving grid, a coarse grid consisting of 0.6 × 10 5 elements, a medium grid with 1.2 × 10 5 elements, and a fine grid comprising 2.6 × 10 5 elements were employed.Secondly, for the background, a coarse grid of 0.3 × 10 5 elements, a medium grid of 0.6 × 10 5 elements, and a fine grid of 1.2 × 10 5 elements were used.For the time step independence study, the medium mesh was employed with time intervals of 500, 2000, and 4000 for each cycle.Table 1 and Table 2 show the power output coefficient and relative variation of power output coefficient over time for different mesh sizes and time steps, respectively.The medium−mesh resolution satisfies the accuracy in space and nearly agrees among different time steps.In validating the numerical model and mesh, a comparison with the results from [50] and [20] was conducted, revealing a favorable agreement, as depicted in Figure 5. Based on these results and to balance between precision and processing time, we opted for a medium of 1.

Output Power and Efficiency
This study investigates the performance of output power from flapping energy harvesters through the application of the proposed hybrid pitching motion, varying the parameter β and pitching angles under  * = 0.14 and h o /c = 1.0.Table 3 summarizes the C pt , C pl , and C pm over one cycle for quantitative comparison.In Table 3, the proposed hybrid motion (β = 1.5) shows the maximum extraction power, C pt = 1.15, and efficiency = 40.78%,respectively, at the optimum pitch angle of θ 0 = 75° [50],  * = 0.14, and h o /c = 1.0.
Table 3 shows the analysis of efficiency across different pitching motions and values of β for a pitch angle of 75º, revealing significant variations in performance.The sinusoidal motion establishes a baseline efficiency of 40.31%.In comparison, the non−sinusoidal and hybrid motions demonstrate varying efficiencies due to an increase in the maximum vertical displacement, d, of the flapping flat plate, which directly impacts the efficiency.As β is 1.5, the hybrid motion of 40.78% is slightly higher than the sinusoidal and non−sinusoidal pitching motions.
Figure 6 shows the variation of C pl , C pm , C pt , and η for the sinusoidal, non−sinusoidal, and hybrid pitching motions by changing the parameter β.The plot shows the hybrid motion surpassing the sinusoidal and non−sinusoidal pitching motions in total power extraction C pt , pushing power, C pl , and efficiency, β.In contrast, the moment power, C pm , in the non−sinusoidal motion is higher than the sinusoidal and hybrid motions.At β = 1.0, the C pt and η of the non−sinusoidal motion is slightly higher than those of the hybrid motion because the non−sinusoidal motion exhibits characteristics at β = 1.0 similar to those of the sinusoidal motion, as shown in Figure 2.However, excessively increasing β could reduce performance, such as β = 3.0 and 4.0.This could be seen in the calculated results that reveal that the efficiency of output power is defined by the contributions from the different combinations of C pl , C pm , and corresponding heaving and pitching velocities, all of which are affected by β values.For a specific nominal pitch angle, there is an optimal value of β at which the power output is maximized.The C pt and η values for the sinusoidal, non−sinusoidal, and hybrid pitching motions, evaluated at various pitch angles ranging from 55° to 80° with a constant reduced frequency  * = 0.14 and a heave amplitude to chord ratio h o /c = 1.0, are summarized and listed in Table 4.Although the full range of β values from 1.0 to 4.0 has been calculated for each pitch angle, only the C pt values resulting from the optimum β values are shown in this table.Table 4 and Figure 7 illustrate that hybrid pitching motions consistently deliver higher power coefficient values C pt compared to their sinusoidal and non−sinusoidal counterparts across identical pitch angles.The incorporation of hybrid motion attributes can substantially boost power generation specifically; at a 55° pitch angle, the hybrid motion (β = 3.0) achieves the highest relative increment of 68.6% in C pt compared to the sinusoidal motion.The output power of the hybrid motion at a pitch angle of 70º demonstrates a 16% enhancement when compared to the sinusoidal motion with β = 2.0.In the case of the optimal pitch angle for a sinusoidal motion of 75° [50], the output power of the hybrid motion shows a 10.57% enhancement.At the pitch angle of 80°, the advantage of the hybrid motion is also evident, where it exhibits an 11.10% increment over the sinusoidal motion, despite both having a β value of 1.0.The trend across the data indicates that increasing the pitch angles in the hybrid pitching motions could enhance the power output of the system.Table 4 shows the C pt values for each type of motion over a range of pitch angles.
Figure 7b shows the efficiency for the sinusoidal, non−sinusoidal, and hybrid pitching motions across varying pitch angles from 55° to 80°.The hybrid motions, combining attributes of both sinusoidal and non−sinusoidal dynamics, consistently achieve higher efficiencies throughout the range, peaking at 75°.The graph clearly demonstrates that the hybrid motions are consistently outperforming the other two types.

Mechanism of Energy Output Enhancement
Figure 8 shows the calculated results of C L , C M , C pl , C pm , and C pt for the sinusoidal motion and the function of β for the hybrid pitching motions.Figure 8a shows an enhancement in the pushing force coefficient C L across the oscillation cycle.C L contributes positively to the pushing power coefficient for every case of the proposed hybrid pitching motion.
Within the time interval (0.0 ≤ t/T ≤ 0.25), the rapid attainment of maximum pitch angles results in a heightened pushing power C pl compared to the sinusoidal motion.As depicted in Figure 8a, following the peaks, the pushing forces coefficient of the hybrid motions drops to lower values compared to the sinusoidal pitching motion.The C L values for the hybrid pitching motions surpass those of the sinusoidal pitching motion until t/T = 0.015, indicating that C pl increases with β.This can be attributed to the fact that as β increases, the flat plate reaches its maximum pitching angle more swiftly compared to the gradual approach of the sinusoidal motion, thereby resulting in an increase in C pl , as shown in Figure 8c.
Figure 8b,d illustrate the pitching moment coefficient, C M (t), and the moment power coefficient, C pm , with distinct peaks and valleys that align with specific phases of the oscillation cycle.Figure 8d shows that the moment power coefficient, C pm , could reach its minimum value by increasing parameter β through the process loop.
Figure 8e presents the variations of C pt , exhibiting its dependency on both pushing power and moment power coefficients across different pitching motions.At specific times t/T = 0.25 and t/T = 0.75, the coefficients related to the pushing force, pushing power, and total power show nearly identical values for both the sinusoidal and hybrid pitching motions.This result, as shown in Figure 8a,c,e, arises because the pitch angles are the same at these specific times for all pitching motion cases.Figure 10 shows the static pressure for three different pitching motions, sinusoidal, non−sinusoidal (β = 1.5), and hybrid (β = 1.5), at  * = 0.14 and θ 0 = 75° in a flapping flat plate.It is evident that the pressure on the windward side of the plate exceeds that on the leeward side when the flat plate is descending between t/T = 0.0 and 0.50, thereby contributing favorably to pushing force for extracting flow energy.The pressure contours of the non−sinusoidal and hybrid cases look similar during 0.0 ≤ t/T ≤ 0.25, and the pressure contours of the sinusoidal and hybrid cases look similar during 0.25 < t/T ≤ 0.50.The detailed surface pressure distribution is plotted in Figure 11.The pressure distribution around the flat plate directly influences the enhancement of aerodynamic forces, such as pushing and moment, which are directly linked to energy extraction performance.Figure 11a shows that the hybrid and non−sinusoidal motions are distinguished by their higher pressure difference than sinusoidal differences at t/T = 0.05. Figure 11b,c show the calculated surface pressures at t/T = 0.10, and the pitch angles for both the hybrid and non−sinusoidal motions exceed those of the sinusoidal motion.In these cases, the increased pitch angles lead to a greater negative pressure on the leeward side of the plate.In the hybrid and non−sinusoidal motions, the projected length along the x−axis is shorter compared to the sinusoidal motion.As illustrated in Figure 11d,e for t/T = 0.20 and t/T = 0.25, respectively, the pitch angles of the sinusoidal, non−sinusoidal, and hybrid motions reach similar values, resulting in similar pressure distributions.This result is also reflected in the pushing force coefficient, pushing power coefficient, and total power output, as depicted in Figure 8a,c,e.The non−sinusoidal motion maintains a higher pitch angle during this time period (0.25 < t/T ≤ 0.50) compared to both the sinusoidal and hybrid motions.Because of the higher pitching angle, it exhibits a shorter projection length along the x−axis than the other cases.This short projection length results in a reduced pushing force, subsequently lowering the power output of the non−sinusoidal motion.This relationship is shown in Figure 11f-i

Conclusions
This study introduced and investigated the proposed hybrid pitching motion for oscillating flat plates within an energy harvester framework, aimed at enhancing the efficiency of power extraction.The hybrid motion effectively combines the advantageous aspects and characteristics of both the sinusoidal and non−sinusoidal motions, yielding an optimized pitch angle variation that improves the performance of the harvester.
The numerical analysis demonstrated that the hybrid motion outperforms the traditional sinusoidal and non−sinusoidal motions across various wing pitch angles (55°, 65°, 70°, 75°, and 80°).Notably, the power output was improved by approximately 10.57% over the optimum condition for the oscillating flat plate at a wing pitch angle of 75°, underscoring the utility of this proposed approach.The maximum energy extraction for the proposed hybrid pitching motion was at a pitch angle of 70°, showing a 16.0% increase in power output compared to the same pitch angle condition in the sinusoidal pitching motion.The proposed hybrid motion presents a notable advancement in the field of energy harvesting, offering a viable path for improving the performance of oscillating flat plate harvesters.

Figure 1 .
Figure 1.Illustrations of an oscillating flat plate energy harvester.

Figure 4 .
Figure 4. Simulation space and grid model: (a) simulation space and (b) moving area and the magnified view of the meshed computational domain near the corner of the wing wall.

2 ×Figure 5 .
Figure 5.Comparison of the numerical results and the present study and the validation data: (a) pushing power coefficient for an oscillating plate at Re 1.1 × 10 3 by Usoh et al. [50]; (b) pushing force coefficient at Re 500,000 by Alam and Sohn [20]; (c) pushing power coefficient at Re = 500,000 by Alam and Sohn [20].

Figure 6 .
Figure 6.Comparison of sinusoidal, non−sinusoidal, and the proposed hybrid pitching motions versus various values of β (a) pushing power coefficient; (b) moment power coefficient; (c) output power coefficient; and (d) efficiency at  * = 0.14 and θ O = 75°.

Figure 7 .
Figure 7. Variation of (a) output power coefficient and (b) efficiency for different pitching amplitude and different values of β at h o /c = 1.0 and  * = 0.14.

Figure 8 .
Figure 8. Variations of (a) pushing force coefficients C L (t); (b) pitching moment coefficients C M (t); (c) instantaneous pushing power coefficient C pl ; (d) instantaneous moment power coefficient C pm ; and (e) instantaneous total power coefficient C pt for sinusoidal and hybrid motions with different values of β at  * = 0.14 and θ O = 75°.

Figure 9
Figure9shows the streamline and vorticity contours in half of the cycle for sinusoidal, non−sinusoidal (β = 1.5), and hybrid (β = 1.5) at θ 0 = 75°.Subsequent to the downward pitch, a leading−edge vortex (LEV) develops on the windward side.This vortex then disappears by the stagnation flow on the upper surface and generates the LEV on the leeward lower surface near t/T = 0.15.After this point and throughout the descent,

Table 2 .
Grid time step independence analysis.

Table 3 .
The instantaneous C pt , C pl , and C pm for sinusoidal motion, non-sinusoidal motion, and hybrid pitching motions at different values of β for θ O = 75° and  * = 0.14.

Table 4 .
Maximum output power coefficient for sinusoidal, non−sinusoidal, and hybrid pitching motions with different pitching angles at  * = 0.14 and h o /c = 1.0.