Experimental Validation of the Capture Chamber Model in Mutriku
MOWC Wave Power Plant
AITOR J. GARRIDO, SALVADOR CAYUELA, AMPARO VILLASANTE, IZASKUN GARRIDO
Automatic Control Group - ACG,
Inst. of Research and Development of Processes - IIDP,
Faculty of Engineering of Bilbao,
University of the Basque Country (UPV/EHU),
Pº Rafael Moreno 3, Bilbao, 48013,
SPAIN
Abstract: - Wave energy holds the potential to fulfill 15% of the EU's energy demand by 2050, thereby
reducing CO2 emissions by 136 million metric tons per megawatt-hour, as outlined in the EU Energy Road
Map. Similarly, the Spanish Renewable Energies Plan underscores the significant marine energy potential in
Spain, particularly emphasizing wave energy. Within this framework, Oscillating Water Column (OWC)
converters currently stand as among the most promising wave energy conversion technologies, offering the
capability to harness ocean energy from various on-shore and floating structures. This paper introduces an
analytical model of the wave capture chamber parameterized for a specific on-shore OWC wave power plant.
The model is specifically adapted and parameterized for the Mutriku Marine Offshore Wave Power Plant
located on the coast of the Spanish Basque Country. Subsequently, validation is conducted using both real wave
entry data measured on-site and experimental output power data generated in the plant.
Key-Words: - Wave Energy, Energy Converter, Modeling, Data Collection, Empirical Verification, Validation,
Marine Energy, Renewable Energy.
Received: April 25, 2023. Revised: February 23, 2024. Accepted: March 29, 2024. Published: May 29, 2024.
1 Introduction
Ocean Energy Europe has reported that wave energy
has the potential to capture over 3000 terawatt-hours
(TWh) annually, a quantity commensurate with the
predominant portion of European energy demand.
Additionally, the United States Department of
Energy has appraised the wave energy resource
potential to range between 1,594 and 2,640 TWh per
year, [1].
Accordingly, estimates for wave energy
production exhibit alignment, indicating an
aggregate potential of approximately 29,500 TWh
annually.
Concurrently, notable initiatives are underway
in the Spanish Basque Country through the Nereida
Marine Ocean Wave Current (MOWC)
experimental/commercial Project. Spearheaded by
the Basque Energy Agency (EVE), this undertaking
seeks to validate the feasibility of Oscillating Water
Column (OWC) technology with Wells turbine
power take-off. Situated within a recently
constructed breakwater in Mutriku on the north
coast of Spain, the project encompasses 16 18.5 kW
turbines, yielding a combined power output of 296
kW (Figure 1). This effort achieved a milestone of
2.4 gigawatt-hours (GWh) in 2021, [2].
Fig. 1: Turbo-generator modules employed within
the Mutriku wave plant
2 System Description
OWC-based converters serve as mechanisms
designed to convert the mechanical energy inherent
in oceanic waves into electrical power. This process
is facilitated through the utilization of a capture
chamber in conjunction with a turbo-generator
module.
WSEAS TRANSACTIONS on SYSTEMS and CONTROL
DOI: 10.37394/23203.2024.19.16
Aitor J. Garrido, Salvador Cayuela,
Amparo Villasante, Izaskun Garrido
E-ISSN: 2224-2856
153
Volume 19, 2024
The capture chamber, a stationary structure,
features a lower section exposed to the sea,
consistently submerged beneath the Still Water
Level (SWL), as depicted in the schematic
representation provided in Figure 2. In this
configuration, the undulating motion of ocean
waves induces an oscillatory airflow within the
chamber, alternately pushing and pulling internal
air, resulting in compression and decompression
cycles. Consequently, a pressure disparity arises
across the turbine, instigating its rotation. The
magnitude of this pressure differential (dp) can be
characterized by the following equation, [3], [4]:
 
