Kailiang Wu

Linfeng Xu, Shengrong Ding, Kailiang Wu

High-order accurate simulations of special relativistic hydrodynamics (RHD) are prone to numerical breakdown if intrinsic physical constraints (positive rest-mass density/pressure and subluminal velocity) are violated near strong discontinuities. In this work, we develop a robust and efficient physical-constraint-preserving (PCP) flux-limiting framework for high-order schemes, using finite-difference WENO as a representative example. By leveraging the geometric quasilinearization (GQL) representation, which equivalently reformulates the nonlinear RHD constraints into a family of linear inequalities, we integrate a Zalesak-type Flux-Corrected Transport (FCT) update into a scalar-style limiter that acts directly on conservative variables. A critical innovation is the explicit, non-iterative determination of limiting parameters via a rational stereographic parameterization of the GQL normal vector. This technique transforms the required worst-case minimization over auxiliary variables into a generalized Rayleigh-quotient formulation, allowing the optimal parameters to be obtained by solving small symmetric eigenvalue problems ($2\times2$ in 1D; $(d+1)\times(d+1)$ in $d$ dimensions). Relaxed variants are further introduced to reduce computational costs in multidimensions while retaining the PCP guarantee. Extensive numerical benchmarks ranging from 1D to 3D, including ultra-relativistic Riemann problems and astrophysical jets, demonstrate that the proposed method robustly enforces physical admissibility, sharply resolves discontinuities, and maintains design-order accuracy for smooth solutions.

Manting Peng, Kailiang Wu

This paper establishes the first rigorous superconvergence theory for semidiscrete and fully discrete central discontinuous Galerkin (CDG) methods for linear hyperbolic equations on overlapping meshes. While the optimal $L^2$ convergence of $\mathbb{Q}^k$ CDG schemes was established on uniform Cartesian meshes by Liu, Shu, and Zhang [ SIAM J. Numer. Anal.}, 56 (2018), pp. 520--541], their observed $\mathcal{O}(h^{k+2})$ pointwise superconvergence has remained unproven, due to the loss of standard single-mesh Galerkin orthogonality inherent in the CDG overlapping structure. To overcome this fundamental barrier, we introduce a projection-correction framework that identifies a hidden superconvergent mechanism: an asymptotic weak residual cancellation in one dimension, and a high-order cancellation-by-aggregation (HOCA) mechanism in multiple dimensions. This HOCA approach overcomes the analytical challenge posed by coupled primal-dual directional residuals, recovering critical error cancellation properties absent from the standard variational formulation. Consequently, we provide the rigorous proof of the conjectured $\mathcal{O}(h^{k+2})$ pointwise superconvergence in the discrete $\ell^{\infty}$ norm across all superconvergent points. Furthermore, we reveal that under a systematically corrected initialization, this framework yields a previously undiscovered, stronger cell-average superconvergence estimate of order $\mathcal{O}(h^{\min\{2k+1,k+3\}})$. The theory is extended to fully discrete explicit Runge--Kutta CDG schemes, where stagewise corrected errors are constructed to preserve spatial superconvergence up to temporal truncation errors, yielding a stable reconstruction-based postprocessing estimate. Numerical experiments in one and two spatial dimensions confirm the sharpness of the theoretical rates.