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

We present in this work a new reconstruction scheme, so-called MUSCL-THINC-BVD scheme, to solve the five-equation model for interfacial two phase flows. This scheme employs the traditional shock capturing MUSCL (Monotone Upstream-centered Schemes for Conservation Law) scheme as well as the interface sharpening THINC (Tangent of Hyperbola for INterface Capturing) scheme as two building-blocks of spatial reconstruction using the BVD (boundary variation diminishing) principle that minimizes the variations (jumps) of the reconstructed variables at cell boundaries, and thus effectively reduces the numerical dissipations in numerical solutions. The MUSCL-THINC-BVD scheme is implemented to all state variables and volume fraction, which realizes the consistency among volume fraction and other physical variables. Read More


Laterally periodic nanostructures were investigated with grazing incidence small angle X-ray scattering (GISAXS) by using the diffraction patterns to reconstruct the surface shape. To model visible light scattering, rigorous calculations of the near and far field by numerically solving Maxwell's equations with a finite-element method are well established. The application of this technique to X-rays is still challenging, due to the discrepancy between incident wavelength and finite-element size. Read More


This paper proposes the use of the Spectral method to simulate diffusive moisture transfer through porous materials, which can be strongly nonlinear and can significantly affect sensible and latent heat transfer. An alternative way for computing solutions by considering a separated representation is presented, which can be applied to both linear and nonlinear diffusive problems, considering highly moisture-dependent properties. The Spectral method is compared with the classical implicit Euler and Crank-Nicolson schemes. Read More


We have developed the first laser damage simulation algorithm capable of determining crater and surface modification morphology from microscopic physics. Rapid progress in the field of high intensity ultrafast lasers has enabled its utility in a myriad of applications. Simulation plays an important role in this research by allowing for closer analysis of the physical mechanisms involved, but current techniques struggle to meet both the spatial scope or resolution requirements for modeling such dynamics, typically specializing in one or the other. Read More


A multi-scale framework was recently proposed for more realistic molecular dynamics simulations in continuum solvent models by coupling a molecular mechanics treatment of solute with a fluid mechanics treatment of solvent, where we formulated the physical model and developed a numerical fluid dynamics integrator. In this study, we incorporated the fluid dynamics integrator with the Amber simulation engine to conduct atomistic simulations of biomolecules. At this stage of the development, only nonelectrostatic interactions, i. Read More


A self-consistent procedure for ion-induced soft error rate calculation in space environment taking into account Error Correcting Codes is proposed. The method is based on partitioning of multiple cell events into groups with different multiplicities. Read More


A procedure of space proton-induced Soft Error Rate calculation based on heavy ion testing data is proposed. The approach relies on Geant4-assisted Monte Carlo simulation of the secondary particle LET spectra produced by proton interactions. Read More


We recalculate the leading relativistic corrections for the ground electronic state of the hydrogen molecule using variational method with explicitly correlated functions which satisfy the interelectronic cusp condition. The new computational approach allowed for the control of the numerical precision which reached about 8 significant digits. More importantly, the updated theoretical energies became discrepant with the known experimental values and we conclude that the yet unknown relativistic recoil corrections might be larger than previously anticipated. Read More


We evaluate the practical performance of the zero-multiple summation method (ZMM), a method for approximately calculating electrostatic interactions in molecular dynamics simulations. The performance of the ZMM is compared with that of the smooth particle mesh Ewald method (SPME). Even though the ZMM uses a larger cutoff distance than the SPME does, the performance of the ZMM is found to be comparable to or better than that of the SPME. Read More


There is considerable interest in collective effects in hybrid systems formed by molecular or atomic ensembles strongly coupled by an electromagnetic resonance. For analyzing such collective effects, we develop an efficient and general theoretical formalism based on the natural modes of the resonator. The main strength of our approach is its generality and the high level of analyticity enabled by modal analysis, which allows one to model complex hybrid systems without any restriction on the resonator shapes or material properties, and to perform statistical computations to predict general properties that are robust to spatial and polarization disorders. Read More