Computer Modeling in Engineering & Sciences |

DOI: 10.32604/cmes.2021.013597

ARTICLE

A Double-Phase High-Frequency Traveling Magnetic Field Developed for Contactless Stirring of Low-Conducting Liquid Materials

1College of Material Science and Opto-Electronic Technology, University of the Chinese Academy of Science, Beijing, China

2Guangwei (Zhongshan) Intelligent Technology Co., Ltd., Huairou Science City, Beijing, China

3SIMAP/EPM, Universite Grenoble Alpes, Grenoble, France

*Corresponding Author: Xiaodong Wang. Email: xiaodong.wang@ucas.edu.cn

Received: 12 August 2020; Accepted: 21 December 2020

Abstract: The use of low electrically conducting liquids is more and more widespread. This is the case for molten glass, salt or slag processing, ionic liquids used in biotechnology, batteries in energy storage and metallurgy. The present paper deals with the design of a new electromagnetic induction device that can heat and stir low electrically-conducting liquids. It consists of a resistance-capacity-inductance circuit coupled with a low-conducting liquid load. The device is supplied by a unique electric power source delivering a single-phase high frequency electric current. The main working principle of the circuit is based on a double oscillating circuit inductor connected to the solid-state transistor generator. This technique, which yields a set of coupled oscillating circuits, consists of coupling a forced phase and an induced phase, neglecting the influence of the electric parameters of the loading part (i.e., the low-conductivity liquid). It is shown that such an inductor is capable to provide a two-phase AC traveling magnetic field at high frequency. To better understand the working principle, the present work improves a previous existing simplified theory by taking into account a complex electrical equivalent diagram due to the different mutual couplings between the two inductors and the two corresponding induced current sets. A more detailed theoretical model is provided, and the key and sensitive elements are elaborated. Based on this theory, equipment is designed to provide a stirring effect on sodium chloride-salted water at 40 S/m. It is shown that such a device fed by several hundred kiloHertz electric currents is able to mimic a linear motor. A set of optimized operating parameters are proposed to guide the experiment. A pure electromagnetic numerical model is presented. Numerical modelling of the load is performed in order to assess the efficiency of the stirrer with a salt water load. Such a device can generate a significant liquid motion with both controlled flow patterns and adjustable amplitude. Based on the magnetohydrodynamic theory, numerical modeling of the salt water flow generated by the stirrer confirms its feasibility.

Keywords: Alternating magnetic field; electromagnetic stirring; low electrically conducting liquids; ionic liquid

The field of magneto-hydro-energetics (MHE) couples fluid mechanics, electromagnetism and heat transfer. MHE can be applied to various electrically conducting fluid media, for example, liquid metals, electrolytes, molten oxides, ionic liquids [1–3], ionized gases or plasma. Such a field concerns various applications from astrophysics to industrial processes on a wide range of scales.

Many industrial processes involve AC magnetic fields. Induction effects are widely used in the area of processing of conducting liquid materials. The most widespread applications are induction heating and stirring of liquid metals such as steel, aluminium during their elaboration. As far as low conducting liquids are concerned, the applications deal with molten salt, electrolyte or molten glass processing. For example fluoride salts are used as electrolytes in aluminium reduction cells as well as coolant fluids in molten salt reactors [1–11]. Low conducting ionic liquids are also involved in energy storage, electrolyte materials for Li/Na batteries and biotechnology applications [2,3].

The most important advantages of induction devices are the fast heating speed, the good efficiency of the process and the absence of electrical/chemical contact and pollution between the melt and the heating device. The latter advantage is crucial when high-temperature liquids are considered.

The electrical conductivity of low-conducting liquid materials is rather low, as shown in Tab. 1. Their electrical conductivity is mainly related to the mobility of ions in the bath. It may strongly depend on the temperature, e.g., molten oxides or salts, and to a lesser extent on the frequency of the applied electric currents. Efficient coupling between an AC magnetic field and a low-conducting load, i.e., heating and stirring by Lorentz forces, requires the use of high-frequency fields. Optimum coupling conditions occur when the electromagnetic skin depth is on the order of the pool diameter [1]. Specifically, the optimum coupling, both for single-phase and multiphase magnetic fields, is determined by the ratio of the electromagnetic skin depth

