Next Article in Journal
Equilibrium Interaction Strategies for Integrated Energy System Incorporating Demand-Side Management Based on Stackelberg Game Approach
Previous Article in Journal
Energy Management of Green Port Multi-Energy Microgrid Based on Fuzzy Logic Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Fluid Dynamics Simulation on Blade Geometry of Novel Axial FlowTurbine for Wave Energy Extraction

Department of Mechanical Engineering, North Carolina A&T State University, Greensboro, NC 27411, USA
*
Author to whom correspondence should be addressed.
Energies 2024, 17(14), 3602; https://doi.org/10.3390/en17143602
Submission received: 12 June 2024 / Revised: 11 July 2024 / Accepted: 16 July 2024 / Published: 22 July 2024
(This article belongs to the Topic Energy from Sea Waves)

Abstract

:
Wave energy converters (WECs) utilizing the Oscillating Water Column (OWC) principle have gained prominence for harnessing kinetic energy from ocean waves. This study explores an innovative approach by transforming the pivoting Denniss–Auld turbine blades into a fixed configuration, offering a simplified alternative. The fixed-blade design emulates the maximum pivot points during the OWC’s exhalation and inhalation phases. Traditional Denniss–Auld turbines rely on complex hub systems to enable controllable blade rotation for performance optimization. This research examines the turbine’s efficiency without mechanical actuation. The simulations were conducted using ANSYS™ CFX 2023 R2 to solve the three-dimensional, incompressible, steady-state Reynolds-Averaged Navier–Stokes (RANS) equations, employing the k-ω SST turbulence model to close the system of equations. A grid convergence study was performed, and the numerical results were validated against available experimental and numerical data. An in-depth analysis of the intricate flow field around the turbine blades was also conducted. The modified Denniss–Auld turbine demonstrated a broad operating range, avoiding stalling at high flow coefficients and exhibiting performance characteristics like an impulse turbine. However, the peak efficiency was 12%, significantly lower than that of conventional Denniss–Auld and impulse turbines. Future research should focus on expanding the design space through parametric studies to enhance turbine efficiency and power output.

1. Introduction

