{"id":"ad5d9f6a-f8b7-444a-a81a-f600e605ea2a","arxiv_id":"2411.18117","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Weakly viscoelastic tube drawing equations show elastic effects typically accelerate hole closure and destabilize the process, except for very large inlet holes or weak stretching.","lead":"This paper works out the math of stretching a stretchy polymer tube into a fibre and finds that, surprisingly, the elasticity usually makes the central hole close faster rather than slower. The results give fibre makers a better picture of when polymer materials will preserve or close internal holes, and how stable the drawing process will be.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'elasticity always destabilizes at Re=0' conclusion rests on a first-order eigenvalue calculation whose equations are omitted, so the manuscript alone cannot support it.","rationale":"The reader's stated weakest assumption is the validity of the O(De) perturbation and the exclusion of hole-closure regimes. That concern is real but explicitly scoped by the authors, and it does not invalidate the hole-closure mechanism inside the stated non-closure regime. The stability claim is different in kind: Section 5 draws a universal conclusion ('always destabilizing' at Re=0) from an omitted calculation. The leading-order Newtonian stability problem is recoverable, but the first-order correction omega_1 requires the adjoint problem and solvability condition, and the manuscript does not give enough information to check their derivation. This is the most load-bearing weakness because the abstract and conclusions advertise the stability result as a central, counterintuitive finding. My read does not move the verdict: the paper should remain conditional on supplying the omitted equations or reproducible code. I did not find a concrete sign error or internal inconsistency in the steady-state O(De) system; the C^(1) decomposition and the physical explanation via the second normal stress difference are coherent, and the q restriction to phi >= 0.5 is stated explicitly. Thus the hole-closure part is credible, while the stability part is unverifiable as submitted.","tokens_in":26876,"tokens_out":11177,"duration_ms":110011,"concrete_test":"Obtain the omitted first-order eigenvalue system and adjoint solvability condition from the authors (for example, the Mathematica notebook referenced in Section 5), or require the full equations in an appendix, and independently recompute Re(omega_1) at D0c for the Re=0 cases in Figures 12-15, including at least phi=0.95, Ca=2, alpha=0, beta=0. If the recomputed sign is not positive at the critical draw ratio, the 'always destabilizing' claim is false; if it is positive, the manuscript should still state the complete equations so the result is reproducible.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 5 states that the linearized first-order equations and the Fredholm/solvability condition for omega_1 can be derived with Mathematica, but the detailed expressions are omitted. The central stability claim—Re(omega_1)>0 at the Newtonian critical draw ratio for Re=0, asserted as 'always', in contrast to solid viscoelastic threads—depends entirely on the sign of this omitted correction. The leading-order operator is standard, but the first-order correction involves a non-self-adjoint eigenvalue problem, an adjoint boundary-value problem, and delicate integration by parts; an error in any boundary term would reverse the conclusion. The paper supplies no code, no tables of eigenvalue data, and no explicit adjoint, and the plotted sweeps cover a finite parameter set. This is not an internal contradiction, but it is a load-bearing verification gap. The steady-state hole-closure mechanism is presented much more transparently and appears coherent; the stability result, however, is not independently checkable from the manuscript as it stands, and if it fails, the abstract's contrast with solid-thread drawing must be withdrawn.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives long-wave asymptotic equations for the drawing of an axisymmetric Giesekus tube at small Deborah number, starting from the full slender geometry and treating De as a perturbation about the Newtonian tube-drawing state of Fitt et al. (2002). The O(De) correction fields are computed from the leading-order solution, and the authors show that the outlet hole correction h(1)(z=1) is typically negative (enhanced hole closure) except for large inlet hole ratios phi or draw ratios close to unity, which they explain through the second normal stress difference and its axial/radial advection. The paper then presents a linear stability analysis of the same base state, reporting that for negligible inertia (Re=0) the elastic correction to the growth rate omega1 has positive real part at the Newtonian critical draw ratio, while for non-zero inertia it can be stabilizing or destabilizing depending on parameters.","tokens_in":27072,"tokens_out":2399,"duration_ms":25025,"significance":"If the results are correct, this is a valuable contribution to the fluid mechanics of microstructured fibre drawing: it gives the first systematic weakly-viscoelastic tube-drawing model, identifies the second normal stress difference as the controlling mechanism, and produces a falsifiable prediction about the sign of the hole-size correction. The steady-state part is transparent and self-contained, and the leading-order reduction correctly matches the known Newtonian equations; the paper also makes a clear contrast with solid-thread viscoelastic drawing, which is mechanistically interesting. The main limitation is that the central stability claim is not independently checkable from the manuscript because the first-order eigenvalue system and the solvability condition for omega1 are omitted, with no code, tables, or explicit adjoint supplied. The steady-state derivation and the physical explanation in Section 4 are strong and, in my assessment, the more reliable part of the paper.","major_comments":[{"comment":"The central claim that \"elastic effects are always destabilizing for negligible inertia\" (abstract and Section 6) rests entirely on the sign of Re(omega1) at the Newtonian critical draw ratio. However, Section 5 states that \"the detailed expressions of these equations are omitted here for the sake of brevity\" and that the eigenvalue problem can be derived using Mathematica. The paper supplies neither the first-order eigenvalue equations, the adjoint problem, the Fredholm/solvability algebra, nor any numerical data (tables, convergence studies, code) that would allow the reader to verify the sign of Re(omega1). Because the first-order correction involves a non-self-adjoint eigenvalue problem and integration by parts with boundary terms, a sign error in any boundary contribution would reverse the conclusion. This is a load-bearing verification gap: the abstract's contrast with solid-thread drawing cannot be supported by the manuscript as it stands. I request that the omitted expressions be written out, or the relevant code/data be made available, at least for the representative parameter sets plotted in Figures 12-15.","section":"§5 (linear stability analysis)"},{"comment":"The paper's abstract and conclusions state that elastic effects \"enhance hole closure for most parameter values,\" but the analysis explicitly excludes the hole-closure regime phi < 0.5 and small capillary numbers (text below Eq. (4.3) and Section 4.1). Within the stated domain the claim is supported by the plotted sweeps, but the phrase \"most parameter values\" is broader than the computed and physically regular regime. The authors should either soften the universality language or, better, state precisely the parameter region over which h(1)(z=1)<0 has been verified, including the De-sufficiency condition needed to prevent O(De^2) terms from changing the sign.","section":"§4.1, Figures 3-6"},{"comment":"The expansion omega = omega0 + De*omega1 assumes that the eigenvalue perturbation is analytic in De and that the leading-order operator has a simple eigenvalue at the critical point. This is plausible for the plotted cases but not demonstrated; if the eigenvalue is defective or if the critical draw ratio is not isolated, the first-order correction to the growth rate would not be given by the stated solvability condition. The manuscript should at least state the non-defectiveness assumption or verify it numerically for the reported eigenvalues.","section":"§5, Eqs. (5.1)-(5.6)"}],"minor_comments":[{"comment":"The scaling for Ca in Eq. (2.25) uses the total viscosity eta0 and the small parameter epsilon; this is fine, but the reader should be reminded that Ca depends on epsilon, so the statement that Ca is O(1) in the asymptotic limit is an ordering assumption; a one-line clarification would help.","section":"§2.1, around Eq. (2.14)"},{"comment":"The labels 1-4 in Eq. (4.7) are very helpful, but the text describes terms 1 and 2 as arising from axial and radial advection without giving the corresponding signs explicitly for term 2. Since the sign of the integral of C(1) is the crux of the hole-size mechanism, a short sign table or a direct inequality would make Section 4.2 easier to follow.","section":"§4.2, Eq. (4.7)"},{"comment":"There are several typographical artifacts, e.g. the LaTeX commands \"/greaterorsimilar\" in Sections 4.1.1 and 5 should be rendered as symbols, and the reference \"Geyling, F. and GM, H. (1980)\" appears to have an incomplete author name. A final proofreading pass is needed.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the steady-state analysis is a solid contribution and could be published after the stability section is made verifiable. The omission of the first-order eigenvalue system is not merely a presentation issue because the paper's headline claim about always-destabilizing elasticity depends on it. I would suggest asking the authors to either include the full equations (perhaps in an appendix) or provide reproducible code/data for the stability results; without that, the stability conclusions should be explicitly labelled as numerical observations over a finite parameter set rather than universal claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about this paper. First, it is the first systematic derivation of long-wave equations for a weakly viscoelastic tube with an internal hole, and the steady-state part is careful and credible: the De=0 limit reduces to Fitt et al.'s Newtonian tube equations, and the O(De) corrections are computed from those leading-order solutions without fitted parameters. Second, the headline stability result — elasticity always destabilizes at Re=0, unlike solid viscoelastic threads — rests on an eigenvalue calculation whose explicit equations are omitted. As it stands, that claim is not independently checkable from the manuscript.\n\nWhat is genuinely new is the hole-closure mechanism. The decomposition of C^(1) into four terms, and the explanation in terms of radial versus axial advection of the second normal stress difference, is illuminating and appears physically sound. The finding that elasticity enhances hole closure for moderate holes but suppresses it for very large holes or weak drawing is a useful, non-obvious result. The steady-state numerics look consistent, and the figures match the described behavior.\n\nThe soft spots are concentrated in Section 5. The text says the first-order eigenvalue problem, adjoint, and Fredholm solvability condition 'can easily be derived using Mathematica' but gives none of the expressions, no code, no tables of omega_1, and no convergence checks. The conclusion 'always destabilizing' at Re=0 is supported only by plotted curves over a finite parameter set. This is not an internal contradiction, but it is load-bearing: if any boundary term in the integration by parts has the wrong sign, the abstract's contrast with solid threads is wrong. A referee cannot verify it from the paper alone. The fix is straightforward in principle — include the omitted equations in an appendix or ship the Mathematica notebook/data — but the current manuscript does not have it.\n\nA minor, explicitly-stated caveat: the 'most parameter values' phrasing covers only the regime phi >= 0.5 and moderate-to-large Ca, where leading-order hole closure does not occur. The paper says this, so it is not a hidden flaw, but the abstract slightly overstates the scope.\n\nBottom line: this is a serious, careful asymptotic paper with a credible steady-state contribution and an interesting physical mechanism. I would send it to peer review, and would accept it after the authors either provide the omitted stability derivations/code or restrict the 'always' claim to the computed cases. I would cite the steady-state equations and the mechanism; I would be cautious about citing the stability conclusion.","headline":"Careful asymptotic derivation of viscoelastic tube drawing with a credible hole-closure mechanism; the 'always destabilizing' stability claim is not checkable from the manuscript without the omitted eigenvalue details or code.","tokens_in":27589,"tokens_out":3851,"would_cite":true,"duration_ms":31121,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76A10","76D45","76E15"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper establishes that weak polymer elasticity in tube drawing mostly hastens hole closure and, for negligible inertia, always destabilizes the process.","keywords":["viscoelastic tube drawing","Giesekus model","weakly viscoelastic asymptotics","second normal stress difference","hole closure","draw resonance","microstructured optical fibres","linear stability"],"falsifier":"Run a full two-dimensional numerical simulation of the Giesekus equations for $De \\approx 0.1$, $D \\approx 1.5$, $Ca \\approx 1.8$, $\\varphi \\approx 0.5$, $Re = 0$, $\\alpha = 0$, $\\beta = 0.1$: if the outlet hole radius is larger than the Newtonian value, the predicted sign $h^{(1)}(z=1) < 0$ is wrong. Separately, at the critical draw ratio for $Re=0$, compute the real part of the elastic growth-rate correction $\\omega_1$: a negative value would disprove the claim that elasticity is always destabilizing when inertia is negligible.","tokens_in":26669,"feed_emoji":"🕳️","tokens_out":10673,"duration_ms":88715,"temperature":0.7,"pith_summary":"The paper asks whether adding polymer elasticity to the fluid used to draw microstructured ('holey') optical fibres helps or hurts the two manufacturing goals: keeping the internal hole open and avoiding oscillatory draw resonance. Working in the weakly viscoelastic limit of the Giesekus constitutive model, it derives one-dimensional long-wave equations as small corrections to the well-studied Newtonian tube-drawing flow. The main result is a reversal of intuition: elastic stresses generally make the outlet hole smaller than the Newtonian prediction, and they are always destabilizing when inertia is negligible, unlike viscoelastic solid threads. A small regime—very large inlet holes or draw ratios close to one—does benefit from elasticity, preserving a larger hole. If the paper is right, polymer selection and draw-protocol design for holey fibres must be based on this sign change, not on the familiar pinching-suppression behaviour of elastic jets.","feed_headline":"Viscoelasticity closes tube holes faster, and destabilizes drawing","feed_subtitle":"Polymer elasticity shrinks internal holes and lowers stability margins in fibre drawing—opposite to intuition.","key_machinery":"The load-bearing machinery is a double asymptotic expansion: the usual slender-tube (long-wave) expansion in the aspect ratio $\\epsilon$, and then a low-Deborah-number expansion $\\psi = \\psi^{(0)} + De\\,\\psi^{(1)}$ around the Newtonian state. The leading-order system reduces to the Newtonian tube equations, and the first-order system determines the corrections $u^{(1)}$, $h^{(1)}$, and $H^{(1)}$. The central diagnostic identity is $h^{(1)}(z=1) = \\frac{1}{h^{(0)}(z=1) D} \\int_0^1 C^{(1)} dz$, which shows that the outlet hole correction is the integrated elastic correction to the radial-flow strength $C$. The paper rewrites $C^{(1)}$ in Eq. (4.7) as four terms: axial advection of the leading-order second normal stress difference, radial advection, a Giesekus mobility term, and a geometric response term; term 2 dominates for moderate holes, while term 1 dominates near the inlet for very large holes. The stability part uses the adjoint (Fredholm solvability) condition on the linearized perturbation equations to compute the elastic growth-rate correction $\\omega_1$.","core_discovery":"To first order in the Deborah number $De$ (the ratio of polymer relaxation time to device transit time), the elastic correction to the outlet hole radius, $h^{(1)}(z=1)$, has the opposite sign from what elasticity-driven pinching suppression would suggest. The sign is set by the elastic correction $C^{(1)}$ to the strength of the surface-tension-driven radial flow, and ultimately by the second normal stress difference induced by the Giesekus stresses. For a tube whose leading-order hole stays open, $C^{(1)}$ is typically negative, so the hole closes faster than in the Newtonian case; the opposite occurs only for large inlet hole ratio $\\varphi$ or draw ratio $D$ close to 1, where axial advection of stress near the inlet produces a strong outward flow. In the linear stability problem, the first-order elastic correction $\\omega_1$ to the growth rate is positive at the Newtonian critical draw ratio when $Re=0$, so elasticity always destabilizes; for non-zero inertia, $\\omega_1$ can be negative or positive depending on $Ca$, $\\varphi$, $\\alpha$, and $Re$.","pith_inferences":["Editorial inference: the paper's mechanism implies that rheological characterisation for holey-fibre polymers should prioritise the second normal stress difference, since hole evolution is independent of the first normal stress difference at this order.","Editorial inference: the sign reversal between solid threads and tubes suggests a crossover as hole size shrinks; testing very small $\\varphi$ (where the leading-order hole remains open) could locate where elastic destabilisation turns into stabilisation.","Editorial inference: direct numerical simulation of the full Giesekus equations for $De$ around 0.05 to 0.1 would test whether the long-wave first-order corrections remain accurate away from the asymptotic limit.","Editorial inference: for manufacturing, the paper implies a regime map in $(\\varphi, D, Re)$ where elasticity is beneficial; one could optimise toward large-hole, low-draw conditions while using inertia to restore stability."],"forward_implications":["For moderate inlet hole ratios, a weakly viscoelastic polymer melt will produce a smaller outlet hole than a Newtonian melt under the same draw conditions, so manufacturing compensation must account for elasticity.","For tubes with very large inlet holes and draw ratios close to one, elasticity can enlarge the outlet hole, offering a possible benefit for extrusion of large-air-hole microstructured fibres.","At negligible inertia, viscoelasticity lowers the stability margin: the critical draw ratio for draw resonance is reached earlier than the Newtonian prediction.","With inertia, the effect is parameter-dependent: elasticity stabilizes for sufficiently large capillary number or sufficiently small inlet holes, and destabilizes otherwise.","As the solvent fraction $\\beta$ approaches one, all elastic corrections vanish and the Newtonian tube-drawing behaviour is recovered, so the results interpolate cleanly between polymer and Newtonian processing."],"supporting_citations":[{"why":"Supplies the constitutive model (Giesekus stress with mobility factor alpha) used throughout the derivation.","marker":"Giesekus (1982)"},{"why":"Provides the Newtonian capillary-drawing long-wave model that the leading-order De=0 system recovers.","marker":"Fitt et al. (2002)"},{"why":"Gives the Newtonian stability analysis for drawing microstructured optical fibres, the baseline for the leading-order eigenvalue problem.","marker":"Wylie et al. (2023)"},{"why":"Establishes the pinching-suppression intuition for elastic liquids that the paper's hole-closure result overturns.","marker":"Entov and Hinch (1997)"},{"why":"Supports the same pinching-suppression intuition for viscoelastic jets, the background against which the sign of h^(1) is surprising.","marker":"Li and Fontelos (2003)"},{"why":"Shows elastic effects stabilize viscoelastic solid-thread spinning, the contrast case for the tube's instability result.","marker":"Zhou and Kumar (2010)"},{"why":"Documents the physical importance of the second normal stress difference, the quantity the paper identifies as the controlling mechanism.","marker":"Maklad and Poole (2021)"},{"why":"Supplies the weakly viscoelastic perturbation methodology (psi = psi^(0) + De psi^(1)) that the derivation follows.","marker":"Boyko and Stone (2022)"}],"fun_headline_variants":["Elasticity accelerates hole closure, destabilizes drawing","Viscoelasticity hastens hole closure and hurts stability","Elastic effects shrink tube holes faster, undermine stability","Elasticity closes tube holes faster, destabilizes fibre drawing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central results are first-order corrections in the Deborah number (relaxation time divided by transit time through the device), so the conclusions hold only when that number is small enough that the neglected quadratic corrections cannot reverse the sign of $h^{(1)}$ or $\\omega_1$, and only for inlet hole ratios large enough that the leading-order hole does not close (the paper sets aside $\\varphi < 0.5$ and strong surface tension).","fun_headline_variants_meta":{"raw":{"variants":["Elasticity accelerates hole closure, destabilizes drawing","Viscoelasticity hastens hole closure and hurts stability","Elastic effects shrink tube holes faster, undermine stability","Elasticity closes tube holes faster, destabilizes fibre drawing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000984,"raw_usage":{"total_tokens":4187,"prompt_tokens":967,"completion_tokens":3220,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":3154}},"tokens_in":583,"tokens_out":3220,"duration_ms":27313,"temperature":1.0,"reasoning_tokens":3154,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:30:32.339003+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a full two-dimensional numerical simulation of the Giesekus equations for $De \\approx 0.1$, $D \\approx 1.5$, $Ca \\approx 1.8$, $\\varphi \\approx 0.5$, $Re = 0$, $\\alpha = 0$, $\\beta = 0.1$: if the outlet hole radius is larger than the Newtonian value, the predicted sign $h^{(1)}(z=1) < 0$ is wrong. Separately, at the critical draw ratio for $Re=0$, compute the real part of the elastic growth-rate correction $\\omega_1$: a negative value would disprove the claim that elasticity is always destabilizing when inertia is negligible.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the constitutive model (Giesekus stress with mobility factor alpha) used throughout the derivation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Newtonian capillary-drawing long-wave model that the leading-order De=0 system recovers."},{"cited_title":"J., Papri, N","cited_arxiv_id":null,"evidence_quote":"Gives the Newtonian stability analysis for drawing microstructured optical fibres, the baseline for the leading-order eigenvalue problem."},{"cited_title":"and Hinch, E","cited_arxiv_id":null,"evidence_quote":"Establishes the pinching-suppression intuition for elastic liquids that the paper's hole-closure result overturns."},{"cited_title":"and Kumar, S","cited_arxiv_id":null,"evidence_quote":"Shows elastic effects stabilize viscoelastic solid-thread spinning, the contrast case for the tube's instability result."},{"cited_title":"and Poole, R","cited_arxiv_id":null,"evidence_quote":"Documents the physical importance of the second normal stress difference, the quantity the paper identifies as the controlling mechanism."},{"cited_title":"and Stone, H","cited_arxiv_id":null,"evidence_quote":"Supplies the weakly viscoelastic perturbation methodology (psi = psi^(0) + De psi^(1)) that the derivation follows."}],"review_version":1}