65N50 Mesh generation and refinement
Refine
Language
- English (7)
Keywords
- AFEM (2)
- eigenvalue problem (2)
- AMFEM (1)
- Adaptive finite element methods (1)
- Krylov subspace method (1)
- adaptive finite element method (1)
- adaptive finite element method (AFEM) (1)
- adaptive mesh refinement (1)
- adaptive mixed finite element method (1)
- compressed sensing (1)
- convergence (1)
- convergence analysis (1)
- dictionary (1)
- eigenvalue problems (1)
- elliptic eigenvalue problem (1)
- elliptic obstacle problems (1)
- error reduction (1)
- finite element method (FEM) (1)
- goal oriented error estimation (1)
- hierarchical basis (1)
- homotopy (1)
- linear programming (1)
- mutual incoherence (1)
- optimal control (1)
- optimal convergence (1)
- partial differential equation (1)
- restricted isometry property (1)
- sparse solution (1)
This paper presents three different adaptive algorithms for eigenvalue problems
associated with non-selfadjoint partial differential
operators. The basis for the developed algorithms is a homotopy method.
The homotopy method starts from a well-understood selfadjoint problem,
for which well-established adaptive methods are available.
Apart from the adaptive grid refinement, the progress of the homotopy as
well as the solution of the iterative method are adapted to balance the contributions
of the different error sources.
The first algorithm balances the homotopy, discretization and approximation errors with respect
to a fixed step-size $\tau$ in the homotopy.
The second algorithm combines the adaptive step-size control for the homotopy
with an adaptation in space that ensures an error below a fixed tolerance $\varepsilon$.
The third algorithm allows the complete adaptivity in space,
homotopy step-size as well as the iterative algebraic eigenvalue solver.
All three algorithms are compared in numerical examples.
The paper proposes goal-oriented error estimation and mesh refinement
for optimal control problems with elliptic PDE constraints using the value
of the reduced cost functional as quantity of interest. Error representation,
hierarchical error estimators, and greedy-style error indicators are derived and
compared to their counterparts when using the all-at-once cost functional as
quantity of interest. Finally, the efficiency of the error estimator and generated
meshes are demonstrated on numerical examples.
Various applications in fluid dynamics and computational continuum mechanics motivate the development of reliable and efficient adaptive algorithms for mixed finite element methods. In order to save degrees of freedom, not all but just some selected set of finite element domains are refined. Hence the fundamental question of convergence as well as the question of optimality require new mathematical arguments. The presented adaptive algorithm for Raviart-Thomas mixed finite element methods solves the Poisson model problem, with optimal convergence rate.
Chen, Holst, and Xu presented "convergence and optimality of adaptive mixed finite element methods" (2008) following arguments of Rob Stevenson for the conforming finite element method. Their algorithm reduces oscillations separately, before approximating the solution by some adaptive algorithm in the spirit of W. Dörfler (1996). The algorithm proposed here appears more natural in switching to either reduction of the edge-error estimator or of the oscillations.
We consider a new adaptive finite element (AFEM) algorithm for elliptic PDE-eigenvalue problems.
In contrast to other approaches we incorporate the iterative solution of the resulting finite dimensional
algebraic eigenvalue problems into the adaptation process.
In this way we can balance the costs of the adaption process
for the mesh with the costs for the iterative eigenvalue method. We present error estimates that incorporate
the discretization errors, approximation errors in the eigenvalue solver and roundoff errors and use
these for the adaptation process. We show that for the adaptation process it is possible to restrict to
very few iterations
of a Krylov subspace solver for the eigenvalue problem on coarse meshes.
We present several examples and show that this new approach achieves
much better complexity than previous AFEM approaches which assume that the algebraic
eigenvalue problem is solved to full accuracy.
We consider an adaptive finite element method (AFEM) for obstacle problems associated with linear second order elliptic boundary value problems and prove a reduction in the energy norm of the discretization error which leads to $R$-linear convergence. This result is shown to hold up to a consistery error due to the extension of the discrete multipliers (point functionals) to $H^{-1}$ and a possible mismatch between the continuous and discrete coincidence and noncoincidence sets. The AFEM is based on a residual-type error estimator consisting of element and edge residuals. The a posteriori error analysis reveals that the significant difference to the unconstrained case lies in the fact that these residuals only have to be taken into account within the discrete noncoincidence set. The proof of the error reduction property uses the reliability and the discrete local efficiency of the estimator as well as a perturbed Galerkin orthogonality. Numerical results are given illustrating the performance of the AFEM.
A new concept is introduced for the adaptive finite element discretization of partial differential equations that have a sparsely
representable solution. Motivated by recent work on compressed sensing, a recursive mesh refinement procedure is presented that uses linear programming to find a good approximation to the sparse solution on a given refinement level. Then only those parts of the mesh are refined that belong to nonzero expansion coefficients. Error estimates for this procedure are refined and the behavior of the procedure is demonstrated via some simple elliptic model problems.
A refined a posteriori error analysis for symmetric eigenvalue problems and the convergence of the first-order adaptive finite element method (AFEM) is presented. The $H^1$ stability of the $L^2$ projection provides reliability and efficiency of the edge-contribution of standard residual-based error estimators for $P_1$ finite element methods. In fact, the volume contributions and even oscillations can be omitted for Courant finite element methods. This allows for a refined averaging scheme and so improves [Dong Mao, Lihua Shen and Aihui Zhou, Adaptive finite element algorithms for eigenvalue problems based on local averaging type a posteriori error estimates, Advanced in Computational Mathematics, 2006, 25: 135-160]. The proposed AFEM
monitors the edge-contributions in a bulk criterion and so enables a contraction property up to higher-order terms and global convergence. Numerical experiments exploit the remaining $L^2$ error contributions and confirm our theoretical findings. The averaging schemes show a high accuracy and the AFEM leads to optimal empirical convergence rates.