Browsing by Subject "Preconditioner"
Now showing 1 - 3 of 3
Results Per Page
Sort Options
Item Fast parallel solution of heterogeneous 3D time-harmonic wave equations(2012-12) Poulson, Jack Lesly; Ying, Lexing; Engquist, Bjorn; Fomel, Sergey; Ghattas, Omar; van de Geijn, RobertSeveral advancements related to the solution of 3D time-harmonic wave equations are presented, especially in the context of a parallel moving-PML sweeping preconditioner for problems without large-scale resonances. The main contribution of this dissertation is the introduction of an efficient parallel sweeping preconditioner and its subsequent application to several challenging velocity models. For instance, 3D seismic problems approaching a billion degrees of freedom have been solved in just a few minutes using several thousand processors. The setup and application costs of the sequential algorithm were also respectively refined to O(γ^2 N^(4/3)) and O(γ N log N), where N denotes the total number of degrees of freedom in the 3D volume and γ(ω) denotes the modestly frequency-dependent number of grid points per Perfectly Matched Layer discretization. Furthermore, high-performance parallel algorithms are proposed for performing multifrontal triangular solves with many right-hand sides, and a custom compression scheme is introduced which builds upon the translation invariance of free-space Green’s functions in order to justify the replacement of each dense matrix within a certain modified multifrontal method with the sum of a small number of Kronecker products. For the sake of reproducibility, every algorithm exercised within this dissertation is made available as part of the open source packages Clique and Parallel Sweeping Preconditioner (PSP).Item Multiscale mortar mixed finite element methods for flow problems in highly heterogeneous porous media(2013-12) Xiao, Hailong; Arbogast, Todd James, 1957-We use Darcy's law and conservation of mass to model the flow of a fluid through a porous medium. It is a second order elliptic system with a heterogeneous coefficient. We consider the equations written in mixed form. In the heterogeneous case, we define a new multiscale mortar space that incorporates purely local information from homogenization theory to better approximate the solution along the interfaces with just a few degrees of freedom. In the case of a locally periodic heterogeneous coefficient of period epsilon, we prove that the new method achieves both optimal order error estimates in the discretization parameters and good approximation when epsilon is small. Moreover, we present numerical examples to assess its performance when the coefficient is not obviously locally periodic. We show that the new mortar method works well, and better than polynomial mortar spaces. On the other hand, we also propose to use multiscale mortars as a coarse component to construct a two-level preconditioner for the saddle point linear system arising from the fine scale discretization of the mixed finite element system. The two-level preconditioners are constructed based on the interfaces. We propose a framework to define the interpolation operators for the face based two-level preconditioners for different combination of coarse and fine scale mortar spaces for matching and nonmatching grids. In this dissertation, we show that for quasi-homogeneous problems and matching grids, the condition number of the preconditioned interface operator is bounded by (log(H/h))², which is the same as the traditional two-level preconditioners, for quasi-homogeneous problems. We show several numerical examples to demonstrate that for the strongly heterogeneous porous media, it is often desirable and even necessary to use a higher dimensional coarse mortar space to construct the coarse preconditioner to achieve convergence. We apply our ideas to study slightly compressible single phase and two-phase flow in a porous medium. We find that for the nonlinear single phase problem, the two-level preconditioners could be successfully applied to the symmetrized linear system. For the two-phase problem, using the fine scale, instead of multiscale, velocity solutions from the flow problem can greatly benefit the transport problem.Item Preconditioning for the mixed formulation of linear plane elasticity(Texas A&M University, 2005-11-01) Wang, YanqiuIn this dissertation, we study the mixed finite element method for the linear plane elasticity problem and iterative solvers for the resulting discrete system. We use the Arnold-Winther Element in the mixed finite element discretization. An overlapping Schwarz preconditioner and a multigrid preconditioner for the discrete system are developed and analyzed. We start by introducing the mixed formulation (stress-displacement formulation) for the linear plane elasticity problem and its discretization. A detailed analysis of the Arnold-Winther Element is given. The finite element discretization of the mixed formulation leads to a symmetric indefinite linear system. Next, we study efficient iterative solvers for the symmetric indefinite linear system which arises from the mixed finite element discretization of the linear plane elasticity problem. The preconditioned Minimum Residual Method is considered. It is shown that the problem of constructing a preconditioner for the indefinite linear system can be reduced to the problem of constructing a preconditioner for the H(div) problem in the Arnold-Winther finite element space. Our main work involves developing an overlapping Schwarz preconditioner and a multigrid preconditioner for the H(div) problem. We give condition number estimates for the preconditioned systems together with supporting numerical results.