Convergent finite element discretizations of the nonstationary incompressible magnetohydrodynamics system

The incompressible MHD equations couple Navier-Stokes equations with Maxwell's equations to describe the flow of a viscous, incompressible, and electrically conducting fluid in a Lipschitz domain $\Omega \subset \mathbb{R}^3$. We verify convergence of iterates of different coupling and decoupling fully discrete schemes towards weak solutions for vanishing discretization parameters. Optimal first order of convergence is shown in the presence of strong solutions for a splitting scheme which decouples the computation of velocity field, pressure, and magnetic fields at every iteration step

Optimal control for the thin-film equation: Convergence of a multi-parameter approach to track state constraints avoiding degeneracies

We consider an optimal control problem subject to the thin-film equation which is deduced from the Navier--Stokes equation. The PDE constraint lacks well-posedness for general right-hand sides due to possible degeneracies; state constraints are used to circumvent this problematic issue and to ensure well-posedness, and the rigorous derivation of necessary optimality conditions for the optimal control problem is performed. A multi-parameter regularization is considered which addresses both, the possibly degenerate term in the equation and the state constraint, and convergence is shown for vanishing regularization parameters by decoupling both effects. The fully regularized optimal control problem allows for practical simulations which are provided, including the control of a dewetting scenario, to evidence the need of the state constraint, and to motivate proper scalings of involved regularization and numerical parameters

Optimal Strong Rates of Convergence for a Space-Time Discretization of the Stochastic Allen-Cahn Equation with multiplicative noise

The stochastic Allen-Cahn equation with multiplicative noise involves the nonlinear drift operator ${\mathscr A}(x) = \Delta x - \bigl(\vert x\vert^2 -1\bigr)x$. We use the fact that ${\mathscr A}(x) = -{\mathcal J}^{\prime}(x)$ satisfies a weak monotonicity property to deduce uniform bounds in strong norms for solutions of the temporal, as well as of the spatio-temporal discretization of the problem. This weak monotonicity property then allows for the estimate $\underset{1 \leq j \leq J}\sup {\mathbb E}\bigl[ \Vert X_{t_j} - Y^j\Vert_{{\mathbb L}^2}^2\bigr] \leq C_{\delta}(k^{1-\delta} + h^2)$ for all small $\delta>0$, where $X$ is the strong variational solution of the stochastic Allen-Cahn equation, while $\big\{Y^j:0\le j\le J\big\}$ solves a structure preserving finite element based space-time discretization of the problem on a temporal mesh $\{ t_j;\, 1 \leq j \leq J\}$ of size $k>0$ which covers $[0,T]$

Multiscale resolution in the computation of crystalline microstructure

Summary.: This paper addresses the numerical approximation of microstructures in crystalline phase transitions without surface energy. It is shown that branching of different variants near interfaces of twinned martensite and austenite phases leads to reduced energies in finite element approximations. Such behavior of minimizing deformations is understood for an extended model that involves surface energies. Moreover, the closely related question of the role of different growth conditions of the employed bulk energy is discussed. By explicit construction of discrete deformations in lowest order finite element spaces we prove upper bounds for the energy and thereby clarify the question of the dependence of the convergence rate upon growth conditions and lamination orders. For first order laminates the estimates are optima

Stabilization methods in relaxed micromagnetism

The magnetization of a ferromagnetic sample solves a non-convex variational problem, where its relaxation by convexifying the energy density resolves relevant macroscopic information. The numerical analysis of the relaxed model has to deal with a constrained convex but degenerated, nonlocal energy functional in mixed formulation for magnetic potential u and magnetization m. In [C. Carstensen and A. Prohl, Numer. Math. 90 (2001) 65–99], the conforming P1 - (P0)d-element in d=2,3 spatial dimensions is shown to lead to an ill-posed discrete problem in relaxed micromagnetism, and suboptimal convergence. This observation motivated a non-conforming finite element method which leads to a well-posed discrete problem, with solutions converging at optimal rate. In this work, we provide both an a priori and a posteriori error analysis for two stabilized conforming methods which account for inter-element jumps of the piecewise constant magnetization. Both methods converge at optimal rate; the new approach is applied to a macroscopic nonstationary ferromagnetic model [M. Kružík and A. Prohl, Adv. Math. Sci. Appl. 14 (2004) 665–681 – M. Kružík and T. Roubíček, Z. Angew. Math. Phys. 55 (2004) 159–182 ]

Error analysis of a mixed finite element method for the Cahn-Hilliard equation

Summary.: We propose and analyze a semi-discrete and a fully discrete mixed finite element method for the Cahn-Hilliard equation ut + Δ(ɛΔu−ɛ−1f(u)) = 0, where ɛ > 0 is a small parameter. Error estimates which are quasi-optimal order in time and optimal order in space are shown for the proposed methods under minimum regularity assumptions on the initial data and the domain. In particular, it is shown that all error bounds depend on only in some lower polynomial order for small ɛ. The cruxes of our analysis are to establish stability estimates for the discrete solutions, to use a spectrum estimate result of Alikakos and Fusco [2], and Chen [15] to prove a discrete counterpart of it for a linearized Cahn-Hilliard operator to handle the nonlinear term on a stretched time grid. The ideas and techniques developed in this paper also enable us to prove convergence of the fully discrete finite element solution to the solution of the Hele-Shaw (Mullins-Sekerka) problem as ɛ → 0 in [29

Rate of convergence of regularization procedures and finite element approximations for the total variation flow

Summary: We derive rates of convergence for regularization procedures (characterized by a parameter ɛ) and finite element approximations of the total variation flow, which arises from image processing, geometric analysis and materials sciences. Practically useful error estimates, which depend on only in low polynomial orders, are established for the proposed fully discrete finite element approximations. As a result, scaling laws which relate mesh parameters to the regularization parameter are also obtained. Numerical experiments are provided to validate the theoretical results and show efficiency of the proposed numerical method

Convergence of a fully discrete finite element method for a degenerate parabolic system modelling nematic liquid crystals with variable degree of orientation

Published versio

A convergent and constraint-preserving finite element method for the p-harmonic flow into spheres

Published versio