pith. sign in

arxiv: 2604.26463 · v5 · submitted 2026-04-29 · ⚛️ physics.flu-dyn

A Hybrid Gas-Kinetic Scheme and Discrete Velocity Method for Continuum and Rarefied Flows

Hangkong Wu , Yuze Zhu , Yajun Zhu , Kun Xu This is my paper

Pith reviewed 2026-05-08 03:11 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords hybrid methodgas-kinetic schemediscrete velocity methodcontinuum and rarefied flowsnumerical collision timeasymptotic preservingshock capturingadaptive strategies
0
0 comments X

The pith

A hybrid gas-kinetic and discrete velocity scheme recovers Navier-Stokes solutions in dense flows while capturing free molecular behavior in rarefied regimes.

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

The paper develops a hybrid method that combines the efficiency of gas-kinetic schemes for continuum flows with the accuracy of discrete velocity methods for rarefied flows. It achieves this by balancing equilibrium distributions from one approach with upwind-reconstructed non-equilibrium distributions from the other, using a numerical collision time as the bridge. Adaptive partitioning based on local Knudsen and Mach numbers further cuts computational cost without sacrificing fidelity. A sympathetic reader would care because many engineering flows, such as those around vehicles or in vacuum systems, span both regimes and currently require switching between separate solvers. Systematic tests on boundary layers, cavities, shocks, and cylinders confirm the hybrid recovers expected limits in each regime.

Core claim

The hybrid approach balances the equilibrium distribution function in GKS with the upwind-reconstructed non-equilibrium distribution function in DVM through a numerical collision time. This balancing strategy ensures to recover Navier-Stokes solutions in the continuum limit (asymptotic preserving), while naturally capturing free molecular flows in the rarefied limit. The numerical collision time also improves robustness for shock capturing, and adaptive domain partitioning based on local Knudsen and Mach numbers reduces overall cost.

What carries the argument

The numerical collision time that blends equilibrium distributions from the gas-kinetic scheme with non-equilibrium distributions from the discrete velocity method.

If this is right

  • In continuum limits the scheme reduces exactly to Navier-Stokes solutions.
  • In rarefied limits it reproduces free-molecular behavior without additional modifications.
  • Shock waves are captured more robustly than in pure GKS due to the added collision-time control.
  • Adaptive switching based on local Knudsen and Mach numbers automatically assigns the cheaper method where appropriate.

Where Pith is reading between the lines

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

  • The same blending idea could be applied to other pairs of kinetic and macroscopic schemes for multi-scale transport problems.
  • Extension to unsteady flows would require checking whether the collision-time definition remains stable during transients.
  • The adaptive criterion might be refined by including higher moments to detect non-equilibrium more precisely.

Load-bearing premise

A single tunable numerical collision time can blend the two distribution functions without creating instability or accuracy loss in any flow regime.

What would settle it

A simulation of a known shock structure or free-molecular flow where the hybrid results deviate systematically from established reference solutions or analytical limits would show the balancing fails.

Figures

Figures reproduced from arXiv: 2604.26463 by Hangkong Wu, Kun Xu, Yajun Zhu, Yuze Zhu.

