# Solving the Cahn-Hilliard variational inequality with a semi-smooth Newton method

Luise Blank; Martin Butz; Harald Garcke

ESAIM: Control, Optimisation and Calculus of Variations (2011)

- Volume: 17, Issue: 4, page 931-954
- ISSN: 1292-8119

## Access Full Article

top## Abstract

top## How to cite

topBlank, Luise, Butz, Martin, and Garcke, Harald. "Solving the Cahn-Hilliard variational inequality with a semi-smooth Newton method." ESAIM: Control, Optimisation and Calculus of Variations 17.4 (2011): 931-954. <http://eudml.org/doc/221916>.

@article{Blank2011,

abstract = {
The Cahn-Hilliard variational inequality is a non-standard
parabolic variational inequality of fourth order for which
straightforward numerical
approaches cannot be applied. We propose a primal-dual active set
method which can be interpreted as a semi-smooth Newton method as
solution technique for the discretized Cahn-Hilliard variational
inequality. A (semi-)implicit Euler discretization is used in time
and a piecewise linear finite element discretization of splitting
type is used in space leading to a discrete variational inequality
of saddle point type in each time step. In each iteration of the
primal-dual active set method a linearized system resulting from the
discretization of two coupled elliptic equations which are defined
on different sets has to be solved. We show local
convergence of the primal-dual active set method and demonstrate its
efficiency with several numerical simulations.
},

author = {Blank, Luise, Butz, Martin, Garcke, Harald},

journal = {ESAIM: Control, Optimisation and Calculus of Variations},

keywords = {Cahn-Hilliard equation; active-set methods;
semi-smooth Newton methods; gradient flows; PDE-constraint
optimization; saddle point structure; primal-dual active set method; discrete variational inequality},

language = {eng},

month = {11},

number = {4},

pages = {931-954},

publisher = {EDP Sciences},

title = {Solving the Cahn-Hilliard variational inequality with a semi-smooth Newton method},

url = {http://eudml.org/doc/221916},

volume = {17},

year = {2011},

}

TY - JOUR

AU - Blank, Luise

AU - Butz, Martin

AU - Garcke, Harald

TI - Solving the Cahn-Hilliard variational inequality with a semi-smooth Newton method

JO - ESAIM: Control, Optimisation and Calculus of Variations

DA - 2011/11//

PB - EDP Sciences

VL - 17

IS - 4

SP - 931

EP - 954

AB -
The Cahn-Hilliard variational inequality is a non-standard
parabolic variational inequality of fourth order for which
straightforward numerical
approaches cannot be applied. We propose a primal-dual active set
method which can be interpreted as a semi-smooth Newton method as
solution technique for the discretized Cahn-Hilliard variational
inequality. A (semi-)implicit Euler discretization is used in time
and a piecewise linear finite element discretization of splitting
type is used in space leading to a discrete variational inequality
of saddle point type in each time step. In each iteration of the
primal-dual active set method a linearized system resulting from the
discretization of two coupled elliptic equations which are defined
on different sets has to be solved. We show local
convergence of the primal-dual active set method and demonstrate its
efficiency with several numerical simulations.

LA - eng

KW - Cahn-Hilliard equation; active-set methods;
semi-smooth Newton methods; gradient flows; PDE-constraint
optimization; saddle point structure; primal-dual active set method; discrete variational inequality

UR - http://eudml.org/doc/221916

ER -

## References

top- R.A. Adams, Sobolev spaces, Pure and Applied Mathematics65. Academic Press, New York-London (1975).
- L. Banas and R. Nürnberg, A multigrid method for the Cahn-Hilliard equation with obstacle potential. Appl. Math. Comput.213 (2009) 290–303.
- J.W. Barrett, J.F. Blowey and H. Garcke, Finite element approximation of the Cahn–Hilliard equation with degenerate mobility. SIAM J. Numer. Anal.37 (1999) 286–318.
- J.W. Barrett, R. Nürnberg and V. Styles, Finite element approximation of a void electromigration model. SIAM J. Numer. Anal.42 (2004) 738–772.
- L. Blank, H. Garcke, L. Sarbu and V. Styles, Primal-dual active set methods for Allen-Cahn variational inequalities with non-local constraints. Preprint SPP1253-09-01 (2009).
- J.F. Blowey and C.M. Elliott, The Cahn-Hilliard gradient theory for phase separation with nonsmooth free energy. I. Mathematical analysis. Eur. J. Appl. Math.2 (1991) 233–280.
- J.F. Blowey and C.M. Elliott, The Cahn-Hilliard gradient theory for phase separation with nonsmooth free energy. II. Numerical analysis. Eur. J. Appl. Math.3 (1992) 147–179.
- J.F. Blowey and C.M. Elliott, Curvature dependent phase boundary motion and parabolic double obstacle problems, in Degenerate Diffusions, W.-M. Ni, L.A. Peletier and J.L. Vazquez Eds., IMA Vol. Math. Appl.47, Springer, New York (1993) 19–60.
- J.F. Blowey and C.M. Elliott, A phase field model with a double obstacle potential, in Motion by mean curvature, G. Buttazzo and A. Visintin Eds., de Gruyter (1994) 1–22.
- J.W. Cahn and J.E. Hilliard, Free energy of a nonuniform system. I. Interfacial energy. J. Chem. Phys.28 (1958) 258–267.
- I. Capuzzo Dolcetta, S.F. Vita and R. March, Area-preserving curve-shortening flows: From phase separation to image processing. Interfaces and Free Boundaries4 (2002) 325–434.
- X. Chen, Global asymptotic limit of solutions of the Cahn-Hilliard equation. J. Differential Geom.44 (1996) 262–311.
- X. Chen, Z. Nashed and L. Qi, Smoothing methods and semismooth methods for nondifferentiable operator equations. SIAM J. Numer. Anal.38 (2000) 1200–1216.
- M. Copetti and C.M. Elliott, Numerical analysis of the Cahn-Hilliard equation with a logarithmic free energy. Numer. Math.63 (1992) 39–65.
- T.A. Davis, Algorithm 832: UMFPACK, an unsymmetric-pattern multifrontal method. ACM Trans. Math. Soft.30 (2003) 196–199.
- T.A. Davis, A column pre-ordering strategy for the unsymmetric-pattern multifrontal method. ACM Trans. Math. Soft.34 (2003) 165–195.
- T.A. Davis and I.S. Duff, An unsymmetric-pattern multifrontal method for sparse LU factorization. SIAM J. Matrix Anal. Appl.18 (1997) 140–158.
- I.S. Duff and J.K. Reid, The multifrontal solution of indefinite sparse symmetric linear. ACM Trans. Math. Soft.9 (1983) 302–325.
- C.M. Elliott, The Cahn-Hilliard model for the kinetics of phase separation, in Mathematical models for phase change problems, Internat. Ser. Numer. Math.88, Birkhäuser, Basel (1989).
- C.M. Elliott and A.R. Gardiner, One dimensional phase field computations, Numerical Analysis 1993, Proceedings of Dundee Conference, D.F. Griffiths and G.A. Watson Eds., Longman Scientific and Technical (1994) 56–74.
- C.M. Elliott and S. Luckhaus, A generalised diffusion equation for phase separation of a multi-component mixture with interfacial free energy. SFB 256, University of Bonn, Preprint 195 (1991).
- C.M. Elliott and J. Ockendon, Weak and Variational Methods for Moving Boundary Problems, Pitman Research Notes in Mathematics59. Pitman (1982).
- L.C. Evans, Partial differential equations, Graduate Studies in Mathematics19. American Mathematical Society, Providence (1998).
- A. Friedman, Variational principles and free-boundary problems – Pure and Applied Mathematics. John Wiley & Sons, Inc., New York (1982).
- H. Garcke, Mechanical effects in the Cahn-Hilliard model: A review on mathematical results, in Mathematical Methods and Models in phase transitions, A. Miranvielle Ed., Nova Science Publ. (2005) 43–77.
- C. Gräser, Analysis und Approximation der Cahn-Hilliard Gleichung mit Hindernispotential. Diplomarbeit, Freie Universität Berlin, Fachbereich Mathematik und Informatik (2004).
- C. Gräser and R. Kornhuber, On preconditioned Uzawa-type iterations for a saddle point problem with inequality constraints, in Domain decomposition methods in science and engineering XVI, Lect. Notes Comput. Sci. Eng.55, Springer, Berlin (2007) 91–102.
- C. Gräser and R. Kornhuber, Nonsmooth Newton methods for set-valued saddle point problems. SIAM J. Numer. Anal.47 (2009) 1251–1273.
- C. Gräser and R. Kornhuber, Multigrid methods for obstacle problems. J. Comput. Math.27 (2009) 1–44.
- M. Hintermüller, K. Ito and K. Kunisch, The primal-dual active set strategy as a semismooth Newton method. SIAM J. Optim.13 (2002) 865–888.
- K. Ito and K. Kunisch, Semi-smooth Newton methods for variational inequalities of the first kind. ESAIM: M2AN37 (2003) 41–62.
- B.M. Irons, A frontal solution scheme for finite element analysis. Int. J. Numer. Methods Eng.2 (1970) 5–32.
- D. Kinderlehrer and G. Stampacchia, An introduction to variational inequalities and their applications, Pure and Applied Mathematics88. Academic Press, Inc., New York-London (1980).
- E. Kuhl and D.W. Schmid, Computational modeling of mineral unmixing and growth: An application of the Cahn-Hilliard equation. Comp. Mech.39 (2007) 439–451.
- P.-L. Lions and B. Mercier, Splitting algorithms for the sum of two nonlinear operators. SIAM J. Numer. Anal.16 (1979) 964–979.
- J.W.H. Liu, The multifrontal method for sparse matrix solution: Theory and practice. SIAM Rev.34 (1992) 82–109.
- J. Lowengrub and L. Truskinovsky, Quasi-incompressible Cahn-Hilliard fluids and topological transitions. R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci.454 (1978) 2617–2654.
- A. Novick-Cohen, The Cahn-Hilliard equation: mathematical and modeling perspectives. Adv. Math. Sci. Appl.8 (1998) 965–985.
- R.L. Pego, Front migration in the nonlinear Cahn–Hilliard equation. Proc. Roy. Soc. London, Ser. A422 (1989) 116–133.
- A. Schmidt and K.G. Siebert, Design of adaptive finite element software: The finite element toolbox ALBERTA, Lect. Notes Comput. Sci. Eng.42. Springer, Berlin (2005).
- B. Stoth, Convergence of the Cahn-Hilliard equation to the Mullins-Sekerka problem in spherical symmetry. J. Diff. Equ.125 (1996) 154–183.
- S. Tremaine, On the origin of irregular structure in Saturn's rings. Ast. J.125 (2003) 894–901.
- F. Tröltzsch, Optimale Steuerung partieller Differentialgleichungen: Theorie, Verfahren und Anwendungen. Vieweg Verlag (2005).
- S. Zhou and M.Y. Wang, Multimaterial structural topology optimization with a generalized Cahn-Hilliard model of multiphase transition. Struct. Multidisc. Optim.33 (2007) 89–111.

## NotesEmbed ?

topTo embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.