pith. sign in

arxiv: 2606.26915 · v1 · pith:XZOY73OVnew · submitted 2026-06-25 · 🧮 math-ph · math.MP

Wave-Particle Decomposition for Kinetic Equations I: Theory and Numerics

Chang Liu , Kun Xu This is my paper

Pith reviewed 2026-06-26 02:30 UTC · model grok-4.3

classification 🧮 math-ph math.MP
keywords wave-particle decompositionkinetic relaxation equationsKnudsen spectrumasymptotic preserving methodsChapman-Enskog expansionmacro-micro discretizationgas kinetic schemes
0
0 comments X

The pith

A wave-particle decomposition of the distribution function produces a unified kinetic system valid across the full Knudsen spectrum.

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

The paper formulates a continuous decomposition of the distribution function in relaxation-type kinetic equations by using the characteristic integral solution around a local evolution timescale. This splits the function into an analytically accumulated wave component inside a kinetic horizon and a surviving particle component outside it. The resulting system consists of a source-free total conservation law, a wave equation, and a particle equation. The wave operator recovers fluid equations through a horizon-dependent Chapman-Enskog expansion, while the particle equation carries the non-equilibrium transport. The decomposition supports a conservative macro-micro discretization whose active kinetic degrees of freedom adapt locally to the regime.

Core claim

Leveraging the characteristic integral solution of the kinetic relaxation equation, the distribution function is decomposed into a wave component that accumulates analytically within the local horizon and a particle component defined by the collisionless transport that survives beyond the horizon. This yields a closed wave-particle system comprising a source-free total conservation law, a wave equation whose moments produce Euler and Navier-Stokes fluxes with horizon-dependent coefficients, and a particle equation for the remaining non-equilibrium kinetic transport. The formulation holds across the entire Knudsen spectrum and is discretized by advancing total conservative variables with a fi

What carries the argument

The wave-particle decomposition defined around a local evolution timescale and its associated kinetic horizon, which uses the characteristic integral solution to separate the distribution function into an analytically accumulated wave component and a purely kinetic particle component.

Load-bearing premise

The decomposition assumes that a local evolution timescale and associated kinetic horizon can be defined such that the characteristic integral solution accurately separates the analytically accumulated wave component from the surviving kinetic particle component.

What would settle it

A direct numerical comparison in which the wave-particle solution deviates systematically from both the full kinetic reference in the rarefied limit and the Navier-Stokes reference in the continuum limit, for any choice of local horizon, would falsify the unification claim.

Figures

Figures reproduced from arXiv: 2606.26915 by Chang Liu, Kun Xu.

