Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAspects of an adaptive finite element method for the fractional Laplacian: a priori and a posteriori error estimates, efficient implementation and multigrid solver
132
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We develop all of the components needed to construct an adaptive finite element code that can be used to approximate fractional partial differential equations, on non-trivial domains in $dgeq 1$ dimensions. Our main approach consists of taking tools that have been shown to be effective for adaptive boundary element methods and, where necessary, modifying them so that they can be applied to the fractional PDE case. Improved a priori error estimates are derived for the case of quasi-uniform meshes which are seen to deliver sub-optimal rates of convergence owing to the presence of singularities. Attention is then turned to the development of an a posteriori error estimate and error indicators which are suitable for driving an adaptive refinement procedure. We assume that the resulting refined meshes are locally quasi-uniform and develop efficient methods for the assembly of the resulting linear algebraic systems and their solution using iterative methods, including the multigrid method. The storage of the dense matrices along with efficient techniques for computing the dense matrix vector products needed for the iterative solution is also considered. The performance and efficiency of the resulting algorithm is illustrated for a variety of examples.
rate research
Read More
A numerical scheme is presented for approximating fractional order Poisson problems in two and three dimensions. The scheme is based on reformulating the original problem posed over $Omega$ on the extruded domain $mathcal{C}=Omegatimes[0,infty)$ following Caffarelli and Silvestre (2007). The resulting degenerate elliptic integer order PDE is then approximated using a hybrid FEM-spectral scheme. Finite elements are used in the direction parallel to the problem domain $Omega$, and an appropriate spectral method is used in the extruded direction. The spectral part of the scheme requires that we approximate the true eigenvalues of the integer order Laplacian over $Omega$. We derive an a priori error estimate which takes account of the error arising from using an approximation in place of the true eigenvalues. We further present a strategy for choosing approximations of the eigenvalues based on Weyls law and finite element discretizations of the eigenvalue problem. The system of linear algebraic equations arising from the hybrid FEM-spectral scheme is decomposed into blocks which can be solved effectively using standard iterative solvers such as multigrid and conjugate gradient. Numerical examples in two and three dimensions show that the approach is quasi-optimal in terms of complexity.
In this work we study a residual based a posteriori error estimation for the CutFEM method applied to an elliptic model problem. We consider the problem with non-polygonal boundary and the analysis takes into account the geometry and data approximation on the boundary. The reliability and efficiency are theoretically proved. Moreover, constants are robust with respect to how the domain boundary cuts the mesh.
For the Stokes equation over 2D and 3D domains, explicit a posteriori and a priori error estimation are novelly developed for the finite element solution. The difficulty in handling the divergence-free condition of the Stokes equation is solved by utilizing the extended hypercircle method along with the Scott-Vogelius finite element scheme. Since all terms in the error estimation have explicit values, by further applying the interval arithmetic and verified computing algorithms, the computed results provide rigorous estimation for the approximation error. As an application of the proposed error estimation, the eigenvalue problem of the Stokes operator is considered and rigorous bounds for the eigenvalues are obtained. The efficiency of proposed error estimation is demonstrated by solving the Stokes equation on both convex and non-convex 3D domains.
We explore the connection between fractional order partial differential equations in two or more spatial dimensions with boundary integral operators to develop techniques that enable one to efficiently tackle the integral fractional Laplacian. In particular, we develop techniques for the treatment of the dense stiffness matrix including the computation of the entries, the efficient assembly and storage of a sparse approximation and the efficient solution of the resulting equations. The main idea consists of generalising proven techniques for the treatment of boundary integral equations to general fractional orders. Importantly, the approximation does not make any strong assumptions on the shape of the underlying domain and does not rely on any special structure of the matrix that could be exploited by fast transforms. We demonstrate the flexibility and performance of this approach in a couple of two-dimensional numerical examples.
Finite element exterior calculus (FEEC) has been developed over the past decade as a framework for constructing and analyzing stable and accurate numerical methods for partial differential equations by employing differential complexes. The recent work of Arnold, Falk and Winther cite{ArFaWi2010} includes a well-developed theory of finite element methods for Hodge Laplace problems, including a priori error estimates. In this work we focus on developing a posteriori error estimates in which the computational error is bounded by some computable functional of the discrete solution and problem data. More precisely, we prove a posteriori error estimates of residual type for Arnold-Falk-Winther mixed finite element methods for Hodge-de Rham Laplace problems. While a number of previous works consider a posteriori error estimation for Maxwells equations and mixed formulations of the scalar Laplacian, the approach we take is distinguished by unified treatment of the various Hodge Laplace problems arising in the de Rham complex, consistent use of the language and analytical framework of differential forms, and the development of a posteriori error estimates for harmonic forms and the effects of their approximation on the resulting numerical method for the Hodge Laplacian.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
