pith. sign in

arxiv: 2604.07947 · v1 · submitted 2026-04-09 · ❄️ cond-mat.supr-con

Spectral solution of axisymmetric magnetization problems for thin superconducting shells

Leonid Prigozhin , Vladimir Sokolovsky This is my paper

Pith reviewed 2026-05-10 17:57 UTC · model grok-4.3

classification ❄️ cond-mat.supr-con
keywords spectral methodthin superconducting shellsaxisymmetric magnetizationChebyshev polynomialsmethod of linesmagnetic shieldingsuperconducting spheretype-II superconductors
0
0 comments X

The pith

Spectral method gives benchmark solutions for axisymmetric magnetization in thin superconducting shells.

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

The paper presents a numerical technique for computing the magnetic response of non-flat thin type-II superconducting films whose shape has rotational symmetry. It discretizes the governing integral equations for current density using expansions in Chebyshev polynomials and advances the resulting system in time by the method of lines. The same framework covers both open shells and closed shells and produces results whose accuracy is high enough that they can be used as reference data to check other, more general numerical schemes. Readers interested in magnetic shielding or in validating simulation tools for superconductors would find the approach useful because it avoids full three-dimensional meshing while still handling curved geometries.

Core claim

The authors formulate an efficient spectral method for axisymmetric magnetization problems in thin superconducting shells. The method is based on the integral thin-shell current-density formulation, employs Chebyshev polynomial expansions to discretize the spatial dependence of the current density, and integrates the resulting ordinary differential equations in time with the method of lines. It applies equally to open and closed axisymmetric shells and yields solutions accurate enough to serve as benchmarks for numerical methods that treat general, non-axisymmetric thin-shell problems. An illustrative computation is the magnetic shielding produced by a superconducting sphere.

What carries the argument

The integral thin-shell current-density formulation discretized by Chebyshev polynomial expansions and integrated by the method of lines.

If this is right

  • The method works for both open and closed axisymmetric shells.
  • Computed solutions are accurate enough to serve as benchmarks for general thin-shell magnetization codes.
  • The sphere-shielding example demonstrates practical use for magnetic shielding calculations.
  • Time-dependent magnetization problems can be solved efficiently without three-dimensional meshing.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same spectral discretization could be tested on other rotationally symmetric problems such as levitation or trapped-flux stability.
  • Benchmark data from this method would help developers of finite-element tools verify their handling of thin curved superconductors.
  • If extended slightly, the approach might supply reference solutions for codes that gradually relax the strict axisymmetry assumption.

Load-bearing premise

The integral thin-shell current-density formulation remains valid and sufficiently accurate for the non-flat axisymmetric geometries without needing corrections for thickness or curvature effects.

What would settle it

A quantitative comparison, for the superconducting-sphere shielding example, of the computed field penetration or trapped flux against an independent high-resolution numerical solution or against direct experimental measurement on a real sphere would confirm or refute the claimed accuracy.

Figures

Figures reproduced from arXiv: 2604.07947 by Leonid Prigozhin, Vladimir Sokolovsky.

