In this paper, an arbitrary Lagrangian–Eulerian local discontinuous Galerkin (ALE-LDG) method for Hamilton–Jacobi equations will be developed, analyzed and numerically tested. This method is based on the time-dependent approximation space defined on the moving mesh. A priori error estimates will be stated with respect to the L∞ 􏰒(0, T ; L2)norm. In particular, the optimal (k + 1) convergence in one dimension and the suboptimal (k + 1 ) convergence in two dimensions will be proven for the semi-discrete method, when a local Lax–Friedrichs flux and piecewise polynomials of degree k on the reference cell are used. Furthermore, the validity of the geometric conservation law will be proven for the fully-discrete method. Also, the link between the piecewise constant ALE-LDG method and the monotone scheme, which converges to the unique viscosity solution, will be shown. The capability of the method will be demonstrated by a variety of one and two dimensional numerical examples with convex and noneconvex Hamiltonian.

We present a uniformly first order accurate numerical method for solving the Klein-Gordon-Zakharov (KGZ) system with two dimensionless parameters 0 < ε ≤ 1 and 0 <γ ≤ 1, which are inversely proportional to the plasma frequency and the acoustic speed, respectively. In the simultaneous high-plasma-frequency and subsonic limit regime, i.e. ε < γ → 0+, the KGZ system collapses to a cubic Schrödinger equation, and the solution propagates waves with O ( ε2)-wavelength in time and meanwhile contains rapid outgoing initial layers with speed O (1/ γ ) in space due to the incompatibility of the initial data. By presenting a multiscale decomposition of the KGZ system, we propose a multiscale time integrator Fourier pseudospectral method which is explicit, efficient and uniformly accurate for solving the KGZ system for all 0 < ε < γ ≤ 1. Numerical results are reported to show the efficiency and accuracy of scheme. Finally, the method is applied to investigate the convergence rates of the KGZ system to its limiting models when ε < γ → 0+.

In this paper, we introduce a general approach for proving optimal $L^2$ error estimates for the semi-discrete local discontinuous Galerkin (LDG) methods solving linear high order wave equations. The optimal order of error estimates hold not only for the solution itself but also for the auxiliary variables in the LDG method approximating the various order derivatives of the solution. Several examples including the one-dimensional third order wave equation, one-dimensional fifth order wave equation, and multi-dimensional Schr\"{o}dinger equation are explored to demonstrate this approach. The main idea is to derive energy stability for the various auxiliary variables in the LDG discretization, via using the scheme and its time derivatives with different test functions. Special projections are utilized to eliminate the jump terms at the cell boundaries in the error estimate in order to achieve the optimal order of accuracy.

No abstract uploaded!

We present a new formulaiton of the incompressible Navier-Stokes equation in terms of an auxiliary field that differs from the velocity by a gauge transformation. The gauge freedom allows us to assign simple and specific boundary conditions for both the auxiliary field and the gauge field, thus eliminating the issue of pressure boundary condition in the usual primitive variable formulation. The resulting dynamic and kinematic equations can then be solved by standard methods for heat and Poisson equations. A normal mode analysis suggests that in contrast to the classical projection method, the gauge method does not suffer from the problem of numerical boundary layers. Thus the subtleties on the spatial discretization for the projection methods are removed. Consequently, the projection step in the gauge method can be accomplished by standard Poissin solves. We demonstrate the efficiency and accuracy of the gauge method by several numerical examples, including the flow past cylinder.