We consider the numerical solution of the time-fractional diffusion-wave equation on two-dimensional and three-dimensional unbounded spatial domains. Introduce an artificial boundary and find the exact and approximate artificial boundary conditions for the given problem, which lead to a bounded computational domain. Using the exact or approximating boundary conditions on the artificial boundary, the original problem is reduced to an initial-boundary-value problem on the bounded computational domain which is respectively equivalent to or approximates the original problem. Finite difference methods are used to solve the reduced problems on the bounded computational domain and the stability of these finite difference methods are proved. The numerical results demonstrate that the method given in this paper is effective and feasible.

No abstract uploaded!

This paper is concerned with an interior penalty discontinuous Galerkin (IPDG) method based on a flexible type of non-polynomial local approximation space for the Helmholtz equation with varying wavenumber. The local approximation space consists of multiple polynomial-modulated phase functions which can be chosen according to the phase information of the solution. We obtain some {approximation} properties for this space and \textit{a prior} $L^2$ error estimates for the \textit{h}-convergence of the IPDG method using duality argument. We also provide ample numerical examples to show that, building phase information into the local spaces often gives more accurate results comparing to using the standard polynomial spaces.

In this paper we consider a discontinuous Galerkin discretization of the ideal magnetohydrodynamics (MHD) equations on unstructured meshes, and the divergence free constraint ($\nabla \cdot \B = 0$) of its magnetic field $\B$. We first present two approaches for maintaining the divergence free constraint, namely the approach of a locally divergence free projection inspired by locally divergence free elements \cite{Li2005}, and another approach of the divergence cleaning technique given by Dedner et al. \cite{Dedner2002}. By combining these two approaches we obtain a efficient method at the almost same numerical cost. Finally, numerical experiments are performed to show the capacity and efficiency of the scheme.

A new high order finite-difference method utilizing the idea of Harten ENO subcell resolution method is proposed for chemical reactive flows and combustion. In reaction problems, when the reaction time scale is very small, e.g., orders of magnitude smaller than the fluid dynamics time scales, the governing equations will become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present proposed method is a modified fractional step method which solves the convection step and reaction step separately. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with the computed flow variables in the shock region modified by the Harten subcell resolution idea. For numerical experiments, a fifth-order finite-difference WENO scheme and its anti-diffusion WENO variant are considered. A wide range of 1D and 2D scalar and Euler system test cases are investigated. Studies indicate that for the considered test cases, the new method maintains high order accuracy in space for smooth flows, and for stiff source terms with discontinuities, it can capture the correct propagation speed of discontinuities in very coarse meshes with reasonable CFL numbers.