Refine
Language
- English (8)
Keywords
- Cahn-Larche system (1)
- Dune (1)
- Schur-Newton multigrid (1)
- abstract interface (1)
- adaptive finite elements (1)
- biomechanics (1)
- c++ (1)
- generic programming (1)
- grid-based methods (1)
- hierarchical grids (1)
We present globally convergent multigrid methods for the nonsymmetric
obstacle problems as arising from the discretization of Black–Scholes models of
American options with local volatilities and discrete data. No tuning or regularization
parameters occur. Our approach relies on symmetrization by transformation
and data recovery by superconvergence.
The purpose of the paper is to apply monotone multigrid methods
to static and dynamic biomechanical contact problems.
In space, a finite element method involving a mortar
discretization of the contact conditions is used.
In time, a new contact--stabilized Newmark scheme is presented.
Numerical experiments for a two body Hertzian contact problem
and a biomechanical knee problem are reported.
We present and analyze novel hierarchical a posteriori error estimates
for self-adjoint elliptic obstacle problems.
Our approach differs from straightforward, but non-reliable estimators~\cite{RHWHoppe_RKornhuber_1994a}
by an additional extra term accounting for the deviation
of the discrete free boundary in the localization step.
We prove efficiency and reliability
on a saturation assumption and a regularity condition on the underlying grid.
Heuristic arguments suggest
that the extra term is of higher order and preserves full locality.
Numerical computations confirm our theoretical findings.
In this review, we intend to clarify the underlying ideas and the relations
between various multigrid methods ranging from subset decomposition,
to projected subspace decomposition and truncated multigrid.
In addition, we present a novel globally convergent inexact active set method
which is closely related to truncated multigrid. The numerical properties
of algorithms are carefully assessed by means of a degenerate problem and
a problem with a complicated coincidence set.
In a companion paper [Matheon-Preprint 403] we introduced an abstract definition of a parallel and adaptive hierarchical grid for scientific computing. Based on this
definition we derive an efficient interface specification as a set of C++ classes.
This interface separates the applications from the grid data structures.
Thus, user implementations become independent of the underlying grid
implementation. Modern C++ template techniques are used to provide an
interface implementation without big performance losses.
The implementation is realized as part of the
software environment DUNE.
Numerical tests demonstrate the flexibility and the efficiency of our approach.
We present a hierarchical a~posteriori error analysis for the
minimum value of the energy functional in symmetric obstacle
problems. The main result is that the energy of the exact solution
is, up to data oscillation, equivalent to an appropriate
hierarchical estimator. The proof of the main result does not
invoke any saturation assumption. Moreover, we prove an a
posteriori error estimate indicating that the estimator from
\cite{RHWHoppe_RKornhuber_1994a} is asymptotically reliable and we
give sufficient conditions for the validity of a saturation
assumption. Finally, we corroborate and complement our theoretical
results with numerical experiments.
We introduce nonsmooth Schur-Newton methods for the solution of the nonlinear discrete saddle-point problems arising from discretized vector-valued Cahn-Hilliard equations with logarithmic and obstacle potentials. The discrete problems are obtained by semi-implicit discretization in time and a first order finite element discretization in space. We incorporate the linear constraints that enforce solutions to stay on the Gibbs simplex using Lagrangian multipliers and prove existence of these multipliers under the assumption of a non-trivial initial condition for the order parameters.
Coarsening in solder alloys is a widely accepted indicator for
possible failure of joints in electronic devices.
Based on the well-established Cahn--Larch\'e model
with logarithmic chemical energy density (see, e.g., Dreyer and Mueller 2001),
we present a numerical environment for the efficient
and reliable simulation of coarsening in binary alloys.
Main features are adaptive mesh refinement based on hierarchical error estimates,
fast and reliable algebraic solution by multigrid and Schur--Newton multigrid methods,
and the quantification of the coarsening speed by the temporal growth of mean phase radii.
We provide a detailed description and a numerical assessment of the algorithm and its different components,
together with a practical application to a eutectic AgCu brazing alloy.