The rapid expansion of the world economy and population is driving a surge in global energy demand. Consequently, greenhouse gas emissions have significantly increased, leading to ecological challenges and negatively impacting the Earth’s climate. Therefore, it is crucial to explore clean energy solutions incorporating renewable energy sources to address the growing energy needs and environmental concerns.
Renewable energy sources have the potential to meet two-thirds of global energy demand and help mitigate global temperature rise by preventing greenhouse gas emissions. Reducing greenhouse gas emissions from now until 2050 is essential to limiting the average global surface temperature increase to 2 °C [1]. Consequently, global leaders and policymakers have implemented measures to adopt and integrate renewable energy technologies in response to climate change.
Ocean wave energy has emerged as one of the most promising eco-friendly alternatives to fossil fuels among the various renewable energy sources. Over the past two decades, it has garnered significant interest due to its potential to substantially supplement the electricity demand of many countries [2]. Furthermore, the high energy density, minimal ecological impacts, and negligible energy losses of waves traveling long distances have intrigued researchers, driving the development of feasible technologies to harness wave energy [3]. Current research indicates that the potential of ocean wave energy is approximately 30,000 TWh per year. If harvested efficiently, this could meet global electricity demands [4]. Wave energy converters (WECs) extract energy from waves, converting it into mechanical energy and electrical energy. The commercialization of wave energy hinges on enhancing the performance of WECs [5]. In recent decades, technologies such as the Oscillating Water Column (OWC), floating buoys, point absorbers, and overtopping devices have been extensively explored for wave energy conversion. However, due to their simplicity, OWC WECs are the most widely used devices for converting wave energy [6].
In an OWC, the rotor is the only moving part, rotating at high rpm and connected to an electrical generator to generate electricity. The OWC utilizes the oscillatory motion of ocean waves, where air becomes trapped above the free water surface and undergoes cycles of inhalation and exhalation. This results in the OWC system experiencing unsteady bi-directional airflow, necessitating specialized Power Take-Off (PTO) systems capable of managing this complex flow. According to recent research, the most widely used turbines in OWCs are the Wells [7] and impulse-type concepts [8].
The Wells turbine stands out among axial flow air turbines as one of the most economical and straightforward methods for harnessing wave energy. Proposed by Professor Wells, this self-rectifying axial flow turbine operates within an OWC using compressed air. Its rotor comprises several symmetrical airfoil blades arranged radially and staggered at 90 degrees around the hub (see Figure 1). Due to the symmetric orientation of the blades, the tangential force acting on them remains consistent regardless of the airflow direction, as depicted in Figure 2. Consequently, the Wells turbine consistently rotates in the same direction, even when exposed to bi-directional airflow.
Figure 1. Schematic diagram of OWC and Wells turbine [9]. (a) OWC. (b) Wells turbine.
Figure 1. Schematic diagram of OWC and Wells turbine [9]. (a) OWC. (b) Wells turbine.
Energies 17 03602 g001
Despite its advantages, the Wells turbine suffers from inherent shortcomings such as a narrow operating range, a low aerodynamic efficiency, high noise levels, a high axial force coefficient, a low tangential force coefficient, and poor starting characteristics compared to conventional turbines (Figure 3). Consequently, several self-rectifying air turbines have been developed, analyzed, and continuously improved to address these limitations.
Wave energy converters utilizing self-rectifying impulse turbines are gaining popularity as viable alternatives to the Wells turbine [10]. Two sets of guide vanes are symmetrically arranged on both sides of the rotor in bi-directional turbines, whereas uni-directional turbines have only one set.
Setoguchi et al. [11] examined the current advancements in wave energy conversion utilizing self-rectifying air turbines. They conducted numerical simulations to evaluate the overall performance of these turbines under irregular sea wave conditions, focusing on their starting and running characteristics. The study found that impulse turbines demonstrated superior running performance compared to the Wells turbine. Due to their rapid start capability, impulse turbines have a longer electricity generation duration than Wells turbines. Moreover, they operate at lower rotational speeds due to higher torque coefficients under load-free conditions than Wells turbines.
Figure 2. Uni-directional rotation of the Wells turbine blade [12].
Figure 2. Uni-directional rotation of the Wells turbine blade [12].
Energies 17 03602 g002
Figure 3. Operating curve of Denniss–Auld turbine [12].
Figure 3. Operating curve of Denniss–Auld turbine [12].
Energies 17 03602 g003
Falcao et al. [13] assessed the performance of different self-rectifying turbines across various ocean wave energy levels and OWC applications. The study revealed that the Wells turbine can achieve a peak efficiency of 75%, a benchmark surpassed only by the most advanced impulse turbines. However, in energetic sea conditions, the average efficiency of the Wells turbine drops significantly compared to that of impulse turbines. It is generally observed that the average efficiency of impulse turbines remains consistent regardless of sea state. However, it typically relies on movable guide vanes for optimal performance in all conditions.
Hu et al. [14] conducted a numerical study on both steady and unsteady flow characteristics in a Wells turbine. They observed a stall near the blade tip under steady flow conditions, while radial flow helped maintain attached boundary layers near the hub. The study investigated two control methods for unsteady flow: constant angular velocity and constant damping moment. Regarding constant angular velocity, the authors noted the hysteresis effect impacting the Wells turbine, particularly at high frequencies of sea waves.
Halder et al. [15] performed a multi-objective optimization of blade sweep for a Wells turbine. As a result, they achieved a 28.28% increase in the peak torque coefficient. However, this improvement came at the cost of a 13.5% decrease in efficiency.
Das et al. [16] experimentally studied a biplane Wells turbine under varying loads and periodic airflows. They observed that during the turbine’s acceleration and deceleration phases, the pressure drop and torque coefficients exhibited hysteresis behavior, attributed to the OWC turbine system’s capacitive characteristics.
Wang et al. [17] introduced an innovative biomimetic design inspired by hawkmoth wings for the Wells turbine, aiming to enhance its efficiency in energy harvesting. The study found that the biomimetic design modifies the flow field and can achieve higher peak efficiency, especially at high angles of attack.
Folley et al. [18] assessed the performance of the contra-rotating Wells turbine installed in the LIMPET wave power station based on a theoretical analysis and model tests. The study suggests that a biplane or monoplane Wells turbine with guide vanes may offer better performance for OWC applications than the contra-rotating Wells turbine, which exhibited reduced efficiency in their evaluations.
In response to the limitations of Wells and impulse turbines, researchers have developed alternative self-rectifying axial flow turbines such as the HydroAir turbine and the Denniss–Auld turbine. The HydroAir turbine addresses these challenges by increasing the distance between the guide vane rows and rotor blades, reducing the excessive flow angle encountered at the second row of guide vanes [19]. Decreasing the flow velocity at the entrance to the second row of guide vanes can reduce kinetic energy losses.
Another type of self-rectifying axial flow turbine is the Denniss–Auld turbine [20,21]. The Denniss–Auld turbine and the variable-pitch Wells turbine share similar characteristics but differ in the angles at which their blades are set (Figure 4). For instance, the rotor blades are positioned in the range of α < γ < π α and θ < γ < θ for the Denniss–Auld turbine and the variable-pitch Wells turbine, respectively, where α 20 ° 35 ° ,   θ 25 ° .
Furthermore, in the Wells turbine, the flow attaches to the blade’s leading edge. In contrast, the Denniss–Auld rotor blades alternate between acting as leading edges and trailing edges depending on the direction of fluid flow. This characteristic necessitates that the rotor blades of the Denniss–Auld turbine are designed with identical edges.
When the OWC system inhales or exhales, the rotor blades of the Denniss–Auld turbine must pivot rapidly between their extreme positions. In contrast, Wells turbine rotor blades require smaller angular ranges of blade pitch to operate effectively.
The blade pitching angles are critical to the aerodynamic performance of a Denniss–Auld turbine. Figure 3 shows the operating curve of the Denniss–Auld turbine for variable blade pitch angles, γ. According to the analytical model prediction using blade element theory, the peak efficiency occurs at around γ = 40° and decreases with an increase in γ. So, the optimum blade angle should be between 20° and 40°. This entails the basis of our design comprising γ = 30° (Figure 5) for improved performance.
This study addresses the difficulties mentioned above by redesigning the Denniss–Auld turbine blades to reduce mechanical complexity by preventing pivoting. In addition, we have investigated the aerodynamic performance characteristics of a modified Denniss–Auld turbine comprising novel airfoil-shaped blades using ANSYS™ CFX 2023 R2 by solving the steady-state, incompressible, three-dimensional Reynolds-Averaged Navier–Stokes (RANS) equations coupled with the k- ω SST turbulence model [22,23,24]. Ultimately, the novelty of this research is the design of a modified Denniss–Auld turbine without pivoting with the goal of mitigating stall at high flow coefficients. Under steady-state flow conditions, we determined the turbine’s performance in non-dimensional coefficients, i.e., torque coefficient, efficiency, and input coefficient, for a wide range of flow coefficients. Numerical accuracy was achieved through a grid independence study [4,22,25,26,27,28,29]. Furthermore, we validated our numerical modeling by comparing the results to available experimental [30] and numerical results [31] for a baseline Wells turbine with a uniform tip clearance.

2. Numerical Methods

Numerical simulations for this study were performed using the pressure-based solver within ANSYS™ CFX, where velocity and pressure equations are solved simultaneously by treating them as a single equation system. In addition, a fully implicit discretization approach was undertaken to discretize the equations for computational modeling. For instance, the advection and diffusion terms in the RANS equations were treated with a second-order scheme and a shape function-based approach, respectively. The details of the equations can be found in reference [32].
When computing turbulent flow, its pressure and velocity field decompose into mean and fluctuating parts. RANS equations, governing mean flow, are derived from averaging the Navier-Stokes equations over time. However, due to the inherently nonlinear nature of the Navier-Stokes equations, the velocity field still experiences fluctuations resulting from the convective acceleration (Equation (1)) [33]. As a result, a nonlinear term, the Reynolds stress, appears, which Boussinesq considers a function of the mean flow component concerning the closure problem (Equation (2)) [33].
This study adopts the k- ω SST turbulence model since it predicts flow separation more accurately than other RANS models. In addition, the model can also resolve complex boundary layer flows with adverse pressure gradients (e.g., Turbomachinery).
ρ U t + . ρ U U = p + . μ U + U T + ρ g 2 3 μ . U . ρ U U ¯ R e y n o l d s S t r e s s ,
ρ U U ¯ = ν t U + U T M e a n V e l o c i t y G r a d i e n t s 2 3 ρ k δ i , j ,
where δ i , j is the Kronecker delta.

2.1. Numerical Modeling

