BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Hehu Xie

Hehu Xie appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath.OC

Catalog footprint

What is connected

37works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2024arXiv

Augmented Subspace Scheme for Eigenvalue Problem by Weak Galerkin Finite Element Method

This study proposes a class of augmented subspace schemes for the weak Galerkin (WG) finite element method used to solve eigenvalue problems. The augmented subspace is built with the conforming linear finite element space defined on the coarse mesh and the eigenfunction approximations in the WG finite element space defined on the fine mesh. Based on this augmented subspace, solving the eigenvalue problem in the fine WG finite element space can be reduced to the solution of the linear boundary value problem in the same WG finite element space and a low dimensional eigenvalue problem in the augmented subspace. The proposed augmented subspace techniques have the second order convergence rate with respect to the coarse mesh size, as demonstrated by the accompanying error estimates. Finally, a few numerical examples are provided to validate the proposed numerical techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A New Subspace Iteration Algorithm for Solving Generalized Eigenvalue Problems

It is needed to solve generalized eigenvalue problems (GEP) in many applications, such as the numerical simulation of vibration analysis, quantum mechanics, electronic structure, etc. The subspace iteration is a kind of widely used algorithm to solve eigenvalue problems. To solve the generalized eigenvalue problem, one kind of subspace iteration method, Chebyshev-Davidson algorithm, is proposed recently. In Chebyshev-Davidson algorithm, the Chebyshev polynomial filter technique is incorporated in the subspace iteration. In this paper, based on Chebyshev-Davidson algorithm, a new subspace iteration algorithm is constructed. In the new algorithm, the Chebyshev filter and inexact Rayleigh quotient iteration techniques are combined together to enlarge the subspace in the iteration. Numerical results of a vibration analysis problem show that the number of iteration and computing time of the proposed algorithm is much less than that of the Chebyshev-Davidson algorithm and some typical GEP solution algorithms. Furthermore, the new algorithm is more stable and reliable than the Chebyshev-Davidson algorithm in the numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An Efficient Adaptive Finite Element Method for Eigenvalue Problems

The aim of this paper is to propose an efficient adaptive finite element method for eigenvalue problems based on the multilevel correction scheme and inverse power method. This method involves solving associated boundary value problems on each adaptive partitions and very low dimensional eigenvalue problems on some special meshes which are controlled by the proposed algorithm. Since we Hence the efficiency of solving eigenvalue problems can be improved to be similar to the adaptive finite element method for the associated boundary value problems. The convergence and optimal complexity is theoretically verified and numerically demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the convergence to local limit of nonlocal models with approximated interaction neighborhoods

Many nonlocal models have adopted Euclidean balls as the nonlocal interaction neighborhoods. When solving them numerically, it is sometimes convenient to adopt polygonal approximations of such balls. A crucial question is, to what extent such approximations affect the nonlocal operators and the corresponding solutions. While recent works have analyzed this issue for a fixed horizon parameter, the question remains open in the case of a small or vanishing horizon parameter, which happens often in many practical applications and has significant impact on the reliability and robustness of nonlocal modeling and simulations. In this work, we are interested in addressing this issue and establishing the convergence of the nonlocal solutions associated with polygonally approximated interaction neighborhoods to the local limit of the original nonlocal solutions. Our finding reveals that the new nonlocal solution does not converge to the correct local limit when the number of sides of polygons is uniformly bounded. On the other hand, if the number of sides tends to infinity, the desired convergence can be established. These results may be used to guide future computational studies of nonlocal models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On accelerating a multilevel correction adaptive finite element method for Kohn-Sham equation

Based on the numerical method proposed in [G. Hu, X. Xie, F. Xu, J. Comput. Phys., 355 (2018), 436-449.] for Kohn-Sham equation, further improvement on the efficiency is obtained in this paper by i). designing a numerical method with the strategy of separately handling the nonlinear Hartree potential and exchange-correlation potential, and ii).parallelizing the algorithm in an eigenpairwise approach. The feasibility of two approaches are analyzed in detail, and the new algorithm is described completely. Compared with previous results, a significant improvement of numerical efficiency can be observed from plenty of numerical experiments, which make the new method more suitable for the practical problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Parallel Augmented Subspace Method for Eigenvalue Problems

