In this paper, we propose a numerical method to solve isotropic elliptic equations on point cloud by generalizing the point integral method. The idea of the point integral method is to approximate the differential operators by integral operators and discretize the corresponding integral equation on point cloud. The key step is to get the integral approximation. In this paper, with proper kernel function, we get an integral approximation for the elliptic operators with isotropic coefficients. Moreover, the integral approximation has been proved to keep the coercivity of the original elliptic operator. The convergence of the point integral method is also proved.

In this paper, we generalize the so-called inverse Lax-Wendroff boundary treatment \cite{ilw} for the inflow boundary of a linear hyperbolic problem discretized by the recently introduced central compact schemes. The outflow boundary is treated by the classical extrapolation and a stability analysis for the resulting scheme is provided. To ensure the stability of the considered schemes provided with the chosen boundaries, the G-K-S theory is used, first in the semidiscrete case then in the fully discrete case with the third-order TVD Runge-Kutta time discretization. Afterwards, due to the high algebraic complexity of the G-K-S theory, the stability is analyzed by visualizing the eigenspectrum of the discretized operators. We show in this paper that the results obtained with these two different approaches are perfectly consistent. We also illustrate the high accuracy of the presented schemes on simple test problems.

This is the second of a series of papers on the subject of projection methods for viscous incompressible flow calculations. The purpose of the present paper is to explain why the accuracy of the velocity approximation is not affected by (1) the numerical boundary layers in the approximation of pressure and the intermediate velocity field and (2) the noncommutativity of the projection operator and the laplacian. This is done by using a Godunov–Ryabenki type of analysis in a rigorous fashion. By doing so, we hope to be able to convey the message that normal mode analysis is basically sufficient for understanding the stability and accuracy of a finite-difference method for the Navier–Stokes equation even in the presence of boundaries. As an example, we analyze the second-order projection method based on pressure increment formulations used by van Kan and Bell, Colella, and Glaz. The leading order error term in this case is of $O(\Delta t)$ and behaves as high frequency oscillations over the whole domain, compared with the $O(\Delta t^{{1 / 2}} )$ numerical boundary layers found in the second-order Kim–Momn method.

Computing stationary states is an important topic for phase field crystal (PFC) models. Great efforts have been made for energy dissipation of the numerical schemes when using gradient flows. However, it is always time-consuming due to the requirement of small effective time steps. In this paper, we propose an adaptive accelerated proximal gradient method for finding the stationary states of PFC models. The energy dissipation is guaranteed and the convergence property is established for the discretized energy functional. Moreover, the connections between generalized proximal operator with classical (semi-)implicit and explicit schemes for gradient flow are given. Extensive numerical experiments, including two three dimensional periodic crystals in Landau-Brazovskii (LB) model and a two dimensional quasicrystal in Lifshitz-Petrich (LP) model, demonstrate that our approach has adaptive time steps which lead to significant acceleration over semi-implicit methods for computing complex structures. Furthermore, our result reveals a deep physical mechanism of the simple LB model via which the sigma phase is first discovered.

A Deuflhard-type exponential integrator sine pseudospectral (DEI-SP) method is proposed and analyzed for solving the generalized improved Boussinesq (GIBq) equation. The numerical scheme is based on a second-order exponential integrator for time integration and a sine pseudospectral discretization in space. Rigorous analysis and abundant experiments show that the method converges quadratically and spectrally in time and space, respectively. Finally the DEI-SP method is applied to investigate the complicated and interesting long-time dynamics of the GIBq equation.