A Numerical Investigation of the Energy Efﬁciency Enhancement of Oscillating Water Column Wave Energy Converter Systems

: This work focuses on the geometry effects over the performance of oscillating water column (OWC)-type wave energy converter (WEC) systems and searches for the OWC geometries that enhance the energy efﬁciencies under the same wave conditions. To analyze the hydrodynamic performances of the WEC systems, an in-house smoothed particle hydrodynamics (SPH) code based on weakly compressible ﬂuid approach is utilized. The energy efﬁciency enhancement studies of the determined OWC device are carried out with a two-step geometry modiﬁcation procedure. The ﬁrst step starts with the validation of the free-surface elevation and orbital velocity time histories. Then, a three-by-three simulation matrix that depends on the geometrical design parameters of chamber length and front wall draft is run at three different wave conditions, and the OWC geometry that produces the maximum energy efﬁciency is determined. In the second step, the corner regions of the obtained optimal geometry are chamfered, and another simulation matrix is tested at the wave condition that yields maximum wave energy. It is observed in this step that the energy efﬁciency index can still be improved by 4.3% by only chamfering the back face of the OWC chamber. To scrutinize the physical grounds of this increase, the correlation between the time-averaged vorticity and energy efﬁciency is presented. Finally, the performance of the best conﬁguration is also examined in three different wave periods, where the suggested geometry shows better performance with respect to base geometry results in all wave conditions.


