Pith. sign in

REVIEW 3 major objections 4 minor 22 references

Curvilinear coordinates and curvature in radiative transport

T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A single curvature formula gives the streaming operator in any curvilinear transport frame

desk verdict The central streaming formula (Eq 27) is not correct—the derivation drops a nonzero off-diagonal shape-operator term and mis-signs the κ_n term—but the geometric framework and the classical examples are solid enough that the paper deserves a salvageable referee round. read the letter →

arxiv 2508.20852 v1 pith:RFM7WZ6Q submitted 2025-08-28 math.NA cs.NAmath.APnucl-th

Johannes Krotz , Ryan G. McClarren This is my paper · ORCID
classification math.NAcs.NAmath.APnucl-th
keywords radiativetransportstreamingoperatorcurvilinearcoordinatesorthonormalframecurvatureshapeintegralcurvesfoliation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper derives a general expression for the streaming term Ω·∇Ψ in radiative-transport and kinetic equations when the direction Ω is parametrized by a projection µ onto a spatially varying unit vector n and an angle ω in the orthogonal plane spanned by t and b. The expression rewrites the angular derivatives in geometric terms: the drift in µ is controlled by the curvature of a surface normal to n and by the curvature of the integral curves of n, while the drift in ω is controlled by how t and b bend and twist. If correct, this gives practitioners one master formula for writing the streaming operator in any curvilinear coordinate system, replacing the case-by-case geometric derivations common in the literature. It also makes term-by-term vanishing visible: flat surfaces, straight integral curves, geodesics, and non-twisting frames all eliminate terms without further computation.

What carries the argument

The machinery is the parametrization Ω(μ,ω)=μn+√(1−μ²)(cos ω t+sin ω b), with (n,t,b) an integrable orthonormal frame allowed to vary in space. Treating μ and ω as spatial fields converts Ω·∇Ψ into a chain rule, and the coefficients are identified with the shape operator of the foliation surface X^n and with curvature vectors κ_n, κ_t, κ_b of the frame's integral curves. Orthonormal-frame identities reduce the μ-coefficient to a squared cosine/sine combination, and the homothetic-scaling remark shows why all curvature terms scale as 1/ρ for self-similar families of surfaces.

What would settle it

For the translating-paraboloid frame with a≠b, compute t·∇_b n + b·∇_t n directly at a generic point; the sum is nonzero, so repeating the derivation from Eq. (17) without discarding the sin(ω)cos(ω) term changes the coefficient of ∂Ψ/∂μ. Direct numerical differentiation of Ω·∇μ along the direction Ω̂_∥ would then show the missing cross term.

Watch

Extended reading notes

Core claim

The paper's central claim is Eq. (27), a general formula for the streaming operator when Ω is written in a local orthonormal frame (n,t,b): Ω·∇Ψ equals the spatial derivative Ω·∇_r Ψ, plus a term multiplying ∂Ψ/∂μ built from the normal curvature C(r,ω) of the surface orthogonal to n and the curvature κ_n of the integral curves of n, plus a term multiplying ∂Ψ/∂ω built from the curvatures κ_t, κ_b and the twisting t·∇_n b of the frame around γ_n. Two interchangeable forms are given, one in terms of curve curvatures and one in terms of shape operators, and the paper demonstrates the formula on cylindrical, spherical, elliptical, and translated-graph coordinates, including a homothetic-scaling

Load-bearing premise

The load-bearing step at Eq. (17) assumes t·∇_b n + b·∇_t n = 0, equivalent to t and b being principal directions of the surface X^n; the paper does not state this, the translating-paraboloid example with a≠b violates it, and if it fails Eq. (27) misses the sin(ω)cos(ω) cross term.

Editorial extensions

