LAUSR

Abstract A 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 finite-volume 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.