Remote Access Mathematics of Computation
Green Open Access

Mathematics of Computation

ISSN 1088-6842(online) ISSN 0025-5718(print)

 
 

 

A necessary condition for the stability of a difference approximation to a hyperbolic system of partial differential equations


Author: Anne M. Burns
Journal: Math. Comp. 32 (1978), 707-724
MSC: Primary 65M10; Secondary 65N10
DOI: https://doi.org/10.1090/S0025-5718-1978-0492034-2
MathSciNet review: 492034
Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: We are interested in the boundary conditions for a difference approximation to a hyperbolic system of partial differential equations $ {u_t} = A{u_x}$, $ u(x,0) = F(x)$, $ Ku(0,t) = 0$ in the quarter plane $ x \geqslant 0$, $ t \geqslant 0$.

We consider approximations of the type:

$\displaystyle {u_j}(t + \Delta t) = \sum\limits_{k = - r}^p {{C_k}{u_{j + k}}(t),\quad j = 1,2, \ldots ,} $

$\displaystyle {u_j} + \sum\limits_{k = 1}^s {{\alpha _{jk}}{u_k}(t + \Delta t) = 0,\quad j = - r + 1, \ldots ,0.} $

If N is the null space of K and E is the "negative" eigenspace of A, then the system of partial differential equations is well-posed if and only if $ K \cap E = \left\{ 0 \right\}$ and Rank $ K = $ the number of negative eigenvalues of A.

In direct analogy to this, we prove that for a difference scheme of the above type with $ r = p = 1$, $ K\prime = I + \Sigma _{k = 1}^s\;{\alpha _k}$ and $ N\prime = $ null space of $ K\prime$, a necessary condition for stability is $ N\prime \cap E = \left\{ 0 \right\}$. If, in addition, a condition proven by S. J. Osher to be sufficient for stability is not satisfied, then Rank $ K = $ the number of negative eigenvalues of A is also necessary for stability. We then generalize this result to the case $ r > 1$, $ p > 1$.

Together these conditions imply that "extrapolation" on "negative" eigenvectors leads to instability; "extrapolation" on "positive" eigenvectors is "almost necessary. "Extrapolation" on "positive" eigenvectors and not on "negative" eigenvectors is sufficient for stability.


References [Enhancements On Off] (What's this?)

  • [1] S. K. GODUNOV & V. S. RYABENKI, Theory of Difference Schemes, An Introduction, North-Holland, Amsterdam, 1964. MR 0181117 (31:5346)
  • [2] REUBEN HERSH, "Mixed problems in several variables," J. Math. Mech., v. 12, 1963. MR 0147790 (26:5304)
  • [3] H. O. KREISS, "Stability theory for difference approximations of mixed initial boundary value problems. I," Math. Comp., v. 22, 1968, pp. 703-714. MR 0241010 (39:2355)
  • [4] H. O. KREISS & J. OLIGER, Methods for the Approximate Solution of Time Dependent Problems, GARP Publications Series no. 10, Feb. 1973.
  • [5] PETER D. LAX & BURTON WENDROFF, "Difference schemes for hyperbolic equations," Comm. Pure Appl. Math., v. 17, 1964. MR 0170484 (30:722)
  • [6] S. J. OSHER, "Systems of difference equations with general homogeneous boundary conditions," Trans. Amer. Math. Soc., v. 137, 1969, pp. 177-201. MR 0237982 (38:6259)
  • [7] S. J. OSHER, "On systems of difference equations with wrong boundary values," Math. Comp., v. 23, 1969, pp. 567-572. MR 0247785 (40:1046)
  • [8] S. J. OSHER, "On certain Toeplitz operators in two variables," Pacific J. Math., v. 34, 1970, pp. 123-129. MR 0267408 (42:2310)

Similar Articles

Retrieve articles in Mathematics of Computation with MSC: 65M10, 65N10

Retrieve articles in all journals with MSC: 65M10, 65N10


Additional Information

DOI: https://doi.org/10.1090/S0025-5718-1978-0492034-2
Article copyright: © Copyright 1978 American Mathematical Society

American Mathematical Society