Electrically Induced Instabilities of Liquid Metal Free Surfaces

Similar documents
Stability of Liquid Metal Interface Affected by a High-Frequency Magnetic Field

51. IWK Internationales Wissenschaftliches Kolloquium International Scientific Colloquium

Experimental and numerical investigation on particle-induced liquid metal flow using Lorentz force velocimetry

Numerical Study of Magnet Systems for Lorentz Force Velocimetry in Electrically Low Conducting Fluids

Suppression of Temperature Fluctuations by Rotating Magnetic Field in a Large Scale Rayleigh-Bénard Cell

AC & DC Magnetic Levitation and Semi-Levitation Modelling

Investigations and Experiments of Sophisticated Magnet Systems for a first Lorentz Force Velocimeter for Electrolytes

Free Surface Controlled by Magnetic Fields

APPLICATION OF A MULTI-DEGREE-OF-FREEDOM SENSOR IN LOCAL LORENTZ FORCE VELOCIMETRY USING A SMALL-SIZE PERMANENT MAGNET SYSTEM

RESAGK, Christian 1 ; DIETHOLD, Christian 2 ; FRÖHLICH, WERNER, Michael 3 ;

Electromagnetic flow rate measurement in molten tin circulating in a closed-loop test system

Lorentz force velocimetry using small-size permanent magnet systems and a multi-degree-of-freedom force/torque sensor

GCSE PHYSICS REVISION LIST

ENGI Multiple Integration Page 8-01

LIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE

MAGNETOHYDRODYNAMICS

Numerical Simulation of Lorentz Force Enhanced Flow Patterns within Glass Melts

Mechanism of atomization in a two-layer Couette flow

Care should be taken to give an appropriate number of significant figures in the final answers to calculations.

THE RIGA DYNAMO EXPERIMENTS

Graduate Preliminary Examination, Thursday, January 6, Part I 2

FORCE MEASUREMENTS BY STRAIN GAUGE SENSORS AS PART OF TIME-OF-FLIGHT FLOW RATE CONTROL

Eddy Current Interaction of a Magnetic Dipole With a Translating Solid Bar

9th European LS-DYNA Conference 2013

St Olave s Physics Department. Year 11 Mock Revision Checklist

INFLUENCE OF GEOMETRICAL ASPECT RATIO ON THE OSCILLATORY MARANGONI CONVECTION IN LIQUID BRIDGES R. MONTI, R. SAVINO, M. LAPPA *

Higher Orders Instability of a Hollow Jet Endowed with Surface Tension

PHYSICS ADVANCED HIGHER. Unit 3 Electromagnetism Homework

Turbulence Model Affect on Heat Exchange Characteristics Through the Beam Window for European Spallation Source

Interfacial waves in steady and oscillatory, two-layer Couette flows

PHYS 432 Physics of Fluids: Instabilities

arxiv:physics/ v1 14 Nov 2006

arxiv:patt-sol/ v1 7 Mar 1997

On the effects of Non-Newtonian fluids above the ribbing instability

COPPER FOR BUSBARS CHAPTER 4: SHORT-CIRCUIT EFFECTS

1 cm b. 4.4 mm c. 2.2 cm d. 4.4 cm v

1.060 Engineering Mechanics II Spring Problem Set 1

Chapter 5. Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice

Cornu s method for the determination of elastic constants of a Perspex beam Surjeet Singh

Liquid Feed Injection in a High Density Riser

Contents. I Introduction 1. Preface. xiii

EFFECT OF WALL JET ON OSCILLATION MODE OF IMPINGING JET

Table of Contents. Foreword... xiii. Preface... xv

Predicting Breakup Characteristics of Liquid Jets Disturbed by Practical Piezoelectric Devices

Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions

Fig. 8.1 illustrates the three measurements. air medium A. ray 1. air medium A. ray 2. air medium A. ray 3. Fig For Examiner s Use

CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer

Answers to questions in each section should be tied together and handed in separately.

Sample paper 1. Question 1. What is the dimensional formula of torque? A. MLT -2 B. MT -2 C. ML 2 T -2 D. MLT -1 E. ML 3 T -2.

Exploration of the Feasibility of Polar Drive on the LMJ. Lindsay M. Mitchel. Spencerport High School. Spencerport, New York

