Matrix structure and convergence behavior of the matched eigenfunction method for computing heave wave forces on generalized concentric bodies

This paper introduces a unified matched eigenfunction expansion method (MEEM) framework for generalized concentric bodies that demonstrates significantly faster convergence and smaller matrix sizes compared to traditional boundary element methods, while maintaining high accuracy for both vertical and slanted geometries.

The Gist

Imagine you are trying to predict how a giant, floating offshore structure (like a wave energy converter) will bob up and down when hit by ocean waves. To do this safely and efficiently, engineers need to calculate the "push" and "pull" forces the water exerts on the structure.

For decades, the standard way to do this has been like trying to map a coastline by taking millions of tiny, individual measurements with a ruler. This method, called the Boundary Element Method (BEM), is accurate but incredibly slow and computationally heavy. It's like trying to solve a puzzle by cutting every single piece into a million smaller fragments just to be sure they fit.

This paper introduces a smarter, faster way to solve the same puzzle using a method called Matched Eigenfunction Expansion (MEEM). Here is how the paper explains it, using simple analogies:

1. The "Lego Tower" vs. The "Pixelated Image"

The standard method (BEM) treats the water around the object like a digital image made of millions of tiny pixels. To get a clear picture, you need a massive number of pixels, which takes a long time to process.

The new method (MEEM) treats the water like a Lego tower built from specific, pre-made shapes. Instead of measuring every tiny point, the math breaks the water down into concentric rings (like tree rings or a target). Inside each ring, the water's movement is described by a known mathematical "recipe" (an eigenfunction). You only need to figure out the "ingredients" (coefficients) for a few of these recipes to get the whole picture.

2. The "Matching Game"

The core trick of this method is matching. Imagine you have a series of nested rings of water. The method ensures that the water pressure and speed flow smoothly from one ring to the next, just like ensuring the water level is the same where two connected buckets meet.

The authors organized these matching rules into a giant matrix (a grid of numbers). They discovered this grid has a very specific, sparse pattern—like a highway with only two lanes of traffic instead of a gridlock of cars. Because the grid is so organized and "sparse," the computer can solve it incredibly fast.

3. Handling "Slanted" Shapes

Real-world objects aren't always perfect cylinders; they often have slanted sides (like a cone or a funnel). The standard way to handle this with MEEM is to approximate the slant by stacking many thin, flat rings on top of each other, like a staircase trying to mimic a ramp.

The paper tested how many "steps" are needed to make the staircase look like a smooth ramp. They found that:

• Gentle slopes need fewer steps.
• Steep slopes need more steps.
• Even with a "staircase" approximation, the method can predict the forces on the object with less than 5% error, even for steep angles, which is accurate enough for engineering.

4. The Speed Demon

The most exciting finding is the speed comparison. The authors pitted their new method against the industry-standard software (Capytaine).

• Accuracy: Both methods can achieve the same level of accuracy (2% error).
• Speed: The new method is 10 times faster (an order of magnitude).
• Size: The new method uses a mathematical "matrix" that is 100 times smaller (two orders of magnitude) than the one used by the standard method.

The Analogy: If the standard method is like driving a heavy truck through a city to deliver a package, the new method is like using a high-speed drone. They both get the package to the same destination, but the drone gets there much faster and with less fuel.

5. Why This Matters

The paper concludes that this method is a powerful tool for optimization. Because it is so fast, engineers can now test thousands of different shapes for offshore structures in the time it used to take to test just one. This allows them to find the "perfect" design much quicker, potentially saving money and improving the safety of marine structures.

In summary: The paper proves that by using a clever mathematical "recipe" approach instead of a brute-force "pixel" approach, we can calculate wave forces on floating structures much faster and with smaller computer requirements, without losing accuracy.

Technical Summary

Technical Summary: Matrix Structure and Convergence Behavior of the Matched Eigenfunction Method for Computing Heave Wave Forces on Generalized Concentric Bodies

Problem Statement
The calculation of hydrodynamic coefficients (added mass, radiation damping, and excitation forces) for offshore structures is a critical component of marine design and optimization. Traditional Boundary Element Method (BEM) solvers, while versatile for arbitrary geometries, often present a computational bottleneck due to their reliance on dense mesh discretization, which scales poorly with the number of panels required for high accuracy. While the Matched Eigenfunction Expansion Method (MEEM) offers a semi-analytical alternative with potential computational benefits, its adoption in offshore engineering has been limited. Previous literature has lacked direct comparisons of MEEM's accuracy and performance against modern BEM solvers, and the method's convergence behavior—particularly for slanted (conical) geometries and general concentric bodies—remains undocumented. Consequently, the practical utility of MEEM for complex, non-vertical geometries is unclear.

