The theory of the SQUID is violating mathematics since it assumes that you can have a differentiable scalar field, take its gradient because it is differentiable, and then take a circular integral within the vector-field obtained by taking the gradient of the scalar field and then end up with a non-zero value for a line intetegral in this field. This violates the most basic tenets of vector calculus. A model; that requires that one must violate mathematical rules is not a model at all.
Once again, everyone but Johan is wrong. But in this case it's not just all physicists that are wrong, it's every mathematician since Cauchy and Riemann.
I'm actually going to try to explain this to you, Johan (a ridiculous endeavor, clearly).
Consider the complex plane. Every point is labelled by a complex number z=x+iy. Consider the function f(z)=z. You can regard it as a complex scalar field - a particularly simple one, and one that is obviously differentiable (it's actually much more than differentiable, it's analytic).
There another way we can write f(z)=z: as [latex]$f(z)=z=r e^{i \phi}$[/latex], where [latex]$r=\sqrt{x^2 + y^2}$[/latex] and [latex]$\phi=\tan^{-1}{y/x}$[/latex] as usual. As you should already know, φ=0 on the positive real axis, φ=π/2 on the positive imaginary axis, φ=π on the negative real axis, and φ=-3π/2 on the negative imaginary axis. At any point except just below the positive real axis, φ increases linearly as you move counterclockwise around any arc of a circle centered at z=0.
So, consider the gradient of φ at a point. The gradient is a vector, so let's look at the gradient component perpendicular to the radial direction from the origin. That gradient component is very simple, it's just 1/r (the radial component is simple too, it's 0).
Now compute the line integral of the gradient of φ around any circle centered on the origin. Since the component of the gradient tangent to the circle is constant and equal to 1/r, the result is 2πr*1/r=2π. Shocking - the integral of the gradient of a "scalar field" isn't zero!
How did that happen? Very simple: φ has a branch-cut discontinuity. But φ is the phase of of a perfectly differentiable scalar field, and it's closely related to its logarithm.
Note the following:
0) All I've done is explain the residue theorem (a foundational result of 19th century mathematics) in a slightly unconventional way, since [latex]$2\pi = (1/i)\int_C z^{-1} dz = (1/i)\int_C \partial_z \ln z ~dz = \int_0^{2\pi} \partial_\phi \phi = 2\pi$[/latex].
1) This is mathematically identical to the analysis of SQUIDs you've just declared "mathematically impossible".
2) Your other statements about the physics of Josephson junctions are just as plainly wrong. For example, your assertion that superconductors cease to superconduct in the DC Josephson effect as soon as the voltage is non-zero is immediately falsified by just about every experiment done since 1962, as is your assertion that the AC J. effect doesn't exist. I linked to a paper that shows the latter, which you of course ignored.
3) This will go way over Johan's head, but for anyone else, note that if the function f(z)=z
n, the integral I did would give 2πn. In general the phase of a smooth function must wind an integer number of times n. The fact that n must be an integer (for well behaved functions f(z)) is essentially the origin of flux quantization and explains much of the behavior of SQUIDs (and Josephson junctions in general).