A second order splitting method for the Cahn-Hilliard equation

From MaRDI portal





The Cahn-Hilliard equation models the phase separation which occurs upon the quenching of a binary mixture into the spinodal region: \[ (1)\quad \partial u/\partial t=\Delta \phi (u)-\gamma \Delta^ 2u\quad x\in \Omega,\quad t>0, \] plus initial and boundary conditions and with \(\phi (u)=\psi '(u);\quad \psi (u)=\gamma_ 2(u^ 4/4)+\gamma_ 1(u^ 4/4)+\gamma_ 1(u^ 3/3)+(\gamma_ 0/2)u^ 2;\) \(\gamma_ 2>0\). Here u(x,t) means the concentration of one component of such a mixture. By introducing the chemical potential \(w=\phi (u)-\gamma \Delta u,\) the author substitutes to (1) a second order system of partial differential equations. Then, he considers a semi-discrete finite element method (in space) which requires only continuous elements. The optimal order error bounds are derived in various norms for an implementation which uses mass lumping. The continuous problem has an energy based Lyapunov functional and it proved that this property still holds for the discrete problem.




Cited in
(only showing first 100 items - show all)








This page was built for publication: A second order splitting method for the Cahn-Hilliard equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1118373)