This paper develops an adaptive moving mesh method to solve a phase field model for the mixture of two incompressible fluids. The projection method is implemented on a half-staggered, moving quadrilateral mesh to keep the velocity field divergence-free, and the conjugate gradient or multigrid method is employed to solve the discrete Poisson equations. The current algorithm is composed by two independent parts: evolution of the governing equations and mesh-redistribution. In the first part, the incompressible Navier-Stokes equations are solved on a fixed half-staggered mesh by the rotational incremental pressure-correction scheme, and the Allen-Cahn type of phase equation is approximated by a conservative, second-order accurate central difference scheme, where the Lagrangian multiplier is used to preserve the mass-conservation of the phase field. The second part is an iteration procedure. During the mesh redistribution, the phase field is remapped onto the newly resulting meshes by the high-resolution conservative interpolation, while the non-conservative interpolation algorithm is applied to the velocity field. The projection technique is used to obtain a divergence-free velocity field at the end of this part. The resultant numerical scheme is stable, mass conservative, highly efficient and fast, and capable of handling variable density and viscosity. Several numerical experiments are presented to demonstrate the efficiency and robustness of the proposed algorithm.

This paper presents an efficient adaptive mesh redistribution method to solve a non-linear Dirac (NLD) equation. Our algorithm is formed by three parts: the NLD evolution, the iterative mesh redistribution of the coarse mesh and the local uniform refinement of the final coarse mesh. At each time level, the equidistribution principle is first employed to iteratively redistribute coarse mesh points, and the scalar monitor function is subsequently interpolated on the coarse mesh in order to do one new iteration and improve the grid adaptivity. After an adaptive coarse mesh is generated ideally and finally, each coarse mesh interval is equally divided into some fine cells to give an adaptive fine mesh of the physical domain, and then the solution vector is remapped on the resulting new fine mesh by an affine method. The NLD equation is finally solved by using a high resolution shock-capturing method on the (fixed) non-uniform fine mesh. Extensive numerical experiments demonstrate that the proposed adaptive mesh method gives the third-order rate of convergence, and yields an efficient and fast NLD solver that tracks and resolves both small, local and large solution gradients automatically.

This note gives a numerical investigation for the popular high resolution conservative schemes when applied to inviscid, compressible, perfect gas flows with an initial high density ratio as well as a high pressure ratio. The results show that they work very inefficiently and may give inaccurate numerical results even over a very fine mesh when applied to such a problem. Numerical tests show that increasing the order of accuracy of the numerical schemes does not help much in improving the numerical results. How to cure this difficulty is still open.

This Letter presents a numerical study of the interaction dynamics for the solitary waves of a nonlinear Dirac field with scalar self-interaction by using a fourth order accurate Runge-Kutta discontinuous Galerkin (RKDG) method. Some new interaction phenomena are observed: (a) a new quasi-stable long-lived oscillating bound state from the binary collisions of a single-humped soliton and a two-humped soliton; (b) collapse in binary and ternary collisions; (c) strongly inelastic interaction in ternary collisions; and (d) bound states with a short or long lifetime from ternary collisions.

Based on a simple projection of the solution increments of the underlying partial differential equations (PDEs) at each local time level, this paper presents a difference scheme for nonlinear Hamilton-Jacobi (H-J) equations with varying time and space grids. The scheme is of good consistency and monotone under a local CFL-type condition. Moreover, one may deduce a conservative local time step scheme similar to Osher and Sanders scheme approximating hyperbolic conservation law (CL) from our scheme according to the close relation between CLs and H-J equations. Second order accurate schemes are constructed by combining the reconstruction technique with a second order accurate Runge-Kutta time discretization scheme or a Lax-Wendroff type method. They keep some good properties of the global time step schemes, including stability and convergence, and can be applied to solve numerically the initial-boundary-value problems of viscous H-J equations. They are also suitable to parallel computing. Numerical errors and the experimental rate of convergence in the Lp-norm, p=1, 2, and ∞, are obtained for several one-and two-dimensional problems. The results show that the present schemes are of higher order accuracy.