pith. sign in

arxiv: 2605.22274 · v1 · pith:AA24DVHBnew · submitted 2026-05-21 · ⚛️ physics.flu-dyn · physics.comp-ph

A unified gas-kinetic wave-particle method for multiscale binary-species gas mixtures

Pith reviewed 2026-05-22 02:54 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn physics.comp-ph
keywords unified gas-kinetic methodwave-particle methodbinary-species mixturesmultiscale flowshypersonic flowDSMC validationrarefied gas dynamicsdiffusion modeling
0
0 comments X

The pith

A unified gas-kinetic wave-particle method models binary-species gas mixtures from continuum to rarefied regimes while capturing inter-species velocity and temperature differences.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops a unified gas-kinetic wave-particle method that automatically splits the gas distribution function into analytical hydrodynamic waves for near-equilibrium behavior and discrete particles for non-equilibrium parts. The approach applies the Groppi model to set local equilibrium velocity and temperature, recovering proper viscosity and diffusion, while adding a Shakhov correction to match the heat conduction coefficient. Diffusion appears in both the source term through operator splitting and in flux calculations via the characteristic solution, with strict consistency enforced between the wave and particle representations. Tests across regimes demonstrate that the method tracks species-specific velocity and temperature differences and produces wall pressure, shear stress, and heat flux coefficients that align with DSMC data in hypersonic cases.

Core claim

The UGKWP method decomposes the distribution function into waves and particles that respectively handle near-equilibrium and non-equilibrium portions. It adopts the Groppi et al. model for the macroscopic velocity and temperature of the target equilibrium distribution to obtain correct viscosity and diffusion coefficients in the continuum limit and incorporates the Shakhov model to correct the Prandtl number for heat conduction. Diffusion is incorporated through both the source term and the flux evolution while preserving consistency between wave and particle descriptions; an improved collision-time model governs the free-transport time of high-speed particles. Numerical experiments confirm

What carries the argument

Automatic decomposition of the gas distribution function into analytical hydrodynamic waves and discrete particles, with consistent source-term and flux modeling based on the Groppi equilibrium and Shakhov correction.

If this is right

  • The method reproduces species-specific velocity and temperature differences throughout the continuum-to-rarefied range.
  • Wall pressure, shear stress, and heat flux coefficients for hypersonic flows match DSMC reference data.
  • Diffusion effects are captured simultaneously in the source term and in the characteristic flux integral without breaking wave-particle consistency.
  • An improved collision-time correction governs free transport of high-speed particles.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same wave-particle split could be applied to plasma or radiation transport problems mentioned as target applications.
  • Because the decomposition is automatic, the scheme may eliminate the need for explicit regime-switching in mixed-scale simulations.
  • The enforced consistency between wave and particle parts suggests the method could serve as a stable platform for adding additional physical models such as chemical reactions.

Load-bearing premise

The Groppi model together with the Shakhov correction recovers the correct viscosity, diffusion, and heat conduction while the wave and particle descriptions remain strictly consistent in both source and flux.

What would settle it

A hypersonic binary-species simulation in which the computed species velocity and temperature profiles or the wall heat-flux coefficient deviates measurably from corresponding DSMC results.

Figures

Figures reproduced from arXiv: 2605.22274 by Chengwen Zhong, Junzhe Cao, Kun Xu, Wenpei Long, Yufeng Wei.

