An important property for finite difference schemes designed on curvilinear meshes is the exact preservation of free-stream solutions. This property is difficult to fulfill for high order conservative essentially non-oscillatory (WENO) finite difference schemes. In this paper we explore an alternative flux formulation for such finite difference schemes which can preserve free-stream solutions, based on the numerical technique for the metric terms, which can be applied to this alternative flux formulation but is difficult to be applied to the standard finite difference formulation. Free-stream and vortex preservation properties are investigated, and comparison with standard finite difference WENO schemes is made. Theoretical derivation and numerical results show that the finite difference WENO schemes based on the alternative flux formulation can preserve free-stream and vortex solutions on both stationary and dynamically generalized coordinate systems, hence giving much better performance than the standard finite difference WENO schemes for such problems.

We study a robust and efficient eigensolver for computing a few smallest positive eigenvalues of the three-dimensional Maxwell’s transmission eigenvalue problem. The discretized governing equations by the Nedelec edge element result in a large-scale quadratic eigenvalue problem (QEP) for which the spectrum contains many zero eigenvalues and the coefficient matrices consist of patterns in the matrix form XY^{−1}Z, both of which prevent existing eigenvalue solvers from being efficient. To remedy these difficulties, we rewrite the QEP as a particular nonlinear eigenvalue problem and develop a secant-type iteration, together with an indefinite locally optimal block preconditioned conjugate gradient method (LOBPCG), to sequentially compute the desired positive eigenvalues. Furthermore, we propose a novel method to solve the linear systems in each iteration of LOBPCG. Intensive numerical experiments show that our proposed method is robust, although the desired real eigenvalues are surrounded by complex eigenvalues.

In this paper, we develop an efficient and energy stable fully-discrete local discontinuous Galerkin (LDG) method for the Cahn–Hilliard–Hele–Shaw (CHHS) system. The semi-discrete energy stability of the LDG method is proved firstly. Due to the strict time step restriction ( t = O( x4)) of explicit time discretization methods for stability, we introduce a semi- implicit time integration scheme which is based on a convex splitting of the discrete Cahn– Hilliard energy. The unconditional energy stability has been proved for this fully-discrete LDG scheme. The fully-discrete equations at the implicit time level are nonlinear. Thus, the nonlinear Full Approximation Scheme (FAS) multigrid method has been applied to solve this system of algebraic equations, which has been shown the nearly optimal complexity numerically. Numerical results are also given to illustrate the accuracy and capability of the LDG method coupled with the multigrid solver.

No abstract uploaded!

No abstract uploaded!