If this is right

  • In any frame where the surfaces X^n are flat and the integral curves of n are straight, the coefficient of ∂Ψ/∂μ vanishes, simplifying or removing the μ-dependence of the streaming operator.
  • For frames whose t and b curves are geodesics of X^n and do not twist around γ_n, the entire ∂Ψ/∂ω term drops; this is the geometric reason the sphere example reduces to a single sin(ω)cot(θ)/ρ term.
  • The formula turns a tedious algebraic chore into a geometric check: each streaming coefficient can be diagnosed as zero from flatness, straightness, geodesy, or non-twisting before any curvature is computed.
  • For homothetic foliations such as spheres, cylinders, and ellipsoids, all curvature terms scale as 1/ρ, so the streaming coefficients can be computed once at one scale and rescaled.
  • When X^n are level sets of Ψ with holonomy SO(2), the appendix argument shows ∂Ψ/∂ω vanishes, giving the standard reduced transport equation in μ and space only.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the frame's t and b directions are not principal directions of X^n, the sin(ω)cos(ω) cross term kept at Eq. (17) is nonzero; Eq. (27) would then need an extra correction term, so the practical domain of the formula is frames adapted to the surface curvature.
  • The curve-curvature and shape-operator forms of the angular coefficients can be seen as two gauges for the same geometric drift; switching between them may help build the conservative discretizations the appendix shows are generically unavailable outside planes and spheres.
  • The translated-graph examples suggest a practical recipe for layered or extruded geometries: since every streaming coefficient depends only on the base graph G0 and not on the translation height z, coefficients can be precomputed once from the cross-section profile.
  • The appendix's holonomy condition reframes dimensional reduction in ω as a parallel-transport property of the orthogonal foliation, which could guide coordinate selection for kinetic models with toroidal or helical symmetry.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The manuscript derives a general expression for the streaming operator Ω·∇Ψ in transport equations when the direction Ω is parameterized by (μ,ω) in a local orthonormal frame (n,t,b). It proposes Eq. (27) as a universal formula involving the surface curvature C(r,ω), the curvature κ_n of the integral curve of n, the curvatures κ_t,κ_b, and a winding term t·∇_n b. The derivation is split into Ω·∇μ and Ω·∇ω, with curve-curvature and surface-curvature forms, and is illustrated on cylindrical, spherical, ellipsoidal, and translating-graph coordinate examples.

Significance. The geometric framework is appealing and, if correct, would be a useful reference for practitioners choosing coordinate systems in kinetic/radiative transfer. The paper also makes a pedagogical contribution by expressing angular derivative terms through curvature of integral curves and shape operators, and the examples are worked in detail. However, the central formula is invalid as stated: the derivation of ∇Ωμ contains a missing off-diagonal shape-operator term and a sign error. These errors change the streaming operator in all but trivial geometries, so the paper's main claim of a general streaming expression is not supported.

major comments (3)
  1. [§2.1, Eq. (17)] The derivation cancels the term sin(ω)cos(ω)(t·∇_b n + b·∇_t n) by invoking Remark 1. This is incorrect: Remark 1 gives t·∇_b n = −n·∇_b t and b·∇_t n = −n·∇_t b, which are not negatives of each other. For a surface with normal n and orthonormal tangents t,b, Weingarten symmetry yields t·∇_b n = b·∇_t n = S_tb, so the cross term is 2(1−μ²)sinω cosω S_tb. It vanishes only if t,b are principal directions. Integrability of the frame does not imply this. The paper's own translating-paraboloid frame has nonzero S_tb and S_bt (Eqs. (56)–(57)), so the cancellation is invalid in general and in the paper's own example. Consequently Eqs. (18), (21), and (27) are missing this term.
  2. [§2.1, Eqs. (15)–(16), (21)] The sign convention in Remark 2 defines κ_u := −∇_u u. In Eq. (16) the term μ(Ω∥·∇_n n) is replaced by μ√(1−μ²)(Ω̂∥·κ_n). This has the wrong sign: ∇_n n = −κ_n, so Eq. (16) should be −μ√(1−μ²)(Ω̂∥·κ_n). The same sign error propagates into Eqs. (21) and (27). A direct check on the translating paraboloid confirms that the sign must be flipped for the formula to reproduce Ω·∇(Ω·n).
  3. [§3.5, Eqs. (56)–(57)] The translating-paraboloid example is advertised as an illustration of the general framework, but it explicitly violates the assumption needed for the Eq. (17) cancellation: the off-diagonal shape-operator entries S_tb and S_bt are nonzero (Eqs. (56)–(57)). At (x,y)=(1,1) with a=1,b=2, the sum t·∇_b n + b·∇_t n is nonzero, so the sinω cosω term does not vanish. Thus the example cannot be used to validate Eq. (27); it instead demonstrates the missing term.
