Stabilized variational formulation of an oldroyd-B fluid flow equations on a Graphic Processing Unit (GPU) architecture

The governing equations of the flow of an oldroyd-B fluid are discretized using the finite element method. To overcome the convective nature of the momentum equation, the Galerkin/Least-Squares Finite Element Method (GLS/FEM) is used while the Discrete Elastic-Viscous Stress-Splitting (DEVSS) method is used to overcome the instability due to the absence of diffusion in the constitutive equations. The discretized equations are implemented on a hybrid system between the Graphics Processing Unit (GPU) architecture using Compute-Unified-Device-Architecture (CUDA) and a multi-core CPU. The implementation is applied successfully to simulate the blood flow in abdominal aortic aneurysm. To accelerate application performance on the GPU several optimized approaches are adopted. The most significant approach is the coloring technique that is used to assemble the global matrix. Numerical experiments show that the hybrid CPU-GPU implementation has a 26 time speedup over the multi-core CPU implementations. (C) 2020 Published by Elsevier B.V.

Automated summary

The governing equations of the flow of an Oldroyd-B fluid are discretized using the finite element method, and the GLS/FEM method is used to handle the convective nature of the momentum equation, while the DEVSS method is used to overcome instability caused by the lack of diffusion in the constitutive equations. The discretized equations are successfully implemented on a hybrid system between GPU architecture and multi-core CPU, accelerating application performance by adopting optimized approaches, leading to a 26 times speedup over multi-core CPU implementations.