We develop efficient moving mesh algorithms for one-and two-dimensional hyperbolic systems of conservation laws.The algorithms are formed by two independent parts: PDE evolution and mesh-redistribution.The first part can be any appropriate high-resolution scheme, and the second part is based on an iterative procedure.In each iteration, meshes are first redistributed by an quidistribution principle, and then on the resulting new grids the underlying numerical solutions are updated by a conservative-interpolation formula proposed in this work.The iteration for the meshredistribution at a given time step is complete when the meshes governed by a nonlinear equation reach the equilibrium state.The main idea of the proposed method is to keep the mass-conservation of the underlying numerical solution at each redistribution step.In one dimension, we can show that the underlying numerical approximation obtained in the mesh-redistribution part satisfies the desired TVD property, which guarantees that the numerical solution at any time level is TVD, provided that the PDE solver in the first part satisfies such a property.Sev eral test problems in one and two dimensions are computed using the proposed moving mesh algorithm.The computations demonstrate that our methods are efficient for solving problems with shock discontinuities, obtaining the same resolution with a much smaller number of grid points than the uniform mesh approach.

It was generally expected that monotone schemes are oscillation-free for hyperbolic conservation laws. However, recently local oscillations were observed and usually understood to be caused by relative phase errors. In order to further explain this, we first investigate the discretization of initial data that trigger the chequerboard mode, the highest frequency mode. Then we proceed to use the discrete Fourier analysis and the modified equation analysis to distinguish the dissipative and dispersive effects of numerical schemes for low frequency and high frequency modes, respectively. It is shown that the relative phase error is of order O(1) for the high frequency modes unj=λnkeiξj, ξ≈π, but of order O(ξ2) for low frequency modes (ξ≈0). In order to avoid numerical oscillations, the relative phase errors should be offset by numerical dissipation of at least the same order. Numerical damping, i.e. the zero order term in the corresponding modified equation, is important to dissipate the oscillations caused by the relative phase errors of high frequency modes. This is in contrast to the role of numerical viscosity, the second order term, which is the lowest order term usually present to suppress the relative phase errors of low frequency modes.

This paper presents an adaptive moving mesh algorithm for two-dimensional (2D) ideal magnetohydrodynamics (MHD) that utilizes a staggered constrained transport technique to keep the magnetic field divergence-free. The algorithm consists of two independent parts: MHD evolution and mesh-redistribution. The first part is a high-resolution, divergencefree, shock-capturing scheme on a fixed quadrangular mesh, while the second part is an iterative procedure. In each iteration, mesh points are first redistributed, and then a conservative-interpolation formula is used to calculate the remapped cell-averages of the mass, momentum, and total energy on the resulting new mesh; the magnetic potential is remapped to the new mesh in a non-conservative way and is reconstructed to give a divergence-free magnetic field on the new mesh. Several numerical examples are given to demonstrate that the proposed method can achieve high numerical accuracy, track and resolve strong shock waves in ideal MHD problems, and preserve divergence-free property of the magnetic field. Numerical examples include the smooth Alfve´n wave problem, 2D and 2.5D shock tube problems, two rotor problems, the stringent blast problem, and the cloud-shock interaction problem.

This paper presents a second-order accurate adaptive generalized Riemann problem (GRP) scheme for one and two dimensional compressible fluid flows. The current scheme consists of two independent parts: Mesh redistribution and PDE evolution. The first part is an iterative procedure. In each iteration, mesh points are first redistributed, and then a conservative-interpolation formula is used to calculate the cell-averages and the slopes of conservative variables on the resulting new mesh. The second part is to evolve the compressible fluid flows on a fixed nonuniform mesh with the Eulerian GRP scheme, which is directly extended to two dimensional arbitrary quadrilateral meshes. Several numerical examples show that the current adaptive GRP scheme does not only improve the resolution as well as accuracy of numerical solutions with a few mesh points, but also reduces possible errors or oscillations effectively.