Search Results

Now showing 1 - 6 of 6
  • Item
    Optimal and robust a posteriori error estimates in L∞(L2) for the approximation of Allen-Cahn equations past singularities
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2009) Bartels, Sören; Müller, Rüdiger
    Optimal a posteriori error estimates in L∞(L2) are derived for the finite element approximation of Allen-Cahn equations. The estimates depend on the inverse of a small parameter only in a low order polynomial and are valid past topological changes of the evolving interface. The error analysis employs an elliptic reconstruction of the approximate solution and applies to a large class of conforming, nonconforming, mixed, and discontinuous Galerkin methods. Numerical experiments illustrate the theoretical results.
  • Item
    Error control for the approximation of Allen-Cahn and Cahn-Hilliard equations with a logarithmic potential
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2010) Bartels, Sören; Müller, Rüdiger
    A fully computable upper bound for the finite element approximation error of Allen-Cahn and Cahn-Hilliard equations with logarithmic potentials is derived. Numerical experiments show that for the sharp interface limit this bound is robust past topological changes. Modifications of the abstract results to derive quasi-optimal error estimates in different norms for lowest order finite element methods are discussed and lead to weaker conditions on the residuals under which the conditional error estimates hold.
  • Item
    Numerical approach to a model for quasistatic damage with spatial BV-regularization
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Bartels, Sören; Milicevic, Marijo; Thomas, Marita
    We address a model for rate-independent, partial, isotropic damage in quasistatic small strain linear elasticity, featuring a damage variable with spatial BV-regularization. Discrete solutions are obtained using an alternate time-discrete scheme and the Variable-ADMM algorithm to solve the constrained nonsmooth optimization problem that determines the damage variable at each time step. We prove convergence of the method and show that discrete solutions approximate a semistable energetic solution of the rate-independent system. Moreover, we present our numerical results for two benchmark problems.
  • Item
    Fully discrete approximation of rate-independent damage models with gradient regularization
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Bartels, Sören; Milicevic, Marijo; Thomas, Marita; Weber, Nico
    This work provides a convergence analysis of a time-discrete scheme coupled with a finite-element approximation in space for a model for partial, rate-independent damage featuring a gradient regularization as well as a non-smooth constraint to account for the unidirectionality of the damage evolution. The numerical algorithm to solve the coupled problem of quasistatic small strain linear elasticity with rate-independent gradient damage is based on a Variable ADMM-method to approximate the nonsmooth contribution. Space-discretization is based on P1 finite elements and the algorithm directly couples the time-step size with the spatial grid size h. For a wide class of gradient regularizations, which allows both for Sobolev functions of integrability exponent r ∈ (1, ∞) and for BV-functions, it is shown that solutions obtained with the algorithm approximate as h → 0 a semistable energetic solution of the original problem. The latter is characterized by a minimality property for the displacements, a semistability inequality for the damage variable and an energy dissipation estimate. Numerical benchmark experiments confirm the stability of the method.
  • Item
    Quasistatic small-strain plasticity in the limit of vanishing hardening and its numerical approximation
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2010) Bartels, Sören; Mielke, Alexander; Roubíček, Thomáš
    The quasistatic rate-independent evolution of the Prager--Ziegler-type model of linearized plasticity with hardening is shown to converge to the rate-independent evolution of the Prandtl--Reuss elastic/perfectly plastic model. Based on the concept of energetic solutions we study the convergence of the solutions in the limit for hardening coefficients converging to 0 by using the abstract method of Gamma-convergence for rate-independent systems. An unconditionally convergent numerical scheme is devised and 2D and 3D numerical experiments are presented. A two-sided energy inequality is a posteriori verified to document experimental convergence rates.
  • Item
    Approximation schemes for materials with discontinuities
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Bartels, Sören; Milicevic, Marijo; Thomas, Marita; Tornquist, Sven; Weber, Nico
    Damage and fracture phenomena are related to the evolution of discontinuities both in space and in time. This contribution deals with methods from mathematical and numerical analysis to handle these: Suitable mathematical formulations and time-discrete schemes for problems with discontinuities in time are presented. For the treatment of problems with discontinuities in space, the focus lies on FE-methods for minimization problems in the space of functions of bounded variation. The developed methods are used to introduce fully discrete schemes for a rate-independent damage model and for the viscous approximation of a model for dynamic phase-field fracture. Convergence of the schemes is discussed.