The performance characteristics of a Wells turbine comprising symmetrical NACA0015 blades, fixed at a 90-degree staggered angle with uniform tip clearance, were assessed numerically and compared to existing data. Because tip leakage is among the most significant secondary fluid flow phenomena found in Turbomachinery, tip clearance has been used to predict turbine performance accurately. The specifications of the simulated turbine are given in Table 1.
When performing CFD analyses, it is necessary to discretize the computational domain to solve the governing equations for fluid flow. Therefore, the CFD results are prone to discretization errors. The quality of the computational grid can often become the primary cause of the error. Figure 6 depicts the unstructured mesh used in this study.
Within ANSYS meshing, we can specify a cell size within a volume. Even though the mesh is unstructured and the cell sizes vary slightly, we selected four grid sizes for the mesh resolution within our control volume: coarse (10 mm), medium (5 mm), fine (2.5 mm), and extra-fine (1.25 mm).
It is necessary to resolve the boundary layer flow adequately around the rotor blades to determine the dynamic stall point of the turbine accurately. For this reason, hexahedron inflation layers were densely placed around the blade surface to resolve the boundary layer. The computational domain extends four and six times the axial chord length in the leading- and trailing-edge directions, respectively. The boundary condition at the inlet is set to a uniform axial velocity with a specified turbulence intensity of 5%. A uniform static pressure boundary condition is imposed at the outlet of the computational domain. No-slip wall boundary conditions are set on the blade surface, hub, and tip.
Considering that the Wells turbine has a symmetrical geometry, the computations are limited to one blade-to-blade passage along with periodic boundary conditions, as shown in Figure 7. The details of the meshing and boundary conditions are provided in Table 2. From the definition, the dimensionless wall distance, y+, is given by the following [31]:
y + = y u T ν ,
Where y , u T , a n d   ν signify the absolute first layer distance from the wall, friction velocity, and kinematic viscosity of the fluid, respectively. The k-ω SST turbulence model, because it is sensitive to the mesh, requires high resolution with a value of y+ < 1 to model the viscous sublayer region. Accounting for the required y+ value, the first prism layer height obtained was 1.1 × 10−5 m. There was a total of 20 inflation layers, with a growth factor of 1.2.
In the fluid zone, the Moving Reference Frame (MRF) was used, which rotates at the same speed as the turbine rotor. Hence, the following conservation equations were solved:
I.
The conservation of mass
ρ t + . ρ v r = 0 ,
II.
The conservation of momentum
ρ v t + . ρ v r v + ρ [ ( ω a ^ ) X v v t ] = p + . τ _ ,
Here, v t , v   , ω , τ _ , and a ^ indicate the transitional velocity, absolute velocity, angular velocity, viscous stress, and a unit vector defining the axis of rotation, respectively.
The performance characteristics of the turbine were assessed for various flow coefficients, ϕ , where
ϕ = v U t i p ,
Since we used pseudo-transient simulations in ANSYS CFXTM, time step values were required even with steady flow. Therefore, for the first few flow configurations ( ϕ < 0.2) far from the blade stall, the time step size, Δt, was set as 0.05 s, selected from the Auto Timescale option [32]. However, the results were oscillatory and unsteady for the pre- and post-stall conditions ( ϕ > 0.2). So, we reduced to Δt = 0.005 s; therefore, the fluctuations decreased, leading to convergence.
ANSYS CFX uses second-order accurate approximations. The second-order backward Euler scheme achieves time discretization with implicit and conservative time stepping. For additional details, please see reference [34].
The variation in the flow coefficient is achieved by changing the axial velocity, v, whereas the circumferential velocity at the tip radius, Utip, remains constant. The root mean square (RMS) residuals of the governing equations are monitored to determine convergence. As part of the convergence verification, torque output and pressure drop across the turbine are also observed. Sometimes, set values of the convergence criteria of each residual are lowered to guarantee that the monitored quantities remain constant. The net mass imbalance in this study was less than 0.058 percent.

2.2. Validation of the Numerical Model

Our research considers seven steady-state flow conditions, each characterized by an associated flow coefficient, ϕ. The turbine operating range is chosen as 0.075 ϕ 0.275 , or equivalently 4 ° α 15 ° , where α is the incident angle of attack.
According to the literature [3], the following non-dimensional coefficients have been used to describe turbine performance:
The torque coefficient, CT,
C T = T ρ ω 2 R 5 ,
The pressure drop coefficient, Δ P 0 * ,
Δ P 0 * = Δ p 0 ρ ω 2 R 2 ,
The input coefficient, CA,
C A = 2 Δ p 0 Q ρ v 2 + U R 2 b c z v ,
The efficiency, η ,
η = T ω Δ p 0 Q ,
where Q , ρ , T , Δ p 0 , ω , z , U R , c , b , a n d   v represent the volumetric flow rate, air density, blade torque, static pressure drops across the turbine, angular velocity, number of blades, circumferential velocity at blade mean radius, blade chord length, blade height, and axial inlet velocity, respectively.
The Reynolds number is defined as follows:
R e = ρ c v 2 + U t i p 2 μ ,
A grid independence study is essential for obtaining accurate results in numerical simulations. Therefore, we picked four different cell sizes to study the effect of grid resolution on the results. Then, we compared the numerical and experimental results of the computed torque coefficients for each mesh at ϕ = 0.225 .
As evident in Figure 8, there is a substantial drop in the torque coefficient below the medium grid. The discrepancy between the medium, fine, and extra-fine grids was ~24%, ~4%, and ~0.5%, respectively. Due to the negligible variation in the torque coefficient between the fine and extra-fine grid, the fine grid was used for subsequent simulations. To simulate the flow, the workstation specifications of Intel (R) Xeon (R) Gold 6258 R CPU @ 2.7 GHz and 2.69 GHz (2 processors) with a 56-core system were used. The equipment was sourced from Dell, based in Round Rock, TX, USA. A table listing the cell size statistics and a figure showing the grid independence study are provided in Table 3 and Figure 8, respectively.
We selected the torque coefficient as the parameter evaluated for each mesh to reach grid independence, consistent with the literature [35,36]. However, some studies [37] also employed efficiency too. In either case, we see a negligible change in the results for the extra-fine grid (See Table 4).

3. Results and Discussion

