pith. machine review for the scientific record. sign in

arxiv: 2605.03479 · v1 · submitted 2026-05-05 · ⚛️ physics.flu-dyn

Recognition: unknown

A Time-Domain Harmonic Balance Unified Gas-Kinetic Scheme for Temporally Periodic Flows Across all Knudsen Regimes

Authors on Pith no claims yet

Pith reviewed 2026-05-07 14:20 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords harmonic balanceunified gas-kinetic schemeperiodic flowsKnudsen regimestime-spectral methodcavity flowsmultiscale simulationgas kinetics
0
0 comments X

The pith

The HB-UGKS reformulates time-periodic flows as a block-coupled quasi-steady system to simulate them efficiently across all Knudsen regimes.

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

This paper introduces a harmonic balance version of the unified gas-kinetic scheme to compute flows that repeat regularly in time, whether the gas behaves like a continuous fluid or as individual particles. It rewrites the repeating problem as a set of linked steady equations joined by a time-spectral source term drawn from Fourier harmonics, then advances the whole block in pseudo-time with local stepping so every instant in the cycle is found together. The UGKS flux keeps transport and collisions coupled inside each evaluation, preserving physical consistency from near-continuum to free-molecular conditions. Validation on two cavity problems shows that the fundamental harmonic suffices for small-amplitude shear oscillations while several harmonics are needed for large thermal modulations, and the approach cuts wall-clock time by factors up to 9 in high-frequency cases compared with explicit time marching.

Core claim

The time-domain harmonic balance approach converts the unsteady periodic problem into a block-coupled quasi-steady system via a time-spectral source term; when this system is solved with the UGKS flux that maintains transport-collision consistency, the resulting scheme reproduces the full periodic flow field across the entire Knudsen range while allowing pseudo-time marching and concurrent resolution of all sub-time levels.

What carries the argument

The time-spectral source term generated by the Fourier decomposition, which couples the harmonic coefficients into one large block system that is marched in pseudo-time.

If this is right

  • For small-amplitude shear-driven oscillatory cavities the fundamental harmonic alone captures the anti-resonance phenomenon and matches hydrodynamic damping from linearized analyses across Knudsen and Strouhal numbers.
  • For thermally driven cavities with large temperature swings, higher-order harmonics are required to recover strong nonlinear waveform distortions and rarefaction effects.
  • Computational wall-clock time drops substantially relative to explicit time-domain methods, with the largest gains in high-frequency regimes reaching speedup factors of 9.0 and 8.26 in the two tested cases.
  • The multiscale fidelity of the underlying UGKS flux is retained without modification, so the same code covers the full range from hydrodynamic to free-molecular regimes.

Where Pith is reading between the lines

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

  • The same block-coupled formulation could be applied to other repeating rarefied-gas problems such as oscillating membranes or pulsed micro-jets, potentially enabling broad parameter sweeps over frequency and Knudsen number that explicit methods cannot afford.
  • Because all time levels are solved simultaneously, the scheme naturally supplies the entire periodic waveform, which could be used directly for stability analysis or for coupling to structural models without post-processing.
  • One testable extension is to examine whether a modest increase in retained harmonics can approximate slowly drifting or quasi-periodic flows without reverting to full unsteady integration.

Load-bearing premise

A modest number of Fourier harmonics plus the UGKS flux coupling will faithfully reproduce the full nonlinear time-periodic dynamics without significant truncation error or loss of stability, particularly when temperature modulations are large.

What would settle it

A side-by-side run of HB-UGKS (with a fixed small number of harmonics) against a fully resolved explicit time-marching UGKS for the thermally driven cavity at high Strouhal number and moderate Knudsen number; if the velocity or temperature fields differ by more than a few percent at corresponding instants, the truncation premise fails.

Figures

Figures reproduced from arXiv: 2605.03479 by Hangkong Wu, Kun Xu, Yufeng Wei, Yuze Zhu.

