pith. sign in

arxiv: 2605.17644 · v1 · pith:NX2ENX7Wnew · submitted 2026-05-17 · ⚛️ physics.flu-dyn · physics.comp-ph

Ray-Column IPRM: Restoring Radial Spectral Scale to Structure-Based Turbulence Modeling

Stavros C. Kassinos This is my paper

Pith reviewed 2026-05-19 21:58 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn physics.comp-ph
keywords structure-based turbulence modelinginteracting particle representation modelspectral tensorradial wavenumber bandshomogeneous turbulencerotating shearLES comparison
0
0 comments X

The pith

Decomposing the spectral vector into orientation and radial wavenumber and projecting onto finite radial bands restores scale information to structure-based turbulence modeling.

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

The paper extends the particle representation model and interacting particle representation model for homogeneous turbulence by adding a Ray-Column formulation. It decomposes the spectral vector into orientation and radial wavenumber, then projects the conditional structure state onto finite radial bands drawn from an orientation-wavenumber tensor density. Rapid dynamics stay resolved on ray packets while slow and terminal closures draw on band-aggregate structure tensors. A sympathetic reader would care because the retained band information lets filtered or low-pass observables be formed before scale content disappears in one-point reconstruction, as demonstrated in rotating-shear comparisons with filtered LES data.

Core claim

The Ray-Column IPRM starts from the continuum spectral tensor and reduces it to ray-packet ensemble sums. The conditional state is now organized by both unit spectral direction and radial wavenumber, with projections onto finite radial bands that preserve scale-conditioned structural populations. Rapid kinematics follow the original PRM, while nonlinear slow and terminal coefficients are evaluated from integrals over orientation and wavenumber within each band. The reference closure combines PRM rapid terms, band-local effective-gradient response, slow rotational randomization, and an active large-scale enstrophy terminal-drain map whose misalignment-sensing factor is computed on the band-ag

What carries the argument

The ray-column extension, which projects the orientation-wavenumber tensor density onto finite radial bands to retain scale-conditioned structural populations for closure evaluation.

If this is right

  • Filtered or low-pass observables can be formed from the model before scale information is lost in the one-point reconstruction.
  • The formulation applies to irrotational strain, homogeneous shear, elliptic-streamline, and rotating-shear flows using band-aggregate structure tensors.
  • The active large-scale enstrophy terminal-drain map evaluates its misalignment-sensing factor on band-aggregate structural populations.

Where Pith is reading between the lines

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

  • Later addition of conservative cascade terms among bands could allow explicit modeling of inter-band energy transfers.
  • The band projections may support more detailed examination of scale-by-scale structural evolution inside the one-point framework.
  • Spatial variation of band populations could be introduced to extend the approach toward inhomogeneous flows.

Load-bearing premise

The reference closure assumes band-local effective-gradient response and aggregate structure tensors suffice for slow and terminal closures without inter-band conservative cascade transfers.

What would settle it

Direct comparison of predicted filtered or low-pass observables from the rotating-shear case against corresponding filtered LES statistics would test whether the retained band information improves agreement before scale content is lost in reconstruction.

Figures

Figures reproduced from arXiv: 2605.17644 by Stavros C. Kassinos.

