## Numerical approximation of the spectra of non-compact operators arising in buckling problems

*
Monique Dauge,
Manil Suri
*

The finite element method approximates the spectrum of
an operator *S* by computing the spectra of a sequence of operators
*S*_{N}
defined in terms of the finite element spaces. For the case that
*S* is compact, convergence of the approximate spectra follows from
the convergence of *S*_{N} to
*S*
in the operator norm.

We consider the case that *S* is non-compact, in
which case such operator norm convergence cannot take place, and the
approximations may be polluted by spurious eigenvalues. Pollution-free
convergence of the eigenvalues can, however, be
guaranteed outside the **essential numerical range** of
*S*, which is
related to the essential spectrum of *S*.

We present results for estimating this essential numerical range and apply them
to an algorithm for the buckling of
three-dimensional bodies (that gives rise to a non-compact *S*).
Our results show, for instance, that for the example of a circular disc,
the algorithm will be free of spurious eigenvalues provided the body is thin enough.
The case that singularities in the stresses can lead to
non-physical spectral values being approximated is also investigated.

The analysis in this paper is made complete in the Next Paper
where we prove that under generic assumptions thin plates produce sequences of buckling eigenvalues
which tend to 0 as the thickness tends to 0, thus escaping the essential spectrum region.

Rennes, Baltimore, November 2001.

Preprint Rennes 01-54

*Journal of Numerical Mathematics * **10**, 3, 193--219, 2002.