| dc.identifier.uri | http://hdl.handle.net/11401/77232 | |
| dc.description.sponsorship | This work is sponsored by the Stony Brook University Graduate School in compliance with the requirements for completion of degree. | en_US |
| dc.format | Monograph | |
| dc.format.medium | Electronic Resource | en_US |
| dc.language.iso | en_US | |
| dc.publisher | The Graduate School, Stony Brook University: Stony Brook, NY. | |
| dc.type | Dissertation | |
| dcterms.abstract | Iterative methods are some of the most important techniques in solving large scale linear systems. Compared to direct methods, iterative solvers have the advantage of lower storage cost and better scalability as the problem size increases. Among the methods, multigrid solvers and multigrid preconditioners are often the optimal choices for solving systems that arise from partial differential equations (PDEs). In this dissertation, we introduce several efficient algorithms for various scientific applications. Our first algorithm is a specialized geometric multigrid solver for ill-conditioned systems from Helmholtz equations. Such equations appear in climate models with pure Neumann boundary condition and small wave numbers. In numerical linear algebra, ill-conditioned and even singular systems are inherently hard to solve. Many standard methods are either slow or non-convergent. We demonstrate that our solver delivers accurate solutions with fast convergence. The second algorithm, HyGA, is a general hybrid geometric+algebraic multigrid framework for elliptic type PDEs. It leverages the rigor, accuracy and efficiency of geometric multigrid (GMG) for hierarchical unstructured meshes, with the flexibility of algebraic multigrid (AMG) at the coarsest level. We conduct numerical experiments using Poisson equations in both 2-D and 3-D, and demonstrate the advantage of HyGA over classical GMG and AMG. Besides the aforementioned algorithms, we introduce an orthogonally projected implicit null-space method (OPINS) for saddle point systems. The traditional null-space method is inefficient because it is expensive to find the null-space explicitly. Some alternatives, notably constraint-preconditioned or projected Krylov methods, are relatively efficient, but they can suffer from numerical instability. OPINS is equivalent to the null-space method with an orthogonal projector, without forming the orthogonal basis of the null space explicitly. Our results show that it is more stable than projected Krylov methods while achieving similar efficiency. | |
| dcterms.available | 2017-09-20T16:52:15Z | |
| dcterms.contributor | Jiao, Xiangmin | en_US |
| dcterms.contributor | Reuter, Matthew | en_US |
| dcterms.contributor | Khairoutdinov, Marat | en_US |
| dcterms.contributor | Samulyak, Roman | en_US |
| dcterms.contributor | Zhang, Minghua | en_US |
| dcterms.contributor | . | en_US |
| dcterms.creator | Lu, Cao | |
| dcterms.dateAccepted | 2017-09-20T16:52:15Z | |
| dcterms.dateSubmitted | 2017-09-20T16:52:15Z | |
| dcterms.description | Department of Applied Mathematics and Statistics | en_US |
| dcterms.extent | 116 pg. | en_US |
| dcterms.format | Application/PDF | en_US |
| dcterms.format | Monograph | |
| dcterms.identifier | http://hdl.handle.net/11401/77232 | |
| dcterms.issued | 2016-12-01 | |
| dcterms.language | en_US | |
| dcterms.provenance | Made available in DSpace on 2017-09-20T16:52:15Z (GMT). No. of bitstreams: 1
Lu_grad.sunysb_0771E_12947.pdf: 1913317 bytes, checksum: 3eb3ce58268b2a1ec131e7918b96315c (MD5)
Previous issue date: 1 | en |
| dcterms.publisher | The Graduate School, Stony Brook University: Stony Brook, NY. | |
| dcterms.subject | Applied mathematics | |
| dcterms.subject | geometric+algebraic multigrid, KKT systems, Krylov subspace methods, multigrid methods, null-space method, singular systems | |
| dcterms.title | Efficient Iterative and Multigrid Solvers with Applications | |
| dcterms.type | Dissertation | |
