This paper presents a mass-lumped Virtual Element Method (VEM) with explicit Strong Stability-Preserving Runge-Kutta (SSP-RK) time integration for two-dimensional parabolic problems on general polygonal meshes. A diagonal mass matrix is constructed via row-sum operations combined with flooring to ensure uniform positivity. Stabilization terms vanish identically under row summation, so the lumped weights derive solely from the $L^2$ projector and are computable through a small polynomial system at cost $\mathcal{O}(N_k^3)$ per element. The resulting lumped bilinear form satisfies $L^2$-equivalence with edge-count-independent constants, yielding a symmetric positive definite discrete inner product. A mesh-robust spectral bound $\lambda_{\max}\big((\hat{\mathbf{M}}_h)^{-1}\mathbf{K}_h\big) \le C_{\mathrm{inv}}^2/\hat{\beta}_* \cdot h^{-2}$ is established with constants depending only on spatial dimension, polynomial order, and mesh regularity. This delivers the classical diffusion-type CFL condition $\Delta t=\mathcal{O}(h^2)$ for forward Euler stability and extends to higher-order SSP-RK schemes, guaranteeing preservation of energy decay, positivity, and discrete maximum principles. Numerical experiments on distorted quadrilaterals, serendipity elements, and Voronoi polygons validate the theoretical predictions: the lumped VEM with $k=1$ achieves optimal convergence rates ($\mathcal{O}(h)$ in $H^1$, $\mathcal{O}(h^2)$ in $L^2$) with no degradation from geometric distortion or mass lumping, while SSP-RK integrators remain stable under the predicted $\Delta t\propto h^{2}$ scaling
翻译:暂无翻译