New block triangular preconditioners for saddle point linear systems with highly singular blocks.
On a multilevel Krylov method for the Helmholtz equation preconditioned by shifted Laplacian.
On adaptive BDDC for the flow in heterogeneous porous media
We study a method based on Balancing Domain Decomposition by Constraints (BDDC) for numerical solution of a single-phase flow in heterogeneous porous media. The method solves for both flux and pressure variables. The fluxes are resolved in three steps: the coarse solve is followed by subdomain solves and last we look for a divergence-free flux correction and pressures using conjugate gradients with the BDDC preconditioner. Our main contribution is an application of the adaptive algorithm for selection...
On finite element approximation of fluid structure interaction by Taylor-Hood and Scott-Vogelius elements
This paper focuses on mathematical modeling and finite element simulation of fluid-structure interaction problems. A simplified problem of two-dimensional incompressible fluid flow interacting with a rigid structure, whose motion is described with one degree of freedom, is considered. The problem is mathematically described and numerically approximated using the finite element method. Two possibilities, namely Taylor-Hood and Scott-Vogelius elements are presented and implemented. Finally, numerical...
On the equivalence of primal and dual substructuring preconditioners.
On the parameter in augmented Lagrangian preconditioning for isogeometric discretizations of the Navier-Stokes equations
In this paper, we deal with the optimal choice of the parameter for augmented Lagrangian preconditioning of GMRES method for efficient solution of linear systems obtained from discretization of the incompressible Navier-Stokes equations. We consider discretization of the equations using the B-spline based isogeometric analysis approach. We are interested in the dependence of the convergence on the parameter for various problem parameters (Reynolds number, mesh refinement) and especially for...
On the relation between the AINV and the FAPINV algorithms.
Operator preconditioning with efficient applications for nonlinear elliptic problems
This paper is devoted to the numerical solution of nonlinear elliptic partial differential equations. Such problems describe various phenomena in science. An approach that exploits Hilbert space theory in the numerical study of elliptic PDEs is the idea of preconditioning operators. In this survey paper we briefly summarize the main lines of this theory with various applications.
-regular splitting iterative methods for non-Hermitian positive definite linear systems.
Parallel solution of elasticity problems using overlapping aggregations
The finite element (FE) solution of geotechnical elasticity problems leads to the solution of a large system of linear equations. For solving the system, we use the preconditioned conjugate gradient (PCG) method with two-level additive Schwarz preconditioner. The preconditioning is realised in parallel. A coarse space is usually constructed using an aggregation technique. If the finite element spaces for coarse and fine problems on structural grids are fully compatible, relations between elements...
Preconditioning of two-by-two block matrix systems with square matrix blocks, with applications
Two-by-two block matrices of special form with square matrix blocks arise in important applications, such as in optimal control of partial differential equations and in high order time integration methods. Two solution methods involving very efficient preconditioned matrices, one based on a Schur complement reduction of the given system and one based on a transformation matrix with a perturbation of one of the given matrix blocks are presented. The first method involves an additional inner solution...
Robust domain decomposition preconditioners for abstract symmetric positive definite bilinear forms
An abstract framework for constructing stable decompositions of the spaces corresponding to general symmetric positive definite problems into “local” subspaces and a global “coarse” space is developed. Particular applications of this abstract framework include practically important problems in porous media applications such as: the scalar elliptic (pressure) equation and the stream function formulation of its mixed form, Stokes’ and Brinkman’s equations. The constant in the corresponding abstract...
Robust domain decomposition preconditioners for abstract symmetric positive definite bilinear forms
An abstract framework for constructing stable decompositions of the spaces corresponding to general symmetric positive definite problems into “local” subspaces and a global “coarse” space is developed. Particular applications of this abstract framework include practically important problems in porous media applications such as: the scalar elliptic (pressure) equation and the stream function formulation of its mixed form, Stokes’ and Brinkman’s equations....
Robust operator estimates and the application to substructuring methods for first-order systems
We discuss a family of discontinuous Petrov–Galerkin (DPG) schemes for quite general partial differential operators. The starting point of our analysis is the DPG method introduced by [Demkowicz et al., SIAM J. Numer. Anal. 49 (2011) 1788–1809; Zitelli et al., J. Comput. Phys. 230 (2011) 2406–2432]. This discretization results in a sparse positive definite linear algebraic system which can be obtained from a saddle point problem by an element-wise Schur complement reduction applied to the test space....
Robust preconditioners for the matrix free truncated Newton method
New positive definite preconditioners for the matrix free truncated Newton method are given. Corresponding algorithms are described in detail. Results of numerical experiments that confirm the efficiency and robustness of the preconditioned truncated Newton method are reported.
Robust semi-coarsening multilevel preconditioning of biquadratic FEM systems
While a large amount of papers are dealing with robust multilevel methods and algorithms for linear FEM elliptic systems, the related higher order FEM problems are much less studied. Moreover, we know that the standard hierarchical basis two-level splittings deteriorate for strongly anisotropic problems. A first robust multilevel preconditioner for higher order FEM systems obtained after discretizations of elliptic problems with an anisotropic diffusion tensor is presented in this paper. We study...
Spectral discretization of Darcy equations coupled with Navier-Stokes equations by vorticity-velocity-pressure formulation
We consider a model coupling the Darcy equations in a porous medium with the Navier-Stokes equations in the cracks, for which the coupling is provided by the pressure's continuity on the interface. We discretize the coupled problem by the spectral element method combined with a nonoverlapping domain decomposition method. We prove the existence of solution for the discrete problem and establish an error estimation. We conclude with some numerical tests confirming the results of our analysis.
Sweeping preconditioners for elastic wave propagation with spectral element methods
We present a parallel preconditioning method for the iterative solution of the time-harmonic elastic wave equation which makes use of higher-order spectral elements to reduce pollution error. In particular, the method leverages perfectly matched layer boundary conditions to efficiently approximate the Schur complement matrices of a block LDLT factorization. Both sequential and parallel versions of the algorithm are discussed and results for large-scale problems from exploration geophysics are presented....
Total overlapping Schwarz' preconditioners for elliptic problems
A variant of the Total Overlapping Schwarz (TOS) method has been introduced in [Ben Belgacem et al., C. R. Acad. Sci., Sér. 1 Math. 336 (2003) 277–282] as an iterative algorithm to approximate the absorbing boundary condition, in unbounded domains. That same method turns to be an efficient tool to make numerical zooms in regions of a particular interest. The TOS method enjoys, then, the ability to compute small structures one wants to capture and the reliability to obtain the behavior of the solution...
Total overlapping Schwarz' preconditioners for elliptic problems
A variant of the Total Overlapping Schwarz (TOS) method has been introduced in [Ben Belgacem et al., C. R. Acad. Sci., Sér. 1 Math.336 (2003) 277–282] as an iterative algorithm to approximate the absorbing boundary condition, in unbounded domains. That same method turns to be an efficient tool to make numerical zooms in regions of a particular interest. The TOS method enjoys, then, the ability to compute small structures one wants to capture and the reliability to obtain the...