A type of parallel augmented subspace scheme for eigenvalue problems is proposed by using coarse space in the multigrid method. With the help of coarse space in multigrid method, solving the eigenvalue problem in the finest space is decomposed into solving the standard linear boundary value problems and very low dimensional eigenvalue problems. The computational efficiency can be improved since there is no direct eigenvalue solving in the finest space and the multigrid method can act as the solver for the deduced linear boundary value problems. Furthermore, for different eigenvalues, the corresponding boundary value problem and low dimensional eigenvalue problem can be solved in the parallel way since they are independent of each other and there exists no data exchanging. This property means that we do not need to do the orthogonalization in the highest dimensional spaces. This is the main aim of this paper since avoiding orthogonalization can improve the scalability of the proposed numerical method. Some numerical examples are provided to validate the proposed parallel augmented subspace method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Algebraic Multigrid Method for Eigenvalue Problems in Some Different Cases

The aim of this paper is to develop an algebraic multigrid method to solve eigenvalue problems based on the combination of the multilevel correction scheme and the algebraic multigrid method for linear equations. Our approach uses the algebraic multigrid method setup procedure to construct the hierarchy and the intergrid transfer operators. In this algebraic multigrid scheme, a large scale eigenvalue problem is solved by some algebraic multigrid smoothing steps in the hierarchy and very small-dimensional eigenvalue problems. To emphasize the efficiency and flexibility of the proposed method, here we consider a set of test eigenvalue problems, discretized on unstructured meshes, with different shape of domain, singularity, and discontinuous parameters. Moreover, global convergence independent of the number of desired eigenvalues is obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Multigrid Method for the Ground State Solution of Bose-Einstein Condensates Based on Newton Iteration

In this paper, a new kind of multigrid method is proposed for the ground state solution of Bose-Einstein condensates based on Newton iteration method. Instead of treating eigenvalue $λ$ and eigenvector $u$ respectively, we regard the eigenpair $(λ, u)$ as one element in the composite space $\R \times H_0^1(Ω)$ and then Newton iteration method is adopted for the nonlinear problem. Thus in this multigrid scheme, we only need to solve a linear discrete boundary value problem in every refined space, which can improve the overall efficiency for the simulation of Bose-Einstein condensations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Multilevel Correction Method for Interior Transmission Eigenvalue Problem

In this paper, we give a numerical analysis for the transmission eigenvalue problem by the finite element method. A type of multilevel correction method is proposed to solve the transmission eigenvalue problem. The multilevel correction method can transform the transmission eigenvalue solving in the finest finite element space to a sequence of linear problems and some transmission eigenvalue solving in a very low dimensional spaces. Since the main computational work is to solve the sequence of linear problems, the multilevel correction method improves the overfull efficiency of the transmission eigenvalue solving. Some numerical examples are provided to validate the theoretical results and the efficiency of the proposed numerical scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Multilevel Correction Scheme for Nonsymmetric Eigenvalue Problems by Finite Element Methods

A multilevel correction scheme is proposed to solve defective and nodefective of nonsymmetric partial differential operators by the finite element method. The method includes multi correction steps in a sequence of finite element spaces. In each correction step, we only need to solve two source problems on a finer finite element space and two eigenvalue problems on the coarsest finite element space. The accuracy of the eigenpair approximation is improved after each correction step. This correction scheme improves overall efficiency of the finite element method in solving nonsymmetric eigenvalue problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Anisotropic meshes and stabilized parameters for the stabilized finite element methods

We propose a numerical strategy to generate the anisotropic meshes and select the appropriate stabilized parameters simultaneously for two dimensional convection-dominated convection-diffusion equations by stabilized continuous linear finite elements. Since the discretized error in a suitable norm can be bounded by the sum of interpolation error and its variants in different norms, we replace them by some terms which contain the Hessian matrix of the true solution, convective fields, and the geometric properties such as directed edges and the area of the triangle. Based on this observation, the shape, size and equidistribution requirements are used to derive the corresponding metric tensor and the stabilized parameters. It is easily found from our derivation that the optimal stabilized parameter is coupled with the optimal metric tensor on each element. Some numerical results are also provided to validate the stability and efficiency of the proposed numerical strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computable Error Estimates for Ground State Solution of Bose-Einstein Condensates

In this paper, we propose a computable error estimate of the Gross-Pitaevskii equation for ground state solution of Bose-Einstein condensates by general conforming finite element methods on general meshes. Based on the proposed error estimate, asymptotic lower bounds of the smallest eigenvalue and ground state energy can be obtained. Several numerical examples are presented to validate our theoretical results in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fully Computable Error Bounds for Eigenvalue Problem

