Recognition: unknown
A Time-Domain Harmonic Balance Unified Gas-Kinetic Scheme for Temporally Periodic Flows Across all Knudsen Regimes
Pith reviewed 2026-05-07 14:20 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [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.
- [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
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
-
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
-
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
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
axioms (2)
- domain assumption Temporally periodic flows admit a convergent Fourier series representation with a finite number of harmonics sufficient for engineering accuracy.
- domain assumption The unified gas-kinetic scheme preserves essential transport-collision coupling in its interface fluxes across the full Knudsen range.
Reference graph
Works this paper leans on
-
[1]
S. J. Lighthill, Acoustic streaming, Journal of Sound and Vibration 61 (3) (1978) 391–
1978
-
[2]
doi:https://doi.org/10.1016/0022-460X(78)90388-7
-
[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
1998
-
[4]
C.Cercignani, Rarefiedgasdynamics: frombasicconceptstoactualcalculations, Vol.21, Cambridge university press, 2000
2000
-
[5]
Sone, Molecular gas dynamics: theory, techniques, and applications, Springer, 2007
Y. Sone, Molecular gas dynamics: theory, techniques, and applications, Springer, 2007
2007
-
[6]
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]
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]
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
1954
-
[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
1963
-
[10]
G. A. Bird, Direct simulation and the boltzmann equation, Physics of Fluids 13 (11) (1970) 2676–2681
1970
-
[11]
G. A. Bird, Molecular gas dynamics and the direct simulation of gas flows, Oxford university press, 1994
1994
-
[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]
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]
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]
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]
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
1965
-
[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]
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]
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]
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]
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]
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
2010
-
[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
2012
-
[24]
Z. Guo, K. Xu, Progress of discrete unified gas-kinetic scheme for multiscale flows, Advances in Aerodynamics 3 (1) (2021) 6
2021
-
[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]
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–
2017
-
[27]
doi:https://doi.org/10.1016/j.jcp.2016.12.022
-
[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]
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]
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]
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]
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]
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]
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
2020
-
[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
2022
-
[36]
K. C. Hall, W. S. Clark, Linearized euler predictions of unsteady aerodynamic loads in cascades, AIAA journal 31 (3) (1993) 540–550
1993
-
[37]
W. S. Clark, K. C. Hall, A time-linearized navier–stokes analysis of stall flutter, J. Turbomach. 122 (3) (2000) 467–476
2000
-
[38]
L. Wu, J. M. Reese, Y. Zhang, Oscillatory rarefied gas flow inside rectangular cavities, Journal of Fluid Mechanics 748 (2014) 350–367
2014
-
[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
2002
-
[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
2013
-
[41]
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]
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]
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
2019
-
[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
2018
-
[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
2025
-
[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
1968
-
[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]
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
1982
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.