A GPU accelerated level set reinitialization for an adaptive discontinuous Galerkin method. (August 2016)

• Record Type: Journal Article
• Authors: Karakus, A.; Warburton, T.; Aksel, M.H.; Sert, C.
• Is Part Of: Computers & mathematics with applications. Volume 72:issue 3(2016)
• Page Start: 755
• Page End: 767
• Publication Date: 2016-08
• Subjects: Discontinuous Galerkin -- Level set -- Local time integration -- Hamilton–Jacobi -- GPU computing; Electronic data processing -- Periodicals; Mathematics -- Data processing -- Periodicals; 510.28541
• DOI: 10.1016/j.camwa.2016.05.025
• Languages: English
• ISSNs: 0898-1221

Abstract

Abstract: GPU accelerated high order reconstruction of signed distance function of the level set method is studied. The flow based reinitialization equation is discretized in space by using a nodal discontinuous Galerkin method on adaptive unstructured grids. Artificial diffusion with a modal decay rate based regularity estimator is used to damp out high frequency solution components near kinks, where mesh adaptivity is applied. A two rate Adams-Bashforth time integrator is developed to avoid time step restrictions resulting from artificial diffusion stabilization and local mesh refinement. Platform independence of the solver is achieved by using an extensible multi-threading programming API that allows runtime selection of different computing devices (GPU and CPU) and threading interfaces (CUDA, OpenCL and OpenMP). Overall, a highly scalable numerical scheme that preserves the simplicity of the original level set method is obtained. Performance and accuracy of the method to construct signed distance function on highly disturbed initial data with smooth and non-smooth interfaces are tested through distinct two- and three-dimensional problems.