5,221 research outputs found
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
Fast finite difference solvers for singular solutions of the elliptic Monge-Amp\`ere equation
The elliptic Monge-Ampere equation is a fully nonlinear Partial Differential
Equation which originated in geometric surface theory, and has been applied in
dynamic meteorology, elasticity, geometric optics, image processing and image
registration. Solutions can be singular, in which case standard numerical
approaches fail. In this article we build a finite difference solver for the
Monge-Ampere equation, which converges even for singular solutions. Regularity
results are used to select a priori between a stable, provably convergent
monotone discretization and an accurate finite difference discretization in
different regions of the computational domain. This allows singular solutions
to be computed using a stable method, and regular solutions to be computed more
accurately. The resulting nonlinear equations are then solved by Newton's
method. Computational results in two and three dimensions validate the claims
of accuracy and solution speed. A computational example is presented which
demonstrates the necessity of the use of the monotone scheme near
singularities.Comment: 23 pages, 4 figures, 4 tables; added arxiv links to references, added
coment
Easy Integral Surfaces: A Fast, Quad-based Stream and Path Surface Algorithm
a fast, quad-based stream and path surface algorith
Multilevel Solvers for Unstructured Surface Meshes
Parameterization of unstructured surface meshes is of fundamental importance in many applications of digital geometry processing. Such parameterization approaches give rise to large and exceedingly ill-conditioned systems which are difficult or impossible to solve without the use of sophisticated multilevel preconditioning strategies. Since the underlying meshes are very fine to begin with, such multilevel preconditioners require mesh coarsening to build an appropriate hierarchy. In this paper we consider several strategies for the construction of hierarchies using ideas from mesh simplification algorithms used in the computer graphics literature. We introduce two novel hierarchy construction schemes and demonstrate their superior performance when used in conjunction with a multigrid preconditioner
- …