382 research outputs found
Solving the Poisson equation on small aspect ratio domains using unstructured meshes
We discuss the ill conditioning of the matrix for the discretised Poisson
equation in the small aspect ratio limit, and motivate this problem in the
context of nonhydrostatic ocean modelling. Efficient iterative solvers for the
Poisson equation in small aspect ratio domains are crucial for the successful
development of nonhydrostatic ocean models on unstructured meshes. We introduce
a new multigrid preconditioner for the Poisson problem which can be used with
finite element discretisations on general unstructured meshes; this
preconditioner is motivated by the fact that the Poisson problem has a
condition number which is independent of aspect ratio when Dirichlet boundary
conditions are imposed on the top surface of the domain. This leads to the
first level in an algebraic multigrid solver (which can be extended by further
conventional algebraic multigrid stages), and an additive smoother. We
illustrate the method with numerical tests on unstructured meshes, which show
that the preconditioner makes a dramatic improvement on a more standard
multigrid preconditioner approach, and also show that the additive smoother
produces better results than standard SOR smoothing. This new solver method
makes it feasible to run nonhydrostatic unstructured mesh ocean models in small
aspect ratio domains.Comment: submitted to Ocean Modellin
A fully-coupled discontinuous Galerkin method for two-phase flow in porous media with discontinuous capillary pressure
In this paper we formulate and test numerically a fully-coupled discontinuous
Galerkin (DG) method for incompressible two-phase flow with discontinuous
capillary pressure. The spatial discretization uses the symmetric interior
penalty DG formulation with weighted averages and is based on a wetting-phase
potential / capillary potential formulation of the two-phase flow system. After
discretizing in time with diagonally implicit Runge-Kutta schemes the resulting
systems of nonlinear algebraic equations are solved with Newton's method and
the arising systems of linear equations are solved efficiently and in parallel
with an algebraic multigrid method. The new scheme is investigated for various
test problems from the literature and is also compared to a cell-centered
finite volume scheme in terms of accuracy and time to solution. We find that
the method is accurate, robust and efficient. In particular no post-processing
of the DG velocity field is necessary in contrast to results reported by
several authors for decoupled schemes. Moreover, the solver scales well in
parallel and three-dimensional problems with up to nearly 100 million degrees
of freedom per time step have been computed on 1000 processors
A robust adaptive algebraic multigrid linear solver for structural mechanics
The numerical simulation of structural mechanics applications via finite
elements usually requires the solution of large-size and ill-conditioned linear
systems, especially when accurate results are sought for derived variables
interpolated with lower order functions, like stress or deformation fields.
Such task represents the most time-consuming kernel in commercial simulators;
thus, it is of significant interest the development of robust and efficient
linear solvers for such applications. In this context, direct solvers, which
are based on LU factorization techniques, are often used due to their
robustness and easy setup; however, they can reach only superlinear complexity,
in the best case, thus, have limited applicability depending on the problem
size. On the other hand, iterative solvers based on algebraic multigrid (AMG)
preconditioners can reach up to linear complexity for sufficiently regular
problems but do not always converge and require more knowledge from the user
for an efficient setup. In this work, we present an adaptive AMG method
specifically designed to improve its usability and efficiency in the solution
of structural problems. We show numerical results for several practical
applications with millions of unknowns and compare our method with two
state-of-the-art linear solvers proving its efficiency and robustness.Comment: 50 pages, 16 figures, submitted to CMAM
Master of Science
thesisThe algebraic multigrid (AMG) method is often used as a preconditioner in Krylov subspace solvers such as the conjugate gradient method. An AMG preconditioner hierarchically aggregates the degrees of freedom during the coarsening phase in order to eciently account for lower-frequency errors. Each degree of freedom in the coarser level corresponds to one of the aggregates in the ner level. The aggregation in each level in the hierarchy has a signicant impact on the eectiveness of AMG as a preconditioner. The aggregation can be formulated as a partitioning problem on the graph induced from the matrix representation of a linear system. The contributions of this work are as follows: rst, a GPU implementation of a \bottom-up" partitioning scheme based on maximal independent sets (MIS), including an ecient conditioning scheme for enforcing partition size constraints; second, three novel topological metrics, convexity, eccentricity, and minimum enclosing ball, for measuring partition quality; third, empirical test results comparing our MIS-Based aggregation methods with the MeTis graph partioning library, showing that the metrics correlate more strongly with AMG performance than the commonly used edge-cut metric, and that for ner aggregations, MIS-based aggregation is better suited for AMG coarsening than is the \top down" MeTis graph partitioning library, but that for coarser aggregations, MeTis performs better
Shifted Laplacian multigrid for the elastic Helmholtz equation
The shifted Laplacian multigrid method is a well known approach for
preconditioning the indefinite linear system arising from the discretization of
the acoustic Helmholtz equation. This equation is used to model wave
propagation in the frequency domain. However, in some cases the acoustic
equation is not sufficient for modeling the physics of the wave propagation,
and one has to consider the elastic Helmholtz equation. Such a case arises in
geophysical seismic imaging applications, where the earth's subsurface is the
elastic medium. The elastic Helmholtz equation is much harder to solve than its
acoustic counterpart, partially because it is three times larger, and partially
because it models more complicated physics. Despite this, there are very few
solvers available for the elastic equation compared to the array of solvers
that are available for the acoustic one. In this work we extend the shifted
Laplacian approach to the elastic Helmholtz equation, by combining the complex
shift idea with approaches for linear elasticity. We demonstrate the efficiency
and properties of our solver using numerical experiments for problems with
heterogeneous media in two and three dimensions
- …