Next Article in Journal
Diverse Decoding for Abstractive Document Summarization
Previous Article in Journal
Study on Optical Properties and Internal Quantum Efficiency Measurement of GaN-based Green LEDs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Study on the Effect of Non-Sinusoidal Motion on the Energy Extraction Performance of Parallel Foils

Shaanxi Engineering Laboratory of Turbomachinery and Power Equipment, School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an, Shaanxi Province 710049, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(3), 384; https://doi.org/10.3390/app9030384
Submission received: 2 January 2019 / Revised: 17 January 2019 / Accepted: 20 January 2019 / Published: 23 January 2019
(This article belongs to the Section Energy Science and Technology)

Abstract

:
The effect of non-sinusoidal motion which influences the energy extraction performance of foil is considered in this paper. Two oscillation motions, the combined non-sinusoidal plunging and sinusoidal pitching motion, as well as the combined non-sinusoidal pitching and sinusoidal plunging motion, are selected to investigate the oscillation process of two-dimensional parallel foils numerically. The optimal oscillation motion and average power coefficient at different combined motions are gained. The effects of the plunging motion and pitching motion at different oscillation motions are analyzed, and the evolution law of the foil lift force and vortex field are obtained. It is indicated that the non-sinusoidal motion has a significant influence on energy extraction. When the motion is combined (non-sinusoidal plunging and sinusoidal pitching motion), the best extraction performance is gained at Kh = −0.5. The maximal CPm is 0.375 and the maximal η is 0.188. When the motion is combined (non-sinusoidal pitching and sinusoidal plunging motion), the maximal CPm is 0.623 and the maximal η is 0.312 which appear at Kθ = 2. For the same frequency, the more the plunging motion is similar to the sinusoidal motion, the more energy is extracted by foils. While the more the pitching motion approximates to the square wave, the worse the achieved extraction performance is.

1. Introduction

As a new energy extraction device using clean energy, foil is thought to have a bright prospect, especially in the current situation of energy shortages and environmental pollution. Since the combined plunging and pitching motion was put forward and realized to extract energy in 1981 [1], plenty of researchers have discussed the energy extraction characteristics of foil by numerical and experimental methods. According to the laminar computation model, Campobasso et al. [2] simulated the energy extraction performance of foil at a low-velocity flow field with Re = 1100. It was found that the energy extraction efficiency can achieve 35% and the deep stall of airfoil has a positive effect on efficiency. Abiru [3] discussed the energy extraction performance and the effect factors (such as the plunging and pitching amplitude, frequency and the free stream velocity) by experiments. The results showed that the power coefficient and efficiency increase with the increase in the amplitude. The experiment system efficiency is observed to be at least 32%. A wind energy extraction device was designed by Shimizu et al. [4], which is based on the flutter mechanism of foil. The device worked under a low frequency and gained a high extraction performance and efficiency. Bryant et al. [5] established a quasi-stable state model and a semi-empirical model of dynamic stall, optimizing the aerodynamic force and energy extraction performance of foil. They observed that the semi-empirical model of the dynamic stall is effective in predicting the aerodynamic force and efficiency. The numerical research on the extraction performance of foil discusses the potential effect factors like the oscillation motion, the motion parameter, the Re number and the profile of foil by the control variable method and the traditional optimization algorithm, aiming at finding a higher power output and extraction efficiency [6,7,8,9,10,11,12,13]. The oscillation process always adopts the sinusoidal motion, however, some researchers tried to use non-sinusoidal motion on foil for further investigations. Platzer et al. [14,15,16] designed a foil device extracting energy in non-sinusoidal motion and conducted numerical simulations on the power output, the extraction efficiency and the feasibility of this device. Xiao et al. [17,18] carried out an analysis of different non-sinusoidal pitching motions which are similar to trapezoidal ones, their research revealed that foil under non-sinusoidal motions like trapezoid motion has a higher power output and extraction efficiency at some St number. Ashraf [19] discussed the effect of different non-sinusoidal motions on the flow field around the foil. They found that the power coefficient of non-sinusoidal motions improved by 17% compared to sinusoidal motion.
Recently, research about the effect factors of foil extraction performance has been becoming more and more perfect and complete. However, most research focused on single foil and only a few of the researchers analyzed the extraction characteristics of foils under parallel configuration. According to our previous research [20], the flow field around two parallel foils will affect each other during the oscillation motions, which results in a better performance than that around the single foil. Therefore, it is necessary to further explore the influence of non-sinusoidal motion to the energy extraction performance of parallel foils.
In this paper, two oscillation motions are considered, which are the combined non-sinusoidal plunging and sinusoidal pitching motion, as well as the combined non-sinusoidal pitching and sinusoidal plunging motion. The numerical analysis of energy extraction characteristics is conducted. The rest of this paper is introduced as follows. Firstly, the computation model, the numerical method and the oscillation control equations of parallel foils are described. Then, the energy extraction characteristics of parallel foils under two different combined motions are analyzed in detail. Finally, the conclusions about the effect from non-sinusoidal motions are drawn.

