Journal of Computational Physics

A numerical solution of a model for a superconductor field problem

View publication


A model of a magnetic field problem occurring in connection with Josephson junction devices is derived, and numerical solutions are obtained. The model is of mathematical interest, because the magnetic vector potential satisfies inhomogeneous Helmholtz equations in part of the region, i.e., the superconductors, and the Laplace equation elsewhere. Moreover, the inhomogeneities are the gauge constants for the potential, which are different for each superconductor, and their magnitudes are proportional to the currents flowing in the superconductors. These constants are directly related to the self and mutual inductances of the superconducting elements in the device. The numerical solution is obtained by the iterative use of a fast Poisson solver. Chebyshev acceleration is used to reduce the number of iterations required to obtain a solution. A typical problem involves solving 100,000 simultaneous equations, which the algorithm used with this model does in 20 iterations, requiring three minutes of CPU time on an IBM YM/370/168. Excellent agreement is obtained between calculated and observed values for the inductances. © 1979.