minor comments (4)
  1. [Keywords] Typos: 'Curvelinear coordiantes' should be 'Curvilinear coordinates'.
  2. [§3.4, text after Eq. (48)] The sentence 'All that is left is to plug these terms back into the formulas for ∇ΩΨ, which we will refrain from due to the length of the terms' is awkwardly phrased; the authors could either provide the final expression or explicitly say it is omitted for brevity.
  3. [§2.2, Eqs. (25)–(26)] The surface-curvature form of ∇Ωω is introduced without a fully detailed derivation, and the notation (1,0)S_Xb(...) is hard to parse. Since this form is not used in the final examples, the authors should either clarify it or remove it.
  4. [General] There are several typographical inconsistencies, e.g., 'Froebenius' for 'Frobenius', and some equation numbers are referenced imprecisely (e.g., 'equation (28)' in Section 3 where Eq. (27) is meant).

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the streaming-operator derivation is a self-contained differential-geometry calculation.

full rationale

The paper's central derivation, Eqs. (10)-(27), proceeds by applying the chain rule, orthonormal-frame identities (Remark 1), the definition of integral-curve curvature (Remark 2), and standard results on foliations and the Weingarten map (Remarks 3-4), each cited to external classical sources (Frobenius, Lee, Kreyszig, Weingarten). No fitted parameter is renamed as a prediction, no target formula is used as an input, and no load-bearing assumption is justified solely by a self-citation. The quantity C(r,ω) introduced in Eq. (20) is not an independent input; it is a shorthand for the already-computed quadratic form (cosω, sinω)·S_X^n·(cosω, sinω)^T. All examples are substitutions into the derived general expression. Even if a reader disputes the algebraic correctness of a step (e.g., the cancellation at Eq. (17) or the sign convention for κ_n), that is a mathematical-correctness concern, not circularity. The derivation is self-contained against external benchmarks, so the circularity score is 0.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No fitted parameters or invented physical entities. The derivation rests on standard differential geometry plus an unstated, false-in-general algebraic assumption that the off-diagonal shape operator component vanishes.

assumptions (4)
  • domain assumption n, t, b form an integrable orthonormal vector field
    Assumed in §1 and used throughout; integrability allows surfaces X^n orthogonal to n to exist via Remark 3.
  • standard math Frobenius theorem characterization of foliations in R3
    Used in Remark 3 to assert V·rot V = 0 is necessary and sufficient for surfaces orthogonal to V.
  • standard math Weingarten theorem and shape operator identities
    Used in Remarks 4 and §2.1 to identify derivative terms with surface curvature.
  • ad hoc to paper The cross term t·∇_b n + b·∇_t n vanishes
    Implicitly assumed in Eq (17) when simplifying (1−μ²)sin(ω)cos(ω)(t·∇_b n + b·∇_t n) to zero. Not stated, not true in general; corresponds to the off-diagonal shape operator term being zero.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Curvilinear coordinates and curvature in radiative transport." pith.science (2026). https://pith.science/paper/RFM7WZ6Q

@misc{pith2026250820852,
  author       = {Pith},
  title        = {Pith review of: Curvilinear coordinates and curvature in radiative transport},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RFM7WZ6Q}},
  note         = {Machine review of arXiv:2508.20852}
}
read the original abstract

We derive a general expression for the streaming term in radiative transport equa- tions and other transport problems when formulated in curvilinear coordinates, emphasizing coordinate systems adapted to the geometry of the domain and the directional dependence of particle transport. By parametrizing the angular vari- able using a local orthonormal frame, we express directional derivatives in terms of curvature-related quantities that reflect the geometry of underlying spatial man- ifolds. Our formulation highlights how the interaction between coordinate choices and curvature influences the streaming operator, offering geometric interpretations of its components. The resulting framework offers intuitive insight into when and how angular dependence can be simplified and may guide the selection of coordinate systems that balance analytical tractability and computational efficiency.

Figures

Figures reproduced from arXiv: 2508.20852 by the authors.

