Bor Plestenjak,Daniel A. Kiefer,Hauke Gravenkamp
Abstract:Eigenvalues of parameter-dependent quadratic eigenvalue problems form eigencurves. The critical points on these curves, where the derivative vanishes, are of practical interest.
A particular example is found in the dispersion curves of elastic waveguides, where such points are called zero-group-velocity (ZGV) points. Recently, it was revealed that the problem of computing ZGV points can be modeled as a multiparameter eigenvalue problem (MEP), and several numerical methods were devised. Due to their complexity, these methods are feasible only for problems involving small matrices. In this paper, we improve the efficiency of these methods by exploiting the link to the Sylvester equation. This approach enables the computation of ZGV points for problems with much larger matrices, such as multi-layered plates and three-dimensional structures of complex cross-sections.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
The paper aims to solve the efficiency problem of calculating Zero - Group - Velocity (ZGV) points in waveguides. Specifically, the paper focuses on the Parameter - Dependent Quadratic Eigenvalue Problem (QEP), which has the form:
\[
W(k, \omega)u := ((ik)^2L_2 + ikL_1 + L_0 + \omega^2M)u = 0
\]
where \(L_0, L_1, L_2\) and \(M\) are real \(n\times n\) matrices, usually obtained by (semi - ) discretization of boundary value problems. The solutions \((k, \omega)\) form characteristic curves \(\omega(k)\), and we are interested in the critical points where the derivative on these curves is zero, that is, \(\omega'(k)=\frac{\partial\omega}{\partial k} = 0\).
In the application of elastic waveguides, these characteristic curves are called dispersion curves, and the slope \(c_g=\omega'\) is called the group velocity, which describes the speed of energy propagation. The points on the dispersion curve where the group velocity is zero are called Zero - Group - Velocity (ZGV) points. For complex multilayer plates and three - dimensional structures, calculating these points is challenging, especially when the matrix scale is large.
### Main contributions
To improve the computational efficiency, the author introduced a method based on the Sylvester equation. This method can significantly accelerate the calculation and reduce the memory requirements, thus making it possible to handle larger - scale matrix problems. Specific improvements include:
1. **Utilizing the Sylvester equation**: By transforming the original problem into the Sylvester equation, the linear system can be solved more efficiently.
2. **Combining the Method of Fixed Relative Distance (MFRD)**: Provide initial approximate values, and then use the locally convergent Gauss - Newton method for refinement.
3. **Scanning algorithm**: Combining MFRD and the Gauss - Newton method, a scanning algorithm was developed to calculate all ZGV points within a given wavenumber interval.
### Application background
The main application background of this research is the Lamb wave calculation in anisotropic elastic waveguides. Through the above methods, ZGV points can be effectively found in larger - scale and more complex structures, which is of great significance for understanding the energy propagation characteristics in waveguides.
In summary, by introducing new mathematical tools and algorithms, this paper solves the efficiency problem of calculating zero - group - velocity points in complex waveguide structures, providing strong support for practical engineering applications.