󰇛
󰇛󰇜󰇜 (1)
where:
: Air density (kg/m3).
: Blade’s height (m).
: Blade’s chord’s length (m).
: Blades number.
: Blade’s section area (m2).
: Turbine’s mean diameter (m).
: Turbine’s angular speed (rad/s).
Fig. 2: System’s scheme
Regardless of the direction of airflow, the
rotational motion of the turbine remains unaffected,
primarily attributable to the specific design
characteristics inherent in the conventional turbine
configuration, notably the Wells turbine, [5], [6]
(Figure 3).
Fig. 3: Unused Wells turbine
The symmetric blade configuration
characteristic of self-rectifying turbines facilitates
unidirectional rotation. Nonetheless, this feature
also engenders an undesirable phenomenon known
as the Stalling effect, wherein the turbine ceases
rotation upon reaching a threshold airflow velocity.
Comprising the power take-off system (PTO),
the turbo-generator module incorporates both the
turbine and an induction generator, typically of the
Doubly Fed Induction Generator (DFIG) type [7],
[8], [9]. Its principal function revolves around
converting the oscillatory pressure differentials into
electrical power. Various control strategies may be
implemented to facilitate this process, [10], [11],
[12].
The mechanical dynamics governing the turbo-
generator block can be described by the following
equation:
󰇗 (2)
3 Theoretical Model Framework
Various system modeling approaches may be
employed to derive an appropriate model, [13], [14],
[15]. In this endeavor, the initial step entails
determining the pressure drop value based on the
characteristics of the input waves. As indicated by
equation (1), the properties of the input waves
directly influence the airflow speed. Consequently,
it is essential to establish a relationship between the
waves and airflow speed, [16], [17], [18].
To achieve this, wave dynamics must be taken
into account, which are governed by different
theories contingent upon the specific attributes of
the wave under consideration. According to Airy
linear theory, a wind wave can be represented as an
ideal sinusoidal wave, [19]. The mathematical
expression delineating the surface profile of such
waves is formulated as follows:
󰇛󰇜 󰇣
󰇛 󰇜󰇤 (3)
Furthermore, the calculation of the volume of
water in the chamber can be performed by
considering the volume of air present within the
Oscillating Water Column (OWC):
󰇛󰇜 󰇛󰇜 (4)
where  denote the volumes of the capture
chamber and water, respectively.
The water volume can also be determined by
integrating the variation of the water level across the
Oscillating Water Column (OWC) area, yielding the
expression:
WSEAS TRANSACTIONS on SYSTEMS and CONTROL
DOI: 10.37394/23203.2024.19.16
Aitor J. Garrido, Salvador Cayuela,
Amparo Villasante, Izaskun Garrido
E-ISSN: 2224-2856
154
Volume 19, 2024
󰇛󰇜 

 (5)
being l the length of the chamber.
Hence, the instantaneous airflow can be represented
as:
󰇛󰇜
 (6)
Now, by considering the geometry of the
chamber, it becomes feasible to derive the airflow
velocity required in equation (1):
󰇛󰇜

 
 (7)
being D the duct’s diameter.
However, as it may be observed from equation
(1), the rotational speed is also a required parameter.
This rotational speed is contingent upon the torque
exerted by the turbo-generator module, which can
be derived from Equation (2).
󰇟
󰇛󰇜󰇠 (8)
The Torque Coefficient (Ct) is, in turn,
interconnected with the Flow Coefficient (ф)
through the characteristic curves of the turbine, as
show in Figure 4, [20].
Fig. 4: Torque Coefficient vs. Flow Coefficient
Where the parameter ф represents a dimensionless
quantity associated with the tangent of the angle of
attack at the blade tip that can be straightforwardly
calculated as:
 (9)
Hence, the rotational speed of the turbo-
generator is derived from the established DFIG
equations.
By taken into account all the aforementioned
relationships outlined from Equation (1), an
expression for the pressure drop across the turbine,
in terms of the wave entry, can be formulated as
follows:
 󰇡
󰇢
󰇡


