OpenCL-based implementation of an unstructured edge-based finite element convection-diffusion solver on graphics hardware
Autor: | Farshid Mossaiby, Riccardo Rossi, Pooyan Dadvand, Sergio Idelsohn |
---|---|
Rok vydání: | 2011 |
Předmět: |
Numerical Analysis
Discretization Computer science business.industry Applied Mathematics Graphics hardware Numerical analysis General Engineering Finite difference Parallel computing Computational fluid dynamics Solver Data structure 01 natural sciences Finite element method 010305 fluids & plasmas Computational science Computer Science::Performance 010101 applied mathematics 0103 physical sciences 0101 mathematics business |
Zdroj: | International Journal for Numerical Methods in Engineering. 89:1635-1651 |
ISSN: | 0029-5981 |
DOI: | 10.1002/nme.3302 |
Popis: | SUMMARY The solution of problems in computational fluid dynamics (CFD) represents a classical field for the application of advanced numerical methods. Many different approaches were developed over the years to address CFD applications. Good examples are finite volumes, finite differences (FD), and finite elements (FE) but also newer approaches such as the lattice-Boltzmann (LB), smooth particle hydrodynamics or the particle finite element method. FD and LB methods on regular grids are known to be superior in terms of raw computing speed, but using such regular discretization represents an important limitation in dealing with complex geometries. Here, we concentrate on unstructured approaches which are less common in the GPU world. We employ a nonstandard FE approach which leverages an optimized edge-based data structure allowing a highly parallel implementation. Such technique is applied to the ‘convection-diffusion’ problem, which is often considered as a first step towards CFD because of similarities to the nonconservative form of the Navier–Stokes equations. In this regard, an existing highly optimized parallel OpenMP solver is ported to graphics hardware based on the OpenCL platform. The optimizations performed are discussed in detail. A number of benchmarks prove that the GPU-accelerated OpenCL code consistently outperforms the OpenMP version. Copyright © 2011 John Wiley & Sons, Ltd. |
Databáze: | OpenAIRE |
Externí odkaz: |