Figure 1
Figure 1. Figure 1: Positioning of the present wave-particle decomposition relative to the unified gas-kinetic (UGKS, [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: illustrates the structure of this local wave operator. The upper integration limit is the local kinetic horizon, the kernel contains the local relaxation target, and the exponential factor measures the collisionless attenuation along the backward molecular path. Local Time Horizon Local Maxwellian Source Collisional Decay Rate [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Macro-micro coupling procedure of the WPD algorithm. The macroscopic conservative state de [PITH_FULL_IMAGE:figures/full_fig_p023_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: L∞ convergence curves for the smooth-wave accuracy test. The computational domain is x ∈ [0, 2], with the initial interface at x = 1. A uniform mesh with Nx = 800 cells is used. The velocity space is discretized by Nξ = 513 uniform points over [−8, 8], and the final time is t = 0.2. The same physical mesh, velocity mesh, gas model, and time step are used for UGKS and WPD-SN ; in all WPD-SN results reported… view at source ↗
Figure 5
Figure 5. Figure 5: Sod shock tube at Kn = 1. X Density 0.5 1 1.5 0.2 0.4 0.6 (a) Density X Velocity 0.5 1 1.5 0 0.2 0.4 0.6 0.8 Ref. UGKS WP_SN (b) Velocity X Temperature 0.5 1 1.5 1.4 1.6 1.8 2 2.2 2.4 Ref. UGKS WP_SN (c) Temperature [PITH_FULL_IMAGE:figures/full_fig_p029_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Sod shock tube at Kn = 10−2 . 5.1.3. Couette flow The second one-dimensional benchmark is the planar Couette flow between two parallel diffuse-reflection walls. The lower wall is stationary and the upper wall moves tangentially, while the initial density and temperature are uniform. This test is sensitive to wall-induced 29 [PITH_FULL_IMAGE:figures/full_fig_p029_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Sod shock tube at Kn = 10−4 . rarefaction and velocity slip, and therefore provides a direct check of the kinetic wall treat￾ment in the deterministic WPD-SN formulation. The physical domain is discretized by Ny = 128 cells in the wall-normal direction. A sequence of Knudsen numbers from Kn = 1 down to 10−4 , including intermediate rarefaction levels, is considered. For rarefied and transition regimes, the… view at source ↗
Figure 8
Figure 8. Figure 8: Couette-flow profiles obtained by WPD-SN and the reference solutions. 5.1.4. Normal shock structure The third one-dimensional benchmark is the steady normal shock structure. This problem provides a more stringent kinetic test than the Sod tube because the solution contains a stationary non-equilibrium layer whose thickness and high-order moments are controlled by the collision model. We compute monatomic n… view at source ↗
Figure 9
Figure 9. Figure 9: Normal shock structure at M = 8 and Kn = 1. X Normalized -5 0 5 0 0.2 0.4 0.6 0.8 1 Ref. Density WP_SN Density Ref. Velocity WP_SN Velocity Ref. Temperature WP_SN Temperature Ma = 10, Kn=1 (a) M = 10, density X Stress Heat Flux -5 0 5 0 2 4 6 8 10 12 14 -100 -80 -60 -40 -20 Ref. Stress 0 WP_SN Stress Ref. Heat flux WP_SN Heat flux Ma = 10, Kn=1 (b) M = 10, stress [PITH_FULL_IMAGE:figures/full_fig_p032_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Normal shock structure at M = 10 and Kn = 1. pressible cavity, where the WPD formulation recovers the standard viscous vortex structure and agrees with the benchmark data of Ghia et al. For the WPD splitting, the local particle fraction is measured by βP = exp(−T /τ ), βW = 1 − βP , where T is the wave-particle horizon and τ is the local relaxation time. With CFLl = 1, the 32 [PITH_FULL_IMAGE:figures/ful… view at source ↗
Figure 11
Figure 11. Figure 11: Velocity distribution functions across the normal shock layer. [PITH_FULL_IMAGE:figures/full_fig_p033_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Rarefied lid-driven cavity flow at Kn = 1. X Y 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 RHO 1.18 1.16 1.14 1.12 1.1 1.08 1.06 1.04 1.02 1 0.98 0.96 0.94 0.92 0.9 0.88 0.86 (a) Density X Y 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 T 1.016 1.014 1.012 1.01 1.008 1.006 1.004 1.002 1 0.998 0.996 0.994 (b) Temperature [PITH_FULL_IMAGE:figures/full_fig_p035_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Rarefied lid-driven cavity flow at Kn = 0.075. 5.2.2. Continuum cavity The continuum limit is assessed by the classical lid-driven cavity at Re = 1000 and Re = 3200. In these cases the wave component dominates and the WPD scheme reduces to the continuum GKS behavior. The Re = 1000 case is computed on a 200 × 200 mesh, while the Re = 3200 case uses a finer 600 × 600 mesh. The velocity profiles are compared… view at source ↗
Figure 14
Figure 14. Figure 14: Centerline velocity profiles for rarefied cavity flows. [PITH_FULL_IMAGE:figures/full_fig_p036_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Continuum cavity flow at Re = 1000 compared with the Ghia et al. benchmark. 36 [PITH_FULL_IMAGE:figures/full_fig_p036_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Continuum cavity flow at Re = 3200 compared with the Ghia et al. benchmark. 5.3. Hypersonic Flow Around a Cylinder The hypersonic cylinder case examines the WPD framework in a non-equilibrium external flow with a curved diffuse-reflection wall. The free-stream Mach number is M∞ = 8. The cylinder radius is Rc = 1, and the far-field radius is 15Rc. A body-fitted O-grid is used with 100 cells in the circumfe… view at source ↗
Figure 17
Figure 17. Figure 17: Cylinder pressure and temperature contours at [PITH_FULL_IMAGE:figures/full_fig_p039_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Cylinder velocity contours at Kn = 1. 39 [PITH_FULL_IMAGE:figures/full_fig_p039_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: Cylinder stagnation-line density and wall pressure at [PITH_FULL_IMAGE:figures/full_fig_p040_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: Cylinder pressure and temperature contours at [PITH_FULL_IMAGE:figures/full_fig_p040_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: Cylinder velocity contours at Kn = 1. R RHO 1.5 2 2.5 3 3.5 4 5 10 15 20 WP_MC Reference (a) WPD-MC (local), ρ θ P 1 2 3 4 5 6 5 10 15 20 25 30 35 40 WP_SN Reference (b) WPD-MC (local), p [PITH_FULL_IMAGE:figures/full_fig_p041_21.png] view at source ↗
Figure 22
Figure 22. Figure 22: Cylinder stagnation-line density and pressure at [PITH_FULL_IMAGE:figures/full_fig_p041_22.png] view at source ↗
Figure 23
Figure 23. Figure 23: Equivalent particle-number fields at Kn = 1. more localized, and the equivalent particle-number field in [PITH_FULL_IMAGE:figures/full_fig_p042_23.png] view at source ↗
Figure 24
Figure 24. Figure 24: Cylinder pressure and temperature contours at [PITH_FULL_IMAGE:figures/full_fig_p042_24.png] view at source ↗
Figure 25
Figure 25. Figure 25: Cylinder velocity contours at Kn = 10−2 . R RHO 1.5 2 2.5 3 5 10 15 20 WP_SN Reference (a) WPD-SN , ρ θ P 1 2 3 4 5 6 5 10 15 20 25 30 35 40 45 WP_SN Reference (b) WPD-SN , p [PITH_FULL_IMAGE:figures/full_fig_p043_25.png] view at source ↗
Figure 26
Figure 26. Figure 26: Cylinder stagnation-line density and wall pressure at [PITH_FULL_IMAGE:figures/full_fig_p043_26.png] view at source ↗
Figure 27
Figure 27. Figure 27: Cylinder pressure and temperature contours at [PITH_FULL_IMAGE:figures/full_fig_p044_27.png] view at source ↗
Figure 28
Figure 28. Figure 28: Cylinder velocity contours at Kn = 10−2 . 44 [PITH_FULL_IMAGE:figures/full_fig_p044_28.png] view at source ↗
Figure 29
Figure 29. Figure 29: Cylinder stagnation-line density and pressure at [PITH_FULL_IMAGE:figures/full_fig_p045_29.png] view at source ↗
Figure 30
Figure 30. Figure 30: Equivalent particle-number fields at Kn = 10−2 . 5.3.3. Cylinder flow at Kn = 10−3 The Kn = 10−3 case further approaches the continuum limit. Figures 31–36 show that the wave component dominates most of the domain, while the particle component is retained where strong gradients and non-equilibrium effects remain. The pressure, temperature, and velocity contours are smooth and consistent between WPD-SN and… view at source ↗
Figure 31
Figure 31. Figure 31: Cylinder pressure and temperature contours at [PITH_FULL_IMAGE:figures/full_fig_p046_31.png] view at source ↗
Figure 32
Figure 32. Figure 32: Cylinder velocity contours at Kn = 10−3 . 46 [PITH_FULL_IMAGE:figures/full_fig_p046_32.png] view at source ↗
Figure 33
Figure 33. Figure 33: Cylinder stagnation-line density and wall pressure at [PITH_FULL_IMAGE:figures/full_fig_p047_33.png] view at source ↗
Figure 34
Figure 34. Figure 34: Cylinder pressure and temperature contours at [PITH_FULL_IMAGE:figures/full_fig_p047_34.png] view at source ↗
Figure 35
Figure 35. Figure 35: Cylinder velocity contours at Kn = 10−3 . R RHO 1.5 2 2.5 3 5 10 15 20 WP_MC Reference (a) WPD-MC (local), ρ θ P 1 2 3 4 5 6 5 10 15 20 25 30 35 40 45 WP_MC Reference (b) WPD-MC (local), p [PITH_FULL_IMAGE:figures/full_fig_p048_35.png] view at source ↗
Figure 36
Figure 36. Figure 36: Cylinder stagnation-line density and pressure at [PITH_FULL_IMAGE:figures/full_fig_p048_36.png] view at source ↗
Figure 37
Figure 37. Figure 37: Equivalent particle-number fields at Kn = 10−3 . The computational saving in the cylinder calculation is measured by both the transported kinetic representation and the wall time. On the 100 × 64 mesh, the deterministic SN solver uses 64 × 64 velocity points and two reduced distributions: NSN = 100 × 64 × 64 × 64 × 2 = 5.24 × 107 . For stochastic cases, the effective kinetic DoF is computed from the parti… view at source ↗
Figure 38
Figure 38. Figure 38: Three-dimensional x38 case at Kn = 10−1 . 50 [PITH_FULL_IMAGE:figures/full_fig_p050_38.png] view at source ↗
Figure 39
Figure 39. Figure 39: Three-dimensional x38 case at Kn = 10−1 . X Y Z p: 0.8 1.2 2 4 6 8 10 12 14 16 18 20 22 24 26 (a) Pressure X Y Z Mach: 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 (b) Mach number [PITH_FULL_IMAGE:figures/full_fig_p051_39.png] view at source ↗
Figure 40
Figure 40. Figure 40: Three-dimensional x38 case at Kn = 10−1 . 6. Conclusion In this work, we have formulated a wave-particle decomposition featuring a local kinetic horizon for kinetic relaxation equations. At the theoretical level, the characteristic integral solution effectively partitions the distribution function into an accumulated wave contribution and a collisionless particle contribution. The resulting unified wave-p… view at source ↗
Figure 41
Figure 41. Figure 41: Three-dimensional x38 case at Kn = 10−2 . X Y Z u: 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 (a) u X Y Z v: -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 (b) v [PITH_FULL_IMAGE:figures/full_fig_p052_41.png] view at source ↗
Figure 42
Figure 42. Figure 42: Three-dimensional x38 case at Kn = 10−2 . weighted Euler and Navier–Stokes fluxes, while the particle equation accurately captures the remaining kinetic transport. At the algorithmic level, we have constructed a conservative macro-micro numerical method based on this coupled system. The total conservative variables are advanced within a finite-volume framework using the sum of the wave and particle fluxes… view at source ↗
Figure 43
Figure 43. Figure 43: Three-dimensional x38 case at Kn = 10−2 . X Y Z rho: 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 (a) Density X Y Z T: 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 (b) Temperature [PITH_FULL_IMAGE:figures/full_fig_p053_43.png] view at source ↗
Figure 44
Figure 44. Figure 44: Three-dimensional x38 case at Kn = 10−4 . Coupling these components through the target state, the local-horizon source, and conser￾vative moment reconstruction ensures that the total macroscopic conservation law is strictly preserved at the discrete level. This formulation distinguishes itself from the global time-step wave-particle splitting of the UGKWP method by defining the decomposition directly at t… view at source ↗
Figure 45
Figure 45. Figure 45: Three-dimensional x38 case at Kn = 10−4 . X Y Z p: 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 (a) Pressure X Y Z Mach: 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 (b) Mach number [PITH_FULL_IMAGE:figures/full_fig_p054_45.png] view at source ↗
Figure 46
Figure 46. Figure 46: Three-dimensional x38 case at Kn = 10−4 . particle component is generated by the collisionless factor of the integral solution and is therefore interpreted as a fractional kinetic population rather than a signed residual. Our analysis demonstrates the discrete asymptotic-preserving continuum limit, rarefied￾regime consistency, and the appropriate scaling of active kinetic degrees of freedom. Com￾prehensiv… view at source ↗
read the original abstract

This paper presents a wave-particle decomposition (WPD) for kinetic relaxation equations, formulated around a local evolution timescale and its associated kinetic horizon. By leveraging the characteristic integral solution, we decompose the distribution function into an analytically accumulated wave component and a purely kinetic particle component. The latter is defined by the collisionless transport that survives beyond a prescribed local domain of influence, termed the horizon. This continuous formulation yields a unified wave-particle system valid across the entire Knudsen spectrum, comprising a source-free total conservation law, a wave equation, and a particle equation. The wave operator admits a Chapman--Enskog expansion, whose moments yield Euler and Navier--Stokes fluxes with horizon-dependent coefficients, while the particle equation governs the remaining non-equilibrium kinetic transport. At the algorithmic level, this system is discretized by a conservative macro-micro method. The total conservative variables are advanced by a finite-volume update using the sum of a Navier--Stokes gas-kinetic wave flux and a particle flux computed by either a deterministic discrete-ordinate $(S_N)$ method or a Monte Carlo representation. Unlike the global time-step splitting in the unified gas-kinetic wave-particle (UGKWP) method, the present partition is defined at the PDE level and governed by the local ratio of evolution timescale to relaxation time. The particle component is therefore a fractional kinetic population generated by the collisionless factor of the integral solution. Formal analysis establishes the asymptotic-preserving continuum limit, rarefied-regime consistency, and regime-adaptive scaling of active kinetic degrees of freedom. Numerical tests in one, two, and three dimensions validate the accuracy, multiscale capability, and efficiency of the framework.

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 wave-particle decomposition (WPD) for kinetic relaxation equations, formulated around a local evolution timescale and associated kinetic horizon. Leveraging the characteristic integral solution, the distribution is partitioned into an analytically accumulated wave component and a collisionless particle component surviving beyond the horizon. This yields a unified continuous system valid across the Knudsen spectrum, consisting of a source-free total conservation law, a wave equation whose moments admit a Chapman-Enskog expansion to horizon-dependent Euler and Navier-Stokes fluxes, and a particle equation for non-equilibrium transport. A conservative macro-micro discretization is introduced, with formal analysis establishing asymptotic-preserving continuum limits, rarefied-regime consistency, and regime-adaptive scaling of kinetic degrees of freedom; numerical tests in 1D-3D are presented.

Significance. If the decomposition is shown to be exact, this PDE-level approach would represent a meaningful extension of prior UGKWP frameworks by eliminating global time-step splitting in favor of a local, horizon-based partition. The claimed formal analysis of the AP property and the conservative macro-micro scheme would be strengths, as would the regime-adaptive particle population that reduces computational cost in near-continuum regions while retaining kinetic fidelity in rarefied zones. Successful validation across dimensions supports potential utility for multiscale problems.

major comments (2)
  1. [Formulation of the wave-particle system (around the local horizon)] The central claim that the decomposition produces a source-free total conservation law and a unified system valid for arbitrary Knudsen numbers rests on the premise that a local kinetic horizon can be defined such that the characteristic integral exactly separates the wave and particle components without residual sources or regime-dependent errors. An explicit formula for the horizon (in terms of local evolution timescale and relaxation time) and a proof that the separation holds for general initial data must be supplied, as this is load-bearing for the asymptotic-preserving limit and the absence of spurious sources.
  2. [Chapman--Enskog expansion of the wave operator] § on Chapman-Enskog expansion: the horizon-dependent coefficients in the resulting Euler and Navier-Stokes fluxes are asserted to arise from moments of the wave operator, but the explicit dependence on the horizon parameter and the reduction to standard fluid fluxes in the appropriate limits (e.g., horizon o au or horizon o au/ au) should be derived in detail to confirm consistency with the claimed continuum and rarefied limits.
minor comments (2)
  1. [Notation and definitions] The notation for the local evolution timescale and the collisionless factor in the integral solution could be clarified with a dedicated definition box or explicit equation reference to aid readers in following the regime-adaptive scaling argument.
  2. [Discretization and numerical tests] In the numerical section, the description of how the particle flux is computed via S_N or Monte Carlo should include a brief statement on conservation enforcement at the discrete level to match the continuous source-free property.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. We address the two major comments below and will revise the manuscript to supply the requested explicit formula, proof, and detailed derivations.

read point-by-point responses

  1. Referee: [Formulation of the wave-particle system (around the local horizon)] The central claim that the decomposition produces a source-free total conservation law and a unified system valid for arbitrary Knudsen numbers rests on the premise that a local kinetic horizon can be defined such that the characteristic integral exactly separates the wave and particle components without residual sources or regime-dependent errors. An explicit formula for the horizon (in terms of local evolution timescale and relaxation time) and a proof that the separation holds for general initial data must be supplied, as this is load-bearing for the asymptotic-preserving limit and the absence of spurious sources.

    Authors: We agree that an explicit formula and a self-contained proof are necessary to make the claims fully rigorous. The current manuscript defines the horizon via the local ratio of evolution timescale to relaxation time and derives the decomposition from the characteristic integral, but we will add the precise formula (horizon = f(τ_e, τ)) together with a proof that the integral solution separates exactly into wave and particle components for arbitrary initial data, with no residual sources. These additions will be placed in the formulation section to directly support the source-free conservation law and the AP property. revision: yes

  2. Referee: [Chapman--Enskog expansion of the wave operator] § on Chapman-Enskog expansion: the horizon-dependent coefficients in the resulting Euler and Navier-Stokes fluxes are asserted to arise from moments of the wave operator, but the explicit dependence on the horizon parameter and the reduction to standard fluid fluxes in the appropriate limits (e.g., horizon o au or horizon o au/au) should be derived in detail to confirm consistency with the claimed continuum and rarefied limits.

    Authors: We concur that the explicit dependence and the limiting cases require a more detailed derivation. In the revised manuscript we will expand the Chapman–Enskog section to display the horizon-dependent coefficients explicitly, derive the resulting Euler and Navier–Stokes fluxes from the moments of the wave operator, and prove the recovery of the standard fluid fluxes when the horizon tends to the relaxation time (continuum limit) and when the horizon tends to zero (rarefied limit). revision: yes

Circularity Check

0 steps flagged

No circularity: PDE-level decomposition is independently formulated

full rationale

The paper defines the wave-particle decomposition directly from the characteristic integral solution of the relaxation equation around a local evolution timescale and kinetic horizon, yielding a source-free conservation law plus separate wave and particle equations. This construction is presented as distinct from the prior UGKWP global time-step splitting, with the new partition governed at the PDE level by the local timescale ratio. Formal analysis of asymptotic-preserving limits and regime-adaptive scaling follows from the stated decomposition and Chapman-Enskog expansion of the wave operator; no step reduces a claimed result to a fitted input, self-citation chain, or definitional renaming. The horizon is an explicit modeling choice in the formulation rather than an output derived from the target properties.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claim rests on the definition of the local horizon and the integral solution decomposition, which introduces a new entity and relies on standard kinetic theory assumptions. No explicit fitted parameters are named in the abstract.

axioms (2)
  • standard math Characteristic integral solution exists for the kinetic equation
    Invoked to decompose the distribution function into wave and particle components.
  • domain assumption Chapman-Enskog expansion applies to the wave operator
    Used to obtain Euler and Navier-Stokes fluxes with horizon-dependent coefficients.
invented entities (1)
  • kinetic horizon no independent evidence
    purpose: To define the boundary between wave and particle components based on local domain of influence
    New concept introduced to separate the analytically accumulated wave from surviving kinetic transport.

pith-pipeline@v0.9.1-grok · 5825 in / 1645 out tokens · 56486 ms · 2026-06-26T02:30:26.225910+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

25 extracted references · 6 canonical work pages

  1. [1]

    Boltzmann, Weitere Studien uber das Warmegleichgewicht unter Gasmolekulen, Sitzungsberichte der Kaiserlichen Akademie der Wissenschaften66 (1872) 275–370

    L. Boltzmann, Weitere Studien uber das Warmegleichgewicht unter Gasmolekulen, Sitzungsberichte der Kaiserlichen Akademie der Wissenschaften66 (1872) 275–370

  2. [2]

    Cercignani,The Boltzmann Equation and Its Applications, Springer, New York, 1988

    C. Cercignani,The Boltzmann Equation and Its Applications, Springer, New York, 1988

    1988

  3. [3]

    Chapman, T.G

    S. Chapman, T.G. Cowling,The Mathematical Theory of Non-Uniform Gases, 3rd ed., Cambridge University Press, Cambridge, 1970

    1970

  4. [4]

    Shakhov, Generalization of the Krook kinetic relaxation equation,Fluid Dyn.3 (1968) 95–96

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

    1968

  5. [5]

    Bird,Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Oxford University Press, Oxford, 1994

    G.A. Bird,Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Oxford University Press, Oxford, 1994

    1994

  6. [6]

    Grad, On the kinetic theory of rarefied gases,Commun

    H. Grad, On the kinetic theory of rarefied gases,Commun. Pure Appl. Math.2 (1949) 331–407

    1949

  7. [7]

    Pareschi, G

    L. Pareschi, G. Russo, Implicit-explicit Runge-Kutta schemes and applications to hy- perbolic systems with relaxation,J. Sci. Comput.25 (2005) 129–155

    2005

  8. [8]

    Filbet, S

    F. Filbet, S. Jin, A class of asymptotic-preserving schemes for kinetic equations and related problems with stiff sources,J. Comput. Phys.229 (2010) 7625–7648

    2010

  9. [9]

    Liu, S.-H

    T.-P. Liu, S.-H. Yu, Boltzmann equation: Micro-macro decompositions and positivity of shock profiles,Commun. Math. Phys.246 (2004) 133–179

    2004

  10. [10]

    Gamba, S

    I.M. Gamba, S. Jin, L. Liu, Micro-macro decomposition based asymptotic-preserving numerical schemes and numerical moments conservation for collisional nonlinear kinetic equations,J. Comput. Phys.382 (2019) 264–290

    2019

  11. [11]

    Xu, J.-C

    K. Xu, J.-C. Huang, A unified gas-kinetic scheme for continuum and rarefied flows,J. Comput. Phys.229 (2010) 7747–7764

    2010

  12. [12]

    Xu,Direct Modeling for Computational Fluid Dynamics: Construction and Applica- tion of Unified Gas-Kinetic Schemes, World Scientific, 2014

    K. Xu,Direct Modeling for Computational Fluid Dynamics: Construction and Applica- tion of Unified Gas-Kinetic Schemes, World Scientific, 2014

    2014

  13. [13]

    Huang, K

    J.-C. Huang, K. Xu, P. Yu, A unified gas-kinetic scheme for continuum and rarefied flows II: multi-dimensional cases,Commun. Comput. Phys.12 (2012) 662–690. 56

    2012

  14. [14]

    C. Liu, K. Xu, A unified gas kinetic scheme for continuum and rarefied flows V: Multi- scale and multi-component plasma transport,Commun. Comput. Phys.22 (2017) 1175– 1223, doi:10.4208/cicp.OA-2017-0102

    work page doi:10.4208/cicp.oa-2017-0102 2017

  15. [15]

    C. Liu, Z. Wang, K. Xu, A unified gas-kinetic scheme for continuum and rarefied flows VI: Dilute disperse gas-particle multiphase system,J. Comput. Phys.386 (2019) 264– 295, doi:10.1016/j.jcp.2018.12.040

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

  16. [16]

    Xu, A gas-kinetic BGK scheme for the Navier-Stokes equations and its connection with artificial dissipation and Godunov method,J

    K. Xu, A gas-kinetic BGK scheme for the Navier-Stokes equations and its connection with artificial dissipation and Godunov method,J. Comput. Phys.171 (2001) 289–335

    2001

  17. [17]

    Z. Guo, K. Xu, R. Wang, Discrete unified gas kinetic scheme for all Knudsen number flows: low-speed isothermal case,Phys. Rev. E88 (2013) 033305

    2013

  18. [18]

    C. Liu, Y. Zhu, K. Xu, Unified gas-kinetic wave-particle methods I: Continuum and rarefied gas flow,J. Comput. Phys.401 (2020) 108977

    2020

  19. [19]

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

    work page doi:10.1186/s42774-021-00062-1 2021

  20. [20]

    Fluids31 (2019) 067105

    Y.Zhu, C.Liu, C.Zhong, K.Xu, Unifiedgas-kineticwave-particlemethodsII:Multiscale simulation on unstructured mesh,Phys. Fluids31 (2019) 067105

    2019

  21. [21]

    W. Li, C. Liu, Y. Zhu, J. Zhang, K. Xu, Unified gas-kinetic wave-particle methods III: Multiscale photon transport,J. Comput. Phys.408 (2020) 109280

    2020

  22. [22]

    X. Yang, K. Xu, Wave-particle based multiscale modeling and simula- tion of non-equilibrium turbulent flows,Computers & Fluids(2026) 107162, doi:10.1016/j.compfluid.2026.107162

    work page doi:10.1016/j.compfluid.2026.107162 2026

  23. [23]

    Z. Guo, K. Xu, Y. Zhu, A unified gas-kinetic framework from Boltzmann to Navier- Stokes scales,Advances in Aerodynamics(2026) 8:8, https://doi.org/10.1186/s42774- 025-00252-1

    work page doi:10.1186/s42774- 2026

  24. [24]

    L. Luo, L. Wu, Multiscale simulation of rarefied gas dynamics via di- rect intermittent GSIS-DSMC coupling,Advances in Aerodynamics(2024) 6:22, https://doi.org/10.1186/s42774-024-00188-y

    work page doi:10.1186/s42774-024-00188-y 2024

  25. [25]

    W. Su, L. Zhu, P. Wang, Y. Zhang, L. Wu, Can we find steady-state solutions to multiscale rarefied gas flows within dozens of iterations?J. Comput. Phys.407 (2020) 109245. 57

    2020