The present publication is a continuation of preliminary works performed by Ernst et al. [4]. Its purpose concerns the design of high-frequency multiphase electromagnetic stirrers. Since multiphase generators working at high frequencies are not available, a special design of the ensemble including generator-capacitor-load is needed.

Ernst et al. [4] provided some elements in order to design a two-phase magnetic field device with several hundred kilohertz mimicking a linear motor. Ernst et al. [4] proposed the main working principle of the double oscillating circuit inductor connected to the solid-state transistor generator to form the equipment delivered by CELES Company. This technique, which yields a set of coupled oscillating circuits, consists of coupling a forced phase and an induced phase and neglects the influence of the electric parameters of the loading part (i.e., the low-conductivity liquid).

To better understand the working principle, the present work improves the previous simplified theory by taking into account a complex electrical equivalent diagram due to the different mutual couplings between the two inductors and the two corresponding induced current sets. A more detailed theoretical model is provided, and the key and sensitive elements are elaborated. Based on this theory, equipment is designed to use sodium chloride-salted water at 40 S/m to provide a stirring effect. A pure electromagnetics numerical model is presented. Numerical modelling of the load is performed in order to assess the efficiency of the stirrer with a salt water load.

The present electromagnetic model includes the two inductors coupling to the salt water-induced charge and all the electrical circuit elements connecting the two inductors to the different capacitors, connection inductances, and generator. This model is applied to a frequency sweeping study, which clearly shows the different resonance peaks corresponding to the different working resonance frequencies caught by the generator during the reception tests. Then, it is explained how the transistor generator can switch from one resonance frequency to another resonance frequency and the consequences of the two inductor currents (module difference, phase shift, etc.). Finally, in Section 4 an analytical theory based on the transformer theory is presented that is inspired by the former simplified analytical theory. It is nevertheless more complex since it takes into account the different couplings between the two inductors and the corresponding two sets of induced currents in the salt water, as in the numerical model. The latest numerical model and improved analytical model provide a better understanding of the experimental behavior of the equipment and, in turn, improved mastery of the whole process during stirring operations.

In Section 5, detailed numerical modeling of the experimental set-up is presented. It demonstrates that the traveling magnetic field is able to generate a liquid motion. This is confirmed by some experimental results obtained in a former identical installation, which consists of a salt water cylindrical tank with a diameter of 0.30 m and a height of 0.35 m (Section 6).

2 The Double-Phase Magnetic Traveling Field Stirring Installation

The installation (Fig. 1) includes a 25 kW high frequency transistor solid state converter generator feeding a set of two parallel oscillating circuits, with each consisting of a capacitor battery and an inductor. Each inductor is made of two turns connected in series-opposition, thus providing a “four-poles” magnetic field configuration. The two identical inductors, corresponding to the two electrical phases producing the magnetic traveling field, as in a linear motor, are situated with one inserted into the other, and one inductor can be shifted vertically with respect to the other so that their mutual inductance can be adjusted to reach the optimal electrical working conditions, producing equality between the two inductor current modules and 90

3 Numerical Model of the Double-Phase Inductor with Connection to the Generator and Resonance Frequency Determination

The simplified electrical diagram of the COMSOLTM model is shown in Fig. 2 along the current component values from the experimental tests. The lower forced phase (in red) is connected to the generator output with a

The main result of frequency sweeping is displayed in Fig. 3, which shows the

In Tab. 2, the main results are summarized for the three resonance frequencies, f1, f0 and f2. According to the 3rd and 5th columns, for the two resonance frequencies f1 and f2 onto which the generator locks, the I2/I1 modules ratio is close to 1, and the phase shifts of I2 with respect to I1 are respectively −20

4 Improved Analytical Model of the Coupled Circuits

In the present model, we will take into account all the possible couplings between the two forced and induced coils and the two corresponding sets of induced currents. These two sets of induced currents can be considered as real forced and induced inductors and as virtual coils that are similar to the “images” of the real forced and induced inductors and have the corresponding resistance, self-inductance and mutual inductance relative to the other elements. This is sketched in Fig. 4, with the different couplings corresponding to the different mutual inductances defined between the inductors and induced current virtual coils (for clarity, the induced charge is set out of the inductor coils). The electrical equivalent circuit of the whole system in Fig. 2 is represented in Fig. 5. As the system is supposed to be symmetric between the forced phase and induced phase, implying that homologous parameters are the same, the magnitudes in Figs. 4 and 5 are defined in Tab. 3 (when necessary and for convenience of notation, the indices i and j are used for the forced and induced phases).

Without taking the C1 capacitor into account in the first state, the rest of the electrical diagram in Fig. 5 can be considered as a multi-secondary transformer whose primary is the forced-phase (1) Coil (connected to the generator) and whose secondaries are the induced-phase (2) Coil (connected to capacitor C2) on the one hand and the two induced virtual coils (connected in short circuit and coupling to all the other circuits) on the other hand.

Thus, the classical transformer equations of the system (the C1 capacitor is not yet taken into account) for the forced phase (1), the induced phase (2) and the two induced-current virtual coils are as follows:

Eqs. (4) and (5) can be simplified because the electromagnetic skin depth in the salt water is on the same order of magnitude as the container radius; thus, one can show that

Using the notations above and assuming that

This is the whole impedance of the forced and induced system seen from the forced-phase coil, showing the resistive part and the reactive part. By putting this impedance in parallel with the C1(= C) capacitor, we obtain the total impedance

with:

The induced-to-forced-phases currents complex ratio I2/I1 can also be written with the former notation as follows:

The three resonance frequencies are given by canceling the imaginary part of (7), which requires the solution of a 3rd-degree equation with the x variable. The solutions, after some simplifications due to small terms that can be neglected (this is shown in the numerical application with experimental values), are

The experimental values of the different components in Fig. 5 which are obtained either from either preliminary measurements (mainly by the capacitive discharge method owing to an impedance-meter) or from the previous numerical COMSOLTM model are then applied. The main results are provided in Tab. 6. These values are for the most part very close to the homologous values in Tab. 2, which shows that the analytical model is sufficiently accurate and greatly aids in understanding the process at work. As explained previously, the generator works by locking on either frequency f1 or f2, which correspond to either −25

5 Modeling of the Magnetohydrodynamic Problem

To assess the capability of the stirring device, full magnetohydrodynamic numerical modeling is performed. Fig. 6 provides a sketch of the computation geometry, which comprises a salt water load surrounded by the two-phase coil. The geometry is 2D-axisymmetric. A widely-used Multiphysics finite-element software COMSOLTM (https://www.comsol.com/release/5.6) is used to calculate in a coupled way the magnetic field, Joule dissipation, electromagnetic forces, temperature field and electromagnetically driven flow in the load. The working parameters and physical properties are reported in Tab. 7. The detailed mathematical description of the model (equations, boundary conditions and main numerical data are summarized in Appendix 2. A transient protocol is used to calculate the flow and temperature fields, respectively from rest and room temperature. Indeed, because of the efficient Joule heating of the bath, its temperature rises quite rapidly, and the computations must be stopped before 300–400 s.

Numerical results are presented in Figs. 7–9. Fig. 7 shows the magnetic field distribution in the load and around the coils, whereas Figs. 8 and 9 illustrate the liquid motions. Fluid flow is generated by both electromagnetic forces and buoyancy. In both cases, the flow pattern consists of a single vortex in a meridian plane. Nevertheless, the velocity magnitudes are different. Natural convection due to Joule heating is localized mainly along the vertical walls because of the skin depth effect. Accordingly, buoyancy leads to an upward motion along the vertical walls, whereas the traveling magnetic field may be upward or downward. In Fig. 8, electromagnetic forces are oriented downward. Thus, the electromagnetic forces are opposed to buoyancy, and the velocity magnitude is less than that in the case of an upward traveling magnetic field presented in Fig. 9. Note that the maximum phase shift was used (I2/I1 phase shift imposed to −90

The latter results are confirmed by an analysis of the order of magnitude of the various forces and flow velocity. The typical computed magnetic flux density amplitude along the pool wall is approximately 10 mT. In the present geometry, the pole pitch is 126 mm; hence, the value of the synchronous velocity Vs is

This a quite a low value compared with the liquid metal stirring forces but can be sufficient to stir an electrolyte. Indeed, balancing the fluid inertia and the electromagnetic body force provides an upper bound on the liquid velocity Uem [12], namely,

where

The above estimate is in accordance with the numerical modeling. The numerical results confirmed by the analysis of the order of magnitude indicate that the electromagnetic device may be efficient in acting on a low-conducting liquid.

A former similar stirring installation has been designed in 2004 [4] with the same double phase inductor connected to a triode generator working at about the same frequency. The salt water filled tank was identical to this presented in this paper. Small hollow silver coated glass spheres were immersed in the water in order to show the stirring path thanks to their tracking. To do this a vertical meridian plane in the salt water was illuminated with a vertical slot projected by a slide projector. This vertical meridian plane was then photographed with an exposure time of a few seconds thus showing the glass spheres trajectory whose length is proportional to the velocity. With the triode generator it was possible to almost reach the optimal double phase stirring conditions (mod(I2/I1) and phase (I2/I1) respectively close to 1

The well-organized flow pattern results from additive vertical both upward natural thermo-convection and electromagnetic stirring forces in the peripheral tank area, as shown by the numerical model in Fig. 9. The longest observed spherical particles trail corresponds to a velocity of about 3.5 cm/s which shows the ability of this double phase inductor to produce an efficient stirring.

A high-frequency traveling magnetic field apparatus has been designed. After optimizing the operating parameters, the output power of the generator can be as high as 25 kW. The power of the generator clearly dynamically changes with the impedance of the circuits. The closer the impedance of the external circuit is to the internal resistance of the high-frequency power source, the greater the output power of the high-frequency power source. The external circuit impedance varies with the distance between the two coils. Within this adjustable position range, the greater the distance is, the greater the impedance of the loop, and the greater the difference from the internal resistance of the high-frequency power source is, the higher the output power of the high-frequency power source. The impedance of the external circuit is also related to the conductivity of the salt water solution: The closer the conductivity is to 40 S/m, the closer the impedance of the circuit is to the internal resistance of the high-frequency power source, and the greater the output power of the high-frequency generator. In short, if we achieve a high output power from the high-frequency generator, then the distance between the coils should not be too large, and the conductivity of the salt water should be close to 40 S/m.

The key element in forming a traveling magnetic field is that the phase shift and amplitude ratio of the two-phase current. The phase shift of the currents is related to the distance between the coils. The greater the distance is, the larger the phase difference. It is also related to the conductivity of the load solution. The larger the conductivity of the load solution is, the larger the phase shift. The amplitude ratio of the two-phase current is related to the distance between the coils. The greater the distance is, the smaller the amplitude ratio of the two-phase current. In short, to form a traveling magnetic field, the distance between the coils should be neither too large nor too small. Note that the electrical system may sometimes be subject to instabilities linked to the tuning of the operating frequency. As far as the stirring effect is concerned, numerical modeling shows that the device is efficient enough to generate significant liquid velocities in salt water. Further experimental investigations are underway to validate the present theoretical results.

Funding Statement: This study was supported by the Instrument and Equipment Development Project of the Chinese Academy of Sciences (YJKYYQ20200053), the “Double First-Class” Construction Fund (111800XX62), and the Mechanical Engineering Discipline Construction Fund (111800M000).

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

- Fautrelle, Y., Ernst, R., & Moreau, R. (2009). Magnetohydrodynamics applied to materials processing.
*International Journal of Materials Research*,*100(10)*, 1389-1398. [Google Scholar] [CrossRef] - Watanabe, M., Thomas, M. L., Zhang, S., Ueno, K., & Yasuda, T. (2017). Application of ionic liquids to energy storage and conversion materials and devices.
*Chemical Reviews*,*117(10)*, 7190-7239. [Google Scholar] [CrossRef] - Egorova, K. S., Gordeev, E. G., & Ananikov, V. (2017). Biological activity of ionic liquids and their application in pharmaceutics and medicine.
*Chemical Reviews*,*117(10)*, 7132-7189. [Google Scholar] [CrossRef] - Ernst, R., Perrier, D., Brun, P., & Lacombe, J. (2005). Multiphase electromagnetic stirring of low conducting liquids.
*COMPEL International Journal of Computations and Mathematics in Electrical*,*24(1)*, 334-343. [Google Scholar] [CrossRef] - Khomich, V. Y., Moshkunov, S. I., Rebrov, I. E. (2016). Multiphase high-voltage high-frequency sawtooth wave generator. 2016 International Symposium on Power Electronics, Electrical Drives, Automation and Motion, Anacapri (Italy). DOI 10.1109/SPEEDAM.2016.7525962. [CrossRef]
- Farahat, R., Eissa, M., Megahed, G., Fathy, A., & Abdel-Gawad, S. (2019). Effect of EAF slag temperature and composition on its electrical conductivity.
*ISIJ International*,*59(2)*, 216-220. [Google Scholar] [CrossRef] - Mills, K. C., Yuan, L., & Jones, R. T. (2011). Estimating the physical properties of slags.
*Journal of the Southern African Institute of Mining and Metallurgy*,*111*, 649-658. [Google Scholar] - Han, Y. U., Min, D. J. (2016). The cationic effect on properties and structure of CaO–MgO–SiO2 melts. Advances in Molten Slags, Fluxes, and Salts: Proceedings of the 10th International Conference on Molten Slags, Fluxes and Salts. pp. 501–509. Springer, Cham. DOI 10.1007/978-3-319-48769-4_53. [CrossRef]
- Courtessole, C., & Etay, J. (2011). Interface, flows and transfer in an electromagnetic process devoted to liquid/liquid extraction.
*Advanced Engineering Materials*,*13(7)*, 556-562. [Google Scholar] [CrossRef] - Zhang, S., Sun, N., He, X., Lu, X., & Zhang, X. (2006). Physical properties of ionic liquids: Database and evaluation.
*Journal of Physical and Chemical Reference Data*,*35(4)*, 1475-1517. [Google Scholar] [CrossRef] - Ernst, R., Mangelinck-Noel, N., Hamburger, J., Garnier, C., & Ramoni, P. (2005). Grain size reduction by electromagnetic stirring inside gold alloys.
*European Physical Journal Applied Physics*,*30(3)*, 215-222. [Google Scholar] [CrossRef] - Davidson, P. A. (2001). An introduction to magnetohydrodynamics. UK: Cambridge University Press.

Appendix 1

Appendix 2: Mathematical Description of the Modeling of the Magnetohydrodynamic Problem

• Physical equations and boundary conditions:

1/Electromagnetics:

— activated in all domains (water+air)/frequency domain solver/quadratic discretization

— Equation:

— with: A: Complex vector potential/u: Velocity/B: Magnetic flux density/

— Boundary conditions: Outer boundaries: Magnetic insulation

2/Heat transfer:

— activated in the water domain/time dependent solver/linear discretization

— Equation:

— with: T: Temperature/u: Velocity/t: Time/

— Boundary conditions:

3/Fluid mechanics:

— activated in the water domain/time dependent solver/incompressible/laminar flow/Boussinesq approximation/P1 + P1 discretization/Streamline diffusion + crosswind diffusion (Navier Stokes equation)

— Equations:

— with: u: Velocity/t: Time/

— Boundary conditions: Lower + side boundaries: No slip wall/upper boundary: Open boundary

• Main model numerical data:

— total number of mesh elements: 2118

— total number of degrees of freedom: 6721

— stored time step of dependent solver: 1 s

— maximum Reynolds number 3800.

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. |