This paper is concerned with the computable error estimates for the eigenvalue problem which is solved by the general conforming finite element methods on the general meshes. Based on the computable error estimate, we can give an asymptotically lower bound of the general eigenvalues. Furthermore, we also give a guaranteed upper bound of the error estimates for the first eigenfunction approximation and a guaranteed lower bound of the first eigenvalue based on computable error estimator. Some numerical examples are presented to validate the theoretical results deduced in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Cascadic Multigrid Method for GPE Problem

A cascadic multigrid method is proposed for the GPE problem based on the multilevel correction scheme. With this new scheme, the ground state eigenvalue problem on the finest space can be solved by smoothing steps on a series of multilevel finite element spaces and some nonlinear eigenvalue problem solving on a very low-dimensional space. Choosing the appropriate sequence of finite element spaces and the number of smoothing steps, the optimal convergence rate with the optimal computational work can be arrived. Some numerical experiments are presented to validate our theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Full Multigrid Method for Nonlinear Eigenvalue Problems

This paper is to introduce a type of full multigrid method for the nonlinear eigenvalue problem. The main idea is to transform the solution of nonlinear eigenvalue problem into a series of solutions of the corresponding linear boundary value problems on the sequence of finite element spaces and nonlinear eigenvalue problems on the coarsest finite element space. The linearized boundary value problems are solved by some multigrid iterations. Besides the multigrid iteration, all other efficient iteration methods for solving boundary value problems can serve as the linear problem solver. We will prove that the computational work of this new scheme is truly optimal, the same as solving the linear corresponding boundary value problem. In this case, this type of iteration scheme certainly improves the overfull efficiency of solving nonlinear eigenvalue problems. Some numerical experiments are presented to validate the efficiency of the new method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Multigrid Method for Nonlinear Eigenvalue Problems

A multigrid method is proposed for solving nonlinear eigenvalue problems by the finite element method. With this new scheme, solving nonlinear eigenvalue problem is decomposed to a series of solutions of linear boundary value problems on multilevel finite element spaces and a series of small scale nonlinear eigenvalue problems. The computational work of this new scheme can reach almost the same as the solution of the corresponding linear boundary value problem. Therefore, this type of multilevel correction scheme improves the overfull efficiency of the nonlinear eigenvalue problem solving.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Multilevel Newton Iteration Method for Eigenvalue Problems

We propose a new type of multilevel method for solving eigenvalue problems based on Newton iteration. With the proposed iteration method, solving eigenvalue problem on the finest finite element space is replaced by solving a small scale eigenvalue problem in a coarse space and solving a series of augmented linear problems, derived by Newton step in the corresponding series of finite element spaces. This iteration scheme improves overall efficiency of the finite element method for solving eigenvalue problems. Finally, some numerical examples are provided to validate the efficiency of the proposed numerical scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sharp Error Bounds for Piecewise Polynomial Approximation: Revisit and Application to Elliptic PDE Eigenvalue Computation

In this paper, we revisit approximation properties of piecewise polynomial spaces, which contain more than ${\cal P}_{r-1}$ but not ${\cal P}_r$. We develop more accurate upper and lower error bounds that are sharper than those used in literature. These new error bounds, especially the lower bounds are particular useful to finite element methods. As an important application, we establish sharp lower bounds of the discretization error for Laplace and $2m$-th order elliptic eigenvalue problems in various finite element spaces under shape regular triangulations, and investigate the asymptotic convergence behavior for large numerical eigenvalue approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The weak Galerkin method for eigenvalue problems

This article is devoted to computing the eigenvalue of the Laplace eigenvalue problem by the weak Galerkin (WG) finite element method with emphasis on obtaining lower bounds. The WG method is on the use of weak functions and their weak derivatives defined as distributions. Weak functions and weak derivatives can be approximated by polynomials with various degrees. Different combination of polynomial spaces leads to different WG finite element methods, which makes WG methods highly flexible and efficient in practical computation. We establish the optimal-order error estimates for the WG finite element approximation for the eigenvalue problem. Comparing with the classical nonconforming finite element method which can just provide lower bound approximation by linear elements with only the second order convergence, the WG methods can naturally provide lower bound approximation with a high order convergence (larger than $2$). Some numerical results are also presented to demonstrate the efficiency of our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Cascadic Multigrid Method for Eigenvalue Problem

A cascadic multigrid method is proposed for eigenvalue problems based on the multilevel correction scheme. With this new scheme, an eigenvalue problem on the finest space can be solved by smoothing steps on a series of multilevel finite element spaces and eigenvalue problem solving on the coarsest finite element space. Choosing the appropriate sequence of finite element spaces and the number of smoothing steps, the optimal convergence rate with the optimal computational work can be arrived. Some numerical experiments are presented to validate our theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Full Multigrid Method for Eigenvalue Problems

