A wave-particle model for reacting mixtures
An arXiv preprint reports that an extended UGKWP method agreed with DSMC reference simulations, the benchmark runs used in the study, across reacting gas-mixture tests. The comparisons covered flow fields, species mole fractions and wall quantities in hypersonic flow around a cylinder, one-dimensional shock structures and a three-dimensional side jet over a blunt cone.
UGKWP combines two ways of representing the gas. Hydrodynamic waves represent near-equilibrium distributions, while numerical particles represent nonequilibrium ones. The study extends this wave-particle split to multiscale gas-mixture flows with a single elementary reaction.
Chemical source terms are applied first at the macroscopic level, meaning the bulk-flow update, and then incorporated into the wave-particle update. In the monatomic setting, particles in free transport remain chemically inactive. The multispecies extension covers more than two species under a non-extreme mass-ratio condition.
The numerical cases use four species and the same Zeldovich-type elementary reaction: O2 + N ⇋ NO + O. The chemical scope is therefore one elementary reaction.
Cylinder profiles differed by reaction-energy condition
The baseline cylinder validation used a freestream Mach number of 9 and a Knudsen number of 0.1. The Knudsen number is used here to distinguish more rarefied and more continuum-like conditions. The study reported agreement with DSMC for the flow field, species mole fractions and wall quantities in this case.
Compared with the ΔE = 0 J reference case, the forward exothermic, or heat-releasing, condition had a larger temperature peak, a farther shock location and a larger heat-flux peak. The forward endothermic, or heat-absorbing, condition had lower temperature and heat-flux peaks and a nearer shock.
The cylinder comparison also included Kn∞ = 1 and Kn∞ = 0.02. At Kn∞ = 0.02, the oxygen mole fraction, χO2, approached zero near stagnation; at Kn∞ = 1, product mole fractions reached about 10 percent. Both simulations were reported as accurate relative to DSMC.
Separate one-dimensional shock-structure tests at Ma = 3 and Ma = 5 were also reported to match DSMC.
The three-dimensional side-jet case
In a separate three-dimensional side-jet case over a blunt cone, the far field was pure O2 and the jet was full nitrogen at Mach 2. The freestream was Mach 10.15 with Kn∞ = 0.002; the jet density was 200 times the far-field density, and the mesh contained 309,240 cells.
The forward exothermic case showed a large high-temperature region and increased product concentrations near the jet exit and wake. The forward endothermic case showed limited reaction rates and small product mole fractions.
Reference particle counts were 800 per cell in the two- and three-dimensional cases and 5,000 per cell in the shock structures. The CFL number was 0.8, and isothermal walls used full accommodation.
Scope of the computational result
The reported agreement applies to a focused computational setup: four monatomic species, one elementary reaction and a treatment in which free-transport particles remain chemically inactive. The multispecies extension beyond two species assumes a non-extreme mass ratio.
The work is an arXiv preprint, version 1, dated 26 August 2026. It acknowledges support from the National Key R&D Program of China, the National Natural Science Foundation of China and the Hong Kong Research Grants Council. Supporting data are available from the corresponding author upon reasonable request.
Paper data and sources
Original title: A unified gas-kinetic wave-particle method for multiscale gas-mixture flow with an elementary chemical reaction
Authors: Cao Junzhe, Wei Yufeng, Long Wenpei et al.
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-26
DOI: Not available
Original paper · Full text