This study examined turbine performance for a constant rotational speed of 2000 rev/min at various inlet velocities. The flow coefficient ϕ captures the operating range of the turbine, where 0.075 ϕ 0.275 . Detailed information on the range of the inlet velocity, flow coefficient, and Reynolds number is shown in Table 5.
In the current CFD simulations, we compare the experimental and numerical data. Comparing the computed results to the experiment [30] and CFD [31] data, Figure 9 illustrates a high level of agreement up to ϕ = 0.2 for the pressure drop coefficient ( Δ P 0 * ), torque coefficient (CT), and efficiency ( η ). But for ϕ > 0.2 , i.e., at a high angle of attack, deep stalling is evident from the rapid drop in blade torque and turbine efficiency. This is a result of leading-edge flow separation. Furthermore, the RANS simulations overestimate the torque. This is typical of even well-resolved RANS computations; they are often inaccurate past the blade stall point [38]. Solving complex three-dimensional separated flow is a principal factor in the resulting deviation in the numerical and experimental results. Yet, under the same geometrical conditions, the present work is consistent with existing CFD results [31].
After validating the numerical procedure, we simulated our modified Denniss–Auld turbine. The simulations were carried out for a range of axial inlet flow velocities and rotational speeds corresponding to a flow coefficient range of 0.41 ≤ ϕ ≤ 2.25. Figure 10, Figure 11 and Figure 12 show that the turbine fails to produce positive torque at the lower end of the flow coefficient range, resulting in negative efficiency, suggesting that it cannot generate power (i.e., torque) corresponding to the input rotational speed (rpm). Thus, the negative efficiency implies that the given rpm and inlet velocity do not physically allow the turbine to rotate. Due to its inability to produce positive torque at a low flow coefficient, ϕ < 0.44 , we simulated the turbine at high flow coefficients, which correspond essentially to flow conditions at a high angle of attack.
This study could have considered reporting the results in terms of different Tip Speed Ratios (TSRs), λ. TSR is defined as follows:
λ = U t i p v ,
By comparing Equations (6) and (12), we can deduce
ϕ = 1 λ ,
Consequently, the results we demonstrated in terms of ϕ should be similar if expressed through λ .
At higher flow coefficients, the modified Denniss–Auld turbine has similarities to the impulse turbine. For instance, the torque coefficient increases with the flow coefficient; this increase is similar in behavior to the impulse turbine. The turbine blades do not reach stall at higher flow coefficients, contrary to Wells turbines, as shown in Figure 10. In addition, the slopes of the torque coefficient curves are milder at low flow coefficients but become steeper at higher flow coefficients.
Figure 11 shows that the modified Denniss–Auld turbine maintains a relatively constant efficiency for a wide range of flow coefficients, like the impulse turbine. Nonetheless, it has a lower efficiency compared to the impulse turbine with fixed guide vanes. Moreover, this figure shows an increase in turbine efficiency with increasing axial inlet velocity.
While Wells turbines have a higher peak efficiency than other turbines, their operating range is generally lower. However, this modified Denniss–Auld turbine has a wider operating range depending on the incident flow velocity (Figure 12). In other words, the turbines can achieve high rotational speeds if the axial flow velocity is high enough. The results are comparable to the article by Finnigan et al. [21], according to Figure 13 below. It shows that the maximum lift coefficient increases with the Reynolds number; the high lift coefficient causes the rise in turbine torque.
Figure 14 illustrates how the Reynolds number affects the turbine’s peak efficiency. Variations in the Reynolds number ranged from 6.5 × 104 to 3.1 × 105. We observed a substantially linear increase in peak efficiency with various flow coefficients.
In Figure 15, the torque coefficient is plotted as a function of turbine rotational speed at a given inlet axial flow velocity. It is evident from this figure that increasing the turbine’s rotational speed for a certain axial velocity reduces the torque coefficient drastically. A steep slope is observed at an incident speed of 4.7 m/s, but it becomes more gradual with increasing rotational speeds and axial velocities. In addition, the curves shift rightward, accommodating the high rotational speed associated with an increase in inlet velocity, which suggests the turbine, in turn, can operate within a broader range. Thus, the turbine can experience higher rotational speeds when subjected to greater axial velocity at the inlet.
Figure 16 shows a plot of the input coefficient versus the flow coefficient. The input coefficients are lower when the flow coefficients are smaller; therefore, the turbine harnesses less power at the lower range of flow coefficients. At a lower coefficient range ( ϕ < 0.9 ), the slope of the input coefficient is steeper, so the rate of change in power production is sharper than at a flow coefficient range of ϕ > 0.9 . Inlets of varying axial velocities show similar trends.
According to Figure 17, the input coefficient drops exponentially with increasing turbine rotational speed. Also, a decline in the slope of the curves is associated with an increase in the incident axial velocity.
The torque coefficient steadily increases with the Reynolds number at peak efficiency, as shown in Figure 18. It is evident that the input coefficients increase until a Reynolds number of 1 × 105 and then remain constant. This results in a gradual increase in peak efficiency.
Figure 19 provides a comparison of the modified Denniss–Auld turbine with the standard Denniss–Auld turbine in terms of non-dimensional coefficients for a wide range of flow conditions. For γ = 20   a n d   40 , the CT and CA increase almost linearly with ϕ , and the slope decreases with increasing γ .
On the contrary, the CT rises exponentially with ϕ , which means that the modified Denniss–Auld turbine produces much higher torque than the standard design at higher ϕ . CA demonstrates logarithmic growth; the sharper increase occurs in the lower range of ϕ , and the gains become flattened out for higher ϕ .
The trend of the efficiency plots is similar for all the γ , except for γ = 20 . The reasoning behind this could be attributed to the shortcomings associated with the existing test rig to resolve the peak efficiency. The modified Denniss–Auld turbine produces the lowest aerodynamic efficiency over the entire range of ϕ . The peak efficiency of the standard Denniss–Auld turbine is 63% and occurs for γ = 40 at ϕ = 1 . However, the modified Denniss–Auld turbine has a peak efficiency of 11.73% that occurs at ϕ = 0.9 .
Figure 20 and Figure 21 show the coefficient of pressure for a varying flow coefficient at an inlet velocity of 7.854 m/s. The results show that with the increasing flow coefficient, the peak pressure point moves downstream on the pressure side of the turbine. The peak efficiency of the turbine is at a flow coefficient of 0.83. This is a result of the flow attachment point no longer being on the leading edge of the turbine. This is evident by the inflection point and decreases in pressure at x/c = 0.042 at a flow coefficient of 1.25 and 2.5. Increasing the flow coefficient to 2.5 increases the inflection point and causes a further decrease in pressure. This results in the decrease in turbine efficiency discussed earlier. Furthermore, the rounded shape of the trailing edge of the turbine causes the pressure side to experience a severe decrease in pressure leading to lower pressure compared to the suction side.
Figure 22 shows a plot of Cp at a flow coefficient of 0.42 at 40% and 80% span of the turbine. The plot at 40% span shows a slightly higher pressure difference indicating better aerodynamic performance compared to the 80% span location. However, the trends of both plots are similar.
Figure 23 and Figure 24 show the pressure with increasing flow coefficient at 40% and 80% span of the turbine. The contours at 40% span illustrate the change in pressure distribution on the turbine at various flow coefficients (Figure 23). These results show that at a flow coefficient below 1.25 the peak pressure occurs on the leading edge. Above ϕ = 1.25, the peak surface pressure is downstream of the leading edge (pressure side of the airfoil), leading to a decrease in turbine efficiency. Similarly, the same behavior is displayed at the 80% span position (Figure 24).
Due to the surface curvature of the modified Denniss–Auld turbine, the torque coefficient increases exponentially throughout the flow coefficient range investigated in this study. However, a penalty occurs in the corresponding pressure drops.
A comparison of the rotational speed and torque versus the flow coefficient and various inlet velocities is shown in Figure 25. The results show that the torque increases sharply at a flow coefficient of between 0.5 and 1.0. From there, the torque increases but less substantially. The intersection between the rotational speed and the torque also nearly represents the peak power output, as shown in Figure 26. As expected, the power output increases the inlet velocity up to a flow coefficient of ~1.0; then, the power decreases significantly with the flow coefficient. Figure 25 and Figure 26 represent a performance map that can be used to determine the power output for a range of inlet velocities, rotational speed, and flow coefficients that can be used for future optimization of the studied turbine blades.
Figure 27 shows the wall shear distribution at the 40% and 80% span locations. Both locations show similar trends where the flow attaches at the leading edge (x/c = 0). The attachment location moves downstream on the pressure side of the turbine with an increasing flow coefficient. In the cases shown, there is a sharp increase in shear stress moving towards the suction side downstream of the incipient leading-edge attach point. The increase in shear stress is reduced with increasing flow coefficient at x/c = 0.1. The boundary layer separates from the suction side until it reattaches close to the trailing edge. On the pressure side, the flow displays different behavior. For flow coefficients φ = 0.125 and φ = 0.83, the boundary layer is attached along the chord of the turbine. However, at φ = 0.063, the flow shows an increasing downstream negative shear stress that is typically representative of an adverse pressure gradient but not so much as to cause separation until very close to the trailing edge at x/c = 0.98.
Figure 28 and Figure 29 show that the flow field streamlines at the 40% and 80% span locations, respectively. The flow field is highly three-dimensional, mostly separated on the suction side. At the 40% span location, beginning at a flow coefficient of 0.42, two large vortices are seen on the suction side of the turbine. As the flow coefficient increases, the vortices convect away from the blade. We observe a similar highly separated flow at the 80% span location. The comparison of the 80% and 40% spans reveals that the vortices have already moved away from the blade at the lowest flow coefficient (φ = 0.42) and advect away from the blade with an increasing flow coefficient. Additionally, we predicted a decrease in the velocity magnitude upstream of the leading edge. This decrease leads to an increase in the static pressure, thus increasing the pressure drop.