Figure 1
Figure 1. Figure 1: The computational mesh of the flat plate boundary layer view at source ↗
Figure 2
Figure 2. Figure 2: Comparison between the numerical results and the Blasius solutions for the flat view at source ↗
Figure 3
Figure 3. Figure 3: The weighting factors in the whole computational domain view at source ↗
Figure 4
Figure 4. Figure 4: Comparison of wall-clock time between the hybrid approaches with and without view at source ↗
Figure 5
Figure 5. Figure 5: The computational mesh of the lid-driven cavity flow view at source ↗
Figure 6
Figure 6. Figure 6: The velocity contours in the whole computational domain (UGKS: the colored view at source ↗
Figure 7
Figure 7. Figure 7: The density residual at Kn = 10: a) convergence histories; b) wall-clock time (a) U (b) V view at source ↗
Figure 8
Figure 8. Figure 8: The velocity contours in the whole computational domain (UGKS: the colored view at source ↗
Figure 9
Figure 9. Figure 9: The density residual at Kn = 0.075: a) convergence histories; b) wall-clock time (a) U (b) V view at source ↗
Figure 10
Figure 10. Figure 10: The velocity profiles along the central line at view at source ↗
Figure 11
Figure 11. Figure 11: The streamlines with a mesh resolution of 61 view at source ↗
Figure 12
Figure 12. Figure 12: The axial distribution of density and temperature for the shock-structure case at view at source ↗
Figure 13
Figure 13. Figure 13: The axial distribution of the weighting factors and comparison of computational view at source ↗
Figure 14
Figure 14. Figure 14: The axial distribution of density and temperature for the shock-structure case at view at source ↗
Figure 15
Figure 15. Figure 15: The axial distribution of the weighting factors and comparison of computational view at source ↗
Figure 16
Figure 16. Figure 16: The axial distribution of density and temperature for the shock-structure case at view at source ↗
Figure 17
Figure 17. Figure 17: The axial distribution of the weighting factors and comparison of computational view at source ↗
Figure 18
Figure 18. Figure 18: The computational mesh of the flow past a semi-cylinder view at source ↗
Figure 19
Figure 19. Figure 19: Comparison of density, axial velocity, temperature contours at three different view at source ↗
Figure 20
Figure 20. Figure 20: Comparison of density, axial velocity, temperature along the axial direction at view at source ↗
Figure 21
Figure 21. Figure 21: Comparison of the residual convergence histories and wall-clock time between the view at source ↗
read the original abstract

The gas-kinetic scheme (GKS) provides high computational efficiency and accuracy for continuum flow simulations but is unable to reliably capture rarefaction effects. In contrast, although the discrete velocity method (DVM) is better suited for rarefied flows, it exhibits reduced accuracy and slow convergence when applied to continuum regimes. To overcome these limitations, this work proposes a hybrid GKS-DVM method that integrates the strengths of both approaches. The hybrid approach balances the equilibrium distribution function in GKS with the upwind-reconstructed non-equilibrium distribution function in DVM through a numerical collision time. This balancing strategy ensures to recover Navier-Stokes solutions in the continuum limit (asymptotic preserving), while naturally capturing free molecular flows in the rarefied limit. Moreover, the introduction of a numerical collision time significantly enhances robustness in shock capturing for continuum flow applications. To further reduce computational cost of the hybrid approach, several adaptive strategies based on the local Knudsen number and Mach number have been proposed. The effectiveness and accuracy of the proposed hybrid method are systematically assessed through four representative test cases: a flat-plate boundary layer, a lid-driven cavity flow, shock structures, and flow past a semi-cylinder. The first case is subjected to continuum conditions, while the latter two span a broad range of Knudsen numbers. The results demonstrate that the proposed method achieves high solution accuracy and computational efficiency across both continuum and rarefied flow regimes.

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

3 major / 3 minor

Summary. The manuscript proposes a hybrid gas-kinetic scheme (GKS) and discrete velocity method (DVM) for continuum and rarefied flows. The hybrid balances the GKS equilibrium distribution function with the DVM upwind-reconstructed non-equilibrium distribution via a numerical collision time, augmented by adaptive domain partitioning based on local Knudsen and Mach numbers. The approach is asserted to be asymptotic-preserving (recovering Navier-Stokes in the Kn→0 limit) while capturing free-molecular behavior at high Kn, with improved shock robustness. Validation is presented on four cases: flat-plate boundary layer (continuum), lid-driven cavity, shock structures, and semi-cylinder flow (spanning Knudsen numbers).

Significance. If the central construction holds, the hybrid offers a practical route to multi-regime kinetic simulations by retaining GKS efficiency in continuum regions and DVM fidelity in rarefied regions. The numerical collision time for shock robustness and the Kn/Ma-based adaptivity are potentially useful engineering features. The test suite covers relevant regimes, but quantitative error metrics and convergence data would strengthen the significance assessment.

