Matrix structure and convergence behavior of the matched eigenfunction method for computing heave wave forces on generalized concentric bodies
Pith reviewed 2026-05-20 01:45 UTC · model grok-4.3
The pith
The matched eigenfunction expansion method computes heave wave forces on slanted concentric cylinders within 5 percent of boundary element results using matrices two orders of magnitude smaller.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the matched eigenfunction expansion method supplies a single framework for an arbitrary number of fixed or heaving surface-piercing annular cylinders with continuous and radially monotonic profiles; numerical experiments then show that this framework approximates hydrodynamic coefficients of slanted geometries to within 5 percent of Capytaine even at 15 degrees from vertical, while reaching 2 percent convergence an order of magnitude faster than Capytaine and with a matrix size two orders of magnitude smaller.
What carries the argument
The block matrix structure obtained by dividing the fluid domain into annular regions and matching eigenfunction expansions of the velocity potential at each interface, which organizes the linear system solved for the added mass, radiation damping, and excitation coefficients.
If this is right
- Hydrodynamic coefficients for a broad range of continuous radially monotonic shapes become computable with far less effort than boundary element methods.
- Optimization loops that iterate over many offshore structure geometries become practical because each evaluation finishes an order of magnitude sooner.
- Both fixed and heaving surface-piercing bodies can be treated inside the same matching procedure.
- Convergence rates are now documented for slanted profiles, giving users a clear rule of thumb for choosing matrix size.
Where Pith is reading between the lines
- The compact matrix size could let the method run inside larger coupled simulations that include multiple structures or time-varying waves.
- A hybrid scheme that switches to boundary element methods only near sharp features might extend the approach to bodies that violate the monotonic-profile assumption.
- Insights into the block structure may suggest similar semi-analytical reductions for related problems such as diffraction by arrays of cylinders or motion in irregular waves.
Load-bearing premise
The body profiles must stay continuous and radially monotonic so that the eigenfunction expansions in neighboring annular regions can be matched directly without extra interface conditions or geometric approximations.
What would settle it
A direct numerical comparison on a body whose profile contains a radial discontinuity or reversal would show whether the hydrodynamic coefficients diverge sharply from a high-resolution boundary element reference solution.
Figures
read the original abstract
Structural survival of offshore structures is crucial for the growing marine economy. Calculating the added mass, radiation damping, and excitation coefficients to quantify wave loads with the traditional boundary element method (BEM) presents a computational bottleneck. The matched eigenfunction expansion method (MEEM), a long-known but rarely-used alternative, offers computational benefits due to its semi-analytical nature. However, previous work fails to directly compare its accuracy and computational performance with BEM, leaving the extent of its utility unknown. Furthermore, the geometry-dependent convergence for cylindrical and slanted geometries has not yet been documented, making the method's practicality for general geometries unclear. This paper presents a unifying MEEM framework for modeling an arbitrary number of fixed or heaving surface-piercing annular cylinders with continuous and radially-monotonic body profiles, and explores the method's block matrix structure, convergence behavior, ability to accurately approximate slanted geometries, and computational advantages over the BEM solver Capytaine. The numerical experiments show that MEEM can compute hydrodynamic coefficients of slanted geometries within 5% of Capytaine, even for angles as steep as 15 degrees from vertical. Finally, MEEM can achieve 2% convergence of its hydrodynamic coefficients an order of magnitude faster than Capytaine with a matrix size two orders of magnitude smaller, making it a computationally effective alternative to traditional BEM solvers. These contributions enable hydrodynamic analysis of a broad range of shapes with increased speed and confidence, paving the way for future optimization studies to yield improved designs.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a unifying matched eigenfunction expansion method (MEEM) framework for an arbitrary number of fixed or heaving surface-piercing annular cylinders possessing continuous and radially-monotonic body profiles. It analyzes the resulting block-matrix structure, documents geometry-dependent convergence for both cylindrical and slanted cases, demonstrates approximation of slanted geometries, and reports computational comparisons against the BEM solver Capytaine, claiming 5% agreement on hydrodynamic coefficients for slants up to 15° from vertical together with order-of-magnitude speed-ups at 2% convergence using matrices two orders of magnitude smaller.
Significance. If the accuracy and performance claims survive clarification of the geometric-proxy issue, the work supplies a semi-analytical, matrix-structured alternative to BEM that could materially accelerate hydrodynamic-coefficient evaluation for generalized offshore geometries and thereby support subsequent optimization studies.
major comments (2)
- [Numerical experiments / slanted-geometry results] Numerical experiments (results section on slanted geometries): the headline claim that MEEM reproduces Capytaine coefficients within 5% for 15° slants rests on replacing the true linear slant with a continuous radially-monotonic profile for the MEEM computation. It is not stated whether the Capytaine reference solutions are obtained on the identical monotonic proxy or on the exact slanted geometry; without this clarification the reported tolerance conflates profile-approximation error with MEEM truncation and interface-matching error, weakening attribution of the observed accuracy to the method itself.
- [Convergence behavior and performance comparison] Convergence and performance subsection: the statements that 2% convergence is reached an order of magnitude faster with a matrix two orders of magnitude smaller lack explicit truncation criteria (number of retained eigenfunctions per region), definition of matrix size (total unknowns versus per-annulus blocks), Capytaine mesh resolution and convergence tolerance, and precise error-bar construction. These omissions make the quantitative speed-up claims difficult to reproduce or compare directly.
minor comments (2)
- [Abstract] Abstract: the phrase 'generalized concentric bodies' is immediately qualified by the continuous radially-monotonic restriction; a single clarifying sentence would prevent readers from over-generalizing the scope.
- [Throughout] Notation: ensure that the block-matrix partitioning and the matching conditions at vertical interfaces are denoted consistently between the theoretical derivation and the numerical implementation sections.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments, which help strengthen the clarity and reproducibility of the work. We address each major comment below and will incorporate revisions to resolve the identified ambiguities.
read point-by-point responses
-
Referee: Numerical experiments (results section on slanted geometries): the headline claim that MEEM reproduces Capytaine coefficients within 5% for 15° slants rests on replacing the true linear slant with a continuous radially-monotonic profile for the MEEM computation. It is not stated whether the Capytaine reference solutions are obtained on the identical monotonic proxy or on the exact slanted geometry; without this clarification the reported tolerance conflates profile-approximation error with MEEM truncation and interface-matching error, weakening attribution of the observed accuracy to the method itself.
Authors: We agree that explicit clarification is required. In the numerical experiments, Capytaine reference solutions were computed on the exact linear-slant geometry, while MEEM was applied to a continuous radially-monotonic profile chosen to approximate the slant. The reported 5% agreement therefore represents the combined effect of profile approximation and MEEM truncation/interface-matching error. Because the manuscript positions MEEM as a method for monotonic profiles that can approximate slanted bodies, this comparison is intentional; however, we will revise the text to state the distinction unambiguously so that readers can correctly attribute the observed accuracy. revision: yes
-
Referee: Convergence and performance subsection: the statements that 2% convergence is reached an order of magnitude faster with a matrix two orders of magnitude smaller lack explicit truncation criteria (number of retained eigenfunctions per region), definition of matrix size (total unknowns versus per-annulus blocks), Capytaine mesh resolution and convergence tolerance, and precise error-bar construction. These omissions make the quantitative speed-up claims difficult to reproduce or compare directly.
Authors: We accept that the current presentation lacks sufficient detail for direct reproduction. In the revised manuscript we will add: (i) the truncation criterion (number of retained eigenfunctions retained in each fluid region), (ii) an explicit definition of matrix size as the total number of unknowns in the assembled block system, (iii) the Capytaine mesh density and solver tolerance employed for the reference solutions, and (iv) the precise procedure used to construct the error bars shown in the convergence plots. These additions will allow readers to verify the reported order-of-magnitude speed-up at 2% convergence. revision: yes
Circularity Check
No significant circularity; quantitative claims rest on external Capytaine benchmarks
full rationale
The paper defines a MEEM framework for continuous radially-monotonic annular bodies, derives the block-matrix structure and interface matching conditions from standard eigenfunction expansions in annular regions, and validates all accuracy (within 5% for 15° slants) and convergence (2% with order-of-magnitude speed-up) claims through direct numerical comparison to the independent Capytaine BEM solver. No fitted parameters are renamed as predictions, no self-citation chain is load-bearing for the central results, and the monotonic-profile assumption is stated explicitly as a modeling precondition rather than derived from the outputs. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The velocity potential in each fluid region can be expanded in a complete set of eigenfunctions that satisfy the linearized free-surface and bottom boundary conditions.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Convergence … modeled as a power-law envelope ϵ ≈ (N_im / β_im )^−α_im (Fig. 6).
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Michelén and Ancellin, Matthieu and Haji, Maha , month = jan, year =
Khanal, Kapil and Ströfer, Carlos A. Michelén and Ancellin, Matthieu and Haji, Maha , month = jan, year =. Fully. doi:10.48550/arXiv.2501.06988 , abstract =
-
[2]
McCabe, Rebecca and Khanal, Kapil and Haji, Maha , month = aug, year =. Open-source toolbox for semi-analytical hydrodynamic coefficients via the matched eigenfunction expansion method , url =. doi:10.5281/zenodo.14504017 , abstract =
-
[3]
Best, Hope and Khanal, Kapil and McCabe, Rebecca and Jiang, Ruiyang and Treacy, Collin and Haji, Maha , year =
-
[4]
Chau, Fun Pang and Yeung, Ronald W. , month = jan, year =. Inertia and. 25th
- [5]
-
[6]
Nguyen, Nhu , month = jun, year =. Theoretical. doi:10.1115/OMAE2024-128018 , abstract =
-
[7]
Applied Ocean Research , author =
Theoretical modeling of a bottom-raised oscillating surge wave energy converter structural loadings and power performances , volume =. Applied Ocean Research , author =. 2024 , keywords =. doi:10.1016/j.apor.2024.104031 , abstract =
-
[8]
Hiptmair, Ralf and Moiola, Andrea and Perugia, Ilaria , year =. A survey of. Building bridges: connections and challenges in modern approaches to numerical
-
[9]
A novel two-objective optimization computational framework for a two-body heaving wave energy converter , volume =. Renewable Energy , author =. 2022 , keywords =. doi:10.1016/j.renene.2022.04.006 , abstract =
-
[10]
Advances in Engineering Software , author =
Trefftz method: an overview , volume =. Advances in Engineering Software , author =. 1995 , keywords =. doi:10.1016/0965-9978(95)00067-4 , abstract =
-
[11]
Seah, Robert K. M. and Yeung, Ronald W. , month = apr, year =. Symmetric
- [12]
-
[13]
Fuchs, R. A. and MacCamy, R. C. , month = dec, year =. Wave forces on piles: a diffraction theory , shorttitle =
-
[14]
Proceedings of the Royal Society of London
The forces on a circular cylinder submerged in a uniform stream , volume =. Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences , author =. 1936 , pages =. doi:10.1098/rspa.1936.0212 , abstract =
-
[15]
Michele, Simone and Sammarco, Paolo and d’Errico, Michele , month = aug, year =. Theory of the synchronous motion of an array of floating flap gates oscillating wave surge converter , volume =. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences , publisher =. doi:10.1098/rspa.2016.0174 , abstract =
-
[16]
Edwards, Emma C. and Yue, Dick K.-P. , month = feb, year =. Optimisation of the geometry of axisymmetric point-absorber wave energy converters , volume =. Journal of Fluid Mechanics , publisher =. doi:10.1017/jfm.2021.993 , abstract =
-
[17]
Garrett, C. J. R. , month = sep, year =. Bottomless harbours , volume =. Journal of Fluid Mechanics , publisher =. doi:10.1017/S0022112070002495 , abstract =
-
[18]
Chatjigeorgiou, Ioannis K. , year =. Analytical. doi:10.1017/9781316838983 , abstract =
-
[19]
Applied Ocean Research , publisher =
Reflection of oblique ocean water waves by a vertical porous structure placed on a multi-step impermeable bottom , volume =. Applied Ocean Research , publisher =. 2014 , pages =. doi:10.1016/j.apor.2014.07.001 , abstract =
-
[20]
Garrett, C. J. R. , month = mar, year =. Wave forces on a circular dock , volume =. Journal of Fluid Mechanics , publisher =. doi:10.1017/S0022112071000430 , abstract =
-
[21]
Second-order hydrodynamic effects on an arrangement of two concentric truncated vertical cylinders , volume =. Marine Structures , author =. 2009 , keywords =. doi:10.1016/j.marstruc.2008.12.003 , abstract =
-
[22]
Hydrodynamic characteristics of floating toroidal bodies , volume =. Ocean Engineering , author =. 1997 , pages =. doi:10.1016/S0029-8018(96)00020-0 , abstract =
-
[23]
China Ocean Engineering , author =
Hydrodynamic analysis and shape optimization for vertical axisymmetric wave energy converters , volume =. China Ocean Engineering , author =. 2016 , keywords =. doi:10.1007/s13344-016-0062-2 , abstract =
-
[24]
Behaviour of vertical bodies of revolution in waves , volume =. Ocean Engineering , author =. 1986 , pages =. doi:10.1016/0029-8018(86)90037-5 , abstract =
-
[25]
Coastal Engineering , author =
Analysis of orthogonal wave reflection by a caisson with open front chamber filled with sloping rubble mound , volume =. Coastal Engineering , author =. 2014 , keywords =. doi:10.1016/j.coastaleng.2014.05.002 , abstract =
-
[26]
Applied Ocean Research , author =
Added mass and damping of a vertical cylinder in finite-depth waters , volume =. Applied Ocean Research , author =. 1981 , pages =. doi:10.1016/0141-1187(81)90101-2 , abstract =
-
[27]
A novel solution to the second-order wave radiation force on an oscillating truncated cylinder based on the application of control surfaces , volume =. Ocean Engineering , author =. 2020 , keywords =. doi:10.1016/j.oceaneng.2020.107278 , abstract =
-
[28]
Applied Ocean Research , author =
The added mass of bodies heaving at low frequency in water of finite depth , volume =. Applied Ocean Research , author =. 1991 , pages =. doi:10.1016/S0141-1187(05)80036-7 , abstract =
-
[29]
IEEE Journal of Oceanic Engineering , author =
Hydrodynamic. IEEE Journal of Oceanic Engineering , author =. 2015 , note =. doi:10.1109/JOE.2014.2344951 , abstract =
-
[30]
Performance assessment of a. Renewable Energy , author =. 2024 , keywords =. doi:10.1016/j.renene.2024.121547 , abstract =
-
[31]
European Journal of Mechanics - B/Fluids , author =
Hydrodynamics of the oscillating wave surge converter in the open ocean , volume =. European Journal of Mechanics - B/Fluids , author =. 2013 , keywords =. doi:10.1016/j.euromechflu.2013.01.007 , abstract =
-
[32]
Applied Ocean Research , author =
Hydrodynamic coefficients in heave of two concentric surface-piercing truncated circular cylinders , volume =. Applied Ocean Research , author =. 2004 , keywords =. doi:10.1016/j.apor.2005.03.002 , abstract =
-
[33]
Chau, Fun Pang and Yeung, Ronald W. , year =. Inertia,. doi:10.1115/OMAE2012-83987 , abstract =
-
[34]
Chau, Fun Pang and Yeung, Ronald W. , month = jan, year =. Inertia and
-
[35]
Theory and applications of ocean surface waves: Part 1: linear aspects , author=. 2005 , publisher=
work page 2005
-
[36]
IET Renewable Power Generation , volume=
Energy-maximising moment-based constrained optimal control of ocean wave energy farms , author=. IET Renewable Power Generation , volume=. 2021 , publisher=
work page 2021
-
[37]
Analytical methods in marine hydrodynamics , author=. 2018 , publisher=
work page 2018
-
[38]
International Conference on Offshore Mechanics and Arctic Engineering , volume=
Inertia, damping, and wave excitation of heaving coaxial cylinders , author=. International Conference on Offshore Mechanics and Arctic Engineering , volume=. 2012 , organization=
work page 2012
- [39]
- [40]
-
[41]
T. Barker and D. G. Schaeffer and P. Boh\'orquez and J. M. N. T. Gray , title =. J. Fluid. Mech. , year =
- [42]
-
[43]
A. V. Briukhanov and S. S. Grigorian and S. M. Miagkov and M. Y. Plam and I. E. Shurova and M. E. Eglit and Y. L. Yakimov , title =. Physics of Snow and Ice, Proceedings of the International Conference on Low Temperature Science , volume =
- [44]
- [45]
-
[46]
S.C.R. Dennis , title =. Ninth Intl Conf. on Numerical Methods in Fluid Dynamics , publisher =. 1985 , editor =
work page 1985
-
[47]
Formation of levees, troughs and elevated channels by avalanches on erodible slopes , author =. J. Fluid Mech. , volume =
- [48]
-
[49]
and Saut, Jean Claude , title =
Joseph, Daniel D. and Saut, Jean Claude , title =. Theoretical and Computational Fluid Dynamics , publisher=
-
[50]
M.G. Worster , title =. Interactive dynamics of convection and solidification , publisher =. 1992 , editor =
work page 1992
- [51]
- [52]
-
[53]
C.M. Linton and D.V. Evans , title =. Phil.\ Trans.\ R. Soc.\ Lond. , year =
- [54]
- [55]
- [56]
-
[57]
On the oscillations Near and at resonance in open pipes , year =
L. On the oscillations Near and at resonance in open pipes , year =. J. Engng Maths , volume =
- [58]
-
[59]
Camera Calibration Toolbox for Matlab , howpublished =
Bouguet, J.-Y , year=. Camera Calibration Toolbox for Matlab , howpublished =
-
[60]
Techno-economic related metrics for a wave energy converters feasibility assessment , author=. Sustainability , volume=. 2016 , publisher=
work page 2016
-
[61]
IEEE Transactions on Sustainable Energy , volume=
Control co-design of power take-off systems for wave energy converters using WecOptTool , author=. IEEE Transactions on Sustainable Energy , volume=. 2023 , publisher=
work page 2023
-
[62]
Super capacitor energy storage system design for wave energy converter demonstration , author=. 2020 International Symposium on Power Electronics, Electrical Drives, Automation and Motion (SPEEDAM) , pages=. 2020 , organization=
work page 2020
-
[63]
International Marine Energy Journal , volume=
Degrees of freedom effects on a laboratory scale WEC point absorber , author=. International Marine Energy Journal , volume=. 2025 , url=
work page 2025
-
[64]
Methodology for design and economic analysis of marine energy conversion (MEC) technologies , author=. 2014 , journal=
work page 2014
-
[65]
Hydrodynamic characteristics of floating toroidal bodies , author=. Ocean engineering , volume=. 1997 , publisher=
work page 1997
-
[66]
International Conference on Offshore Mechanics and Arctic Engineering , volume=
Hydrodynamic analysis of a vertical axisymmetric oscillating water column device floating in finite depth waters , author=. International Conference on Offshore Mechanics and Arctic Engineering , volume=. 2012 , organization=
work page 2012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.