In this paper, a new discontinuous Galerkin method is developed for the parabolic equation with jump coefficients satisfying the continuous flow condition. Theoretical analysis shows that this method is L2 stable. When the finite element space consists of interpolative polynomials of degrees k, the convergent rate of the semi-discrete discontinuous Galerkin scheme has an order of . Numerical examples for both 1-dimensional and 2-dimensional problems demonstrate the validity of the new method.

The inverse problem considered in this paper is to determine the shape and the impedance of crack from a knowledge of the time-harmonic incident field and the corresponding far field pattern of the scattered waves in two-dimension. The combined single- and double-layer potential is used to approach the scattered waves. As an important feature, this method does not require the solution of u and ∂u/∂v at each iteration. An approximate method is presented and the convergence of this method is proven. Numerical examples are given to show that this method is both accurate and simple to use.

Two algorithms for constructing a class of compactly supported conjugate symmetric complex tight wavelet frames ψ = {ψ1, ψ2} are derived. Firstly, a necessary and sufficient condition for constructing the conjugate symmetric complex tight wavelet frames is established. Secondly, based on a given conjugate symmetric low pass filter, a description of a family of complex wavelet frame solutions is provided when the low pass filter is of even length. When one wavelet is conjugate symmetric and the other is conjugate antisymmetric, the two wavelet filters can be obtained by matching the roots of associated polynomials. Finally, two examples are given to illustrate how to use our method to construct conjugate symmetric complex tight wavelet frames which have some vanishing moments.

The numerical solution of the harmonic heat map flow problems with blowup in finite or infinite time is considered using an adaptive moving mesh method. A properly chosen monitor function is derived so that the moving mesh method can be used to simulate blowup and produce accurate blowup profiles which agree with formal asymptotic analysis. Moreover, the moving mesh method has finite time blowup when the underlying continuous problem does. In situations where the continuous problem has infinite time blowup, the moving mesh method exhibits finite time blowup with a blowup time tending to infinity as the number of mesh points increases. The inadequacy of a uniform mesh solution is clearly demonstrated.

In this paper, we propose a robust finite volume scheme to numerically solve the shallow water equations on complex rough topography. The major difficulty of this problem is introduced by the stiff friction force term and the wet/dry interface tracking. An analytical integration method is presented for the friction force term to remove the stiffness. In the vicinity of wet/dry interface, the numerical stability can be attained by introducing an empirical parameter, the water depth tolerance, as extensively adopted in literatures. We propose a problem independent formulation for this parameter, which provides a stable scheme and preserves the overall truncation error of . The method is applied to solve problems with complex rough topography, coupled with h-adaptive mesh techniques to demonstrate its robustness and efficiency.

Aiming at the isoparametric bilinear finite volume element scheme, we initially derive an asymptotic expansion and a high accuracy combination formula of the derivatives in the sense of pointwise by employing the energy-embedded method on uniform grids. Furthermore, we prove that the approximate derivatives are convergent of order two. Finally, numerical examples verify the theoretical results.

This paper is concerned with obtaining an approximate solution and an approximate derivative of the solution for neutral Volterra integro-differential equation with a weakly singular kernel. The solution of this equation, even for analytic data, is not smooth on the entire interval of integration. The Jacobi collocation discretization is proposed for the given equation. A rigorous analysis of error bound is also provided which theoretically justifies that both the error of approximate solution and the error of approximate derivative of the solution decay exponentially in L∞ norm and weighted L2 norm. Numerical results are presented to demonstrate the effectiveness of the spectral method.