4. Conclusions

This study aimed to investigate a novel bi-directional axial flow turbine for the OWC wave energy converter through three-dimensional CFD simulations. Modifying the Denniss–Auld turbine blades, we examined how blade geometry impacts aerodynamic performance. The results demonstrated that the modified Denniss–Auld turbine behaves similarly to impulse turbines at higher flow coefficients. Unlike Wells turbines, the modified Denniss–Auld turbine does not experience blade stall at higher flow coefficients. It maintains fairly constant efficiency across a wide range of flow coefficients, akin to the impulse turbine.
Although increasing axial inlet velocity can enhance efficiency, the modified Denniss–Auld turbine remains significantly less efficient than an impulse turbine. The proposed turbine has the advantage of requiring less mechanical complexity, as it does not need a variable pitching mechanism for the turbine blades. However, its lower efficiency suggests that further optimization is necessary.
Future work should focus on optimizing the turbine blades’ geometry, potentially incorporating active flow control to improve aerodynamic performance.

Author Contributions

Conceptualization, M.N.U., Y.G. and P.M.A.; methodology, M.N.U. and Y.G.; software, M.N.U.; validation, M.N.U., Y.G. and P.M.A.; formal analysis, M.N.U., Y.G. and P.M.A.; investigation, M.N.U.; resources, P.M.A.; data curation, M.N.U. and Y.G.; writing—original draft preparation, M.N.U.; writing—review and editing, Y.G. and P.M.A.; visualization, M.N.U.; supervision, P.M.A.; project administration, P.M.A.; funding acquisition, P.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the North Carolina Renewable Ocean Energy Program (NCROEP), grant number 111110-740106-72000-0000.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Gielen, D.; Boshell, F.; Saygin, D.; Bazilian, M.D.; Wagner, N.; Gorini, R. The Role of Renewable Energy in the Global Energy Transformation. Energy Strategy Rev. 2019, 24, 38–50. [Google Scholar] [CrossRef]
  2. Reguero, B.G.; Losada, I.J.; Méndez, F.J. A Global Wave Power Resource and Its Seasonal, Interannual and Long-Term Variability. Appl. Energy 2015, 148, 366–380. [Google Scholar] [CrossRef]
  3. Marco, T.; Pranzo, D.; Camporeale, S.; Pascazio, G. Improved Design of High Solidity Wells Turbine. In Proceedings of the The Ninth European Wave and Tidal Energy Conference (EWTEC2011), Southampton, UK, 5–9 September 2011. [Google Scholar]
  4. Uddin, M.N.; Atkinson, M.; Opoku, F. A Computational Fluid Dynamics Investigation of a Numerically Simulated Wave Tank. Am. J. Mech. Eng. 2020, 8, 40–49. [Google Scholar] [CrossRef]
  5. Opoku, F.; Uddin, M.N.; Atkinson, M. A Review of Computational Methods for Studying Oscillating Water Columns—The Navier-Stokes Based Equation Approach. Renew. Sustain. Energy Rev. 2023, 174, 113124. [Google Scholar] [CrossRef]
  6. Alkhalifa, A.S.; Uddin, M.N.; Atkinson, M. Aerodynamic Performance Analysis of Trailing Edge Serrations on a Wells Turbine. Energies 2022, 15, 9075. [Google Scholar] [CrossRef]
  7. Das, T.K.; Samad, A. Wells Turbine as a Power Take-Off Mechanism for Wave Energy Converters. In Ocean Wave Energy Systems: Hydrodynamics, Power Takeoff and Control Systems; Samad, A., Sannasiraj, S.A., Sundar, V., Halder, P., Eds.; Ocean Engineering & Oceanography; Springer International Publishing: Cham, Switerlands, 2022; pp. 365–396. ISBN 978-3-030-78716-5. [Google Scholar]
  8. Wang, X.; Yan, Y.; Wang, W.-Q. Effects of Blade Parameters on the Hydrodynamic Performance of an Impulse Turbine of Oscillation Water Column Wave Energy. Ocean Eng. 2021, 224, 108760. [Google Scholar] [CrossRef]
  9. Kumar, P.M.; Samad, A. Nature-Inspired Design of a Turbine Blade Harnessing Wave Energy. Proc. Inst. Mech. Eng. Part A J. Power Energy 2020, 234, 670–689. [Google Scholar] [CrossRef]
  10. Gato, L.M.C.; Carrelhas, A.A.D.; Cunha, A.F.A. Performance Improvement of the Axial Self-Rectifying Impulse Air-Turbine for Wave Energy Conversion by Multi-Row Guide Vanes: Design and Experimental Results. Energy Convers. Manag. 2021, 243, 114305. [Google Scholar] [CrossRef]
  11. Setoguchi, T.; Takao, M. Current Status of Self Rectifying Air Turbines for Wave Energy Conversion. Energy Convers. Manag. 2006, 47, 2382–2396. [Google Scholar] [CrossRef]
  12. Bhattacharyya, R.; McCormick, M.E. Chapter 4—Power Transfer Systems. In Elsevier Ocean Engineering Series; Bhattacharyya, R., McCormick, M.E., Eds.; Wave Energy Conversion; Elsevier: Amsterdam, The Netherlands, 2003; Volume 6, pp. 37–46. [Google Scholar]
  13. Falcão, A.F.O.; Henriques, J.C.C.; Gato, L.M.C. Self-Rectifying Air Turbines for Wave Energy Conversion: A Comparative Analysis. Renew. Sustain. Energy Rev. 2018, 91, 1231–1241. [Google Scholar] [CrossRef]
  14. Hu, Q.; Li, Y. Unsteady RANS Simulations of Wells Turbine Under Transient Flow Conditions. J. Offshore Mech. Arct. Eng. 2017, 140, 011901. [Google Scholar] [CrossRef]
  15. Halder, P.; Rhee, S.H.; Samad, A. Numerical Optimization of Wells Turbine for Wave Energy Extraction. Int. J. Nav. Archit. Ocean Eng. 2017, 9, 11–24. [Google Scholar] [CrossRef]
  16. Das, T.K.; Kumar, K.; Samad, A. Experimental Analysis of a Biplane Wells Turbine under Different Load Conditions. Energy 2020, 206, 118205. [Google Scholar] [CrossRef]
  17. Wang, H.; Li, Y.; Pröbsting, S.; Zheng, X.; Hu, Q.; Chen, D.; Liu, H.; Wu, X. Aerodynamics of the Wells Turbine with a Hawkmoth-Inspired Blade Design. Bioinspir. Biomim. 2020, 15, 066001. [Google Scholar] [CrossRef]
  18. Folley, M.; Curran, R.; Whittaker, T. Comparison of LIMPET Contra-Rotating Wells Turbine with Theoretical and Model Test Predictions. Ocean Eng. 2006, 33, 1056–1069. [Google Scholar] [CrossRef]
  19. Falcão, A.F.O.; Henriques, J.C.C. Oscillating-Water-Column Wave Energy Converters and Air Turbines: A Review. Renew. Energy 2016, 85, 1391–1424. [Google Scholar] [CrossRef]
  20. Curran, R.; Denniss, T.; Boake, C. Multidisciplinary Design for Performance: Ocean Wave Energy Conversion. Proc. Int. Offshore Polar Eng. Conf. 2000, 4, 434–441. [Google Scholar]
  21. Finnigan, T.; Auld, D. Model Testing of a Variable-Pitch Aerodynamic Turbine. In ISOPE International Ocean and Polar Engineering Conference; ISOPE: Mountain View, CA, USA, 2003. [Google Scholar]
  22. Khan, A.; Irfan, M.; Niazi, U.M.; Shah, I.; Legutko, S.; Rahman, S.; Alwadie, A.S.; Jalalah, M.; Glowacz, A.; Khan, M.K.A. Centrifugal Compressor Stall Control by the Application of Engineered Surface Roughness on Diffuser Shroud Using Numerical Simulations. Materials 2021, 14, 2033. [Google Scholar] [CrossRef] [PubMed]
  23. Teso-Fz-Betoño, D.; Juica, M.; Portal-Porras, K.; Fernandez-Gamiz, U.; Zulueta, E. Estimating the Reattachment Length by Realizing a Comparison between URANS K-Omega SST and LES WALE Models on a Symmetric Geometry. Symmetry 2021, 13, 1555. [Google Scholar] [CrossRef]
  24. Zhao, M.; Wan, D.; Gao, Y. Comparative Study of Different Turbulence Models for Cavitational Flows around NACA0012 Hydrofoil. J. Mar. Sci. Eng. 2021, 9, 742. [Google Scholar] [CrossRef]
  25. Opoku, F.; Atkinson, M.D.; Uddin, M.N. Computational Fluid Dynamics Investigation of a Novel Oscillating Water Column. In AIAA SCITECH 2022 Forum; American Institute of Aeronautics and Astronautics: San Diego, CA, USA, 2022. [Google Scholar]
  26. Opoku, F.; Atkinson, M.; Uddin, M.N. Numerical Investigation of an Offshore Oscillating Water Column. Am. J. Mech. Eng. 2020, 8, 88–105. [Google Scholar] [CrossRef]
  27. Uddin, M.N.; Atkinson, M.; Das, T.K.; Esau, S. Computational Fluid Dynamics Investigation of a Novel Bi-Directional Axial Flow Turbine for Wave Energy Conversion. In ASME Power Conference; American Society of Mechanical Engineers Digital Collection; American Society of Mechanical Engineers: Pittsburgh, PA, USA, 2022. [Google Scholar]
  28. Uddin, M.N.; Atkinson, M.D.; Opoku, F. Numerical Investigation of a Modified Denniss-Auld Turbine. In AIAA AVIATION 2022 Forum; AIAA AVIATION Forum; American Institute of Aeronautics and Astronautics: Chicago, IL, USA, 2022. [Google Scholar]
  29. Karim, S.M.S.; Lin, Y.-L.; Kaplan, M.L. Formation Mechanisms of the Mesoscale Environment Conducive to a Downslope Windstorm over the Cuyamaca Mountains Associated with Santa Ana Wind during the Cedar Fire (2003). J. Appl. Meteorol. Climatol. 2022, 61, 1797–1818. [Google Scholar] [CrossRef]
  30. Curran, R.; Gato, L.M.C. The Energy Conversion Performance of Several Types of Wells Turbine Designs. Proc. Inst. Mech. Eng. Part A J. Power Energy 1997, 211, 133–145. [Google Scholar] [CrossRef]
  31. Das, T.K.; Samad, A. Influence of Stall Fences on the Performance of Wells Turbine. Energy 2020, 194, 116864. [Google Scholar] [CrossRef]
  32. CFX-Solver, A. Theory Guide. Release Ll. 2006. Available online: https://www.academia.edu/ (accessed on 15 July 2024).
  33. Alfonsi, G. Reynolds-Averaged Navier-Stokes Equations for Turbulence Modeling. Appl. Mech. Rev. 2009, 62, 040802. [Google Scholar] [CrossRef]
  34. Bayraktar, E.; Mierka, O.; Turek, S. Benchmark Computations of 3D Laminar Flow around a Cylinder with CFX, OpenFOAM and FeatFlow. IJCSE 2012, 7, 253. [Google Scholar] [CrossRef]
  35. Halder, P.; Samad, A.; Thévenin, D. Improved Design of a Wells Turbine for Higher Operating Range. Renew. Energy 2017, 106, 122–134. [Google Scholar] [CrossRef]
  36. Ghisu, T.; Puddu, P.; Cambuli, F. A Detailed Analysis of the Unsteady Flow within a Wells Turbine. Proc. Inst. Mech. Eng. Part A J. Power Energy 2017, 231, 197–214. [Google Scholar] [CrossRef]
  37. Shaaban, S.; Abdel Hafiz, A. Effect of Duct Geometry on Wells Turbine Performance. Energy Convers. Manag. 2012, 61, 51–58. [Google Scholar] [CrossRef]
  38. Schobeiri, M.T.; Abdelfattah, S. On the Reliability of RANS and URANS Numerical Results for High-Pressure Turbine Simulations: A Benchmark Experimental and Numerical Study on Performance and Interstage Flow Behavior of High-Pressure Turbines at Design and Off-Design Conditions Using Two Different Turbine Designs. J. Turbomach. 2013, 135. [Google Scholar] [CrossRef]