Figure 1
Figure 1. Given a parametrization, i.e., a set of coordinates for r ∈ R 3 , say r = r(u 1 , u2 , u3 ), a natural choice for n, t, and b would be the canonical tangent vectors n = c1 ∂r ∂u1 , t = c2 ∂r ∂u2 , and b = n × t, where the coefficients c1, c2 ensure normalization. Note that this canonical choice is not necessarily orthogonal; if it is not, additional modifications are required. Orthogonalization via Gram-Schmidt is o… view at source ↗
Figure 1
Figure 1. Visualization of the parametrization of [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. 2D visualization of γ n in spherical coordinates and elliptical coordinates. The integral curves γ n following the normal vectors of circles are straight lines, while those of ellipses are curved. Setting u = n in the previous remark, we find that the second term in equation (13) is µ(Ω∥ · ∇nn) = µ p 1 − µ2(Ωˆ ∥ · ∇nn) = µ p 1 − µ2  Ωˆ ∥ · κ n  (16) expressing the term entirely in terms of µ and the curvature of t… view at source ↗
Figures from the paper (7 more)
Figure 3
Figure 3. Figure 3: Example manifold Xn with normal vector n and tangent vectors t and Ω∥ . The curvature at r in direction ω relative to t, denoted by C(r, ω) is the normal curvature of red curve starting at r. By Weingarten’s theorem, the coordinate representation of SX can also be calc…
Figure 4
Figure 4. Figure 4: : Example of non-geodesic integral curve [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]
Figure 5
Figure 5. Figure 5: : Example of a twisting, but flat manifold (helicoid). Diagonal entries of [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]
Figure 6
Figure 6. Figure 6: Visualization of Ω-basis choice n, t, b in cylindrical symmetry. 3.1. Cylindrical I Let us begin with a classical example useful in cases with cylindrical symmetry. As￾sume that r is expressed as r = (ρ cos(ϕ), ρ sin(ϕ), z) through the coordinates ρ, ϕ and z. Common ch…
Figure 7
Figure 7. Figure 7: Visualization of Ω-basis choice n, t, b in spherical symmetry. curves circling the poles parallel to the equator. The choice of n, t and b is depicted in figure 7 for clarity. We express equation (5) here is as Ω · ∇Ψ = Ω · ∇rΨ + h (1 − µ 2 )C(r, ω) + µ p 1 − µ2  Ωˆ ∥…
Figure 8
Figure 8. Figure 8: : 2D example of homothetically rescaled manifolds. Base manifold [PITH_FULL_IMAGE:figures/full_fig_p018_8.png]
Figure 9
Figure 9. Figure 9: : 2D example of translating graphs. Base graph [PITH_FULL_IMAGE:figures/full_fig_p020_9.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

22 extracted references · 19 canonical work pages

  1. [1]

    Simon R. Arridge. Optical tomography in medical imaging. Inverse Problems, 15 0 (2): 0 R41--R93, 1999

    work page 1999

  2. [3]

    On the modeling of social crowds: A survey of models, speculations, and perspectives

    Nicola Bellomo, Abdelghani Bellouquid, Yulong Tao, and Michael Winkler. On the modeling of social crowds: A survey of models, speculations, and perspectives. Mathematical Models and Methods in Applied Sciences, 23 0 (04): 0 1230004, 2013

    work page 2013

  3. [4]

    Sur les groupes d'holonomie homog\`ene des vari\'et\'es \`a connexion affine et des vari\'et\'es riemanniennes

    Marcel Berger. Sur les groupes d'holonomie homog\`ene des vari\'et\'es \`a connexion affine et des vari\'et\'es riemanniennes. Bull. Soc. Math. France, 83: 0 279--330, 1955. ISSN 0037-9484. URL http://www.numdam.org/item?id=BSMF_1955__83__279_0

    work page 1955

  4. [5]

    The Boltzmann Equation and Its Applications, volume 67 of Applied Mathematical Sciences

    Carlo Cercignani. The Boltzmann Equation and Its Applications, volume 67 of Applied Mathematical Sciences. Springer, 1988

    work page 1988

  5. [6]

    Radiative Transfer

    Subrahmanyan Chandrasekhar. Radiative Transfer. Dover Publications, 1960

    work page 1960

  6. [7]

    A survey of closed self-shrinkers with symmetry

    Gregory Drugan, Hojoo Lee, and Xuan Hien Nguyen. A survey of closed self-shrinkers with symmetry. Results in Mathematics, 73 0 (1): 0 32, 2018. ISSN 1420-9012. doi:10.1007/s00025-018-0763-3. URL https://doi.org/10.1007/s00025-018-0763-3

    work page doi:10.1007/s00025-018-0763-3 2018

  7. [8]

    Freimanis

    J. Freimanis. On vector radiative transfer equation in curvilinear coordinate systems. Journal of Quantitative Spectroscopy and Radiative Transfer, 112 0 (13): 0 2134--2148, 2011. ISSN 0022-4073. doi:https://doi.org/10.1016/j.jqsrt.2011.04.007. URL https://www.sciencedirect.com/science/article/pii/S0022407311001646. Polarimetric Detection, Characterizatio...

    work page doi:10.1016/j.jqsrt.2011.04.007 2011

  8. [9]

    Frobenius

    G. Frobenius. Ueber das pfaffsche problem. Journal für die reine und angewandte Mathematik, 1877 0 (82): 0 230--315, 1877. doi:doi:10.1515/crll.1877.82.230. URL https://doi.org/10.1515/crll.1877.82.230

    work page doi:10.1515/crll.1877.82.230

Show all 22 references
  1. [10]

    Self-similar shrinkers in R ^3 , January 2008

    Tobias Kohn. Self-similar shrinkers in R ^3 , January 2008

    2008

  2. [11]

    Kreyszig

    E. Kreyszig. Differential Geometry. Dover Books on Mathematics. Dover Publications, 2013. ISBN 9780486318622. URL https://books.google.com/books?id=Vw3CAgAAQBAJ

    2013

  3. [12]

    Edward W. Larsen. The description of particle transport problems with helical symmetry. Nuclear Science and Engineering, 156 0 (1): 0 68--73, 2007. doi:10.13182/NSE07-A2685

    2007 doi

  4. [13]

    John M. Lee. Introduction to Smooth Manifolds. 2000

    2000

  5. [14]

    John M. Lee. Introduction to Riemannian Manifolds. Springer, Cham, Switzerland, second edition, 2018. doi:10.1007/978-3-319-91755-9. URL https://doi.org/10.1007/978-3-319-91755-9

    2018 doi

  6. [15]

    Microscale heat conduction in dielectric thin films

    Arun Majumdar. Microscale heat conduction in dielectric thin films. Journal of Heat Transfer, 115 0 (1): 0 7--16, 1993

    1993

  7. [16]

    Nicholson

    Dwight R. Nicholson. Introduction to Plasma Theory. John Wiley & Sons, 1983

    1983

  8. [17]

    Cosmic Ray Astrophysics

    Reinhard Schlickeiser. Cosmic Ray Astrophysics. Springer, 2002

    2002

  9. [18]

    Weingarten

    J. Weingarten. Ueber eine klasse auf einander abwickelbarer flächen. Journal für die reine und angewandte Mathematik, 59: 0 382--393, 1861. URL http://eudml.org/doc/147844

  10. [19]

    @esa (Ref

    \@ifxundefined[1] #1\@undefined \@firstoftwo \@secondoftwo \@ifnum[1] #1 \@firstoftwo \@secondoftwo \@ifx[1] #1 \@firstoftwo \@secondoftwo [2] @ #1 \@temptokena #2 #1 @ \@temptokena \@ifclassloaded agu2001 natbib The agu2001 class already includes natbib coding, so you should ...

  11. [20]

    \@lbibitem[] @bibitem@first@sw\@secondoftwo \@lbibitem[#1]#2 \@extra@b@citeb \@ifundefined br@#2\@extra@b@citeb \@namedef br@#2 \@nameuse br@#2\@extra@b@citeb \@ifundefined b@#2\@extra@b@citeb @num @parse #2 @tmp #1 NAT@b@open@#2 NAT@b@shut@#2 \@ifnum @merge>\@ne @bibitem@firs...

  12. [21]

    !1A Qa

    @open @close @open @close and [1] URL: #1 \@ifundefined chapter * \@mkboth \@ifxundefined @sectionbib * \@mkboth * \@mkboth\@gobbletwo \@ifclassloaded amsart * \@ifclassloaded amsbook * \@ifxundefined @heading @heading NAT@ctr thebibliography [1] @ \@biblabel @NAT@ctr \@bibset...

    2011

  13. [22]

    , " * write output.state after.block = add.period

    ENTRY address archive author booktitle chapter collaboration doi edition editor eid howpublished institution journal key lastchecked month note number numpages organization pages publisher school series title translator type url urldate volume year label extra.label sort.label...

  14. [23]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 gl...

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.