Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHigh accuracy methods for eigenvalues of elliptic operators by nonconforming elements
93
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this paper, three high-accuracy methods for eigenvalues of second order elliptic operators are proposed by using the nonconforming Crouzeix-Raviart(CR for short) element and the nonconforming enriched Crouzeix-Raviart(ECR for short) element. They are based on a crucial full one order superconvergence of the first order mixed Raviart-Thomas(RT for short) element. The main ingredient of such a superconvergence analysis is to employ a discrete Helmholtz decomposition of the difference between the canonical interpolation and the finite element solution of the RT element. In particular, it allows for some vital cancellation between terms in one key sum of boundary terms. Consequently, a full one order superconvergence follows from a special relation between the CR element and the RT element, and the equivalence between the ECR element and the RT element for these two nonconforming elements. These superconvergence results improve those in literature from a half order to a full one order for the RT element, the CR element and the ECR element. Based on the aforementioned superconvergence of the RT element, asymptotic expansions of eigenvalues are established and employed to achieve high accuracy extrapolation methods for these two nonconforming elements. In contrast to a classic analysis in literature, the novelty herein is to use not only the canonical interpolations of these nonconforming elements but also that of the RT element to analyze such asymptotic expansions. Based on the superconvergence of these nonconforming elements, asymptotically exact a posteriori error estimators of eigenvalues are constructed and analyzed for them. Finally, two post-processing methods are proposed to improve accuracy of approximate eigenvalues by employing these a posteriori error estimators.Numerical tests are provided to justify and compare the performance of the aforementioned methods.
rate research
Read More
In this paper, an important discovery has been found for nonconforming immersed finite element (IFE) methods using integral-value degrees of freedom for solving elliptic interface problems. We show that those IFE methods can only achieve suboptimal convergence rates (i.e., $O(h^{1/2})$ in the $H^1$ norm and $O(h)$ in the $L^2$ norm) if the tangential derivative of the exact solution and the jump of the coefficient are not zero on the interface. A nontrivial counter example is also provided to support our theoretical analysis. To recover the optimal convergence rates, we develop a new nonconforming IFE method with additional terms locally on interface edges. The unisolvence of IFE basis functions is proved on arbitrary triangles. Furthermore, we derive the optimal approximation capabilities of both the Crouzeix-Raviart and the rotated-$Q_1$ IFE spaces for interface problems with variable coefficients via a unified approach different from multipoint Taylor expansions. Finally, optimal error estimates in both $H^1$- and $L^2$- norms are proved and confirmed with numerical experiments.
We propose two nonconforming finite elements to approximate a boundary value problem arising from strain gradient elasticity, which is a higher-order perturbation of the linearized elastic system. Our elements are H$^2-$nonconforming while H$^1-$conforming. We show both elements converges in the energy norm uniformly with respect to the perturbation parameter.
This paper is devoted to the computation of transmission eigenvalues in the inverse acoustic scattering theory. This problem is first reformulated as a two by two boundary system of boundary integral equations. Next, utilizing the Schur complement technique, we develop a Schur complement operator with regularization to obtain a reduced system of boundary integral equations. The Nystr{o}m discretization is then used to obtain an eigenvalue problem for a matrix. We employ the recursive integral method for the numerical computation of the matrix eigenvalue. Numerical results show that the proposed method is efficient and reduces computational costs.
We examine some numerical iterative methods for computing the eigenvalues and eigenvectors of real matrices. The five methods examined here range from the simple power iteration method to the more complicated QR iteration method. The derivations, procedure, and advantages of each method are briefly discussed.
We establish a new H2 Korns inequality and its discrete analog, which greatly simplify the construction of nonconforming elements for a linear strain gradient elastic model. The Specht triangle [41] and the NZT tetrahedron [45] are analyzed as two typical representatives for robust nonconforming elements in the sense that the rate of convergence is independent of the small material parameter. We construct new regularized interpolation estimate and the enriching operator for both elements, and prove the error estimates under minimal smoothness assumption on the solution. Numerical results are consistent with the theoretical prediction.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
