1 UVOD Leta 1832 je Michael Faraday poizkusil izmeriti hitrost reke Temze blizu mostu Waterloo [1]. Njegova ideja je bila, da je pretok reke izpostavljen zemeljskemu magnetnemu polju, pri čemer se inducira napetost v dveh elektrodah, postavljenih pravokotno skozi pretok. Danes njegov izum elektromagnetnega merilnika pretoka prinaša širok uspeh v kemijski in živilski industriji in ni omejen z vročimi in agresivnimi pretoki, vendar pa je neuspešen pri največjem izzivu merjenja pretoka [2, 3]. To so taline, ki dosežejo visoko temperaturo (železo, aluminij,steklo) ali za lokalne meritve, kjer ovire ne omogočajo dostopa do kanala, saj je metoda z elektromagnetnim merilnikom pretoka zasnovana z velikimi magneti. Problem je tudi vstavljanje elektrod v talino. Te nam omejujejo temperaturno območje uporabe veliko pod točko tališča elektrod. Merjenje pretoka z Lorenzovo silo je že v 50. letih prejšnega stolejta odkril profesor A. Shercliff [2]. Praktična aplikacija te metode se je zaradi tehničnega razvoja v proizvajanju močnih permanentnih magnetih, natančnih principih merjenja sile in procesne simulacije razvila šele v zadnjih letih [4]. Merjenje pretoka z Lorenzovo silo je brezkontaktna metoda, ki deluje na principu interakcije električno prevodne tekočine z magnetnim poljem permanentnega magneta. Deluje na principu vrtinčnih tokov, ki povzročijo, da zaradi tretjega Newtonovega zakona magnet občuti silo, ki je merljiva in odvisna od hitrosti tekočine. Trenutno se razvija v metalurgiji in na drugih področjih. 2 FIZIKALNO OZADJE Brezkontaktno merjenje pretoka električno prevodnih tekočin s pomočjo Lorentzove sile je izvedeno z merjenjem Lorentzove elektromagnetne sile, ki nastane zaradi toka prevodne tekočine pod vplivom izmeničnega magnetnega polja. Faradayev zakon indukcije namreč pravi, da se v kovini ali prevodni tekočini, ki se premika v magnetnem polju, inducirajo vrtinčni tokovi, in sicer na območjih maksimalnega gradienta magnetnega polja. Ti vrtinčni tokovi na podlagi Amperovega zakona ustvarijo lastno magnetno polje. Medsebojna akcija med vrtinčnimi tokovi in primarnim magnetnim poljem ustvari Lorentzovo elektromagnetno silo, ki nasprotuje gibanju oz. toku tekočine. Tudi tukaj velja tretji Newtonov zakon, ki pravi, da akciji vedno nasprotuje reakcija, zato na vir primarnega magnetnega polja Prejet 30. julij, 2015 Odobren 12. avgust, 2015 BREZKONTAKTNO MERJENJE PRETOKA PREVODNIH TEKOČIN S POMOČJO LORENTZOVE SILE 213 deluje nasprotno enaka sila, kot je tista, ki nasprotuje gibanju tekočine. Z meritvijo te reakcijske sile lahko izmerimo hitrost tekočine, saj je ta sila proporcionalna hitrosti tekočine. Ko se torej prevodna tekočina premika skozi magnetno polje, ki je ustvarjeno bodisi s trajnim magnetom bodisi s tokovno vzbujano tuljavo, v njej nastanejo vrtinčni tokovi, ki ustvarijo Lorentzovo silo, ki ima gostoto 𝑓 = 𝑗 × B⃗⃗⃗, in nasprotuje pretoku tekočine. Gostoto te sile lahko aproksimiramo kot 𝐹 ~ 𝜎 𝑣 𝐵2 (1), kjer je 𝜎 električna prevodnost tekočine, 𝑣 njena hitrost in 𝐵 jakost magnetnega polja [2]. Lorentzova sila je torej proporcionalna električni prevodnosti merjene tekočine, relativni hitrosti tekočine glede na magnet in glede na samo jakost električnega polja. Primarno magnetno polje �⃗⃗�(𝑟) lahko ustvarimo s trajnim magnetom ali s primarnim tokom skozi tuljavo 𝐽(𝑟) (glej sliko 1). Gibanje tekočine pod tako vzbujanim magnetnim poljem inducira vrtinčne tokove, ki so prikazani na sliki 2, in jih večkrat poimenujemo sekundarni tokovi ter jih označimo z 𝑗(𝑟). Interakcija med temi vrtinčnimi tokovi in primarnim magnetnim poljem ustvari Lorentzovo elektromagnetno silo v premikajoči se tekočini, ki je enaka �⃗�𝑓 = ∫ 𝑗 × �⃗⃗�𝑑 3𝑟 𝑓 (2) in nasprotuje toku tekočine. Sekundarni tokovi ustvarijo svoje lastno magnetno polje �⃗⃗�(𝑟), sekundarno magnetno polje, ki pa deluje na primarni električni sili in vzbudi Lorentzovo silo �⃗�𝑚 na trajni magnet ali tuljavo �⃗�𝑚 = ∫ 𝐽 × �⃗⃗�𝑑 3𝑟 𝑚 (3). V skladu s tretjim Newtonovim zakonom imata ti dve sili, sila na tekočino �⃗�𝑓 in sila na magnetni sistem �⃗�𝑚, enako velikost, vendar nasprotno smer, torej �⃗�𝑚 = −�⃗�𝑓 (4). Slika 1: Princip merjenja. a) Postavitev tuljave in b) struktura primarnega magnetnega polja za longitudinalni fluksmeter. c) in d) Isto za transverzalni fluksmeter. Povzeto po [5]. Slika 2: Princip merjenja. a) Prevodno tekočino izpostavimo magnetnemu polju. b) Magnetno polje inducira vrtinčne tokove v tekočini. c) Vrtinčni tokovi vzbudijo sekundarno magnetno polje. d) Pojavi se zaviralna sila na tekočino, nasprotna na silo na magnet. Povzeto po [6]. Kvantitativno lahko velikost merjene sile pridobimo tudi tako, da si pomagamo z razmerami v našem sistemu (glej sliko 3). Majhen trajni magnet z dipolnim momentom 𝑚 postavimo na razdaljo 𝐿 nad neskončno tekočino, ki teče z enakomerno hitrostjo 𝑣, ki poteka vzporedno glede na njeno površino. Slika 1: Ploskovna porazdelitev magnetnega polja. a) Primarno magnetno polje B in vrtinčni tokovi J, ustvarjeni z delovanjem magnetnega dipola na enakomerno premikajočo se prevodno tekočino. b) Sekundarno magnetno polje b, vzporedno na vrtinčne tokove J. Povzeto po [7]. 214 KOLAR, KRIVEC, SIMONIČ Slika 2: Graf odvisnosti sile, ki jo povzročijo vrtinčni tokovi od hitrosti tekočine, prikazan s tremi različno močnimi magneti. Povzeto po [5]. Slika 3: Princip merjenja sile z interferometrom. (1) Laser, (2) delilnik žarkov, (3) referenčni odbijalec žarkov (4) merilni odbijalec žarkov, (5) foto-detektor. Povzeto po [8]. Pri analizi takega sistema predpostavimo, da je magnet dipol z dipolnim momentom �⃗⃗⃗� = 𝑚�̂�𝑧, katerega magnetno polje je podano z �⃗⃗�(�⃗⃗�) = 𝜇0 4𝜋 {3 (�⃗⃗⃗⃗�∙�⃗⃗�)�⃗⃗� 𝑅5 − �⃗⃗⃗⃗� 𝑅3 } (5), kjer je �⃗⃗� = 𝑟 − 𝐿�̂�𝑧 in 𝑅 = |�⃗⃗�|. S predpostavko, da je �⃗� = 𝑣�̂�𝑥, za vse 𝑧 < 0, lahko izračunamo vse vrtinčne tokove s pomočjo Ohmovega zakona za premikajoče se električno prevodne tekočine 𝐽 = 𝜎(−∇𝜙 + �⃗� × �⃗⃗�) (6) z robnimi pogoji 𝐽𝑧 = 0, kjer je 𝑧 = 0, in 𝐽𝑧 → 0, ko gre 𝑧 → 1 [3]. Najprej pridobimo skalarni električni potencial kot 𝜙(𝑟) = − 𝜇0𝑣𝑚 4𝜋 𝑥 𝑅3 (7), od koder lahko preračunamo gostoto električnega toka. Ko imamo to znano, lahko uporabimo Biot-Savartov zakon za izračun sekundarnega magnetnega polja �⃗⃗�(𝑟). Sila je nato podana kot �⃗� = (�⃗⃗⃗� ∙ ∇)�⃗⃗� (8), kjer moramo gradient �⃗⃗� izračunati na poziciji dipola. Za dani problem lahko tako analitično pridobimo rešitev 𝐹 = 𝜇0 2𝜎𝑣𝑚2 128𝜋𝐿3 �̂�𝑧 (9), ki nam da približek 𝐹~𝜇0 2𝜎𝑣𝑚2𝐿−3 (10). 3 DIREKTNO MERJENJE SILE Permanentni magnet je postavljen v bližino pretoka, ki ga želimo izmeriti. Lorentzova sila generira vrtinčne tokove. Ti po istoimenskem pravilu inducirajo magnetno polje, ki nasprotuje spremembi magnetnega polja, ki ga je ustvarilo. Pripadajoča sila magnetnega polja vrtinčnih tokov zavira tok tekočine. Po Newtonovem tretjem zakonu mora na permanentni magnet delovati nasprotno enaka sila, ki je šibka, vendar merljiva, in se uporablja za merjenje pretoka elektrolitske snovi. Merilni sistem je sestavljen iz dveh vzporednih vzmeti, na katerih sta pritrjena permanentna magneta, katerih magnetno polje je pravokotno na smer pretoka. Silo, ki jo povzroči Lorentzova sila, lahko merimo s pomočjo odklona, tega pa s torzijo vzmeti. Torzija vzmetnega elementa je izmerjena z merilnikom torzije, žico, kateri se s spreminjanjem dolžine spreminja upornost. Približek sile lahko zapišemo z enačbo (10) 𝐹~𝜇0 2𝜎𝑣𝑚2𝐿−3 [5]. Sila je proporcionalna hitrosti merjene tekočine in tudi prevodnosti, kar vidimo tudi s slike 4. Pri taki meritvi predpostavljamo, da je prevodnost snovi poznana v vsakem trenutku meritve oz. je konstantna, preostali elementi enačbe pa nam določajo občutljivost. To je zelo uporabna ugotovitev, saj vidimo, da je sila odvisna od kvadrata magnetizacije (ali kvadrata električnega toka v primeru elektromagneta) [5]. Občutljivost lahko torej povečamo z močjo magnetnega polja. Iz tega sledi, da bi bila lahko brezkontaktna metoda uporabna tudi za slabo prevodne snovi, kot je npr. steklo, ki je nedosegljivo drugim brezkontaktnim metodam merjenja pretoka. Pri merjenju pretoka slabo prevodnih snovi se odklon meri z interferometrom, kot je prikazano na sliki 5. Magneta drži vzmet iz kvarčnega stekla, ki ima izjemne mehanske lastnosti, kot so npr. nizka elastičnost in nizko temperaturno raztezanje, kar pripelje do linearnosti deformacije vzmeti [8, 9]. Z interferometrom lahko deformacijo zaznamo že na 0.1 nm natančno. BREZKONTAKTNO MERJENJE PRETOKA PREVODNIH TEKOČIN S POMOČJO LORENTZOVE SILE 215 4 ROTACIJSKI MERILNIKI PRETOKA Rotacijski merilnik pretoka, ki deluje na brezkontaktnem principu, je sestavljen iz permanentnih magnetov enakomerno razporejenih po obodu diska, nameščenega na osi v bližini cevi, po kateri se pretaka električno prevodna tekočina. Takšen merilnik imenujemo rotacijski merilnik hitrosti pretoka na osnovi Lorentzove sile (»Lorentz force velocimetry«). Magneti so postavljeni tako, da je njihovo magnetno polje pravokotno na os vrtenja. Magnetno polje permanentnih magnetov povzroči v pretoku vrtinčne tokove, ki po Amperovem zakonu ustvarijo lastno magnetno polje. To sekundarno polje z Lorentzovo silo zaviralno vpliva na pretok tekočine. Po tretjem Newtonovem zakonu, recipročnosti sil, deluje sila z isto absolutno vrednostjo in nasprotnim predznakom na magnet. Elektromagnetna sila na magnete na disku je šibkejša na bolj oddaljeni strani in tako povzroči vrtenje.[6] Zdaj lahko merimo kotno hitrost rotorja, ki je premo sorazmerna s hitrostjo pretoka in obratnosorazmerno z razdaljo med magnetom in cevjo. Druga geometrijska oblika rotacijskega merilnika je en sam cilindričen magnet, katerega magnetno polje je vzporedno z osjo vrtenja.[10] V nasprotju s prvo različico se tukaj magnet vrti zaradi navora, ki ga povzročajo magnetne sile. Ko se ta magnet začne vrteti, nastanejo dodatni vrtinčni tokovi, ki zavirajo rotacijo. Ravnovesno rotacijo dobimo, ko se zaviralni navor izenači z gonilnim. Prednost takega načina je, da je neodvisen od električne prevodnosti tekoče kovine (torej tudi temperature) in moči magneta. Rotacijski merilnik lahko uporabimo tudi kot elektromagnetno črpalko, če na rotor apliciramo zunanji vir rotacije.[11] Slika 4: Prikaz rotacijskega merilnika z uporabo več magnetov. Povzeto po [12]. 5 POSKUS Poskus, katerega namen je dokazati obravnavani princip merjenja, je bil narejen z uporabo pretočne zanke, prikazane na sliki 8. Zanka je sestavljena iz polietilenske cevi, skozi katero se pretaka elektrolitska raztopina natrijevega tetrabonata (Na2B4O7·10H2O) s prevodnostjo ρ = 20 mS/m pri sobni temperaturi. Za pretok tekočine skrbi pretočna črpalka s hitrostjo pretoka 𝑣 = 1,7 m/s. Presek cevi je S = 10 -4 m 2 in je nameščena na minimalni razdalji 1 mm od magneta. Uporabljen je magnet v obliki kocke, stranica je dolga d = 25 mm. Razred magnetizacije magneta je N38, magnetno polje na površini magneta B = 1220 mT. Sila na magnet je merjena s pomočjo teže, in sicer s tehtnico Mettler Toledo hx204, ki meri na 0,1 mg natančno oziroma zazna silo 9,8*10 -7 N. Ker je pretok usmerjen vertikalno ima v neposredni bližini magneta sila samo komponento Z in lahko le-to merimo kot spremembo teže na tehtnici. Največja mogoč sila z uporabljenim sistemom bi bila 30 µN. Ker je tehtnica zelo občutljiva, smo morali najprej izločiti šum okolice tako, da smo naš sistem testirali brez magneta. Pri končnih meritvah smo izmerili povprečno silo 19 µN. Pri uporabi enačbe (9) to pomeni hitrost pretoka 1,05 m/s. Pri meritvah so bila velika odstopanja, vendar je bilo pri vsakem vklopu črpalke očitno, da se generira sila, ki je veliko večja od merilnega pogreška tehtnice. Slika 8: Sestava merilnega okolja. (1) Črpalka, (2) cevi, (3) magnet, (4) tehtnica. 6 SKLEP Opisano brezkontaktno merjenje pretoka prevodnih tekočin obeta velik napredek v metalurgiji, saj imajo takšni merilniki daljšo življenjsko dobo in so natančnejši, hkrati pa cenejši, ker so večinoma izdelani iz standardnih materialov. Z opisanim postopkom, lahko izmerimo tudi hitrost gibanja kovin v trdnem stanju, z njim lahko detektiramo mikrodefekte v materialu, pomaga nam pa tudi pri izdelavi bolj homogenih kovin. V Sloveniji je ta tematika slabo poznana, medtem ko v tujini zbuja čedalje večjo pozornost znanstvene srenje. Prav tako se predvideva, da bo z bolj preciznimi merilniki mogoče merjenje pretoka elektrolitskih tekočin [13]. 216 KOLAR, KRIVEC, SIMONIČ