A matrix-free parallel solution method for the three-dimensional heterogeneous Helmholtz equation

The Helmholtz equation is related to seismic exploration, sonar, antennas, and medical imaging applications. It is one of the most challenging problems to solve in terms of accuracy and convergence due to the scalability issues of the numerical solvers. For 3D large-scale applications, high-performance parallel solvers are also needed. In this...

Combining p-multigrid and Multigrid Reduction in Time methods to obtain a scalable solver for Isogeometric Analysis

The use of sequential time integration schemes becomes more and more the bottleneck within large-scale computations due to a stagnation of processor’s clock speeds. In this study, we combine the parallel-in-time Multigrid Reduction in Time method with a p-multigrid method to obtain a scalable solver specifically designed for Isogeometric...

Multigrid Reduced in Time for Isogeometric Analysis

Isogeometric Analysis (IgA) can be seen as the natural extension of the Finite Element Method (FEM) to high-order B-spline basis functions. Combined with a time inte- gration scheme within the method of lines, IgA has become a viable alternative to FEM for time-dependent problems. However, as processors' clock speeds are no longer increasing but...

p-multigrid methods and their comparison to h-multigrid methods within Isogeometric Analysis

Over the years, Isogeometric Analysis has shown to be a successful alternative to the Finite Element Method (FEM). However, solving the resulting linear systems of equations efficiently remains a challenging task. In this paper, we consider a p-multigrid method, in which coarsening is applied in the spline degree p instead of the mesh width h...

Efficient multigrid based solvers for B-spline MPM

The Material Point Method (MPM) has been applied successfully to problems in engineering which involve large deformations and history-dependent material behavior. However, the classical method suffers from some shortcomings which influence the quality of the numerical solution significantly. High-order B-spline basis functions solve the problem...

Efficient Multigrid based solvers for isogeometric analysis

Introduced in [1], Isogeometric Analysis (IgA) has become widely accepted in academia and industry. However, solving the resulting linear systems remains a challenging task. For instance, the condition number of the Poisson operator scales quadratically with the mesh width h, but, in contrast to standard Finite Elements, exponentially with the...

Evaluation of multilevel sequentially semiseparable preconditioners on computational fluid dynamics benchmark problems using Incompressible Flow and Iterative Solver Software

This paper studies a new preconditioning technique for sparse systems arising from discretized partial differential equations in computational fluid dynamics problems. This preconditioning technique exploits the multilevel sequentially semiseparable (MSSS) structure of the system matrix. MSSS matrix computations give a data-sparse way to...

Reduction of computing time for least-squares migration based on the Helmholtz equation by graphics processing units

In geophysical applications, the interest in least-squares migration (LSM) as an imaging algorithm is increasing due to the demand for more accurate solutions and the development of high-performance computing. The computational engine of LSM in this work is the numerical solution of the 3D Helmholtz equation in the frequency domain. The...

Reduction of computing time for least-squares migration based on the Helmholtz equation by graphics processing units

In geophysical applications, the interest in least-squares migration (LSM) as an imaging algorithm is increasing due to the demand for more accurate solutions and the development of high-performance computing. The computational engine of LSM in this work is the numerical solution of the 3D Helmholtz equation in the frequency domain. The...

Evaluation of Multilevel Sequentially Semiseparable Preconditioners on CFD Benchmark Problems Using IFISS

This paper studies a new preconditioning technique for sparse systems arising from discretized partial differential equations (PDEs) in computational fluid dynamics (CFD), which exploit the multilevel sequentially semiseparable (MSSS) structure of the system matrix. MSSS matrix computations give a data-sparse way to approximate the LU...

Comparison of the deflated preconditioned conjugate gradient method and algebraic multigrid for composite materials

Many applications in computational science and engineering concern composite materials, which are characterized by large discontinuities in the material properties. Such applications require fine-scale finite-element meshes, which lead to large linear systems that are challenging to solve with current direct and iterative solutions algorithms....

Towards faster solution of large power flow problems

Current and future developments in the power system industry demand fast power flow solvers for larger power flow problems. The established methods are no longer viable for such problems, as they are not scalable in the problem size. In this paper, the use of Newton-Krylov power flow methods is proposed, and a multitude of preconditioning...

A Comparison of Two-Level Preconditioners Based on Multigrid and Deflation

It is well known that two-level and multilevel preconditioned conjugate gradient (PCG) methods provide efficient techniques for solving large and sparse linear systems whose coefficient matrices are symmetric and positive definite. A two-level PCG method combines a traditional (one-level) preconditioner, such as incomplete Cholesky, with a...

Comparison of Two-Level Preconditioners Derived from Deflation, Domain Decomposition and Multigrid Methods

For various applications, it is well-known that a multi-level, in particular two-level, preconditioned CG (PCG) method is an efficient method for solving large and sparse linear systems with a coefficient matrix that is symmetric positive definite. The corresponding two-level preconditioner combines traditional and projection-type...

A novel multigrid based preconditioner for heterogeneous Helmholtz problems

An iterative solution method, in the form of a preconditioner for a Krylov subspace method, is presented for the Helmholtz equation. The preconditioner is based on a Helmholtz type differential operator with a complex term. A multigrid iteration is used for approximately inverting the preconditioner. The choice of multigrid components for the...