1,867 research outputs found
Preconditioning for Allen-Cahn variational inequalities with non-local constraints
The solution of Allen-Cahn variational inequalities with mass constraints is of interest in many applications. This problem can be solved both in its scalar and vector-valued form as a PDE-constrained optimization problem by means of a primal-dual active set method. At the heart of this method lies the solution of linear systems in saddle point form. In this paper we propose the use of Krylov-subspace solvers and suitable preconditioners for the saddle point systems. Numerical results illustrate the competitiveness of this approach
A Multigrid Optimization Algorithm for the Numerical Solution of Quasilinear Variational Inequalities Involving the -Laplacian
In this paper we propose a multigrid optimization algorithm (MG/OPT) for the
numerical solution of a class of quasilinear variational inequalities of the
second kind. This approach is enabled by the fact that the solution of the
variational inequality is given by the minimizer of a nonsmooth energy
functional, involving the -Laplace operator. We propose a Huber
regularization of the functional and a finite element discretization for the
problem. Further, we analyze the regularity of the discretized energy
functional, and we are able to prove that its Jacobian is slantly
differentiable. This regularity property is useful to analyze the convergence
of the MG/OPT algorithm. In fact, we demostrate that the algorithm is globally
convergent by using a mean value theorem for semismooth functions. Finally, we
apply the MG/OPT algorithm to the numerical simulation of the viscoplastic flow
of Bingham, Casson and Herschel-Bulkley fluids in a pipe. Several experiments
are carried out to show the efficiency of the proposed algorithm when solving
this kind of fluid mechanics problems
A direct algorithm in some free boundary problems
In this paper we propose a new algorithm for the well known elliptic obstacle problem and for parabolic variational inequalities like one and two phase Stefan problem and of obstacle type. Our approach enters the category of fixed domain methods and solves just linear elliptic or parabolic equations and their discretization at each iteration. We prove stability and convergence properties. The approximating coincidence set is explicitly computed and it converges in the Hausdorff-Pompeiu sense to the searched geometry. In the numerical examples, the algorithm has a very fast convergence and the obtained solutions (including the free boundaries) are accurate
- …