We introduce two exponential-type integrators for the “good” Boussinesq equation. They are of orders one and two, respectively, and they require lower spatial regularity of the solution compared to classical exponential integrators. For the first integrator, we prove first-order convergence in H^r for solutions in H^{r+1} with r>1/2. This new integrator even converges (with lower order) in H^r for solutions in H^r, i.e., without any additional smoothness assumptions. For the second integrator, we prove second-order convergence in H^r for solutions in H^{r+3} with r>1/2 and convergence in L^2 for solutions in H^3. Numerical results are reported to illustrate the established error estimates. The experiments clearly demonstrate that the new exponential-type integrators are favorable over classical exponential integrators for initial data with low regularity.

The surveys in the field of nonlinear filtering (NLF) are enumerous. Most of them are application-oriented and served as the tutorials for the practioners. The local approaches, including Kalman filter and its invariants, have already been studied from various point of views, due to its off-the-shelf implementation and wide applications. However, it cannot give good estimation of the states in highly nonlinear system or with non-Gaussian initial conditional density functions. Moreover, while the local methods only approximate the mean and variance, the global ones seek the way to directly obtain the conditional density function of the states. Consequently, all the statistical information is acquired. In this survey, we shall briefly go through the local approaches and put emphases on the existing three major global approaches: finite-dimensional NLF, sequential Monte Carlo methods (particle filter) and the {Yau-Yau's} on- and off-line solver of Duncan-Mortensen-Zakai's equation {\cite{YY_08}}. The discussions are mainly from the mathematical point of view.

No abstract uploaded!

No abstract uploaded!

The numerical approximation of high-frequency wave propagation in inhomogeneous media is a challenging problem. In particular, computing high-frequency solutions by direct simulations requires several points per wavelength for stability and usually requires many points per wavelength for a satisfactory accuracy. In this paper, we propose a new method for the acoustic wave equation in inhomogeneous media in the time domain to achieve superior accuracy and stability without using a large number of unknowns. The method is based on a discontinuous Galerkin discretization together with carefully chosen basis functions. To obtain the basis functions, we use the idea from geometrical optics and construct the basis functions by using the leading order term in the asymptotic expansion. Also, we use a wavefront tracking method and a dimension reduction procedure to obtain dominant rays in each cell. We show