A GPU accelerated level set reinitialization for an adaptive discontinuous Galerkin method

2016-08-01
KARAKUS, A.
WARBURTON, T.
AKSEL, MEHMET HALUK
Sert, Cüneyt
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.
COMPUTERS & MATHEMATICS WITH APPLICATIONS

Suggestions

A finite element variational multiscale method for the Navier-Stokes equations
Volker, John; Kaya Merdan, Songül (Society for Industrial & Applied Mathematics (SIAM), 2005-01-01)
This paper presents a variational multiscale method (VMS) for the incompressible Navier-Stokes equations which is defined by a large scale space L-H for the velocity deformation tensor and a turbulent viscosity nu(T). The connection of this method to the standard formulation of a VMS is explained. The conditions on L-H under which the VMS can be implemented easily and efficiently into an existing finite element code for solving the Navier - Stokes equations are studied. Numerical tests with the Smagorinsky ...
Interval oscillation criteria for second-order forced delay dynamic equations with mixed nonlinearities
Agarwal, Ravi P.; Anderson, Douglas R.; Zafer, Ağacık (Elsevier BV, 2010-01-01)
Interval oscillation criteria are established for second-order forced delay dynamic equations on time scales containing mixed nonlinearities of the form
A parallel multithreaded sparse triangular linear system solver
Cugu, Ilke; Manguoğlu, Murat (Elsevier BV, 2020-07-15)
We propose a parallel sparse triangular linear system solver based on the Spike algorithm. Sparse triangular systems are required to be solved in many applications. Often, they are a bottleneck due to their inherently sequential nature. Furthermore, typically many successive systems with the same coefficient matrix and with different right hand side vectors are required to be solved. The proposed solver decouples the problem at the cost of extra arithmetic operations as in the banded case. Compared to the b...
A nested iterative scheme for computation of incompressible flows in long domains
Manguoğlu, Murat; Tezduyar, Tayfun E.; Sathe, Sunil (Springer Science and Business Media LLC, 2008-12-01)
We present an effective preconditioning technique for solving the nonsymmetric linear systems encountered in computation of incompressible flows in long domains. The application category we focus on is arterial fluid mechanics. These linear systems are solved using a nested iterative scheme with an outer Richardson scheme and an inner iteration that is handled via a Krylov subspace method. Test computations that demonstrate the robustness of our nested scheme are presented.
The differential quadrature solution of nonlinear reaction-diffusion and wave equations using several time-integration schemes
Meral, Gulnihal; Tezer, Münevver (Wiley, 2011-04-01)
Three different time-integration schemes, namely the finite difference method (FDM) with a relaxation parameter, the least-squares method (LSM) and the finite element method (FEM), are applied to the differential quadrature (DQM) solution of one-dimensional nonlinear reaction-diffusion and wave equations. In the solution procedure, the space derivatives are discretized using DQM, which may also be used without the need of boundary conditions. The aim of the paper is to find computationally more efficient ti...
Citation Formats
A. KARAKUS, T. WARBURTON, M. H. AKSEL, and C. Sert, “A GPU accelerated level set reinitialization for an adaptive discontinuous Galerkin method,” COMPUTERS & MATHEMATICS WITH APPLICATIONS, pp. 755–767, 2016, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/41839.