• In this paper, we study the phase field models with fractional-order in time. The phase field models have been widely used to study coarsening dynamics of material systems with microstructures. It is known that phase field models are usually derived from energy variation so that they obey some energy dissipation laws intrinsically. Recently, many works have been published on investigating fractional-order phase field models, but little is known of the corresponding energy dissipation laws. We focus on the time-fractional phase field models and report that the effective free energy and roughness obey a universal power-law scaling dynamics during coarsening. Mainly, the effective free energy and roughness in the time-fractional phase field models scale by following a similar power law as the integer phase field models, where the power is linearly proportional to the fractional order. This universal scaling law is verified numerically against several phase field models, including the Cahn-Hilliard equations with different variable mobilities and molecular beam epitaxy models. This new finding sheds light on potential applications of time fractional phase field models in studying coarsening dynamics and crystal growths.
  • In this paper, we propose a time-fractional molecular beam epitaxy (MBE) model with slope selection and its efficient, accurate, full discrete, linear numerical approximation. The numerical scheme utilizes the fast algorithm for the Caputo fractional derivative operator in time discretization and Fourier spectral method in spatial discretization. Refinement tests are conducted to verify the $2-\alpha$ order of time convergence, with $\alpha \in (0, 1]$ the fractional order of derivative. Several numerical simulations are presented to demonstrate the accuracy and efficiency of our newly proposed scheme. By exploring the fast algorithm calculating the Caputo fractional derivative, our numerical scheme makes it practice for long time simulation of MBE coarsening, which is essential for MBE model in practice. With the proposed fractional MBE model, we observe that the scaling law for the energy decays as $ O(t^{-\frac{\alpha}{3}})$ and the roughness increases as $O(t^{\frac{\alpha}{3}})$, during the coarsening dynamics with random initial condition. That is to say, the coarsening rate of MBE model could be manipulated by the fractional order $\alpha$, and it is linearly proportional to $\alpha$. This is the first time in literature to report/discover such scaling correlation. It provides a potential application field for fractional differential equations. Besides, the numerical approximation strategy proposed in this paper can be readily applied to study many classes of time-fractional and high dimensional phase field models.
  • In this paper, we consider the numerical approximations for the fourth order Cahn-Hilliard equation with concentration dependent mobility, and the logarithmic Flory-Huggins potential. One challenge in solving such a diffusive system numerically is how to develop proper temporal discretization for nonlinear terms in order to preserve the energy stability at the time-discrete level. We resolve this issue by developing a set of the first and second order time marching schemes based on a novel, called "Invariant Energy Quadratization" approach. Its novelty is that the proposed scheme is linear and symmetric positive definite because all nonlinear terms are treated semi-explicitly. We further prove all proposed schemes are unconditionally energy stable rigorously. Various 2D and 3D numerical simulations are presented to demonstrate the stability, accuracy and efficiency of the proposed schemes thereafter.
  • How to develop efficient numerical schemes while preserving the energy stability at the discrete level is a challenging issue for the three component Cahn-Hilliard phase-field model. In this paper, we develop first and second order temporal approximation schemes based on the "Invariant Energy Quadratization" approach, where all nonlinear terms are treated semi-explicitly. Consequently, the resulting numerical schemes lead to a well-posed linear system with the symmetric positive definite operator to be solved at each time step. We rigorously prove that the proposed schemes are unconditionally energy stable. Various 2D and 3D numerical simulations are presented to demonstrate the stability and the accuracy of the schemes.
  • In this paper, we consider the numerical approximations for the fourth order viscous Cahn-Hilliard equation with the hyperbolic relaxation. The main challenge in solving such a diffusive system numerically is how to develop high order temporal discretization for the hyperbolic and nonlinear terms that allows large time step while preserving the unconditional energy stability, i.e., the energy dissipative structure at the time-discrete level. We resolve this issue by developing two second order time marching schemes using the recently developed "Invariant Energy Quadratization" approach where all nonlinear terms are treated semi-explicitly. In each time step, one only needs to solve a linear system that is symmetric positive definite. We rigorously prove all proposed schemes are unconditionally energy stable, which is further verified by time step refinement test numerically. Moreover, a number of 2D and 3D numerical simulations are presented to demonstrate the stability, accuracy and efficiency of the proposed schemes.
  • Let $G=(V,E)$ and $H$ be two graphs. Packing problem is to find in $G$ the largest number of independent subgraphs each of which is isomorphic to $H$. Let $U\subset{V}$. If the graph $G-U$ has no subgraph isomorphic to $H$, $U$ is a cover of $G$. Covering problem is to find the smallest set $U$. The vertex-disjoint tree packing was not sufficiently discussed in literature but has its applications in data encryption and in communication networks such as multi-cast routing protocol design. In this paper, we give the kind of $(k+1)$-connected graph $G'$ into which we can pack independently the subgraphs that are each isomorphic to the $(2^{k+1}-1)$-order perfect binary tree $T_k$. We prove that in $G'$ the largest number of vertex-disjoint subgraphs isomorphic to $T_k$ is equal to the smallest number of vertices that cover all subgraphs isomorphic to $T_k$. Then, we propose that $T_k$ does not have the \emph{Erd\H{o}s-P\'{o}sa} property. We also prove that the $T_k$ packing problem in an arbitrary graph is NP-hard, and propose the distributed approximation algorithms.
  • We say that a cyclotomic polynomial \Phi_{n}(x) has order three if n is the product of three distinct primes, p<q<r. Let A(n) be the largest absolute value of a coefficient of \Phi_{n}(x) and M(p) be the maximum of A(pqr). In 1968, Sister Marion Beiter conjectured that A(pqr)<=(p+1)/2. In 2008, Yves Gallot and Pieter Moree showed that the conjecture is false for every p>=11, and they proposed the Corrected Beiter conjecture: A(pqr)<=2p/3. Here we will give a proof of this conjecture.
  • We say that a cyclotomic polynomial Phi_{n}(x) has order three if n is the product of three distinct primes, p<q<r. Let A(n) be the largest absolute value of a coefficient of Phi_{n}(x) and M(p) be the maximum of A(pqr). In 1968, Sister Marion Beiter conjectured that A(pqr)<=(p+1)/2. In 2008, Yves Gallot and Pieter Moree showed that the conjecture is false for every p>=11, and they proposed the Corrected Beiter conjecture: M(p)<=2p/3. Here we will give a sufficient condition for the Corrected Beiter conjecture and prove it when p=7.