major comments (3)
  1. [§3.2] §3.2 (numerical collision time definition): the balancing of equilibrium (GKS) and non-equilibrium (DVM) parts through τ_num is load-bearing for both the asymptotic-preserving claim and the rarefied-capture claim, yet the manuscript does not derive the resulting moments or show that the hybrid distribution recovers the correct Chapman-Enskog expansion without residual errors in the continuum limit.
  2. [§4.3] §4.3 (shock-structure results): the paper reports qualitative agreement but supplies no L2 or L∞ error norms against reference Boltzmann or DSMC solutions across the tested Knudsen range; without these, the claim that the method “naturally captures” rarefied behavior cannot be assessed quantitatively.
  3. [§3.4] §3.4 (adaptive Kn/Ma partitioning): the switching thresholds and interface treatment between GKS and DVM regions are not accompanied by a sensitivity study or proof of solution continuity; small variations in the local Kn/Ma criteria could alter the hybrid interface and undermine the reported efficiency gains.
minor comments (3)
  1. [§2] Notation for the equilibrium and non-equilibrium distribution functions is introduced inconsistently between the abstract and §2; a single, early definition would improve readability.
  2. Figure captions for the semi-cylinder and cavity cases omit the specific Knudsen and Mach numbers used in each panel, forcing the reader to cross-reference the text.
  3. [§1] The literature review omits several recent hybrid kinetic schemes (e.g., those combining BGK and DVM with similar adaptivity); adding these would better situate the novelty.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive comments. We address each major comment below and indicate planned revisions.

read point-by-point responses

  1. Referee: [§3.2] §3.2 (numerical collision time definition): the balancing of equilibrium (GKS) and non-equilibrium (DVM) parts through τ_num is load-bearing for both the asymptotic-preserving claim and the rarefied-capture claim, yet the manuscript does not derive the resulting moments or show that the hybrid distribution recovers the correct Chapman-Enskog expansion without residual errors in the continuum limit.

    Authors: We acknowledge the value of an explicit moment derivation. The numerical collision time is constructed so that small local Kn forces the hybrid distribution to the GKS equilibrium (known to recover NS via Chapman-Enskog), while large Kn recovers the DVM non-equilibrium part. We will add a dedicated derivation of the hybrid moments in the revised manuscript to confirm absence of residual errors in the continuum limit. revision: yes

  2. Referee: [§4.3] §4.3 (shock-structure results): the paper reports qualitative agreement but supplies no L2 or L∞ error norms against reference Boltzmann or DSMC solutions across the tested Knudsen range; without these, the claim that the method “naturally captures” rarefied behavior cannot be assessed quantitatively.

    Authors: We agree that quantitative norms strengthen the validation. In the revision we will compute and tabulate L2 and L∞ errors for the shock-structure cases against available Boltzmann/DSMC references over the tested Kn range. revision: yes

  3. Referee: [§3.4] §3.4 (adaptive Kn/Ma partitioning): the switching thresholds and interface treatment between GKS and DVM regions are not accompanied by a sensitivity study or proof of solution continuity; small variations in the local Kn/Ma criteria could alter the hybrid interface and undermine the reported efficiency gains.

    Authors: The Kn/Ma criteria are chosen from local flow physics to maintain accuracy while improving efficiency. We will add a sensitivity study showing the effect of modest threshold variations on both solution fields and CPU time. The interface blending is formulated to ensure continuous distribution functions; while a general analytic proof is difficult for adaptive partitions, the numerical results already demonstrate smooth transitions without spurious oscillations. revision: partial

Circularity Check

0 steps flagged

No significant circularity; hybrid construction is self-contained synthesis of established schemes

full rationale

