A p-Multigrid Accelerated Nodal Spectral Element Method for Free-Surface Incompressible Navier-Stokes Model of Nonlinear Water Waves

Anders Melander,Wojciech Laskowski,Spencer J. Sherwin,Allan P. Engsig-Karup

2024-11-22

Abstract:We present a spectral element model for general-purpose simulation of non-overturning nonlinear water waves using the incompressible Navier-Stokes equations (INSE) with a free surface. The numerical implementation of the spectral element method is inspired by the related work by Engsig-Karup et al. (2016) and is based on nodal Lagrange basis functions, mass matrix-based integration and gradient recovery using global $L^2$ projections. The resulting model leverages the high-order accurate -- possibly exponential -- error convergence and has support for geometric flexibility allowing for computationally efficient simulations of nonlinear wave propagation. An explicit fourth-order accurate Runge-Kutta scheme is employed for the temporal integration, and a mixed-stage numerical discretization is the basis for a pressure-velocity coupling that makes it possible to maintain high-order accuracy in both the temporal and spatial discretizations while preserving mass conservation. Furthermore, the numerical scheme is accelerated by solving the discrete Poisson problem using an iterative solver strategy based on a geometric $p$-multigrid method. This problem constitutes the main computational bottleneck in INSE models. It is shown through numerical experiments, that the model achieves spectral convergence in the velocity fields for highly nonlinear waves, and there is excellent agreement with experimental data for the simulation of the classical benchmark of harmonic wave generation over a submerged bar. The geometric $p$-multigrid solver demonstrates $O(n)$ computational scalability simulations, making it a suitable efficient solver strategy as a candidate for extensions to more complex, real-world scenarios.

Numerical Analysis,Fluid Dynamics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop an efficient and accurate numerical method for simulating nonlinear water wave phenomena. Specifically, the authors proposed a free - surface incompressible Navier - Stokes equation model based on the high - order nodal spectral element method, aiming to achieve the following goals:

1. **Efficient and accurate simulation**: By using the high - order spectral element method, this model can achieve rapid convergence (spectral convergence), thus obtaining high - precision results at a lower computational cost.
2. **Handling complex geometries**: This method supports geometric flexibility and allows for efficient simulation of domains with complex boundaries.
3. **Accelerating the solution of the Poisson problem**: To improve computational efficiency, the authors introduced the geometric p - multigrid method to accelerate the solution of the discretized Poisson problem, which is the main computational bottleneck in the INSE model.

### Specific problem description

The main challenges mentioned in the paper include:

- **Nonlinear water wave simulation**: Accurately simulating nonlinear water waves is very important in fields such as oceanography, coastal and offshore engineering, and environmental science. Although the existing fully nonlinear potential flow (FNPF) models have relatively high computational efficiency, they cannot consider viscous and rotational effects, resulting in over - prediction of wave heights in wave - structure interaction applications.
- **Selection of numerical methods**: Traditional implicit surface - capturing techniques (such as the level - set method and the VOF method) can handle breaking waves, but they have disadvantages in terms of computational cost and overall accuracy. Therefore, it is necessary to develop an explicit tracking method that can maintain high accuracy and handle free surfaces.
- **Solution of the Poisson problem**: In CFD solvers, the solution of the Poisson problem is usually a significant computational bottleneck, affecting runtime performance. For this reason, efficient iterative solver strategies are required to effectively solve the linear systems generated by numerical discretization.

### Paper contributions

To solve the above problems, the authors made the following contributions:

1. **New numerical free - surface INSE model**: Developed a new numerical free - surface incompressible Navier - Stokes solver based on the high - order nodal spectral element method, which is suitable for wave propagation and wave - structure interaction applications.
2. **Geometric p - multigrid acceleration**: Proposed a geometric p - multigrid method to accelerate the solution of the Poisson problem, achieving O(n) computational scalability, making this method suitable for extension to more complex real - world scenarios.
3. **Verification and experimental comparison**: Verified through numerical experiments that this model achieves spectral convergence in the velocity field of highly nonlinear waves and has good agreement with the experimental data in classical benchmark tests.

Through these improvements, this model can significantly reduce the computational cost while ensuring high accuracy, providing new tools and methods for nonlinear water wave simulation.