Eliminating Gibbs phenomena: A non-linear Petrov–Galerkin method for the convection–diffusion–reaction equation. (1st September 2020)

• Record Type: Journal Article
• Authors: Houston, Paul; Roggendorf, Sarah; van der Zee, Kristoffer G.
• Is Part Of: Computers & mathematics with applications. Volume 80:issue 5(2020)
• Page Start: 851
• Page End: 873
• Publication Date: 2020-09-01
• Subjects: Convection–diffusion -- Petrov–Galerkin -- Gibbs phenomenon -- Finite element methods -- Banach spaces; Electronic data processing -- Periodicals; Mathematics -- Data processing -- Periodicals; 510.28541
• DOI: 10.1016/j.camwa.2020.03.025
• Languages: English
• ISSNs: 0898-1221

Abstract

Abstract: In this article we consider the numerical approximation of the convection-diffusion-reaction equation. One of the main challenges of designing a numerical method for this problem is that boundary layers occurring in the convection-dominated case can lead to non-physical oscillations in the numerical approximation, often referred to as Gibbs phenomena. The idea of this article is to consider the approximation problem as a residual minimization in dual norms in L q -type Sobolev spaces, with 1 < q < ∞ . We then apply a non-standard, non-linear Petrov–Galerkin discretization, that is applicable to reflexive Banach spaces such that the space itself and its dual are strictly convex. Similar to discontinuous Petrov–Galerkin methods, this method is based on minimizing the residual in a dual norm. Replacing the intractable dual norm by a suitable discrete dual norm gives rise to a non-linear inexact mixed method. This generalizes the Petrov–Galerkin framework developed in the context of discontinuous Petrov–Galerkin methods to more general Banach spaces. For the convection–diffusion–reaction equation, this yields a generalization of a similar approach from the L 2 -setting to the L q -setting. A key advantage of considering a more general Banach space setting is that, in certain cases, the oscillations in the numerical approximation vanish as q tends to 1, as we will demonstrate using a few simple numerical examples.