Introduction
Because of global warming due to the dramatic increase in CO 2 levels mainly associated with the massive fossil fuel consumption, energy production relying on renewable sources (i.e., solar, wind, and wave energies) has increased tremendously in the last two decades [1,2]. In 2017, it was reported that almost 20% of the total energy generation of humankind is based upon the renewable energy sources [3]. Among renewable energy resources, the high theoretical potential of the ocean energy [4] makes free-surface waves a prominent candidate for generating green electricity [5].
The most prominent barrier for the commercialization of wave energy converters (WECs) is the high initial capital and operational expenses associated with harvesting energy from waves in hazardous offshore and nearshore conditions [6]. Thus far, multiple efforts have been made to solve these hurdles by introducing various WEC concepts. To improve the efficiency of the WECs, various concepts and designs have been developed, where the number of patents has reached more than a thousand as of 2002 [7]. Despite the considerable variation in the suggested concepts, WECs are generally categorized by types of harvesting technology such as oscillating water columns, oscillating wave surge converters, submerged pressure differential, and overtopping devices [5]. These systems can be placed along the shoreline as stationary structures or employed as floating structures on the open seas. General overviews of available concepts and designs for wave energy converters can be found in [8][9][10][11].
The oscillating water column (OWC) is one of the most promising concepts that have reached the prototype development stage in some regions, e.g., the LIMPET experimental plant in Islay Island, Scotland [12], the Mutriku power plant in Spain [13], the Pico power plant in Portugal, and the REWEC3 experimental plant in Italy [12]. OWC systems consist of a hollow and partially immersed storage that is open to water with an air-trapped chamber. Ocean waves cause the free surface within the tank to oscillate vertically, which pressurizes and depressurizes the air inside the chamber. The pneumatic force of this trapped air enters and leaves the bidirectional turbine located on the device chamber and turns into electrical energy through power take-off systems [14].
In recent years, several experimental and numerical studies have been conducted to investigate and characterize the hydrodynamic behavior of WECs from different aspects. Morris-Tomas et al. [15] performed experiments to study the effect of the front wall configuration on the hydrodynamic efficiency of OWCs. Chang et al. [16] conducted an experimental investigation to determine the best geometrical design for the OWC device to achieve optimum hydrodynamic efficiency. Various geometric parameters were tested. The results showed that chamber geometry, especially the backplate slope, substantially influenced the OWC performance.
Although experimental methods capture and reveal such devices' hydrodynamic behavior more realistically, they are frequently regarded to be costly and time-consuming, and may also suffer from the scale effects at the design stage. As an alternative to experimental approaches, computational investigation of the hydrodynamic performance of OWCs using numerical wave tanks is a reliable and less expensive method. Simonetti et al. [17] used a mesh-based finite volume method in the OpenFOAM environment to analyze a fixed asymmetric OWC device. The results showed that chamber length, front wall draft, and relative OWC PTO damping had a strong effect on OWC efficiency. As shown in [18], chamber length has a substantial influence on OWC performance when compared to other OWC geometric characteristics. In addition, the angle of inclination in the bottom profile is an important parameter for shoreline OWC efficiency. Kamath et al. [19] investigated the hydrodynamic behavior of a 2D OWC under various wave conditions; they used the REEF3D CFD code in their numerical study. In the simulation, lower wave steepness caused the inner free surface to move firmly similar to a piston, but higher wave steepness caused nonuniform motion.
One of the important engineering problems that can significantly affect the performance of OWC devices and cause strong nonlinearity behavior inside OWC chambers is the resonant phenomenon (gap and harbor resonance). The modifications of the OWC geometry and topographical slope can alter the resonant frequency of the water body inside the chamber. Gao et al. [20] used OpenFOAM to study the impact of topographical variation on the fluid gap resonance between a stationary box and a vertical wall using incident regular waves. They discovered that the fluid resonant frequency decreased monotonically as the topographical slope increased. Gao et al. [21] studied the effects of Bragg reflection on harbor oscillations, which are caused by two types of incoming steady-state waves. During their simulations, they thoroughly examined the impact of topographic characteristics on oscillations.
In addition to widely implemented mesh-based computational solutions [19,[22][23][24], which may have difficulties in capturing the nonlinear dynamics of high free-surface deformations in the problem area, and the instantaneous response of WEC devices in severe wave conditions, in the last decade, meshless techniques have attained significant interest in providing high-fidelity results [25][26][27][28]. Luo et al. [29] reviewed the state of the art in the application of meshless techniques to solve hydrodynamic issues in marine and coastal engineering such as wave formation, propagation, and refraction, as well as related turbulence generation and dissipation, air entrainment, and mass transfer. Fu et al. [26] used a semi-Lagrangian meshless frame to investigate the sloshing phenomenon in a two-dimensional numerical tank. Their approach was shown to provide more stable and accurate results than one with an explicit mesh-based scheme.
Smoothed particle hydrodynamics (SPH) is one of the most widely utilized meshfree computational approaches for modeling diverse physical fluid flow conditions [30][31][32][33][34]. He et al. [32] developed a coupled weakly compressible and total Lagrangian SPH (WC-TL SPH) approach to model the interactions of elastic entities with free-surface flows. Khayyer et al. [35] briefly reviewed recent progress in the development of totally Lagrangian, mesh-free computational techniques for hydroelastic fluid-structure interactions in marine engineering and highlighted some important issues in this context. Almasi et al. [30] presented a multiphase incompressible smoothed particle hydrodynamics (ISPH) method to examine complicated multiphysics electrohydrodynamics (EHD) issues. Furthermore, the SPH technique is commonly employed to research hydrodynamic problems because of its particular benefits in modeling highly nonlinear free-surface problems [36][37][38][39]. Gotoh et el. [40] provided an update on recent advances in particle methods. They discussed recent advances in accuracy, fluid-structure interactions, and computational efficiency. Lyu et al. [41] provided a detailed overview of SPH-based hydrodynamic simulations for ocean energy systems. The main focus is on three topics: SPH-based numerical fluid tanks, multiphysics SPH methods for modeling ocean energy systems, and, finally, computing efficiency and capacity. Ozbulut et al. [36] utilized the SPH technique in order to model a numerical wave tank. A wide variety of test cases with different wavelengths and steepness ratios for regular and irregular ocean waves were investigated using piston-and flap-type wavemakers.
The application of the SPH method to model the wave energy converters has also increased over the past few years [38,[42][43][44]. Crespo et al. [42] simulated wave interaction with an offshore OWC device in a numerical wave tank by using an open-source SPH code (DualSPHysics). It was proven that their models can accurately create the free-surface evolution inside the chamber when the air pressure is neglected. Wen et al. [44] used an improved SPH model to examine the hydrodynamic performance of onshore OWC and considered the turbulence effect during their work. It was shown that a sloshing phenomena could be detected in the OWC chamber when the front wall depth was reduced. More recently, Ropero-Giralda et al. [38] utilized the SPH method to investigate the efficiency and survivability of a point-absorber-type WEC device under different regular wave conditions. Recently, Quartier et al. [38] used the SPH approach to numerically model an OWC WEC. In order to numerically simulate the power take-off (PTO) system, a force is applied to a plate that is floating on the free surface of the OWC chamber. This eliminates the air phase simulation, which is computationally costly in SPH approaches.

Objectives of This Paper
In this study, the weakly compressible SPH (WCSPH) approximation [45] is utilized to investigate the hydrodynamic characteristics and the energy efficiencies of OWC-type wave energy converter system. The current geometrical model of the present study is assumed to be the simplified version of the 3D case due to its longitudinal symmetry, where the validity of the results can be generalized for the 3D cases by ignoring the effects of higher-order terms related to the eddy-vortices. The proposed in-house SPH code is validated by an experimental study [46] consisting of a stationary OWC with an open chamber. Additionally, a particle resolution and convergence test is also performed to determine the optimum particle resolution on the problem domain. After achieving accurate results by utilizing the proposed SPH scheme on an OWC system, a two-step efficiency enhancement study is performed based on the geometrical modifications of the chamber part in the OWC device to obtain the configuration that provides maximum available wave energy. In the first step, the geometrical parameters (including length of chamber and front wall draft) are varied to obtain the maximum harvested energy in three different wave periods. Using the best geometrical configuration that leads to maximum wave energy production among all wave period conditions, the second step is applied by chamfering the corner parts of the chamber region, considering all possible configurations. As a final scrutiny, an investigation to obtain the relationship between the time-averaged magnitude of vorticity in the OWC inner channel and the hydrodynamic efficiency is performed through a quantitative comparison among all test cases.
The novelty of this study lies in the systematic and comprehensive analysis of the OWC system geometric configurations to increase the amount of extracted energy under equal wave conditions. It can be said as one of the major outcomes of this study that the proposed SPH scheme can accurately capture the energy level changes even after applying slight modifications to the chamber corner geometries, where the correlation between vorticity magnitudes in the flow domain and harvested energy is verified in all geometry configurations. In conclusion, due to the highly robust and accurate capabilities of the utilized in-house code of the proposed work, it becomes possible to claim that minor modifications to the chamber geometry may further improve the overall efficiency of the OWC system under the same wave conditions. The paper's organization is as follows: In the second section, the governing equations, and the discretized form of these equations based on the proposed WCSPH approach, are mentioned. Following the definition of the physical and geometrical parameters of the problem, hydrodynamic efficiency equations of the OWC device are presented in the third section. In the fourth section, the robustness, consistency, and accuracy of our numerical method are examined through comparing the experimental and theoretical results found in the literature [46]. Then, the effect of geometrical modifications on the OWC system is scrutinized with quantitative and qualitative measures in the fifth section. Finally, the discussions on the obtained results and concluding remarks are highlighted in the last section.

Governing Equations and SPH Modeling
In the present study, the continuity and conservation of linear momentum equations are employed to govern the motions of barotropic and isothermal fluid domain: where D/Dt is the material time derivative, u represents the velocity, ρ denotes the density, p is the pressure, µ is the dynamic viscosity, and g is the gravitational acceleration. In passing, it would be explanatory to state that turbulent fluctuations and associated dissipation within the computational domain are prudently ignored with the assumption that the smallest length-scale of turbulent movement does not affect the free-surface deformation, notably. The governing equations are discretized by the well-known WCSPH approach [47] which has been widely utilized in many engineering applications with free-surface hydrodynamics modeling [38,39,48,49]. To close the governing equation system, the WCSPH approach utilizes an equation of state that couples pressure and density variations. The following equation of state is used during the simulations of this study: where c 0 is the reference speed of sound, ρ 0 is the reference density for water equal to 1000 kg/m 3 , and γ is the specific heat ratio of water, taken as 7. The reference speed of sound is determined by keeping the Mach number (M) below 0.1, which satisfies the density variations less than 1% [50]. To guarantee the limited density variations among all simulated wave conditions of the present study, the value of c 0 is set at 25 times the maximum wave group velocity, which gives 200 m/s, and satisfies the incompressibility condition, especially in the regions close to the moving boundary and beach. It should be remembered here that assigning the physical value of speed of sound in water leads to undesirable time step values in terms of high computational costs and also leads to a negligible contribution in the enforcement of incompressibility constraints.
In the SPH method, the fluid domain is represented by a collection of free particles, all of which contain information about field variables such as density, velocity, and pressure. The SPH formulation relies essentially on the integral representation of a function, which is discretizable using a particle approximation [51]. Given particle i, the values of an arbitrary field variable, which is either a scalar f i or a vector f s i , are computed by an interpolation procedure that uses the amounts of the neighboring particles, represented as f j or f s j . To clarify the notation used in the rest of the article, note that mixed notation (both direct and index notations) is used for either the vector or tensor fields. In direct notation, vectors are denoted by bold lowercase letters, while tensor fields are denoted by bold uppercase letters. In index notation, italic Latin indices are used to denote vector or tensor components, and they are always superscripted unless otherwise indicated. Einstein's summation convention applies, where repeating indices imply summation over them. According to the SPH technique, concisely presented here, the value of any field function is computed by the following procedure: where subscripts i and j represent the particle index, N is equal to the total number of neighbor particles for particle i. W ij , or in full form W ij (r ij , h), is the piecewise quintic kernel function [27], which depends on the dimensionless parameter q = ||r i − r j || 2 /h, h is the smoothing length equivalent to 1.33 times the initial particle distance (δx), and α d is a coefficient depending on the size of the problem domain. In a two-dimension domain, α d is equal to 7/(478πh 2 ). V j denotes the particle volume, which is computed as Before showing the discretization methodology of the governing equations through WCSPH approximation, it will be convenient to define the applied treatments for function derivatives. To enhance the accuracy of gradient and divergence operations for any scalar or vector valued field variables, a corrective SPH (CSPH) formulation [52,53] is included into the numerical scheme of this work: In Equation (8), α sl i is the second rank correction tensor, which helps to avoid potential inconsistencies such as particle clustering and truncated support domains at boundaries. By representing the inverse of α sl i as B pl i , if both sides of Equation (8) are multiplied by B pl i , the gradient correction of kernel matrix can be found as B pl i ∂W ij /∂x l i , and this can be described as B i · ∇ i W ij in direct notation.
By imposing these CSPH treatments, the described governing equations of fluid motion can be discretized through WCSPH approximation: ∇ i is the differentiation taken with respect to the position of particle i. In Equation (10), K = 2(n + 2), where n is the dimension of the problem domain [54] and ν is the kinematic viscosity, which is taken as 10 −6 (m 2 /s) for water.
In the weakly compressible SPH method, accurate density computation is very critical because the pressure values are linked to the density values through the equation of state, and small density disturbances can cause large fluctuations in the pressure field, which in turn reduces the precision and stability of numerical approach [55]. To mitigate the adverse effects of oscillatory pressure field distributions during evolution of fluid flow, the density smoothing treatment, which is frequently used in the SPH literature [56,57], is incorporated into the present study's numerical scheme: where ρ i represents the corrected density and is computed by interpolation between neighboring particle densities. As a final correction algorithm, hybrid velocity-variance-based free surface (VFS) and artificial particle displacement (APD) treatment is employed in the simulations of this study. The VFS algorithm is applied only to free-surface particles, which aims to ensure an averaged velocity and avoid excessive scattering on the free surface. The free-surface particles are defined as the particles with a neighboring particle number less than 65% of the average number of neighboring particles in the problem area. The velocities of the free-surface particles are calculated as follows: where u i indicates the corrected particle velocity, and is a dimensionless constant which was suggested to be taken between 0.05-0.1 times the initial particle distance (δx) in our previous studies [36,37,55], and accordingly it is taken as 0.075 δx in all simulations of this study. On the other hand, the APD algorithm is used for the fully populated regions of the fluid domain which considerably improves the precision of the interpolation processes of the proposed SPH method by providing a uniform particle distribution throughout flow evolution. The formulation of the APD algorithm is defined as [55]: where r i denotes the updated position vector of the particles that are inside the supported domain.
In this study, the wall particle velocities are set to the velocity of the boundary and no additional formulation is applied to fluid particles that are close to the wall boundaries. The wall boundaries are represented by fixed four-layer solid particles. Due to the inherent advantages of the SPH approach in the modeling of free-surface problems, the kinematic and dynamic free-surface boundary conditions are captured inherently without applying any treatment. The time integration process is carried through a predictor-corrector scheme, which is defined in detail in our previous study [27], and the time-step value defined by the Courant-Freidrichs-Lewy (CFL) condition, described as: where C CFL = 0.4, and c w is the produced wave system's celerity [36].

Numerical Wave Tank
In this work, the energy efficiency of the OWC-type wave energy converter device is examined through generating regular waves by using flap-type moving boundaries. The simulation cases are chosen to provide intermediate water conditions that avoid shallow water effects in the wave tank and the high computational costs associated with deep water conditions. In addition, all test cases are simulated at a constant depth to reduce the problem parameters in the simulation matrix and focus directly on the geometric parameters of the OWC chamber. Finally, the value of the wave height is chosen to produce a comparable variation of the chamber pre-draft (y) heights. Considering all of these design criteria, the wave properties of this study remain at the upper limit of second-order Stokes waves according to the Le Méhauté [58] abacus given in Figure 1. The relation between wave amplitude (ζ a ) and full stroke of the flap (S 0 = 2S) is obtained through the Biesel transfer function [59]: where k represents the wave number and d is the initial water depth at the flap boundary.
To avoid the sudden initial movement of the wavemaker, the wavemaker is set into a motion with gradually increasing displacement through utilizing a simple mathematical relation [36]: where x 0 and z 0 represent the primary positions of the wavemaker particles. The change in the vertical locations of the wavemaker particles are deemed trivial and are not updated. In Figure 2, two different two-dimensional numerical wave tanks are schematically presented. The first wave tank (top) is employed for generating the numerical results that are compared with the experimental study of Iturrioz et al. [46], while the second one (bottom) is utilized for producing simulation results for all oscillating wave energy converter test cases. The Stokes second-order regular waves with a period of 3.2 s and wave height of 0.08 m are generated in the simulations when comparing the SPH results with the experimental measurements provided in [46]. The geometrical details of the baseline OWC device are described as follows: The overall length of the channel is L 3 = 120 m, the distance between the flap and the OWC front wall is L 2 = 110 m, and the distance between the flap and the toe of the sloping wall is L 1 = 102.5 m. The still-water depth at the flap boundary (d) is equal to 10 m. θ is the 4/10 inclined beach's angle, and w and y represent the optimization variables, namely, the chamber length and draft of the front wall, respectively.

OWC Hydrodynamic Efficiency
The present study mainly concentrates on the investigation of increasing the harvested total wave energy in OWC systems by modifying the chamber geometries. Therefore, the energy losses due to the mechanical efficiency of turbine and energy conversion efficiency of the PTO system are not considered. To pursue an objective comparison between the efficiency levels of different OWC geometries, an energy efficiency index formulation is identified according to the ratio of calculated mean pneumatic power absorbed by the OWC (without any contribution from the damping effects of PTO system) to the incident wave power per length multiplied by the chamber width b: where the chamber width is taken as unity. The average pneumatic power harvested by the OWC can be expressed as follows: where T s indicates the total duration of the simulation, P(t) air represents the instantaneous differential air column pressure obtained from numerical simulation, Q t is the air volume flux inside the OWC chamber, and dt represents the time step of the simulation: In Equation (21), w is the length of the chamber and v(t) is the instantaneous vertical velocity of the free surface in the OWC chamber, which is calculated by averaging the vertical velocity of the free-surface particles inside the OWC chamber at each time step.
The relationship between the air pressure and air flow rate at the turbine inlet can be expressed with the following parabolic expression [60][61][62]: where q t is the volume flow rate of air per second per unit width (q t = Q t /b), and k dm is the damping term in the OWC chamber, which is defined as: Here, ρ air denotes the air density (taken as 1.225 kg/m 3 ) and B * is a dimensionless damping coefficient representing the relation between the pressure drop in the OWC chamber and the turbine's flow rate. The damping coefficient depends on the utilized turbine characteristics which may significantly affect the overall performance of the OWC device [60]. To keep the simplicity and for the sake of quantitative comparison in the representation of the energy efficiencies of each OWC geometries, B * is considered as unity in all simulation cases. One can assign the specific value of B * to determine the total power harvested by the OWC for the given wave conditions and OWC device geometries.
The average incident wave power per unit width for the Stokes second-order wave is indicated as [59]: where H i is the incident wave height, ω denotes the angular frequency of the incident wave, and d is the depth of water. By substituting all of the expressions into the OWC energy efficiency index, it can be simply calculated as follows:

Verification and Validation Studies for the OWC Simulations
In this section, the verification and validation of the proposed numerical scheme is realized through comparing the free-surface deformations and velocity components within the OWC chamber with the available theoretical and experimental results of the literature. In the former test case, the time series of the free-surface elevation and velocity components close to the free surface are verified with the theoretical results of second-order Stokes waves. For the latter case, the numerical performance of the proposed SPH algorithm is validated through the comparison of free-surface elevation in the center of the OWC device with experimental and time-domain results of Iturrioz et al. [46]. Additionally, a resolution test was also performed to achieve the convergence of the numerical scheme and obtain the optimum number of particle distribution in the fluid domain.
To verify the characteristics of generated waves in the numerical OWC tank, the wave kinematic quantities, namely, the free-surface elevation and orbital velocities at x = 55 m and z = −3 m, are compared with the Stokes second-order wave theory [59]. Considering the second-order wave theory, the theoretical time series of water surface elevation and orbital velocities of water particles for a defined position (x and z) are calculated as follows [59]: As can be seen from the time series given in Figure 4, the generated waves in the numerical wave tank are in good agreement with the theoretical wave characteristics.  The results of the numerical approach proposed in this study are validated and verified with the experimental and theoretical findings of Iturrioz et al. [46] in terms of the free-surface time series plotted for the center location of the OWC device given. Iturrioz et al. [46] presented a simplified time-domain model for a stationary OWC. In this model, they used the Cummins integro-differential equation to express the motion of a floating body.
To examine the particle size independence of the proposed algorithm, a convergence study is also carried out with five different particle resolutions. Furthermore, the accuracy of the numerical simulation results is measured using the root mean square error (RMSE) to make a quantitative comparison between the performances of particle resolutions, defined as: where ζ n,i and ζ e,i represent the water surface elevation within the OWC chamber for any i th sample of the numerical and experimental results, respectively, and N is the number of samples for each simulation. Figure 5 displays the comparative graph of the numerical simulations with variable particle resolutions, while Table 1 indicates the RMSE values, computational performances of each cases. One can say that SPH results of the present study can capture the overall wave characteristics inside the OWC chamber in all resolutions and produces compatible free-surface profiles with experimental results with the increasing particle resolutions. It can also be stated that SPH predicts the nonlinear characteristics of the wave trough deformations inside the chamber better than the time-domain model solution.
In light of the obtained RMSE values, it can be said that there is no significant difference between d/40 and d/50 resolutions. On the other hand, using d/60 resolutions results in a significant increase in computational costs. The d/40 resolution is still accurate enough to compare with experimental data without increasing the particle number and computation time. Following the observations and outcomes of the convergence study, d/40 particle resolution is utilized in all simulations of this work.

Results and Discussion
This section presents the simulation results for the OWC-type wave energy converter device. Initially, a simulation matrix is generated by altering the chamber length and front wall draft values of the OWC system for three distinct wave conditions, leading to twenty-seven different configurations. Then, the configuration with the largest EEI value is chosen as the reference or baseline case for which the effect of corner chamfering on the wave energy output capacity of the OWC device is investigated.

Performance Assessment of Chamber Length and Front Wall Draft Configurations
Recall that investigating the impact of geometrical combinations and wave conditions on the efficiency of the stationary OWC-type wave energy converters constitutes the main objective of this study. In this context, a three-by-three simulation matrix was created based on the variation of the characteristic lengths of the system, namely, chamber length (w) and front wall draft (y). Additionally, whole simulation matrix was run for three different wave period conditions at each geometrical configuration. It may be convenient to recall that second-order Stokes waves were generated in all OWC device simulations and wave heights were set to a constant value of 2.5 m to reduce the complexity of the comparison. In all cases, the simulations were carried out for 150 s.
The obtained EEI results for all geometry configurations and wave conditions are presented in Table 2. It can be seen that the anterior wall immersion depth may have a significant effect on OWC efficiency for the wave period T = 6 s where the EEI values drop dramatically with the increasing draft values. In the cases of wave period T = 8 s, the EEI values highly depend on both geometry parameters, where it can be deduced that there is a tendency of decreasing efficiency with higher draft sizes. In contrast, in the simulation results of the wave conditions with T = 10 s, the increase in chamber length and draft values generally leads to a positive effect on the wave energy efficiencies except the maximum values of both parameters. The comparative analysis on the results of all test cases indicates that the amount of the harvested wave energy highly depends on the OWC geometry for each particular wave condition and the collected energy may be increased by the appropriate design of OWC geometry according to the dominating wave characteristics of the region. If the utilized regular wave periods are considered, it can be clearly inferred that these OWC geometries can be a good candidate for a coastal region with the peak wave period of T = 8 s and average wave height of H = 2.5 m, where the maximum EEI values are achieved for nearly all geometrical configurations. To provide further details for the EEI calculations of all geometries, the time series of free-surface elevation, velocity, and air pressure inside the chamber are plotted in Figure 6. As can be seen in Figure 6, an increase in the front wall draft of the OWC for a fixed chamber length generally results in a decrease in the EEI index due to a decrease in the surface velocity of the water and the air pressure in the chamber. For a fixed front wall draft, a decrease in chamber length generally results in an increase in surface velocity and air pressure in the OWC chamber, leading to an increase in EEI index.

Effect of the Corner Chamfering on the Energy Effeciency
Having determined the chamber geometry configuration and wave condition that yields the maximum EEI value, a systematical investigation of this chamber geometry and wave characteristic was performed to seek for further improvement in the total efficiency by only chamfering the corners inside the OWC channel. Thus, seven additional geometry configurations (see Figure 3) were created, wherefore the dimensional details are described in Table 3 together with the best geometry of the previous analyses. As can be seen in Table 3, all simulation cases except for cases 3 and 7 cause a drop in the OWC efficiency compared to case 1 (base model). Modifying the channel walls' geometry in cases 3 and 7 leads to an increase in efficiency by 3.4 and 4.3 percent, respectively. The possible reason that lies behind this improvement is considered to be related to the amount of vorticity levels. To extract the vortical characteristics in all geometry configurations, a strategy for the objective quantification of vorticities is proposed. According to this strategy, a correlation between the amount of vorticity generated inside the common region for all geometries (shown as blue in Figure 3) and the EEI was established. The absolute amount of dimensionless vorticity (ω * = ω(λ/g) 0.5 with λ denoting the wavelength) is averaged over the total simulation time. To extract this correlation and represent the results in a succinct manner, the dimensionless vorticity and EEI values of each test cases are plotted in sequential order in Figure 7. It is numerically proven that EEI values are improving with the decrease in vorticity levels, which is actually an achievement of proper chamfers applied in the OWC chamber geometry.  Table 3. Figure 8 presents the vorticity fields of test cases with maximum and minimum EEI values whereby one can clearly see that the magnitude of the vorticity field is smaller for the case with the maximum EEI value. This observation can be attributed to the fact that the vorticity inherently leads to energy dissipation. The chamfering hampers the formation of vorticity through providing a streamlined flow pattern and larger space for fluid flow. The comparison of test cases 3 and 7 reveals that they have nearly identical EEI and vorticity values, which suggests that the chamfer formed by the dimensional factors of A3-B3 pairs does not have a notable influence on the EEI and vorticity magnitude. In all simulations, it is observed that the existence of the chamfer defined by the geometrical factor of the A1-B1 pair reduces the EEI value.
One can easily observe that there is an inverse relationship between the energy efficiency and the amount of vorticity inside the determined region. Using the chamfer configuration of cases 3 and 7 results in lower vortex flow characteristics in the blue region (Figure 3), and hence leads to lower energy losses in the OWC system. Particularly, the chamfer on the back wall of the chamber encourages more streamlined flow along the vertical direction and enhances the free-surface velocity and the air pressure inside the chamber which can be seen in the time-series graphs given in Figure 9.    The working performance of the best geometrical configuration (test case 3) was also evaluated for the wave periods of T = 6, 8, and 10 s. Achieved efficiencies are tabulated and compared with the findings of base geometry in Table 4. As can be seen from the table, modifying the channel walls' geometry to case number 3 leads to an increase in efficiency in all wave conditions, where the amount of increase depends on the incident wave characteristics. Figure 10 displays the time series of free-surface elevations, velocities, and air pressure in the chamber for all simulated wave periods.
When considering the whole OWC device simulation results, it can be claimed that utilizing the obtained chamber length and draft dimensions together with implying the chamfer modifications of test case 3 can be a good candidate for a coastal region with the peak wave period of T = 8 s and characteristic wave height of H = 2.5 m.  Figure 10. Time series of free-surface elevation (top), velocity (center), and air pressure drop (bottom) inside the OWC chamber for the best performance geometry configuration (case 3) at each wave period.

Conclusions
The search for possible energy efficiency enhancement in OWC-type WECs under the same wave conditions comprises the main objective of this study. A weakly compressible SPH approach was employed to investigate the hydrodynamic characteristics and the energy efficiency of an OWC-type wave energy converter system. Numerical simulations were performed in a two-dimensional numerical wave tank that generates regular waves using a flap-type moving boundary.
As an initial validation study, the numerical performance of the proposed SPH scheme was examined by comparing the free-surface elevation in the middle of the OWC device with experimental and time-domain method results of Iturrioz et al. [46]. To investigate the particle size independence of the numerical scheme and obtain the optimum particle size, a convergence study was conducted, where the accuracy of the numerical results was assessed by comparing the root mean square error with experimental results. The outcomes of these simulations proved that the proposed SPH scheme can capture the free-surface time series given by the experimental and time-domain method study of [46]. Secondly, the free-surface elevation and the orbital velocity time series of the proposed SPH scheme were compared with the theoretical second-order Stokes wave theory findings, where highly compatible results were observed.
Following the validation of the numerical method for the utilized OWC geometry, a two-step efficiency enhancement study was conducted through systematical geometry modifications. In the first step, the chamber length and front wall draft were varied to obtain the best geometry that leads to maximum harvested energy in three different wave periods. After determining the geometrical configuration that gave the maximum efficiency, the second-step geometry modification was carried out on the optimal geometry of the first step. A simulation matrix that consisted of the combinations of chamfered corner regions of the OWC chamber was run and it was observed that the energy efficiency of the OWC geometry can be still improved by 4.3%. Furthermore, the enhanced OWC geometry configuration also harvested more wave energy in different wave conditions with respect to the base geometry, which had no chamfering at the corners. As a final investigation, a quantified relationship between the time-averaged vorticity inside the chamber region and the energy efficiency was built. The vorticity analyses inside the chamber region indicated that reducing the vorticity magnitudes results in the higher total energy efficiency of the device, which is an expected outcome of energy conservation laws in physics.
To summarize and make an overall assessment of all test cases, the potential of the harvested wave energy highly depends on the OWC geometry for each particular wave condition, and the collected energy may be increased by the appropriate design of OWC geometry according to the dominating wave climates of the region. Nonetheless, the inspection of the effect of the water depth is not included in the limelight of this study. As trials at different water depths will lead to different oscillation characteristics of the free surface inside the chamber, it may be beneficial to extend the comprehensiveness of this study to various water depths and wave characteristics in future studies.