Figure 4. Air turbines: Wells turbine, Denniss–Auld turbine, and impulse turbine [12]. (a) Wells turbine. (b) Denniss-Auld turbine. (c) Impulse turbine.
Figure 4. Air turbines: Wells turbine, Denniss–Auld turbine, and impulse turbine [12]. (a) Wells turbine. (b) Denniss-Auld turbine. (c) Impulse turbine.
Energies 17 03602 g004
Figure 5. (a) Denniss–Auld turbine and (b) modified Denniss–Auld turbine.
Figure 5. (a) Denniss–Auld turbine and (b) modified Denniss–Auld turbine.
Energies 17 03602 g005
Figure 6. Computational mesh generation around the blade: (a) mesh near leading edge and (b) mesh near trailing edge.
Figure 6. Computational mesh generation around the blade: (a) mesh near leading edge and (b) mesh near trailing edge.
Energies 17 03602 g006
Figure 7. Computational domain and boundary conditions.
Figure 7. Computational domain and boundary conditions.
Energies 17 03602 g007
Figure 8. Grid independence study on turbine torque coefficient at ϕ = 0.225 .
Figure 8. Grid independence study on turbine torque coefficient at ϕ = 0.225 .
Energies 17 03602 g008
Figure 9. Validation of current CFD results with experimental [30] and other CFD [31] results in non-dimensional parameters: (a) torque coefficient; (b) pressure drop coefficient; and (c) efficiency.
Figure 9. Validation of current CFD results with experimental [30] and other CFD [31] results in non-dimensional parameters: (a) torque coefficient; (b) pressure drop coefficient; and (c) efficiency.
Energies 17 03602 g009
Figure 10. Torque coefficient vs. flow coefficient for various inlet axial flow velocities.
Figure 10. Torque coefficient vs. flow coefficient for various inlet axial flow velocities.
Energies 17 03602 g010
Figure 11. Efficiency vs. flow coefficient for various inlet axial flow velocities.
Figure 11. Efficiency vs. flow coefficient for various inlet axial flow velocities.
Energies 17 03602 g011
Figure 12. Efficiency vs. number of revolutions (RPM) for various inlet axial flow velocities.
Figure 12. Efficiency vs. number of revolutions (RPM) for various inlet axial flow velocities.
Energies 17 03602 g012
Figure 13. Angle of attack at maximum lift coefficient and maximum lift coefficient vs. Reynolds number.
Figure 13. Angle of attack at maximum lift coefficient and maximum lift coefficient vs. Reynolds number.
Energies 17 03602 g013
Figure 14. Peak efficiency vs. Reynolds number.
Figure 14. Peak efficiency vs. Reynolds number.
Energies 17 03602 g014
Figure 15. Torque coefficient vs. number of revolutions (RPM) for various inlet axial flow velocities.
Figure 15. Torque coefficient vs. number of revolutions (RPM) for various inlet axial flow velocities.
Energies 17 03602 g015
Figure 16. Input coefficient vs. flow coefficient for various inlet axial flow velocities.
Figure 16. Input coefficient vs. flow coefficient for various inlet axial flow velocities.
Energies 17 03602 g016
Figure 17. Input coefficient vs. number of revolutions (RPM) for various inlet axial flow velocities.
Figure 17. Input coefficient vs. number of revolutions (RPM) for various inlet axial flow velocities.
Energies 17 03602 g017
Figure 18. Torque coefficient and input coefficient vs. Reynolds number at peak efficiency.
Figure 18. Torque coefficient and input coefficient vs. Reynolds number at peak efficiency.
Energies 17 03602 g018
Figure 19. Comparison with standard Denniss–Auld turbine [21]: (a) torque coefficient; (b) input coefficient; (c) efficiency; and (d) peak efficiency at various blade pitching angles.
Figure 19. Comparison with standard Denniss–Auld turbine [21]: (a) torque coefficient; (b) input coefficient; (c) efficiency; and (d) peak efficiency at various blade pitching angles.
Energies 17 03602 g019
Figure 20. Cp comparison at a plane through 40% from hub at different flow coefficients.
Figure 20. Cp comparison at a plane through 40% from hub at different flow coefficients.
Energies 17 03602 g020
Figure 21. Cp comparison at a plane through 80% from hub at different flow coefficients.
Figure 21. Cp comparison at a plane through 80% from hub at different flow coefficients.
Energies 17 03602 g021
Figure 22. Cp comparison at a plane through 80% and 40% from hub at flow coefficient of 0.42.
Figure 22. Cp comparison at a plane through 80% and 40% from hub at flow coefficient of 0.42.
Energies 17 03602 g022
Figure 23. Cp at a plane through 40% from hub at various flow coefficients.
Figure 23. Cp at a plane through 40% from hub at various flow coefficients.
Energies 17 03602 g023
Figure 24. Cp at a plane through 80% from hub at various flow coefficients.
Figure 24. Cp at a plane through 80% from hub at various flow coefficients.
Energies 17 03602 g024
Figure 25. Rotational speed and torque at various flow coefficients.
Figure 25. Rotational speed and torque at various flow coefficients.
Energies 17 03602 g025
Figure 26. Power and rotational speed at various flow coefficients.
Figure 26. Power and rotational speed at various flow coefficients.
Energies 17 03602 g026
Figure 27. Wall shear for various flow coefficients at 40% and 80% span locations.
Figure 27. Wall shear for various flow coefficients at 40% and 80% span locations.
Energies 17 03602 g027
Figure 28. Streamlines at a plane through 40% from hub at various flow coefficients.
Figure 28. Streamlines at a plane through 40% from hub at various flow coefficients.
Energies 17 03602 g028aEnergies 17 03602 g028b
Figure 29. Streamlines at a plane through 80% from hub at various flow coefficients.
Figure 29. Streamlines at a plane through 80% from hub at various flow coefficients.
Energies 17 03602 g029aEnergies 17 03602 g029b
Table 1. Main features of the turbine models.
Table 1. Main features of the turbine models.
FeaturesParameters
Number of blades8
Blade chord length125 mm
Hub diameter400 mm
Tip diameter600 mm
Solidity at mean radius0.64
Tip clearance1% of the chord length
Hub-to-tip ratio0.67
Aspect ratio0.75
Table 2. Meshing and appropriate boundary conditions.
Table 2. Meshing and appropriate boundary conditions.
ParameterSingle Turbine
Mesh/natureUnstructured
No. of inflation layers20
y+<1
FluidAir at 25 °C
Turbulence modelk- ω SST
InletUniform inlet velocity
OutletPressure outlet
Hub, tip, and bladeNo-slip wall
Number of revolutions2000 rpm
Convergence criterion1 × 10−5
Mass imbalance0.001
Table 3. Grid sizes.
Table 3. Grid sizes.
MeshNo. of Cells (Millions)Cell Size (mm)Turbulence ModelSimulation Time
Coarse3.610.0k- ω SST8 h
Medium7.65.014 h
Fine18.12.520 h
Extra-fine321.2526 h
Table 4. Uncertainty estimation.
Table 4. Uncertainty estimation.
ParametersNumbers
Grid size (mm)Extra-fine, fine, medium, and coarse1.25, 2.5, 5.0, 10.0
Average cell size (h)h1, h2, h3, h48 × 10−4 m, 0.0017 m, 0.0033 m, 0.0057 m
Grid refinement factor (r)r21, r32, r432.07, 1.9681, 1.7318
Performance parameter ( ϕ ) ϕ 1 , ϕ 2 ,   ϕ 3 ,   ϕ 4 0.1273, 0.1267, 0.1216, 0.0925
Apparent orderp2.83
Approximate relative error e a 21 ,   e a 32 ,   e a 43 0.0050, 0.0403, 0.2393
Extrapolated values ϕ e x t 21 ,   ϕ e x t 32 ,   ϕ e x t 43 0.1279, 0.1319, 0.1613
Extrapolated relative error e e x t 21 ,   e e x t 32 ,   e e x t 43 0.0096, 0.0787, 0.4268
GCI G C I e x t r a f i n e 21 ,   G C I f i n e 32 ,
G C I m e d i u m 43
0.0058, 0.0520, 0.4088
Table 5. Inlet axial velocity and flow coefficient ranges.
Table 5. Inlet axial velocity and flow coefficient ranges.
Flow   Coefficient ,   ϕ Inlet Axial Velocity, v (m/s)Reynolds Number, Re
0.0754.715.25 × 105
0.1257.855.27 × 105
0.17510.995.31 × 105
0.2012.945.34 × 105
0.22514.145.36 × 105
0.2515.525.39 × 105
0.27517.285.43 × 105
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Uddin, M.N.; Gao, Y.; Akangah, P.M. Computational Fluid Dynamics Simulation on Blade Geometry of Novel Axial FlowTurbine for Wave Energy Extraction. Energies 2024, 17, 3602. https://doi.org/10.3390/en17143602

AMA Style

Uddin MN, Gao Y, Akangah PM. Computational Fluid Dynamics Simulation on Blade Geometry of Novel Axial FlowTurbine for Wave Energy Extraction. Energies. 2024; 17(14):3602. https://doi.org/10.3390/en17143602

Chicago/Turabian Style

Uddin, Mohammad Nasim, Yang Gao, and Paul M. Akangah. 2024. "Computational Fluid Dynamics Simulation on Blade Geometry of Novel Axial FlowTurbine for Wave Energy Extraction" Energies 17, no. 14: 3602. https://doi.org/10.3390/en17143602

APA Style

Uddin, M. N., Gao, Y., & Akangah, P. M. (2024). Computational Fluid Dynamics Simulation on Blade Geometry of Novel Axial FlowTurbine for Wave Energy Extraction. Energies, 17(14), 3602. https://doi.org/10.3390/en17143602

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