“…As a well-known approach, the convex splitting approach [8] can guarantee the numerical energy stability. For solving the PFC model, the convex splitting approach combined with different spatial discretization methods is studied, including spectral method [4,19], finite difference method [13,16,29], finite element method [9,26], and local discontinuous Galerkin method [11,12]. However, all second-order schemes in the above study are nonlinear, and their implementations are too complex and costly owing to nonlinear equations to be solved at each time step.…”