Rigorously justified local time stepping in the unified gas-kinetic wave-particle method for steady multiscale flow simulation
Abstract
Local time stepping (LTS) can accelerate convergence to steady states in kinetic simulations with large variations in the local time steps across the computational domain. When neighboring cells advance with unequal time steps, the time-averaged particle flux must be balanced across their common interface. For particle-based or hybrid wave-particle methods under finite volume method (FVM) framework, we rigorously establish a sufficient condition for time-averaged interfacial particle-flux balance: fixed positive local time steps together with proportional particle-mass scaling. When a particle crosses from cell L to cell R, its mass is scaled by ΔtR/ΔtL. In the unified gas-kinetic wave-particle (UGKWP) implementation, the same ratio is applied to the remaining free-transport time of the crossing particle. LTS also affects the wave-particle decomposition and the time integration of the wave fluxes in UGKWP. The wave-particle decomposition in cell i is determined by the local ratio Δti/τi, which better reflects the local relation between the observation scale and relaxation time for multiscale cases. The equilibrium and analytic free transport wave fluxes are integrated over and normalized by the corresponding cell-side time steps to obtain the interfacial time-averaged wave fluxes. The UGKWP-LTS method is used to simulate the hypersonic flow past a cylinder at Kn=0.01 and $0.1$, and a flat plate at Kn=0.0169. In all three cases, the surface quantities obtained with UGKWP-LTS agree well with the reference data. Relative to global time stepping (GTS), UGKWP-LTS achieves step-count speedups of 6.6×, 3.8×, and 20× for the three cases, respectively. The corresponding wall-clock speedups are 7.1×, 4.5×, and approximately 21×.
Related
- Unified Gas-Kinetic Wave-Particle Method for Multiscale Simulation of Vlasov-Poisson-Fokker-Planck System
- A wave-particle decomposition framework for multiscale kinetic transport: continuous-spectrum equations and coupled wave-particle iteration
- Dual Laser and Shadow Particle Image Velocimetry for Multiscale Fluid-Structure Interactions
- An implicit multiscale solver for steady kinetic equations based on wave-particle decomposition
- Multiscale Quantum-Classical Simulation of Protein Autoproteolysis in a Million-Atom Membrane Environment
BibTeX
@misc{guo2026rigorously,
title = {{Rigorously justified local time stepping in the unified gas-kinetic wave-particle method for steady multiscale flow simulation}},
author = {Wenzhi Guo and Junzhe Cao and Wenpei Long and Kun Xu},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2607.09552},
eprint = {2607.09552},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2607.09552},
}
RIS
TY - UNPB
TI - Rigorously justified local time stepping in the unified gas-kinetic wave-particle method for steady multiscale flow simulation
AU - Wenzhi Guo
AU - Junzhe Cao
AU - Wenpei Long
AU - Kun Xu
PY - 2026
DA - 2026/09/01
PB - arXiv
DO - 10.48550/arxiv.2607.09552
UR - https://arxiv.org/abs/2607.09552
AB - Local time stepping (LTS) can accelerate convergence to steady states in kinetic simulations with large variations in the local time steps across the computational domain. When neighboring cells advance with unequal time steps, the time-averaged particle flux must be balanced across their common interface. For particle-based or hybrid wave-particle methods under finite volume method (FVM) framework, we rigorously establish a sufficient condition for time-averaged interfacial particle-flux balance: fixed positive local time steps together with proportional particle-mass scaling. When a particle crosses from cell L to cell R, its mass is scaled by Δt_R/Δt_L. In the unified gas-kinetic wave-particle (UGKWP) implementation, the same ratio is applied to the remaining free-transport time of the crossing particle. LTS also affects the wave-particle decomposition and the time integration of the wave fluxes in UGKWP. The wave-particle decomposition in cell i is determined by the local ratio Δt_i/τ_i, which better reflects the local relation between the observation scale and relaxation time for multiscale cases. The equilibrium and analytic free transport wave fluxes are integrated over and normalized by the corresponding cell-side time steps to obtain the interfacial time-averaged wave fluxes. The UGKWP-LTS method is used to simulate the hypersonic flow past a cylinder at Kn=0.01 and $0.1$, and a flat plate at Kn=0.0169. In all three cases, the surface quantities obtained with UGKWP-LTS agree well with the reference data. Relative to global time stepping (GTS), UGKWP-LTS achieves step-count speedups of 6.6×, 3.8×, and 20× for the three cases, respectively. The corresponding wall-clock speedups are 7.1×, 4.5×, and approximately 21×.
ER -