709 research outputs found
An odyssey into local refinement and multilevel preconditioning III: Implementation and numerical experiments
In this paper, we examine a number of additive and multiplicative multilevel iterative methods and preconditioners in the setting of two-dimensional local mesh refinement. While standard multilevel methods are effective for uniform refinement-based discretizations of elliptic equations, they tend to be less effective for algebraic systems, which arise from discretizations on locally refined meshes, losing their optimal behavior in both storage and computational complexity. Our primary focus here is on Bramble, Pasciak, and Xu (BPX)-style additive and multiplicative multilevel preconditioners, and on various stabilizations of the additive and multiplicative hierarchical basis (HB) method, and their use in the local mesh refinement setting. In parts I and II of this trilogy, it was shown that both BPX and wavelet stabilizations of HB have uniformly bounded condition numbers on several classes of locally refined two- and three-dimensional meshes based on fairly standard (and easily implementable) red and red-green mesh refinement algorithms. In this third part of the trilogy, we describe in detail the implementation of these types of algorithms, including detailed discussions of the data structures and traversal algorithms we employ for obtaining optimal storage and computational complexity in our implementations. We show how each of the algorithms can be implemented using standard data types, available in languages such as C and FORTRAN, so that the resulting algorithms have optimal (linear) storage requirements, and so that the resulting multilevel method or preconditioner can be applied with optimal (linear) computational costs. We have successfully used these data structure ideas for both MATLAB and C implementations using the FEtk, an open source finite element software package. We finish the paper with a sequence of numerical experiments illustrating the effectiveness of a number of BPX and stabilized HB variants for several examples requiring local refinement
A Two-Level Method for Mimetic Finite Difference Discretizations of Elliptic Problems
We propose and analyze a two-level method for mimetic finite difference
approximations of second order elliptic boundary value problems. We prove that
the two-level algorithm is uniformly convergent, i.e., the number of iterations
needed to achieve convergence is uniformly bounded independently of the
characteristic size of the underling partition. We also show that the resulting
scheme provides a uniform preconditioner with respect to the number of degrees
of freedom. Numerical results that validate the theory are also presented
Multigrid Preconditioning for a Space-Time Spectral-Element Discontinuous-Galerkin Solver
In this work we examine a multigrid preconditioning approach in the context of a high- order tensor-product discontinuous-Galerkin spectral-element solver. We couple multigrid ideas together with memory lean and efficient tensor-product preconditioned matrix-free smoothers. Block ILU(0)-preconditioned GMRES smoothers are employed on the coarsest spaces. The performance is evaluated on nonlinear problems arising from unsteady scale- resolving solutions of the Navier-Stokes equations: separated low-Mach unsteady ow over an airfoil from laminar to turbulent ow. A reduction in the number of ne space iterations is observed, which proves the efficiency of the approach in terms of preconditioning the linear systems, however this gain was not reflected in the CPU time. Finally, the preconditioner is successfully applied to problems characterized by stiff source terms such as the set of RANS equations, where the simple tensor product preconditioner fails. Theoretical justification about the findings is reported and future work is outlined
Multilevel preconditioning based on discrete symmetrization for convection-diffusion equations
AbstractThe subject of this paper is an additive multilevel preconditioning approach for convection-diffusion problems. Our particular interest is in the convergence behavior for convection-dominated problems which are discretized by the streamline diffusion method. The multilevel preconditioner is based on a transformation of the discrete problem which reduces the relative size of the skew-symmetric part of the operator. For the constant coefficient case, an analysis of the convergence properties of this multilevel preconditioner is given in terms of its dependence on the convection size. Moreover, the results of computational experiments for more general convection-diffusion problems are presented and our new preconditioner is compared to standard multilevel preconditioning
Matrix-equation-based strategies for convection-diffusion equations
We are interested in the numerical solution of nonsymmetric linear systems
arising from the discretization of convection-diffusion partial differential
equations with separable coefficients and dominant convection. Preconditioners
based on the matrix equation formulation of the problem are proposed, which
naturally approximate the original discretized problem. For certain types of
convection coefficients, we show that the explicit solution of the matrix
equation can effectively replace the linear system solution. Numerical
experiments with data stemming from two and three dimensional problems are
reported, illustrating the potential of the proposed methodology
Multilevel Markov Chain Monte Carlo Method for High-Contrast Single-Phase Flow Problems
In this paper we propose a general framework for the uncertainty
quantification of quantities of interest for high-contrast single-phase flow
problems. It is based on the generalized multiscale finite element method
(GMsFEM) and multilevel Monte Carlo (MLMC) methods. The former provides a
hierarchy of approximations of different resolution, whereas the latter gives
an efficient way to estimate quantities of interest using samples on different
levels. The number of basis functions in the online GMsFEM stage can be varied
to determine the solution resolution and the computational cost, and to
efficiently generate samples at different levels. In particular, it is cheap to
generate samples on coarse grids but with low resolution, and it is expensive
to generate samples on fine grids with high accuracy. By suitably choosing the
number of samples at different levels, one can leverage the expensive
computation in larger fine-grid spaces toward smaller coarse-grid spaces, while
retaining the accuracy of the final Monte Carlo estimate. Further, we describe
a multilevel Markov chain Monte Carlo method, which sequentially screens the
proposal with different levels of approximations and reduces the number of
evaluations required on fine grids, while combining the samples at different
levels to arrive at an accurate estimate. The framework seamlessly integrates
the multiscale features of the GMsFEM with the multilevel feature of the MLMC
methods following the work in \cite{ketelson2013}, and our numerical
experiments illustrate its efficiency and accuracy in comparison with standard
Monte Carlo estimates.Comment: 29 pages, 6 figure
- …