2. Numerical Model and Method

2.1. Model of Parallel Foils

The parallel foils considered in this paper are two-dimensional NACA 0012 airfoils with the chord length c = 1 m. As shown in Figure 1, two foils are arranged next to each other along the vertical direction of the chord length, and the space between the two foils is δ. Parallel foils move in the combined plunging and pitching motion, the foil A and foil B (see Figure 1) are the starting positions of foil oscillation. As the oscillation starts, the foils move outwards during t/T = 0–0.5. The transient pitching angle α(t), plunging displacement h(t), pitching amplitude θ0 and plunging amplitude H0 are shown in Figure 1. While the parallel foils move to their plunging limited positions, the foils start the reserve motions. The foils get closer and closer during t/T = 0.5–1, and finally, they get back to the starting position. An oscillation cycle of the foils has thus finished. The pitching axis of the foil is one-third chord length to the leading edge of the foil and δ is selected as 1.2 m throughout the whole numerical analysis.

2.2. Numerical Method

According to the finite volume method, the numerical simulations of the parallel foils oscillating in different non-sinusoidal motions are conducted by the commercial software Fluent 16.0 and the incompressible viscous flow field around foils is analyzed. The spatial term is discretized in the form of the second-order upwind type, while the transient term is discretized by the second-order-accurate backward implicit scheme [21]. All cases are investigated at Re = 1100, and the laminar model is selected to discuss the evolution law of the flow field around parallel foils.

2.3. Design of Computation Domain

In order to save computation resources and time while keeping a high simulation accuracy, the computation domain of the parallel foils is divided into several parts. The region close to the parallel foils is the inner zone, and the region outside the inner zone is the outer zone, as shown in Figure 2. Along the parallel foils’ surface, the no-slip boundary condition is adopted. The interface between the inner and outer zone is set as the interface which connects the structured mesh and unstructured mesh. The velocity inlet is defined as U = 1 m/s and the pressure outlet is set. Based on the C programming language, the plunging and pitching motion characteristics in all simulations are controlled by the user-defined function (UDF) in Fluent. The dynamic mesh technique is applied to realize the reasonable transformation and reconstruction in the inner zone during the oscillation process, the quality of grids around the parallel foils is kept at a proper accuracy.
As displayed in Figure 3, large transformation due to the oscillation of foils occurs in the inner zone grid. Thus, the unstructured mesh is selected for the inner zone, which is denser and where it is easier to improve the grid quality. The structured mesh is adopted at the foil surface to ensure the computation accuracy of the boundary layer near the foil surface. For the outer zone, there is a much smaller transformation during oscillation, so the sparse structured mesh is used here to save the computation time. The mesh contains 7.5 × 104 cells, the height of the first cell in the boundary layer is set as 10−5c to satisfy the computation accuracy.

2.4. Kinematics

