{"id":"28c17cac-b25e-4916-b02f-a5b4291141a1","arxiv_id":"1908.04713","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The authors derive closed-form eddy-current formulas for a spinning brass disk or ring in an anti-symmetric magnetic field and match a measured spin-down time to within 8.61% after calibrating friction.","lead":"This paper derives a closed-form formula for the eddy currents that slow a spinning metal gyroscope above a pair of magnets, and it checks the predicted slowdown against a strobe-light experiment. The method reduces the electromagnetic problem to a simple Poisson equation, which makes it useful for teaching and for quick engineering estimates of eddy-current damping.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (30) as printed violates charge conservation: ∇·j1 = σωk1 a²b²/(4r³) sinθ ≠ 0, so the central closed-form current is not a solution of the stated boundary value problem.","rationale":"The paper's quasi-static reduction to a Poisson equation is a standard and plausible approach, and the overall logic of computing eddy currents from ∇·j = 0 with the j·n = 0 boundary condition is sound. The strongest_claim, however, rests on Equations (30) and (31) being correct closed-form solutions. A direct check shows that Equation (30) does not satisfy charge conservation in the ring geometry: the printed θ-component has the wrong sign for the a²b²/r² term. This is a concrete, text-level inconsistency in the central analytical result, not merely a disagreement with a numerical benchmark. Because the spin torque depends only on j_r, the headline experimental comparison (8.61%) may be unaffected by this particular sign error, and the error may be a typographical slip rather than a conceptual failure. That is why I recommend CONDITIONAL rather than REJECT: the paper should not be accepted until the formulas are corrected and independently verified. A separate concern about solving the disk and rings as isolated regions without interface continuity is also present and reinforces the need for an independent numerical or analytical check of the assembled rotor. The reader's weakest_assumption focused on the fitted magnetic-field polynomial and the neglected in-plane field; that is a valid concern about validation strength, but the more load-bearing issue for the central claim is the internal inconsistency of the printed current formula. Hence partial agreement with the reader. The paper deserves credit for a clear derivation outline, explicit boundary conditions, and an experimental comparison, but with no machine-checked proof, no shipped code, and displayed formulas that fail a basic consistency test, the closed-form result cannot be accepted as-is.","tokens_in":10089,"tokens_out":16348,"duration_ms":168990,"concrete_test":"Recompute the divergence of the printed j1 in Equation (30) in cylindrical coordinates. The calculation gives ∇·j1 = σωk1 a²b²/(4 r³) sinθ, which is nonzero. Then re-solve the Poisson/Neumann problem of Section V for the potential φ(r,θ) = R(r) sinθ with the stated source and boundary conditions, and form j = σ(−∇φ + f); the resulting θ component is σωk1/8 (a²+b² + a²b²/r² − 3r²) cosθ. If this re-derivation confirms the sign discrepancy, Equation (30) must be corrected. Repeat the same consistency check for Equation (31): if its divergence is also nonzero, the closed-form current claim fails in full. For the interface concern, solve the same Poisson problem on the connected disk-plus-ring geometry with continuity conditions at the shared radius and compare the total torque with the sum of the separately solved regions.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim is a closed-form eddy-current solution, Equations (30) and (31), from which the damping torque is computed. As printed, Equation (30) is internally inconsistent. Substituting the displayed j1(r,θ) into ∇·j = 0 in cylindrical coordinates gives ∇·j1 = σωk1 a²b²/(4 r³) sinθ, which is nonzero wherever the ring has a finite inner radius. The θ-component's a²b²/r² term carries the wrong sign: independent solution of the Poisson/Neumann problem stated in Section V yields j1θ = σωk1/8 (a²+b² + a²b²/r² − 3r²) cosθ, with a plus sign before a²b²/r², not the minus sign printed. The radial component, and hence the spin torque M = −λω which depends only on j_r, appears to be correct, so the reported 8.61% comparison may survive a sign correction; but the abstract's claim of a closed-form current distribution is not established as written. A second, independent concern is that the disk and rings are solved as separate boundary-value problems with j·n = 0 at r = a and r = b, without imposing continuity of potential or normal current at the interfaces of a single electrically connected brass rotor. If these regions are in electrical contact, the assembled conductor's current distribution differs from the sum of independent solutions, so even a corrected (30)–(31) may not describe the actual gyroscope.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents an analytical treatment of eddy-current damping in a spinning brass gyroscope. Working in the lab frame and assuming a static magnetic field and quasi-static currents, the authors reduce the problem to a Poisson equation with Neumann boundary conditions, approximate the vertical field component over the rotor as B_z = k1 z + k3 z^3, solve the boundary-value problem for a disk/ring geometry, and obtain closed-form eddy-current expressions. They then compute the magnetic torque as −λω, predict an exponential decay of the rotation frequency, and compare the decay time with strobe-light measurements, reporting a relative error of 8.61%.","tokens_in":10388,"tokens_out":19963,"duration_ms":163627,"significance":"If correct, the paper would offer a compact, pedagogical alternative to series solutions for eddy currents in thin rotating conductors, and the experimental comparison would provide a useful validation. The core quasi-static reduction (Eq. 14) and the divergence calculation in Eq. (25) are standard and check out. However, the printed current formula (30) fails charge conservation as written, the treatment of the composite rotor as independent regions is not physically justified, and the experimental comparison involves calibrated or fitted inputs; these issues currently limit the significance of the paper.","major_comments":[{"comment":"The printed θ-component of j1 contains a minus sign in front of a²b²/r². Direct solution of the stated boundary-value problem (25)–(28) gives j1θ = (σωk1/8)(a² + b² + a²b²/r² − 3r²) cosθ, with a plus sign before the a²b²/r² term. With the published minus sign, ∇·j1 = σωk1 a²b²/(4r³) sinθ, which is nonzero wherever a > 0, so Eq. (30) is not a solution of the charge-conserving problem. The radial component appears correct, so the torque may survive a sign correction, but the central claim of a closed-form current distribution is not correct as printed. Eq. (31) should be rechecked with the same divergence test.","section":"Section V, Eq. (30)"},{"comment":"The plate and the rings are solved as independent boundary-value problems, each with j_r = 0 at r = a and r = b. If these are parts of a single electrically connected brass rotor, as the description around Fig. 5 suggests, those surfaces are internal conducting interfaces rather than insulating boundaries; the true current distribution must satisfy continuity of the potential and of the normal current across them. The manuscript provides no justification for treating the regions as isolated, so even after correcting Eq. (30), the closed-form currents may not describe the actual assembled rotor. The authors should state whether the regions are electrically isolated or solve the coupled problem.","section":"Section V, Eq. (25) and Fig. 5"},{"comment":"The torque integral leading to τ = 7.839 s is not shown. The paper says 'we can easily calculate' the integral in Eq. (32), but it gives neither an expression for λ nor the integration steps, so the theoretical decay time cannot be reproduced or checked. For a paper whose central quantitative claim is the 8.61% agreement in τ, the derivation of λ as a function of a, b, σ, k1, and k3 should be provided explicitly.","section":"Section VI, Eqs. (32)–(36)"},{"comment":"The reported 8.61% agreement is not a parameter-free prediction. In Section IV, Br is described as a fitting parameter and k1, k3 are fitted to the field model; in Section VI, a friction multiplier is fitted so that the theoretical constant term matches the experimental curve. Because the empirical friction torque (36) contains both a constant term M0 and a linear term αω, the fitted multiplier also modifies the effective decay rate (λ + mα)/I, so τ is partly determined by the fit. The paper should state explicitly which quantities are predicted and which are calibrated, and should quantify the sensitivity of τ to the friction multiplier. In addition, the final paragraph of Section VII estimates the neglected in-plane field contribution as about 10% in torque, which is larger than the reported 8.61% discrepancy; this error budget needs a quantitative treatment before the agreement can be interpreted as a validation of the model.","section":"Section VI and Section VII"}],"minor_comments":[{"comment":"The same vertical component of the magnetic field is denoted both By and Bz in Eq. (23); please fix the coordinate notation so the axes and the field-component subscripts are consistent.","section":"Eq. (23)"},{"comment":"The sentence 'Here we make use of cylindrical coordinates (or polar coordinates since the problem is 2D) them' is ungrammatical and should be rewritten.","section":"Section V"},{"comment":"The radial component of Eq. (30) is written as a sum of two terms; combining them into a single fraction would make it easier to see that the boundary conditions j_r(a) = j_r(b) = 0 are satisfied.","section":"Section V, Eq. (30)"},{"comment":"The text states 'Entering the dimensions and the conductivity of our gyroscope, we get the result that τ = 7.839s' but does not give the conductivity value or the rotor dimensions used in the calculation; please include these inputs.","section":"Section VI"},{"comment":"The title of reference [9] contains a typo: 'Three-Dimentional' should be 'Three-Dimensional'.","section":"Reference [9]"},{"comment":"The phrase 'top-down anti-symmetry' is unclear; presumably 'top-bottom anti-symmetry' is intended.","section":"Fig. 6 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper appears to be a well-intentioned student project with a standard and mostly correct core derivation, but the sign error in the printed current formula, the unaddressed interface condition for the composite rotor, and the fitted nature of the experimental comparison are serious issues that prevent acceptance in the current form. The journal should require a corrected derivation, an explicit expression for λ, and an honest account of which inputs are calibrated before considering publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The paper actually does something new: it solves the eddy-current problem in the lab frame for a thin ring or disk in an anti-symmetric polynomial vertical field, giving closed-form currents and a damping time. Smythe and Schieber worked in the rotating frame; Nurge's sphere is a series. A closed-form lab-frame answer for this geometry is useful and teachable. The quasi-static reduction to ∇²φ = ∇·[(ω×r)×B] is clean and standard, and the radial current—the component that produces the magnetic torque—checks out.\n\nThe problem is that the printed current distribution is not a solution of the stated BVP. In Eq(30), the θ-component has the wrong sign on the a²b²/r² term. Substitution gives ∇·j1 = σωk1 a²b²/(4r³) sinθ ≠ 0, so charge is not conserved. The correct coefficient has a plus sign, and with that the divergence vanishes. Since the torque depends only on j_r, the 8.61% comparison may survive, but the abstract's claim of a closed-form current distribution is false as written.\n\nThe experimental validation is also weaker than it looks. Br, k1, k3, and a friction multiplier are fitted, so the 8.61% agreement is calibrated, not a free prediction. The uncalibrated decay time is off by roughly half. Credit where due: the strobe-light measurement is clever, and the authors explicitly estimate the neglected in-plane field at about 10% error.\n\nA second modeling concern: the rotor is a plate plus two rings, solved as independent BVPs with j·n=0 at the interfaces. If the brass parts are in electrical contact, the assembled current distribution is not the sum of independent solutions; you need continuity of potential and normal current. The paper never addresses this.\n\nOverall, the derivation is mostly sound, the writing is clear, and the pedagogical value is real. But the central formula needs a sign correction and the interface condition needs justification. I would send it to peer review—it deserves a serious referee—but I would demand a corrected version. As it stands, I would not cite Eq(30).","headline":"A clean quasi-static derivation of eddy-current damping that prints the wrong sign in one current component; the torque and the 8.61% comparison survive, but the central formula is not a solution as written.","tokens_in":10938,"tokens_out":10261,"would_cite":false,"duration_ms":90888,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Eddy-current damping in a spinning gyroscope is captured by a Poisson equation whose closed-form solution matches the measured slowdown to within 8.61%.","keywords":["eddy current","electromagnetic damping","gyroscope teslameter","Poisson equation","closed-form solution","magnetic braking torque","quasi-static approximation"],"falsifier":"Measure the horizontal field component in the rotor plane with a calibrated Hall probe and include its torque in the equation of motion; if the decay-time gap closes from 8.61 percent to near zero, the neglect is confirmed, whereas a sizeable remaining gap would implicate the cubic polynomial fit or the quasi-static assumption.","tokens_in":9838,"feed_emoji":"🧲","tokens_out":4883,"duration_ms":50682,"temperature":0.7,"pith_summary":"Eddy currents in a spinning conductor are usually hard to compute because they couple electric and magnetic fields in time. This paper shows that, for a thin rotating brass gyroscope under quasi-static conditions, the problem collapses to a static boundary-value problem for a scalar potential, namely a Poisson equation whose source term is set by the motion of the rotor through the applied field. For the anti-symmetric vertical field $B_y=k_1 z + k_3 z^3$ produced by two cuboid magnets, the induced current distribution has a closed form, and the magnetic braking torque is exactly proportional to the angular velocity. The resulting exponential slowdown, supplemented by a linear-plus-constant frictional torque, reproduces the strobe-light measurement, with the characteristic decay time off by 8.61 percent. The significance is practical: a problem normally requiring numerical simulation becomes a calculation that students and engineers can do with a desktop computer.","feed_headline":"Closed-form eddy currents match gyroscope damping to 8.61%","feed_subtitle":"A Poisson-equation solution for a spinning brass disk in a cubic magnetic field yields the measured slowdown curve.","key_machinery":"The mechanism carrying the argument is the reduction of Faraday's law plus charge conservation to a Poisson equation for a scalar electric potential $\\phi$, with the boundary condition that no current leaves the conductor. Because the rotator is thin and the applied field is taken to have only a $y$ component that depends on the in-plane coordinate $z$, the problem becomes two-dimensional in polar coordinates. Separation of variables through angular eigenfunctions $\\phi_n=A_n(r)\\cos(n\\theta)+B_n(r)\\sin(n\\theta)$ turns the equation into Euler ordinary differential equations, and the particular-plus-harmonic superposition produces the closed-form currents in Eqs. (30) and (31). The same machinery delivers a torque linear in angular velocity and hence the exponential decay law that is compared with experiment.","core_discovery":"The paper's central claim is that the eddy-current problem in a thin spinning conductor in a static magnetic field can be reduced, without solving the time-dependent Maxwell equations, to the boundary value problem $\\nabla^2\\phi = \\nabla\\cdot[(\\boldsymbol{\\omega}\\times\\mathbf{r})\\times\\mathbf{B}]$ with $\\mathbf{j}\\cdot\\hat{\\mathbf{n}}=0$ on the surface. Solving this Poisson equation in polar coordinates, by splitting the solution into a particular part and a harmonic part, yields closed-form current densities $\\mathbf{j}_1$ and $\\mathbf{j}_3$ for the linear and cubic terms of the vertical field. The magnetic torque from these currents is $M=-\\lambda\\omega$, so the rotational equation $I\\,d\\omega/dt=-\\lambda\\omega$ predicts $\\omega(t)=\\omega_0 e^{-t/\\tau}$. After calibrating the magnet remanence and one frictional multiplier, the computed decay time differs from the experiment by 8.61 percent, and the authors trace most of this residual to their neglect of the horizontal field component.","pith_inferences":["A natural follow-up experiment is to place the magnets so that the cubic term in $B_y$ is negligible, leaving only the linear term; the decay should then be governed entirely by Eq. (30), which would isolate that formula's accuracy.","The same scalar-potential reduction should extend to other thin rotors such as annular plates or disks of nonuniform thickness, provided the boundary conditions on each edge are imposed separately.","If the horizontal-field torque were measured or computed independently, the reported 8.61 percent discrepancy should shrink to nearly zero, making the neglect of that component a testable quantitative claim rather than a qualitative error estimate."],"forward_implications":["Any static vertical field in the rotor region can be expanded as $k_1 z + k_3 z^3$ or decomposed into symmetric and anti-symmetric parts, so the method transfers to other magnet geometries without numerical simulation.","The magnetic torque is exactly linear in $\\omega$, so the deceleration curve is exponential; this is a sharp prediction that high-precision measurements can confirm or reject.","The closed-form currents give a quantitative laboratory check for eddy-current demonstrations and a quick estimate of eddy-current losses in thin conducting disks.","The dominant identified error comes from a single omission, the horizontal field component, so including that component in the torque integral is a direct route toward closing the 8.61 percent gap."],"supporting_citations":[{"why":"Earlier rotating-disk eddy-current calculation that this paper generalizes from cylindrical poles to a cubic vertical-field profile.","marker":"[4]"},{"why":"Alternative rotating-frame derivation whose agreement with [4] anchors the existing eddy-current results this work extends to the lab frame.","marker":"[5]"},{"why":"Rigorous lab-frame treatment of a rotating sphere that yields a Legendre-series solution, the complexity of which motivates the closed-form approach here.","marker":"[6]"},{"why":"Maxwell equations, the starting point from which the paper derives the reduced Poisson equation.","marker":"[7]"},{"why":"Ohm's law in the form $\\mathbf{j}=\\sigma(\\mathbf{E}+\\mathbf{f})$, the constitutive relation linking current to electric field and motional force.","marker":"[8]"},{"why":"Magnetostatic formula for rectangular permanent magnets, used to obtain the field distribution and to fit the coefficients $k_1$ and $k_3$.","marker":"[9]"},{"why":"Provides the quasi-static condition $\\epsilon\\,\\partial\\mathbf{E}/\\partial t \\ll \\sigma\\mathbf{E}$ used to justify treating the current as nearly time-independent.","marker":"[12]"}],"fun_headline_variants":["Eddy current solved analytically, matches gyroscope to 8.61% error","Poisson equation yields closed-form eddy current, 8.61% error","Analytical eddy current solution predicts gyroscope spin-down to 8.61%","Closed-form eddy currents from Poisson equation, verified to 8.61%","Damping prediction hits 8.61% error with analytical eddy current"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that, in the rotor region, the vertical magnetic field depends only on the in-plane coordinate $z$ and is described by $B_z=k_1 z + k_3 z^3$, while the comparable horizontal field can be neglected; the authors estimate that neglect contributes about 10 percent of the torque, nearly the entire reported 8.61 percent discrepancy.","fun_headline_variants_meta":{"raw":{"variants":["Eddy current solved analytically, matches gyroscope to 8.61% error","Poisson equation yields closed-form eddy current, 8.61% error","Analytical eddy current solution predicts gyroscope spin-down to 8.61%","Closed-form eddy currents from Poisson equation, verified to 8.61%","Damping prediction hits 8.61% error with analytical eddy current"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000643,"raw_usage":{"total_tokens":2891,"prompt_tokens":815,"completion_tokens":2076,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":431,"completion_tokens_details":{"reasoning_tokens":1982}},"tokens_in":431,"tokens_out":2076,"duration_ms":15105,"temperature":1.0,"reasoning_tokens":1982,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:34:30.075461+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the horizontal field component in the rotor plane with a calibrated Hall probe and include its torque in the equation of motion; if the decay-time gap closes from 8.61 percent to near zero, the neglect is confirmed, whereas a sizeable remaining gap would implicate the cubic polynomial fit or the quasi-static assumption.","supporting_citations":[{"cited_title":"Gyroscope Teslameter","cited_arxiv_id":null,"evidence_quote":"Earlier rotating-disk eddy-current calculation that this paper generalizes from cylindrical poles to a cubic vertical-field profile."},{"cited_title":"Magnet / Gyroscope MYSTERY! Solve this unseen video","cited_arxiv_id":null,"evidence_quote":"Alternative rotating-frame derivation whose agreement with [4] anchors the existing eddy-current results this work extends to the lab frame."},{"cited_title":"Feynman, Robert B","cited_arxiv_id":null,"evidence_quote":"Rigorous lab-frame treatment of a rotating sphere that yields a Legendre-series solution, the complexity of which motivates the closed-form approach here."},{"cited_title":"On Eddy Currents in a Rotating Disk","cited_arxiv_id":null,"evidence_quote":"Maxwell equations, the starting point from which the paper derives the reduced Poisson equation."},{"cited_title":"Braking torque on rotating sheet in stationary magnetic field","cited_arxiv_id":null,"evidence_quote":"Ohm's law in the form $\\mathbf{j}=\\sigma(\\mathbf{E}+\\mathbf{f})$, the constitutive relation linking current to electric field and motional force."},{"cited_title":"Drag and lift forces between a rotating conductive sphere and a cylindrical magnet","cited_arxiv_id":null,"evidence_quote":"Magnetostatic formula for rectangular permanent magnets, used to obtain the field distribution and to fit the coefficients $k_1$ and $k_3$."},{"cited_title":"Three-Dimentional Pure Permanent Magnet Undulator Design Theory","cited_arxiv_id":null,"evidence_quote":"Provides the quasi-static condition $\\epsilon\\,\\partial\\mathbf{E}/\\partial t \\ll \\sigma\\mathbf{E}$ used to justify treating the current as nearly time-independent."}],"review_version":1}