Scattering of internal gravity waves

Exercise: concepts from chapter 10

Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection

Fluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow

BRITISH PHYSICS OLYMPIAD BPhO Round 1 Section 2 18 th November 2016

Simple examples of MHD equilibria

Gesellschaft für Schwerionenforschung mbh (GSI), Planckstrasse 1, D Darmstadt, Germany

Numerical modeling of magnetic induction and heating in injection molding tools

Simultaneous Induction Heating and Electromagnetic Stirring of a Molten Glass Bath

Physics 42 Exam 3 Spring 2016 Name: M T W

PHYA2. (JUN15PHYA201) WMP/Jun15/PHYA2/E4. General Certificate of Education Advanced Subsidiary Examination June Mechanics, Materials and Waves

VERIFICATION TEST OF A THREE-SURFACE-MULTI-LAYERED CHANNEL TO REDUCE MHD PRESSURE DROP

MAGNETOHYDRODYNAMICS Vol. 51 (2015), No. 3, pp

Experiments on capillary-gravity waves of solitary type on deep water

Electromagnetic interaction between a rising spherical particle in a conducting liquid and a. localized magnetic field for

Magnetic Levitation of a Large Mass of Liquid Metal

Intermittency in spiral Poiseuille flow

TORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR

Modelling of Rotating Permanent Magnet Induced Liquid Metal Stirring

# Ans Workings / Remarks

Temperature, Pressure and Velocity measurements on the Ranque-Hilsch Vortex Tube

True/False. Circle the correct answer. (1pt each, 7pts total) 3. Radiation doesn t occur in materials that are transparent such as gases.

Exam 4 (Final) Solutions

IIT JEE PAPER The question paper consists of 3 parts (Chemistry, Mathematics and Physics). Each part has 4 sections.

Equilibrium Evolution in the ZaP Flow Z-Pinch

Bridging by Light. Robert Schroll Wendy Zhang. University of Chicago. 9/17/2004 Brown Bag Talk p. 1

Self-Excited Vibration in Hydraulic Ball Check Valve

The aims of this experiment were to obtain values for Young s modulus and Poisson s ratio for

AC : A STUDENT PROJECT ON RAYLEIGH-BENARD CONVECTION

TEST 2. This test is on the final sections of this session's syllabus and. should be attempted by all students.

PH 102 Exam II. 1. Solve 8 problems out of the 12 below. All problems have equal weight.

Vibration of submillimeter-size supported droplets

Magnetic Levitation and Noncoalescence of Liquid Helium

10.52 Mechanics of Fluids Spring 2006 Problem Set 3

Density Field Measurement by Digital Laser Speckle Photography

APRIL 2015 EXAMINATION version A PHY 132H1S Duration - 2 hours

Welcome A STUDY ON UNIFORMITY OF A MAGNET. Excerpt from the Proceedings of the 2014 COMSOL Conference in Bangalore

Induction thermography on CFRP and the role of anisotropy

High-performance Forehearth Coloring using Lorentz Forces

1. In Young s double slit experiment, when the illumination is white light, the higherorder fringes are in color.

Numerical 2D Modelling of Turbulent Melt Flow in CZ System with AC Magnetic Fields

Dynamics of Air-Blown Dimples

5. You may use the following values of physical constants wherever necessary. Class XII Physics (042) Sample Question Paper

Impact of the forces due to CLIQ discharges on the MQXF Beam Screen. Marco Morrone, Cedric Garion TE-VSC-DLM

Calculus Relationships in AP Physics C: Electricity and Magnetism

Chapter 9. Reflection, Refraction and Polarization

Fig. 2.1 I =... A [2] Suggest why it would be impossible for overhead cables carrying an alternating current to float in the Earth s magnetic field.

ASTRA simulations of the slice longitudinal momentum spread along the beamline for PITZ

Last Name: First Name:. UNIVERSITY OF TORONTO Faculty of Arts and Science. APRIL 2013 EXAMINATION version 2 PHY 205H1S

Measurement of the Liquid Film Thickness in. Micro Tube Slug Flow

Transcription:

International Scientific Colloquium Modelling for Material Processing Riga, June 8-9, 2006 Electrically Induced Instabilities of Liquid Metal Free Surfaces D. Schulze, Ch. Karcher, V. Kocourek, J.U. Mohring Abstract The dynamic behaviour of liquid metal submitted to a high-frequency magnetic field is investigated experimentally. We study both a liquid metal drop placed on a curved glass plate and a liquid metal surface confined in a thin annulus gap. We observe several drop instabilities depending on drop volume and magnitude and frequency of the applied field as well as various deformations of the free surface depending on frequency and field magnitude. Introduction High-frequency magnetic fields can be effectively used to shape and to control free surfaces of electrically conducting liquids like liquid metals, see Moreau [1] and Davidson [2] for an overview. Fautrelle et al. [3] study the dynamic behaviour of liquid metal drops in lowfrequency magnetic fields. They show that for low values of the magnetic field, the wave pattern is axisymmetric and the oscillation frequency is identical to that of the electromagnetic forces. For larger magnetic field amplitudes, the observed waves are structured, but no longer axisymmetric. For even larger magnetic field strengths, a peculiar chaotic pattern appears, characterized by fingers and cavity formation. They observe, that under different conditions (electromagnetic fields), the amplitude of 2D azimuthal modes increases upon increasing the frequency of the field. In a theoretical study, Fautrelle and Sneyd [4] find that the instability of a liquid metal free surface submitted to a parallel alternating magnetic field sets in when the magnitude of the magnetic field exceeds a certain critical value B c. Within this context we study experimentally the effects of a high-frequency magnetic field on the dynamical behaviour of liquid metal surfaces. The results of the experimental observations are in good agreement with predictions of an analytical model. 1. Experimental Methods In both experimental configurations the electromagnetic field is generated by an inductor supplied with alternating electrical current I cos(2πt). We increase the feeding current within the range 0 < I < 350 A and vary its frequency within the range 20 khz < f < 50 khz. As a test liquid we use the alloy Indium-Gallium-Tin called Galinstan. This metal shows a melting temperature of 19 C allowing precise measurements at room temperature. To avoid the oxidation of the free surface, the surface is fully covered by diluted hydrochloric acid or diluted sodium hydroxide. The liquid metal shows the following properties at room temperature: electrical conductivity σ = 3.46 10 6 S/m, density ρ = 6440 kg/m 3, surface tension σ S = 0.718 N/m, and a static contact angle of Θ = 168. The eddy currents induced in the liquid metal lead to an overall rise in temperature, thereby changing material properties like electrical conductivity and surface tension. These unwelcome Joule heat losses can only partly be removed through intensive water-cooling, 50

limiting each experimental run to only a few seconds. All the measurements take place in the range between 30 C and 60 C. The used generator allows the frequency to be adjusted only in discrete steps of some khz. The inductor current is measured by using a shunt and a spectrum analyzer device. During the experiments we observe the free surfaces shapes from above via a highspeed camera taking up to 307 frames per second. The maximum resolution is 1024 x 1024 pixels, with a colour depth of 8 bit. Due to the finite camera memory of 512 MB the recording length is limited to t = 6.65 s. b = 33 Hydrochlorid acid r = 48 2 R High speed camera Galinstan volume V with and without magnetic field Water - cooled bottom 2 R 0 Lamp Inductor with N = 10 Fig. 1. Experimental set-up for the drop 2. Drop Experiment In this model experiments, Kocourek et al. [5] study the static deformation of a circular liquid metal drop due to an electromagnetic pressure that is generated by a highfrequency magnetic field. Applications of this investigations include electron beam evaporation, laser welding and cold crucible technologies. The test facility for the liquid metal drop experiment is show in Fig. 1. The used experimental techniques allow us to measure static drop deformations R/R 0 as well as the critical current I c for different frequencies and drop volumes V. For relatively small value of the inductor current I, the highfrequency field results in a stable symmetrical drop squeezing. Initially flat drops (due to gravity) show a nearly hemispherical shape when an inductor current of a frequency of 20 Fig. 2. Static drop deformations (squeezing) at 20 khz khz is switched on. Fig. 2 51