The frequency of the parallel foils applied in this paper is the reduced frequency, which is the non-dimensional frequency. It is defined as follows:
k = 2 π f c U
where f represents the frequency, c represents the chord length of foil, U represents the velocity of the free stream.
Based on the research of Xiao et al. [17], the power coefficient shows a similar trend under the different nominal angle of attack α0. α0 is set at 15° in this paper, the corresponding pitching amplitude θ0 can be calculated as follows:
θ 0 = arctan ( 2 π f H 0 c U ) + α 0
where H0 is the dimensionless plunging amplitude and it is set at 1.
Our previous study focused on parallel foils with combined sinusoidal plunging and sinusoidal pitching motion [20]. To discuss the effect of non-sinusoidal motion, the energy extraction characteristics of parallel foils under the combined non-sinusoidal plunging and sinusoidal pitching motion, as well as the combined non-sinusoidal pitching and sinusoidal plunging motion, are both investigated. Sinusoidal pitching oscillation and sinusoidal plunging oscillation are defined as Equations (3) and (4):
α ( t ) = θ 0 sin ( 2 π f t + β )
h ( t ) = H 0 c sin ( 2 π f t + ϕ + β )
where β for the upper foil is set at 180° and β for the lower foil is set at 0°. φ represents the phase difference between the plunging motion and pitching motion; φ is set as 90° [10].
The non-sinusoidal pitching motion is calculated based on the control parameter of pitching motion Kθ:
α ( t ) = { θ 0 sin 1 [ K θ sin ( 2 π f t + β ) ] sin 1 ( K θ ) ,     1 K θ < 0 θ 0 sin ( 2 π f t + β ) ,     K θ = 0 θ 0 tanh [ K θ sin ( 2 π f t + β ) ] tanh K θ ,     K θ > 0
In a similar way, non-sinusoidal plunging motion is calculated based on the control parameter of plunging motion Kh:
h ( t ) = { H 0 c sin 1 [ K h sin ( 2 π f t + ϕ + β ) ] sin 1 ( K h ) ,     1 K h < 0 H 0 c sin ( 2 π f t + ϕ + β ) ,     K h = 0 H 0 c tanh [ K h sin ( 2 π f t + ϕ + β ) ] tanh K h ,     K h > 0
In this paper, different non-sinusoidal plunging motions and non-sinusoidal pitching motions are considered, including the cases under Kh = −1, −0.5, 0, 1 and 2 as well as the cases under Kθ = −1, −0.5, 0, 1 and 2. The curves of the five non-sinusoidal plunging motions and non-sinusoidal pitching motions at k = 0.8 are presented in Figure 4a,b, respectively.
The energy extraction process of the foil over one cycle consists of two parts, the plunging oscillation and the pitching oscillation, the power output P is defined as follows:
P = F y ( t ) V y ( t ) + M ( t ) ω ( t )
where Fy(t) is the lift force in the direction of the y-axis, M(t) is the torque, ω(t) is the velocity of the pitching oscillation.
The transient power coefficient CP can be computed by Equation (8):
C P = P 1 2 ρ U 3 c
The average power coefficient is composed by the pitching motion part and the plunging motion part, the definition is shown as follows:
C P m = C P h m + C P θ m = 1 T U 0 T [ C l ( t ) V y ( t ) + C m ( t ) ω ( t ) ] d t
The lift force coefficient Cl and the torque coefficient Cm are defined as Equations (10) and (11):
C l = F y 1 2 ρ U 2 c
C m = M 1 2 ρ U 2 c 2
The extraction efficiency η is computed as the ratio of the extracted power to the total power of the free stream flow in the swept area, the definition is as Equation (12):
η = P m 1 2 ρ U 3 A = C P m c A
where A represents the area of the foil’s swept region, which is selected as 2H0c [17].

2.5. Validation study

The grid independent and time-independent validation are conducted to ensure that the mesh strategy adopted and the time step set in the simulation have no influence on the analysis accuracy. As shown in Figure 5, three kinds of grids are selected, which are 5.2 × 104 cells, 7.5 × 104 cells and 1.12 × 105 cells, respectively. At the same time, three kinds of time steps are adopted, which are 600, 1200 and 2400 time steps in one oscillation cycle, respectively. The variations of the transient power coefficient CP in a cycle are obtained. The difference between the three kinds of mesh strategies is small enough to ignore. Additionally, the results among the tested time steps show little difference. Thus, the mesh strategy with 7.5 × 104 cells and 1200-time steps per cycle are adopted in all numerical simulations.

3. Results and Discussions

3.1. Combined Non-Sinusoidal Plunging and Sinusoidal Pitching Motion

In this research, the energy extraction performance of parallel foils at combined non-sinusoidal plunging and sinusoidal pitching motion is investigated firstly. Five different plunging motions, which are controlled by the non-sinusoidal plunging parameter Kh (Kh is set as −1, −0.5, 0, 1 and 2 respectively), are simulated. As shown in Figure 6, when Kh is −0.5, 0 and 1, the average power coefficient CPm ascends firstly and then decreases with the increasing k, showing a stable trend in general. At these three motions, all peak values of CPm appear at k = 0.8. Additionally, the maximal CPm is achieved when Kh = −0.5, which is 0.375. When k < 1.0, the best energy extraction performance of foil is observed at Kh = −0.5 under different k. Additionally, the optimal performance is gained at Kh = 1 when k > 1.0.
It is found from Figure 6 that CPm decreases obviously and becomes a negative value gradually with the increase in k when Kh = −1 and Kh = 2. This means these two oscillation motions are in the energy consumption situation and not the energy extraction situation in the range of frequency considered in this paper. It can be concluded from Figure 4 and Figure 6 that, at the same frequency, the more the plunging motion is similar to the sinusoidal motion, the more the energy that can be extracted by parallel foils. At the same time, the extraction efficiency η shows a similar trend with the CPm curves at the non-sinusoidal plunging motion. When k = 0.8 and Kh = −0.5, the maximal extraction efficiency η is observed, which is 0.188.
Parallel foils in this paper oscillate as the combined plunging and pitching motion. Thus, to further investigate the effect factors and change law of energy extraction performance, the average power coefficient CPm is divided into the one caused by the plunging motion CPhm and the one caused by pitching motion CPθm. In general, the plunging motion is the process that extracting energy and the pitching motion is the process that consuming energy. However, as shown in Figure 7, when Kh = −1, the plunging motion converts the energy extraction situation into the energy consumption situation with increasing k. While the CPhm increases gradually with the increasing k at the other four non-sinusoidal plunging motions. It indicates that the plunging motion with Kh = −1 is not appropriate to the energy extraction of the foil. For the pitching motion, they are all sinusoidal motion in different cases. In the figure, the little difference of CPθm at the same frequency is from the difference of the flow field caused by non-sinusoidal plunging motion. With the rise in k, CPθm declines obviously, which means the energy consumption increases. This is because the foil oscillates faster at high frequencies and the angular velocity of the adjusting process of the foil near the plunging limited position is higher, consuming more energy. In summary, the difference of CPm curves in Figure 6 is mainly from the different non-sinusoidal plunging motion.
The amount of energy extracted during the plunging oscillation process is related to the synchronicity between the lift coefficient Cl and the plunging velocity Vy. The Cl curves of the upper foil at Kh = 0 and Kh = 2 with k = 0.8 are shown in Figure 8. By comparing these two cases, it is concluded that the synchronicity of Cl and Vy is better at Kh = 0. While the synchronicity is not well at Kh = 2, even the totally opposite trend of Cl and Vy during t/T = 0.17–0.3 and t/T = 0.67–0.81 is observed. Although the peak value of Cl at Kh = 2 is obviously higher than that at Kh = 0, the CPm at Kh = 2 is relatively low because of the bad synchronicity of Cl and Vy. Thus, the energy extraction characteristics of foil are relatively poor when Kh = 2.

3.2. Combined Non-Sinusoidal Pitching and Sinusoidal Plunging Motion

The effects of the five different non-sinusoidal pitching motions on the extraction performance of parallel foils are considered. The motion of the foil is combined non-sinusoidal pitching and sinusoidal plunging motion. The non-sinusoidal pitching motion is controlled by the parameter Kθ. The energy extraction performance values at different Kθ and k are simulated, as shown in Figure 9. It can be found that the CPm at four kinds of motions ascends firstly and then decreases with increasing k, except the case at Kθ = −1. The peak values of CPm at these four motions all appear at k = 0.8. This means k = 0.8 is the optimal frequency which could allow for the maximal power output in general. When Kθ = 2, the maximal CPm is obtained, which is 0.623.
When the frequency is low, the effects of different oscillation motions on the energy extraction performance are less. With the increase in k, the influence from the non-sinusoidal pitching motion on CPm increases. Foils go into the energy consumption situation earlier and faster in the case with Kθ = −1 and Kθ = 2 at high frequencies. The extraction efficiency η has the same trend with CPm. The optimal extraction efficiency η achieves 0.312 at k = 0.8 with Kθ = 2.
To investigate the influence of plunging motion and pitching motion on energy extraction, we separate CPm into the plunging part CPhm and pitching part CPθm, as shown in Figure 10. In the plunging oscillation process, the energy extracted by the plunging motion under four motions (except Kθ = −1) ascends with the increasing k. This is due to the increasing frequency lifts the plunging velocity Vy of foil, more energy is extracted by plunging motion as a result.
When the frequency is relatively low, the different Kθ has little influence on the power consumption of the pitching motion. Some non-sinusoidal pitching motions obviously cause the energy consumption with the increase in k, like Kθ = 1 and Kθ = 2. This indicates that the more the non-sinusoidal pitching motion is similar to the square wave, the worse the extraction performance that is obtained. This is because the more the pitching motion is like a square wave, the higher the average pitching velocity in an oscillation cycle will be, causing higher energy consumption. That is the reason that causes the sharp decline of CPm at high frequencies in the case of Kθ = 2 (see Figure 9).
Figure 11 shows the Cl curves of the upper foil under two typical non-sinusoidal pitching motions Kθ = 1 and Kθ = 2 at k = 0.8. In these two cases, the plunging oscillation of the foil moves as with the same sinusoidal plunging motion. However, because of the effect from the non-sinusoidal pitching motion, the flow field around the foil is different, resulting in different forces in the y-direction. That is to say, the different pitching motions that cause different Cl values in the plunging motions are the same. This explains the reason for different CPhm at the same k in Figure 10. From Figure 11, it is observed that the synchronicity of Cl and Vy at these two motions is similar. When Kθ is 2, Cl fluctuates obviously near the plunging limited position, the peak value of Cl is higher. At the same time, Cl in an oscillation cycle at Kθ = 2 is higher than that at Kθ = 1 in general. Thus, from Figure 10 and Figure 11, the case with Kθ = 2 extracts more energy.
To analyze the difference of the Cl curve at different non-sinusoidal pitching motions in Figure 11, we discuss the evolution of the vortex field around parallel foils in the cases of Kθ = 1 and Kθ = 2 during t/T = 0–0.5. Comparing the vortex field in Figure 12a,b, it is found that the vortex structures at Kθ = 2 are bigger than those at Kθ = 1 during every moment. The bigger low-pressure area is formed because of the bigger vortex structure, increasing the lift force. So Cl at Kθ = 2 is higher than that at Kθ = 1 in Figure 11. The Cl curve in Figure 11b has a significant fluctuation near t/T = 0.05, combining with Figure 12b. It is concluded that a big trailing edge vortex sheds at this time, causing a sharp change of Cl. Meanwhile, the leading edge vortex is formed earlier at Kθ = 2.

4. Conclusions

In this paper, we considered the effect of non-sinusoidal motion on the energy extraction performance of parallel foils. The energy extraction processes of two oscillation motions (the combined non-sinusoidal plunging and sinusoidal pitching motion, as well as the combined non-sinusoidal pitching and sinusoidal plunging motion) with Re = 1100 are simulated. Five non-sinusoidal motions are discussed and the optimal motion, the optimal frequency and the average power coefficient for these two combined oscillation motions are obtained. The contribution of plunging oscillation and pitching oscillation to energy extraction performance of foil is analyzed and the relationship between the lift force and the vortex structure is investigated.
The results show that both the non-sinusoidal plunging motion and the non-sinusoidal pitching motion have an obvious effect on the energy extraction characteristics of foil. When the foils oscillate as the combined non-sinusoidal plunging and sinusoidal pitching motion, the optimal extraction performance of foil appears when Kh is −0.5; the highest power output CPm and extraction efficiency η are 0.375 and 0.188, respectively. When k < 1.0, the combined motion with Kh = −0.5 shows a better extraction performance, while the better performance appears at the combined motion with Kh = 1 when k > 1.0. At the same k, the more the plunging motion approximates to the sinusoidal motion, the more energy is extracted by the parallel foils.
When the motion is the combined non-sinusoidal pitching and sinusoidal plunging, the peak values of CPm under four different motions (except Kθ = −1) are gained at k = 0.8. This means that the optimal frequency which can gain the best extraction performance of foil in most motions is 0.8. The optimal CPm and η are achieved at Kθ = 2, which are 0.623 and 0.312, respectively. At a high frequency, the effect of the non-sinusoidal pitching motion on the extraction performance is obvious. With the increase in k, the non-sinusoidal pitching motion at Kθ = 1 and Kθ = 2 cause a high consumption of energy. That is to say, the more the non-sinusoidal pitching motion is similar to the square wave, the worse extraction performance gained.
According to the vortex fields under different non-sinusoidal pitching motions, it was found that the vortex structure at Kθ = 2 is much bigger than that at Kθ = 1 at every moment. The lift force of foil is higher at Kθ = 2, and the synchronicity between Cl and Vy at these two motions are similar. As a result, the oscillation process of the foils extracts more energy when Kθ is 2.

Author Contributions

Y.W. and Y.X. established the analysis method of the energy extraction performance of parallel foils; Y.W. performed the numerical investigation and obtained the conclusions. F.Z. established the computational model, discussed the results and wrote the paper.

Funding

This research was funded by the National Natural Science Foundation of China, project number is 11872294.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Afoil’s oscillation area, m2
clength of foil, m
Cllift coefficient
Cmtorque coefficient
CPpower coefficient
ffrequency, Hz
Fyaerodynamic lift force in y direction, N
h(t)plunging motion
H0non-dimensional plunging amplitude
kreduced frequency
Khcontrol parameter of plunging motion
Kθcontrol parameter of pitching motion
Mtorque, N·m
Ppower, W
ReReynolds number
StStrouhal number
ttime, s
Toscillation cycle, s
Ufree stream velocity, m·s−1
Vy(t)plunging velocity, m·s−1
α(t)pitching motion
α0nominal angle of attack
δspace between parallel foils, m
ηefficiency
θ0pitching amplitude, °
ρfluid density, kg·m−3
φphase difference between plunging and pitching motion, °
ωpitching velocity, rad·s−1

Superscript

hplunging motion
maverage value

References

  1. McKinney, W.; DeLaurier, J. The wingmill: An oscillating-wing windmill. J. Energy 1981, 5, 109–115. [Google Scholar] [CrossRef]
  2. Campobasso, M.S.; Drofelnik, J. Compressible Navier-Stokes analysis of an oscillating wing in a power-extraction regime using efficient low-speed preconditioning. Comput. Fluids 2012, 67, 26–40. [Google Scholar] [CrossRef]
  3. Abiru, H.; Yoshitake, A. Study on a flapping wing hydroelectric power generation system. J. Environ. Eng. 2011, 6, 178–186. [Google Scholar] [CrossRef]
  4. Shimizu, E.; Isogai, K.; Obayashi, S. Multiobjective design study of a flapping wing power generator. J. Fluids Eng. 2008, 130, 021104. [Google Scholar] [CrossRef]
  5. Han, R.; Cherry, J.; Kallenberg, R. Modeling an oscillating water foil for hydro-kinetic power generator using COMSOL 3.5a. In Proceedings of the COMSOL Conference, Boston, MA, USA, 7–9 October 2010. [Google Scholar]
  6. Jones, K.D.; Platzer, M.F. Numerical computation of flapping-wing propulsion and power extraction. In Proceedings of the 35th Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 6–9 January 1997. [Google Scholar]
  7. Jones, K.D.; Lindsey, K.; Platzer, M.F. An investigation of the fluid-structure interaction in an oscillating-wing micro-hydropower generator. Phys. Atom. Nucl. 2003, 65, 325–328. [Google Scholar]
  8. Kim, J.; Quang, L.T.; Hwan, K.J.; Sitorus, P.E.; Tambunan, I.H.; Kang, T. Experimental and numerical study of a dual configuration for a flapping tidal current generator. Bioinspir. Biomim. 2015, 10, 046015. [Google Scholar] [CrossRef] [PubMed]
  9. Dumas, G.; Kinsey, T. Eulerian simulations of oscillating airfoils in power extraction regime. Adv. Fluid Mech. VI 2006, 245–254. [Google Scholar] [CrossRef]
  10. Kinsey, T.; Dumas, G. Parametric study of an oscillating airfoil in a power-extraction regime. AIAA J. 2008, 46, 1318–1330. [Google Scholar] [CrossRef]
  11. Tuncer, I.H.; Kaya, M. Optimization of flapping airfoils for maximum thrust and propulsive efficiency. AIAA J. 2005, 43, 2329–2336. [Google Scholar] [CrossRef]
  12. Bryant, M.; Shafer, M.W.; Garcia, E. Power and efficiency analysis of a flapping wing wind energy harvester. Proc. SPIE 2012, 8341, 83410E. [Google Scholar] [CrossRef]
  13. Sitorus, P.E.; Le, T.Q.; Ko, J.H.; Truong, Q.T.; Tambunan, I.H.; Kang, T. Progress on development of a lab-scale flapping-type tidal energy harvesting system in KIOST. In Proceedings of the 2013 IEEE Conference on Clean Energy and Technology, Lankgkawi, Malaysia, 18–20 November 2013. [Google Scholar]
  14. Platzer, M.F.; Sarigul-Klijn, N. A new oscillating-foil power generator for sailingship-based renewable energy generation. In Proceedings of the ASME 2010 4th International Conference on Energy Sustainability, Phoenix, AZ, USA, 17–22 May 2010. [Google Scholar]
  15. Platzer, M.F.; Ashraf, M.A.; Young, J. Extracting power in jet streams: Pushing the performance of flapping-wing technology. In Proceedings of the 27th International Congress of the Aeronautical Sciences, Nice, France, 19–24 September 2010. [Google Scholar]
  16. Platzer, M.F.; Young, J.; Lai, J.C.S. Flapping-wing technology: The potential for air vehicle propulsion and airborne power generation. In Proceedings of the 26th International Congress of the Aeronautical Sciences, Anchorge, AK, USA, 14–19 September 2008. [Google Scholar]
  17. Xiao, Q.; Liao, W.; Yang, S. How motion trajectory affects the energy extraction performance of an oscillating foil. In Proceedings of the 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 6 January 2010. [Google Scholar]
  18. Xiao, Q.; Liao, W.; Yang, S.; Yan, P. How motion trajectory affects energy extraction performance of a biomimic energy generator with an oscillating foil. Renew. Energy 2012, 37, 61–75. [Google Scholar] [CrossRef]
  19. Ashraf, M.A.; Young, J.; Lai, J.C.S.; Platzer, M.F. Numerical analysis of an oscillating-wing wind and hydropower generator. AIAA J. 2011, 49, 1374–1386. [Google Scholar] [CrossRef]
  20. Wang, Y.L.; Jiang, W.; Xie, Y.H. Numerical investigation into the effects of motion parameters on energy extraction of the parallel foils. J. Fluids Eng. 2018, 141, 061104. [Google Scholar] [CrossRef]
  21. Young, J.; Lai, J.C.S. Mechanisms influencing the efficiency of oscillating airfoil propulsion. AIAA J. 2007, 45, 1695–1702. [Google Scholar] [CrossRef]
Figure 1. The structure and oscillation process of parallel foils: (a) t/T = 0–0.5; (b) t/T = 0.5–1.
Figure 1. The structure and oscillation process of parallel foils: (a) t/T = 0–0.5; (b) t/T = 0.5–1.
Applsci 09 00384 g001
Figure 2. The computation domains arrangement.
Figure 2. The computation domains arrangement.
Applsci 09 00384 g002
Figure 3. The mesh strategy.
Figure 3. The mesh strategy.
Applsci 09 00384 g003
Figure 4. The curves under the different non-sinusoidal motions at k = 0.8: (a) non-sinusoidal plunging motion; (b) non-sinusoidal pitching motion.
Figure 4. The curves under the different non-sinusoidal motions at k = 0.8: (a) non-sinusoidal plunging motion; (b) non-sinusoidal pitching motion.
Applsci 09 00384 g004
Figure 5. The variations of CP under different mesh strategies and time steps.
Figure 5. The variations of CP under different mesh strategies and time steps.
Applsci 09 00384 g005
Figure 6. The CPm and η of parallel foils at different Kh: (a) CPm curves under different k; (b) η curves under different k.
Figure 6. The CPm and η of parallel foils at different Kh: (a) CPm curves under different k; (b) η curves under different k.
Applsci 09 00384 g006
Figure 7. The CPhm and CPθm of parallel foils at different Kh: (a) CPhm curves under different k; (b) CPθm curves under different k.
Figure 7. The CPhm and CPθm of parallel foils at different Kh: (a) CPhm curves under different k; (b) CPθm curves under different k.
Applsci 09 00384 g007
Figure 8. The Cl curves of the upper foil under different Kh at k = 0.8: (a) Kh = 0; (b) Kh = 2.
Figure 8. The Cl curves of the upper foil under different Kh at k = 0.8: (a) Kh = 0; (b) Kh = 2.
Applsci 09 00384 g008
Figure 9. The CPm and η of parallel foils at different Kθ: (a) CPm curves under different k; (b) η curves under different k.
Figure 9. The CPm and η of parallel foils at different Kθ: (a) CPm curves under different k; (b) η curves under different k.
Applsci 09 00384 g009
Figure 10. The CPhm and CPθm of parallel foils at different Kθ: (a) CPhm curves under different k; (b) CPθm curves under different k.
Figure 10. The CPhm and CPθm of parallel foils at different Kθ: (a) CPhm curves under different k; (b) CPθm curves under different k.
Applsci 09 00384 g010
Figure 11. Cl curves of upper foil under different Kθ at k = 0.8: (a) Kθ = 1; (b) Kθ = 2.
Figure 11. Cl curves of upper foil under different Kθ at k = 0.8: (a) Kθ = 1; (b) Kθ = 2.
Applsci 09 00384 g011
Figure 12. The vortex field during t/T = 0–0.5 under different Kθ at k = 0.8: (a) Kθ = 1; (b) Kθ = 2.
Figure 12. The vortex field during t/T = 0–0.5 under different Kθ at k = 0.8: (a) Kθ = 1; (b) Kθ = 2.
Applsci 09 00384 g012

Share and Cite

MDPI and ACS Style

Wang, Y.; Zhu, F.; Xie, Y. Numerical Study on the Effect of Non-Sinusoidal Motion on the Energy Extraction Performance of Parallel Foils. Appl. Sci. 2019, 9, 384. https://doi.org/10.3390/app9030384

AMA Style

Wang Y, Zhu F, Xie Y. Numerical Study on the Effect of Non-Sinusoidal Motion on the Energy Extraction Performance of Parallel Foils. Applied Sciences. 2019; 9(3):384. https://doi.org/10.3390/app9030384

Chicago/Turabian Style

Wang, Yulu, Fahui Zhu, and Yonghui Xie. 2019. "Numerical Study on the Effect of Non-Sinusoidal Motion on the Energy Extraction Performance of Parallel Foils" Applied Sciences 9, no. 3: 384. https://doi.org/10.3390/app9030384

APA Style

Wang, Y., Zhu, F., & Xie, Y. (2019). Numerical Study on the Effect of Non-Sinusoidal Motion on the Energy Extraction Performance of Parallel Foils. Applied Sciences, 9(3), 384. https://doi.org/10.3390/app9030384

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop