We describe the University of Florida Sparse Matrix Collection, a large and actively growing set of sparse matrices that arise in real applications The Collection is widely used by the numerical linear algebra community for the development and performance evaluation of sparse matrix algorithms It allows for robust and repeatable experiments: robust because performance results with artificially generated matrices can be misleading, and repeatable because matrices are curated and made publicly available in many formats Its matrices cover a wide spectrum of domains, include those arising from problems with underlying 2D or 3D geometry (as structural engineering, computational fluid dynamics, model reduction, electromagnetics, semiconductor devices, thermodynamics, materials, acoustics, computer graphics/vision, robotics/kinematics, and other discretizations) and those that typically do not have such geometry (optimization, circuit simulation, economic and financial modeling, theoretical and quantum chemistry, chemical process simulation, mathematics and statistics, power networks, and other networks and graphs) We provide software for accessing and managing the Collection, from MATLAB™, Mathematica™, Fortran, and C, as well as an online search capability Graph visualization of the matrices is provided, and a new multilevel coarsening scheme is proposed to facilitate this task

The university of Florida sparse matrix collection

Large linear systems of saddle point type arise in a wide variety of applications throughout computational science and engineering. Due to their indefiniteness and often poor spectral properties, such linear systems represent a significant challenge for solver developers. In recent years there has been a surge of interest in saddle point problems, and numerous solution techniques have been proposed for this type of system. The aim of this paper is to present and discuss a large selection of solution methods for linear systems in saddle point form, with an emphasis on iterative methods for large and sparse problems.

/pdf/numerical-solution-of-saddle-point-problems-l583iob5c2.pdf

Numerical solution of saddle point problems

5. M. Green, J. Schwarz, and E. Witten, Superstring theory, Cambridge Univ. Press, 1987. 6. J. Isenberg, P. Yasskin, and P. Green, Non-self-dual gauge fields, Phys. Lett. 78B (1978), 462-464. 7. B. Kostant, Graded manifolds, graded Lie theory, and prequantization, Differential Geometric Methods in Mathematicas Physics, Lecture Notes in Math., vol. 570, SpringerVerlag, Berlin and New York, 1977. 8. C. LeBrun, Thickenings and gauge fields, Class. Quantum Grav. 3 (1986), 1039-1059. 9. , Thickenings and conformai gravity, preprint, 1989. 10. C. LeBrun and M. Rothstein, Moduli of super Riemann surfaces, Commun. Math. Phys. 117(1988), 159-176. 11. Y. Manin, Critical dimensions of string theories and the dualizing sheaf on the moduli space of (super) curves, Funct. Anal. Appl. 20 (1987), 244-245. 12. R. Penrose and W. Rindler, Spinors and space-time, V.2, spinor and twistor methods in space-time geometry, Cambridge Univ. Press, 1986. 13. R. Ward, On self-dual gauge fields, Phys. Lett. 61A (1977), 81-82. 14. E. Witten, An interpretation of classical Yang-Mills theory, Phys. Lett. 77NB (1978), 394-398. 15. , Twistor-like transform in ten dimensions, Nucl. Phys. B266 (1986), 245-264. 16. , Physics and geometry, Proc. Internat. Congr. Math., Berkeley, 1986, pp. 267302, Amer. Math. Soc, Providence, R.I., 1987.

Infinite-Dimensional Dynamical Systems in Mechanics and Physics (Roger Temam)

An overview of the extended/generalized finite element method (GEFM/XFEM) with emphasis on methodological issues is presented. This method enables the accurate approximation of solutions that involve jumps, kinks, singularities, and other locally non-smooth features within elements. This is achieved by enriching the polynomial approximation space of the classical finite element method. The GEFM/XFEM has shown its potential in a variety of applications that involve non-smooth solutions near interfaces: Among them are the simulation of cracks, shear bands, dislocations, solidification, and multi-field problems. Copyright © 2010 John Wiley & Sons, Ltd.

The extended/generalized finite element method: An overview of the method and its applications

http://www.mathcs.emory.edu/%7Ebenzi/Web_papers/survey.pdf

Preconditioning techniques for large linear systems: a survey

https://www.igpm.rwth-aachen.de/Download/reports/pdf/IGPM264.pdf

An extended pressure finite element space for two-phase incompressible flows with surface tension

In this paper a new finite element approach for the discretization of elliptic partial differential equations on surfaces is treated. The main idea is to use finite element spaces that are induced by triangulations of an “outer” domain to discretize the partial differential equation on the surface. The method is particularly suitable for problems in which there is a coupling with a flow problem in an outer domain that contains the surface. We give an analysis that shows that the method has optimal order of convergence both in the $H^1$- and in the $L^2$-norm. Results of numerical experiments are included that confirm this optimality.

/pdf/a-finite-element-method-for-elliptic-equations-on-surfaces-3vxx0v5zzt.pdf

A Finite Element Method for Elliptic Equations on Surfaces

Introduction.- Part I One-phase incompressible flows.- Mathematical models.- Finite element discretization.- Time integration.-

Numerical Methods for Two-phase Incompressible Flows

In this paper a stabilizing augmented Lagrangian technique for the Stokes equations is studied. The method is consistent and hence does not change the continuous solution. We show that this stabilization improves the well-posedness of the continuous problem for small values of the viscosity coefficient. We analyze the influence of this stabilization on the accuracy of the finite element solution and on the convergence properties of the inexact Uzawa method.

/pdf/grad-div-stablilization-for-stokes-equations-tyyuugxiip.pdf

Grad-div stablilization for Stokes equations

We present a method that has been developed for the efficient numerical simulation of two-phase incompressible flows. For capturing the interface between the phases the level set technique is applied. The continuous model consists of the incompressible Navier–Stokes equations coupled with an advection equation for the level set function. The effect of surface tension is modeled by a localized force term at the interface (so-called continuum surface force approach). For spatial discretization of velocity, pressure and the level set function conforming finite elements on a hierarchy of nested tetrahedral grids are used. In the finite element setting we can apply a special technique to the localized force term, which is based on a partial integration rule for the Laplace–Beltrami operator. Due to this approach the second order derivatives coming from the curvature can be eliminated. For the time discretization we apply a variant of the fractional step θ-scheme. The discrete saddle point problems that occur in each time step are solved using an inexact Uzawa method combined with multigrid techniques. For reparametrization of the level set function a new variant of the fast marching method is introduced. A special feature of the solver is that it combines the level set method with finite element discretization, Laplace–Beltrami partial integration, multilevel local refinement and multigrid solution techniques. All these components of the solver are described. Results of numerical experiments are presented.

/pdf/a-finite-element-based-level-set-method-for-two-phase-ralws0ygav.pdf

Arnold Reusken

Papers

An extended pressure finite element space for two-phase incompressible flows with surface tension

A Finite Element Method for Elliptic Equations on Surfaces

Numerical Methods for Two-phase Incompressible Flows

Grad-div stablilization for Stokes equations

A finite element based level set method for two-phase incompressible flows