summarizes our experimental results on drop squeezing. It shows the normalized drop radius R/R 0 as function of the drop volume V, where R 0 denotes the drop radius without magnetic field. At first, small drops are inherently stabilized by surface tension. Upon increasing the inductor current beyond a certain value, the static shape of the drop becomes unstable. We observe the onset of drop oscillations (waves in the r ϕ plane). The experimental findings are illustrated in Fig. 3. Upon increasing the Volume V, we observe that the mode number m of the exited axisymmetric waves increase likewise. For V = 1 ml we observe a cigar-type instability with mode number m = 2 (see Fig. 3 left) while for V = 5 ml the first unstable mode shows a triangular symmetry with m= 3 (see Fig. 3 right). As mentioned before, for V < 0.5 ml we are unable to detect any instability. In fact, we cannot rule out that this Fig.3. Drop deformations at current I = 100 A, frequency f = 20 khz at mode 2 and 3 state becomes unstable at some higher value of the current I that lies beyond our experimental feasibility. Fig. 4 shows the corresponding stability diagram. Here the critical current I c is plotted against the volume V. Beneath the curves, i.e. I < I c, we observe stable circular drops squeezed by the induced electromagnetic pressure. Above the curves, i.e. I > I c, these static states are unstable and waves of a particular mode number m and oscillation frequencies f c are observed. Upon increasing the current the first unstable mode is characterized by m = 2. This mode corresponds to drop oscillations of cigar-type shape as shown in Fig. 3 (left). This mode is present in the range 0.5 ml < V < 5 ml. At small drop volume, the critical current I c decreases rapidly with increasing V but at higher values of V it starts to increase. In this range a transition to mode number m = 3 with triangular Fig. 4. Stability diagram for 20 khz at mode m = 2 and 3 symmetry takes place, cf. Fig. 4 (right). This scenario is repeated at a volume of about V = 10 ml where the transition from m = 3 to m = 4 occurs, i.e. the transition from modes with triangular symmetry to modes with square symmetry. 52

3. Annulus Gap In the second model experiment, Mohring et al. [6] study the electromagnetically induced deformation of a free surface confined in an annual gap. The motivation for these investigations comes from electromagnetic slit-sealing. Inductor with N = 2 * 5 b = 33 r = 48 High speed camera L = 12...22 Dilated sodium hydroxide Ring-like-lamp Free surface Mirror Galinstan H = 20...33 a = 1...6 Water-cooled bottom (copper) D = 150 Fig. 5. Investigated arrangement for the annual gap Perspex All length in mm! The apparatus used is shown schematically in Fig. 5. The characteristic dimensions are the annulus width a, the outer annulus diameter D, the fluid height H and the coupling gap L. The set-up consists of a water-cooled bottom made of copper. The inner and outer walls of the annulus are formed by circular rings of transparent perspex of thickness 10 mm. The width of the annulus is varied by using inner cylinders of various diameters. The gap of the annulus is filled with Galinstan. The coupling gap L can be adjusted by a special device sustaining the bottom. The deformations of liquid metal surface are observed from above by a high-speed camera using a special mirror in the centre. During the measurements we vary the electrical parameters like the frequency f and the inductor current I. At a first value I of the inductor current the initially undisturbed flat surface becomes unstable to small-amplitude waves characterized by a certain wave number that depends strongly on inductor frequency. The amplitude of the waves increases upon increasing I. Upon another slight increase of the inductor current, these surface oscillations are superseded by high-amplitude static deformation of sine wave shape. This shape, shown in Fig. 6 (left), represents a frozen state, i.e. the position of the maximum is nearly timeindependent. The observed wavelength λ is dependent on both geometrical and electrical parameters. The static shape is the result of minimizing the total energy of the free surface consisting of: (i) the electromagnetic energy, (ii) the potential energy, and (iii) the energy of surface tension. 53

