104 research outputs found
Discretizations & Efficient Linear Solvers for Problems Related to Fluid Flow
Numerical solutions to fluid flow problems involve solving the linear systems arising from the discretization of the Stokes equation or a variant of it, which often have a saddle point structure and are difficult to solve. Geometric multigrid is a parallelizable method that can efficiently solve these linear systems especially for a large number of unknowns. We consider two approaches to solve these linear systems using geometric multigrid: First, we use a block preconditioner and apply geometric multigrid as in inner solver to the velocity block only. We develop deal.II tutorial step-56 to compare the use of geometric multigrid to other popular alternatives. This method is found to be competitive in serial computations in terms of performance and memory usage. Second, we design a special smoother to apply multigrid to the whole linear system. This smoother is analyzed as a Schwarz method using conforming and inf-sup stable discretization spaces. The resulting method is found to be competitive to a similar multigrid method using non-conforming finite elements that were studied by Kanschat and Mao. This approach has the potential to be superior to the first approach. Finally, expanding on the research done by Dannberg and Heister, we explore the analysis of a three-field Stokes formulation that is used to describe melt migration in the earth\u27s mantle. Multiple discretizations were studied to find the best one to use in the ASPECT software package. We also explore improvements to ASPECT\u27s linear solvers for this formulation utilizing block preconditioners and algebraic multigrid
Iterative solution of saddle point problems using divergence-free finite elements with applications to groundwater flow
Available from British Library Document Supply Centre-DSC:DXN041433 / BLDSC - British Library Document Supply CentreSIGLEGBUnited Kingdo
Iterative methods for heterogeneous media
EThOS - Electronic Theses Online ServiceGBUnited Kingdo
Hybridised multigrid preconditioners for a compatible finite element dynamical core
Compatible finite element discretisations for the atmospheric equations of
motion have recently attracted considerable interest. Semi-implicit
timestepping methods require the repeated solution of a large saddle-point
system of linear equations. Preconditioning this system is challenging since
the velocity mass matrix is non-diagonal, leading to a dense Schur complement.
Hybridisable discretisations overcome this issue: weakly enforcing continuity
of the velocity field with Lagrange multipliers leads to a sparse system of
equations, which has a similar structure to the pressure Schur complement in
traditional approaches. We describe how the hybridised sparse system can be
preconditioned with a non-nested two-level preconditioner. To solve the coarse
system, we use the multigrid pressure solver that is employed in the
approximate Schur complement method previously proposed by the some of the
authors. Our approach significantly reduces the number of solver iterations.
The method shows excellent performance and scales to large numbers of cores in
the Met Office next-generation climate- and weather prediction model LFRic.Comment: 24 pages, 13 figures, 5 tables; accepted for publication in Quarterly
Journal of the Royal Meteorological Societ
KSPHPDDM and PCHPDDM: Extending PETSc with advanced Krylov methods and robust multilevel overlapping Schwarz preconditioners
[EN] Contemporary applications in computational science and engineering often require the solution of linear systems which may be of different sizes, shapes, and structures. The goal of this paper is to explain how two libraries, PETSc and HPDDM, have been interfaced in order to offer end-users robust overlapping Schwarz preconditioners and advanced Krylov methods featuring recycling and the ability to deal with multiple right-hand sides. The flexibility of the implementation is showcased and explained with minimalist, easy-to-run, and reproducible examples, to ease the integration of these algorithms into more advanced frameworks. The examples provided cover applications from eigenanalysis, elasticity, combustion, and electromagnetism.Jose E. Roman was supported by the Spanish Agencia Estatal de Investigacion (AEI) under project SLEPc-DA (PID2019-107379RB-I00)Jolivet, P.; Roman, JE.; Zampini, S. (2021). KSPHPDDM and PCHPDDM: Extending PETSc with advanced Krylov methods and robust multilevel overlapping Schwarz preconditioners. Computers & Mathematics with Applications. 84:277-295. https://doi.org/10.1016/j.camwa.2021.01.0032772958
Schnelle Löser für partielle Differentialgleichungen
[no abstract available
- …