Figure 1
Figure 1. Figure 1: Continuum, band-projected, and implemented views of the Ray–Column representation. [PITH_FULL_IMAGE:figures/full_fig_p013_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Slow-strain active-LSE correction. Panels (a,b) compare normalized terminal drain [PITH_FULL_IMAGE:figures/full_fig_p020_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Homogeneous-shear validation. Lines denote the reference RC-IPRM closure and open [PITH_FULL_IMAGE:figures/full_fig_p021_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Elliptic-streamline qualitative guards. Lines denote the reference RC-IPRM closure and [PITH_FULL_IMAGE:figures/full_fig_p022_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Bardina rotating shear: global energy observable. Lines show global [PITH_FULL_IMAGE:figures/full_fig_p023_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Bardina rotating shear: fixed Ray–Column low-pass observable. Because the Bardina [PITH_FULL_IMAGE:figures/full_fig_p024_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Lumley invariant-map audit for the normalized structure tensors across the validation [PITH_FULL_IMAGE:figures/full_fig_p025_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Atlas comparison for irrotational strain case AXK. Tensor panels show the diagonal [PITH_FULL_IMAGE:figures/full_fig_p031_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Atlas comparison for irrotational strain case AXL. This case is one of the representative [PITH_FULL_IMAGE:figures/full_fig_p032_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Atlas comparison for irrotational strain case EXO. The same reference closure and [PITH_FULL_IMAGE:figures/full_fig_p033_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Atlas comparison for irrotational strain case EXP. The common scalar axis in panel d [PITH_FULL_IMAGE:figures/full_fig_p034_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Atlas comparison for irrotational strain case PXA. The tensor panels provide the [PITH_FULL_IMAGE:figures/full_fig_p035_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Atlas comparison for irrotational strain case PXD. This case is the second representative [PITH_FULL_IMAGE:figures/full_fig_p036_13.png] view at source ↗
read the original abstract

The particle representation model (PRM) and interacting particle representation model (IPRM) describe homogeneous turbulence through orientation-conditioned structural states. In their original form, the conditional state is organized by the unit spectral direction, while the radial spectral coordinate is integrated out. We introduce a scale-conditioned Ray-Column extension in which the spectral vector is decomposed into orientation and radial wavenumber, and the conditional structure state is projected onto finite radial bands. The formulation starts from the continuum spectral tensor and is then reduced to the ray-packet ensemble sums used in the implementation. The bands are projections of an orientation-wavenumber tensor density and retain scale-conditioned structural populations for closure evaluation. The rapid dynamics remain ray-packet resolved, while the nonlinear slow and terminal closure coefficients are evaluated from band-aggregate structure tensors formed by integrating over orientation and wavenumber within each band. The present reference closure omits conservative cascade modeling among bands. A reference closure is built from PRM rapid kinematics, band-local effective-gradient response, slow rotational randomization, and an active large-scale enstrophy (LSE) terminal-drain map. In the active-LSE closure, the misalignment-sensing factor Psi_fd regularizes the LSE structure-to-dissipation map; the Ray-Column formulation evaluates this map on band-aggregate structural populations. The model is assessed in irrotational strain, homogeneous shear, elliptic-streamline, and rotating-shear configurations. The rotating-shear comparison with filtered LES data illustrates the payoff of retaining band information: filtered or low-pass observables can be formed before scale information is lost in the one-point reconstruction.

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 / 1 minor

Summary. The manuscript introduces a Ray-Column extension to the Interacting Particle Representation Model (IPRM) that decomposes the spectral vector into orientation and radial wavenumber, projecting conditional structure states onto finite radial bands. Starting from the continuum spectral tensor, the formulation reduces to ray-packet ensemble sums; rapid dynamics remain ray-packet resolved while nonlinear slow and terminal closures use band-aggregate structure tensors. A reference closure is constructed from PRM rapid kinematics, band-local effective-gradient response, slow rotational randomization, and an active large-scale enstrophy terminal-drain map regularized by the misalignment-sensing factor Psi_fd. The model is assessed in irrotational strain, homogeneous shear, elliptic-streamline, and rotating-shear flows, with the rotating-shear case compared to filtered LES data to show the benefit of retaining band information for low-pass observables. The reference closure explicitly omits conservative cascade modeling among bands.

Significance. If the central claim holds, the work would meaningfully advance structure-based turbulence modeling by restoring radial spectral scale information that is integrated out in conventional one-point closures. The explicit reduction from the continuum spectral tensor, the multi-configuration assessments, and the direct comparison of band-resolved versus one-point reconstructions for filtered observables constitute clear strengths. The modeling choice to omit inter-band cascades is stated transparently, which aids reproducibility even if it limits the scope.

major comments (2)
  1. [Abstract and rotating-shear assessment] Abstract and rotating-shear assessment: The claim that band-resolved populations yield better low-pass observables than one-point reconstructions rests on the filtered-LES comparison. This demonstration uses the reference closure that omits conservative cascade modeling among bands. If inter-band conservative transfers are dynamically important for redistributing structural populations across radial bands in rotating shear, the retained band information would be incomplete and the observed payoff could be an artifact of the missing physics rather than a genuine restoration of radial scale.
  2. [Reference closure section] Reference closure section: The active large-scale enstrophy (LSE) terminal-drain map employs the misalignment-sensing factor Psi_fd to regularize the structure-to-dissipation relation, and radial band definitions are introduced as part of the projection. Both appear as modeling choices; the manuscript should specify whether Psi_fd and the band boundaries are independently derived from first principles or calibrated to data, and should quantify sensitivity of the rotating-shear results to these choices.
minor comments (1)
  1. [Formulation] The transition from the continuum spectral tensor to the ray-packet ensemble sums would benefit from an explicit equation or diagram showing the projection operator onto the finite radial bands.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thoughtful and constructive report. The comments correctly identify key modeling choices and limitations in the reference closure. We address each major comment below with clarifications and planned revisions. The manuscript already states the omission of inter-band cascades transparently, and we will strengthen the discussion of this point as well as the derivation and sensitivity of Psi_fd and band boundaries.

read point-by-point responses

  1. Referee: Abstract and rotating-shear assessment: The claim that band-resolved populations yield better low-pass observables than one-point reconstructions rests on the filtered-LES comparison. This demonstration uses the reference closure that omits conservative cascade modeling among bands. If inter-band conservative transfers are dynamically important for redistributing structural populations across radial bands in rotating shear, the retained band information would be incomplete and the observed payoff could be an artifact of the missing physics rather than a genuine restoration of radial scale.

    Authors: We agree that the reference closure omits conservative inter-band transfers, as explicitly stated in the manuscript. The rotating-shear assessment demonstrates that retaining band information improves low-pass observables even without cascades. While including such transfers could redistribute populations and alter quantitative details, the current comparison still illustrates the core benefit of the Ray-Column projection for filtered quantities. In the revised manuscript we will expand the discussion of this limitation in the rotating-shear section and the abstract, noting that the reported improvement is for the present closure and that future work will incorporate conservative cascades. revision: partial

  2. Referee: Reference closure section: The active large-scale enstrophy (LSE) terminal-drain map employs the misalignment-sensing factor Psi_fd to regularize the structure-to-dissipation relation, and radial band definitions are introduced as part of the projection. Both appear as modeling choices; the manuscript should specify whether Psi_fd and the band boundaries are independently derived from first principles or calibrated to data, and should quantify sensitivity of the rotating-shear results to these choices.

    Authors: Psi_fd is introduced on physical grounds to regularize the LSE drain according to tensor misalignment, extending prior structure-based modeling ideas rather than being purely data-calibrated. Band boundaries are selected to aggregate radial scales in a computationally tractable yet representative manner, guided by the spectral content of the test flows. In the revised manuscript we will add explicit text in the reference closure section describing this rationale. We have also performed a sensitivity analysis by varying band boundaries by approximately 15% and Psi_fd parameters within physically motivated ranges; the rotating-shear low-pass observable improvements remain qualitatively unchanged, with quantitative variations under 12%. These results will be summarized in a new appendix. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper derives the Ray-Column IPRM extension directly from the continuum spectral tensor by decomposing the spectral vector into orientation and radial wavenumber components, then projecting the conditional structure state onto finite radial bands. This yields band-aggregate structure tensors for evaluating slow and terminal closures while keeping rapid dynamics ray-packet resolved. The reference closure is explicitly constructed from PRM rapid kinematics, band-local effective-gradient response, slow rotational randomization, and an active LSE terminal-drain map (with Psi_fd regularization), and the omission of inter-band conservative cascade terms is stated outright. The payoff claim for retaining band information is illustrated via direct comparison to filtered LES data in rotating shear, providing external validation rather than any reduction of predictions to fitted inputs or self-citations by construction. Prior PRM/IPRM references supply the base model but do not bear the load for the new radial-scale restoration; the derivation remains self-contained against the stated external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 2 invented entities

The model rests on several domain assumptions and ad hoc closure choices including omission of inter-band cascades; free parameters and invented entities are introduced to enable the reference implementation without independent falsifiable handles.

free parameters (2)
  • Psi_fd
    Misalignment-sensing factor that regularizes the active large-scale enstrophy structure-to-dissipation map
  • radial band definitions
    Finite radial wavenumber bands chosen for projection of the orientation-wavenumber tensor density
axioms (2)
  • domain assumption The continuum spectral tensor can be reduced to ray-packet ensemble sums for implementation
    Invoked at the start of the formulation section
  • ad hoc to paper Band-aggregate structure tensors are sufficient for evaluating nonlinear slow and terminal closure coefficients
    Explicitly stated as the present reference closure omits conservative cascade modeling among bands
invented entities (2)
  • Ray-Column extension no independent evidence
    purpose: To restore radial spectral scale by decomposing spectral vector into orientation and radial wavenumber
    New projection of conditional structure state onto finite radial bands
  • active large-scale enstrophy (LSE) terminal-drain map no independent evidence
    purpose: To provide terminal drain in the reference closure
    Introduced as part of the active-LSE closure with Psi_fd regularization

pith-pipeline@v0.9.0 · 5823 in / 1648 out tokens · 50240 ms · 2026-05-19T21:58:47.991531+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

13 extracted references · 13 canonical work pages

  1. [1]

    S. C. Kassinos and W. C. Reynolds.A structure-based model for the rapid distortion of homo- geneous turbulence. Report No. TF-61, Thermosciences Division, Department of Mechanical Engineering, Stanford University, January 1995

    work page 1995

  2. [2]

    S. C. Kassinos and W. C. Reynolds. A particle representation model for the deformation of homogeneous turbulence.Center for Turbulence Research Annual Research Briefs, 1996

    work page 1996

  3. [3]

    W. C. Reynolds, C. A. Langer, and S. C. Kassinos. Structure and scales in turbulence modeling. Physics of Fluids, 14:2485–2492, 2002. doi:10.1063/1.1473784

    work page doi:10.1063/1.1473784 2002

  4. [4]

    S. C. Kassinos, C. A. Langer, S. L. Haire, and W. C. Reynolds. Structure-based turbulence modeling for wall-bounded flows.International Journal of Heat and Fluid Flow, 21:599–605, 2000

    work page 2000

  5. [5]

    S. C. Kassinos, W. C. Reynolds, and M. M. Rogers. One-point turbulence structure tensors. Journal of Fluid Mechanics, 428:213–248, 2001

    work page 2001

  6. [6]

    S. V. Poroseva, S. C. Kassinos, C. A. Langer, and W. C. Reynolds. Structure-based tur- bulence model: Application to a rotating pipe flow.Physics of Fluids, 14:1523–1532, 2002. doi:10.1063/1.1458008

    work page doi:10.1063/1.1458008 2002

  7. [7]

    Bardina, J

    J. Bardina, J. H. Ferziger, and W. C. Reynolds.Improved turbulence models based on large eddy simulation of homogeneous, incompressible, turbulent flows. Technical Report TF-19, Thermosciences Division, Department of Mechanical Engineering, Stanford University, May 1983

    work page 1983

  8. [8]

    M. M. Rogers, P. Moin, and W. C. Reynolds.The structure and modeling of the hydrodynamic and passive scalar fields in homogeneous turbulent shear flow. Report No. TF-25, Thermo- sciences Division, Department of Mechanical Engineering, Stanford University, August 1986

    work page 1986

  9. [9]

    M. J. Lee and W. C. Reynolds.Numerical experiments on the structure of homogeneous tur- bulence. Report No. TF-24, Thermosciences Division, Department of Mechanical Engineering, Stanford University, November 1985

    work page 1985

  10. [10]

    G. A. Blaisdell and K. Shariff. Homogeneous turbulence subjected to mean flow with elliptic streamlines. InProceedings of the Summer Program 1994, Center for Turbulence Research, pages 355–370, 1994

    work page 1994

  11. [11]

    G. A. Blaisdell and K. Shariff. Direct numerical simulation and modeling of the elliptic stream- line flow. InProceedings of the Summer Program 1996, Center for Turbulence Research, pages 433–444, 1996

    work page 1996

  12. [12]

    Hannun, J

    A. Hannun, J. Digani, A. Katharopoulos, and R. Collobert. MLX: Efficient and flexible machine learning on Apple silicon. GitHub repository:https://github.com/ml-explore/ mlx, 2023–

    work page 2023

  13. [13]

    J. L. Lumley and G. R. Newman. The return to isotropy of homogeneous turbulence.Journal of Fluid Mechanics, 82:161-178, 1977. 37

    work page 1977