Survey
* Your assessment is very important for improving the work of artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the work of artificial intelligence, which forms the content of this project
Physical and Numerical Simulation of Geotechnical Engineering 9th Issue, Dec. 2012 Analysis on Coupled Fields Inside Passage by Finite Element Method and Its Applications WANG Hao, LIU Lily College of Science, Liaoning Technical University, Fuxin, Liaoning 123000, P.R.China [email protected] ABSTRACT: By analyzing the dynamic flow driven by the electrolyte solution in the magnetic field and electric field that generated respectively by solenoid and parallel plate electrodes, this article did the electromagnetic field calculation and flow field calculation via ANSYS software, and the results from the electromagnetic coupling analysis were obtained thereby. KEYWORDS: ANSYS, Solenoid, Electromagnetic force, Electrolyte solution I NTRODUCTION The electrolyte solution inside passage driven by electromagnetic force is a coupled field composed of electric field, magnetic field and flow field, the electric field will be generated if the powered plates of both positive and negative polarity are parallelized on both sides of the passage. This passage is also inside the long straight solenoid and the magnetic field can be generated by electrifying the solenoid. The electromagnetic interaction will create Lorentz force, which will change the distribution of the flow field and boost forward the electrolyte solution. According to the finite element solution process of the flow field and electromagnetic field obtained by analyzing the flow field and electromagnetic field via finite element method, the coupled field inside the passage can be calculated. The general steps of calculating the coupled field are showed as Figure 1. Begin Determine the finite element analysis model Initialize the flow field areas, parameters and variables Calculate the distribution of the electromagnetic field at t0 Calculate the electromagnetic force upon the electrolyte solution Put the electromagnetic force into the fluid momentum equation; calculate the flow field distribution at ti Put Vti into the Maxwell Equation; calculate the electromagnetic field distribution at ti Y ▏Vti-Vti-T ▏≤δ N End Figure 1 Program design flow of calculating the coupled field 2 ANSYS ANALYSIS MODEL The two-dimensional calculation model used in the ANSYS analysis is showed in Figure 2. The passage is in rectangular shape with an interior magnetic field generated © ST. PLUM-BLOSSOM PRESS PTY LTD by solenoid or multiphase alternative superconducting magnet. Two sets of phases with a phase number of 20 are on the upper and lower sides respectively. Under the effect of magnetic field and actuated by it, the electrolyte solution in the passage will flow. Analysis on Coupled Fields Inside Passage by Finite Element Method and Its Applications DOI: 10.5503/J.PNSGE.2012.09.017 Figure 2 Calculation model for ANSYS analysis Preferences: (1) Passage of the electrolyte solution: length=1200mm; width=200mm; height=140mm (2) Solenoid: length=1500mm; diameter=1.5mm; number of windings=1000/2000; conducting wire: enameled copper wire with a electrical resistivity of 1.7× 10-8Ω m; conducted current: 0-15A. (3) Set the electric conductivity of the electrolyte solution as a constant number σ =10.0s/m, with a density of ρ =1100kg/m3 and a viscosity coefficient of μ =0.0017Ns/m2. (4) DC magnetic field follows the positive direction of Y axis, while AC travelling wave magnetic field follows the positive direction of X axis. The exciting current in the coil is J=100A. The end effect at inlet and outlet of the passage is neglected 3 MESH ANALYSIS GENERATION FOR ANSYS For finite element simulation via computation model, the first thing is to generate grid mesh, see Figure 3. As for plane-type two-dimension simulation analysis, the quadrangular cell is used to divide the whole model. Figure 3 Mesh generation for ANSYS analysis when DC electric pressure is 30V. The distribution of electric field vectors in electrolyte solution between the polar plates is showed in Figure 4. 4 CALCULATION OF ELECTRIC FIELD According to the above model, the electric conductivity of the electrolyte solution can be calculated as 10.0s/m Figure 4 Distribution of electric field vectors in electrolyte solution 88 Physical and Numerical Simulation of Geotechnical Engineering 9th Issue, Dec. 2012 map for magnetic intensity vector and magnetic field intensity is showed in Figure 6. Due to the internal diameter of the solenoid is smaller than its length, the magnetic field in the solenoid is short of homogeneity, and the edge effect is relatively evident. The magnetic field distribution follows the solenoid axial line (Z axis) direction can be observed in the below pictures. 5 CALCULATION OF MAGNETIC FIELD The length of the solenoid of the enameled copper wire is 1.5m with a diameter of 0.0015m. The conducted current is 0.5-30A while the number of windings is 2000. The vector graph of magnetic induction density generated in electrified solenoid and distribution map of magnetic induction density are showed in Figure 5 respectively. The Figure 5 Distribution map of magnetic induction density B generated in electrified solenoid Figure 6 Distribution map of magnetic field intensity H generated in electrified solenoid are showed in Figure 7. Due to the differences of electric pressures between the polar plates and the differences of the current intensity passing the solenoid, different vector graph analyzing results for electromagnetic force are obtained and showed in Figure 8. 6 CALCULATION OF ELECTROMAGNETIC FORCE The electric field generated between the polar plates and the magnetic field generated in the solenoid act on the electrolyte solution, driving it by electromagnetic force. The calculation vector graphs of the electromagnetic force 89 Analysis on Coupled Fields Inside Passage by Finite Element Method and Its Applications DOI: 10.5503/J.PNSGE.2012.09.017 Figure 7 Vector analysis for electrolyte solution under electromagnetic force Figure 8 Analysis on electromagnetic force when change the electric pressure between the polar plates and electric current in the solenoid specific time can be calculated. When the electrified electric current inside coils of the solenoid is 5A, the actuating speed distribution of the flow field can be observed in figure 9. 7 CALCULATION OF FLOW FIELD DRIVING By coupling the electromagnetic field inside solenoid with the passage of the electrolyte solution, the actuating speed distribution of the flow field inside passage at a Figure 9 Vector distribution of the actuating speed of the electrolyte solution flow field inside passage 90 Physical and Numerical Simulation of Geotechnical Engineering 9th Issue, Dec. 2012 8 CALCULATION OF FLOW FIELD The resulting actuating speed acts upon the electrolyte solution in the passages that outside electromagnetic field, and via calculation the liquid flow situation is shown in Figure 10. Putting circular cylinder into the electrolyte solution outside passage, along with the adding of the actuating speeds the possibly generable circumfluence phenomena on the circular cylinder can be observed in Figure 11. Different actuating speeds usually mean different size and distribution of the eddy-making upon the circular cylinder. While the flow field distribution of the liquid of same electrolyte solution density at different electric current that generating different actuating speeds can be observed in Figure 12. Figure 10 Flow field speed distribution of the electrolyte solution inside passage Figure 11 Circumfluence distribution of the circular cylinder in electrolyte solution of the exterior passage Figure 12 Circumfluence distribution of the circular cylinder in electrolyte solution at different electric current 91 Analysis on Coupled Fields Inside Passage by Finite Element Method and Its Applications DOI: 10.5503/J.PNSGE.2012.09.017 [1]. 9 CONCLUSIONS Based on gradually perfected finite element ANSYS calculation method, this article mainly analyzed the obtained display, image and flow speed field for the conducting fluid driven flow and image display for Circumfluence of the circular cylinder etc. by combining the equation of electrolyte solution flow driven by electromagnetic force with discretization process, coupling method, model building and special tridimensional calculation. The analytical results for electric field, magnetic field, electromagnetic force and coupled field etc. were obtained thereby. [2]. AUTHOR IN BRIEF [6]. [3]. [4]. [5]. [7]. Wang Hao (1964-), male, Liaoning Linghai, professor and doctor, mainly works on theory and application of electric field and magnetic field. [8]. [9]. REFERENCES 92 HAN Jiang. Fundamental research on alternative magnetic fluid propulsion, Electrical Engineering Institute of Chinese Academy of Sciences, Magisterial Thesis, 2002. ZHAO Gang, LIN Jinfu, JU Zixiang. General situation of hydrokinetics analytical model for magnetic fluid propulsion passage, New Technology for Electrical Engineering and Energy, 1996, 3, 16-20. WANG Hua, CAO Gang. Analysis of fluid and solid coupled module for liquid containers based on ANSYS, Journal of Chongqing science and technology college (Natural Science Edition), 02, 2006. YU Yating, DU Pingan, WANG Zhenwei. Current studies on applications of finite element method, Mechanical Design, 03, 2005. FANG Sheng, ZHANG Peiqiang. Simulation of particle movement for magnetic rheological fluid and its structural evolution under rotating magnetic field, Journal of Chemical Physics, 5, 562-566, 2001. TAN Zuowu, YUN Jialing. Mathematical model for magnetic fluid thruster, Journal of Watercraft, 1997. 37-43. WANG Hao. Study on flowing of electrolyte under electromagnetic field, Liaoning Technical University, 2007. Kim, S J, Lee, Choung M. Investigation of the flow around a circular cylinder under the influence of an electromagnetic field. Exp. Fluids. 2000, 28, 252-260. Seong-jae Kim, Choung Mook Lee. Control of flows around a circular cylinder: suppression of oscillatory lift force. Fluid Dynamics Research. 2001, 29, 47-63.