Zhangpeng Sun

Tiao Lu, Zhangpeng Sun

For the stationary Wigner equation with inflow boundary conditions, its numerical convergence with respect to the velocity mesh size are deteriorated due to the singularity at velocity zero. In this paper, using the fact that the solution of the stationary Wigner equation is subject to an algebraic constraint, we prove that the Wigner equation can be written into a form with a bounded operator $\mathcal{B}[V]$, which is equivalent to the operator $\mathcal{A}[V]=Θ[V]/v$ in the original Wigner equation under some conditions. Then the discrete operators discretizing $\mathcal{B}[V]$ are proved to be uniformly bounded with respect to the mesh size. Based on the therectical findings, a signularity-free numerical method is proposed. Numerical reuslts are proivded to show our improved numerical scheme performs much better in numerical convergence than the original scheme based on discretizing $\mathcal{A}[V]$.

Weifeng Hou, Zhangpeng Sun, Wenqi Yao et al.

The compact finite difference method is a powerful tool for discretizing conservation laws, owing to its inherent flexibility in developing high-resolution and highly stable schemes. In this paper, we propose a framework for the design of genuine globally conservative compact finite difference schemes, which addresses a critical requirement in conservation laws. Within our framework, we rigorously establish that the discrete conservation law maintains strict conservation for flux functions in polynomial spaces with optimal algebraic order, i.e., the discrete scheme achieves an optimal algebraic precision.Our work advances the existing conservative compact finite difference schemes, which rely on approaches to maintaining global conservation that are fundamentally consistent with the method proposed by Lele [Lele, J. Comput. Phys., 1992]. As an application, we propose an algorithm for designing globally conservative fourth-order schemes, aimed at optimizing resolution and asymptotic stability. Three schemes are generated using the algorithm, with their excellent performance across multiple aspects validated through numerical experiments.