Comparison of the values found in different articles with that measured in this paper. The values and relative errors have been obtained directly from graphs, at 77 K. Available values for H (a,b) and H ‖ (c) in monocrystals are shown. H ‖ (c) is always greater than H (a,b)
It has long been known that a repulsive force arises between a magnetic field (generated, for instance, by a permanent magnet - PM) and a superconductor –Sc (Arkadiev, 1947). This force is due to the repulsion of the magnetic field away from the superconductor – the Meissner effect. Type I superconductors only can be in the Meissner state, which means that a magnetic field will be always expelled from the superconductor, independently of its poles orientation. Nevertheless, type II superconductors may be in two different states: first, provided the magnetic field is low enough, they are at a Meissner state similar to type I superconductors. In this Meissner state they absolutely expel the magnetic field and prevalent repulsive forces appear. Second, for magnetic fields larger than the so-called First Critical Field HC1, the magnetic flux penetrates the superconductor creating a magnetization which contributes to an attractive resulting force. This second state is known as mixed state.
In 1953 Simon first tried to make a superconducting bearing (Simon, 1953) using superconductors in the mixed state.The first engine using a superconducting bearing was made in 1958 (Buchhold, 1960). After the discovery of high critical temperature superconductors (Bednorz & Müller, 1986), the Meissner repulsive force has become a popular way of demonstrating superconducting properties (Early et al., 1988).For calculating forces between a magnet and a superconductor it is necessary to have models that describe both the flux penetration state and the Meissner state repulsion. The first one can be solved by using conventional methods to compute forces between magnetic elements and magnetized volumes. However, for the Meissner state the question has remained open until these last years.
Several models using the method of images to calculate superconducting repulsion forces (Lin, 2006; Yang & Zheng, 2007) have been proposed. However, this method of images is limited to a few geometrical configurations that can be solved exactly, and the physical interpretation of the method is under discussion (Giaro et al., 1990; Perez-Diaz & Garcia-Prada, 2007). Furthermore, some discrepancies within experiments still exist (Hull, 2000).
A general local model based on London’s and Maxwell’s equations has been developed to describe the mechanics of the superconductor-permanent magnet system (Perez-Diaz et al., 2008). Due to its differential form, this expression can be easily implemented in a finite elements analysis (FEA) and is consequently appliable to any shape of superconductor in pure Meissner state (Diez-Jimenez et al. 2010).
In this chapter, we present the demonstration of the model, the implementation into a finite elements program, the experimental verification and its limit of application. To finish, we show an example of usage.
2. Magneto-mechanics of a superconductor in Meissner state
A superconductor is in a pure Meissner state when it is exposed to an externally applied magnetic field, , lower than a certain value, HC1. HC1 is a characteristic of the material (Alario & Vincent, 1991), which depends on temperature. In this case, it is assumed that both and are equal to zero inside the superconductor. When a magnetic field is then applied –for example by moving a permanent magnet close to the superconductor - a surface current density is generated on the outermost surfaces. According to the London equation, this current is confined only to a depth of λ(T). Type II superconductors, such as the rare earth oxide high temperature superconductors, have the highest values for λ, reaching typical values of thousands of Amstrongs (Umezawa & Crabtree, 1998). Therefore, as this paper deals with macroscopic elements, it can be approximated that current density has an infinitely localized surface current
where is a surface current density tangent to the surface vector field and is a Dirac delta function on
In the static limit it can be assumed that. Therefore it may simply be written as:
Note that is continuous and fulfils:, provided the permanent magnet does not touch the superconductor surface. On the contrary, both and are discontinuous at the superconducting surface. In particular, bothand are discontinuous.
By using the divergence theorem (Jackson, 1975) on a small parallelepiped with volume V, a face just above the superconductor surface and another parallel face under it, it can be written that:
where S is the surface defining this parallelepiped andits normal vector. By using Maxwell law (3) it can be reduced to:
But, taking (1) into account, and considering under the superconducting surface it can be written that:
where is evaluated at z=0+ (limit above the superconductor surface).
As this result is independent of the small parallelepiped previously chosen, the integrands must equal:
Therefore, an expression for the superconducting current as a function of the applied magnetic field may be written:
All expressions shown use the MKS unit system.
Applying the divergence theorem clearly shows that the total charge is always conserved, for whichever surface shape the superconductor has, provided the source of the applied field is outside the superconductor:
which is consistent with the previous expression for the superconducting current.
The total field thus generated for a semi-infinite plane is equivalent to that generated by the method of images (Cansiz & Hull, 2005; Hellman et al. 1988; Hull & Cansiz, 1999; Moon, 1994). As the expression (11) has been derived using only local arguments, it may be used for any shape of superconductor. It does not depend on the curvature of the surface.
2.1. Force calculation
The external force (by unit surface) experienced by the superconductor can be calculated by using Lorentz force.
it can be written that:
This is a local and exact expression for the “pressure” or more precisely “stress” or force per unit surface on the superconductor, which depends only on the applied magnetic field. It is useful for any shape of superconductor. This differs fundamentally from the general expression used to calculate the magnetic stress between magnetic materials as given by Moon.
According to Newton’s law, the force exerted by the superconductor on the magnet is simply the opposite one. Therefore, for any shape of superconductor, the force exerted by the superconductor on the magnet can be written as:
where the integration extends over the whole surface of the superconductor.
2.2. Torque calculation
The torque suffered by the superconductor can easily be deduced as :
where is the position vector between the differential surface element and the center of mass of the superconductor bulk. Again, the integration extends over the whole surface of the superconductor.
In order to calculate the moment applied over the magnet, the must be the position vector between the differential surface element and the center of mass of the permanent magnet. As noted previously, the force exerted by the superconductor on the external magnetic field (in this case a single permanent magnet) is simply the opposite one, and the same applies for the torque.
3. Finite elements implementation
Due to this differential form equation (16) can be easily implemented in a finite element program. A FEM algorithm has been adapted for the commercial software ANSYS. The SURF154 element of ANSYS was used insofar as it has defined a set of useful attributes e.g. the surface normal direction. The algorithm is valid in the context of a common electromechanical simulation. The steps for the simulation were:
Select Element Type: SOLID98 (with a maximum of one degree of freedom MAG) and SURF154.
Create the different materials to be used. For the superconductor bulk, air properties were used.
Generate the geometries of the volumes for the electromagnetic system.
Assign materials’ properties to each volume, selecting air for the superconductor.
Mesh the whole system with the SOLID98 element (as fine as is considered adequate - discussed further below).
Mesh the superconductor surface with the SURF154 element.
Apply the electromechanical loads to the system.
Solve the electromagnetic equation system.
Once the system has been solved, the algorithm can be applied using a Command List. Fig. 1 shows a flow-diagram of the procedure.
This procedure has to be performed for each piece of superconductor in the system. Should there be more than one piece, a different internal SURF154 element must be created and accordingly, the number of SURF154 elements in the first step (ESEL) must be changed.
These steps provide the three components of the force vector. The torque applied on the superconductor can also be calculated, from which the torques values can then be derived.
3.1. Results provide by the post-processing
The algorithm has been tested using one of the most common experiments found in relevant literature: a permanent magnet oriented vertically over a superconductor pile in any arbitrary position.
Firstly, an electromagnetic system composed of a small magnet suspended over a superconducting cylinder was designed, as shown in Fig. 2.
The dimensions of the superconductor were: 20 mm diameter and 7.5 mm height, and the small magnet was characterized by a coercivity of 875 kA/m and a remanence of 1.18 T, with a 3.5 mm diameter and a 2 mm height. The magnet’s centroid was placed 10 mm over the surface of the superconductor. The entire system was placed in surroundings measuring 100x100x100 mm.
The results that can be obtained are the distribution of forces, torque and current densities per surface element. In fig 3, these distributions for an arbitrary position of the magnet over the superconductor are shown.
The same simulation was repeated several times with different meshes, increasing the number of elements for the whole simulation. Using the parameter α, the fineness of the mesh can be defined as the ratio between the maximum of the area of the elements and the total area of the superconductor multiplied by 100.
Different meshes along with their respective α parameters are shown in fig. 4. The number of SURF154 elements and the values of the solution are also displayed.
Fig. 5 shows the relative error of the calculations in relation to the analytical solution (with a magnetic moment of 0.016 A m2). The higher the number of surface elements, the smaller the relative error of the result. For example, where α is smaller than 0.1 %, the resulting relative error is less than 3 %.
In Fig. force versus z are shown for different values of α. The FEM results tend towards the analytical values as α decreases. It is noted that the magnetic dipole approximation made for the analytical calculation only remains valid where there is a large distance between the permanent magnet and the superconductor.
The convergence of the algorithm has been checked in relation to the finite elements’ size, and compared to analytical solutions for simple geometries. An α parameter has been proposed to assess the relative error in the results. The results showed good accuracy, whilst not requiring high specification computing technology.
4. Experimental verification
Different experiments were carried out in order to check the validity of the model. Some of them will be summarized in the following.
4.1. Force measurement
The following methodology was used to measure the forces: a cylindrical superconductor made of polycrystalline YBa2Cu3O7-x,manufactured by CAN superconductors (Kamenice 25168, Czech Republic) was immersed in a bath of liquid nitrogen N2 (77 K) at ambient pressure. The cylinder had a diameter of 45 mm and a height of 13 mm. It was fixed to a nitrogen vessel. The vessel, containing the superconductor, was placed on a lab jack stand to adjust the height. A small cylindrical permanent magnet was used, which had a coercivity of 875 kA/m, a remanence of 1.18 T, and had a diameter of 5 mm and a height of 5 mm. All experimental measurements followed the same coordinate system shown in Fig. 7. The origin of the coordinates was set at the center of the upper surface of the superconductor.
The permanent magnet was placed over the superconductor (Z coordinate),and fixed vertically to a PVC cantilever according to its magnetization direction (θ = 90º). The cantilever had 2 pairs of strain gauges to measure vertical forces at its extremes. This strain gauge configuration is not sensitive to the lateral and axial forces. The torques were neglected due to the size of the magnet. The PVC cantilever was joined to a 3D positioning table. The position of the magnet was then fixed in relation to the superconductor surface with a precision of 0.1 mm. The strain gauges were calibrated using a dynamometer and a set of 12 references forces. The calibration constant was established by least squares fitting in K = (3.87±0.14)×10-4 N/με, with a correlation coefficient of R2 = 0.997.
The measurement for every position was made in zero field cooling conditions (ZFC). The vertical forces were recorded where X = 0.0, 5.0, 10.0, 15.0, 17.5, 20.0, 22.5 and 25.0 ± 0.1 mm; at 3 different heights from the surface of the superconductor: 12.0, 10.0, and 8.0 ± 0.1 mm. Furthermore, measurements were taken in the center of the upper superconductor face, X = 0 mm from Z = 7.0 to 14.0 ± 0.1 mm. The Y position was always fixed at 0 mm.
These positions were chosen in order to avoid exceeding a limit of 3.5 mT of magnetic flux density at any point of the superconductor surface. Using this limit ensures the Meissner state is retained. Regardless of this, after every measurement the remanent magnetization of the superconductor bulk was checked and in most cases no measurable magnetization was found.
The FEA errors are estimated to be less than 3 %.
The figures show positive agreement between experimental and theoretical values. Only Fig. 9 shows an appreciable difference for values approaching X = 20 mm. However, it must be pointed that the radius of the superconductor is 22.5 mm. It is only in these surroundings that a very low remanent magnetization was recorded, which indicates a non complete Meissner state. This explains why some experimental values were lower than those of a complete Meissner state.
4.2. Equilibrium angle measurement
In addition to previous experiments, the mechanical behavior of a magnet which has the ability to tilt over the superconductor in the Meissner state was also studied in this paper. In the present experiment only one degree of freedom was permitted in the tilt angle of the magnet (θ coordinate). The equilibrium angle of the permanent magnet over the cylindrical superconductor was measured for different relative positions. The results can be used to understand not only how the permanent magnet is repelled, but also how it turns when it is released over a superconductor.
A cylindrical permanent magnet (made of NdFeB with a coercivity of 875 kA/m and a remanence of 1.29 T) was placed over the superconductor. Their dimensions were 6.3 mm in diameter and 25.4 mm in length and it had a magnetization direction parallel to its axis of revolution. A rigid plastic circular rod was fixed in the center of mass, perpendicular to the axis of revolution. This rod was used as the shaft in a plastic bearing, which was lubricated with oil. The whole bearing system was joined to a 3D displacement table. This arrangement ensured it was possible to control the position of the permanent magnet with an accuracy of 0.1 mm, and the only permitted degree of freedom was the rotation around the Y axis. Concentric to the bearing, a graduate goniometer measured the angle of rotation of the magnet. The whole experiment design is shown in Fig 12.
Fig. 13 shows the comparison between the equilibrium angles measured and those calculated by expression (19).
Again, there was a good agreement between the calculus made according to our model and the experiments. These experiments were carried out in Zero Field cooling condition (ZFC), and consequently there is no remanent magnetization.
5. Limits of application
The lower critical field,
Differences between the predicted Meissner forces and the experimentally measured ones indicate that a part of the sample is in the mixed-state. Establishing with precision the instant when the differences begin will permit us to determine the
The method also determines the zone at the sample where transition from Meissner to mixed state occurs.
For a position of the magnet with respect to the superconductor we define the Meissner Efficacy as
The measurement for every position was made in zero field cooling conditions (ZFC). The origin of coordinates was set at the center of the upper surface of the superconductor. The reference point of the magnet was placed in the center of the lower surface of the magnet. Therefore, the Z coordinate is the distance between the faces of the magnet and the superconductor. X is the distance of the center of the magnet to the axis of the superconductor cylinder (radial position). We have recorded the vertical forces for X = 0.0, 5.0, 10.0, 15.0, 17.5, 20.0, 22.5 and 25.0 ± 0.1 mm; at 4 different heights from the surface of the superconductor: 12.0, 10.0, 8.0 and 6.0 ± 0.1 mm.
Fig. 14 shows the Meissner Efficacy versus the maximum of the surface current density distribution
We observe that for low values of the maximum surface current density, the Meissner Efficacy is just 1.
From a certain value, the Meissner Efficacy decays linearly. From this data we can derive a weighted mean value of
In Table 1
Now, if we use a value of
|C. Bömer et al (2007)|
|Umewaza et al (2007)|
|Kaiser et al (1991)|
|Wu et al (1990)|
|Results for Hc1 (A/m)||2900 ± 250|
6000 ± 2300
|4500 ± 450||3580|
|3226 ± 176|
The uncertainty in the determination of
In fact, the values far from the Meissner state contribute to improve the accuracy of the
This results are in according to the border and thickness effects and border magnetization that have been already described by other authors in an uniform magnetic field ( Brandt, 2000; Morozov et al., 1996; Li et al., 2004; Schmidt et al., 1997):
6. Example of application - permanent magnet over a superconducting torus
We calculate the torque exerted between a superconducting torus and a permanent magnet by using this model. We find that there is a flip effect on the stablest direction of the magnet depending on its position. This could be easily used as a digital detector for proximity.
We consider a full superconducting torus and a cylindrical permanent NdFeB magnet over the superconductor axis (Z axis). In figure 15 we can observe the geometrical configuration of both components. Every calculation is referenced with respect of a Cartesian coordinate system placed in the center of mass of the torus which Z axis is coincident with the axis of the torus.
The superconducting torus has an internal radius RINT = 6 mm and a diameter of the section ØSECTION = 10 mm. The cylindrical permanent magnet has a length LPM = 5 mm and a diameter ØPM = 5 mm. When calculating the magnetic field generated by the magnet we define its magnetic properties as: Coercive magnetic field HCOERCIVITY = 875 kA/m and remanent magnetic flux density BREMANENT = 1.18 T. We assume that the direction of magnetization of the permanent magnet coincides with its axis of revolution.
The variables θ and z are the coordinates we modify in order to analyze the mechanical behavior of the magnet over the superconductor. z is the distance along the Z axis between the center of mass of the torus and the one of the cylindrical permanent magnet. θ is the angle between the axis of the magnet and the vertical Z axis.
The equilibrium angle (θeq) as a function of z can be determined as follows. For a certain z we calculate the Y component of the torque (My) exerted on the magnet by the superconductor as a function of θ and we find the equilibrium angle as the value for which My(θeq)=0. The sign of the slope dMy/dθ at that point determines the stability or instability of the equilibrium point.
In figure 16 the torque (My) exerted on the magnet by the superconductor as a function of θ is shown for z = 0, 3, 6, 9, 12 and 15 mm. The maximum values for the torque exerted to the permanent magnet appear at θ = 45º and θ = 135º for every z. The remarkable fact is that the sign suddenly changes when moving from z = 3 mm to z = 6 mm. The equilibrium points are always at θ = 0º and θ = 90º, but θ = 0º is a stable equilibrium point for z = 0 mm and z = 3 mm, while it is unstable for the rest of the positions. On the other hand θ = 90º is unstable for z = 0 mm and z = 3 mm, but it is stable for the rest of the positions. That means that if you approach a magnet along the Z axis and it is able to rotate, it will be perpendicular to the Z axis while it is at z ≥ 6 mm, but it will suddenly rotate to be parallel to the Z axis when you pass from z = 6 to z ≤ 3.
In figure 17 the variation of the torque at θ = 45º as a function of z. The torque changes its sign between z =3 mm and z =4 mm.
Finally, figure 18 shows the stable equilibrium angle as a function of z. It is evident that, at a certain position between z =+ 3 and z =+ 4 mm we found that the stable equilibrium angle switches from a vertical orientation of the magnet to an horizontal one describing the flip effect claimed in this work.
Therefore, it can be concluded that if you approach a magnet along the Z axis and it is able to rotate, it will be perpendicular to the Z axis while it is at a certain distance (z ≥ 4 mm in our example) and it will change to be parallel to the Z axis for closer positions (z ≤ 3 mm in our example). As the equilibrium angle does not depend on the magnetic moment, the magnet can be much smaller. As a flip in the orientation of a permanent magnet can be easily instrumented, this effect can be easily used as a binary detector for proximity.
Magnet-superconductor forces both in Meissner and mixed states can be calculated with the accuracy required to engineer useful levitating devices.
The implementation of a local differential expression in a finite elements program opens new perspectives to the use of magnet-superconductor devices for engineering. This can be used to calculate forces whatever the size, shape and geometry of the system, for both permanent magnets and electromagnets.
Accuracy and convergence, in addition to the experimental verification for different cases have been tested. There is a good agreement between experimental results and calculation, even with very low-cost computing resources involved.
Moreover, the expression can be used to determine the point when the mixed state arises in a superconductor piece.