Abstract:In this paper, we present a fast divergence-free spectral algorithm (FDSA) for the curl-curl problem. Divergence-free bases in two and three dimensions are constructed by using the generalized Jacobi polynomials. An accurate spectral method with exact preservation of the divergence-free constraint point-wisely is then proposed, and its corresponding error estimate is established. We then present a highly efficient solution algorithm based on a combination of matrix-free preconditioned Krylov subspace iterative method and a fully diagonalizable auxiliary problem, which is derived from the spectral discretisations of generalized eigenvalue problems of Laplace and biharmonic operators. We rigorously prove that the dimensions of the invariant subspace of the preconditioned linear system resulting from the divergence-free spectral method with respect to the dominate eigenvalue $1$, are $(N-3)^2$ and $2(N-3)^3$ for two- and three-dimensional problems with $(N-1)^2$ and $2(N-1)^3$ unknowns, respectively. Thus, the proposed method usually takes only several iterations to converge, and astonishingly, as the problem size (polynomial order) increases, the number of iterations will decrease, even for highly indefinite system and oscillatory solutions. As a result, the computational cost of the solution algorithm is only a small multiple of $N^3$ and $N^4$ floating number operations for 2D and 3D problems, respectively. Plenty of numerical examples for solving the curl-curl problem with both constant and variable coefficients in two and three dimensions are presented to demonstrate the accuracy and efficiency of the proposed method.

Gradient recovery based finite element methods for the two-dimensional quad-curl problem

A Finite Element Method by Patch Reconstruction for the Quad-Curl Problem Using Mixed Formulations

Recovery based finite element method for biharmonic equation in two dimensional

Unfitted finite element method for the quad-curl interface problem

Finite Element Solution in Fluid Mechanics Using Exponential Function Based Interpolation

Distributional Finite Element curl div Complexes and Application to Quad Curl Problems

An $H^2$(curl)-conforming finite element in 2D and its applications to the quad-curl problem

H($\text{curl}^2$)-conforming quadrilateral spectral element method for quad-curl problems

Superconvergent recovery of edge finite element approximation for Maxwell’s equations

A new recovery based C 0 element method for fourth-order equations

Robust mixed finite element methods for a quad-curl singular perturbation problem

Recovery Techniques for Finite Element Methods

Gradient recovery-based adaptive stabilized mixed FEM for the convection–diffusion–reaction equation on surfaces

A highly efficient and accurate divergence-free spectral method for curl-curl equation in two and three dimensions

A least-squares Galerkin approach to gradient recovery for Hamilton-Jacobi-Bellman equation with Cordes coefficients

Gradient recovery for elliptic interface problem: I. body-fitted mesh

A recovery-based linear C0 finite element method for a fourth-order singularly perturbed Monge-Ampère equation

Simple curl-curl-conforming finite elements in two dimensions

Cubic superconvergent finite volume element method for one-dimensional elliptic and parabolic equations

Continuous Linear Finite Element Method for Biharmonic Problems on Surfaces

Biquadratic element discrete duality finite volume method for solving elliptic equations on quadrilateral mesh