Methodology
This work presents a unifying MEEM framework capable of modeling an arbitrary number of fixed or heaving, surface-piercing, concentric annular cylinders with continuous, radially-monotonic profiles. The methodology involves:

1. Mathematical Formulation: The fluid domain is divided into internal regions beneath each cylindrical ring and a single external region. Linear potential flow theory is applied, separating the velocity potential into incident, diffracted, and radiated components. The radiated potential in each region is expressed as a series of eigenfunctions (involving Bessel functions and trigonometric/hyperbolic functions) with unknown coefficients.
2. Matrix Structure: By enforcing continuity of potential and radial velocity at the interfaces between regions, and zero radial velocity at body boundaries, a system of linear algebraic equations ($A\mathbf{x} = \mathbf{b}$) is derived. The authors demonstrate that the system matrix $A$ possesses a block bi-diagonal structure. This structure arises because matching conditions at a specific boundary only couple the unknowns of the two adjacent regions.
3. Numerical Robustness: The implementation addresses several numerical challenges, including overflow in Bessel functions for high frequencies or deep water (mitigated via exponentially scaled Bessel functions and analytical limits) and high condition numbers in the linear system.
4. Slanted Geometry Approximation: To model bodies with slanted or conical surfaces (e.g., CorPower, WaveBot), the method approximates the continuous slope by discretizing it into a finite number of stepped annular rings. The hydrodynamic forces are then integrated over this discretized wetted surface.
5. Validation and Benchmarking: The framework, implemented in the open-source Python package OpenFLASH, is validated against the BEM solver Capytaine. A comprehensive convergence study is conducted to determine the number of eigenfunction terms required for a target accuracy based on dimensionless geometric parameters.

Key Contributions

• Unified Framework: A generalized MEEM formulation for $M$ concentric regions with arbitrary fixed or heaving degrees of freedom, organized into a sparse block matrix structure.
• Convergence Analysis: A systematic study identifying dimensionless parameters (such as the ratio of fluid height to radial width) that govern the convergence rate of hydrodynamic coefficients. The authors develop predictive models to estimate the number of terms needed to achieve a specific error tolerance (e.g., 1% or 2%).
• Slanted Geometry Handling: An evaluation of the step-discretization approach for slanted bodies, quantifying the relationship between slant steepness, the number of subdivisions, and the resulting error in hydrodynamic coefficients.
• Performance Benchmarking: A direct comparison of MEEM and Capytaine regarding accuracy, matrix size, and computational runtime.

Results

• Accuracy: MEEM computes hydrodynamic coefficients for slanted geometries within 5% of Capytaine results, even for slant angles as steep as 15° from the vertical, provided sufficient subdivisions are used.
• Convergence: The study establishes that damping coefficients generally converge faster than added mass. The developed predictive models can estimate the required number of terms to achieve less than 2% error in hydrodynamic coefficients with 95% confidence for random geometric configurations.
• Computational Efficiency: MEEM demonstrates significant speed advantages over Capytaine. To achieve 2% convergence in hydrodynamic coefficients, MEEM requires a matrix size two orders of magnitude smaller (e.g., side length ~60 vs. >1000 for Capytaine) and runs approximately one order of magnitude faster (0.01s vs. 0.1s in the tested scenario).
• Numerical Stability: Despite high condition numbers in the system matrix for certain geometries, the method maintains accuracy, and the use of exponentially scaled Bessel functions extends the valid range of input parameters.

Significance
The paper claims that MEEM serves as a computationally effective alternative to traditional BEM solvers for axisymmetric marine structures. By leveraging its semi-analytical nature and sparse block matrix structure, MEEM enables the rapid calculation of hydrodynamic coefficients with high accuracy. This efficiency is particularly significant for design optimization studies, where hydrodynamic models are often the computational bottleneck. The ability to accurately model slanted geometries and predict convergence requirements a priori allows engineers to perform optimization with increased speed and confidence, potentially yielding improved designs for offshore renewable energy devices and other marine structures. The authors position the OpenFLASH package as a tool to facilitate these future optimization studies and further generalizations of the method.