Calculation of equipotentials and flux lines in axially symmetrical permanent magnet assemblies by computer

Two equivalent theoretical models of permanent magnets are used to develop algorithms for numerically computing the magnetic scalar potential and the magnetic vector potential in the vicinity of an axially symmetric array of pole pieces and permanent magnets. A computer program based on these algorithms calculates equipotential surfaces and flux lines in and around the magnets and pole pieces. In deriving the algorithm for numerically calculating the vector potential a relationship between the magnetic scalar potential and the vector potential was found which enables the program to calculate the vector potential from the scalar potential distribution and thus generate equipotentials and flux lines with only one iterative calculation. An algorithm which calculates the scalar potential of a "floating" pole piece, that is, one on which the scalar potential has not been specified, is developed. The vector potential around the pole piece is determined from the scalar potential calculation, and this information is used to calculate the vector potential and the flux lines within the pole piece. The computer program calculates the coordinates of all points at which the equipotential lines and flux lines cross the Liebmann net. This information is fed to a cathode ray tube plotter which generates a field plot. To deal with systems in which macroscopic currents are present as well as permanent magnets, the iterative Liebmann net calculation of the vector potential is developed, and a method of applying Neumann boundary conditions to the vector potential at high-permeability surfaces is described.