Süleyman Cengizci, Ömür Uğur, Srinivasan Natesan, SUPG-Based Stabilized Finite Element Computations of Convection-Dominated 3D Elliptic PDEs using Shock-Capturing, Journal of Computational and Applied Mathematics, 451: 116022 (December 2024, May 2024).

PlumX

Abstract

This study is interested in stabilized finite element computations of advection-dominated convection–diffusion-type partial differential equations. The governing equations are assumed to be defined on the 3D unit cube and stationary. Towards that end, we stabilize the standard (Bubnov–) Galerkin finite element method (GFEM), which typically suffers from yielding numerical approximations polluted with node-to-node nonphysical oscillations in convection dominance, by employing the streamline-upwind/Petrov–Galerkin (SUPG) formulation. In order to achieve better shock representations near sharp layers, the SUPG-stabilized formulation is also enhanced with a simple and residual-based shock-capturing mechanism, the so-called YZβ technique. Several test computations are performed to assess the efficiency and robustness of the proposed formulation. The test problems are examined under more challenging conditions than those previously studied in the literature, i.e., for flows substantially dominated by convection phenomena. The numerical results and comparisons demonstrate that the SUPG-YZβ combination significantly eliminates both globally spread and localized numerical instabilities. Beyond that, the proposed formulation accomplishes that by using only linear elements on relatively coarser meshes and without making use of any adaptive mesh strategies.

Keywords: Convection-Diffusion, Convection-Dominated, Finite Elements, SUPG, YZβ Shock-Capturing

Orta Doğu Teknik Üniversitesi, Uygulamalı Matematik Enstitüsü, Üniversiteler Mahallesi, Dumlupınar Bulvarı No:1, 06800 Çankaya/Ankara