arxiv: v1 [physics.flu-dyn] 11 Mar 2011

Transcription

1 arxiv: v1 [physics.fl-dyn 11 Mar 011 Interaction of a magnetic dipole with a slowly moving electrically condcting plate Evgeny V. Votyakov Comptational Science Laboratory UCY-CompSci, Department of Mechanical and Manfactring Engineering, University of Cyprs, 75 Kallipoleos, Nicosia 1678, Cyprs karal@gmail.com André Thess Institte of Thermodynamics and Flid Mechanics, Ilmena University of Technology, P. O. Box , Ilmena, Germany thess@t-ilmena.de Abstract. We report an analytical investigation of the force and torqe acting pon a magnetic dipole placed in the vicinity of a moving electrically condcting nonmagnetic plate. This problem is relevant to contactless electromagnetic flow measrement in metallrgy and extends previos theoretical works (Thess et al. Phys. ev. Lett., 96(006), ; New J. Phys. 9(007), 99) to the case where the orientation of the magnetic dipole relative to the plate is arbitrary. It is demonstrated that for the case of low magnetic eynolds nmber the three-dimensional distribtions of the indced electric potential, of the eddy crrents and of the indced magnetic field can be rigorosly derived. It is also shown that all components of the force and torqe can be compted withot any frther approximation. The reslts of the present work serve as a benchmark problem that can be sed to verify nmerical simlations of more complex magnetic field distribtions.

2 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 1. Introdction When an electrically condctive macroscopic body moves in the nonniform magnetic field B created by an external sorce, for instance a permanent magnet, then eddy crrents j are indced in the body. These eddy crrents create the Lorentz force F = j B. The force is directed in the direction opposite to the direction of the body s translation. Moreover, the eddy crrents indce an additional magnetic field b which interacts with the permanent in sch a way as to create a force and a torqe acting pon the magnet. Ths, by measring these forces and torqes one can determine the parameters of the movement sch as the velocity or the direction. This principle is embodied in a contactless electromagnetic flow measrement techniqe called Lorentz force velocimetry [8, [10, [9, [1, [, [6, which permits flow measrement in hot and aggressive flids like liqid alminim or molten steel. When developing measrement systems embodying Lorentz force velocimetry, so-called Lorentz force flowmeters, it is necessary to predict the force and torqe acting pon a complex-shaped magnet system by a trblent liqid metal flowing in pipes, dcts and open channels. Sch predictions are sally performed by nmerically solving the fll set of three-dimensional eqations of magnetohydrodynamics [3. In order to be able to assess the reliability of sch simlations it is necessary to have simple models that are amenable to rigoros analytic treatment. The goal of the present paper is to formlate and solve sch a model. The model to be stdied in the present work is a generalization of a previosly stdied problem [9 to the case when the orientation of a magnetic dipole is arbitrary. More precisely, the athors of [9 investigated the interaction between a moving plate and a magnetic dipole whose orientation is perpendiclar to the srface of the plate. Here we relax the assmption abot the orientation of the dipole and allow the dipole to be arbitrarily oriented. As will be shown in Section, this problem can still be solved exactly and all components of the torqe and the force can be expressed analytically along with the electric and magnetic field. We illstrate or reslts in Section 3 sing some representative plots of the three-dimensional strctre of the eddy crrents and the indced magnetic fields. In Section 4 we smmarize or conclsions and indicate a possible illstrative application of the theory that wold be interesting to investigate.. The problem and its exact soltion Consider a single magnetic dipole with dipole moment m located at a distance h above a flid or solid layer with electrical condctivity σ and thickness w as shown in Fig. 1. The dipole is assmed to be arbitrary oriented in space. For brevity, if it is not specified otherwise, we assme that the index n denotes one of the Cartesian axes, m = me n. The general case can be considered as a linear combination of the three generic cases, n = x,y,z, any general orientation a can be expressed as a = m {a (n) } T, where {a (n) } = (a (x),a (y),a (z) ) is the tensor containing components for each n = x,y,z dipole orientation. The layer moves with an niform velocity = e x and extends in the z

3 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 3 Figre 1. Sketch of the problem: A plate with thickness w moves with constant velocity and interacts with an arbitrarily oriented magnetic dipole. direction from z d to z. The thickness of the layer is then given by w = z z d. Here, the indices d and stand for downward side and pward side, respectively. The borders of the layer in the x and y direction are spposed to be at infinity. This corresponds to the assmption that the distance h is mch less than the horizontal size of the moving solid layer. We frther assme that plate or the flid moves slowly which is eqivalent to the assmption that magnetic eynolds nmber m = µ 0 σw is small. In magnetohydrodynamics (see e.g. [7, [3) this case is referred to as the qasistatic approximation and implies that the magnetic field associated with the eddy crrents is mch smaller than the applied magnetic field. The assmption m 1 makes the present problem amenable to rigoros analytic treatment. For instance, the magnetic eynolds nmber of liqid alminim (σ ) flowing with a velocity = 1 m/s in a pipe with a diameter of 0.1 m is approximately 0.3. In what follows we will se mainly Cartesian coordinates and will not apply an integration of Bessel fnctions in cylindrical coordinates as we did before in [9. Or goal is to derive explicit compact algebraic expressions, which are sefl for fast analytic calclations, for the primary magnetic field B(r), the electric potential φ(r), the eddy crrents j(r) and the secondary magnetic field b(r). Then, all other qantities sch as forces and torqes acting pon the dipole can be readily written down..1. Primary magnetic field The first step in the soltion of or problem is the specification of the magnetic field of the dipole which is referred to as the primary magnetic field. Let the magnetic dipole be at the point r 0 = (x 0,y 0,z 0 ) and we wish to calclate the magnetic field B(r) at a point r = (x, y, z). Before presenting formlae for the primary magnetic field, we introdce the following short notations of the vector = r r 0 : x = x x 0, y = y y 0, z = z z 0, r = x + y, = r + z,

4 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 4 1 = [ [ tanh 1 z 1, = z 3 (1) Now, the magnetic field of a magnetic dipole can be represented as [5: B(r) = µ 0 [m 4π = µ [ ( 0 3 4π m 1 ) = µ [ ( 0 m 4π ) 1 m = µ [ ( ) 0 1 4π m, () These expressions appear more complicated than the formla normally sed. However, the present notation will trn ot to be convenient for the frther analysis and it is the keytoacompact representation ofallfields. Thenitfollows, bytaking(1)andm = me n into accont, that the components of the primary magnetic field can be conveniently given as: B k = µ k n [ 1 = µ k n z [ tanh 1 z, k = x,y,z, (3) where we introdced the abbreviation µ = mµ 0 /4π, and n (.) is the partial derivative with respect to the coordinate specified by the direction n of the dipole. If the dipole m is arbitrary oriented in space, m = m x e x +m y e y +m z e z, then n m. For instance, for the dipole oriented parallel to the z-axis, n = e z and n = z, and the components of the primary magnetic field are:.. Electric potential B x = 3µ x z 5, B y = 3µ y z 5, B z = µ [ 3 z (4) The derivation of the electric potential is an extension of the method developed earlier in [9, and it is presented also in [6. The electric potential φ and eddy crrents j are governed by Ohm s law for an electrically condcting material moving with velocity, which can be written as: j = σ( φ+ B), j = 0. (5) (The second expression is not Ohm s law and does imposes the reqirement that there is no sorce or sink of electric crrents in the moving plate.) In general, the magnetic field in (5) is the sm of the primary magnetic field and the magnetic field associated with the eddy crrents, referred to as the secondary magnetic field. As explained earlier, we assme that the magnetic eynolds nmber m = µ 0 σvh, which characterizes the ratio between the secondary and primary field, is small, so we can approximate the magnetic field in Ohm s law by the primary field given by (3). The primary magnetic field is solenoidal, B = 0, and the velocity is constant. By applying the divergence operator to Ohm s law (5) we therefore obtain the following eqation for the electric potential: φ = [ B = B [ [ B = 0. (6)

5 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 5 Hence, φ is a harmonic fnction obeying Laplace s eqation, φ = 0. Any harmonic fnction in three-dimensional space can be expressed throgh derivatives or integrals of the fnction 1/ defined in (1), where the specific expression for φ depends on bondary conditions. In or case, the domain is infinitively large in x and y directions, and finite in the z direction. At the pper (z = z ) and lower (z = z d ) srfaces of the plate, the vertical electric crrent mst be zero, which is expressed as j z z=zs = 0, and hence, σ( z φ+b y ) z=zs = 0, where the index s stands for either s = or s = d. Hence, we can directly express z φ by setting j z = 0 everywhere in the layer: z φ = B y = µ y n z [ tanh 1 z. (7) where the last eqality directly follows from (3). Now, if we change the order of differentiation and pt z as the first derivative, then the integration is trivial, and we come immediately to the final reslt: φ(x,y,z) = µ { y n [ tanh 1 z } +φ (x,y), (8) where the second term φ (x,y) appears as constant of the integration and so it does not depend on z. This term is determined from the additional bondary condition specifying that the electric potential vanishes when the dipole is at infinite distance from the layer. The reslt is: [ φ (x,y) = lim 1 z y n tanh z 0 { [ } = n lim y tanh 1 z z 0 { = n sgn( z) y } = sgn( z) y n lnr, (9) r where sgn(z) = z/ z is the sign fnction with the property that sgn( z) = 1 for z > z 0 and sgn( z) = 1 for z < z 0. One can readily verify that φ given by (8) is the soltion of (6). To do this, we notice that φ = φ(1/), where 1/ is given by (1) and [1/ = [tanh 1 ( z/) = 0, hence, φ = 0 too. In Appendix A it is shown that the components of electric crrent can be expressed throgh a stream fnction ψ defined as j = (ψe z ). This stream fnction is given by: [ ψ(x,y,z) = µσ { x n ψ (x,y) = lim z 0 x n } +ψ (x,y), (10) = sgn( z) x n lnr. 1 z tanh [ tanh 1 z In the short notation both φ and ψ can be conveniently represented by means of the differentiation of the same axiliary fnction Φ as: φ (n) = µ y n Φ, ψ (n) = µσ x n Φ, (11) { } Φ = sgn( z) tanh 1 z +lnr,

6 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 6 where an pper index n specifies explicitly the direction of the magnetic dipole. Becase the electric crrents are two-dimensional, the soltion (11) describes an electric potential and crrent stream fnction for any horizontal plane of the moving layer. Explicit formlae for ψ and j are given in Appendix A. These formlae can be immediately sed for calclation withot performing cmbersome differentiation..3. Secondary magnetic field and scalar potential The eddy crrents j of the moving plate indce a magnetic field b which we refer to as the secondary magnetic field. When j is time-independent as in or case, the secondary magnetic field b(x,y,z) at a point r = (x,y,z) is niqely determined by Ampere s law j(x,y,z) = 1 µ 0 b(x,y,z). (1) On the other hand, j can be determined with the aid of the stream fnction which has been already above calclated, Eq.(11): j(x,y,z) = Θ(z) [e z ψ(x,y,z). Here the axiliary step fnction Θ(z) is introdced in order to ensre that eddy crrents otside of the plate are zero. For instance, in the case of the layer of finite thickness, the fnction Θ(z) is defined as: Θ(z) = { 1 for z d z z ; 0 for z < z d and z < z, where z d and z are the lower and pper borders of the plate. In the general case, the secondary field can be represented completely as a sm of the stream fnction and a gradient of a scalar fnction a in the form: b(x,y,z) = Θ(z)µ 0 ψ(x,y,z)e z +µ 0 a(x,y,z), (13) Fnction Θ(z) enforces the condition that otside the plate the secondary field b is solely a gradient of the scalar potential a. This atomatically ensres that b = 0, and then the problem how to calclate b is redced into the problem how to calclate the scalar potential a. In order to accomplish this task one has to employ the constraint that b is divergence-free. By taking the divergence of (13), the divergence-free property, b = 0, reslts in the following eqation: [ 1 a = Θ(z) z ψ = Θ(z)µσ x n, (14) where the second eqality directly follows from Eq. (11). So one has to solve now eqation(14) with the bondary conditions a 0 which describes that the secondary magnetic field vanishes at infinity. It is possible to solve Eq. (14) in a straightforward way for any specific n = x,y,z. This reqires explicit partial derivatives with respect to x and n and then cmbersome symbolic calclations for each separate index n by recriting series expansion with Bessel fnctions and frther integration. Sch procedre has been performed in [9 for the

7 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 7 particlar case n = z, see also [6. However, there exists another techniqe to provide a general formla for arbitrary general n. For this we notice that a derivative of any fnction f(), which depends on solely, with respect to n = x,y,z can be taken as a derivative, with the opposite sign, with respect to n 0 = x 0,y 0,z 0, correspondingly, n f() = n0 f(), where is defined in (1). In order to se this property of we represent the scalar potential a by means of a generating (or primitive) fnction A twice differentiated with respect to the coordinates of the dipole, a = µσ x0 n0 A. We insert the above formla for a into (14) and obtain the simpler Poisson s eqation independent of n to determine A as: A = Θ(z) 1. (15) This eqation can be readily solved withot cmbersome integration and series expansion as we shall demonstrate below. Once we obtain the primitive fnction A, the secondary magnetic field can be represented by means of triple differentiation in the form: b (n) = µ 0 µσ [ x0 n0 A. (16).4. Scalar potential and secondary magnetic field for an infinitely thin layer Becase the problem nder consideration is linear, we treat the case of an infinitely thin layer separately. Once a soltion for this case is fond, all other cases can be presented throgh a sperposition of thin layers. Let ẑ be a location inside a thin moving layer with electric crrents indced by an external magnetic dipole. Let s assme that the thickness of the layer is small compared to the distance between the dipole and the layer which is expressed by the condition w = z z d h. Finally let s introdce ẑ = z +z d /. In this case, the fnction Θ(z), which is needed to set the position of the layer, mst be taken in the limit z d z. It gives Θ(z) = wδ(z ẑ), where w is the thickness of the layer which is inserted for the proper physical dimensionality. (Fnction Θ(z) is dimensional in the case of the infinitely thin layer, bt w will be omitted for other cases becase there the proper dimensionality follows directly from an integration over the plate thickness.) Let A 0 be the primitive fnction for the thin layer. Here the index 0 is inserted to distingish the present case from other cases. As follows from Eq.(15), everywhere in the space A 0 = 0 except for z = ẑ. This means that A 0 is a continos harmonic fnction showing a discontinity in the first derivative with respect to z at z = ẑ. As shown in Appendix B, the soltion of (15) is: A 0 (ẑ) = w { 1 ˆ z(ẑ) tanh +lnr ˆ(ẑ) ˆ(ẑ) = r + ˆ z(ẑ), ˆ z(ẑ) = ẑ z 0 + ẑ z. }, (17) Notation ˆ z(ẑ) looks gly bt this is becase we want to stress an analogy in between fnctions and ˆ as explained below. In the given case of a infinitely thin layer, it

8 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 8 is sperflos to pt argment ẑ in the rond brackets of the fnction ˆ z, bt this argment will be necessary afterwards, when we shall extend these reslts to the case of a layer with finite thickness. Note, thata 0 in(17)andφin(11)arealmostidentical, thatis, (w/)φcorresponds to sgn(ˆ z)a 0, with the mapping ˆ and z ˆ z. Moreover, frther analysis shows that this correspondence is valid not only for the entire fnctions A 0 and Φ, bt also for all their derivatives with respect to n = x,y,z and n 0 = x 0,y 0,z 0, that is, (w/) n Φ sgn(ˆ z) n A 0, n Φ = n0 Φ, and n A 0 = n0 A 0. The only exception is z A 0 = z0 A 0 becase (z z 0 ) is taken in A 0 as an absolte vale. By taking this property into accont and recalling that both the eddy crrents j and secondary magnetic field b can be represented throgh a triple differentiation as: j µσ = [ x n Φe z, b µ 0 µσ = [ x 0 n0 A 0. Therefore withot compting derivatives, it is possible to bild p the following correspondence between the secondary magnetic field and eddy crrents: b x (n) sgn(ˆ z) j y (n), b (n) y sgn(ˆ z) j x (n) for n = x,y (18) b x (z) sgn(ˆ z) j y (z), b(z) y sgn(ˆ z) j x (z) for n = z. (19) (The difference in sign for n = z is becase z Φ = z0 Φ while z A 0 = z0 A 0.) The remaining z-components of the secondary field are: b (x) z sgn(ˆ z) b z (z) sgn(ˆ z) j (z) y, b (y) { j (x) y j x (y) z sgn(ˆ z) j x (z), }. (0) The last relationship for b z (z) is the conjectre of the fact that the Φ is a harmonic fnction, i.e. Φ = 0, hence, (j y (x) j x (y) )/(µσ) = x ( x x + y y )Φ = [ x z z Φ and (w/) x z z Φ sgn(ˆ z) x z z A 0. According to the mapping (18-0), the tensor {b (n) } of the secondary magnetic field can be now written down by means of the tensor {j (n) } given in (A.10) for eddy crrents. This tensor is presented in Appendix C and can be straightforwardly sed in algebraic calclations in contrast to the implicit reslts pblished by Priede [6..5. Secondary magnetic field for half-space and for the layer with finite thickness Becase the problem nder consideration is linear, the secondary field of any moving plate can be represented by smming all the fields originating from the layers of the plate derived in the preceding section. Mathematically, this means an integration of the secondary magnetic field tensor, cf. Eq. (C.1), over the thickness of the plate. The limits of the integration in the case of half-space are z and z d. Then, the integration shold be carried ot for the P k coefficients only becase the fnctions f k and g k do not

9 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 9 depend on z. By performing this procedre formally, we obtain P,k = z P 0,k (ẑ) dẑ, (1) where P 0,k and P,k are tensor coefficients of an infinitely thin layer and half-space, correspondingly. Explicit formlae for P,k are presented in Appendix D. When the plate is of finite thickness, z d ẑ z, the coefficients P k are obviosly expressed throgh the coefficients for half-space in the form: z z zd P k (z,z d ) = P 0,k (ẑ) dẑ = P 0,k (ẑ) dẑ P 0,k (ẑ) dẑ z d = P,k (z ) P,k (z d ), () where P,k (z ) and P,k (z d ) are compted according to (D.1) with ˆ z and ˆ taken at points z and z d..6. Magnetic potential energy, force and torqe In the general case, the potential energy U of a magnetic dipole m in a magnetic field B is U = m B, see e.g. [4. In or case, the magnetic field is the indced secondary magnetic field b = µ 0 µσ [ m {b (n) } T m, hence, the potential energy of the dipole is written down as: U = µ 0 µσ [ m {b (n) } T m. = µ 0 µσ [ m x b(x) x +m y b(y) y +m z b(z) z +m x m y b x (y), (3) where the last relationship is obtained by taking into accont a symmetry of the secondary field tensor b (x) y = b x (y), b y (z) = b z (y), and b x (z) = b (x) z. Force F acting pon the dipole can be calclated as a negative gradient of the potential energy: [ F = U, F i = µ 0 µσ m {f (n) i } T m, (4) where a force field tensor {f (n) i } is a partial derivative of the secondary magnetic field tensor with respect to the coordinate axis i = x,y,z, i.e.: {f (n) i } = i {b (n) }. (5) Becase the tensor {b (n) } has been fond above, algebraic formlae for the tensor {f (n) i } can be compted as well. They are presented by Eq.(E.1-E.3). Notice, that these tensors contain 7 terms (3 times 9), while de to symmetry and algebraic simplifications they can be given by means of linear combinations of five coefficients Q k, k = 0,...,4. The torqe T acting pon the dipole is: T = m b = µ 0 µσ ( m [ {b (n) } T m ). (6) It can be compted also also by means of the tensor {b (n) } presented explicitly in (C.1). Above reslts for the potential energy, force and torqe are of the same strctre when the secondary magnetic field is indced by a discrete ensemble of magnetic dipoles

10 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 10 or by a continos macroscopic magnet. (The latter can be obtained by an integration over the space occpied by the magnet.) However in or simplest case, in order to calclate the force and torqe acting pon a single dipole, we jst pt r = 0 and z = z 0 in all the tensors fond already. Mathematically, these tensors are sms of the terms, which have fnctions f k and g k as factors with coefficients P k or Q k. The fnctions f k and g k are eqivalent to coskφ and sinkφ, correspondingly, here, tanφ = y/ x, and, hence, f k and g k are zero at r = 0 for k 1. Hence, at r = 0, all the tensor terms having factors P k or Q k with k 1 mst vanish de to their prodcts with f k and g k. This greatly simplifies all the formlae, and we obtain finally: b = µ 0 µσ P 0 r=0,z=z0 { m z e x +m x e z }, (7) F = µ 0 µσ Q 0 r=0,z=z0 { (3m x +m y +4m z)e x +m x m y e y }, (8) T = µ 0 µσ P 0 r=0,z=z0 { mx m y e x (m x +m z )e y +m y m z e z }, (9) where the non-zero tensor coefficients P 0 and Q 0 at r = 0 and z = z 0 are presented in Appendix F in formlae (F.1, F.) for a thin layer, (F.3, F.4) for a half space, and (F.5, F.6) for the layer with finite thickness. 3. Spatial strctre of eddy crrents and secondary magnetic field In order to illstrate the analytical reslts of previos sections, we plotted eddy crrents and the secondary magnetic field by sing the formlae derived above, Eq. A.6) and (C.1). The eddy crrents below are shown as contor lines of the eddy-crrent stream fnction ψ given by Eq. (A.6). Solid lines correspond to positive ψ and indicate conterclockwise direction of the eddy crrent; whereas dotted lines are for negative contor levels and imply clockwise direction of the eddy crrent. The secondary magnetic field is given in the form of selected three-dimensional magnetic field lines compted by sing Eq. (C.1). According to the right-hand rle, the secondary magnetic field lines emanate from the solid positive loops and penetrate the negative dotted loops. Figre shows behavior at m y = 0 and continosly changing m x and m z provided that m x + m y + m z = 1. This corresponds to the case when the dipole is confined to the plane spanned by the velocity vector and the direction is normal to the srface of the plate. The reslts of Fig. (a) have already been given in [9. In this case, m x = m y = 0,m z = 1, and there are two loops of eddy crrents whose centers are arranged parallel to the direction of the plate movement. The secondary magnetic field lines form a cage. When m x contribtion increases, both loops deform and the positive loop grows in its size and shifts nder the dipole as shown in Fig. (b). Then, the case of Fig. (c) corresponds to a magnetic dipole, which is oriented along the x-direction, i.e. m x = 1, m y = 0, m z = 0, where the positive eddy crrent loop is located completely nder the dipole, and two negative loops are at the periphery. As a reslt, the threedimensional lines of the secondary magnetic field are either completely in the y z plane, if they start at x = 0, or the lines are closed at the peripheral negative loops if they start at y = 0. This behavior is emphasized by Fig. (d) which refers to the same

11 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 11 parameters as Fig. (c) bt plotted from a different perspective. Then, Fig. (e) is a mirror of Fig. (b) with the difference that the directions of the eddy crrents and of the secondary magnetic field are opposite. The same is tre for Fig. (f) as compared to Fig. (a). Figre 3 shows the behavior at m z = 0 and continosly changing m x and m z provided that m x + m y + m z = 1. This corresponds to the case when the magnetic dipole is located in a plane parallel to the srface of the layer. The case of Fig. 3(a), m x = 1, m y = 0, m z = 0 is the same as Fig. (c) and is repeated here in order to gide the eye. By rotating the dipole, the central loop trns together with the lines of the secondary magnetic field, Fig. 3(b). Then, when the main axis of the central loop is parallel to the diagonal of the plate, it splits into two eqal positive loops, and ths one observes for alternating eddy crrents loops for m x = 0, m y = 1 m z = 0, see Fig. 3(c). The magnetic lines emanate from the positive loop and go into the neighboring negative loop. If one looks at the system from above, as in Fig. 3(d) being the same as Fig. 3(c) bt differently oriented, one can see that the line r = 0, where the dipole is located, is the place where the magnetic field lines meet and trn ot, holding the secondary field eqal always zero. Then, the Fig. 3(e) is a mirror of Fig. 3(b) and Fig. 3(f) is a mirror of Fig. 3(a) with the difference that directions of eddy crrents and secondary field are opposite in respect to each other. 4. Conclsion An analytic soltion is obtained for the electromagnetic interaction between an arbitrarily oriented magnetic dipole and a moving electrically condcting solid layer. The soltion incldes the electric potential and eddy electric crrents inside the layer as well as the scalar potential for the secondary magnetic field otside the layer. The formlae obtained are sefl for a Lorentz force velocimetry, becase a moving solid plate can be considered as the limiting case of the mean velocity profile of a trblent liqid metal flow when the eynolds nmber (e) tends to infinity. It has been shown in [9 that the signal of a Lorentz force velocimetry system interacting with a trblent pipe flow converges towards the signal generated by a moving solid cylinder when e. Hence the present model provides a benchmark for nmerical simlations of Lorentz force velocimetry in channel flows at high eynolds nmbers. It shold be finally mentioned that the present model cold be applied in a straightforward way to describe the behavior of an edcational experiment which consists of a spherical permanent magnet rolling down an inclined alminim plate. Depending on the initial orientation of the magnetization the sphere ndergoes different kinds of motion inclding straight motion with periodically changing velocity, wiggling motion and irreglar motion. With minor changes the present model cold be extended to cover this case.

12 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 1 Acknowledgements We are gratefl to the German esearch Fondation(Detsche Forschngsgemeinschaft) for spporting EV s visit to Ilmena in the framework of the research training grop (Gradiertenkolleg) Lorentz force velocimetry and Lorentz force eddy crrent testing.

13 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 13 (m x m y m z ) = ( ) (m x m y m z ) = ( ) (a) (b) (m x m y m z ) = ( ) (c) (d) (m x m y m z ) = ( ) (m x m y m z ) = ( ) (e) (f) Figre. Eddy crrents and secondary magnetic field for the case m y = 0, i.e. when the magnetic dipole is located in the plane spanned by the velocity vector and the direction normal to the srface of the plate. The points where the secondary magnetic field lines emanate are selected manally inside and nearby the eddy-crrent loops in order to gide the eye.

14 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 14 (m x m y m z ) = ( ) (m x m y m z ) = ( ) (a) (b) (m x m y m z ) = ( ) (c) (d) (m x m y m z ) = ( ) (m x m y m z ) = ( ) (e) (f) Figre 3. Eddy crrents and secondary magnetic field for the case m z = 0, i.e. when the magnetic dipole is located in a plane parallel to the srface of the moving plate. The points where the secondary magnetic field lines emanate are selected manally inside and nearby the eddy-crrent loops in order to gide the eye.

15 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 15 Appendix A. Eddy crrents expressed throgh a stream fnction By inserting φ(x,y,z) from (8) into (5), we find the components of electric crrent as: [ } 1 z j x = σ{ x φ} = µσ { x y n tanh + x φ, (A.1) [ } 1 z j y = σ{ y φ B z } = µσ { x x n tanh y φ, (A.) j z = 0, (A.3) To simplify the above expression for j y, we sed B z from (3) and the fact that [ [ 1 z 1 z xx tanh = ( yy + zz ) tanh. (A.4) The obtained eddy crrents are prely horizontal. De to their two-dimensionality they can be expressed as j = (ψe z ) where the stream fnction ψ is [ ψ(x,y,z) = µσ { x n tanh ψ (x,y) = lim z 0 x n } +ψ (x,y), (A.5) = sgn( z) x n lnr. 1 z [ tanh 1 z The fnction ψ (x,y) is a conjectre of the fnction φ (x,y) and it is needed to nllify eddy crrents at the infinite distance from the dipole. When the dipole is oriented arbitrarily, the general expression for ψ can be given explicitly with the aid of the tensor {ψ (n) } = (ψ (x),ψ (y),ψ (z) ) as follows: ψ = µσ ( m {ψ (n) } T) ( ), {ψ (n) } = F 0 +F f F g F 1 f 1, (A.6) where and f 1 = x r, f k = f 1 f k 1 g 1 g k 1, (A.7) g 1 = y r, g k = g 1 f k 1 +f 1 f k 1 ; (A.8) F 0 = z 3, F 1 = r 3, F = z + sgn( z) [ 1 z. (A.9) 3 r The eddy crrents j = µσ ( m {j (n) } T) can be represented with the aid of the tensor {j (n) } which is: C 1 g 1 +C 3 g 3 C 1 f 1 C 3 f 3 C g {j (n) } = 3C 1 f 1 C 3 f 3 C 1 g 1 C 3 g 3 C 0 C f, (A.10) 0 0 0

16 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 16 where C 0 = 3r 1 5 3, C 1 = sgn( z) 3r z, C 4 5 = 3r 5, (A.11) { [ z 3r C 3 = sgn( z) r [ } z r r 3 1. (A.1) ThecoefficientsC k tendtozerowhen z goestoinfinity. Tofindtheeddycrrents nder the dipole, r = 0, one has to make a series expansion arond r = 0. It gives: C 0 r 0 = 1 } { 1+ 3r +O ( r 4), (A.13) z 3 z C 1 r 0 = sgn( z)3r } {1 5r +O ( r 5), (A.14) 4 z 4 z } C r 0 = { 1+ 3r 5r +O ( r 5), (A.15) z 5 z } C 3 r 0 = { sgn( z)5r3 +O ( r 5) (A.16) 8 z 6 Ths at r = 0, C 0 is not zero. It is negative, hence, nder the dipole, the direction of the eddy crrent is opposite to the y-axis, when the dipole is oriented along the z-direction in agreement with the right-hand rle. These series expansions reslts will also be sefl afterwards in the analysis of the secondary magnetic field and forces acting pon the dipole. It is important to notice that the fnctions f k and g k above depend on x and y only and do not depend on z. This will be be sefl below when we compte secondary magnetic field by integrating over the plate thickness z d ẑ z. Moreover, these fnctions f k and g k are specially determined is sch a way, that if the origin of the Cartesian coordinates is selected on the line passing throgh the magnetic dipole perpendiclar totheplate, x 0 = y 0 = 0, theninthecylindrical coordinates, tanϕ = y/x, f k = coskϕ and g k = sinkϕ, i.e. the formlae above represent Forier decomposition. There mst be no anglar dependence in a series expansion arond r = 0 precisely at the point r = 0, therefore, C 0 is the only non vanishing coefficient in (A.13). Appendix B. Primitive fnction for a thin layer Here, we prove that the soltion of (15) for an infinitely thin layer is given by: { } A 0 (ẑ) = w tanh 1 ˆ z(ẑ) +lnr, (B.1) ˆ(ẑ) where ˆ(ẑ) = r + ˆ z(ẑ), ˆ z(ẑ) = ẑ z 0 + ẑ z. Fnction ˆ z(ẑ) can be nderstood from the following geometric interpretation. It is the sm of two terms, the first term, ẑ z 0, is a distance between z 0 and ẑ, and the second term ẑ z is a distance between ẑ and z. Hence, the layer located at z = ẑ

17 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 17 sbdivides the whole space into two half-spaces, z > ẑ and z < ẑ, therefore the fnction ˆ z(ẑ) in straight brackets takes the following form: { z +z 0 ẑ if (ẑ z)(ẑ z 0 ) > 0; ˆ z(ẑ) = (B.) z z 0 = z if (ẑ z)(ẑ z 0 ) < 0. The first line is applicable when z and z 0 are both in the same half-space, hence, ˆ z(ẑ) = z + z 0 ẑ in this case; the second line is applicable when z and z 0 are in different half-spaces, hence, ˆ z(ẑ) = z z0 in this case. Then, it is important that ˆ z(ẑ) shows a discontinity in the first derivative at z = ẑ, i. e.: d dz ˆ z(ẑ) = { +1 if z > ẑ, 1 if z < ẑ, hence 1 d dz ˆ z(ẑ) = δ(z ẑ). Now it is easy to ascertain by means of straightforward calclations that A 0 = 0 everywhere except for the point at z = ẑ while at this point the formla (15) is valid becase: A 0 x + A 0 = w ˆ z(ẑ), y { ˆ 3 A 0 = w ˆ z(ẑ) [ ˆ z(ẑ) } d + 1ˆ d z ˆ 3 dz dz ˆ z(ẑ), A 0 = A 0 x + A 0 y + A 0 z = wδ(z ẑ) 1, where the eqality ˆ = at z = ẑ is taken into accont. Appendix C. Tensor for the secondary field of a thin layer As follows from the mapping (18-0) between indced field and eddy crrents, and the tensor (A.10) for the eddy crrents, the tensor of the indced field can be written down as: {b (n) } = 3f 1 P 1 f 3 P 3 g 1 P 1 g 3 P 3 f P P 0 g 1 P 1 g 3 P 3 f 3 P 3 f 1 P 1 g P P 0 f P g P 4f 1 P 1 (C.1) where the coefficients P k wp 0,k for an infinitely thin layer are obtained from the mapping P 0,k sgn(ˆ z)c k /, k = 1,...,3, by keeping in mind that ˆ and ˆ z are to be sed instead of and z. Specifically, P 0,0 = sgn(ˆ z) { 3r 4ˆ 1 ˆ 3 P 0, = sgn(ˆ z) 3r 4ˆ 5, P 0,3 = ˆ z r 3 }, P 0,1 = 3r ˆ z, 8ˆ 5 [ 3r 4ˆ + 1 r + 1r 3 [ ˆ z ˆ 1. (C.) The series expansion (A.13)-(A.16) is also applicable for (C.). One can see that at r = 0 P 0,0 is not vanishing only, therefore, as follows from (C.1), the nonzero tensor components at r = 0 are b x (z) = b (x) z only.

18 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 18 Appendix D. Secondary magnetic field tensor for a half-space and for a layer with finite thickness For a half-space and for a layer with finite thickness, the secondary field tensor is of the same strctre as for an infinitely thin layer, Eq. (C.1), with the following coefficients for the half space: P,0 = ˆ z 8ˆ 3, P,1 = r 16ˆ 3, P, = sgn(ˆ z) { [ ˆ z 8ˆ + 1 ˆ z 3 4r P,3 = r 16ˆ ˆr + sgn(ˆ z)ˆ r 3 where ˆ z and ˆ shold be compted at z. For the layer with finite thickness, the tensor coefficients are: P k (z,z d ) = P,k (z ) P,k (z d ), k = 0...4; } ˆ 1 [ ˆ z ˆ 1, (D.1), (D.) where P,k (z ) and P,k (z d ) are compted according to (D.1) with ˆ z and ˆ taken at points z and z d. Appendix E. Force field tensor The force field tensor, {f (n) i } = i {b (n) } can be written down in the following way: 3Q 0 +f Q +f 4 Q 4 g Q +g 4 Q 4 3f 1 Q 1 +f 3 Q 3 {f x (n) } = g Q +g 4 Q 4 Q 0 f 4 Q 4 g 1 Q 1 +g 3 Q 3,(E.1) 3f 1 Q 1 f 3 Q 3 g 1 Q 1 g 3 Q 3 4Q 0 +f Q {f y (n) } = {f z (n) } = g Q +g 4 Q 4 Q 0 f 4 Q 4 g 1 Q 1 +g 3 Q 3 Q 0 f 4 Q 4 g Q g 4 Q 4 f 1 Q 1 f 3 Q 3 g 1 Q 1 g 3 Q 3 f 3 Q 3 f 1 Q 1 g Q 3f 1 Q 1 f 3 Q 3 g 1 Q 1 g 3 Q 3 4Q 0 +f Q g 1 Q 1 g 3 Q 3 f 3 Q 3 f 1 Q 1 g Q 4Q 0 f Q g Q 4f 1 Q 1 The coefficients Q k wq 0,k for a infinitely thin layer are: Q 0,0 = 3 ˆ z { } 5r 8ˆ 5 ˆ 1, Q 0,1 = sgn(ˆ z) 3r { } 5r ˆ 5 4ˆ 1, Q 0, = 15r ˆ z, Q 0,3 = sgn(ˆ z) 15r3 8ˆ 7 8ˆ 7, Q 0,4 = 3 ˆ z [ [ 5r ˆ 3 8ˆ ˆ r ˆ z r 4 ˆ 1, (E.). (E.3). (E.4) For a half-space, the coefficients Q k Q,k are obtained by sing the same method as done before for the coefficients P,k. This involves an integration over the thickness of

19 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 19 the plate, Q,k = z Q 0,k(ẑ) dẑ, and the reslts of the integration are: Q,0 = 1 { } 3r 16ˆ 3 ˆ 1, (E.5) Q,1 = sgn(ˆ z) 3r ˆ z 16ˆ, Q 5, = 3r { 16ˆ 5, ˆ z Q,3 = sgn(ˆ z) [1+ 3r + 1 4rˆ 3 4ˆ r 3 Q,4 = 3 [1+ r 4ˆr 4ˆ + r4 + 3ˆ 8ˆ 4 r 4 For a plate with finite thickness we have: [ ˆ z [ ˆ z ˆ 1 Q k = Q,k (z ) Q,k (z d ), k = ˆ 1 } (E.6), (E.7). (E.8) Appendix F. Non-zero tensor coefficients at r = 0 and z = z 0 The coefficients for a thin layer of thickness w z 0 ẑ located at ẑ are: P 0 r=0,z=z0 w P 0,0 r=0,z=z0 = w sgn(z 0 ẑ) 16 z 0 ẑ, (F.1) 3 3 Q 0 r=0,z=z0 w Q 0,0 r=0,z=z0 = w 18 z 0 ẑ 4. (F.) The same coefficients for a half-space are obtained by integration over ẑ from z p to infinity 1 P 0 r=0,z=z0 P,0 r=0,z=z0 = 3 z 0 z, (F.3) 1 Q 0 r=0,z=z0 Q,0 r=0,z=z0 = 18 z 0 z 3. (F.4) Finally, for a plate with finite thickness, z d ẑ z : P 0 r=0,z=z0 = 1 { } 1 3 z 0 z d 1, (F.5) z 0 z Q 0 r=0,z=z0 = 1 { } 1 18 z 0 z d 1. (F.6) 3 z 0 z 3 eferences [1 I. Bcenieks. Electromagnetic indction flowmeter on permanent magnets. In Proceedings of the 8th Int. Pamir Conference on Fndamental and Applied MHD Vol.1, pages , amatelle, France, 00. [ I. Bcenieks. Modelling of rotary indctive electromagnetic flowmeter for liqid metal flow control. In Proceedings of the 5th Int. Symposim on Magnetic Sspension Technology, pages 04 08, Dresden, Germany, 005. [3 P. A. Davidson. An introdction to Magnetohydrodynamics. Cambridge University Press, 001. [4 D. J. Griffiths. Introdction to Electrodynamics, 3rd ed. Englewood Cliffs,, NJ: Prentice Hall, 1989.

20 Interaction of a magnetic dipole with a slowly moving electrically condcting plate 0 [5 J.D. Jackson. Classical Electrodynamics, Third Edition. Wiley, New York, NY, U.S.A., [6 J. Priede. Single-magnet rotary flowmeter for liqid metals. J. Appl. Phys, sbmitted, [7 P. H. oberts. An introdction to Magnetohydrodynamics. Longmans, Green, New York, [8 J. A. Shercliff. The theory of electromagnetic flow-measrement. Cambridge University Press, 196. [9 A. Thess, E. V. Votyakov, Knaepen B., and Zikanov O. Theory of the Lorentz force flowmeter. New J. Phys., 9:99, 007. [10 A. Thess, E. V. Votyakov, and Y. Kolesnikov. Lorentz Force Velocimetry. Phys. ev. Lett., 96:164501, 006.