Document Type


Publication Date



Department of Mechanical Engineering-Engineering Mechanics


The motivation and focus of this work is to generate passive transfer function matrices that model the radiation forces for an array of WECs. Multivariable control design is often based on linear time-invariant (LTI) systems such as state-space or transfer function matrix models. The intended use is for designing real-time control strategies where knowledge of the model’s poles and zeros is helpful. This work presents a passivity-based approach to estimate radiation force transfer functions that accurately replace the convolution operation in the Cummins equation while preserving the physical properties of the radiation function. A two-stage numerical optimization approach is used, the first stage uses readily available algorithms for fitting a radiation damping transfer function matrix to the system’s radiation frequency response. The second stage enforces additional constraints on the form of the transfer function matrix to increase its passivity index. After introducing the passivity-based algorithm to estimate radiation force transfer functions for a single WEC, the algorithm was extended to a WEC array. The proposed approach ensures a high degree of match with the radiation function without degrading its passivity characteristics. The figures of merit that will be assessed are (i) the accuracy of the LTI systems in approximating the radiation function, as measured by the normalized root mean squared error (NRMSE), and (ii) the stability of the overall system, quantified by the input passivity index, (Formula presented.), of the radiation force transfer function matrix.

Publisher's Statement

Copyright: © 2023 by the authors. Licensee MDPI, Basel, Switzerland. Publisher’s version of record:

Publication Title


Creative Commons License

Creative Commons Attribution 4.0 International License
This work is licensed under a Creative Commons Attribution 4.0 International License.


Publisher's PDF