Journal
COMPUTATIONAL MATERIALS SCIENCE
Volume 140, Issue -, Pages 22-31Publisher
ELSEVIER SCIENCE BV
DOI: 10.1016/j.commatsci.2017.08.025
Keywords
Immersed interface method; Cahn-Hilliard equation; Surface wetting; Fourier integration
Categories
Ask authors/readers for more resources
fast immersed boundary method for the Cahn-Hilliard equation is introduced. The decomposition of the fourth-order non-linear Cahn-Hilliard equation into a system of linear parabolic second-order equations allows to pose arbitrary Neumann or surface wetting conditions on the boundary. In space a finitevolume discretization on a regular Cartesian voxel grid allows the use of fast parabolic solvers via Fourier transform of arbitrary convergence order. For the time discretization, a second-order Runge-Kutta scheme is applied. The polynomial approximation of the chemical potential results in a numerical scheme that is unconditionally gradient-stable and allows large time steps. With an additional pre-conditioner for the linear system, the condition of linear system is minimized. By this the convergence is independent of both spatial discretization and time step size. This allows for the simulation of phase separation in large porous complex domains with three dimensions and over hundred million degrees of freedom while applying arbitrary boundary conditions and using large time steps. (C) 2017 Elsevier B.V. All rights reserved.
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