Efficient management of parallelism in object-oriented numerical software libraries
Modern software tools for scientific computing
Finite element approximation of the Cahn-Hilliard equation with concentration dependent mobility
Mathematics of Computation
Finite-difference schemes for nonlinear wave equation that inherit energy conservation property
Journal of Computational and Applied Mathematics
A discontinuous Galerkin method for the Cahn-Hilliard equation
Journal of Computational Physics
libMesh: a C++ library for parallel adaptive mesh refinement/coarsening simulations
Engineering with Computers
A nonstiff, adaptive mesh refinement-based method for the Cahn-Hilliard equation
Journal of Computational Physics
Local discontinuous Galerkin methods for the Cahn-Hilliard type equations
Journal of Computational Physics
A time-adaptive finite volume method for the Cahn-Hilliard and Kuramoto-Sivashinsky equations
Journal of Computational Physics
An efficient moving mesh spectral method for the phase-field model of two-phase flows
Journal of Computational Physics
Hi-index | 31.45 |
The Cahn-Hilliard (CH) equation is a time-dependent fourth-order partial differential equation (PDE). When solving the CH equation via the finite element method (FEM), the domain is discretized by C^1-continuous basis functions or the equation is split into a pair of second-order PDEs, and discretized via C^0-continuous basis functions. In the current work, a quantitative comparison between C^1 Hermite and C^0 Lagrange elements is carried out using a continuous Galerkin FEM formulation. The different discretizations are evaluated using the method of manufactured solutions solved with Newton's method and Jacobian-Free Newton Krylov. It is found that the use of linear Lagrange elements provides the fastest computation time for a given number of elements, while the use of cubic Hermite elements provides the lowest error. The results offer a set of benchmarks to consider when choosing basis functions to solve the CH equation. In addition, an example of microstructure evolution demonstrates the different types of elements for a traditional phase-field model.