In this paper, a full (nested) multigrid scheme is proposed to solve eigenvalue problems. The idea here is to use the multilevel correction method to transform the solution of eigenvalue problem to a series of solutions of the corresponding boundary value problems and eigenvalue problems defined on the coarsest finite element space. The boundary value problems which are define on a sequence of multilevel finite element space can be solved by some multigrid iteration steps. Besides the multigrid iteration, all other efficient iteration methods for solving boundary value problems can serve as linear problem solver. The computational work of this new scheme can reach optimal order the same as solving the corresponding source problem. Therefore, this type of iteration scheme improves the efficiency of eigenvalue problem solving.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Multigrid Method Based On Shifted-Inverse Power Technique for Eigenvalue Problems

A multigrid method is proposed in this paper to solve eigenvalue problems by the finite element method based on the shifted-inverse power iteration technique. With this scheme, solving eigenvalue problem is transformed to a series of nonsingular solutions of boundary value problems on multilevel meshes. Since replacing the difficult eigenvalue solving by the easier solution of boundary value problems, the multigrid way can improve the overall efficiency of the eigenvalue problem solving. Some numerical experiments are presented to validate the efficiency of this new method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Multigrid Method for the Ground State Solution of Bose-Einstein Condensates

A multigrid method is proposed to compute the ground state solution of Bose-Einstein condensations by the finite element method based on the multilevel correction for eigenvalue problems and the multigrid method for linear boundary value problems. In this scheme, obtaining the optimal approximation for the ground state solution of Bose-Einstein condensates includes a sequence of solutions of the linear boundary value problems by the multigrid method on the multilevel meshes and a series of solutions of nonlinear eigenvalue problems on the coarsest finite element space. The total computational work of this scheme can reach almost the optimal order as same as solving the corresponding linear boundary value problem. Therefore, this type of multigrid scheme can improve the overall efficiency for the simulation of Bose-Einstein condensations. Some numerical experiments are provided to validate the efficiency of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A multilevel correction method for optimal controls of elliptic equation

We propose in this paper a multilevel correction method to solve optimal control problems constrained by elliptic equations with the finite element method. In this scheme, solving optimization problem on the finest finite element space is transformed to a series of solutions of linear boundary value problems by the multigrid method on multilevel meshes and a series of solutions of optimization problems on the coarsest finite element space. Our proposed scheme, instead of solving a large scale optimization problem in the finest finite element space, solves only a series of linear boundary value problems and the optimization problems in a very low dimensional finite element space, and thus can improve the overall efficiency for the solution of optimal control problems governed by PDEs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local a Priori Estimate on the General Scale Subdomains

The local a priori estimate for the finite element approximation is essential for underlying the local and parallel technique. It is well known that the constant coefficients in the inequality is independent of the mesh size. But it is not so clear whether the constant depends on the scale of the local subdomains. The aim of this note is to derive a new local a priori estimate on the general scale domains. We also show that the dependence of the constant appearing in the local a priori estimate on the scale of the subdomains.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local and Parallel Finite Element Algorithm Based On Multilevel Discretization for Eigenvalue Problem

A local and parallel algorithm based on the multilevel discretization is proposed in this paper to solve the eigenvalue problem by the finite element method. With this new scheme, solving the eigenvalue problem in the finest grid is transferred to solutions of the eigenvalue problems on the coarsest mesh and a series of solutions of boundary value problems by using the local and parallel algorithm. The computational work in each processor can reach the optimal order. Therefore, this type of multilevel local and parallel method improves the overall efficiency of solving the eigenvalue problem. Some numerical experiments are presented to validate the efficiency of the new method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Multilevel Correction Type of Adaptive Finite Element Method for Eigenvalue Problems

A type of adaptive finite element method for the eigenvalue problems is proposed based on the multilevel correction scheme. In this method, adaptive finite element method to solve eigenvalue problems involves solving associated boundary value problems on the adaptive partitions and small scale eigenvalue problems on the coarsest partitions. Hence the efficiency of solving eigenvalue problems can be improved to be similar to the adaptive finite element method for the associated boundary value problems. The convergence and optimal complexity is theoretically verified and numerically demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Parallel Method for Population Balance Equations Based on the Method of Characteristics

In this paper, we present a parallel scheme to solve the population balance equations based on the method of characteristics and the finite element discretization. The application of the method of characteristics transform the higher dimensional population balance equation into a series of lower dimensional convection-diffusion-reaction equations which can be solved in a parallel way.Some numerical results are presented to show the accuracy and efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Metric tensors for the interpolation error and its gradient in $L^p$ norm

