“…In Sections 5 and 6 we report the numerical results for the cases with two-sided and one-sided degenerate diffusion mobility, respectively, with emphasis on the coarsening mechanisms and rates. Even though there have been numerous computational simulations for CH equations in the literature [3,4,5,7,8,14,16,17,18,19,20,21,23,24,25,27,32,39,40,42,43,44,46,50,51,54,55] and the references cited therein, the results given in this work represent the most comprehensive so far: we take small time steps in our simulations not only to ensure stability but most importantly to ascertain time accuracy; we run the simulations for exceedingly long time intervals so that the solution behaviors reach well into appropriate asymptotic regimes; we conduct numerical studies of a test case involving only an elliptical domain and another test case involving only two spherical domains to demonstrate the consistency of numerical solutions with the physical Gibbs-Thomson effect and to illustrate the communication between disjoint domains and the mechanism behind the coarsening; we also carefully perform large scale simulations on systems involving convoluted spatial patterns to demonstrate the dependence of coarsening rates on the diffusion mobilities and volume fractions. For each case of two-sided and one-sided degenerate mobilities, we analyze three typical scenarios of 75%, 50% and 25% volume fractions corresponding to having large, equal and small positive phases and compare their coarsening rates.…”