Fourier spectral/hp element method: investigation of time-stepping and parallelisation strategies
Fourier spectral/hp element method: investigation of time-stepping and parallelisation strategies
复制标题
DOI：
10.25560/25140
复制
发表时间：
2013-09
期刊：
影响因子：
--
通讯作者：
A. Bolis
中科院分区：
文献类型：
--
作者：
A. Bolis
As computer hardware has evolved, the time required to perform numerical simulations has reduced, allowing investigations of a wide range of new problems. This thesis focuses on algorithm optimisation, to minimise run-time, when solving the incompressible Navier-Stokes equations. Aspects affecting performance related to the discretisation and algorithm parallelisation are investigated in the context of high-order methods. The roles played by numerical approximations and computational strategies are highlighted and it is recognised that a versatile implementation provides additional benefits, allowing an ad-hoc selection of techniques to fit the needs of heterogeneous computing environments. We initially describe the building blocks of a spectral/hp element and pure spectral method and how they can be encapsulated and combined to create a 3D discretisation, the Fourier spectral/hp element method. Time-stepping strategies are also described and encapsulated in a flexible framework based on the General Linear Method. After implementing and validating an incompressible Navier-Stokes solver, two canonical turbulent flows are analysed. Afterward a 2D hyperbolic equation is considered to investigate the efficiency of lowand high-order methods when discretising the spatial and temporal derivatives. We perform parametric studies, monitoring accuracy and CPU-time for different numerical approximations. We identify optimal discretisations, demonstrating that high-order methods are the computationally fastest approach to attain a desired accuracy for this problem. Following the same philosophy, we investigate the benefits of using a hybrid parallel implementation. The message passing model is introduced to parallelise different kernels of an incompressible Navier-Stokes solver. Monitoring the parallel performance of these strategies the most efficient approach is highlighted. We also demonstrate that hybrid parallel solutions can be used to significantly extend the strong scalability limit and support greater parallelism.
