Stefano Giani Homesite

Convergence for Eigenvalue Problems with Mesh Adaptivity

Back to convergence page | Previous | Next

Numerics:

In this example we investigate how our method copes with discontinuous coefficients. We inserted a square subdomain of side 0.5 in the center of the unit square domain. We also choose the function $\mathcal{A}$ to be piecewise constant and to assume the value 100 inside the subdomain and the value 1 outside it.

The jump in the value of $\mathcal{A}$ could produce a jump in the gradient of the eigenfunctions all along the boundary of the subdomain. So the regularity of the eigenfunctions in each subdomain, in the sense of Assumption 1.1, is now between $3/2\leq s+1<2$ . Using uniform refinement, the rate of convergence for eigenvalues should be at least $\mathcal{O}(H^\mathrm{max}_n)^{2s}$ or equivalently $\mathcal{O}(N^{-s})$ , where is the number of DOFs. Instead, using our method we obtain greater orders of convergence for big enough value of $\theta$ and $\tilde{\theta}$ , as can be seen from Table 3. We measure the rate of convergence computing the value of $\beta$ as before. In fact the rate of convergence for $\theta=\tilde{\theta}=0.5$ or is close to the rate of convergence for smooth problems. In this case the exact eigenvalue $\lambda$ is unknown, but we approximate it by computing the eigenvalue on a very fine mesh involving about half a million of DOFs.

In Figure 3 we depict the mesh coming from the fourth iteration of Algorithm 1 with $\theta=\tilde{\theta}=0.8$ . This mesh is the result of multiple refinements using both marking strategies 1 and 2 each time. As can be seen the corners of the subdomain are much more refined than the rest of the mesh. This is clearly the effect of the first marking strategy, since the edge residuals have detected the discontinuity in the gradient of the eigenfunction along the interface.

Finally in Figure 4 we depict the eigenfunction corresponding to the smallest eigenvalue of the problem with discontinuous coefficients. This eigenfunction is the one used to refine the mesh in Figure 3.

$ \theta=\tilde{\theta}=0.2$ — **Table 3:** Comparison of the reduction of the error and DOFs of the adaptive method for the smallest eigenvalue for the problem with discontinuous coefficients.

$ \theta=\tilde{\theta}=0.5$ — **Table 3:** Comparison of the reduction of the error and DOFs of the adaptive method for the smallest eigenvalue for the problem with discontinuous coefficients.

	$\theta=\tilde{\theta}=0.2$	$\theta=\tilde{\theta}=0.5$	$\theta=\tilde{\theta}=0.8$
Iteration	$\vert\lambda-\lambda_n\vert$	DOFs	$\beta$	$\vert\lambda-\lambda_n\vert$	DOFs	$\beta$	$\vert\lambda-\lambda_n\vert$	DOFs	$\beta$
1	1.1071	81	–	1.1071	81	–	1.1071	81	–
2	1.0200	103	0.3410	0.8738	199	0.2632	0.4834	356	0.5597
3	1.0105	129	0.0416	0.5848	314	0.8805	0.2244	799	0.9494
4	1.0039	147	0.0498	0.3983	491	0.8591	0.0990	2235	0.7957
5	0.8968	167	0.8843	0.2766	673	1.1564	0.0401	4764	1.1932
6	0.8076	194	0.6996	0.1933	975	0.9665	0.0180	12375	0.8372
7	0.8008	217	0.0747	0.1346	1476	0.8722	0.0065	12375	1.1888
8	0.7502	237	0.7401	0.0948	2080	1.0237	0.0020	65387	1.4482