Journal
COMPUTATIONAL MECHANICS
Volume 46, Issue 3, Pages 471-493Publisher
SPRINGER
DOI: 10.1007/s00466-010-0490-4
Keywords
Cahn-Hilliard equation; Fourth-order diffusion; Phase separation; Higher-order continuity; Natural element method
Ask authors/readers for more resources
We present a natural element method to treat higher-order spatial derivatives in the Cahn-Hilliard equation. The Cahn-Hilliard equation is a fourth-order nonlinear partial differential equation that allows to model phase separation in binary mixtures. Standard classical C-0-continuous finite element solutions are not suitable because primal variational formulations of fourth-order operators are well-defined and integrable only if the finite element basis functions are piecewise smooth and globally C-1-continuous. To ensure C-1-continuity, we develop a natural-element-based spatial discretization scheme. The C-1-continuous natural element shape functions are achieved by a transformation of the classical Farin interpolant, which is basically obtained by embedding Sibsons natural element coordinates in a Bernstein-Bezier surface representation of a cubic simplex. For the temporal discretization, we apply the (second-order accurate) trapezoidal time integration scheme supplemented with an adaptively adjustable time step size. Numerical examples are presented to demonstrate the efficiency of the computational algorithm in two dimensions. Both periodic Dirichlet and homogeneous Neumann boundary conditions are applied. Also constant and degenerate mobilities are considered. We demonstrate that the use of C-1-continuous natural element shape functions enables the computation of topologically correct solutions on arbitrarily shaped domains.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available