Fourier spectral methods achieve exponential accuracy both on the approximation level and for solving partial differential equations (PDEs), if the solution is analytic. If the solution is discontinuous but piecewise analytic up to the discontinuities, Fourier spectral methods produce poor pointwise accuracy, but still maintains exponential accuracy after post-processing. In earlier work, an extended technique is provided to recover exponential accuracy for functions which have end-point singularities, from the knowledge of point values on standard collocation points. In this paper, we develop a technique to recover exponential accuracy from the first $N$ Fourier coefficients of functions which are analytic in the open interval but have unbounded derivative singularities at end points. With this post-processing method, we are able to obtain exponential accuracy of spectral methods applied to linear transport equations involving such functions.

This paper is concerned with superconvergence properties of discontinuous Galerkin (DG) methods for 2-D linear hyperbolic conservation laws over rectangular meshes when upwind fluxes are used. We prove, under some suitable initial and boundary discretizations, the ($2k+1$)-th order superconvergence rate of the DG approximation at the downwind points and for the cell averages, when piecewise tensor-product polynomials of degree $k$ are used. Moreover, we prove that the gradient of the DG solution is superconvergent with a rate of ($k+1$)-th order at all interior left Radau points; and the function value approximation is superconvergent at all right Radau points with a rate of ($k+2$)-th order. Numerical experiments indicate that the aforementioned superconvergence rates are sharp.

For numerical simulation of detonation, computational cost using uniform meshes is large due to the vast separation in both time and space scales. Adaptive mesh refinement (AMR) is advantageous for problems with vastly different scales. This paper aims to propose an AMR method with high order accuracy for numerical investigation of multi-dimensional detonation. A well-designed AMR method based on finite difference weighted essentially non-oscillatory (WENO) scheme, named as AMR\&WENO is proposed. A new cell-based data structure is used to organize the adaptive meshes. The new data structure makes it possible for cells to communicate with each other quickly and easily. In order to develop an AMR method with high order accuracy, high order prolongations in both space and time are utilized in the data prolongation procedure. Based on the message passing interface (MPI) platform, we have developed a workload balancing parallel AMR\&WENO code using the Hilbert space-filling curve algorithm. Our numerical experiments with detonation simulations indicate that the AMR\&WENO is accurate and have a high resolution. Moreover, we evaluate and compare the performance between the uniform mesh WENO scheme and the parallel AMR\&WENO method. The comparison results provide us further insight into the high performance of the parallel AMR\&WENO method.

Hoogendoorn and Bovy (Transportation Research Part B, 2004, 38(7), 571--592) developed an approach for a pedestrian user-optimal dynamic assignment in continuous time and space. Although their model was proposed for pedestrian traffic, it can also be applied to urban cities. The model is very general, and consists of a conservation law (CL) and a Hamilton-Jacobi-Bellman (HJB) equation that contains a minimum value problem. However, only an isotropic application example was given in their paper. We claim that the HJB equation is difficult to compute numerically in an anisotropic case. To overcome this, we reformulate their model for a dense urban city that is arbitrary in shape and has a single central business district (CBD). In our model, the minimum value problem is only used in the CL portion, and the HJB equation reduces to a Hamilton-Jacobi (HJ) equation for easier computation. The dynamic path equilibrium of our model is proven in a different way from theirs, and a numerical algorithm is also provided to solve the model. Finally, we show a numerical example under the anisotropic case and compare the results with those of the isotropic case.

We analyze discontinuous Galerkin methods using upwind-biased numerical fluxes for time-dependent linear conservation laws. In one dimension, optimal a priori error estimates of order $k+1$ are obtained for the semidiscrete scheme when piecewise polynomials of degree at most $k$ $(k \ge 0)$ are used. Our analysis is valid for arbitrary nonuniform regular meshes and for both periodic boundary conditions and for initial-boundary value problems. We extend the analysis to the multidimensional case on Cartesian meshes when piecewise tensor product polynomials are used, and to the fully discrete scheme with explicit Runge--Kutta time discretization. Numerical experiments are shown to demonstrate the theoretical results.