Magnetic Field Gradient Optimization for Electronic Anti-Fouling Effect in Heat Exchanger
Magnetic Field Gradient Optimization for Electronic Anti-Fouling Effect in Heat Exchanger
Journal of Electrical Engineering and Technology. 2014. Nov, 9(6): 1921-1927
Copyright © 2014, The Korean Institute of Electrical Engineers
• Received : January 08, 2014
• Accepted : April 21, 2014
• Published : November 01, 2014 PDF e-PUB PubReader PPT Export by style
Article
Author
Metrics
Cited by
TagCloud
Yong Han
Detp. of Measurement Technology and Instrumentation Key Laboratory of Hebei Province, School of Electrical Engineering, Yanshan University. Qinhuangdao, China. (hanyong_hit@163.com)
Shu-tao Wang
Corresponding Author: Detp. of Measurement Technology and Instrumentation Key Laboratory of Hebei Province, School of Electrical Engineering, Yanshan University. Qinhuangdao, China. (wangshutao@ysu.edu.cn)

Abstract
A new method for optimizing the magnetic field gradient in the exciting coil of electronic anti-fouling (EAF) system is presented based on changing exciting coil size. In the proposed method, two optimization expressions are deduced based on biot-savart law. The optimization expressions, which can describe the distribution of the magnetic field gradient in the coil, are the function of coil radius and coil length. These optimization expressions can be used to obtain an accurate coil size if the magnetic field gradient on a certain point on the coil’s axis of symmetry is needed to be the maximum value. Comparing with the experimental results and the computation results using Finite Element Method simulation to the magnetic field gradient on the coil’s axis of symmetry, the computation results obtained by the optimization expression in this article can fit the experimental results and the Finite Element Method results very well. This new method can optimize the EAF system’s anti-fouling performance based on improving the magnetic field gradient distribution in the exciting coil.
Keywords
1. Introduction
The scale problem is a big problem to the devices which have the heat exchange process when they work, such as cooling tower and boiler. The scale problem roots on using hard water. Chemicals were used to prevent scale buildup at first [1 - 2] . But using chemicals can bring many problems such as environmental contamination. Then the electromagnetic anti-fouling method was researched and used as a substitute. The electronic anti-fouling (EAF) system is a very useful electromagnetic field water treatment system to the scale problem. The structure diagram of EAF system is shown as Fig. 1 . When EAF system works, EAF Control Unit will load a high frequency alternate signal to the exciting coil which is winded on the tube. The alternate signal takes a certain power, and the induced electromagnetic field generated by exciting coil treats the mineral solution. But the antifouling mechanism of EAF system is still not all clear as yet. So it is scarce of theoretical direction to make a good use of EAF system. PPT Slide
Lager Image
Structure sketch of the EAF system.
Tombacz E et al [3 , 4] studied the effect of a weak magnetic field on hematite sol in stationary and flowing systems. They found that when the flow direction was vertical to the magnetic field direction vector, the cluster of electrolyte particle was formed in the solution. But the cluster of electrolyte particle wasn’t observed in the solution which was placed in the static magnetic field. They believed that the cluster of electrolyte particle phenomenon was caused by the lorentz force. According to the Tombacz E’s standpoint, there are Ca 2+ , CO 3 2− , Mg 2+ , SO 4 2− , HCO 3 and some other electrolyte particle in the hard water, EAF system can make magnetic field in its exciting coil, therefore all kinds of electrolyte particle can be influenced by lorentz force when the hard water is flowing throw the coil. The lorentz force is a very important reason for magnetic field anti-fouling but not the only reason. The magnetic field gradient (MFG) is another important reason for magnetic field anti-fouling. According to the gradient magnetic separation theory, diamagnetism electrolyte particle or paramagnetism electrolyte particle will be acted on by gradient magnetic force when they are in gradient magnetic field. The gradient magnetic force can be expressed as [5 , 6] : PPT Slide
Lager Image
Where c 0 is an environment constant, K m is magnetic susceptibility, H is magnetic field intensity, d H /dz is the gradient of magnetic field. In the hard water, it is known that Ca 2+ and Mg 2+ are diamagnetism, but CO 3 2− , HCO 3 and SO 4 2− are paramagnetism according to Van Vleck’s theory [7 , 8] . Therefore the direction of gradient magnetic force acting on Ca 2+ and Mg 2+ is opposite with the direction of gradient magnetic force acting on CO 3 2− , HCO 3 and SO 4 2− . Therefore the negative ions and positive ions are formed into clusters respectively, and the clusters in the solution can become many seeds of the scale crystal. Then the scale crystal can grow up in the hard water solution instead of growing on the water tube wall. Some researchers found that the anti-fouling effect caused by gradient magnetic field was much better than the antifouling effect caused by uniform magnetic field [9 , 10] . Therefore the anti-fouling performance of EAF system can be optimized by making a biggest magnetic field gradient in the exciting coil.
The magnetic field in the coil is not uniform, and the ratio of coil radius to coil length can influence the magnetic field gradient in the exciting coil. On the same cross section which is vertical with the axis of symmetry of the coil, magnetic field gradient is almost the same. In addition, the axial component of the magnetic field is much bigger than the radial component of the magnetic field, and the tangential component of the magnetic field is zero. Therefore in the present study the axial magnetic field gradient on the axis of symmetry of exciting coil will be optimized by changing the ratio of coil radius to coil length.
2. Optimization Method Modeling
In this paper the exciting coil is single layer for a tightly wound. The cross profile schematic diagram of the coil is shown as Fig. 2 . In Fig. 2 , 2 l is the coil length, r is coil radius, z is the distance between a certain point on the coil’s axis of symmetry and the central point of the coil. PPT Slide
Lager Image
Cross-section diagram of single layer for a tightly wound solenoid.
The computational formula of the magnetic field intensity on the coil’s axis of symmetry is as  : PPT Slide
Lager Image
where Hz is the axial magnetic field intensity, Bz is the axial magnetic induction intensity, μ0 is vacuum permeability, I is the exciting electric current, n is the number of turns per unit length of the coil. Then the Hz ’s derivative with respect to z is as: PPT Slide
Lager Image
It can be known from Eq. (3) that d H z /d z is the function of r, l and z . When z =0, d H z /d z is 0. Therefore the magnetic field gradient on the midpoint of the coil’s axis of symmetry is always 0. It also can be known from Eq. (3) that d H z /d z ’s maximum value can be calculated only when z is a certain constant. In the present article we set that z < l . If the ratio of coil radius to coil length is wanted when magnetic field gradient on a certain constant z is maximum value, the optimization expression should be deduced under two conditions: when the coil length 2 l is a constant or coil radius r is a constant.
- 2.1 When the coil length 2lis a constant
When the coil length 2 l is a constant, assume PPT Slide
Lager Image
then Eq. (3) can be changed into: PPT Slide
Lager Image
Eq. (4) is the function of α , assume PPT Slide
Lager Image
derivative with respect to α is as: PPT Slide
Lager Image PPT Slide
Lager Image
Then we can obtain the ratio of coil radius to coil length by solving Eq. (6) when d Hz /d z is maximum value.
- 2.2 When the coil radiusris a constant
When the coil radius r is a constant, assume β = PPT Slide
Lager Image
then Eq. (3) can be changed into: PPT Slide
Lager Image
Eq. (4) is the function of β , assume PPT Slide
Lager Image
derivative with respect to β is as: PPT Slide
Lager Image PPT Slide
Lager Image
Then we can obtain the ratio of coil radius to coil length by solving Eq. (9) when d Hz /d z is maximum value.
3. Result Simulation Verification and Discussion
- 3.1 Simulation verification for Eq. (6)
We set 2 l =0.2m, electric current density of the coil is 1250000 A/m 3 , n =500/m, and z is 0.02m, 0.05m and 0.08m. In order to verifying the correctness of Eq. (6), first we use Eq. (6) to calculate the maximum value of d H z /d z and ratio of coil radius to coil length when d Hz /d z is maximum value. Then we use Finite Element Method to calculate the value of d Hz /d z when the ratio of coil radius to coil length is different, and find the maximum value of d Hz /d z and find the ratio of coil radius to coil length when d Hz /d z is maximum.
The Finite Element calculation is done by using Ansys 10.0 software, which has powerful computing capability based on the FEM (finite element method). It can do an accurate calculation of eddy current, inductance, magnetic field intensity and so on  . PLANE53 element is used in building the model. PLANE53 is defined by 8 nodes and has up to 4 degrees of freedom (DOF) per node. These DOFs are viz. the z-component of the magnetic vector potential (AZ), the time-integrated electric scalar potential (VOLT), the electric current (CURR), and the electromotive force (emf)  .
When using Eq. (6), the calculating result is shown as Table 1 .
MFG maximum and the corresponding ratio ofrtolwhenz=0.02m, 0.05m and 0.08m. PPT Slide
Lager Image
MFG maximum and the corresponding ratio of r to l when z=0.02m, 0.05m and 0.08m.
When using Finite Element Method, the calculating result is shown as Tables 2 - Table 4 , and the fitting curves of the results is shown as Figs. 3 - Fig. 5 .
MFG maximum with different ratio ofrtolwhenz=0.02m PPT Slide
Lager Image
MFG maximum with different ratio of r to l when z=0.02m
MFG maximum with different ratio ofrtolwhenz=0.05m PPT Slide
Lager Image
MFG maximum with different ratio of r to l when z=0.05m
MFG maximum with different ratio ofrtolwhenz=0.08m PPT Slide
Lager Image
MFG maximum with different ratio of r to l when z=0.08m PPT Slide
Lager Image
Fitting curve of MFG maximum with different ratio of r to l when z=0.02m PPT Slide
Lager Image
Fitting curve of MFG maximum with different ratio of r to l when z=0.05m PPT Slide
Lager Image
Fitting curve of MFG maximum with different ratio of r to l when z=0.08m
It can be seen from Table 1 and Figs (3 - 5) that the results calculated by Eq. (6) can fit well with finite element calculation results. Therefore the correctness of Eq. (6) is verified.
- 3.2 Simulation verification for Eq. (9)
We set r =0.1m, electric current density of the coil is 1250000 A/m 3 , n =500/m, and z is 0.02m, 0.05m and 0.08m. In order to verifying the correctness of Eq. (6), first we use Eq. (9) to calculate the maximum value of d Hz /d z and ratio of coil radius to coil length when d Hz /d z is maximum value. Then we use Finite Element Method to calculate the value of d Hz /d z when the ratio of coil radius to coil length is different, and find the maximum value of d Hz /d z and find the ratio of coil radius to coil length when d Hz /d z is maximum. When using Eq. (9), the calculating result is shown as Table 5 .
MFG maximum and the corresponding ratio ofltorwhenz=0.02m, 0.05m and 0.08m. PPT Slide
Lager Image
MFG maximum and the corresponding ratio of l to r when z=0.02m, 0.05m and 0.08m.
When using Finite Element Method, the calculating result is shown as Table 6 - Table 8 , and the fitting curves of the results is shown as Fig. 6 - Fig. 8 .
MFG maximum with different ratio ofltorwhenz=0.02m PPT Slide
Lager Image
MFG maximum with different ratio of l to r when z=0.02m
MFG maximum with different ratio ofltorwhenz=0.05m PPT Slide
Lager Image
MFG maximum with different ratio of l to r when z=0.05m
MFG maximum with different ratio ofltorwhenz=0.08m PPT Slide
Lager Image
MFG maximum with different ratio of l to r when z=0.08m PPT Slide
Lager Image
Fitting curve of MFG maximum with different ratio of l to r when z=0.02m PPT Slide
Lager Image
Fitting curve of MFG maximum with different ratio of l to r when z=0.05m PPT Slide
Lager Image
Fitting curve of MFG maximum with different ratio of l to r when z=0.08m
It can be seen from Table 5 and Figs (6 - 8) that the results calculated by Eq. (9) can fit well with finite element calculation results. Therefore the correctness of Eq. (9) is verified.
- 3.3 Experimentation Verification for Eqs. (6) and Eq. (9)
In the present study the validity of Eqs. (6) and Eq. (9) was also verified experimentally
We made the coils with different size according to the Finite Element models which were built in the section A and section B of part III above. The exciting coil was made up by copper wire, and the resistivity of the copper wire was 1.7×10 −8 (Ω · m) when ambient temperature was 25℃. The diameter of the copper wire was 0.002m. The exciting electric current was 1.25 A. The magnetic field gradient on the axis of exciting was measured with a weak magnetic field gradient meter (produced by Beijing Zhong Hui Tian Cheng technology co., LTD). The weak magnetic field gradient meter’s highest resolution was 0.01μT, and its recording sensitivity and registration accuracy were 0.004 μT/cm and ±0.25% respectively. The experimental setup is shown in Fig. 9 . PPT Slide
Lager Image
Schematic diagram of experimental setup
The experiment was divided into two stages.
3.3.1 Stage 1
First we let the length of exciting coil as a constant: 2 l = 0.2m, and we changed the coil’s radius from 0.01m to 0.1m. The magnetic field gradient on the axis point of z=0.02m (z is the distance from the medial point of axis) was measured under different coil radius, and the corresponding coil radius was found out when the magnetic field gradient on the axis point of z was maximum. At last the maximum magnetic field gradient value was recorded.
Then we changed the axis points of z=0.05m and z=0.08m and repeated the experimental steps above.
3.3.2 Stage 2
First we let the radius of exciting coil as a constant: r = 0.1m, and we changed the coil’s length from 0.01m to 0.1m. The magnetic field gradient on the axis point of z=0.02m was measured under different coil length, and the corresponding coil length was found out when the magnetic field gradient on the axis point of z was maximum. At last the maximum magnetic field gradient value was recorded.
Then we changed the axis points of z=0.05m and z=0.08m and repeated the experimental steps above.
The unit of maximum magnetic field gradient value measured by weak magnetic field gradient meter was transformed from μT /cm into A/m 2 according to Eq. (10). PPT Slide
Lager Image
where B is magnetic induction intensity (the SI unit is T ) and H is magnetic field intensity (the SI unit is A/m). μ 0 is vacuum permeability. The experimental results are shown in Tables 9 and Table 10 .
The experiment results of stage 1 PPT Slide
Lager Image
The experiment results of stage 1
The experiment results of stage 2 PPT Slide
Lager Image
The experiment results of stage 2
It can be seen from Tables 9 and Table 10 that experiment results can fit well with the results calculated by Eqs. (6) and Eq. (9) which are shown in Tables 1 and Table 5 .
- 3.4 Discussion
1) According to Eq. (1), the magnetic field gradient is inversely proportional to magnetic field intensity when PPT Slide
Lager Image
is a constant. The magnetic field gradient is bigger when the magnetic field intensity is smaller. The magnetic field intensity is proportional to the exciting current. In the present paper we obtain a mathematical method to optimize the magnetic field gradient in exciting coil. It means that a small exciting current can obtain a big magnetic field gradient, which means a good anti-fouling effect.
2) The method to optimize the magnetic field gradient in exciting coil in the present paper can be used to make a good anti-fouling effect in the wanted area in the water tube or the heat exchanger. It is helpful for some area in the water circulation system which is hard to clear.
3) In this paper the exciting coil is single layer for a tightly wound. This kind of coil is widespread used in many field. Therefore the optimization method can be used in some field that the magnetic field gradient in the coil is needed to be controlled accurately.
4. Conclusion
In this paper a mathematical method to optimize the magnetic field gradient in exciting coil is presented for giving a theoretical basis to the optimization design of exciting coil of EAF system, and the correctness of the expressions, which are used to calculate the magnetic field gradient maximum, is verified by Finite Element Method and experimentation. This optimization method can be used to make a biggest magnetic field gradient in the certain area in the coil and improve the EAF system’s antifouling performance by calculating the appropriate ratio of coil radius to coil length. This optimization method can also be used in other coil-using field which the magnetic field gradient is needed to be optimized.
Acknowledgements
This work was financially supported by the National Natural Science Foundation of China (Grant No. 51408525), National Natural Science Foundation of China (Grant No. 61201110), the Specialized Research Foundation for the doctoral program of higher education of China (Grant No. 20101333120004) and Natural Science Foundation of Hebei Province (Grant No.F2014203224).
BIO Yong Han. He received his B.S. and M.S. degrees from Harbin Institute of Technology, Harbin, China, in 2003 and 2008, respectively and the Ph.D. degree in Instrument Science and Technology also from Institute of Technology in 2013. His research interests are applied electromagnetics and electronic anti-fouling technology. Shu-tao Wang. He received his B.S. and M.S. degrees from Yanshan University, Qinhuangdao, China, in 1999 and 2002, respectively. He received his Ph.D. degree in Instrument Science and Technology from Institute of Technology, Harbin, China, in 2006. His research interests are weak signal detection and processing and mechanical fault diagnosis.
References