Multiscale Partially Explicit Splitting with Mass Lumping for High-Contrast Wave Equations

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Li, Shu Fan, Leung, Wing Tat
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866912242171117568
author Li, Shu Fan
Leung, Wing Tat
author_facet Li, Shu Fan
Leung, Wing Tat
contents In this paper, contrast-independent partially explicit time discretization for wave equations in heterogeneous high-contrast media via mass lumping is concerned. By employing a mass lumping scheme to diagonalize the mass matrix, the matrix inversion procedures can be avoided, thereby significantly enhancing computational efficiency especially in the explicit part. In addition, after decoupling the resulting system, higher order time discretization techniques can be applied to get better accuracy within the same time step size. Furthermore, the spatial space is divided into two components: contrast-dependent ("fast") and contrast-independent ("slow") subspaces. Using this decomposition, our objective is to introduce an appropriate time splitting method that ensures stability and guarantees contrast-independent discretization under suitable conditions. We analyze the stability and convergence of the proposed algorithm. In particular, we discuss the second order central difference and higher order Runge-Kutta method for a wave equation. Several numerical examples are presented to confirm our theoretical results and to demonstrate that our proposed algorithm achieves high accuracy while reducing computational costs for high-contrast problems.
format Preprint
id arxiv_https___arxiv_org_abs_2502_16403
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Multiscale Partially Explicit Splitting with Mass Lumping for High-Contrast Wave Equations
Li, Shu Fan
Leung, Wing Tat
Numerical Analysis
In this paper, contrast-independent partially explicit time discretization for wave equations in heterogeneous high-contrast media via mass lumping is concerned. By employing a mass lumping scheme to diagonalize the mass matrix, the matrix inversion procedures can be avoided, thereby significantly enhancing computational efficiency especially in the explicit part. In addition, after decoupling the resulting system, higher order time discretization techniques can be applied to get better accuracy within the same time step size. Furthermore, the spatial space is divided into two components: contrast-dependent ("fast") and contrast-independent ("slow") subspaces. Using this decomposition, our objective is to introduce an appropriate time splitting method that ensures stability and guarantees contrast-independent discretization under suitable conditions. We analyze the stability and convergence of the proposed algorithm. In particular, we discuss the second order central difference and higher order Runge-Kutta method for a wave equation. Several numerical examples are presented to confirm our theoretical results and to demonstrate that our proposed algorithm achieves high accuracy while reducing computational costs for high-contrast problems.
title Multiscale Partially Explicit Splitting with Mass Lumping for High-Contrast Wave Equations
topic Numerical Analysis
url https://arxiv.org/abs/2502.16403