Figure 1
Figure 1. Figure 1: Flowchart of the pseudo-time marching procedure for the HB-UGKS solver. 3. Results and discussion In this section, the HB-UGKS is assessed for two representative oscillatory cavity prob￾lems. The shear-driven case is used to verify the linear response, the high-frequency limit, 11 view at source ↗
Figure 2
Figure 2. Figure 2: Harmonic independence verification: ℜ(Ub1/Ua) and ℑ(Ub1/Ua) along x = 0.5L for NH = 1 and NH = 2 (Kn = 0.1, St = 2). τexy = view at source ↗
Figure 3
Figure 3. Figure 3: Normalized fundamental amplitudes along the top lid at view at source ↗
Figure 4
Figure 4. Figure 4: Spatially averaged normalized amplitudes vs view at source ↗
Figure 5
Figure 5. Figure 5: Normalized density variation (ρn − ρ0)/(ρ0Ua) along the top lid at different phase angles for Kn = 0.1, A = 0.5. 17 view at source ↗
Figure 6
Figure 6. Figure 6: Normalized fundamental harmonic amplitudes of horizontal velocity and shear stress at view at source ↗
Figure 7
Figure 7. Figure 7: Normalized fundamental harmonic amplitudes of horizontal velocity and shear stress at view at source ↗
Figure 8
Figure 8. Figure 8: All-Kn evolution of fundamental harmonic amplitudes at St = 2, A = 1 for different Knudsen numbers. Finally, to demonstrate that the HB-UGKS inherits the multiscale asymptotic preserving (AP) property of the standard UGKS, it is essential to verify its capability to recover macro￾scopic hydrodynamic behaviors in the continuum limit. Although the Harmonic Balance method is inherently formulated for unsteady… view at source ↗
Figure 9
Figure 9. Figure 9: Continuum limit: Ub0/U w and Vb0/U w along the centerlines compared with the benchmark data of Ghia et al. 3.1.4. Evaluation of computational efficiency of the HB-UGKS The harmonic balance equation system is formulated not only to obtain the periodic solution accurately but also to significantly reduce the computational cost. To ensure a rigorous and fair comparison, both the calculations of the HB-UGKS an… view at source ↗
Figure 10
Figure 10. Figure 10: Convergence comparison between HB-UGKS and explicit time-domain UGKS for the shear view at source ↗
Figure 11
Figure 11. Figure 11: Harmonic independence verification: spatial distribution of the vertical velocity view at source ↗
Figure 12
Figure 12. Figure 12: Dynamic temperature profiles along the vertical centerline ( view at source ↗
Figure 13
Figure 13. Figure 13: Vertical velocity (V ) profiles along the cavity centerline (x = 0.5L) at different phase angles for various Knudsen numbers. Markers denote the reference data of Yang et al. [43]. 25 view at source ↗
Figure 14
Figure 14. Figure 14: Phase-resolved contours of the vertical velocity view at source ↗
Figure 15
Figure 15. Figure 15: Temporal evolution of temperature and vertical velocity at the probe point ( view at source ↗
Figure 16
Figure 16. Figure 16: Convergence comparison between HB-UGKS and explicit time-domain UGKS for the thermally view at source ↗
read the original abstract

This paper introduces a time-domain harmonic balance unified gas-kinetic scheme (HB-UGKS) designed to simulate temporally periodic flows across all Knudsen regimes. The harmonic balance approach reformulates the periodic problem into a block-coupled, quasi-steady system via a time-spectral source term. This allows for pseudo-time marching, local time-stepping, and the concurrent resolution of all sub-time levels, drastically reducing wall-clock time. Coupled with the UGKS-which maintains essential transport-collision coupling in its flux evaluations--the framework ensures multiscale validity across the entire Knudsen number range. The method is validated against two representative cavity flows. For a shear-driven oscillatory cavity under small-amplitude excitation, the fundamental harmonic alone accurately resolves the flow dynamics across various Knudsen and Strouhal numbers, successfully capturing the anti-resonance phenomenon and matching hydrodynamic damping predictions from linearized Boltzmann analyses. For a thermally driven cavity with large temperature modulations, higher-order harmonics prove essential to capture strong nonlinear waveform distortions and rarefaction effects. Beyond its physical fidelity, the HB-UGKS demonstrates substantial computational efficiency over explicit time-domain methods. This advantage peaks in high-frequency regimes, achieving speedup factors of 9.0 and 8.26 for the shear-driven and thermally driven cases, respectively.

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 introduces a time-domain harmonic balance unified gas-kinetic scheme (HB-UGKS) for simulating temporally periodic flows across all Knudsen regimes. The harmonic balance reformulation converts the periodic problem into a block-coupled quasi-steady system via a time-spectral source term, enabling pseudo-time marching, local time-stepping, and concurrent resolution of all sub-time levels. Coupled with UGKS flux evaluations that preserve transport-collision coupling, the method is validated on two cavity flows: a shear-driven oscillatory cavity (small-amplitude excitation, where the fundamental harmonic suffices to capture anti-resonance and hydrodynamic damping) and a thermally driven cavity (large temperature modulations, where higher-order harmonics are required for nonlinear waveform distortions and rarefaction effects). The paper reports substantial efficiency gains over explicit time-domain methods, with speedups peaking at 9.0 and 8.26 in high-frequency regimes for the two cases.

Significance. If the accuracy and stability claims are substantiated, the HB-UGKS would represent a useful advance for efficient simulation of periodic rarefied gas flows, extending the multiscale capabilities of UGKS to time-periodic problems without introducing free parameters or ad-hoc fitting. The approach leverages established UGKS properties and standard Fourier harmonic balance, and the reported matching of known hydrodynamic behaviors (anti-resonance, damping) in the shear-driven case provides a solid baseline; the efficiency advantage in high-frequency regimes could be particularly impactful for engineering applications in micro/nano flows if truncation errors are controlled.

major comments (2)
  1. [Validation (thermally driven cavity)] Validation section (thermally driven cavity): The claim that higher-order harmonics are essential to capture strong nonlinear waveform distortions and rarefaction effects, and that the method thereby reproduces the full nonlinear time-periodic dynamics across Knudsen regimes, lacks quantitative support. No L2 or pointwise error norms versus harmonic count, no convergence study against a long-time explicit UGKS reference solution, and no stability analysis of the pseudo-time marching for the block-coupled system are reported. This is load-bearing for the central efficiency claim, as the reported speedup of 8.26 relies on the assumption that a modest number of harmonics plus UGKS flux coupling faithfully reproduces the dynamics without significant truncation error for large temperature modulations.
  2. [Abstract and Validation] Abstract and validation (both cases): No grid-convergence data, error bars, or direct quantitative comparisons (e.g., pointwise or integrated errors) against full time-domain UGKS runs are provided, even for the shear-driven case where the fundamental harmonic is stated to accurately resolve dynamics. This weakens the justification for the speedup factors (9.0 and 8.26) and the assertion of multiscale validity across all Knudsen numbers.
minor comments (2)
  1. [Abstract] The abstract states that the method 'successfully capturing the anti-resonance phenomenon and matching hydrodynamic damping predictions from linearized Boltzmann analyses' but does not cite the specific analyses or provide direct numerical comparison tables/figures.
  2. [Methods] Notation for the time-spectral source term and the block-coupled system could be clarified with an explicit equation reference early in the methods section to aid readability for readers unfamiliar with harmonic balance formulations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We agree that the validation sections would benefit from additional quantitative metrics, error norms, and convergence studies to better support the efficiency claims and multiscale accuracy. We will revise the manuscript to incorporate these elements as outlined in our point-by-point responses below.

read point-by-point responses
  1. Referee: Validation section (thermally driven cavity): The claim that higher-order harmonics are essential to capture strong nonlinear waveform distortions and rarefaction effects, and that the method thereby reproduces the full nonlinear time-periodic dynamics across Knudsen regimes, lacks quantitative support. No L2 or pointwise error norms versus harmonic count, no convergence study against a long-time explicit UGKS reference solution, and no stability analysis of the pseudo-time marching for the block-coupled system are reported. This is load-bearing for the central efficiency claim, as the reported speedup of 8.26 relies on the assumption that a modest number of harmonics plus UGKS flux coupling faithfully reproduces the dynamics without significant truncation error for large temperature modulations.

    Authors: We agree that quantitative support is necessary to substantiate the claims for the thermally driven case. In the revised manuscript, we will add L2 and pointwise error norms as functions of harmonic count, computed against a long-time explicit time-domain UGKS reference solution. We will also include a convergence study with respect to the number of harmonics and a discussion of the stability and convergence behavior of the pseudo-time marching for the block-coupled system. These additions will directly address the truncation error concerns and strengthen the justification for the reported speedup of 8.26. revision: yes

  2. Referee: Abstract and validation (both cases): No grid-convergence data, error bars, or direct quantitative comparisons (e.g., pointwise or integrated errors) against full time-domain UGKS runs are provided, even for the shear-driven case where the fundamental harmonic is stated to accurately resolve dynamics. This weakens the justification for the speedup factors (9.0 and 8.26) and the assertion of multiscale validity across all Knudsen numbers.

    Authors: We acknowledge that the current manuscript lacks explicit grid-convergence data, error bars, and direct quantitative error comparisons against time-domain UGKS. In the revision, we will report grid-convergence studies for both cavity flows, including error norms on successively refined meshes. We will also provide direct L2 and pointwise error comparisons between HB-UGKS and full time-domain UGKS results for representative Knudsen and Strouhal numbers in both cases. Where relevant, we will include error estimates to support the speedup factors and the multiscale validity claims. revision: yes

Circularity Check

0 steps flagged

No circularity: HB-UGKS combines established UGKS flux and standard harmonic balance; performance metrics are measured outcomes, not self-referential

full rationale

The paper formulates HB-UGKS by applying a time-spectral source term to convert the periodic problem into a block-coupled quasi-steady system, then marches it with the existing UGKS flux evaluator. Validation is performed by direct comparison to explicit time-domain UGKS runs and to linearized Boltzmann results on two cavity flows; the reported speedups (9.0 and 8.26) are wall-clock timings from those runs. No parameter is fitted inside the same equations used to define the result, no uniqueness theorem is invoked from the authors' prior work to force the formulation, and the central claims rest on external benchmarks rather than on a closed loop of self-citation or self-definition. The modest number of harmonics is an explicit modeling choice whose truncation error is assessed by the numerical experiments themselves.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The approach relies on standard Fourier representation of periodic functions and the pre-existing transport-collision coupling property of UGKS; no new free parameters, ad-hoc constants, or postulated entities are introduced in the abstract.

axioms (2)
  • domain assumption Temporally periodic flows admit a convergent Fourier series representation with a finite number of harmonics sufficient for engineering accuracy.
    Invoked when the time-dependent problem is recast as a block-coupled quasi-steady system via the time-spectral source term.
  • domain assumption The unified gas-kinetic scheme preserves essential transport-collision coupling in its interface fluxes across the full Knudsen range.
    Stated as the reason the combined HB-UGKS remains valid from continuum to free-molecular regimes.

pith-pipeline@v0.9.0 · 5541 in / 1551 out tokens · 78382 ms · 2026-05-07T14:20:47.334654+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

48 extracted references · 23 canonical work pages

  1. [1]

    S. J. Lighthill, Acoustic streaming, Journal of Sound and Vibration 61 (3) (1978) 391–

  2. [2]

    doi:https://doi.org/10.1016/0022-460X(78)90388-7

  3. [3]

    Ho, Y.-C

    C.-M. Ho, Y.-C. Tai, Micro-electro-mechanical-systems (mems) and fluid flows, Annual Review of Fluid Mechanics 30 (Volume 30, 1998) (1998) 579–612

  4. [4]

    C.Cercignani, Rarefiedgasdynamics: frombasicconceptstoactualcalculations, Vol.21, Cambridge university press, 2000

  5. [5]

    Sone, Molecular gas dynamics: theory, techniques, and applications, Springer, 2007

    Y. Sone, Molecular gas dynamics: theory, techniques, and applications, Springer, 2007

  6. [6]

    Beskok, G

    A. Beskok, G. E. Karniadakis, Report: A model for flows in channels, pipes, and ducts at micro and nano scales, Microscale Thermophysical Engineering 3 (1) (1999) 43–77. doi:10.1080/108939599199864

  7. [7]

    N. G. Hadjiconstantinou, Oscillatory shear-driven gas flows in the transition and free- molecular-flow regimes, Physics of Fluids 17 (10) (2005) 100611. doi:10.1063/1.1874193

  8. [8]

    P. L. Bhatnagar, E. P. Gross, M. Krook, A model for collision processes in gases. i. small amplitude processes in charged and neutral one-component systems, Physical review 94 (3) (1954) 511

  9. [9]

    Bird, Approach to translational equilibrium in a rigid sphere gas, Phys

    G. Bird, Approach to translational equilibrium in a rigid sphere gas, Phys. fluids 6 (1963) 1518–1519

  10. [10]

    G. A. Bird, Direct simulation and the boltzmann equation, Physics of Fluids 13 (11) (1970) 2676–2681

  11. [11]

    G. A. Bird, Molecular gas dynamics and the direct simulation of gas flows, Oxford university press, 1994

  12. [12]

    N. G. Hadjiconstantinou, The limits of navier-stokes theory and kinetic extensions for describing small-scale gaseous hydrodynamics, Physics of Fluids 18 (11) (2006) 111301. doi:10.1063/1.2393436

  13. [13]

    J. Fan, C. Shen, Statistical simulation of low-speed rarefied gas flows, Journal of Com- putational Physics 167 (2) (2001) 393–412. doi:https://doi.org/10.1006/jcph.2000.6681

  14. [14]

    T. M. Homolle, N. G. Hadjiconstantinou, A low-variance deviational simulation monte carlo for the boltzmann equation, Journal of Computational Physics 226 (2) (2007) 2341–2358. doi:https://doi.org/10.1016/j.jcp.2007.07.006

  15. [15]

    Z. Cai, R. Li, Numerical regularized moment method of arbitrary order for boltzmann- bgk equation, SIAM Journal on Scientific Computing 32 (5) (2010) 2875–2907. doi:10.1137/100785466

  16. [16]

    Chu, Kinetic-theoretic description of the formation of a shock wave, The Physics of Fluids 8 (1) (1965) 12–22

    C. Chu, Kinetic-theoretic description of the formation of a shock wave, The Physics of Fluids 8 (1) (1965) 12–22. 30

  17. [17]

    J. E. Broadwell, Study of rarefied shear flow by the discrete velocity method, Journal of Fluid Mechanics 19 (3) (1964) 401–414. doi:10.1017/S0022112064000817

  18. [18]

    L. Mieussens, Discrete-velocity models and numerical schemes for the boltzmann-bgk equation in plane and axisymmetric geometries, Journal of Computational Physics 162 (2) (2000) 429–466. doi:https://doi.org/10.1006/jcph.2000.6548

  19. [19]

    L. Wu, J. Zhang, J. M. Reese, Y. Zhang, A fast spectral method for the boltzmann equation for monatomic gas mixtures, Journal of Computational Physics 298 (2015) 602–621. doi:https://doi.org/10.1016/j.jcp.2015.06.019

  20. [20]

    L. Wu, A fast spectral method for the uehling-uhlenbeck equation for quantum gas mixtures: Homogeneous relaxation and transport coefficients, Journal of Computational Physics 399 (2019) 108924. doi:https://doi.org/10.1016/j.jcp.2019.108924

  21. [21]

    L. Zhu, X. Pi, W. Su, Z.-H. Li, Y. Zhang, L. Wu, General synthetic iterative scheme for nonlinear gas kinetic simulation of multi-scale rarefied gas flows, Journal of Computa- tional Physics 430 (2021) 110091. doi:https://doi.org/10.1016/j.jcp.2020.110091

  22. [22]

    Xu, J.-C

    K. Xu, J.-C. Huang, A unified gas-kinetic scheme for continuum and rarefied flows, Journal of Computational Physics 229 (20) (2010) 7747–7764

  23. [23]

    Huang, K

    J.-C. Huang, K. Xu, P. Yu, A unified gas-kinetic scheme for continuum and rarefied flows ii: multi-dimensional cases, Communications in Computational Physics 12 (3) (2012) 662–690

  24. [24]

    Z. Guo, K. Xu, Progress of discrete unified gas-kinetic scheme for multiscale flows, Advances in Aerodynamics 3 (1) (2021) 6

  25. [25]

    C. Liu, K. Xu, Q. Sun, Q. Cai, A unified gas-kinetic scheme for continuum and rar- efied flows v: Multiscale and multi-component plasma transport, Communications in Computational Physics 22 (5) (2017) 1175–1223. doi:10.4208/cicp.OA-2017-0056

  26. [26]

    T. Xiao, Q. Cai, K. Xu, A well-balanced unified gas-kinetic scheme for multiscale flow transport under gravitational field, Journal of Computational Physics 332 (2017) 475–

  27. [27]

    doi:https://doi.org/10.1016/j.jcp.2016.12.022

  28. [28]

    W. Sun, S. Jiang, K. Xu, S. Li, An asymptotic preserving unified gas kinetic scheme for frequency-dependent radiative transfer equations, Journal of Computational Physics 302 (2015) 222–238. doi:https://doi.org/10.1016/j.jcp.2015.09.002

  29. [29]

    C. Liu, Z. Wang, K. Xu, A unified gas-kinetic scheme for continuum and rarefied flows vi: Dilute disperse gas-particle multiphase system, Journal of Computational Physics 386 (2019) 264–295. doi:https://doi.org/10.1016/j.jcp.2018.12.040

  30. [30]

    P. Du, F. Ning, Simulating periodic unsteady flows using a cubic-spline-based time collocation method, Journal of Turbomachinery 136 (4) (2013) 041014. doi:10.1115/1.4025203. 31

  31. [31]

    Blazek (Ed.), Computational Fluid Dynamics: Principles and Applications (Second Edition), second edition Edition, Elsevier Science, Oxford, 2005, pp

    List of symbols, in: J. Blazek (Ed.), Computational Fluid Dynamics: Principles and Applications (Second Edition), second edition Edition, Elsevier Science, Oxford, 2005, pp. xiii–xvii. doi:https://doi.org/10.1016/B978-008044506-9/50001-1

  32. [32]

    Y. Zhu, C. Zhong, K. Xu, Unified gas-kinetic scheme with multigrid convergence for rarefied flow study, Physics of Fluids 29 (9) (2017) 096102. doi:10.1063/1.4994020

  33. [33]

    Y. Zhu, C. Zhong, K. Xu, An implicit unified gas-kinetic scheme for unsteady flow in all knudsen regimes, Journal of Computational Physics 386 (2019) 190–217. doi:https://doi.org/10.1016/j.jcp.2019.01.033

  34. [34]

    S. Tan, W. Sun, K. Xu, J. Wei, G. Ni, Time implicit unified gas kinetic scheme for 3d multi-group neutron transport simulation, Communications in Computational Physics 28 (3) (2020) 1189–1218

  35. [35]

    X. Xu, Y. Zhu, C. Liu, K. Xu, Ugks-based implicit iterative method for multiscale nonequilibrium flow simulations, SIAM Journal on Scientific Computing 44 (4) (2022) B996–B1017

  36. [36]

    K. C. Hall, W. S. Clark, Linearized euler predictions of unsteady aerodynamic loads in cascades, AIAA journal 31 (3) (1993) 540–550

  37. [37]

    W. S. Clark, K. C. Hall, A time-linearized navier–stokes analysis of stall flutter, J. Turbomach. 122 (3) (2000) 467–476

  38. [38]

    L. Wu, J. M. Reese, Y. Zhang, Oscillatory rarefied gas flow inside rectangular cavities, Journal of Fluid Mechanics 748 (2014) 350–367

  39. [39]

    K. C. Hall, J. P. Thomas, W. S. Clark, Computation of unsteady nonlinear flows in cascades using a harmonic balance technique, AIAA journal 40 (5) (2002) 879–886

  40. [40]

    K. C. Hall, K. Ekici, J. P. Thomas, E. H. Dowell, Harmonic balance methods applied to computational fluid dynamics problems, International Journal of Computational Fluid Dynamics 27 (2) (2013) 52–67

  41. [41]

    Zhang, D

    Y. Zhang, D. Wang, S. Zhang, Y. Zhu, Numerical stability analysis of implicit solution methods for harmonic balance equations, Journal of Turbomachinery 147 (12) (2025) 121010. doi:10.1115/1.4068947

  42. [42]

    Y. Zhu, Y. Zhang, S. Zhang, D. Wang, Investigation of the sigma approximation tech- nique for the solution of the time spectral equation system, Journal of Engineering for Gas Turbines and Power 147 (1) (2024) 011027. doi:10.1115/1.4066586

  43. [43]

    Huang, D

    X. Huang, D. X. Wang, Time-space spectral method for rotor–rotor/stator–stator inter- actions, Journal of Turbomachinery 141 (11) (2019) 111006

  44. [44]

    Huang, H

    X. Huang, H. Wu, D. Wang, Implicit solution of harmonic balance equation system using the lu-sgs method and one-step jacobi/gauss-seidel iteration, International Journal of Computational Fluid Dynamics 32 (4-5) (2018) 218–232. 32

  45. [45]

    W. Yang, Y. Niu, Oscillatory thermally driven rarefied flows using a coupled macro- scopic/mesoscopic method, Journal of National University of Defense Technology 47 (6) (2025) 168–177

  46. [46]

    Shakhov, Generalization of the krook kinetic relaxation equation, Fluid dynamics 3 (5) (1968) 95–96

    E. Shakhov, Generalization of the krook kinetic relaxation equation, Fluid dynamics 3 (5) (1968) 95–96

  47. [47]

    S. Chen, K. Xu, A comparative study of an asymptotic preserving scheme and unified gas-kinetic scheme in continuum flow limit, Journal of computational physics 288 (2015) 52–65. doi:http://dx.doi.org/10.1016/j.jcp.2015.02.014

  48. [48]

    U. Ghia, K. N. Ghia, C. Shin, High-re solutions for incompressible flow using the navier- stokes equations and a multigrid method, Journal of computational physics 48 (3) (1982) 387–411. 33