󰇢
󰇛 󰇜 (10)
Here, the function f(·) denotes the characteristic
curve specific to the turbine at hand.
4 Model Evaluation
In this section, the OWC plant model is computed
and subsequently validated through the use of
experimental data. To achieve this, the model output
is calculated using actual wave surface data series
gathered at the Mutriku plant breakwater using an
acoustic Doppler current profiler RDI 600, with the
aim of computing the theoretically predicted
pressure drop. The obtained results are subsequently
juxtaposed with in situ measured experimental dP
data furnished by the Mutriku OWC plant (Basque
Energy Board/EVE - BIMEP) for the corresponding
time frame. All procedures were conducted in
adherence to the requisite safety protocols and
occupational health guidelines within the production
facilities:
Fig. 5. dP comparison
As illustrated in Figure 5, the results present a
substantial level of agreement between the values
predicted by the model and the experimental data
acquired.
5 Conclusion
This paper outlines a thorough methodology for
modeling and simulating on-shore Oscillating Water
Column (OWC) systems. Building upon prior
relevant knowledge, the proposed approach marks a
notable progression in engineering proficiency, with
WSEAS TRANSACTIONS on SYSTEMS and CONTROL
DOI: 10.37394/23203.2024.19.16
Aitor J. Garrido, Salvador Cayuela,
Amparo Villasante, Izaskun Garrido
E-ISSN: 2224-2856
155
Volume 19, 2024
a particular emphasis on addressing control aspects.
The investigation focuses on the Mutriku MOWC
power plant and proceeds to validate its findings
utilizing experimental data.
The model incorporates the wave, chamber, and
turbo-generator modules, laying the groundwork for
validation via experimental data. The outcomes
reveal a significant congruence between model
projections and experimental observations, thereby
enabling deeper exploration into digital twin
implementation and advanced control strategies.
Acknowledgement:
The authors would like to thank the collaboration of
the Basque Energy Board/EVE - BIMEP through
Agreement UPV/EHUEVE23/6/2011. They would
also like to thank the colleagues of WAKE
European project and the BIMEP for their
collaboration and help.
References:
[1] IPPC. Renewable Energy Sources and
Climate Change Mitigation. IPPC Report.
Cambridge University Press. 2011, [Online].
https://www.ipcc.ch/report/renewable-energy-
sources-and-climate-change-mitigation/
(Accessed Day: May 2, 2024).
[2] BiMEP, [Online].
https://www.bimep.com/en/mutriku-
area/technical-characteristics/ (Accessed Day:
May 2, 2024).
[3] Alberdi, M., Amundarain, M., Barambones,
O., Romero, J.A. “Robust control of
oscillating water column (OWC) devices:
Power generation improvement”. OCEANS
2013 MTS/IEEE - San Diego, USA: An Ocean
in Common, art. no. 6740982. 2013.
[4] Amundarain, M, Alberdi, M. and De la Sen,
M., “Sliding-Mode Control of Wave Power
Generation Plants”, IEEE Transactions on
Industry Applications, vol. 48, no. 6, pp.
2372-2381, 2012. DOI:
10.1109/TIA.2012.2227096.
[5] L. Gato and V. Warfield, “Performance of a
high-solidity Wells turbine for an OWC wave
power plant”, in Proc.Euro Wave Energy
Symp., 1994, Edinburgh, UK, pp. 181-190.
[6] Lekube, J., Otaola, E., Maseda, J. “Flow
control in wells turbines for harnessing
maximum wave power” (2018), Sensors, 18
(2), art. no. 535, DOI: 10.3390/s18020535.
[7] Ahmad I., M'zoughi F., Aboutalebi P.,
Garrido I., Garrido A.J., Fuzzy logic control
of an artificial neural network-based floating
offshore wind turbine model integrated with
four oscillating water columns, Ocean
Engineering, 269, art. no. 113578, 2023, DOI:
10.1016/j.oceaneng.2022.113578.
[8] Romero, J.A., Carrascal, E., Sevillano-
Berasategui, G., Barambones, O. “Low effort
LI nuclear fusion plasma control using model
predictive control laws” in Mathematical
Problems in Engineering, art. no. 527420,
2015. DOI: 10.1155/2015/527420.
[9] Sevillano, M.G., Romero, J.A. “Robust
sliding mode control for tokamaks” in
Mathematical Problems in Engineering, art.
no. 341405, 2012. DOI:
10.1155/2012/341405.
[10] Garrido, A.J., De La Sen, M., Soto, J.C.,
Barambones, O., Garrido, I. “Suboptimal
regulation of a class of bilinear interconnected
systems with finite-time sliding planning
horizons” in Mathematical Problems in
Engineering, art. no. 817063, 2008. DOI:
10.1155/2008/817063.
[11] De la Sen M., Miñambres J.J., Almansa A.,
Soto J.C. “Basic theoretical results for expert
systems. Application to the supervision of
adaptation transients in planar robots”. (2004)
Artificial Intelligence, 152 (2), pp. 173 - 211,
DOI: 10.1016/S0004-3702(03)00136-X.
[12] Goretti Sevillano, M., “Sliding-mode loop
voltage control using ASTRA-Matlab
integration in Tokamak reactors” (2012)
International Journal of Innovative
Computing, Information and Control, 8 (9),
pp. 6473-6489.
[13] Lekube, J., “Rotational Speed Optimization in
Oscillating Water Column Wave Power Plants
Based on Maximum Power Point Tracking”
(2017), IEEE Transactions on Automation
Science and Engineering, 14 (2), art. no.
7552497. DOI: 10.1109/TASE.2016.2596579.
[14] Lekube, J., Otaola, E., Maseda, J. “Flow
control in wells turbines for harnessing
maximum wave power”, Sensors, 18 (2), art.
no. 535, 2018. DOI: 10.3390/s18020535.
[15] Barambones, O., Alkorta, P., Maseda, F.J.
“Simple linear models for plasma control in
Tokamak Reactors” 2008 International
Conference on Control, Automation and
Systems, ICCAS 2008. Fukuoka, Japan, art.
no. 4694261, pp. 2429-2432. 2008. DOI:
10.1109/ICCAS.2008.4694261.
[16] De la Sen, M., “On the robust adaptive
stabilization of a class of nominally first-order
hybrid systems”, IEEE Transactions on
WSEAS TRANSACTIONS on SYSTEMS and CONTROL
DOI: 10.37394/23203.2024.19.16
Aitor J. Garrido, Salvador Cayuela,
Amparo Villasante, Izaskun Garrido
E-ISSN: 2224-2856
156
Volume 19, 2024
Automatic Control, Vol. 44, No. 3, pp. 597-
602, 1999.
[17] Odysseas Chalatsakos, Ander Viñas, Ainhoa
Villaverde, Vicente Queral and Jesús A.
Romero, “Modelling and control of the
UPV/EHU Stellarator”, Proceedings of the
IEEE 23rd Mediterranean Conference on
Control and Automation (MED), June 16-19,
2015. Torremolinos, Spain (2015). DOI:
10.1109/MED.2015.7158785.
[18] Sevillano, M.G., “Control-oriented Automatic
System for Transport Analysis (ASTRA)-
Matlab integration for Tokamaks” (2011)
Energy, 36 (5), pp. 2812-2819. DOI:
10.1016/j.energy.2011.02.022.
[19] M. Alberdi, M. Amundarain and O.
Barambones “Performance of an Ocean
Energy Conversion System with DFIG
Sensorless Control” in Mathematical
Problems in Engineering, vol. 2013, no. 1, pp.
14, 2013. DOI: 10.1155/2013/260514.
[20] M'zoughi F., De La Sen M. “Self-adaptive
global-best harmony search algorithm-based
airflow control of a wells-turbine-based
oscillating-water column”, Applied Sciences,
10 (13), art. no. 4628, 2020. DOI:
10.3390/app10134628.
Contribution of Individual Authors to the
Creation of a Scientific Article (Ghostwriting
Policy)
The authors equally contributed in the present
research, at all stages from the formulation of the
problem to the final findings and solution.
Sources of Funding for Research Presented in a
Scientific Article or Scientific Article Itself
The authors would like to thank the Basque
Government for partially funding their research
work through Grant IT1555-22 and they thank
MICIU/AEI/ 10.13039/501100011033 and ERDF/E
for partially funding their research work through
Grants PID2021-123543OB-C21 and PID2021-
123543OB-C22.
Conflict of Interest
The authors have no conflicts of interest to declare.
Creative Commons Attribution License 4.0
(Attribution 4.0 International, CC BY 4.0)
This article is published under the terms of the
Creative Commons Attribution License 4.0
https://creativecommons.org/licenses/by/4.0/deed.en
_US
WSEAS TRANSACTIONS on SYSTEMS and CONTROL
DOI: 10.37394/23203.2024.19.16
Aitor J. Garrido, Salvador Cayuela,
Amparo Villasante, Izaskun Garrido
E-ISSN: 2224-2856
157
Volume 19, 2024