The paper constructs a hybrid GKS-DVM scheme by introducing a numerical collision time to balance equilibrium (GKS) and upwind non-equilibrium (DVM) distribution functions, with adaptive Kn/Ma partitioning. The asymptotic-preserving property in the continuum limit follows directly from reducing to the established GKS (which recovers NS) when the collision time is appropriately scaled, and free-molecular recovery at high Kn follows from DVM dominance; neither is a fitted prediction renamed as output nor a self-referential definition. No load-bearing step reduces by the paper's equations to its own inputs. Validation occurs via external test cases (boundary layer, cavity, shocks, semi-cylinder) rather than internal consistency loops. Minor self-citations to prior GKS/DVM work exist but are not invoked as uniqueness theorems or ansatzes that close the central claim.

Axiom & Free-Parameter Ledger

1 free parameters · 0 axioms · 0 invented entities

The central claim rests on the effectiveness of an introduced numerical collision time whose precise definition, functional form, and tuning procedure are not supplied in the abstract.

free parameters (1)
  • numerical collision time
    Used to balance equilibrium GKS and non-equilibrium DVM distributions; value and dependence on local flow variables are not detailed.

pith-pipeline@v0.9.0 · 5559 in / 1219 out tokens · 69219 ms · 2026-05-08T03:11:14.033074+00:00 · methodology

Review history (2 revisions) →

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

