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.

No abstract uploaded!

In 1962 Sprind\v{z}uk proved Mahler's conjecture in both the real and complex cases. Baker gave a generalized result by using a modified version of Sprind\v{z}uk's method. Later, Baker and Schmidt derived the Hausdorff dimensions of sets which are defiend in terms of approximation by algebraic numbers of bounded degrees by using Baker's theorem. In this article we will prove two analogue theorems in the fields of formal power series over finite finite fields.

The study of the present paper is inspired by Gangl, Kaneko and Zagier's result of the connection with double zeta values and modular forms. We introduce double Eisenstein series $E_{r,s}$ in positive characteristic with double zeta values $\zeta_A(r,s)$ as their constant term and compute the t-expansions of the double Eisenstein series. Moreover, we derive the shuffle relations of double Eisenstein series which match the shuffle relations of double zeta values in \cite{Huei2015}.

Multizeta values in positive characteristic were first introduced and studied by Thakur. He and Lara Rodr\'{\i}guez \cite{LRT} discovered and conjectured certain zeta-like families. Kuan, Lin and Yu stated more conjectures about zeta-like multizeta values in \cite{LK}. In the present paper we study and give the transparent formula for certain Anderson-Thakur polynomials. This enables us to confirm the conjectured zeta-like families.

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.

The main purpose of this paper is to analyze the stability and error estimates of the local discontinuous Galerkin (LDG) methods coupled with implicit-explicit (IMEX) time discretization schemes, for solving one-dimensional convection-diffusion equations with a nonlinear convection. Both Runge-Kutta and multi-step IMEX methods are considered. By the aid of the energy method, we show that the IMEX LDG schemes are unconditionally stable for the nonlinear problems, in the sense that the time-step $\dt$ is only required to be upper-bounded by a positive constant which depends on the flow velocity and the diffusion coefficient, but is independent of the mesh size $h$. We also give optimal error estimates for the IMEX LDG schemes, under the same temporal condition, if a monotone numerical flux is adopted for the convection. Numerical experiments are given to verify our main results.