Eliminating Gibbs phenomena: A non-linear Petrov-Galerkin method for the convection-diffusion-reaction equation

Houston, Paul; Roggendorf, Sarah; van der Zee, Kristoffer G. (2020). Eliminating Gibbs phenomena: A non-linear Petrov-Galerkin method for the convection-diffusion-reaction equation. Computers and mathematics with applications, 80(5), pp. 851-873. Elsevier 10.1016/j.camwa.2020.03.025

Preview	Text 1908.00996.pdf - Accepted Version Available under License Creative Commons: Attribution-Noncommercial-No Derivative Works (CC-BY-NC-ND). Download (5MB) | Preview
	Text 1-s2.0-S0898122120301292-main.pdf - Published Version Restricted to registered users only Available under License Publisher holds Copyright. Download (2MB) | Request a copy

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.

Interest & Impact

Downloads

50 since deposited on 10 Feb 2021
29 in the past 12 months

Citations

5 Citations in Web of Science ®
6 Citations in Scopus

Search

in Google Scholar™

Services

Actions (login required)

Edit item

Actions (login required)

Edit item