A high‐order PIC method for advection‐dominated flow with application to shallow water waves

Abstract

Summary A hybrid Eulerian‐Lagrangian particle‐in‐cell–type numerical method is developed for the solution of advection‐dominated flow problems. Particular attention is given over to the high‐order transfer of flow properties from the particles to the grid. For smooth flows, the method presented is of formal high‐order accuracy in space. The method is applied to solve the nonlinear shallow water equations resulting in a new, and novel, shock capturing shallow water solver. The approach is able to simulate complex shallow water flows, which can contain an arbitrary number of discontinuities. Both trivial and nontrivial bottom topography is considered, and it is shown that the new scheme is inherently well balanced, exactly satisfying the scriptC‐property. The scheme is verified against several one‐dimensional benchmark shallow water problems. These include cases that involve transcritical flow regimes, shock waves, and nontrivial bathymetry. In all the test cases presented, very good results are obtained.

On the hydrodynamic performance of a vertical pile-restrained WEC-type floating breakwater

“…Further developments can be found in, for example, Brackbill and Ruppel [10] and Brackbill et al [11], which significantly reduce the numerical dissipation of PIC method. Recently, variations of the PIC method have achieved high-order accuracy (see Edwards and Bridson [12], Maljaars et al [13] and Wang and Kelly [14]).…”

A 3D parallel particle-in-cell solver for extreme wave interaction with floating bodies

“…Later, many attempts have been made to improve this method (see Brackbill and Ruppel (1986); Brackbill et al (1988); Yabe (1982, 1983)). More recently, high-order PIC variations have become possible (see Edwards and Bridson (2012), Maljaars et al (2018) and Wang and Kelly (2018)). However, so far this hybrid method has not been very well exploited for use in the ocean engineering field, where modelling complex wavestructure interaction with computational efficiency still remains an important challenge.…”