A uniform strategy to derive metric tensors in two spatial dimension for interpolation errors and their gradients in $L^p$ norm is presented. It generates anisotropic adaptive meshes as quasi-uniform ones in corresponding metric space, with the metric tensor being computed based on a posteriori error estimates in different norms. Numerical results show that the corresponding convergence rates are always optimal.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Multi-level Correction Scheme for Eigenvalue Problems

In this paper, a new type of multi-level correction scheme is proposed for solving eigenvalue problems by finite element method. With this new scheme, the accuracy of eigenpair approximations can be improved after each correction step which only needs to solve a source problem on finer finite element space and an eigenvalue problem on the coarsest finite element space. This correction scheme can improve the efficiency of solving eigenvalue problems by finite element method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Computing the lower and upper bounds of Laplace eigenvalue problem: by combining conforming and nonconforming finite element methods

This article is devoted to computing the lower and upper bounds of the Laplace eigenvalue problem. By using the special nonconforming finite elements, i.e., enriched Crouzeix-Raviart element and extension $Q_1^{\rm rot}$, we get the lower bound of the eigenvalue. Additionally, we also use conforming finite elements to do the postprocessing to get the upper bound of the eigenvalue. The postprocessing method need only to solve the corresponding source problems and a small eigenvalue problem if higher order postprocessing method is implemented. Thus, we can obtain the lower and upper bounds of the eigenvalues simultaneously by solving eigenvalue problem only once. Some numerical results are also presented to validate our theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lower Bounds of the Discretization for Piecewise Polynomials

Assume that $V_h$ is a space of piecewise polynomials of degree less than $r\geq 1$ on a family of quasi-uniform triangulation of size $h$. Then the following well-known upper bound holds for a sufficiently smooth function $u$ and $p\in [1, \infty]$ $$ \inf_{v_h\in V_h}\|u-v_h\|_{j,p,Ω,h} \le C h^{r-j} |u|_{r,p,Ω},\quad 0\le j\le r. $$ In this paper, we prove that, roughly speaking, if $u\not\in V_h$, the above estimate is sharp. Namely, $$ \inf_{v_h\in V_h}\|u-v_h\|_{j,p,Ω,h} \ge c h^{r-j},\quad 0\le j\le r, \ \ 1\leq p\leq \infty, $$ for some $c>0$. The above result is further extended to various situations including more general Sobolev space norms, general shape regular grids and many different types of finite element spaces. As an application, the sharpness of finite element approximation of elliptic problems and the corresponding eigenvalue problems is established.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

New metric tensors for anisotropic mesh generation

A new anisotropic mesh adaptation strategy for finite element solution of elliptic differential equations is presented. It generates anisotropic adaptive meshes as quasi-uniform ones in some metric space, with the metric tensor being computed based on a posteriori error estimates proposed in \cite{YinXie}. The new metric tensor explores more comprehensive information of anisotropy for the true solution than those existing ones. Numerical results show that this approach can be successfully applied to deal with poisson and steady convection-dominated problems. The superior accuracy and efficiency of the new metric tensor to others is illustrated on various numerical examples of complex two-dimensional simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Postprocessing and Higher Order Convergence of Stabilized Finite Element Discretizations of the Stokes Eigenvalue Problem

In this paper, the stabilized finite element method based on local projection is applied to discretize the Stokes eigenvalue problems and the corresponding convergence analysis is given. Furthermore, we also use a method to improve the convergence rate for the eigenpair approximations of the Stokes eigenvalue problem. It is based on a postprocessing strategy that contains solving an additional Stokes source problem on an augmented finite element space which can be constructed either by refining the mesh or by using the same mesh but increasing the order of mixed finite element space. Numerical examples are given to confirm the theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Asymptotic Expansions and Extrapolation of Approximate Eigenvalues for Second Order Elliptic Problems by Mixed Finite Element Methods

In this paper, we derive an asymptotic error expansion for the eigenvalue approximations by the lowest order Raviart-Thomas mixed finite element method for the general second order elliptic eigenvalue problems. Extrapolation based on such an expansion is applied to improve the accuracy of the eigenvalue approximations. Furthermore, we also prove the superclose property between the finite element projection with the finite element approximation of the eigenvalue problems by mixed finite element methods. In order to prove the full order of the eigenvalue extrapolation, we first propose "the auxiliary equation method". The result of this paper provides a general procedure to produce an asymptotic expansions for eigenvalue approximations by mixed finite elements.

0 citations0 reviewsNo code linkedOpen access