Eigenvector decomposition to determine the existence, shape, and location of numerical oscillations in Parabolic PDEs

by   Ethan Jensen, et al.

In this paper, we employed linear algebra and functional analysis to determine necessary and sufficient conditions for oscillation-free and stable solutions to linear and nonlinear parabolic partial differential equations. We applied singular value decomposition and Fourier analysis to various finite difference schemes to extract patterns in the eigenfunctions (sampled by the eigenvectors) and the shape of their eigenspectrum. With these, we determined how the initial and boundary conditions affect the frequency and long term behavior of numerical oscillations, as well as the location of solution regions most sensitive to them.


page 9

page 11


Bayesian Numerical Methods for Nonlinear Partial Differential Equations

The numerical solution of differential equations can be formulated as an...

Spectral methods for nonlinear functionals and functional differential equations

We develop a rigorous convergence analysis for finite-dimensional approx...

The numerical solution of semidiscrete linear evolution problems on the finite interval using the Unified Transform Method

We study a semidiscrete analogue of the Unified Transform Method introdu...

A well conditioned Method of Fundamental Solutions

The method of fundamental solutions (MFS) is a numerical method for solv...

On boundary conditions parametrized by analytic functions

Computer algebra can answer various questions about partial differential...

Extended Dynamic Mode Decomposition for Inhomogeneous Problems

Dynamic mode decomposition (DMD) is a powerful data-driven technique for...

Turnpike in optimal shape design

We introduce and study the turnpike property for time-varying shapes, wi...