A Stabilised Nodal Spectral Element Method for Fully Nonlinear Water Waves

We present an arbitrary-order spectral element method for general-purpose simulation of non-overturning water waves, described by fully nonlinear potential theory. The method can be viewed as a high-order extension of the classical finite element method proposed by Cai et al (1998) \cite{CaiEtAl1998}, although the numerical implementation differs greatly. Features of the proposed spectral element method include: nodal Lagrange basis functions, a general quadrature-free approach and gradient recovery using global $L^2$ projections. The quartic nonlinear terms present in the Zakharov form of the free surface conditions can cause severe aliasing problems and consequently numerical instability for marginally resolved or very steep waves. We show how the scheme can be stabilised through a combination of over-integration of the Galerkin projections and a mild spectral filtering on a per element basis. This effectively removes any aliasing driven instabilities while retaining the high-order accuracy of the numerical scheme. The additional computational cost of the over-integration is found insignificant compared to the cost of solving the Laplace problem. The model is applied to several benchmark cases in two dimensions. The results confirm the high order accuracy of the model (exponential convergence), and demonstrate the potential for accuracy and speedup. The results of numerical experiments are in excellent agreement with both analytical and experimental results for strongly nonlinear and irregular dispersive wave propagation. The benefit of using a high-order -- possibly adapted -- spatial discretization for accurate water wave propagation over long times and distances is particularly attractive for marine hydrodynamics applications.

Comments: Accepted for publication in Journal of Computational Physics April 29, 2016

Similar Publications

In the standard SPH method, the interaction between two particles might be not pairwise when the support domain varies, which can result in a reduction of accuracy. To deal with this problem, a modified SPH approach is presented in this paper. First of all, a Lagrangian kernel is introduced to eliminate spurious distortions of the domain of material stability, and the gradient is corrected by a linear transformation so that linear completeness is satisfied. Read More


It is generally accepted that all models are wrong -- the difficulty is determining which are useful. Here, a useful model is considered as one that is capable of combining data and expert knowledge, through an inversion or calibration process, to adequately characterize the uncertainty in predictions of interest. This paper derives conditions that specify which simplified models are useful and how they should be calibrated. Read More


We investigate the transition to a Landau-Levich-Derjaguin film in forced dewetting using a quadtree adaptive solution to the Navier-Stokes equations with surface tension. A discretization of the capillary forces near the receding contact line is used that yields an equilibrium for a specified contact angle $\theta_\Delta$ called the numerical contact angle. Despite the well-known contact line singularity, dynamic simulations can proceed without any explicit additional numerical procedure, yielding an implicitly dynamic contact angle model. Read More


Ab initio quantum chemistry calculations for systems with large active spaces are notoriously difficult and cannot be successfully tackled by standard methods. In this letter, we generalize a Green's function QM/QM embedding method called self-energy embedding theory (SEET) that has the potential to be successfully employed to treat large active spaces. In generalized SEET, active orbitals are grouped into intersecting groups of few orbitals allowing us to perform multiple parallel calculations yielding results comparable to the full active space treatment. Read More


Drag laws for particles in fluids are often expressed in terms of the undisturbed fluid velocity, defined as the fluid velocity a particle sees before the disturbance develops in the fluid. In two-way coupled point-particle simulations the information from the undisturbed state is not available and must be approximated using the disturbed velocity field. Horwitz and Mani (2016) recently developed a procedure to estimate the undisturbed velocity for particles moving at low Reynolds number and obeying the Stokes drag law. Read More


The Wannier localization problem in quantum physics is mathematically analogous to finding a localized representation of a subspace corresponding to a nonlinear eigenvalue problem. While Wannier localization is well understood for insulating materials with isolated eigenvalues, less is known for metallic systems with entangled eigenvalues. Currently, the most widely used method for treating systems with entangled eigenvalues is to first obtain a reduced subspace (often referred to as disentanglement) and then to solve the Wannier localization problem by treating the reduced subspace as an isolated system. Read More


We provide the necessary framework for carrying out stochastic positive-P and gauge-P simulations of bosonic systems with long range interactions. In these approaches, the quantum evolution is sampled by trajectories in phase space, allowing calculation of correlations without truncation of the Hilbert space or other approximations to the quantum state. The main drawback is that the simulation time is limited by noise arising from interactions. Read More


A new highly efficient method is developed for computation of traveling periodic waves (Stokes waves) on the free surface of deep water. A convergence of numerical approximation is determined by the complex singularites above the free surface for the analytical continuation of the travelling wave into the complex plane. An auxiliary conformal mapping is introduced which moves singularities away from the free surface thus dramatically speeding up numerical convergence by adapting the numerical grid for resolving singularities while being consistent with the fluid dynamics. Read More


Machine learning has proven to be a valuable tool to approximate functions in high-dimensional spaces. Unfortunately, analysis of these models to extract the relevant physics is never as easy as applying machine learning to a large dataset in the first place. Here we present a description of atomic systems that generates machine learning representations with a direct path to physical interpretation. Read More


We developed a combined atomistic-continuum hierarchical multiscale approach to explore the effective thermal conductivity of graphene laminates. To this aim, we first performed molecular dynamics simulations in order to study the heat conduction at atomistic level. Using the non-equilibrium molecular dynamics method, we evaluated the length dependent thermal conductivity of graphene as well as the thermal contact conductance between two individual graphene sheets. Read More