首页 | 本学科首页   官方微博 | 高级检索  
     检索      


Numerical Stability in the Presence of Variable Coefficients
Authors:Ernst Hairer  Arieh Iserles
Institution:1.Section de mathématiques,Université de Genève,Geneva,Switzerland;2.Department of Applied Mathematics and Theoretical Physics, Centre for Mathematical Sciences,University of Cambridge,Cambridge,UK
Abstract:The main concern of this paper is with the stable discretisation of linear partial differential equations of evolution with time-varying coefficients. We commence by demonstrating that an approximation of the first derivative by a skew-symmetric matrix is fundamental in ensuring stability for many differential equations of evolution. This motivates our detailed study of skew-symmetric differentiation matrices for univariate finite-difference methods. We prove that, in order to sustain a skew-symmetric differentiation matrix of order \(p\ge 2\), a grid must satisfy \(2p-3\) polynomial conditions. Moreover, once it satisfies these conditions, it supports a banded skew-symmetric differentiation matrix of this order and of the bandwidth \(2p-1\), which can be derived in a constructive manner. Some applications require not just skew-symmetry, but also that the growth in the elements of the differentiation matrix is at most linear in the number of unknowns. This is always true for our tridiagonal matrices of order 2 but need not be true otherwise, a subject which we explore further. Another subject which we examine is the existence and practical construction of grids that support skew-symmetric differentiation matrices of a given order. We resolve this issue completely for order-two methods. We conclude the paper with a list of open problems and their discussion.
Keywords:
本文献已被 SpringerLink 等数据库收录!
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号