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 investigate the Chern–Ricci flow, an evolution equation of Hermitian metrics generalizing the Kähler–Ricci flow, on elliptic bundles over a Riemann surface of genus greater than one. We show that, starting at any Gauduchon metric, the flow collapses the elliptic fibers and the metrics converge to the pullback of a Kähler–Einstein metric from the base. Some of our estimates are new even for the Kähler–Ricci flow. A consequence of our result is that, on every minimal non-Kähler surface of Kodaira dimension one, the Chern–Ricci flow converges in the sense of Gromov–Hausdorff to an orbifold Kähler–Einstein metric on a Riemann surface.

We study a robust and efficient eigensolver for computing the positive dense spectrum of the two-dimensional transmission eigenvalue problem (TEP) which is derived from the Maxwell’s equation with complex media in pseudo-chiral model and the transverse magnetic mode. The discretized governing equations by the N ́ed ́elec edge element result in a large-scale quadratic eigenvalue problem (QEP). We estimate half of the positive eigenvalues of the QEP are on some interval which forms a dense spectrum of the QEP. The quadratic Jacobi-Davidson method with a so-called non-equivalence deflation technique is proposed to compute the dense spectrum of the QEP. Intensive numerical experiments show that our proposed method makes the convergence efficiently and robustly even it needs to compute more than 5000 desired eigenpairs. Numerical results also illustrate that the computed eigenvalue curves can be approximated by the non- linear functions which can be applied to estimate the density of the eigenvalues for the TEP.

No abstract uploaded!

No abstract uploaded!