Comparison of Five Methods of Computing the Dirichlet–Neumann Operator for the Water Wave Problem

Abstract

Abstract. We compare the effectiveness of solving Dirichlet-Neumann problems via the Craig-Sulem (CS) expansion, the Ablowitz-Fokas-Musslimani (AFM) implicit formulation, the dual AFM formulation (AFM * ), a boundary integral collocation method (BIM), and the transformed field expansion (TFE) method. The first three methods involve highly ill-conditioned intermediate calculations that we show can be overcome using multiple-precision arithmetic. The latter two methods avoid catastrophic cancellation of digits in intermediate results, and are much better suited to numerical computation.For the Craig-Sulem expansion, we explore the cancellation of terms at each order (up to 150th) for three types of wave profiles, namely band-limited, realanalytic, or smooth. For the AFM and AFM * methods, we present an example in which representing the Dirichlet or Neumann data as a series using the AFM basis functions is impossible, causing the methods to fail. The example involves band-limited wave profiles of arbitrarily small amplitude, with analytic Dirichlet data. We then show how to regularize the AFM and AFM * methods by over-sampling the basis functions and using the singular value decomposition or QR-factorization to orthogonalize them. Two additional examples are used to compare all five methods in the context of water waves, namely a largeamplitude standing wave in deep water, and a pair of interacting traveling waves in finite depth.

Stability of periodic traveling flexural‐gravity waves in two dimensions

“…There are several reformulations that exist for the given set of Euler's equations (for a review, see Ref. 27). In this work, the focus is on the implicit formulation that we will refer to as AFM formulation.…”

A three-dimensional Dirichlet-to-Neumann operator for water waves over topography

“…A very detailed comparison between five different Dirichlet-to-Neumann formulations is provided by Wilkening & Vasan [29]. They compare expansion methods against methods for solving the non-local equations directly.…”

Um Método Integral de Contorno para a Modelagem da Propagação de Ondas Internas em um Sistema de Dois Fluidos

“…A aproximação eficiente e precisa desse operador tem representado um dos maiores desafios nesses trabalhos. No trabalho recente [16]é apresentada uma comparação de vários métodos desenvolvidos com esse objetivo.…”

“…Os métodos baseados na formulação do operador Dirichlet-Neumann através da integral de contorno possuem boa precisão e estabilidade numérica [16]. Uma aplicação recente desse método foi apresentada em [9], para fazer um estudo numérico dasórbitas de partículas associadas com ondas de Stokes progressivas sobre a influência de uma correnteza uniforme.…”