Numerical Stability and Catalan Numbers
Abstract
To predict allowable time-step size for the fully discretized nonlinear differential equations, a stability theory is developed using exact determination of an infinite perturbation series. Mathematical induction is used to determine the coefficients of the series. It is discovered that the closed-form equation for the nonlinear shift of generic polynomial non-linearity can be written as a series expansion where the coefficients are the Pfaff-Fuss-Catalan numbers in Combinatorics. This reveals criteria which can be used to analytically determine the allowable time step. It is shown that stability region decreases when the nonlinearity of the differential equation increases. Therefore, the maximum allowable time step is severely limited by the nonlinearity even if an unconditionally stable scheme (in a linear sense) is used. The theory is applied to general system of time-dependent nonlinear Partial Differential Equations.
- Publication:
-
arXiv e-prints
- Pub Date:
- September 2013
- DOI:
- 10.48550/arXiv.1309.4820
- arXiv:
- arXiv:1309.4820
- Bibcode:
- 2013arXiv1309.4820G
- Keywords:
-
- Mathematics - Numerical Analysis;
- Mathematics - Analysis of PDEs;
- Mathematics - Combinatorics;
- 65M12 65L20 39B82 34K20 35B35
- E-Print:
- 18 pages, 4 figures - Details of derivations are moved to appendices. Several typos are fixed. A comprehensive example for a nonlinear PDE is added to the end of the paper