Parallel implementation of a non-hydrostatic model for free surface flows with semi-Lagrangian advection treatment

Jacek A. Jankowski

Abstract. The parallel implementation of an unstructured-grid, three-dimensional, semi-implicit finite difference and finite volume model for the free surface Navier-Stokes equations (UnTRIM) is presented and discussed. The new developments are aimed to make the code available for high performance computing in order to address larger, complex problems in environmental free surface flows. The parallelisation is based on the mesh partitioning method and message passing and has been achieved without negatively affecting any of the advantageous properties of the serial code, like its robustness, accuracy and efficiency. The key issue is a new, autonomous parallel streamline backtracking algorithm, which allows using semi-Lagrangian methods in decomposed meshes without compromising the scalability of the code. The implementation has been carefully verified not only with simple, abstract test cases illustrating the application domain of the code, but also with advanced, high resolution models presently applied for research and engineering projects. The scheme performance and accuracy aspects are researched and discussed.

The paper is available from the journal site (links below) and locally (password).

Reference: Jankowski, J.A. (2009) Parallel implementation of a non-hydrostatic model for free surface flows with semi-Lagrangian advection treatment. International Journal for Numerical Methods in Fluids, 59 (10), 1157–1179, 10 April 2009, published online 5 Aug 2008. doi:10.1002/fld.1859.