• Media type: Text; Report; E-Book
  • Title: Positivity preservation of implicit discretizations of the advection equation
  • Contributor: Hadjimichael, Yiannis [Author]; Ketcheson, David I. [Author]; Lóczi, Lajos [Author]
  • Published: Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2021
  • Issue: published Version
  • Language: English
  • DOI: https://doi.org/10.34657/8602; https://doi.org/10.20347/WIAS.PREPRINT.2846
  • ISSN: 2198-5855
  • Keywords: finite difference ; implicit time-discretization ; spectral collocation ; positivity preservation ; linear partial differential equations
  • Origination:
  • Footnote: Diese Datenquelle enthält auch Bestandsnachweise, die nicht zu einem Volltext führen.
  • Description: We analyze, from the viewpoint of positivity preservation, certain discretizations of a fundamental partial differential equation, the one-dimensional advection equation with periodic boundary condition. The full discretization is obtained by coupling a finite difference spatial semidiscretization (the second- and some higher-order centered difference schemes, or the Fourier spectral collocation method) with an arbitrary _x0012_θ-method in time (including the forward and backward Euler methods, and a second-order method by choosing _x0012_ θ ∈ [0, 1] suitably). The full discretization generates a two-parameter family of circulant matrices M ∈ ℝ m_x0002_xm , where each matrix entry is a rational function in θ and _x0017_ν . Here, _x0017_ν denotes the CFL number, being proportional to the ratio between the temporal and spatial discretization step sizes. The entrywise non-negativity of the matrix M---which is equivalent to the positivity preservation of the fully discrete scheme---is investigated via discrete Fourier analysis and also by solving some low-order parametric linear recursions. We find that positivity preservation of the fully discrete system is impossible if the number of spatial grid points m is even. However, it turns out that positivity preservation of the fully discrete system is recovered for odd values of m provided that θ ≥ 1/2 and ν are chosen suitably. These results are interesting since the systems of ordinary differential equations obtained via the spatial semi-discretizations studied are not positivity preserving.
  • Access State: Open Access