{"id":"389823fe-50b1-4e42-9242-139e37cf5123","arxiv_id":"2502.07930","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Derives approximate spherical and cylindrical versions of the Hammer-Rosen Marshak-wave model, with front positions matching numerical diffusion solutions to about 10 percent for inner radii of 0.1 cm and larger.","lead":"The authors extend the Hammer-Rosen radiation diffusion model from planar slabs to spherical and cylindrical shells, producing a fast ODE for the heat-front position. The model matches numerical diffusion simulations to roughly 10 percent for inner radii as small as 0.1 cm, relevant to inertial confinement fusion capsules.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (32) mis-evaluates the logarithmic integral from Eq. (29), so the printed profile Eq. (36) is not an O(epsilon) solution of the first-order equation; the consistency claim fails even for r0 >= 0.1 cm.","rationale":"The reader's weakest assumption concerned the small-r0 validity of setting the left side of Eq. (21) to zero. That is a genuine limitation and is acknowledged by the authors. However, the more load-bearing defect is algebraic and occurs even in the claimed regime r0 >= 0.1 cm: the step from Eq. (29) to Eq. (32) replaces an integral of a logarithm by a function whose derivative is not that logarithm. As a result, Eq. (36), which is explicitly part of the claimed consistent solution, does not satisfy the first-order equation it is supposed to solve. The front-position ODE Eq. (37) can be recovered from Eq. (29) by the y=1 boundary condition and may still be numerically useful, which explains why the front comparisons look reasonable. But the central claim as worded includes the profile Eq. (36) as part of a consistent solution, and that claim is false as printed. Because the error is independent of resolution, code availability, or the small-r0 asymptotics, the current manuscript cannot be accepted even under the reader's stated conditions. A corrected derivation of the profile (or an explicit statement that Eq. (36) is a separate heuristic profile) would be needed before the consistency claim can be credited; hence the recommendation moves from CONDITIONAL to REJECT.","tokens_in":11429,"tokens_out":30313,"duration_ms":256455,"concrete_test":"Set epsilon=0, constant drive, d=1, gamma=1 cm, r0=0.1 cm. From Eq. (35) compute a=(1/gamma)*d rf/ds=(r0/gamma)*ln(1+gamma/r0)=0.240. Using the printed Eq. (36), form zeta and evaluate the residual R(y)=d*ln((gamma*y+r0)/(gamma+r0)) - d zeta/dy - a at y=0.1: R ~ -1.20. Recompute R using the exact antiderivative F(y)=(r0/gamma)*ln((r0+gamma*y)/r0)+y*ln((r0+gamma*y)/(r0+gamma))-y in place of the psi term of Eq. (33); then R vanishes to machine precision. If the residual is not O(epsilon), Eq. (36) is not the first-order Picard iterate and the profile portion of the central claim is refuted.","verdict_should_be":"REJECT","load_bearing_attack":"The central claim is that Eq. (37) and Eq. (36) together form a consistent approximate solution. The front ODE can be defensible: imposing zeta(1)=0 on the correctly integrated first-order equation reproduces Eq. (35). The profile, however, does not. In passing from Eq. (31) to Eq. (32), the integral of d*ln((gamma*yhat+r0)/(gamma+r0)) from 0 to y is replaced by d*[-y+(r0/gamma)*ln(1+gamma*y/r0)], but the derivative of that replacement is -d*gamma*y/(gamma*y+r0), not the logarithm. The exact antiderivative is F(y)=(r0/gamma)*ln((r0+gamma*y)/r0)+y*ln((r0+gamma*y)/(r0+gamma))-y. Thus Eq. (36) satisfies neither Eq. (29) nor the original Eq. (21) to O(epsilon), independently of how small r0 is. Quantitatively, for epsilon=0, d=1, gamma=1, r0=0.1, y=0.1, Eq. (29) gives ln(0.2/1.1) ~ -1.70, while d zeta/dy from Eq. (36) plus a=(r0/gamma)*ln(1+gamma/r0) ~ 0.240 is -0.50; the residual is ~ -1.20, an O(1) violation. The sign inconsistency noted by the reader in Eq. (24)/(28) is related: Eq. (28) does not follow from Eq. (26). The numerical agreement of the front position cannot validate the profile claim; a corrected profile must be derived and the profile comparisons redone.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript extends the planar, arbitrary-drive radiation diffusion model of Hammer and Rosen to diverging spherical and cylindrical geometries. The authors nondimensionalize the nonlinear diffusion equation with power-law material properties, introduce a moving heat-front coordinate, and use a perturbation/Picard iteration in ε = β/(4+α) to obtain an ODE for the front position rf(t) (Eq. 37) and a temperature profile (Eq. 36). They compare front position and profiles against finite-volume numerical solutions for four materials and inner radii from 1000 cm down to 0.001 cm, reporting roughly 10% agreement for r0 ≥ 0.1 cm and documenting breakdown at smaller radii.","tokens_in":11830,"tokens_out":16889,"duration_ms":132521,"significance":"If the derivation is corrected, the front-position ODE is a useful extension: it is free of fitted parameters, uses C and ε from tabulated material properties, and preserves the arbitrary-drive feature of the planar model. The numerical comparisons in Figures 3–4 provide external support for the ODE. However, the profile portion is currently not a consistent O(ε) solution because of an integration error in Eq. (32); the claimed 'consistent approximate solution' is therefore not yet established, and the temperature-profile comparisons in Figures 5–8 need to be redone. The contribution is potentially significant for semi-analytic modeling of ICF-scale curved systems, but the revision is nontrivial.","major_comments":[{"comment":"The step from Eq. (29) to Eq. (32) is not correct. The exact antiderivative of d ln((γy+r0)/(γ+r0)) on [0,y] is d[ y ln((γy+r0)/(γ+r0)) + (r0/γ) ln(1+γy/r0) − y ], but Eq. (32) retains only the last two terms. Consequently Eq. (36) does not satisfy the first-order equation Eq. (29) to O(ε). The missing term is not negligible: for r0 = 0.1 cm, γ = 1, y = 0.1 it contributes about −0.17 to ζ^{1−ε}, the same size as the retained geometric correction. Because the missing term vanishes at y = 1, the front-position ODE obtained by imposing ζ(1) = 0 may still be correct, but the interior profile is not. Please correct Eq. (32), derive the corresponding profile, and repeat the profile comparisons in Figures 5–8.","section":"II, Eqs. (29)–(36)"},{"comment":"There is a sign inconsistency in the zeroth-order solution. Equation (24) as printed, 1 = (ε−1)(1/γ)∂rf/∂s, gives a negative front speed for ε < 1, and it does not lead to Eq. (25). The subsequent use of (1/γ)∂rf/∂s = 1 + O(ε) in Eq. (31) and the positive front speeds reported in Section III require 1 = (1−ε)(1/γ)∂rf/∂s. Please correct Eq. (23)/(24) and the text around them.","section":"II, Eqs. (23)–(24) and (31)"},{"comment":"The treatment of the geometric term as a small correction is not quantified. The left-hand side of Eq. (21) contains dγ/(yγ+r0)∂ζ/∂y, which for r0 = 0.1 cm and γ of order 0.1–0.25 cm is O(1), not O(ε). The assertion that this side is small 'away from the wavefront' is especially delicate because the term is largest near y = 0. The manuscript should state a smallness criterion (e.g., a bound on γ/r0) and verify it for the cases where the 10% agreement is claimed; otherwise the 'ranges of applicability' are only empirically motivated, not supported by the asymptotic ordering.","section":"II, around Eq. (21)"}],"minor_comments":[{"comment":"The dummy variable of integration should not be denoted y; use a different symbol such as ŷ or y′.","section":"II, Eq. (28)"},{"comment":"The solution of Eq. (22) displayed in Eq. (23) should be shown to satisfy the boundary conditions; this will also clarify the sign correction needed in Eq. (24).","section":"II, Eq. (23)"},{"comment":"References 4 and 16 are the same paper (Milovich et al.); consolidate them.","section":"References"},{"comment":"Reference 1 contains a typo: 'Univeristy' should be 'University'.","section":"References"},{"comment":"The claim of 'about 10% or less' for r0 ≥ 0.1 cm would be better supported by a table of relative errors at 3 ns for each material and r0.","section":"III.A"},{"comment":"Figures 5–8 should state in the captions which analytic profile (equation number) is plotted, so that the corrected profile can be checked against them.","section":"III.B"}],"recommendation":"major_revision","confidential_remarks":"The manuscript contains a useful front-position ODE that appears to be supported by the numerical comparisons, but the profile derivation is currently invalid because of the mis-evaluated logarithmic integral in Eq. (32). The error is local and fixable, so major revision is appropriate rather than rejection. I would ask the editor to require a corrected profile and quantitative profile comparisons, not just a statement that the ODE is unchanged."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper extends Hammer-Rosen's planar Marshak-wave model to spherical and cylindrical shells, giving an ODE for the heat-front position and a closed-form temperature profile. The front ODE looks genuinely useful: it reduces to HR in the planar limit, and the front-position plots agree with numerical diffusion runs to about 10% for r0 ≥ 0.1 cm. That part is worth something.\n\nBut the paper as printed is not consistent. The stress-test note is right: Eq. (32) mishandles the integral of the logarithmic term in Eq. (29). The exact antiderivative has an extra y log((r0+γy)/(r0+γ)) term; the paper drops it. As a result, Eq. (36) is not an O(ε) solution of the first-order equation. I checked the derivative: for ε=0, d=1, γ=1, r0=0.1, y=0.1, Eq. (29) gives -1.70 while Eq. (36) gives -0.50. That's an O(1) residual, not a small correction. So the profile part of the \"consistent solutions\" claim fails, independent of how small r0 is.\n\nThere's also the sign slip in Eq. (24): ζ(1)=0 forces (ε-1)(1/γ)∂rf/∂s = 1, i.e., negative speed for ε<1. Later they use 1+O(ε). That's probably a typo, but it's sloppy in a key spot.\n\nThe good news: the front ODE, Eq. (35)/(37), survives. The erroneous log term vanishes at y=1, so the boundary condition still gives the same ODE. I rederived it with the correct antiderivative and got the same expression. That's why the front-position plots look okay. The temperature profiles, however, are plotted from an incorrect formula, so they cannot be trusted.\n\nWhat's genuinely new is the curvilinear correction and the resulting ODE, which is absent from HR and the other cited work. The paper is honest about the small-r0 limitation and doesn't fit any constants to numerics. But the derivation errors and the lack of code/data mean it needs major revision before publication.\n\nThis deserves a serious referee, but the referee should ask for a corrected profile derivation and replotted figures, plus a reproducibility artifact. I wouldn't cite it in its current form.","headline":"Useful curvilinear extension with a correct-looking front ODE, but the printed profile is mathematically wrong and the consistency claim fails as written.","tokens_in":12299,"tokens_out":7569,"would_cite":false,"duration_ms":53707,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that the planar quasi-analytic model for radiation-driven heat waves can be extended to spherical and cylindrical geometries, yielding a heat-front position that matches numerical diffusion solutions to about 10 percent…","keywords":["radiation diffusion","heat front","spherical geometry","cylindrical geometry","quasi-analytic model","inertial confinement fusion","fixed-point iteration"],"falsifier":"Run a high-resolution numerical diffusion solve for a small-$\\epsilon$ material at $r_0 = 0.05$ cm and compare the heat-front trajectory and the inner-boundary temperature gradient to Eqs. (36)-(37); if the front position lies within 10 percent of the PDE result, the claimed applicability range would extend below 0.1 cm, while a large slowdown and a sharp gradient at $r_0$ would confirm the stated failure mode. Equivalently, evaluate the dropped left-hand side of Eq. (21) on the numerical solution and check whether it is small compared with the right-hand side.","tokens_in":11243,"feed_emoji":"🔥","tokens_out":11612,"duration_ms":95539,"temperature":0.7,"pith_summary":"Radiation-driven heat waves in a spherical or cylindrical shell—heated from an inner surface and expanding outward—are usually studied with full numerical diffusion simulations. This paper claims that a quasi-analytic model, built on the planar radiation-diffusion solution, can be extended to these curved geometries while keeping its biggest practical advantage: the drive temperature at the inner boundary can be an arbitrary function of time. The extension reduces the problem to one ordinary differential equation for the heat-front position $r_f(t)$ and a formula for the temperature profile, and the paper reports that for inner radii $r_0 \\ge 0.1$ cm the predicted front position agrees with numerical diffusion solutions to about 10 percent. The agreement fails as $r_0$ shrinks, because curvature produces a steep temperature gradient near the inner boundary that the perturbative construction does not capture; the model then moves the front too slowly. If correct, the model makes fast iteration on small curved systems such as inertial-confinement-fusion capsules possible without running full simulations for every design variation.","feed_headline":"Heat waves in spheres and cylinders matched to ~10 percent","feed_subtitle":"A quasi-analytic curved-geometry heat-wave model matches numerical diffusion down to 0.1 cm inner radii.","key_machinery":"The load-bearing structure is the moving-coordinate transformation $y = (r-r_0)/(r_f(t)-r_0)$, which turns the expanding wave into a fixed-boundary problem, and the smallness of $\\epsilon = \\beta/(4+\\alpha)$, the ratio formed from the power-law exponents of internal energy and opacity. Setting the left-hand side of the transformed equation to zero away from the front gives the planar-like profile $\\zeta_0 = 1-y$; feeding that profile back through the full equation yields the logarithmic geometric correction $\\psi$ and a consistency condition that becomes the front-position ODE. This ODE is the one piece that must be integrated numerically; once $r_f(t)$ is known, the temperature profile and albedo follow from closed-form expressions. The same machinery reduces cleanly to the planar case as $r_0 \\to \\infty$, and it produces a geometric correction twice as large in spherical as in cylindrical geometry.","core_discovery":"At its core, the paper establishes that the consistency condition at the moving heat front—the requirement that the temperature and its gradient vanish together at the front—can be written as an ODE for the front position in spherical and cylindrical geometry, not just in planar geometry. The derivation transforms the nonlinear radiation-diffusion equation into a moving coordinate $y=(r-r_0)/(r_f-r_0)$, expands in the small material exponent $\\epsilon = \\beta/(4+\\alpha)$, and uses one round of fixed-point iteration around the planar profile $\\zeta_0 = 1-y$ to build in a logarithmic geometric correction. The final front-position ODE, Eq. (37), together with the profile, Eq. (36), is the claimed 'consistent solution'. The paper validates it against finite-volume numerical solutions for four materials spanning a range of $\\epsilon$ and for inner radii from 0.01 to 1000 cm, finding roughly 10 percent agreement in front position for $r_0 \\ge 0.1$ cm and systematic slowdown at smaller radii.","pith_inferences":["The same perturbative strategy—drop the geometric correction at lowest order, then iterate—might be applied to converging waves by reversing the sign of the correction, but the semi-infinite domain assumption would need to be replaced; the paper only identifies converging waves as future work.","Because the failure at small $r_0$ stems from the inner-boundary gradient, a boundary-layer treatment near $r_0$ matched to the outer profile could extend the model below 0.1 cm; this is a testable modification.","The accuracy threshold $r_0 \\ge 0.1$ cm is likely material-dependent through $\\epsilon$; a systematic mapping of the failure radius versus $\\epsilon$ and geometry would turn the reported examples into a design rule.","If adopted as a reduced-order model in multi-dimensional codes, the front-position ODE could serve as a fast subgrid model for radiation flow in curved geometry, with error estimates derived from the $\\epsilon$ and $1/r_0$ expansion."],"forward_implications":["For inner radii down to 0.1 cm, the model gives heat-front positions accurate to about 10 percent for the four materials tested, so design scans can replace expensive PDE solves in that regime.","Because the drive temperature remains arbitrary, the model can handle shaped or time-varying radiation drives, not just constant or power-law boundary conditions.","In cylindrical geometry the geometric correction is smaller than in spherical geometry, so the model is accurate at even smaller inner radii in cylinders.","The model reproduces the qualitative waiting behavior of nonlinear diffusion fronts when the drive changes from constant to ramped, although it exaggerates the front slowdown.","The temperature profile and material albedo are available semi-analytically once the front ODE is integrated, which is what enables rapid iteration for design studies."],"supporting_citations":[{"why":"Supplies the planar quasi-analytic model being extended, including the arbitrary drive temperature and the approximate solution procedure.","marker":"14"},{"why":"Shows that the planar model disagrees with curvilinear simulations at capsule scale, motivating the extension.","marker":"15"},{"why":"Provides the finite-volume spatial discretization used to produce the numerical reference solutions.","marker":"19"},{"why":"Provides the adaptive time-integration solver used for the numerical reference solutions.","marker":"20"},{"why":"Supplies the material parameter sets for the four materials used in the comparison.","marker":"21"}],"fun_headline_variants":["Curved geometry heat waves: analytic model nails 10%","Heat wave model extends to small spheres and cylinders","Spherical and cylindrical heat waves match numerics within 10%","New quasi-analytic heat wave model for tiny capsules","Radiation diffusion in curved shells: model validated to 10%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes that the left-hand side of the transformed diffusion equation—the drive term and the geometric correction—is small away from the heat front, so it can be dropped before the profile iteration; this becomes false as the inner radius approaches zero, where a steep temperature gradient appears near the boundary.","fun_headline_variants_meta":{"raw":{"variants":["Curved geometry heat waves: analytic model nails 10%","Heat wave model extends to small spheres and cylinders","Spherical and cylindrical heat waves match numerics within 10%","New quasi-analytic heat wave model for tiny capsules","Radiation diffusion in curved shells: model validated to 10%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000278,"raw_usage":{"total_tokens":1611,"prompt_tokens":860,"completion_tokens":751,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":476,"completion_tokens_details":{"reasoning_tokens":668}},"tokens_in":476,"tokens_out":751,"duration_ms":7732,"temperature":1.0,"reasoning_tokens":668,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T11:23:49.098795+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a high-resolution numerical diffusion solve for a small-$\\epsilon$ material at $r_0 = 0.05$ cm and compare the heat-front trajectory and the inner-boundary temperature gradient to Eqs. (36)-(37); if the front position lies within 10 percent of the PDE result, the claimed applicability range would extend below 0.1 cm, while a large slowdown and a sharp gradient at $r_0$ would confirm the stated failure mode. Equivalently, evaluate the dropped left-hand side of Eq. (21) on the numerical solution and check whether it is small compared with the right-hand side.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the planar quasi-analytic model being extended, including the arbitrary drive temperature and the approximate solution procedure."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows that the planar model disagrees with curvilinear simulations at capsule scale, motivating the extension."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the finite-volume spatial discretization used to produce the numerical reference solutions."},{"cited_title":"Krief \\ and\\ author R","cited_arxiv_id":null,"evidence_quote":"Provides the adaptive time-integration solver used for the numerical reference solutions."},{"cited_title":"Malka \\ and\\ author S","cited_arxiv_id":null,"evidence_quote":"Supplies the material parameter sets for the four materials used in the comparison."}],"review_version":1}