Premium
Avoiding spurious modes of time discretized operators in transport problems
Author(s) -
Arruda Natália C.,
Almeida R. C.,
Silva R. S.,
Malta Sandra M. C.
Publication year - 2010
Publication title -
international journal for numerical methods in biomedical engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.741
H-Index - 63
eISSN - 2040-7947
pISSN - 2040-7939
DOI - 10.1002/cnm.1159
Subject(s) - spurious relationship , discretization , galerkin method , singularity , nonlinear system , mathematics , stability (learning theory) , petrov–galerkin method , finite element method , mathematical analysis , range (aeronautics) , physics , computer science , engineering , statistics , quantum mechanics , machine learning , thermodynamics , aerospace engineering
The solution of the Crank–Nicolson (CN) method may present spurious modes, which mostly happens in systems with an initial singularity. In this work, we analyze the onset of spurious oscillations for the CN method for large time steps. This analysis is performed on the general one‐dimensional linear advection–diffusion–reaction transport equation with initial singularity that models a broad range of mass transfer phenomena in natural and engineering sciences. The discrete problem is obtained using a stable finite element method in space and the generalized trapezoidal family of methods in time. Depending on the range of parameters, either the Galerkin or the streamline upwind Petrov–Galerkin methods are used to guarantee stability in space. We derive a stability analysis of the fully discrete method and investigate the techniques proposed in the literature to damp oscillations. We propose a new threshold of time oscillations to overcome the spurious modes, which is also extended to deal with nonlinear problems. Copyright © 2008 John Wiley & Sons, Ltd.