High precision solutions of two fourth order eigenvalue problems

We solve the biharmonic eigenvalue problem Delta(2)u = lambda u and the buckling plate problem Delta(2)u = -lambda Delta u on the unit square using a highly accurate spectral Legendre-Galerkin method. We study the nodal lines for the first eigenfunction near a corner for the two problems. Five sign changes are computed and the results show that the eigenfunction exhibits a self similar pattern as one approaches the corner. The amplitudes of the extremal values and the coordinates of their location as measured from the corner are reduced by constant factors. These results are compared with the known asymptotic expansion of the solution near a corner. This comparison shows that the asymptotic expansion is highly accurate already from the first sign change as we have complete agreement between the numerical and the analytical results. Thus, we have an accurate description of the eigenfunction in the entire domain.