Figure 1
Figure 1. Figure 1: Diagram to illustrate the algorithm of UGKWP metho [PITH_FULL_IMAGE:figures/full_fig_p013_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Shock structure in binary gas mixture (Ma− = 1.5, mB/mA = 0.5, χ − B = 0.9): (a) Number density and (b) temperature. λ∞ (a) λ∞ (b) [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Shock structure in binary gas mixture (Ma− = 1.5, mB/mA = 0.5, χ − B = 0.1): (a) Number density and (b) temperature. 15 [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Shock structure in binary gas mixture (Ma− = 3, mB/mA = 0.5, χ − B = 0.9): (a) Number density and (b) temperature. λ∞ (a) λ∞ (b) [PITH_FULL_IMAGE:figures/full_fig_p016_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Shock structure in binary gas mixture (Ma− = 3, mB/mA = 0.5, χ − B = 0.1): (a) Number density and (b) temperature. λ∞ τ τ (a) λ∞ τ τ (b) [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Shock structure in binary gas mixture (Ma− = 1.5, mB/mA = 0.5, χ − B = 0.9): (a) Number density and (b) temperature. 16 [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Shock structure in binary gas mixture (Ma− = 1.5, mB/mA = 0.5, χ − B = 0.1): (a) Number density and (b) temperature. λ∞ τ τ (a) λ∞ τ τ (b) [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Shock structure in binary gas mixture (Ma− = 3, mB/mA = 0.5, χ − B = 0.9): (a) Number density and (b) temperature. λ∞ τ τ (a) λ∞ τ τ (b) [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Shock structure in binary gas mixture (Ma− = 3, mB/mA = 0.5, χ − B = 0.1): (a) Number density and (b) temperature. 17 [PITH_FULL_IMAGE:figures/full_fig_p017_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Mole fraction profile of mass diffusion case in an Ar [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Velocity profile of mass diffusion case in an Ar-Ne m [PITH_FULL_IMAGE:figures/full_fig_p019_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Velocity profile of Couette flow in an Ar-Ne mixture [PITH_FULL_IMAGE:figures/full_fig_p020_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Hypersonic flow around a cylinder in an Ar-Ne mixtu [PITH_FULL_IMAGE:figures/full_fig_p021_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Hypersonic flow around a cylinder in an Ar-Ne mixtu [PITH_FULL_IMAGE:figures/full_fig_p022_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Hypersonic flow around a cylinder in an Ar-Ne mixtu [PITH_FULL_IMAGE:figures/full_fig_p022_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Hypersonic flow around a cylinder in an Ar-Ne mixtu [PITH_FULL_IMAGE:figures/full_fig_p023_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Hypersonic flow around a cylinder in an Ar-Ne mixtu [PITH_FULL_IMAGE:figures/full_fig_p023_17.png] view at source ↗
read the original abstract

This paper presents a unified gas-kinetic wave-particle (UGKWP) method for simulating multiscale binary-species gas mixtures. Benefiting from direct modeling in a discretized space, the UGKWP method enables the automatic decomposition of the gas distribution function into analytical hydrodynamic waves and discrete particles, which respectively describe its near-equilibrium and non-equilibrium parts. This approach offers significant advantages for simulating various multiscale physical phenomena, such as hypersonic flows, plasma transport, and radiation transport. In this study, we employ the model proposed by Groppi et al. [EPL, 96 (2011) 64002] to calculate the macroscopic velocity and temperature of the local target equilibrium distribution function, thereby recovering the correct viscosity and diffusion coefficients in the continuum flow regime. To address the heat conduction coefficient, the Shakhov model is incorporated to correct the Prandtl number. Diffusion effects are accounted for not only in the source term via an operator-splitting method, but also in the flux evolution through the characteristic integral solution, while strictly maintaining consistency between the wave and particle descriptions. Furthermore, the microscopic model for high-speed particles is improved by utilizing a physically corrected collision time to determine their free-transport time. Through a series of numerical tests spanning the continuum to rarefied regimes, the proposed UGKWP method is shown to accurately capture the differences in velocity and temperature between different species. Notably, for hypersonic flows, the predicted wall pressure, shear stress, and heat flux coefficients agree well with DSMC results.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a unified gas-kinetic wave-particle (UGKWP) method for multiscale binary-species gas mixtures. It extends single-species UGKWP by adopting the Groppi et al. model to set the macroscopic velocity and temperature of the local target equilibrium distribution function, thereby recovering viscosity and diffusion coefficients, and incorporates the Shakhov correction to adjust the Prandtl number for heat conduction. Diffusion is treated both via operator splitting in the source term and through the characteristic integral in the flux, with the claim of strict consistency between the analytical wave and discrete particle descriptions. A physically corrected collision time is used for high-speed particles. Numerical tests across continuum to rarefied regimes, including hypersonic flows, show the method captures species velocity and temperature differences and yields wall pressure, shear stress, and heat flux coefficients in good agreement with DSMC.

Significance. If the wave-particle consistency and simultaneous recovery of all three transport coefficients are rigorously established for binary mixtures, the method would provide an efficient, automatic multiscale solver for gas-mixture problems in hypersonic aerodynamics and plasma transport without explicit regime switching. The direct use of established Groppi and Shakhov models together with DSMC validation in challenging flows indicates practical utility for engineering applications.

major comments (2)
  1. [§3.2] §3.2 (Target equilibrium distribution): The Groppi et al. model is used to match viscosity and diffusion, after which the Shakhov correction is applied to the heat-flux moment. For binary mixtures the species collision integrals couple momentum and energy exchange; no explicit first-order Chapman-Enskog expansion or moment calculation is shown to confirm that the diffusion flux remains unaltered by the Prandtl-number correction while wave-particle consistency is preserved.
  2. [§4.3] §4.3 (Flux evolution and consistency): Diffusion is stated to appear in both the operator-split source term and the characteristic flux integral. The manuscript does not provide a direct verification that the species diffusion moment extracted from the wave part equals that from the particle part at O(Kn) after the Shakhov modification is introduced.
minor comments (2)
  1. [§2.4] The definition of the physically corrected collision time for high-speed particles is introduced in §2.4 but its explicit functional form is not written out; adding the formula would improve reproducibility.
  2. [Figure 7] Figure 7 (hypersonic cylinder): the comparison with DSMC would be strengthened by reporting quantitative L2 errors or relative differences for the wall coefficients rather than qualitative visual agreement alone.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments correctly identify points where additional explicit verification would strengthen the presentation. We address each major comment below and will incorporate the requested clarifications in the revised manuscript.

read point-by-point responses
  1. Referee: [§3.2] §3.2 (Target equilibrium distribution): The Groppi et al. model is used to match viscosity and diffusion, after which the Shakhov correction is applied to the heat-flux moment. For binary mixtures the species collision integrals couple momentum and energy exchange; no explicit first-order Chapman-Enskog expansion or moment calculation is shown to confirm that the diffusion flux remains unaltered by the Prandtl-number correction while wave-particle consistency is preserved.

    Authors: We agree that an explicit moment calculation would improve clarity. The Groppi model is constructed so that the first-order Chapman-Enskog expansion of the target equilibrium recovers the correct species diffusion flux and mixture viscosity; the Shakhov correction is applied only to the heat-flux term and does not alter the lower-order moments that determine diffusion. Wave-particle consistency is maintained because both the analytical wave and the sampled particles are generated from the identical modified target distribution and collision time. In the revision we will add a short derivation of the relevant moments in §3.2 to make this explicit. revision: yes

  2. Referee: [§4.3] §4.3 (Flux evolution and consistency): Diffusion is stated to appear in both the operator-split source term and the characteristic flux integral. The manuscript does not provide a direct verification that the species diffusion moment extracted from the wave part equals that from the particle part at O(Kn) after the Shakhov modification is introduced.

    Authors: We appreciate the request for direct verification. Because the wave and particle components are obtained from the same integral solution of the kinetic equation with the Groppi-Shakhov target, their O(Kn) diffusion moments are identical by construction. To address the concern explicitly, we will include a brief analytical comparison of the extracted species diffusion moments from the wave and particle contributions, together with a numerical check on a simple binary relaxation test, in the revised §4.3 or an appendix. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation adopts external models and maintains consistency by construction of the unified scheme

full rationale

The paper explicitly adopts the Groppi et al. model for species velocity/temperature (to recover viscosity and diffusion) and the Shakhov correction for Prandtl number from independent literature citations. These are not derived internally or fitted to the target outputs. The claimed wave-particle consistency in source term and flux evolution follows directly from the UGKWP discretization framework (prior independent work), with no reduction of predictions to self-defined inputs or self-citation chains that bear the central load. Numerical validation against DSMC provides external falsifiability. No self-definitional loops, fitted inputs renamed as predictions, or ansatz smuggling appear in the provided derivation steps.

Axiom & Free-Parameter Ledger

1 free parameters · 3 axioms · 0 invented entities

The central claim depends on the accuracy of the Groppi model for recovering viscosity and diffusion, the Shakhov correction for heat conduction, and the assumption that wave and particle descriptions remain consistent under operator splitting and characteristic integration.

free parameters (1)
  • physically corrected collision time
    Used to set free-transport time for high-speed particles; introduced to improve the microscopic model.
axioms (3)
  • domain assumption Groppi et al. model recovers correct viscosity and diffusion coefficients when used for local target equilibrium velocity and temperature
    Invoked to ensure proper continuum transport properties in the wave description.
  • domain assumption Shakhov model corrects the Prandtl number to obtain the proper heat conduction coefficient
    Added specifically to address heat conduction while keeping other transport coefficients intact.
  • domain assumption Consistency between wave and particle descriptions is maintained when diffusion is included in both source term and flux evolution
    Required for the unified treatment to remain valid across regimes.

pith-pipeline@v0.9.0 · 5819 in / 1705 out tokens · 45668 ms · 2026-05-22T02:54:29.534149+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

65 extracted references · 65 canonical work pages

  1. [1]

    C. Liu, Y. Zhu, K. Xu, Unified gas-kinetic wave-particle m ethods I: Contin- uum and rarefied gas flow, Journal of Computational Physics 40 1 (2020) 108977. doi:https://doi.org/10.1016/j.jcp.2019.108977. 24

  2. [2]

    Y. Zhu, C. Liu, C. Zhong, K. Xu, Unified gas-kinetic wave-p article methods. II. Multiscale simulation on unstructured mesh, Physics of Flu ids 31 (2019) 067105. doi:https://doi.org/10.1063/1.5097645

  3. [3]

    X. Xu, Y. Chen, C. Liu, Z. Li, K. Xu, Unified gas-kinetic wav e-particle meth- ods V: Diatomic molecular flow, Journal of Computational Phy sics 442 (2021) 110496. doi:https://doi.org/10.1016/j.jcp.2021.110496

  4. [4]

    Y. Wei, Y. Zhu, K. Xu, Unified gas-kinetic wave-particle m ethods VII: Diatomic gas with rota- tional and vibrational nonequilibrium, Journal of Computa tional Physics 497 (2024) 112610. doi:https://doi.org/10.1016/j.jcp.2023.112610

  5. [5]

    Y. Chen, Y. Zhu, K. Xu, A three-dimensional unified gas-ki netic wave-particle solver for flow computation in all regimes, Physics of Fluids 32 (9) (2020) 096108. doi:https://doi.org/10.1063/5.0021199

  6. [6]

    X. Yang, C. Liu, X. Ji, W. Shyy, K. Xu, Unified gas-kinetic w ave-particle methods VI: Disperse dilute gas-particle multiphase flow, Communications in Com putational Physics 31 (3) (2022) 669–706. doi:https://doi.org/10.4208/cicp.OA-2021-0153

  7. [7]

    X. Yang, W. Shyy, K. Xu, Unified gas-kinetic wave-particl e method for polydisperse gas-solid particle multiphase flow, Journal of Fluid Mechan ics 983 (3) (2024) A37. doi:https://doi.org/10.1017/jfm.2024.80

  8. [8]

    W. Li, C. Liu, Y. Zhu, J. Zhang, K. Xu, Unified gas-kinetic w ave-particle methods III: Multiscale photon transport, Journal of Computationa l Physics 408 (2020) 109280. doi:https://doi.org/10.1016/j.jcp.2020.109280

  9. [9]

    C. Liu, W. Li, Y. Wang, P. Song, K. Xu, An implicit unified ga s-kinetic wave- particle method for radiative transport process, Physics o f Fluids 35 (11) (2023) 112013. doi:https://doi.org/10.1063/5.0174774

  10. [10]

    X. Yang, Y. Zhu, C. Liu, K. Xu, Unified gas-kinetic wave-p article method for frequency- dependent radiation transport equation, Journal of Comput ational Physics 522 (2025) 113587. doi:https://doi.org/10.1016/j.jcp.2024.113587

  11. [11]

    C. Liu, K. Xu, Unified gas-kinetic wave-particle method s IV: multi-species gas mixture and plasma transport, Advances in Aerodynamics 3 (2021) 9. doi:https://doi.org/10.1186/s42774-021-00062-1

  12. [12]

    Z. Pu, K. Xu, Unified gas-kinetic wave-particle method f or multiscale flow simula- tion of partially ionized plasma, Journal of Computational Physics 530 (2025) 113918. doi:https://doi.org/10.1016/j.jcp.2025.113918

  13. [13]

    Z. Pu, K. Xu, Electromagnetic flow control in hypersonic rarefied environment, Journal of Fluid Mechanics 1033 (2026) A37. doi:https://doi.org/10.1017/jfm.2026.11462

  14. [14]

    X. Yang, K. Xu, Wave-particle based multiscale modelin g and simulation of non-equilibrium turbulent flows (2025). arXiv:2503.07207. 25

  15. [15]

    X. Yang, K. Xu, Wave-particle turbulence simulation of spatially developing round jets: Turbulent flow modeling and method validation, Physics of Fl uids 37 (10) (2025) 101706. doi:https://doi.org/10.1063/5.0293211

  16. [16]

    Xu, Direct modeling for computational fluid dynamics : Construction and application of unified gas-kinetic schemes, World Scientific, 2015

    K. Xu, Direct modeling for computational fluid dynamics : Construction and application of unified gas-kinetic schemes, World Scientific, 2015

  17. [17]

    Xu, A unified computational fluid dynamics framework f rom rarefied to continuum regimes, Cambridge University Press, 2021

    K. Xu, A unified computational fluid dynamics framework f rom rarefied to continuum regimes, Cambridge University Press, 2021

  18. [18]

    K. Xu, A gas-kinetic BGK scheme for the Navier-Stokes eq uations and its connection with artificial dissipation and Godunov method, Journal of Compu tational Physics 171 (20) (2001) 289–335. doi:https://doi.org/10.1006/jcph.2001.6790

  19. [19]

    S. Liu, C. Zhong, M. Fang, Simplified unified wave-partic le method with quantified model- competition mechanism for numerical calculation of multis cale flows, Physical Review E 102 (2020) 013304. doi:https://doi.org/10.1103/PhysRevE.102.013304

  20. [20]

    L. M. Yang, Z. H. Li, C. Shu, Y. Y. Liu, J. Wu, Discrete unifi ed gas-kinetic wave- particle method for flows in all flow regimes, Physical Review E 108 (2023) 015302. doi:https://doi.org/10.1103/PhysRevE.108.015302

  21. [21]

    W. Liu, C. Shu, C. J. Teo, Z. Y. Yuan, Y. Y. Liu, A simple hyd rodynamic- particle method for supersonic rarefied flows, Physics of Flu ids 34 (5) (2022) 057101. doi:https://doi.org/10.1063/5.0088946

  22. [22]

    Z. Guo, Y. Zhu, K. Xu, Kinetic representation of the unifi ed gas-kinetic wave-particle method and beyond, Communications in Computational Physic s 39 (5) (2026) 1512–1535. doi:https://doi.org/10.4208/cicp.OA-2024-0192

  23. [23]

    Z. Guo, K. Xu, Y. Zhu, A unified gas-kinetic framework fro m Boltz- mann to Navier-Stokes scales, Advances in Aerodynamics 8 (2 026) 8. doi:https://doi.org/10.1186/s42774-025-00252-1

  24. [24]

    Xu, J.-C

    K. Xu, J.-C. Huang, A unified gas-kinetic scheme for cont inuum and rar- efied flows, Journal of Computational Physics 229 (20) (2010) 7747–7764. doi:https://doi.org/10.1016/j.jcp.2010.06.032

  25. [25]

    W. Long, Y. Wei, K. Xu, Nonequilibrium flow simulations u sing unified gas-kinetic wave-particle method, AIAA Journal 62 (4) (202 4) 1411–1433. doi:https://doi.org/10.2514/1.J063641

  26. [26]

    Y. Wei, J. Cao, X. Ji, K. Xu, Adaptive wave-particle deco mposition in UGKWP method for high-speed flow simulations, Advances in Aerodyn amics 5 (2023) 25. doi:https://doi.org/10.1186/s42774-023-00156-y

  27. [27]

    J. Cao, Y. Wei, W. Long, C. Zhong, K. Xu, Adaptive criteri on and modification of wave- particle decomposition in UGKWP method for high-speed flow s imulation, Computers & Fluids 305 (2026) 106896. doi:https://doi.org/10.1016/j.compfluid.2025.106896. 26

  28. [28]

    E. H. Hirschel, C. Weiland, Selected aerothermodynami c design problems of hypersonic flight vehicles, Springer Berlin, 2009

  29. [29]

    Andries, K

    P. Andries, K. Aoki, B. Perthame, A consistent BGK-type model for gas mixtures, Journal of Statistical Physics 106 (2002) 993–1018. doi:https://doi.org/10.1023/A:1014033703134

  30. [30]

    Groppi, S

    M. Groppi, S. Monica, G. Spiga, A kinetic ellipsoidal BG K model for a binary gas mixture, Europhysics Letters 96 (6) (2011) 6 4002. doi:https://doi.org/10.1209/0295-5075/96/64002

  31. [31]

    Brull, V

    S. Brull, V. Pavan, J. Schneider, Derivation of a BGK mod el for mix- tures, European Journal of Mechanics - B/Fluids 33 (2012) 74 –86. doi:https://doi.org/10.1016/j.euromechflu.2011.12.003

  32. [32]

    Brull, An ellipsoidal statistical model for gas mixt ures, Communications in Mathematical Sciences 13 (1) (2015) 1–13

    S. Brull, An ellipsoidal statistical model for gas mixt ures, Communications in Mathematical Sciences 13 (1) (2015) 1–13. doi:https://dx.doi.org/10.4310/CMS.2015.v13.n1.a1

  33. [33]

    B. N. Todorova, R. Steijl, Derivation and numerical com parison of Shakhov and Ellipsoidal Statistical kinetic models for a monoatomic gas mixture, Eu ropean Journal of Mechanics - B/Fluids 76 (2019) 390–402. doi:https://doi.org/10.1016/j.euromechflu.2019.04.001

  34. [34]

    F. Hild, M. Pfeiffer, Multi-species modeling in the part icle-based ellipsoidal statistical Bhatnagar-Gross-Krook method including internal degrees of freedom, Journal of Compu- tational Physics 514 (2024) 113226. doi:https://doi.org/10.1016/j.jcp.2024.113226

  35. [35]

    F. Fei, D. Liu, L. Xie, Z. Ren, Y. Hu, A unified stochastic p article method for polyatomic gas mixtures, Computer Physics Communicati ons 321 (2026) 110029. doi:https://doi.org/10.1016/j.cpc.2026.110029

  36. [36]

    F. J. McCormack, Construction of linearized kinetic mo dels for gaseous mix- tures and molecular gases, Physics of Fluids 16 (12) (1973) 2 095–2105. doi:https://doi.org/10.1063/1.1694272

  37. [37]

    J. R. Haack, C. D. Hauck, M. S. Murillo, A conservative, e ntropic mul- tispecies BGK model, Journal of Statistical Physics 168 (20 17) 826–856. doi:https://doi.org/10.1007/s10955-017-1824-9

  38. [38]

    A. V. Bobylev, M. Bisi, M. Groppi, G. Spiga, I. F. Potapen ko, A general consistent BGK model for gas mixtures, Kinetic and Related Models 11 (6) (2018) 1377–1393. doi:https://doi.org/10.3934/krm.2018054

  39. [39]

    Q. Li, J. Zeng, L. Wu, Kinetic modelling of rarefied gas mi xtures with disparate mass in strong non-equilibrium flows, Journal of Fluid Mecha nics 1001 (2024) A5. doi:https://doi.org/10.1017/jfm.2024.1047

  40. [40]

    Zhang, L

    Y. Zhang, L. Zhu, P. Wang, Z. Guo, Discrete unified gas kin etic scheme for flows of bi- nary gas mixture based on the McCormack model, Physics of Flu ids 31 (1) (2019) 017101. doi:https://doi.org/10.1063/1.5063846. 27

  41. [41]

    J. Zeng, Q. Li, L. Wu, General synthetic iterative schem e for rarefied gas mixture flows, Journal of Computational Physics 519 (202 4) 113420. doi:https://doi.org/10.1016/j.jcp.2024.113420

  42. [42]

    Wang, Unified gas-kinetic scheme for the study of non- equilibrium flows, Ph.D

    R. Wang, Unified gas-kinetic scheme for the study of non- equilibrium flows, Ph.D. thesis, The Hong Kong University of Science and Technology (2015)

  43. [43]

    Zhang, L

    Y. Zhang, L. Zhu, R. Wang, Z. Guo, Discrete unified gas kin etic scheme for all Knudsen number flows. III. Binary gas mixtures of maxwell molecules, Physical Review E 97 (2018) 053306. doi:https://doi.org/10.1103/PhysRevE.97.053306

  44. [44]

    Z. Xin, Y. Zhang, Z. Guo, A discrete unified gas-kinetic s cheme for multi-species rarefied flows, Advances in Aerodynamics 5 (20 23) 5. doi:https://doi.org/10.1186/s42774-022-00135-9

  45. [45]

    C. Liu, K. Xu, A unified gas kinetic scheme for continuum a nd rarefied flows V: Multiscale and multi-component plasma transport, Communications in C omputational Physics 22 (5) (2017) 1175–1223. doi:https://doi.org/10.4208/cicp.OA-2017-0102

  46. [46]

    M. Quan, X. Yang, Y. Wei, K. Xu, Radiative hydrodynamic e quations with nonequilibrium radiative transfer, Physical Review D 111 (2025) 103038. doi:https://doi.org/10.1103/PhysRevD.111.103038

  47. [47]

    B. N. Todorova, C. White, R. Steijl, Numerical evaluati on of novel kinetic models for binary gas mixture flows, Physics of Fluids 32 (1) ( 2020) 016102. doi:https://doi.org/10.1063/1.5134040

  48. [48]

    E. M. Shakhov, Generalization of the Krook kinetic rela xation equation, Fluid Dynamics 3 (1968) 95–96. doi:https://doi.org/10.1007/BF01029546

  49. [49]

    X. Xu, Y. Chen, K. Xu, Modeling and computation for non-e quilibrium gas dynamics: Beyond single relaxation time kinetic models, Physics of Fl uids 33 (1) (2021) 011703. doi:https://doi.org/10.1063/5.0036203

  50. [50]

    G. A. Bird, Molecular gas dynamics and the direct simula tion of gas flows, Clarendon Press, 1994

  51. [51]

    Fei, A Navier-Stokes asymptotic preserving direct s imulation Monte Carlo method for multi-species gas flows, Journal of Computational Physi cs 538 (2025) 114196

    F. Fei, A Navier-Stokes asymptotic preserving direct s imulation Monte Carlo method for multi-species gas flows, Journal of Computational Physi cs 538 (2025) 114196. doi:https://doi.org/10.1016/j.jcp.2025.114196

  52. [52]

    L. Luo, J. Zeng, Y. Zhang, W. Li, Q. Li, L. Wu, Enhancing DS MC sim- ulations of rarefied gas mixtures using a fast-converging an d asymptotic-preserving scheme, Computer Methods in Applied Mechanics and Engineer ing 449 (2026) 118508. doi:https://doi.org/10.1016/j.cma.2025.118508

  53. [53]

    P. L. Bhatnagar, E. P. Gross, M. Krook, A model for collis ion processes in gases. I. Small amplitude processes in charged and neutral one-component s ystems, Physical Review Journals Archive 94 (1954) 511–525. doi:https://doi.org/10.1103/PhysRev.94.511. 28

  54. [54]

    Weissman, E

    S. Weissman, E. A. Mason, Determination of gaseous-diff usion coefficients from vis- cosity measurements, The Journal of Chemical Physics 37 (6) (1962) 1289–1300. doi:https://doi.org/10.1063/1.1733277

  55. [55]

    C. R. Wilke, A viscosity equation for gas mixtures, The J ournal of Chemical Physics 18 (4) (1950) 517–519. doi:https://doi.org/10.1063/1.1747673

  56. [56]

    V. A. Wassiljewa, Wärmeleitung in gasgemischen, Physi kalische Zeitschrift 5 (1904) 737–742

  57. [57]

    E. A. Mason, S. C. Saxena, Approximate formula for the th ermal conductivity of gas mixtures, The Physics of Fluids 1 (5) (1958) 361–369. doi:https://doi.org/10.1063/1.1724352

  58. [58]

    B. E. Poling, J. M. Prausnitz, J. P. O’Connell, The prope rties of gases and liquids, McGraw- Hill Education, 2001

  59. [59]

    G. May, B. Srinivasan, A. Jameson, An improved gas-kine tic BGK finite-volume method for three-dimensional transonic flow, Journal of Computationa l Physics 220 (2) (2007) 856–878. doi:https://doi.org/10.1016/j.jcp.2006.05.027

  60. [60]

    Xu, Gas-kinetic schemes for unsteady compressible fl ow simulations, in: 29th CFD Lecture Series 1998-03 von Kárman Institute for Fluid Dynamics, Bel gium, 1998

    K. Xu, Gas-kinetic schemes for unsteady compressible fl ow simulations, in: 29th CFD Lecture Series 1998-03 von Kárman Institute for Fluid Dynamics, Bel gium, 1998

  61. [61]

    S. Liu, Y. Yang, C. Zhong, An extended gas-kinetic schem e for shock structure calculations, Journal of Computational Physics 390 (2019) 1–24. doi:https://doi.org/10.1016/j.jcp.2019.04.016

  62. [62]

    S. Liu, C. Zhong, Investigation of the kinetic model equ ations, Physical Review E 89 (2014) 033306. doi:https://doi.org/10.1103/PhysRevE.89.033306

  63. [63]

    Xu, J.-C

    K. Xu, J.-C. Huang, An improved unified gas-kinetic sche me and the study of shock structures, IMA Journal of Applied Mathematics 76 (5) (2011) 698–711. doi:https://doi.org/10.1093/imamat/hxr002

  64. [64]

    Kosuge, K

    S. Kosuge, K. Aoki, S. Takata, Shock-wave structure for a binary gas mixture: finite-difference analysis of the Boltzmann equation for hard-sphere molecul es, European Journal of Mechanics - B/Fluids 20 (1) (2001) 87–126. doi:https://doi.org/10.1016/S0997-7546(00)00133-3

  65. [65]

    M. T. Ho, L. Wu, I. Graur, Y. Zhang, J. M. Reese, Comparati ve study of the Boltzmann and McCormack equations for Couette and Fouri er flows of binary gaseous mixtures, International Journal of Heat and Mass Tr ansfer 96 (2016) 29–41. doi:https://doi.org/10.1016/j.ijheatmasstransfer.2015.12.068. 29