Auxiliary iterative schemes for the discrete operators on de Rham complex

05/05/2021
by   Zhongjie Lu, et al.
0

The main difficulty in solving the discrete source or eigenvalue problems of the operator d^*d with iterative methods is to deal with its huge kernel, for example, the ∇×∇× and - ∇ ( ∇· ) operator. In this paper, we construct a kind of auxiliary schemes for their discrete systems based on Hodge Laplacian on de Rahm complex. The spectra of the new schemes are Laplace-like. Then many efficient iterative methods and preconditioning techniques can be applied to them. After getting the solutions of the auxiliary schemes, the desired solutions of the original systems can be recovered or recognized through some simple operations. We sum these up as a new framework to compute the discrete source and eigenvalue problems of the operator d^*d using iterative method. We also investigate two preconditioners for the auxiliary schemes, ILU-type method and Multigrid method. Finally, we present plenty of numerical experiments to verify the efficiency of the auxiliary schemes.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
04/10/2021

Convergence analysis of numerical schemes for the Darcy-Forchheimer problem

This paper deals with the Darcy-Forchheimer problem with two kinds of bo...
research
07/14/2021

Solving discrete constrained problems on de Rham complex

The main difficulty in solving the discrete constrained problem is its p...
research
07/16/2021

Nodal auxiliary space preconditioning for the surface de Rham complex

This work develops optimal preconditioners for the discrete H(curl) and ...
research
11/27/2020

Fast auxiliary space preconditioners on surfaces

This work presents new uniform preconditioners for the discrete Laplace-...
research
05/12/2022

An Efficient Operator-Splitting Method for the Eigenvalue Problem of the Monge-Ampère Equation

We develop an efficient operator-splitting method for the eigenvalue pro...
research
01/04/2022

On an eigenvalue property of Summation-By-Parts operators

Summation-By-Parts (SBP) methods provide a systematic way of constructin...
research
10/06/2020

Iterative Methods for Computing Eigenvectors of Nonlinear Operators

In this chapter we are examining several iterative methods for solving n...

Please sign up or login with your details

Forgot password? Click here to reset