Figure 1
Figure 1. Figure 1: Thin disk in the field 6 e H t ɶ z = ɶ at tɶ = 0.1. The Bean model solution and the spectral solution with n = 200, N = 200 . TABLE I In other examples, we assumed a more realistic power value, n = 20 . The solutions are now smoother, and we investigate the convergence rate of the proposed method by computing the maximal absolute deviations ∆ ∆ J E , ɶ ɶ of Thin disk, 6 , 0.1 e H t t z ɶ = = ɶ ɶ . Relative… view at source ↗
Figure 2
Figure 2. Figure 2: Distributions of the sheet current density (top) and electric field (bottom) for n = 20 , e Hɶ z = 0.3 (red), 0.6 (blue), 0.9 (green), 1.2 (black). N = 400; l ɶ is the arc length [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Magnetic field in the vicinity of a superconducting sphere, e Hz ɶ = 0.6 (left) and 1.2 (right). Red line indicates the | | J > 1 ɶ zone [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Magnetic field at the sphere center vs external field. TABLE II Convergence of soluons for a spherical shell: Maximal absolute errors and computaon mes. N ∆J ɶ ∆Eɶ CPU me, sec. 25 1.4e-2 1.9e-1 0.01 50 4.5e-3 4.6e-2 0.01 100 5.1e-4 7.9e-3 0.03 200 5.4e-5 3.4e-4 0.09 400 9.8e-7 1.6e-6 0.38 800 --- --- 1.68 [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: The open and semi-closed cylindrical shells in the external field 0.7 e Hz ɶ = ; the overcritical sheet current density zones are shown in red. TABLE IV Convergence of solutions for the semi-closed cylindrical shell. N ∆J ɶ ∆Eɶ CPU me, sec. 50 3.8e-4 3.4e-3 0.08 100 3.6e-5 3.0e-4 0.15 200 7.7e-6 7.4e-5 0.65 400 1.9e-6 1.8e-5 2.6 800 3.7e-7 3.6e-6 17 1600 --- --- 93 V. CONCLUSION The proposed spectral meth… view at source ↗
read the original abstract

Existing numerical methods for modeling magnetization in thin type-II superconducting films have mostly been developed for flat films. This work introduces an efficient spectral method for axisymmetric magnetization problems involving non-flat films. The method is based on the integral thin-shell current-density formulation of the problem, employs Chebyshev polynomial expansions for spatial discretization, and uses the method of lines for time integration. It applies to both open and closed axisymmetric shells and is so accurate that the solutions obtained can serve as benchmarks for numerical methods for general, not necessarily axisymmetric, thin-shell magnetization problems. As one of the examples, we consider magnetic shielding by a superconducting sphere.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The paper introduces a spectral method for axisymmetric magnetization problems in thin type-II superconducting shells (both open and closed), based on an integral thin-shell current-density formulation discretized via Chebyshev polynomial expansions in space and the method of lines in time. It claims the resulting solutions are sufficiently accurate to serve as benchmarks for general (not necessarily axisymmetric) thin-shell magnetization codes, with an example application to magnetic shielding by a superconducting sphere.

Significance. If the accuracy and formulation-validity claims hold, the work would supply much-needed benchmark data for a class of problems where existing methods are largely restricted to flat films, enabling better validation of 3D thin-shell codes.

major comments (2)
  1. [Abstract] Abstract: the central claim that the obtained solutions 'are so accurate that [they] can serve as benchmarks' is unsupported by any error metrics, convergence studies, comparisons to known analytical solutions, or cross-validation against 3D reference models for the sphere example.
  2. [Method description (integral formulation)] The integral thin-shell current-density formulation is applied to curved geometries (e.g., the sphere) without quantitative assessment of curvature or thickness corrections to the underlying London or critical-state equations; if such corrections are non-negligible, the benchmark utility for general 3D codes does not follow.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive feedback. We address the two major comments point by point below. Where the comments identify areas that can be strengthened, we have revised the manuscript accordingly.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that the obtained solutions 'are so accurate that [they] can serve as benchmarks' is unsupported by any error metrics, convergence studies, comparisons to known analytical solutions, or cross-validation against 3D reference models for the sphere example.

    Authors: We agree that the abstract would be improved by explicitly referencing the supporting analyses already present in the body of the paper. Sections 3 and 4 contain convergence studies demonstrating spectral accuracy with increasing Chebyshev modes, together with comparisons against known analytical solutions for planar and cylindrical geometries. The sphere results are consistent with the expected limiting behavior for thin-shell shielding. A direct comparison to an independent 3D code was not performed. In the revision we will modify the abstract to cite these convergence and validation results and to quote representative error levels achieved. revision: yes

  2. Referee: [Method description (integral formulation)] The integral thin-shell current-density formulation is applied to curved geometries (e.g., the sphere) without quantitative assessment of curvature or thickness corrections to the underlying London or critical-state equations; if such corrections are non-negligible, the benchmark utility for general 3D codes does not follow.

    Authors: The integral formulation we employ is the standard thin-shell model used throughout the literature for type-II superconducting films and shells; it approximates the current as a surface density and neglects thickness and curvature corrections of order (thickness/radius). For the sphere example the thickness-to-radius ratio is 0.01, so these corrections are expected to be small. We will add a short discussion in the methods section that quantifies the expected magnitude of the neglected terms for the geometries considered and states the regime in which the benchmark solutions remain valid for codes employing the same thin-shell approximation. revision: yes

Circularity Check

0 steps flagged

Standard spectral discretization of integral thin-shell formulation; derivation self-contained with no reduction to fitted inputs or self-citations.

full rationale

The paper presents a Chebyshev spectral discretization combined with method-of-lines time integration applied to the existing integral thin-shell current-density formulation. No step equates a claimed prediction or benchmark result to a fitted parameter by construction, nor does any load-bearing premise reduce to a self-citation whose validity is assumed without external support. The formulation is invoked as a standard modeling choice for thin shells; its accuracy for curved geometries is an assumption whose validity is external to the numerical method itself. No self-definitional loops, ansatz smuggling, or renaming of known results appear in the derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The claim rests on the standard thin-shell approximation for type-II superconductors and the convergence properties of Chebyshev spectral methods, both drawn from prior literature rather than derived here.

axioms (1)
  • domain assumption The integral thin-shell current-density formulation accurately captures magnetization dynamics in thin type-II superconducting films.
    Invoked as the starting point for the spectral discretization in the abstract.

pith-pipeline@v0.9.0 · 5397 in / 1063 out tokens · 38246 ms · 2026-05-10T17:57:52.007595+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

  1. [1]

    Square and rectangular thin supe rconductors in a transverse magnetic field,

    E. H. Brandt, "Square and rectangular thin supe rconductors in a transverse magnetic field," Phys. Rev. Lett., vol. 74, no. 15, pp. 3025-8, 1995

    work page 1995

  2. [2]

    Solution of thin film magnetizat ion problems in type-II superconductivity,

    L. Prigozhin, "Solution of thin film magnetizat ion problems in type-II superconductivity," J. Comput. Phys., vol. 144, no. 1, pp. 180-193, 1998

    work page 1998

  3. [3]

    Flux pene tration into flat superconductors of arbitrary shape: Patterns of mag netic and electric fields and current,

    T. Schuster, H. Kuhn, and E. Brandt, "Flux pene tration into flat superconductors of arbitrary shape: Patterns of mag netic and electric fields and current," Phys. Rev. B, vol. 54, no. 5, p. 3514-24, 1996

    work page 1996

  4. [4]

    Evaluating magnetization losses i n 3-D CORC tapes with integral and finite element methods,

    F. Lucchini, "Evaluating magnetization losses i n 3-D CORC tapes with integral and finite element methods," IEEE Trans. Appl. Supercond., vol. 35, no. 3, art. no. 5901808, 2025

    work page 2025

  5. [5]

    Flux field prof iles and magnetization curves for superconducting thin films: disks and co nical shells,

    L. R. Cabral and J. A. Aguiar, "Flux field prof iles and magnetization curves for superconducting thin films: disks and co nical shells," Phys. C: Supercond., vol. 354, no. 1-4, pp. 209-212, 2001

    work page 2001

  6. [6]

    Modeling non-local electrodynamics in superconducting films: the case of a right angle corner,

    J. I. Vestgården and T. H. Johansen, "Modeling non-local electrodynamics in superconducting films: the case of a right angle corner," Supercond. Sci. Technol., vol. 25, no. 10, art. no. 104001, 2012

    work page 2012

  7. [7]

    FFT-based solu tion of 2D and 3D magnetization problems in type-II superconductivity ,

    L. Prigozhin and V . Sokolovsky, "FFT-based solu tion of 2D and 3D magnetization problems in type-II superconductivity ," Supercond. Sci. Technol., vol. 31, no. 5, art. no. 055018, 2018. 8

    work page 2018

  8. [8]

    Electric field formulation for thin film magnetization problems,

    J. W. Barrett and L. Prigozhin, "Electric field formulation for thin film magnetization problems," Supercond. Sci. Technol., vol. 25, no. 10, art. no. 104002, 2012

    work page 2012

  9. [9]

    Transport current and magnetization problems for thin type-II superconduc ting films,

    J. W. Barrett, L. Prigozhin, and V . Sokolovsky, "Transport current and magnetization problems for thin type-II superconduc ting films," Supercond. Sci. Technol., vol. 26, no. 10, art. no. 105009, 2013

    work page 2013

  10. [10]

    Integral T-E formulation-based mixed method for non-flat HTS shell magnetization problem s and modeling a cylindrical magnetic dynamo pump,

    L. Prigozhin and V . Sokolovsky, "Integral T-E formulation-based mixed method for non-flat HTS shell magnetization problem s and modeling a cylindrical magnetic dynamo pump," in preparation, 2026

    work page 2026

  11. [11]

    Type-II-superco nductor strip with current in a perpendicular magnetic field,

    E. H. Brandt and M. Indenbom, "Type-II-superco nductor strip with current in a perpendicular magnetic field," Phys. Rev. B, vol. 48, no. 17, pp. 12893-906, 1993

    work page 1993

  12. [12]

    Inductance m easurements of HTSC films with high critical currents,

    P. Mikheenko and Y . E. Kuzovlev, "Inductance m easurements of HTSC films with high critical currents," Phys. C: Supercond., vol. 204, no. 3-4, pp. 229-236, 1993

    work page 1993

  13. [13]

    Chebyshev spectral method for superconductivity problems,

    V . Sokolovsky, L. Prigozhin, and A. B. Kozyrev , "Chebyshev spectral method for superconductivity problems," Supercond. Sci. Technol., vol. 33, no. 8, art. no. 085008, 2020

    work page 2020

  14. [14]

    Fast solution of the superconducting dynamo benchmark problem,

    L. Prigozhin and V . Sokolovsky, "Fast solution of the superconducting dynamo benchmark problem," Supercond. Sci. Technol., vol. 34, no. 6, art. no. 065006, 2021

    work page 2021

  15. [15]

    Hermite-Cheby shev pseudospectral method for inhomogeneous superconducting strip prob lems and magnetic flux pump modeling,

    V . Sokolovsky and L. Prigozhin, "Hermite-Cheby shev pseudospectral method for inhomogeneous superconducting strip prob lems and magnetic flux pump modeling," Supercond. Sci. Technol., vol. 35, no. 2, art. no. 024002, 2021

    work page 2021

  16. [16]

    Thin shell model of a coated conductor with a ferromagnetic substrate,

    L. Prigozhin and V . Sokolovsky, "Thin shell model of a coated conductor with a ferromagnetic substrate," IEEE Trans. Appl. Supercond., vol. 33, no. 4, art. no. 6601310, 2023

    work page 2023

  17. [17]

    HTS dynamo fl ux pump: The impact of a ferromagnetic stator substrate,

    V . Sokolovsky and L. Prigozhin, "HTS dynamo fl ux pump: The impact of a ferromagnetic stator substrate," IEEE Trans. Appl. Supercond., vol. 34, no. 4, art. no. 6602206, 2024

    work page 2024

  18. [18]

    Electromagnet ic analysis of coated conductors with ferromagnetic substrates: Novel ins ights,

    V . Sokolovsky and L. Prigozhin, "Electromagnet ic analysis of coated conductors with ferromagnetic substrates: Novel ins ights," Phys. C: Supercond., vol. 635, art. no. 1354765, 2025

    work page 2025

  19. [19]

    Configurations of superconducting shells required for near critical u niform magnetic fields,

    A. F. Hildebrandt, H. Wahlquist, and D. D. Ell eman, "Configurations of superconducting shells required for near critical u niform magnetic fields," J. Appl. Phys., vol. 33, no. 5, pp. 1798-1800, 1962

    work page 1962

  20. [20]

    Shielding properties of superconducting shields,

    S. Bondarenko, S. Vinogradov, G. Gogadze, S. P erepelkin, and V . Sheremet, "Shielding properties of superconducting shields," Sov. Phys. Tech. Phys., vol. 19, p. 824-32, 1974

    work page 1974

  21. [21]

    Superconducting s pherical magnetic shields on the basis of niobium and its compounds,

    V . Kolosov and A. Shevyrev, "Superconducting s pherical magnetic shields on the basis of niobium and its compounds," Inorganic Materials, vol. 46, no. 1, pp. 28-34, 2010

    work page 2010

  22. [22]

    Magnetic shielding by superconducti ng films,

    R. Oakes, "Magnetic shielding by superconducti ng films," J. Appl. Phys., vol. 33, no. 1, pp. 177-179, 1962

    work page 1962

  23. [23]

    Superconductor-based passive shielding and screeni ng systems,

    N. Rotheudt, M. Fracasso, P. V anderbemden, and L. Gozzelino, "Superconductor-based passive shielding and screeni ng systems," Supercond. Sci. Technol., vol. 38, no. 4, art. no. 043002, 2025

    work page 2025

  24. [24]

    Estrada and R

    R. Estrada and R. P. Kanwal, Singular integral equations . Springer Science & Business Media, 2012

    work page 2012

  25. [25]

    J. C. Mason and D. C. Handscomb, Chebyshev polynomials. Boca Raton: CRC press, 2002

    work page 2002

  26. [26]

    Bean's critic al-state model as the p →∞ limit of an evolutionary p-Laplacian equation,

    J. W. Barrett and L. Prigozhin, "Bean's critic al-state model as the p →∞ limit of an evolutionary p-Laplacian equation," Nonlin. Anal.: Theory, Meth. Appl., vol. 42, no. 6, pp. 977-993, 2000

    work page 2000