38 extracted references · 32 canonical work pages

  1. [1]

    R. C. Blanchard, R. G. Wilmoth, J. N. Moss, Aerodynamic flight mea- surements and rarefied-flow simulations of Mars entry vehicles, Journal of Spacecraft and Rockets 34 (5) (1997) 687–690.doi:10.2514/2.3271

    work page doi:10.2514/2.3271 1997

  2. [2]

    Roe, Approximate Riemann solvers, parameter vectors, and difference schemes, Journal of Computational Physics 135 (2) (1997) 250–258.doi: 10.1006/jcph.1997.5705

    P. Roe, Approximate Riemann solvers, parameter vectors, and difference schemes, Journal of Computational Physics 135 (2) (1997) 250–258.doi: 10.1006/jcph.1997.5705

    work page doi:10.1006/jcph.1997.5705 1997

  3. [3]

    K. Xu, K. H. Prendergast, Numerical Navier-Stokes solutions from gas kinetic theory, Journal of Computational Physics 114 (1) (1994) 9–17.doi: 10.1006/jcph.1994.1145. 30

    work page doi:10.1006/jcph.1994.1145 1994

  4. [4]

    K. Xu, A gas-kinetic BGK scheme for the Navier–Stokes equations and its connection with artificial dissipation and Godunov method, Journal of Computational Physics 171 (1) (2001) 289–335.doi:10.1006/jcph.2001. 6790

    work page doi:10.1006/jcph.2001 2001

  5. [5]

    P. Wang, M. T. Ho, L. Wu, Z. Guo, Y. Zhang, A comparative study of discrete velocity methods for low-speed rarefied gas flows, Computers & Fluids 161 (2018) 33–46.doi:10.1016/j.compfluid.2017.11.006

    work page doi:10.1016/j.compfluid.2017.11.006 2018

  6. [6]

    L. MIEUSSENS, Discrete velocity model and implicit scheme for the BGK equation of rarefied gas dynamics, Mathematical Models and Meth- ods in Applied Sciences 10 (08) (2000) 1121–1149.doi:10.1142/ S0218202500000562

    2000

  7. [7]

    R. Yuan, C. Zhong, A conservative implicit scheme for steady state solu- tions of diatomic gas flow in all flow regimes, Computer Physics Commu- nications 247 (2020) 106972.doi:10.1016/j.cpc.2019.106972

    work page doi:10.1016/j.cpc.2019.106972 2020

  8. [8]

    Xu, J.-C

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

    work page doi:10.1016/j.jcp.2010.06.032 2010

  9. [9]

    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 Com- putational Physics 12 (3) (2012) 662–690.doi:10.4208/cicp.030511. 220911a

    work page doi:10.4208/cicp.030511 2012

  10. [10]

    Huang, K

    J.-C. Huang, K. Xu, P. Yu, A unified gas-kinetic scheme for continuum and rarefied flows III: Microflow simulations, Communications in Com- putational Physics 14 (5) (2013) 1147–1173.doi:10.4208/cicp.190912. 080213a

    work page doi:10.4208/cicp.190912 2013

  11. [11]

    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, Phys. Rev. 94 (1954) 511–525.doi:10.1103/PhysRev.94.511. 31

    work page doi:10.1103/physrev.94.511 1954

  12. [12]

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

    1968

  13. [13]

    Xu, J.-C

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

    work page doi:10.1093/imamat/hxr002 2011

  14. [14]

    C. Liu, K. Xu, Q. Sun, Q. Cai, A unified gas-kinetic scheme for continuum and rarefied flows iv: Full boltzmann and model equations, Journal of Computational Physics 314 (2016) 305–340.doi:10.1016/j.jcp.2016. 03.014

    work page doi:10.1016/j.jcp.2016 2016

  15. [15]

    S. Liu, C. Zhong, Modified unified kinetic scheme for all flow regimes, Phys. Rev. E 85 (6) (2014) 066705.doi:10.1103/PhysRevE.85.066705

    work page doi:10.1103/physreve.85.066705 2014

  16. [16]

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

    work page doi:10.1016/j.jcp.2015.02.014 2015

  17. [17]

    S. Liu, P. Yu, K. Xu, C. Zhong, Unified gas-kinetic scheme for diatomic molecular simulations in all flow regimes, Journal of Computational Physics 259 (2014) 96–113.doi:10.1016/j.jcp.2013.11.030

    work page doi:10.1016/j.jcp.2013.11.030 2014

  18. [18]

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

    2021

  19. [19]

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

    work page doi:10.1016/j.jcp.2019.01.033 2019

  20. [20]

    Y. Zhu, C. Zhong, K. Xu, Implicit unified gas-kinetic scheme for steady state solutions in all flow regimes, Journal of Computational Physics 315 (2016) 16–38.doi:doi.org/10.1016/j.jcp.2016.03.038. 32

    work page doi:10.1016/j.jcp.2016.03.038 2016

  21. [21]

    W. Long, Y. Wei, K. Xu, An implicit adaptive unified gas-kinetic scheme for steady-state solutions of nonequilibrium flows, Physics of Fluids 36 (10) (2024) 106114.doi:10.1063/5.0232275

    work page doi:10.1063/5.0232275 2024

  22. [22]

    S. Chen, K. Xu, C. Lee, Q. Cai, A unified gas kinetic scheme with moving mesh and velocity space adaptation, Journal of Computational Physics 231 (20) (2012) 6643–6664.doi:https://doi.org/10.1016/j.jcp.2012. 05.019

    work page doi:10.1016/j.jcp.2012 2012

  23. [23]

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

    work page doi:10.1063/1.4994020 2017

  24. [24]

    Zhang, Y

    Y. Zhang, Y. Wei, W. Long, K. Xu, An efficiency and memory-saving pro- gramming paradigm for the unified gas-kinetic scheme, Computer Physics Communications 314 (2025) 109684.doi:doi.org/10.1016/j.cpc.2025. 109684

    work page doi:10.1016/j.cpc.2025 2025

  25. [25]

    Z. Guo, K. Xu, R. Wang, Discrete unified gas kinetic scheme for all knudsen number flows: low-speed isothermal case, Phys. Rev. E 88 (2013) 033305. doi:10.1103/PHYSREVE.88.033305

    work page doi:10.1103/physreve.88.033305 2013

  26. [26]

    B. Shan, P. Wang, Y. Zhang, Z. Guo, Discrete unified gas kinetic scheme for all knudsen number flows. IV. strongly inhomogeneous fluids, Phys. Rev. E 101 (2020) 043303.doi:10.1103/PhysRevE.101.043303

    work page doi:10.1103/physreve.101.043303 2020

  27. [27]

    S. Chen, Z. Guo, K. Xu, Simplification of the unified gas kinetic scheme, Phys. Rev. E 94 (2016) 023313.doi:10.1103/PhysRevE.94.023313

    work page doi:10.1103/physreve.94.023313 2016

  28. [28]

    Zhang, S

    R. Zhang, S. Liu, J. Chen, H. Jin, C. Zhuo, C. Zhong, Im- plicit unified gas-kinetic scheme for steady state solution of hy- personic thermodynamic non-equilibrium flows, Communications in Nonlinear Science and Numerical Simulation 140 (2025) 108367. doi:https://doi.org/10.1016/j.cnsns.2024.108367. 33 URLhttps://www.sciencedirect.com/science/article/pii/ ...

    work page doi:10.1016/j.cnsns.2024.108367 2025

  29. [30]

    W. Su, L. Zhu, L. Wu, Fast convergence and asymptotic preserving of the general synthetic iterative scheme, SIAM Journal on Scientific Computing 42 (6) (2020) B1517–B1540.doi:10.1137/20M132691X

    work page doi:10.1137/20m132691x 2020

  30. [31]

    L. Yang, C. Shu, W. Yang, Z. Chen, H. Dong, An improved discrete velocity method (DVM) for efficient simulation of flows in all flow regimes, Physics of Fluids 30 (6) (2018).doi:10.1063/1.5039479

    work page doi:10.1063/1.5039479 2018

  31. [32]

    L. Yang, C. Shu, W. Yang, J. Wu, An improved three-dimensional implicit discrete velocity method on unstructured meshes for all knudsen number flows, Journal of Computational Physics 396 (2019) 738–760.doi:10. 1016/j.jcp.2019.07.002

    2019

  32. [33]

    R. Yuan, S. Liu, C. Zhong, A novel multiscale discrete velocity method for model kinetic equations, Communications in Nonlinear Science and Nu- merical Simulation 92 (2021) 105473.doi:https://doi.org/10.1016/j. cnsns.2020.105473

    work page doi:10.1016/j 2021

  33. [34]

    Zhang, C

    Q. Zhang, C. Zhuo, J. Mu, C. Zhong, S. Liu, A multiscale discrete velocity method for diatomic molecular gas, Physics of Fluids 35 (7) (2023) 076118. doi:10.1063/5.0153834

    work page doi:10.1063/5.0153834 2023

  34. [35]

    V. Venkatakrishnan, Convergence to steady state solutions of the Euler equations on unstructured grids with limiters, Journal of Computational Physics 118 (1) (1995) 120–130.doi:doi.org/10.1006/jcph.1995.1084

    work page doi:10.1006/jcph.1995.1084 1995

  35. [36]

    X. Ji, W. Shyy, K. Xu, A gradient compression-based compact high-order gas-kinetic scheme on 3D hybrid unstructured meshes, International Jour- 34 nal of Computational Fluid Dynamics 35 (2021) 485–509.doi:doi.org/ 10.1080/10618562.2021.1991329

    work page doi:10.1080/10618562.2021.1991329 2021

  36. [37]

    Zhang, X

    Y. Zhang, X. Ji, K. Xu, A high-order compact gas-kinetic scheme in a rotating coordinates frame and on sliding mesh, International Journal of Computational Fluid Dynamics 37 (2022) 181 – 200.doi:doi:10.1080/ 10618562.2023.2178647

    work page arXiv 2022

  37. [38]

    John, X.-J

    B. John, X.-J. Gu, D. R. Emerson, Effects of incomplete surface accom- modation on non-equilibrium heat transfer in cavity flow: A parallel dsmc study, Computers & Fluids 45 (1) (2011) 197–201.doi:doi.org/10.1016/ j.compfluid.2011.01.044

    2011

  38. [39]

    U. Ghia, K. 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.doi:doi.org/10.1016/ 0021-9991(82)90058-4. 35

    1982