High accuracy periodic solutions to the Sivashinsky equation

Authors:
V. Karlin;V. Maz'ya;G. Schmidt
Affiliations:
Centre for Research in Fire and Explosions, University of Central Lancashire, Preston PR1 2HE, UK;Department of Mathematics, Linköping University, Linköping S-581 83, Sweden;Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstrasse 39, Berlin D-10117, Germany
Venue:
Journal of Computational Physics
Year:
2003

Citing 7
Cited 0

The cellular nature of hydrodynamic flame instability

Physica D
Computational frameworks for the fast Fourier transform

Computational frameworks for the fast Fourier transform
Time-Marching Algorithms for Nonlocal Evolution Equations Based Upon "Approximate Approximations"

SIAM Journal on Scientific Computing
Pseudospectra of Linear Operators

SIAM Review
Numerical solution of nonlinear hypersingular integral equations of the Peierls type in dislocation theory

SIAM Journal on Applied Mathematics
Stability of pole solutions for planar propagating flames: I. exact eigenvalues and eigenfunctions

SIAM Journal on Applied Mathematics
Handbook of Mathematical Functions, With Formulas, Graphs, and Mathematical Tables,

Handbook of Mathematical Functions, With Formulas, Graphs, and Mathematical Tables,

Quantified Score

Hi-index	31.45

Visualization

Abstract

The aim of this work is the accurate calculation of periodic solutions to the Sivashinsky equation, which models dynamics of the long wave instability of laminar premixed flame. A highly accurate computational algorithm was developed in both one and two spatial dimensions and its crucial implementation details are presented. The algorithm is based on the concept of saturated asymptotic approximations and can be straightforwardly extended to a wide variety of nonlinear integro-differential equations. The development of such an algorithm was motivated by difficulties in interpretation of the results of numerical experiments with the Sivashinsky equation using spectral methods. The computations carried out by the algorithm in question are in good agreement with the results obtained earlier by spectral methods. Analysis of the accuracy of obtained numerical solutions and of their stabilization to steady states supports the idea of the instability of the steady coalescent pole solutions (with maximal possible number of poles) to the Sivashinsky equation in large domains through huge linear transient amplification of nonmodal perturbations of small but finite amplitudes.