A numerical model for the optimization of the performance of an innovative overtopping breakwater for wave energy conversion is proposed. The model is based on the stochastic description of the overtopping phenomenon based on the results of extensive laboratory tests, and it is able to simulate the behavior of the device operating under any assigned sequence of sea states, thus allowing to easily obtain results that would otherwise require time consuming and costly physical model tests. The model is used here to identify the main geometrical parameters affecting the performance of the device and to optimize such parameters in order to maximize the average yearly output power. An application to a device embedded in the breakwater of Pantelleria Port (Sicily, Italy) is presented. The model is also proved to be useful to verify the possibility of further increasing the output power through the implementation of specific control strategies concerning the operation of the turbines. The work provides a better understanding of the influence that such a system could have on the energy system of small Mediterranean islands, for example in terms of contribution to the CO2 emission reduction.

The choice of the values of the friction parameters may strongly influence the numerical modeling of the interaction between waves and porous media. Here, an assessment of such an influence is carried out using the OpenFOAM solver IhFoamV1 to simulate the response of the Catania harbor breakwater under extreme wave attack. The numerical model was validated by comparison with an experimental dataset, and a sensitivity analysis of the overtopping discharge estimate to van Gent parameter β was carried out testing values suggested by previous studies. A discussion on the importance of a careful estimate of such a parameter when dealing with the numerical modeling of porous coastal structures is presented. Indeed, variations in the non-dimensional overtopping discharge higher than 150% were observed as a consequence of a small variation (10–20%) in the absolute value of β.