{"id":"c6f1fa62-6f5c-4fbe-8d9d-30993f8942ba","arxiv_id":"2608.10951","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The paper derives a reduced two-coupled Ginzburg-Landau system that governs the leading-order nonlinear dynamics of the non-axisymmetric primary instability in small-gap Taylor-Couette flow for strong counter-rotation, and uses it to characterize helicoidal and ribbon waves.","lead":"This mathematics paper derives a pair of coupled Ginzburg-Landau equations that describe how Taylor-Couette flow between strongly counter-rotating cylinders becomes unstable near the onset of non-axisymmetric waves. The result offers a reduced model for a regime that previous analytic approaches did not cover, potentially helping to interpret experiments and numerical simulations in this small-gap, high-Reynolds limit.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's central claim that the Ginzburg-Landau system describes the original Navier-Stokes dynamics rests on an omitted Step 3, delegated to reference [9], and the paper does not verify that the spectral and nonlinear hypotheses of that reduction theory hold for the limit system (7).","rationale":"The reader's verdict is CONDITIONAL, and my stress-test supports keeping that verdict. The paper is a serious formal derivation: the limit-system expansion in Appendix 4.1 is coherent, the linear stability analysis is standard, and the formulas (38)-(41) for the cubic coefficients are explicit enough to be reproduced. The main theorem (Theorem 1) about solutions of the GL system is likely correct as algebra. What is not established is the logical bridge from the GL system to genuine solutions of the Navier-Stokes equations. The paper flags this itself by omitting Step 3 and citing [9]; under the review rule that self-flagged limitations are in-scope evidence, this omission is the weakest load-bearing point. A referee can accept the reduction as a formality only if the hypotheses of [9] are obviously met; here they are not checked, and the unbounded y-direction makes the spectral-gap question nontrivial. I agree with the reader's weakest_assumption. I do not see an internal algebraic contradiction in the coefficient formulas, so I would not move to REJECT. The numerical convergence issue, also noted by the reader, is real but secondary; it can be addressed by a convergence table. Therefore the verdict stays CONDITIONAL: accept only after Step 3 is either justified by checking [9]'s hypotheses or replaced by a direct numerical validation of the GL prediction against the limit system.","tokens_in":19714,"tokens_out":33998,"duration_ms":324539,"concrete_test":"As a check, carry out the omitted Step 3 for the specific operator in Section 4.4: verify the hypotheses of the reduction theorem of [9] for the linearized operator L0 defined in (31) at T=T_c, α=α_c, B=B_c. Concretely, prove that L0-iω0 and L0-2iω0 are invertible with bounded inverses on the relevant space of functions periodic in z with period 2π/α_c and slowly varying in y, and that Re σ(L0) ≤ -γ < 0 except for the neutral pair ±iω0. If the spectral gap or the nonlinear estimates fail because of the continuous y-spectrum, then the GL system is not proven to govern the original system; if they succeed, the central claim is anchored. A complementary numerical test would be a direct simulation of (7) at μ=-0.9 with T slightly above T_c and a small A/B-mode initial condition, comparing the envelope evolution with (5)-(6).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper states in Section 1: 'This last step can be done by repeating the proofs in [9], up to minor changes. We will omit this step here and refer to [9] for more details.' This is the load-bearing premise of the entire paper: without Step 3, the solutions of the coupled Ginzburg-Landau system (5)-(6) are only solutions of a formally derived amplitude equation, not of the limit Navier-Stokes system (7), let alone of the original Navier-Stokes equations. The paper does not check the hypotheses of the reduction theory in [9] for the present problem. In particular, the linearized operator around the critical point has neutral modes at ±iω0 with wavenumbers (±α_c, β_c), and because the problem is translation invariant in the unbounded y-direction, the spectrum near B=B_c accumulates on the imaginary axis when modulations in y are allowed; a standard center-manifold theorem requires a spectral gap, and the paper does not explain how the spatial-dynamics formulation of [9] supplies that gap. It also does not verify that the nonlinear operator (34) satisfies the required smoothness and growth estimates in the chosen function space. The numerical evaluation of the coefficients b and c, while necessary for the stability conclusions, is secondary: even perfectly accurate coefficients cannot turn a formally derived GL system into a rigorous description of the original PDE if the reduction step is not justified. The qualitative stability regions in Sections 3.2-3.3 would be groundless if Step 3 fails, so this omitted verification is the weakest load-bearing point.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Couette-Taylor instability in the small-gap, slow-rotation, high-Reynolds limit for strongly counter-rotating cylinders (mu < mu_c ~ -0.8), where the primary instability is non-axisymmetric. From the Navier-Stokes equations the authors recall the limit system (7), perform its linear stability analysis, and derive a system of two coupled complex Ginzburg-Landau equations, (5)-(6), for the amplitudes of the left- and right-travelling azimuthal waves. All linear and cubic coefficients are expressed through eigenfunctions and the adjoint problem, and the cubic coefficients b and c are evaluated numerically by Chebyshev collocation. The reduced equations are then used to discuss existence and stability of helicoidal and ribbon waves and to derive a third-order ODE system for spatially modulated solutions. The main new content is the two-amplitude Ginzburg-Landau description, the explicit coefficient formulas in Appendix 4.4, and the resulting stability diagram based on the signs of b_r, c_r and b_r+c_r.","tokens_in":1656,"tokens_out":2729,"duration_ms":201937,"significance":"If the reduction is valid, the paper provides a systematic weakly nonlinear description in a regime that previous small-gap analyses did not cover, and it gives explicit, checkable formulas for the amplitude-equation coefficients in terms of the linearized eigenfunctions and adjoint problem. The stability criteria for helicoidal and ribbon waves are presented transparently, and the reduction to a third-order ODE for spatially modulated states is a useful starting point. The paper is nevertheless conditional in an essential way: the step that lifts solutions of the Ginzburg-Landau system to genuine solutions of the limit Navier-Stokes system is delegated to reference [9] and not proved, and the numerical coefficients that determine the stability thresholds are reported without a convergence study. These two issues are load-bearing for the paper's central claims.","major_comments":[{"comment":"The central claim that the Ginzburg-Landau system (5)-(6) describes the limit Navier-Stokes system (7) rests on the omitted Step 3. The sentence 'This last step can be done by repeating the proofs in [9], up to minor changes' is not sufficient as it stands: the manuscript does not verify that the hypotheses of the reduction theory in [9] hold for the present problem, which is time-dependent, has two coupled complex amplitudes, and is posed on the unbounded y-axis where the linearized spectrum is not gapped when modulations are allowed. Please either provide a reduction theorem for (7) with the spectral-gap and nonlinearity hypotheses stated and verified, or explicitly restrict the claims to the formally derived amplitude system. As written, the existence and stability statements for helicoidal and ribbon waves in Sections 3.2 and 3.3 are theorems only about the reduced system unless Step 3 is supplied.","section":"Section 1, Step 3"},{"comment":"The stability thresholds mu_1, mu_2, mu_3 and the signs of b_r, c_r and b_r+c_r are computed numerically, but the paper contains no convergence study in the number N of Chebyshev collocation points, no error estimates, and only one reported numerical value (b and c at mu=-1). The statement that extended precision guarantees well-converged results is not a substitute. Please include a table of b and c for increasing N at representative values of mu, and preferably an independent check of the crossing points, so that the stability diagram in Figure 1 is reproducible and the thresholds are credible.","section":"Section 4.4.3 and Figure 1"},{"comment":"The perturbation expansion of U lists U_0, U_100, U_010, U_001, U_110, U_011 and U_020, but the coefficient U_200 is never introduced; nevertheless U_200 appears in the displayed equation for a_4. This makes the formula for a_4 unverifiable as written. Please include the (B-B_c)^2 term in the expansion and define U_200, or restructure the calculation so that every coefficient used is defined.","section":"Section 4.3, expansion near (29)"},{"comment":"The text states Phi_0011 = S Phi_1100 and then concludes Phi_0011 = (0, phi_1100_y, -phi_1100_z)^t = Phi_1100. Since phi_1100_y is real and phi_1100_z is pure imaginary, the vector (0, phi_1100_y, -phi_1100_z) is the complex conjugate of Phi_1100, not Phi_1100 itself. If the intended identity is Phi_0011 = conjugate of Phi_1100, this should be stated explicitly, because the sign convention enters the c-coefficient formula (41).","section":"Section 4.4.1, computation of Phi_0011"}],"minor_comments":[{"comment":"The notation for the azimuthal wavenumber is overloaded: B is the rescaled variable B=beta R, B_c is the critical value of B, while zeta_1 contains e^{i beta_c y}; later beta is used for a modulation wavenumber in either y or tilde y. Please make the relation between beta, B, R, y and tilde y consistent, especially in Theorem 1 and equation (14), where y and tilde y appear to be used interchangeably.","section":"Sections 2.2, 3.1, 3.2 and Theorem 1"},{"comment":"The value of mu_c is given as approximately -0.8 in the abstract and introduction and -0.785 in Section 3.2; please unify the values and state the numerical precision.","section":"Throughout"},{"comment":"The caption describes crossings at mu_3, mu_1 and mu_2, but the figure does not show markers or a table of the computed b_r(mu), c_r(mu) values. Adding a small table or markers would substantially improve verifiability.","section":"Figure 1"},{"comment":"The scalar product (22) involves only the x and y components; the paper should state explicitly that the z-component is eliminated through the divergence-free condition, and that the nonzero normalization <zeta_1, zeta_1^*> is checked numerically.","section":"Section 4.3 and Appendix 4.4"},{"comment":"There are typos such as 'exist on the super-critical size' which should read 'side'; similar small typographical errors occur elsewhere and should be corrected.","section":"Section 3.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is a reasonable contribution to the Taylor-Couette amplitude-equation literature, and the explicit coefficient formulas are valuable. The main risk is the omitted rigorous reduction step: the authors should either prove the reduction for the limit system (7) or clearly state the result conditionally on the hypotheses of [9]. The numerical convergence data are essential because the stability diagram depends quantitatively on the computed signs of b_r, c_r and b_r+c_r. With those additions the paper would be acceptable; without them the central claim remains a formal asymptotic statement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the first amplitude-equation treatment of the strongly counter-rotating small-gap Couette-Taylor regime, where the primary instability is non-axisymmetric. The authors derive a two-coupled Ginzburg-Landau system, compute its coefficients from the linear eigenfunctions and adjoint problem via Chebyshev collocation, and use the reduced system to map out helicoidal-wave and ribbon-wave stability, plus some exotic solutions governed by a third-order ODE. Appendix 4.4 is a serious piece of coefficient bookkeeping; the symmetries (zeta2 = S zeta1, the quadratic interactions) are handled cleanly, and the sample values at mu = -1 are concrete. The paper is also honest that the exotic solutions are not fully classified.\n\nThe soft spots are real but not fatal. The load-bearing issue is the omitted Step 3. The paper says the lift from solutions of (5)-(6) to genuine solutions of the limit Navier-Stokes system can be obtained by repeating [9] 'up to minor changes.' That might be true, but it is not demonstrated. The original result in [9] does not obviously cover two coupled GL equations with both axial and azimuthal modes and y-modulation; the spectral-gap and nonlinear-estimate hypotheses are neither checked nor cited as a ready-made theorem covering this exact case. Until that step is supplied, (5)-(6) is a formally derived envelope system, not a proven reduction. This is exactly what a referee should ask for.\n\nSecond, the stability regions in Sections 3.2-3.3 depend on numerically computed coefficients b_r and c_r, and on the crossing values mu1, mu2, mu3. There is no convergence study in N and no error estimate, only a claim of extended precision. That is a moderate issue, not a fatal one; the formulas are explicit and reproducible, but the paper should ship data or a convergence table.\n\nThe citation pattern is fine. [1,2] are the same group's prior related results and are cited as such; [9] is a sensible source for the reduction strategy; the numerical work of Deguchi and Nagata is acknowledged.\n\nWho should read it: anyone working on Taylor-Couette amplitude equations or on formal center-manifold reductions in hydrodynamic stability. It deserves a serious referee, but the referee should press on the reduction step and the numerics. My verdict would be conditional acceptance: require the authors to verify that [9] applies, or state precisely which theorem does, and add a convergence study for the coefficients.","headline":"A serious first amplitude-equation treatment of the strongly counter-rotating small-gap Couette-Taylor regime, with a real gap: the center-manifold step is delegated, not verified, and the GL coefficients lack convergence data.","tokens_in":20582,"tokens_out":2436,"would_cite":true,"duration_ms":24855,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B32","76E07","76D05"],"pacs":[],"model":"deepseek-v4-flash","headline":"For strongly counter-rotating cylinders, the first Taylor-Couette instability becomes non-axisymmetric and is governed by two coupled Ginzburg-Landau equations.","keywords":["Couette-Taylor instability","small-gap limit","counter-rotating cylinders","non-axisymmetric instability","coupled Ginzburg-Landau equations","helicoidal waves","ribbon waves","center-manifold reduction"],"falsifier":"Compute the coefficients $b$ and $c$ (or their real parts) for a fixed $\\mu$ in $(-0.85, -0.8)$ using an independent numerical method, or run direct Navier-Stokes simulations in the SGSRRHR limit to measure the critical Taylor number and the stability of helicoidal versus ribbon waves; if the signs of $b_r$ and $c_r$ (or the predicted thresholds $\\mu_1$, $\\mu_2$, $\\mu_3$) do not match the paper's predictions, the central claim fails.","tokens_in":19434,"feed_emoji":"🌀","tokens_out":8707,"duration_ms":71975,"temperature":0.7,"pith_summary":"For a viscous flow between two cylinders that counter-rotate sufficiently strongly (rotation-rate ratio $\\mu$ below $\\mu_c \\approx -0.8$), the first Couette-Taylor instability is no longer axisymmetric but takes the form of azimuthal travelling waves. The paper shows that in the small-gap, slow-rescaled-rate, high-Reynolds limit, the slow modulation of these left- and right-travelling waves is governed by a system of two coupled complex Ginzburg-Landau equations. It computes every coefficient of this amplitude system, including the cubic nonlinearities, numerically from the linearised eigenfunctions and the adjoint eigenvector. It then determines when helicoidal waves (travelling in both axial and azimuthal directions) and ribbon waves (standing axially, travelling azimuthally) exist and are linearly stable, and identifies a third-order ODE system for more exotic spatially modulated solutions whose classification is left open.","feed_headline":"Strong counter-rotation switches Couette-Taylor onset to travelling waves","feed_subtitle":"Near the critical Taylor number, two coupled Ginzburg-Landau equations decide between helicoidal and ribbon wave states.","key_machinery":"The central object is the limit Navier-Stokes system (7), obtained in the SGSRRHR limit, together with its two-dimensional critical eigenspace spanned by the eigenfunctions $\\zeta_1 = e^{i(\\alpha_c z + \\beta_c y)} U_1(x)$ and its reflected counterpart $\\zeta_2 = S\\zeta_1$. The mechanism carrying the argument is the center-manifold reduction of this limit system to the two coupled complex Ginzburg-Landau equations, whose symmetry group forces the normal form (5)-(6). The coefficients $b$ and $c$ are the load-bearing numbers: they are computed from the quadratic interaction operator $B$ and the adjoint eigenvector via the Fredholm-type formulas (40)-(41), and their real parts determine the existence and stability of helicoidal waves (one amplitude nonzero) and ribbon waves (equal amplitudes, standing axially).","core_discovery":"The paper's central claim is that, in the SGSRRHR regime with $\\mu < \\mu_c$, small-amplitude solutions of the limit Navier-Stokes system near the critical Taylor number $T_c(\\mu)$ are described at leading order by the amplitude equations (5)-(6), namely two coupled complex Ginzburg-Landau equations for the amplitudes $A$ and $B$ of the two non-axisymmetric critical modes. All coefficients of this system—the linear dispersion terms and the cubic coupling coefficients $b$ and $c$—are obtained from the eigenfunctions $\\zeta_1, \\zeta_2$ and the adjoint eigenfunction via explicit formulas (40)-(41), evaluated by Chebyshev collocation. The reduced system admits helicoidal-wave and ribbon-wave solutions, and their existence and stability are governed by the signs of the real parts $b_r, c_r$, $b_r + c_r$, and $b_r - c_r$, with thresholds at $\\mu_1 \\approx -0.814$, $\\mu_2 \\approx -0.848$, and $\\mu_3 \\approx -0.8$. The paper also reduces the search for stationary spatial modulations to a third-order ODE system (18). If these claims hold, they give a systematic weakly nonlinear description of the non-axisymmetric onset in the small-gap limit.","pith_inferences":["If the reduction holds, the same two-mode Ginzburg-Landau structure should appear in other hydrodynamic systems with a reflection symmetry coupling two counter-propagating waves, suggesting the stability criteria derived here are generic rather than specific to the Taylor-Couette geometry.","The threshold values $\\mu_1 \\approx -0.814$, $\\mu_2 \\approx -0.848$, and $\\mu_3 \\approx -0.8$ could be compared with finite-gap experiments or direct numerical simulations at moderate radius ratios; agreement would support the small-gap limit as a quantitative model beyond its formal derivation.","The third-order ODE system (18) is likely to admit spatially periodic or homoclinic orbits that correspond to modulated wave packets; proving or disproving their existence would settle the open classification problem raised in the paper.","A direct check of the paper's coefficients using a different discretization (for instance finite elements or shooting methods) would test the reliability of the Chebyshev-collocation results, since no convergence study in the number of collocation points is reported."],"forward_implications":["For $\\mu < \\mu_c \\approx -0.8$, the first linear instability is non-axisymmetric, with critical axial wavenumber $\\alpha_c(\\mu)$ and azimuthal wavenumber $B_c(\\mu) > 0$, so the onset branch differs fundamentally from the axisymmetric Taylor-vortex branch.","The coupled Ginzburg-Landau system predicts that helicoidal waves exist and are stable for $\\mu_2 < \\mu < \\mu_1$ (with $\\mu_1 \\approx -0.814$ and $\\mu_2 \\approx -0.848$), while ribbon waves are stable for $\\mu < \\mu_2$, with the intervals determined by the signs of $b_r$, $c_r$, $b_r + c_r$, and $b_r - c_r$.","If the coefficients are correct, the weakly nonlinear behaviour near onset for any fixed $\\mu$ in the strongly counter-rotating regime is fully captured by the two-equation amplitude system, so the competition between helicoidal and ribbon states is decided by the ratio $c_r/b_r$.","Provided the lifting to Navier-Stokes is justified via the cited reduction theory, each periodic solution of the reduced system corresponds to a genuine periodic flow in the small-gap limit, giving explicit travelling-wave and standing-wave patterns."],"supporting_citations":[{"why":"Establishes the small-gap limit system for the co-rotating case, whose derivation is recalled here.","marker":"[1]"},{"why":"Extends the small-gap analysis to counter-rotating cylinders and identifies the threshold $\\mu_c$ where the primary instability becomes non-axisymmetric.","marker":"[2]"},{"why":"Supplies the general coupled Ginzburg-Landau framework and the stability analysis of helicoidal and ribbon waves used as the template.","marker":"[3]"},{"why":"Provides the center-manifold reduction theory that justifies lifting solutions of the reduced system to solutions of the Navier-Stokes equations; the paper omits this step and refers to this theory.","marker":"[9]"},{"why":"Provides the narrow-gap limit asymptotics for Couette flow that the small-gap limit builds on.","marker":"[13]"},{"why":"Studies the narrow-gap limit at $\\mu = -1$ without considering non-axisymmetric bifurcations, giving the comparison point for the present analysis.","marker":"[14]"}],"fun_headline_variants":["Very counter-rotating cylinders shift Couette-Taylor instability to waves","Non-axisymmetric onset explained: coupled Ginzburg-Landau equations","Helicoidal and ribbon waves take over in strong counter-rotation","How strong counter-rotation changes Couette-Taylor flow's first instability","Couette-Taylor: when counter-rotation turns onset into travelling waves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the center-manifold reduction to the two coupled Ginzburg-Landau equations faithfully represents the full Navier-Stokes system near onset in the SGSRRHR limit, so that the numerically computed coefficients are the true reduction coefficients and the reduced solutions lift to genuine solutions of the original equations.","fun_headline_variants_meta":{"raw":{"variants":["Very counter-rotating cylinders shift Couette-Taylor instability to waves","Non-axisymmetric onset explained: coupled Ginzburg-Landau equations","Helicoidal and ribbon waves take over in strong counter-rotation","How strong counter-rotation changes Couette-Taylor flow's first instability","Couette-Taylor: when counter-rotation turns onset into travelling waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000457,"raw_usage":{"total_tokens":2649,"prompt_tokens":1033,"completion_tokens":1616,"prompt_tokens_details":{"cached_tokens":1024},"prompt_cache_hit_tokens":1024,"prompt_cache_miss_tokens":9,"completion_tokens_details":{"reasoning_tokens":1521}},"tokens_in":9,"tokens_out":1616,"duration_ms":36165,"temperature":1.0,"reasoning_tokens":1521,"cache_read_input_tokens":1024,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:32:19.761682+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the coefficients $b$ and $c$ (or their real parts) for a fixed $\\mu$ in $(-0.85, -0.8)$ using an independent numerical method, or run direct Navier-Stokes simulations in the SGSRRHR limit to measure the critical Taylor number and the stability of helicoidal versus ribbon waves; if the signs of $b_r$ and $c_r$ (or the predicted thresholds $\\mu_1$, $\\mu_2$, $\\mu_3$) do not match the paper's predictions, the central claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the small-gap limit system for the co-rotating case, whose derivation is recalled here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Extends the small-gap analysis to counter-rotating cylinders and identifies the threshold $\\mu_c$ where the primary instability becomes non-axisymmetric."},{"cited_title":"Chossat, G","cited_arxiv_id":null,"evidence_quote":"Supplies the general coupled Ginzburg-Landau framework and the stability analysis of helicoidal and ribbon waves used as the template."},{"cited_title":"Iooss, A","cited_arxiv_id":null,"evidence_quote":"Provides the center-manifold reduction theory that justifies lifting solutions of the reduced system to solutions of the Navier-Stokes equations; the paper omits this step and refers to this theory."},{"cited_title":"Nagata, Taylor-Couette flow in the narrow-gap limit,Phil","cited_arxiv_id":null,"evidence_quote":"Provides the narrow-gap limit asymptotics for Couette flow that the small-gap limit builds on."},{"cited_title":"Nagata, Taylor-Couette system in the narrow-gap limit, revisited, with a corrigendum to the paper by Nagata (2023, Philosophical Transaction A.), Phil","cited_arxiv_id":null,"evidence_quote":"Studies the narrow-gap limit at $\\mu = -1$ without considering non-axisymmetric bifurcations, giving the comparison point for the present analysis."}],"review_version":1}