At slightly higher values of the critical inductor current I c, the state of static surface deformation switches over into a further kind of instability, the electromagnetic pinch. At this critical value, the electromagnetic pressure rises rapidly in the troughs and exceeds the stabilizing effects of hydrostatic pressure and surface tension pressure. A thin channel growing downwards is formed in the trough of the wave. Fig. 6 (right) shows such a pinch channel in the liquid metal. A drop-like cavity filled with sodium hydroxide is formed at the Sodium hydroxide 10 mm Galinstan Pinch channel Fig. 6. Free surface for a = 3 mm, H = 24 mm, L = 17 mm and f = 21,9 khz with static deformation at I = 127 A (left) and with pinch effect at I = 130 A end of the channel. Due to the proximity effect, the electromagnetic forces on the channel walls act repulsively to each other and thus stabilize the thin channel for a short moment. However, the cavity is very unstable. We repeat the experiments to determine the dependence of the critical inductor current I c on the geometrical parameters a and L, and on the current frequency f. Fig. 7 shows the critical current as a function of the frequency at various gap widths for a constant fluid height ( H = 24 mm ) and coupling gap ( L = 17 mm ). The measured points are interpolated by the dashed lines using a function according to I c f f. ( ) α We obtain higher values of α for lower gap widths. The dotted line is a theoretical curve which was calculated by Karcher and Mohring [7]. They perform a linear stability analysis using both the skin depth and the Hele-Shaw approximation. For the onset of a longwavelength instability this 200 model predicts an 180 exponent of α = 1/4. The 160 measured exponent 140 α ( a = 1mm) = 0.2323 is 120 100 80 only a little smaller then this theoretical value. The critical current for the onset of the pinch is also investigated by varying the geometrical 60 a = 1 mm parameters L, a, and H. It a = 3 mm is found that these a = 6 mm parameters only show an Theory electromagnetic effect as 25 30 35 40 45 both the hydrodynamic 50 and the surface tension I c in A f in khz Fig. 7. Critical inductor current in dependence on frequency f width and frequency f for H = 24 mm and L = 17 mm pressure unchanged. remain 54

Conclusions We have studied experimentally the dynamics of both liquid metal drops and a free liquid metal surfaces in an annulus affected by a high-frequency electromagnetic field. In the drop experiment, at a small magnetic field, we observe static symmetric drop squeezing. However, after increasing the field (I > I c ) these static states become unstable against waves of a particular mode number m. In the experiment with the annulus, upon increasing the inductor current we observe different kinds of instabilities. At low values of the inductor current we observe oscillating capillary surface waves with a wavelength corresponding to the electromagnetic skin depth. At higher values of the current, we find a static surface deformation with a wavelength depending on geometrical and electrical parameters. Eventually, the onset of an electromagnetic pinch is observed. We find a good qualitative and quantitative agreement of the data with predictions of an analytical stability analysis. Acknowledgements The This work was sponsored by the Deutsche Forschungsgemeinschaft (DFG) under grant FOR 421/1-2 Ka-A3 within the frame of the Forschergruppe Magnetofluiddynamik. References [1] Moreau, R.: Magnetohydrodynamisc. Kluver, Dortrecht, 1990. [2] Davidson, P.A.: An Introduction to Magnetohydrodynamisc. Cambridge University Press, 2001. [3] Fautrelle, Y., Etay, J., Daugan, S.: Free-surface horizontal waves generated by low-frequency alternating magnetic analysis. J. Fluid Mech. 527, 285-301, 2005. [4] Fautrelle, Y., Sneyd, A.D.: Instability of a plane conducting free surface submitted to alternating magnetic field. J. Fluid Mech. 375, 65 83, 1998. [5] Kocourek, V., Karcher, Ch., Conrath, M., Schulze, D.: Stability of liquid metal drops affected by a highfrequency magnetic field. submitted to Phys. Rev. E (2006). [6] Mohring, J.U., Karcher, Ch., Schulze, D.: Dynamic behavior of a liquid metal interface under the influence of a high-frequency magnetic field. Phys. Rev. E 71, 301 306, 2005. [7] Karcher, Ch., Mohring, J.U.: Stability of a liquid metal interface affected by a high-frequency magnetic field. Magnetohydrodynamics 39, 267 276, 2003. Authors Prof. Dr.-Ing. Schulze, Dietmar Dipl.-Ing. Kocourek, Vaclav Dipl.-Ing. Mohring, Jens-Uwe Faculty of Electrical Engineering and Information Technology Dr.-Ing. Karcher, Christian Faculty of Mechanical Engineering Technische Universität Ilmenau Postfach 10 05 65 D-98684 Ilmenau, Germany E-mail: dietmar.schulze@tu-ilmenau.de 55