flow rate measurements of liquid metals are required in various technological processes ranging from the cooling of nuclear reactors to the dosing and casting of molten metals.@xcite electromagnetic flowmeters are essential in the diagnostics and automatic control of such processes . a variety of electromagnetic flowmeters have been developed starting from the late 1940s and described by shercliff.@xcite the standard approach is to determine the flow rate by measuring the potential difference induced between a pair of electrodes by a flow of conducting liquid in the magnetic field.@xcite this approach is now well developed and works reliably for common liquids like water,@xcite but not so for liquid metals .
major problem in molten metals , especially at elevated temperatures , is the electrode corrosion and other interfacial effects , which can cause a spurious potential difference between the electrodes .
the electrode problem is avoided by contactless eddy - current flowmeters , which determine the flow rate by sensing the flow - induced perturbation in an applied magnetic field.@xcite the main problem with this type of flowmeters is the weak field perturbation which may be caused not only by the flow .
we showed recently that the flow - induced phase shift of ac magnetic field is more reliable for flow rate measurements than the amplitude perturbation.@xcite    another contactless techniques for flow rate measurements in liquid metals is the so - called magnetic flywheel invented by shercliff,@xcite who prescribes _ a `` plurality '' _ of permanent magnets distributed equidistantly along the circumference of a disk , which is mounted on an axle and placed close to a tube carrying the liquid metal flow .
the eddy currents induced by the flow across the magnetic field interact with the magnets by entraining them , which makes the disk rotate with a rate proportional to that of the flow .
this type of flowmeter , described also in the textbook by shercliff@xcite and extensively used by bucenieks,@xcite was recently successfully reembodied under the name of the lorentz force velocimetry ( lfv).@xcite    recently , we suggested an alternative and much more compact design of such a flowmeter , which conversely to shercliff s flywheel uses just a single magnet mounted on the axle it can freely rotate around and magnetized perpendicularly to it.@xcite we also introduced a basic mathematical model and presented first experimental implementation of this type of flowmeter.@xcite when such a magnet is placed properly at a tube with the liquid metal flow , it starts to revolve similarly to shercliff s flywheel .
but in contrast to the latter , which is driven by the electromagnetic force acting on separate magnets , the single magnet is set into rotation only by the torque . this driving torque is due to the eddy currents induced by the flow across the magnetic field . as the magnet starts to rotate ,
additional eddy currents are induced , which brake the rotation .
an equilibrium rotation rate is attained when the braking torque balances the driving one , and this rate depends only on the flow velocity and the flowmeter arrangement , whereas it is independent of the electromagnetic torque itself .
thus , the equilibrium rotation rate is affected neither by the magnet strength nor by the electrical conductivity of the liquid metal provided that the friction on the magnet is negligible .
this a major advantage of the single - magnet rotary flowmeter over the lfv approach , which relies on direct force measurements.@xcite    in this paper , we present an extended theory of the single - magnet rotary flowmeter and compare it with experimental results .
two limiting cases of long and short magnets are analyzed using linear - dipole and single - dipole approximations .
we obtain simple analytic solutions for the force and torque on slowly moving and rotating magnets due to eddy currents in the layers of infinite horizontal extent and arbitrary depth .
this allows us to find the equilibrium rotation rate of the magnet at which the torques due to the translation and rotation balance each other .
we also consider an active approach , where the force on the magnet is used to control its rotation rate so that the resulting force vanishes .
this rotation rate , similarly to the equilibrium one , is proportional to the layer velocity and independent of its conductivity and the magnet strength .
the torque on a magnetic dipole rotating about an axis normal to a thin sheet has been calculated by smythe using an original receding image method.@xcite reitz uses this method to calculate the lift and drag forces on the coils of various geometries moving with constant velocity above a conducting thin plate.@xcite the lift force , which at high speeds approaches the force between the coil and its image located directly below it , varies as the velocity squared in the low - speed limit considered in this paper .
the drag force is found to vary inversely and directly with the velocity at low and high speeds , respectively .
palmer finds analytical expressions for the eddy current forces on a circular current loop moving with a constant velocity parallel to a thin conducting sheet.@xcite the force on a rectangular coil moving above a conducting slab has been calculated numerically by reitz and davis using the fourier transform method.@xcite the same problem for the magnetic dipole of arbitrary orientation placed next to a thin slowly moving slab is addressed by kirpo _
@xcite the force and torque on a transversely oriented dipole above a slowly moving plane layer of arbitrary thickness has been found analytically in the context of the lfv.@xcite fast computation of forces on moving magnets are of interested also for the eddy current force testing techniques.@xcite    this paper is organized as follows .
the following section presents two simple mathematical models of the single - magnet rotary flowmeter , which are used to calculate analytically the force and torque on the magnets moving and rotating slowly above a layer of infinite lateral extent .
the limits of long and short magnets , which are approximated by linear and point dipoles , are considered in secs .
[ sub:2d ] and [ sub:3d ] , respectively .
section [ sec : exp ] presents the flowmeter implementation details and test results .
the paper is concluded by a summary in sec .
[ sec : con ] .
schematic view of the single - magnet rotary flowmeter with the @xmath0-axis directed out of the plane of figure . ]    consider a horizontally unbounded planar layer of electrical conductivity @xmath1 occupying the lower half - space and moving as a solid body with a constant velocity @xmath2 parallel to a permanent magnet placed at a distance @xmath3 above its surface and rotating with a constant angular velocity @xmath4 around an axis parallel to the surface and perpendicular to @xmath5 velocities are assumed sufficiently low for the magnetic field of induced currents to be negligible compared to the field of the magnet .
the origin of cartesian coordinates is set at the center of the magnet with the @xmath6 @xmath0 and @xmath7 axis directed along @xmath8 @xmath9 and downward normally to the surface , respectively , as shown in fig .
[ cap : sketch ] . in the following
, two limiting cases will be considered in which the magnet will be assumed either much longer or much shorter than @xmath10      we start with a long cylinder magnetized perpendicularly to its axis about which it can freely rotate.@xcite in this case , the magnetic field is approximated by that of a two - dimensional ( linear ) dipole with the vector potential @xmath11 which has only the @xmath0-component @xmath12 where @xmath13 is the vacuum permeability and @xmath14 is the linear dipole moment , which is perpendicular to @xmath15 and directed at the angle @xmath16 from the positive @xmath17 axis ; @xmath18 is the radius vector from the magnet axis , and @xmath19 and @xmath20 are the cylindrical radius and the polar angle in the cylindrical coordinates around the @xmath0-axis .
the magnetic field of linear dipole is @xmath21 eddy currents are induced by two effects : the translation of the layer and the temporal variation of the magnetic field due to its rotation .
the latter vanishes in the co - rotating frame of reference , where the magnetic field is stationary , while the layer appears to move with the resulting velocity @xmath22 which contains also an apparent rotational motion of the layer @xmath23 opposite to that of the magnet .
the density of eddy currents is given by ohm s law for a moving medium @xmath24 where @xmath25 is the electric potential .
in this case , no electric potential is induced because the e.m.f . , @xmath26 is both solenoidal and tangential to the surface .
if the induced magnetic field is negligible as originally assumed , eddy currents can be represented as the superposition @xmath27 where @xmath28 and @xmath29 are the currents induced by the translation and rotation , respectively .
it is important to note that in the approximation under consideration with a fixed magnetic field distribution , eddy currents are determined only by the instantaneous velocities . besides that eddy currents are coplanar to the surface and , thus , mutually independent over the depth of the layer
consequently , a layer of finite thickness may be represented as a semi - infinite one with zero velocity at @xmath30 where @xmath31 is the distance of the lower boundary of finite - thickness layer from the magnet .
this , in turn , is equivalent to the superposition of two semi - infinite layers with the second layer at @xmath32 moving oppositely to the first one at @xmath33 so that the resulting velocity vanishes at @xmath34 in the following , this approach allows us to construct the solution for a finite - thickness layer by taking the difference of two half - space solutions , which is subsequently denoted by @xmath35_{h_{1}}^{h_{2}},$ ] where @xmath36 stands either for the force or the torque due to the eddy currents in a half - space .
moreover , by the same arguments , this approach can easily be extended to @xmath7-dependent velocity distributions @xmath37 , for which general solution can be constructed as a superposition of solutions for thin layers moving with various velocities given by @xmath38 where @xmath39_{h_{1}}^{h_{2}}/\left[z\right]_{h_{1}}^{h_{2}}$ ] for @xmath40 and @xmath41    the linear force density experienced by an infinitely long magnet due to the layer translation , which according to the momentum conservation law is opposite to that acting upon the layer itself , can be written as @xmath42 where the integral is taken over the @xmath43-cross - section of the layer .
the @xmath0-component of force is absent due to the @xmath44 reflection symmetry . in the low - speed limit under consideration , when force varies linearly with the velocity , there is also no @xmath7-component of force .
this is the case because according to the linearity @xmath45 when @xmath46 while the latter transformation is equivalent to the rotation of the coordinate system by @xmath47 around the @xmath7-axis , which leaves the @xmath7-components invariant . in polar
coordinates with the surface defined by @xmath48 we obtain @xmath49 the linear torque density , which because of the aforementioned symmetries has only the @xmath0-component , can be found as @xmath50 where @xmath51 the linear force density due to the rotation , defined by @xmath52 is found as @xmath53 the linear torque density due to the rotation , as that due to the translation above , has only the @xmath0-component @xmath54 which is not defined for a semi - infinite layer because the last integral diverges .
nevertheless , expression ( [ eq : m1-gen ] ) can still be evaluated for the layer of finite depth by substituting the infinite limit in the last integral by @xmath55 which results in @xmath56_{h_{2}}^{h_{1}}.\label{eq : m1-h2}\ ] ]    the solution above becomes unbounded as @xmath57 and , thus , inapplicable to thick layers .
this implies that for a half - space layer the induced magnetic field can not be neglected however slow the rotation of magnet .
the induced magnetic field is related to eddy currents by ampere s law @xmath58 which combined with expressions ( [ eq : b2d ] ) and ( [ eq : ohm ] ) leads to @xmath59 as suggested by the external magnetic field ( [ eq : a2d ] ) , we search for the vector potential in the complex form @xmath60.\label{eq : a - cmp}\ ] ] then eq .
( [ eq : a ] ) for the complex amplitude @xmath61 takes the form @xmath62 where the prime stands for the derivative with respect @xmath63 the general solution of eq .
( [ eq : a - rot ] ) can be written as @xmath64 where @xmath65 is an unknown constant to be determined by matching the induced and externally imposed magnetic fields ; @xmath66 is the modified bessel function of the second kind,@xcite and @xmath67 is the skin depth due to the rotation . at distances @xmath68 smaller than the skin depth @xmath69 @xmath70
is expected to approach the imposed field @xmath71 which yields @xmath72 ) and @xmath73 substituting expressions ( [ eq : a - cmp ] ) and ( [ eq : a - hat ] ) into integral ( [ eq : m1-gen ] ) , the torque on the magnet can be represented as @xmath74 where @xmath75r\,\d r$ ] and the asterisk denotes the complex conjugate . using eq .
( [ eq : a - rot ] ) , after some algebra we obtain @xmath76=\left(\frac{\mu_{0}}{2\pi}\right)^{2}\re\left[k_{0}((1+i)r / d)k_{1}((1-i)r / d)(1+i)r / d\right].\ ] ] for the low rotation rates that satisfy @xmath77 , we obtain @xmath78=\left(\frac{\mu_{0}}{2\pi}\right)^{2}(\ln(2d / r)-\gamma),\ ] ] where @xmath79 is euler s constant . finally , taking into account that @xmath80 and @xmath81 we obtain@xmath82 where @xmath83 is the effective skin depth , which owing to the low velocities under consideration is supposed to be large relative to @xmath84 .
now we can use the results above to find the magnet rotation rate depending on the velocity of layer .
first , if the magnet rotates steadily without a significant friction , the torques due to the translation and rotation are at equilibrium : @xmath85_{h_{2}}^{h_{1}}=0,$ ] which yields @xmath86 note that this equilibrium rotation rate @xmath87 depends neither on the magnet strength nor on the layer conductivity unless @xmath88 in the latter case , the skin effect becomes important and , thus , @xmath31 has to be substituted by @xmath89 in the expression above .
the velocity of layer can be determined also in another way by measuring the force on the magnet , as in the lfv , to control the magnet rotation rate so that the resulting force vanishes , @xmath90_{h_{2}}^{h_{1}}=0.$ ] this results in@xmath91 which again depends linearly on the layer velocity , but does not depend on its conductivity or the magnet strength .      in the other limiting case of a short magnet ,
when the distance to the surface is large or at least comparable to the size of magnet , the latter can be considered as a dipole with the scalar magnetic potential @xmath92 where @xmath93 is the dipole moment and @xmath94 is the fundamental solution of laplace s equation , which satisfies @xmath95 with the dirac delta function on the r.h.s . in the following , we use simplified notation @xmath96 where @xmath97 is the unit vector .
then the dipole magnetic field is given by @xmath98 where the last term is added to ensure the solenoidality of @xmath99 also at @xmath100 the associated dipole current distribution is @xmath101 which easily leads to the classical expressions for the force and torque used later on . in the following
, we will be using also the spherical and cylindrical coordinates associated with the cartesian ones in the usual way . as in the previous section ,
we change to the co - rotating frame of reference and consider the eddy currents due to translation and rotation separately .      for the translation with
@xmath8 the charge conservation @xmath102 applied to eq .
( [ eq : ohm ] ) results in laplace s equation for @xmath103@xmath104 at the surface @xmath105 the normal component of electric current vanishes : @xmath106 in order to find the induced electric potential , firstly , it is important to notice that @xmath99 being a free - space magnetic field satisfies the laplace equation itself .
consequently , the cartesian components of @xmath99 satisfy this equation , too . secondly ,
if @xmath107 in bc ( [ bc : phi0 ] ) satisfies laplace s equation , then @xmath108 satisfies not only bc ( [ bc : phi0 ] ) but also eq .
( [ eq : phi0 ] ) because the integration along a straight line , similarly to the differentiation , are interchangeable with the laplacian . by the same argument , we can interchange integration and differentiation in expression ( [ sol : phi0-gen ] ) , which yields @xmath109 where @xmath110 and @xmath111 is the zeroth degree associated legendre function of the second kind;@xcite @xmath112 is the spherical polar angle from the positive @xmath7-axis ( see fig .
[ cap : sketch ] ) and @xmath68 is the corresponding cylindrical radius .
the `` constant '' of integration @xmath113 which similarly to the first term @xmath114 is supposed to be axisymmetric and also to satisfy the laplace equation , is chosen to regularize @xmath115 at @xmath116 by removing the logarithmic singularity @xmath117 for @xmath118 this results in @xmath119 and@xmath120 where @xmath121 is the spherical radius . for a transversal dipole @xmath122 considered also by thess _ _ at al.__@xcite , solution ( [ sol : phi0-h ] ) reduces to @xmath123 where the last term represents the magnetic potential of a dipole aligned with the @xmath0-axis . for a general dipole orientation ,
expression ( [ eq : h ] ) substituted into solution ( [ sol : phi0-h ] ) after some algebra yields@xmath124,\label{sol : phi0}\ ] ] where @xmath20 is the azimuthal angle from the positive @xmath17-axis in the @xmath125-plane .    the electric potential distribution ( [ sol : phi0 ] )
allows us to calculate the force acting upon the magnet , which , as noted above , is opposite to that acting upon the layer , i.e. , @xmath126 where the integral of the lorentz force density @xmath127 is taken over the layer volume @xmath128 for the longitudinal force component , we obtain@xmath129+\int_{h}^{\infty}(\partial_{xm}^{2}g)^{2}\,\d z\right\ } \,\d s\nonumber \\   & = & \frac{\mu_{0}^{2}\sigma v}{512\pi h^{3}}(3m_{x}^{2}+m_{y}^{2}+4m_{z}^{2}),\label{eq : fx}\end{aligned}\ ] ] where the first integral is taken over the surface @xmath130 at @xmath131 and can be swapped with the second one over the layer depth . for a transversal dipole @xmath132 , the expression above coincides with that of thess _ _ et al.__@xcite as well as with the result of reitz@xcite for a thin sheet in the limit of a slowly moving dipole . as seen
, the force is the strongest on a transversal dipole and reduces on longitudinal and spanwise dipoles by factors of @xmath133 and @xmath134 , respectively .
thus , the force ( [ eq : fx ] ) , in contrast to force ( [ eq : f0-gen ] ) for a long magnet , varies with the dipole orientation in the @xmath43-plane . on a horizontally inclined dipole @xmath135
, there is also a spanwise force component@xmath136 but there is no vertical force whatever the dipole orientation:@xmath137 this is consistent with the result of reitz@xcite stating that in the low - speed limit the lift force is proportional to the velocity squared , while the approximation under consideration takes into account only the part of force proportional to the velocity .
alternatively , the force and torque acting on dipole can be found using the associated current distribution ( [ eq : j ] ) and the induced magnetic field @xmath138 which lead straightforwardly to the classical expressions@xcite @xmath139 where the integrals are taken over the space @xmath140 above the layer .
this is a bit longer but algebraically more straightforward approach , which will be pursued in the following .    in order to find the induced magnetic field @xmath138 it is important to notice that solution ( [ sol : phi0-gen ] ) satisfies condition ( [ bc : phi0 ] ) not only at @xmath131 but at any @xmath141 thus
, the @xmath7-component of the current is absent not only at the surface but throughout the whole layer .
then @xmath102 and @xmath142 imply @xmath143 where @xmath144 is the electric stream function , whose isolines coincide with the eddy current lines . substituting this expression into eq .
( [ eq : ohm ] ) and taking the @xmath7-component of the _ curl _ of the resulting equation , we obtain@xmath145 note that in contrast to the electric potential in eq .
( [ eq : phi0 ] ) , no boundary conditions are required for @xmath144 at the surface .
this is because , firstly , eq . ( [ eq : psi - gen ] ) contains no derivatives in @xmath7 and , secondly , the absence of the @xmath7-component of current at the surface is explicitly ensured by expression ( [ eq : j - psi ] ) .    for translational motion with @xmath8 the r.h.s of eq .
( [ eq : psi - gen ] ) takes the form @xmath146 which by the same arguments as above satisfies laplace s equation and equals to @xmath147 with @xmath148 which satisfies the same equation and , thus , represents the solution of eq .
( [ eq : psi - gen ] ) .
for a transversal dipole @xmath149 ) , the general solution ( [ eq : psi0-h ] ) simplifies to @xmath150 where the last term represents the magnetic potential of a longitudinal dipole .
thus , in this case , the eddy current lines coincide with the isolines of electric potential ( [ sol : phi0 ] ) rotated by @xmath151 about the @xmath7-axis . for a longitudinal dipole @xmath152
the comparison of solutions ( [ eq : psi0-h ] ) and ( [ sol : phi0-h ] ) shows that the eddy current lines and the electric potential isolines are swapped with the corresponding distributions induced by a spanwise dipole @xmath153 and rotated by @xmath151 about the @xmath7-axis . for an arbitrarily oriented dipole , expression ( [ eq : h ] )
substituted into the general solution ( [ eq : psi0-h ] ) after some algebra yields@xmath154.\label{sol : psi0}\ ] ] figure [ fig : mov - xyz ] shows the isolines of the electric potential ( [ sol : phi0 ] ) and those of @xmath155 which represent the eddy current lines , in the @xmath125-plane @xmath156 for three basic dipole orientations along the @xmath17- , @xmath0- and @xmath7-axis .
note that the patterns are self - similar in the @xmath125-plane with the characteristic length scaling directly with @xmath141 namely , these distributions are functions of spatial angles only , while according to ( [ sol : phi0 ] ) and ( [ sol : psi0 ] ) the magnitude of both @xmath103 and @xmath157 falls off as @xmath158    isolines of the electric potential @xmath103 ( dashed ) and of the stream function @xmath157 ( solid ) in the @xmath125-plane , where the latter represent the eddy current lines , induced by translation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c ) .
, title="fig:"](-30,0)(a)isolines of the electric potential @xmath103 ( dashed ) and of the stream function @xmath157 ( solid ) in the @xmath125-plane , where the latter represent the eddy current lines , induced by translation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c ) . , title="fig:"](-30,0)(b )    isolines of the electric potential @xmath103 ( dashed ) and of the stream function @xmath157 ( solid ) in the @xmath125-plane , where the latter represent the eddy current lines , induced by translation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c ) . , title="fig:"](-30,0)(c )    in order to satisfy the solenoidality condition @xmath159 the induced magnetic field is sought as @xmath160 where @xmath161 is the vector potential , which is supposed to satisfy the coulomb gauge @xmath162 then ampere s law leads to @xmath163 which , in turn , results in @xmath164 substituting expression ( [ eq : j - psi ] ) into the last integral , after some algebra we obtain @xmath165 where @xmath166 is the same as used by thess _ _ et al.__@xcite the previous expression implies that @xmath167 similarly to its source @xmath168 has no @xmath7-component . then the induced magnetic field can be written as @xmath169 where @xmath170 represents the scalar magnetic potential , which completely defines @xmath171 outside the layer , where @xmath172 further , using the tensor notation with the einstein summation convention , expressions ( [ eq : f ] ) and ( [ eq : m ] ) can be represented in terms of @xmath173 as @xmath174 where @xmath175 is the anti - symmetric tensor and the subscript after the semicolon denotes the differential with respect to the corresponding coordinate with indices @xmath176 @xmath177 and @xmath178 standing for the @xmath6@xmath0 and @xmath179 directions , respectively . taking into account the symmetry of @xmath180 with respect to the interchange of observation and integration points , the derivatives of @xmath173 at the dipole location @xmath181 in expressions ( [ eq : fi ] ) and ( [ eq : mi ] ) are found as @xmath182 where @xmath183 integrals ( [ eq : x - ij ] ) and ( [ eq : x - ijk ] ) can be evaluated analytically using , for example , the computer algebra system mathematica,@xcite which also allows us to carry out all other analytical transformations . in such a way , we firstly verify that eq .
( [ eq : fi ] ) with expression ( [ eq : psi0-h ] ) substituted into integral ( [ eq : x - ijk ] ) indeed reproduces previous results ( [ eq : fx ] ) , ( [ eq : fy ] ) and ( [ eq : fz ] ) . secondly , eq .
( [ eq : mi ] ) with the same expression for @xmath157 substituted in integral ( [ eq : x - ij ] ) results in@xmath184 which for a transversal dipole again coincides with the results of thess _ _ et al.__@xcite      for a solid - body rotation with the velocity @xmath185 which appears in the co - rotating frame of reference , the charge conservation @xmath102 applied to eq .
( [ eq : ohm ] ) results in@xmath186 while the vanishing of the normal current component at the surface @xmath131 requires@xmath187 in order to solve this problem , firstly , it is important to notice that since @xmath188 satisfies the laplace equation , @xmath189 is a particular solution to eq .
( [ eq : phi1 ] ) .
comparing expressions ( [ eq : phi1-b ] ) and ( [ sol : phi0 ] ) shows that @xmath190 where @xmath103 is given by solution ( [ sol : phi0-h ] ) .
searching for the solution as @xmath191 reduces eq .
( [ eq : phi1 ] ) to the laplace equation for @xmath192 while bc ( [ bc : phi1 ] ) takes the form @xmath193 by the usual arguments , we obtain @xmath194 where @xmath195 and @xmath196 is the first order associated legendre function of the second kind.@xcite using the expressions above , after some transformations we obtain@xmath197.\label{sol : phi1}\end{aligned}\ ] ] which is the electric potential in the co - rotating frame of reference .
instantaneous isolines of the electric potential @xmath198 ( dashed ) and of the stream function @xmath199 ( solid ) , which represent the eddy current lines , induced by a solid - body rotation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c).,title="fig:"](-30,0)(a)instantaneous isolines of the electric potential @xmath198 ( dashed ) and of the stream function @xmath199 ( solid ) , which represent the eddy current lines , induced by a solid - body rotation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c).,title="fig:"](-30,0)(b )    instantaneous isolines of the electric potential @xmath198 ( dashed ) and of the stream function @xmath199 ( solid ) , which represent the eddy current lines , induced by a solid - body rotation in the magnetic field of the dipole aligned with the @xmath17- ( a ) , @xmath0- ( b ) , and @xmath7-axis ( c).,title="fig:"](-30,0)(c )    as before , condition ( [ bc : phi1 ] ) is satisfied not only at @xmath131 but any @xmath7 and , thus , the @xmath7-component of the current vanishes throughout the layer .
consequently , the electric current can again be expressed as ( [ eq : j - psi ] ) with @xmath200 for which eq .
( [ eq : psi - gen ] ) takes the form@xmath201 since the r.h.s . of eq .
( [ eq : psi1 ] ) satisfies the laplace equation , after some algebra we obtain@xmath202 the eddy current lines in the @xmath125-plane , which are the isolines of @xmath200 along with the isolines of the electric potential ( [ sol : phi1 ] ) are shown in fig .
[ fig : rot - xyz ] .
again , the patterns are self - similar and scale directly with the distance @xmath7 from the dipole , while the magnitude of both @xmath199 and @xmath198 falls off as @xmath203 in contrast to the translation considered above , the eddy currents induced by rotation are time - dependent and vary periodically as the dipole orientation changes from the @xmath17- to the @xmath7-axis , which are shown in figs . [
fig : rot - xyz](a ) and ( c ) , respectively . only the electric potential , but no current , is induced by a solid - body rotation when the dipole is aligned with the axis of rotation ( @xmath0 ) . in his case , when the magnetic field is symmetric about the axis of rotation , the induced e.m.f . is irrotational and , thus , compensated by the electric potential gradient , which equals to @xmath204 where @xmath205 is the cylindrical radius from the axis of rotation , @xmath206 is the azimuthal unity vector , and @xmath207 is the azimuthal component of the vector potential , which can be used to describe a general axially symmetric poloidal magnetic field . for the dipole field , we have@xmath208 where @xmath209 is the associated dipole current ( [ eq : j ] ) and @xmath210 substituting this into expression ( [ eq : v1xb ] ) and equating it to the gradient of the electric potential @xmath211 we obtain @xmath212 which is the free - space electric potential in the co - rotating frame of reference shown in fig .
[ fig : rot - xyz](b ) and coinciding with ( [ sol : phi1 ] ) when @xmath213 this potential vanishes in the laboratory frame of reference , where the magnetic field is invariant with respect to the rotation around the symmetry axis .
further , using expression ( [ sol : psi1 ] ) for @xmath199 in integrals ( [ eq : x - ij ] ) and ( [ eq : x - ijk ] ) , which can be evaluated together with eqs .
( [ eq : fi ] ) and ( [ eq : mi ] ) in the same way as for the translation in the previous section , we obtain @xmath214    the first point to note is that both the force and torque vanishes when dipole is aligned with the axis of rotation @xmath215 which , as discussed above , is due to the absence of eddy currents in this case .
secondly , the longitudinal @xmath216 component of force ( [ eq : f1 ] ) , in contrast to that due to the translation given by expression ( [ eq : fx ] ) , is independent of the dipole orientation in the @xmath43-plane and depends only on the magnitude of the dipole moment in this plane : @xmath217 thirdly , there is also a non - zero transversal @xmath218 force , which in contrast to the longitudinal one is purely oscillatory and varies periodically with the dipole orientation in the @xmath43-plane as @xmath219 where @xmath16 is the poloidal angle of the dipole orientation in the @xmath43-plane from the positive @xmath17-axis , which is the same as that for the linear dipole in expression ( [ eq : a2d ] ) ( see fig .
[ cap : sketch ] ) .
moreover , when dipole is inclined to the axis of rotation ( @xmath220 also a spanwise @xmath221 force appears , which similarly to the transversal one alternates periodically around zero as @xmath222 with the dipole rotation . as seen from expression ( [ eq : m1 ] ) , inclination also gives rise to alternating torque components around both the @xmath17- and @xmath7-axis .
the @xmath0-component of torque , which is negative and , thus , opposing the rotation , has not only a constant but also an alternating part , which varies with the double frequency of rotation as @xmath223    instantaneous rotation rate ( [ eq : omgt - trq ] ) versus time ( solid ) and orientation ( dashed ) for a negligible inertia of the magnet.,scaledwidth=50.0% ]    thus , in contrast to the long magnet in @xmath224[sub:2d ] , there is no steady rotation rate , which could balance the @xmath0-component of the constant torque ( @xmath225 ) due to the translation .
the variation of the magnet rotation rate in response to oscillatory torque is constrained by the inertia of magnet . for the oscillatory part of the rotation rate @xmath226 to be negligible compared to the mean one @xmath227 we need @xmath228 where @xmath229 is the moment of inertia of magnet . if this condition is satisfied , the oscillatory component may be neglected in the balance of torques ( [ eq : m0 ] ) and ( [ eq : m1 ] ) around the axis of rotation , which is the @xmath0-axis . for a layer of finite thickness
, this results in @xmath230 the expression above differs only by a factor of @xmath231 from result ( [ eq : omg - frc ] ) for a long magnet in the case of a vanishing longitudinal force . in the opposite limit of a negligible inertia , the torque balance yields@xmath232 this instantaneous rotation rate @xmath233 which is shown in fig . [ fig : omg - t ] versus both the time @xmath234 and orientation @xmath235 alternates between @xmath236 and @xmath237 of the mean value ( [ eq : omgb - trq ] ) @xmath238 the last relation follows from the integration of eq .
( [ eq : omgt - trq ] ) as @xmath239 which also defines @xmath87 parametrically versus @xmath234 in fig .
[ fig : omg - t ] .
note that the oscillations of the rotation rate do not affect its mean value @xmath240 which is defined by ( [ eq : omgb - trq ] ) independently of condition ( [ eq : omeg - osc ] ) .
alternatively , when the magnet rotation rate is actively controlled so that to balance the longitudinal force ( [ eq : fx ] ) with the @xmath17-component of force ( [ eq : f1 ] ) , the instantaneous rotation rate is @xmath241 where the temporal mean value @xmath242 follows from the solution of differential equation ( [ eq : omgt - frc ] ) , which can be written as @xmath243
single - magnet rotary flowmeter with a _ smco_-type magnet ( 1 ) of @xmath244 and @xmath245 in diameter and height , respectively , installed at the side of a rectangular stainless steel duct ( 2 ) of @xmath246 in cross - section carrying the liquid sodium flow .
the rotation rate is measured with an inductive magnetic proximity sensor ( 3 ) and also monitored with a photoelectric barrier sensor ( 4 ) .
, scaledwidth=50.0% ]    the laboratory model of the single - magnet flowmeter shown in fig .
[ exp - pic ] consists of a cylindrical _ smco_-type permanent magnet with diameter @xmath247 and height @xmath248 , magnetized perpendicularly to its axis with @xmath249 surface induction , which holds up to @xmath250 the magnet is mounted on a stainless steel axle held by a housing made of `` macor '' ceramic , which ensures a low mechanical friction and can withstand up to @xmath251    the flowmeter was tested on a sodium loop using about 90 liters of molten _ na _ at @xmath252 with the electrical conductivity of @xmath253 the flow was driven by a linear electromagnetic pump , which provided a maximum velocity of @xmath254 at the flow rate of @xmath255 in a stainless steel duct of @xmath246 cross - section .    the frequency of magnet rotation was determined with an inductive magnetic proximity sensor ( sick mm12 - 60aps - zu0 ) , which was fixed at the distance of @xmath256 perpendicularly to the side surface of magnet at its midheight as shown in fig .
[ exp - pic ] .
the sensor briefly switched off and then on again as the magnetic field along the axis of the sensor changed its direction .
the frequency of sensor transition from the off to the on state , which happens twice per revolution , was measured with keithley 2000 multimeter using the reciprocal frequency counting techniques with the measurement gate time set to @xmath257 the frequency was also monitored with tektronix tds 210 oscilloscope . after changing the pump current ,
the flow was allowed to develop for @xmath258 and then six measurements were taken with the intervals of approximately @xmath259 the standard deviation in the measured rotation frequencies was typically a few tenths of @xmath260    the rotation rates measured at two gap widths between the magnet and the duct are shown in fig .
[ fig : exp](a ) depending on the average sodium velocity , which was determined using ultrasonic doppler velocimetry@xcite with the resolution of @xmath261@xcite although the rotation rate is seen to increase as the magnet is approached to the duct and to vary nearly linearly with the flow velocity in agreement with the theory , the linear fit shows a certain zero offset .
obviously , there is some opposing force , which has to be overcome for the rotation to start .
this requires the flow velocity of at least @xmath262 first , such an opposing force is caused by the dry friction in bearings .
second , the magnet is oriented by the earth s magnetic field , which has a stronger effect than the friction and also has to be overcome for the rotation to start .
third , even stronger effect is caused by the stainless steel walls of the duct , which are weakly ferromagnetic.@xcite the magnet , when approached to the duct , was observed to turn with the dipole moment perpendicular to the wall .
it is also important to note that the magnet , when turned by hand away from the equilibrium orientation , returned to it monotonically without oscillations .
this implies the inertia of magnet to be small relative to the electromagnetic drag force .
thus , the inertia can not significantly contribute to the overcoming of the orienting torque due to the wall magnetization .
however , when magnet rotates slowly , the orienting torque causes the rotation rate to fluctuate , which shows up as an increased scatter in the measured values seen in fig .
[ fig : exp ] at low flow velocities .    when the zero offset is removed and the remaining rotation rate is multiplied by the distance @xmath263 between the magnet axis and the liquid metal , where @xmath264 is the magnet radius and @xmath265 is the thickness of stainless steel wall , we obtain the relative rotation rate @xmath266 which represents a rescaled slope coefficient for fig .
[ fig : exp](a ) and is plotted in fig . [
fig : exp](b ) versus the velocity .
the short - magnet solution ( [ eq : omgb - trq ] ) in two limiting cases of a semi - infinite @xmath267 and a thin ( @xmath268 layer yields @xmath269 and @xmath270 respectively , which are smaller than the measured values @xmath271 seen in fig .
[ fig : exp](b ) . on the other hand ,
the long - magnet solution ( [ eq : omg - trq ] ) , which appears more adequate for the experimental setup , yields @xmath272 for @xmath273 such quantitative differences between the experiment and theory are not surprising given the simplifications underlying the latter .
first , theoretical model does not take into account the finite width of the duct .
but this alone could hardly explain the high rotation rate of the magnet at which its field travels faster than the layer .
this apparently being the case in the experiment implies the presence of significant velocity gradients , which are also ignored by the theory , but could be taken into account as outlined in sec .
[ sub:2d ] . note that strong velocity gradients can be caused by the magnet itself , which due to its large size and strength may act as a magnetic obstacle partially blocking the flow and so increasing its local velocity .    measured ( a ) and rescaled ( b ) magnet rotation rates along with the best linear fits versus the average velocity for @xmath274 and @xmath275 gaps between the magnet and the duct with the distance between the magnet axis and liquid metal @xmath276 where @xmath264 is the magnet radius and @xmath265 is the thickness of stainless steel wall.,title="fig:",scaledwidth=50.0%](-30,30)(a)measured ( a ) and rescaled ( b ) magnet rotation rates along with the best linear fits versus the average velocity for @xmath274 and @xmath275 gaps between the magnet and the duct with the distance between the magnet axis and liquid metal @xmath276 where @xmath264 is the magnet radius and @xmath265 is the thickness of stainless steel wall.,title="fig:",scaledwidth=50.0%](-30,30)(b )
we have presented a theory of single - magnet rotary flowmeter for two limiting cases of long and short magnets , which were modeled as linear and single dipoles , respectively .
simple analytical solutions were obtained for the force and torque on slowly translating and rotating magnets due to eddy currents in layers of arbitrary depth and infinite horizontal extent .
the velocity was assumed to be constant and the motion so slow that the induced magnetic field could be neglected .
the latter assumption was not applicable to the long magnet rotating above a conducting half - space . in this case , to obtain a finite braking torque , the skin effect due to the induced magnetic field had to be taken into account however slow the rotation . for a single dipole of arbitrary orientation ,
compact analytical solutions were obtained in terms of both the electric potential and stream function induced by the layer translation and rotation in the co - rotating frame of reference .
the electric stream function was used further to find the scalar potential of induced magnetic field at the dipole location , which resulted in simple expressions for the force and torque on the dipole .
eventually , we found the equilibrium rotation rate at which the driving torque due to the layer translation is balanced by the braking torque due to the magnet rotation .
an alternative approach was also considered , where the force on the magnet could be used to control its rotation rate so that the resulting force vanishes . in either case
, the resulting rotation rate is directly proportional to the layer velocity and inversely proportional to the distance between the magnet and the liquid metal .
these results were found in a qualitative agreement with the measurements on the liquid sodium flow . a more accurate quantitative agreement with experiment
is limited due to the substantial approximations underlying the theoretical model , which neglects the finite lateral extension of the layer as well as the spatial and temporal variations of the velocity distribution .    in conclusion ,
note that the resulting rotation rate is independent of the magnet strength and the electrical conductivity of the liquid metal provided that the mechanical friction or other external effects are negligible compared to the driving torque .
this is the main advantage of rotary flowmeter over the lfv.@xcite              x. fu , l. hu , k. m. lee , j. zou , x. d. ruan , and h. y. yang , dry calibration of electromagnetic flowmeters based on numerical models combining multiple physical phenomena ( multiphysics ) , _ j. appl .
phys . _ * 108 * , 083908 - 11 ( 2010 ) .
i. bucenieks , modelling of rotary inductive electromagnetic flowmeter for liquid metal flow control , in : _ proceedings of the 5th int .
symposium on magnetic suspension technology , _ 204208 ( dresden , germany , 2005 ) .
g. gerbeth , j. priede , i. bucenieks , a. bojarevics , j. gelfgat , and d. buchenau , verfahren und anordnung zur kontaklosen messung des durchflusses elektrisch leitfhiger medien , _
german patent application _ de 102007046881a1 ( 2009 ) .
j. priede , d. buchenau , and g. gerbeth , force - free and contactless sensor for electromagnetic flowrate measurements , in : _
pamir conference _ , 815820 ( presqule de giens , france , 2008 ) ; _ magnetohydrodynamics _ * 45 * , 451458 ( 2009 ) .
s. eckert , d. buchenau , g. gerbeth , f. stefani and f .-
weiss , some recent developments in the field of measuring techniques and instrumentation for liquid metal flows , _ j. nucl .
tech . _ * 48 * , 490498 ( 2011 ) .