Forward Problem in 3D Magnetic Induction Tomography (MIT)
|
|
- Osborn Turner
- 5 years ago
- Views:
Transcription
1 Forward Problem in 3D Magnetic Induction Tomography (MIT) Manuchehr Soleimani 1, William R B Lionheart 1, Claudia H Riedel 2 and Olaf Dössel 2 1 Department of Mathematics, UMIST, Manchester, UK, m.soleimani@umist.ac.uk 2 Insitut für Biomedizinische Technik, Universität Karlsruhe (TH), Germany Claudia.Riedel@ibt.uni-karlsruhe.de ABSTRACT Simulation of the forward problem in Magnetic Induction Tomography (MIT) is one of the key issues in development of the technology of MIT. We simulate the process of measurement in MIT using an edge based finite element technique. We also present the simulation of this forward problem by a program, MAFIA based on finite integration technique. In both methods the same geometrical arrangements are simulated and results are compared. Keywords Forward Problem, Magnetic Induction Tomography, Edge Finite Element, Finite Integration Technique 1 INTRODUCTION Magnetic Induction Tomography (MIT) is a relatively new imaging modality in which internal electric and magnetic material property distributions are reconstructed from the measured impedance data from sensing and exciting coils. In medical applications the non-invasive measurement of the impedance of biological tissue can yield data of diagnostic relevance. For example a change in impedance can give an indication of the healing process of wounds and of skin irritations. The traditional way is to apply the current directly to the tissue and measure the voltage with electrodes. This leads to stray capacitances between the electrodes as well as between the ground and the patient especially at frequencies above 500 khz (Netz 1998). If the patient has acute pain it is not always possible to touch the skin, so a non-contact method must be used. In some applications (e.g. impedance measurement in the brain) impedance can hardly be measured with surface electrodes (Netz 1998). Such a sensor could be an instrument, which is moved around on the skin of a patient, in order to detect areas with impedances that deviate from normal. The non-contact measurement is based on the idea that a time varying magnetic field induces eddy currents in the conductive tissue. These eddy currents will create a field by themselves and a change in the signal can be detected. This change is expected to be very small, so a high resolution of the measuring system is necessary (Scharfetter 2000). Before attempting to reconstruct the unknown material property, it is necessary to solve the forward problem. A solution of the eddy current forward problem in medical MIT has been published by Merwa (2003) and Hollhaus (2002), which considers the low conductivity case and uses the A r, A r V formulation (Biro, 1999) along with a second order edge finite element for reduced magnetic vector potential (A r ) and a second order nodal finite element for electric scalar potential (V). The electric and magnetic field in the interior will be needed for the image reconstruction. Also for an accurate reconstruction it is necessary to be able to predict the measurement at the sensing coils for an arbitrary inhomogeneous material. The forward problem of MIT is modelled in terms of the electromagnetic fields. The forward models are three-dimensional problems and one cannot constrain electromagnetic field to a plane. 2 ELECTROMAGNETIC LAWS Forward solution, which means simulation of the exciting and sensing procedure where the materials and exciting currents are known. There are several ways to describe the problem. The first are used in geophysics and use a dipole method, later this method is used (Tarjan 1968; Gencer 1998) for a medical application. Griffiths (1999) and Gencer (1998) are using more analytical methods. Also, numerical methods for the forward calculation have been reported by several authors (e.g. Gencer 1999; Scharfetter 2000; Merwa 2003). The objective of the forward solver is the evaluation of the electromagnetic field under known material distributions. In fact we need to evaluate induced voltage on the sensing coils due to induced current, which is the evaluation of the trans-impedance between sensing and exciting coils. 275
2 Therefore the magnetic field of the excitation coil generates a primary solenoidal magnetic field. If conducting material is in the field, eddy currents are produced. These eddy currents themselves create a secondary field, which is approximately in quadrature to the primary field in medical application (Griffiths 1999). From the calculated electromagnetic field we can calculate the induced voltage, which is measured in practice. Using the eddy current approximation and a formulation in terms of electric field, Maxwell s equations can be written 1 µ E + iω ( σ + iωε ) E = iω J e (1) where E is electric field, J e is excitation current density, σ is conductivity,µ is permeability, ω is angular frequency and ε is permittivity. 3 SOLVING THE DISCRETISIED EDDY CURRENT PROBLEM 3.1 Finite Integration Technique (FIT) We use two numerical methods, first a Finite Integration Technique. The simulations are done using MAFIA (Weiland 1996). The Maxwell equations are written in discrete form, the field quantities E and B on different grids. These two grids are dual to each other. This results in a discrete representation of the curl, div and grad operators. This preserves the analytical properties of curl, div, and grad. i.e. div curl = 0 and curl grad = 0 are fulfilled also by their discrete counterparts (Hanhne 1992). Since we deal in with a time harmonic problem, for linear material the equations can be written in frequency space. Materials can be anisotropic with real ε, µ andσ. Therefore the material properties are described by the diagonal matrices D ε, D µ and D σ in the grid and with the matrices D ε, D µ and D σ in the dual grid. Using equation (1) creates a current density j e in a discrete form D 1 A C D S D 1 µ D 1 A CD S ω 2 D ε + 1 iω D σ e = iωj e (2) The discrete curl and div operators are described by the matrices C and S in the grid and by C and S in the dual grid. The geometrical dataset of the grid and the dual grid are described in diagonal matrices D S and D S that determine the length of one mesh segment, D A and D A are also diagonal matrices that describe the surface (Hanhne 1992; Hanhne and Weiland. 1992). 3.2 Finite Element Method (FEM) The second numerical technique is on Finite Element Method. We chose edge finite element over nodal elements, which has well-established advantages for vector fields (Biro 1999). The electromagnetic field in the eddy current problem can be defined either by field or by potentials. There are different possible combinations of these quantities as state variables (Biro 1999). In different regions different approximations of Maxwell s equations can be used. To exclude the wave propagation effects, we neglect the displacement current, which is a low frequency approximation. In the excitation coil region, we assumed there is no eddy current so we can use an electrostatic approximation. With the same assumption, (no eddy current) the sensing coil is considered as a magnetostatic area. In the software we developed, coils with any shape can be simulated and the result of the coil modelling comes explicitly into the eddy current model. We use a formulation in terms of magnetic vector potential. The space of shape function of the first order edge element (one form Withney elements) includes the space of the gradient of the zero form Withney element which is in fact the nodal FEM (Biro 1999), so we can define a modified magnetic vector potential in conductive region * in which E = iωa, in non conductive region A is used. Galerkin s approximation to equation (1) using the A, A* formulation of the eddy current problem is Ω 1 N A dv + ( i ω ( σ + i ωε ) N A ) dv = ( NJ ) e dv. (3) µ Ω e Ω c 276
3 where N is the first order vector based shape function, Ω, Ω e, Ω c are whole region, eddy current region and coil region respectively, and A is the magnetic vector potential. The linear system of equations arising from this edge FEM formulation is solved using a preconditioned BiCSTAB method. 4.1 Simulation of Measuring Values 4 SIMULATION MODEL Using MAFIA the induced voltage that is measured in practise is calculated. The complex measured voltage U meas can be calculated evaluating the line integral for FIT field results and a volume integral for FEM field results where E is the electric field and C is the integration path, j is the current density for an artificially unit current passing the coil and Ω is the coil s volume. 4.2 Simulation of the Excitation Coil U meas = Edl or U meas = E. jdv (4) C Ω In FIT formulation the excitation coil is also simulated with line segments. These line segments are carrying a current with amplitude and a phase. For this formulation we used only one line segment instead of 22. The current of one winding is multiplied by number of turns. The resulting current density is used in equation (5). In FEM formulation we calculate the current density in the coil analytically in simple coil shapes or by finite element solution for a complex shape. For example for a cylindrical coil shape the current density is J = J me ˆ θ, in which J m is the amplitude of the current density and ê θ is the unit vector in tangential direction. The modelling of the current source is important for the accuracy and speed of FEM simulation of the eddy current problem (Ren 1996). 4.3 Boundary Condition In FIT formulation we use open volume boundary condition. For the simulation, 500 iterations are used, which gives us a residual of In FEM software we define far field boundary conditions so we need to mesh area surrounding the coils and conductive object and in the far field we define normal component of the magnetic flux density to be zero. 4.4 Dimension of the Sensor and the Geometrical Arrangement Figure 1 shows the geometrical arrangements for the numerical simulation. The upper coil presents the excitation coil, a coils with 22 windings and an excitation current of 0.2 A. The measurement coil is between the excitation coil and the tissue and represents also a coil with 22 windings. The coils are rectangular coils with length of 10 mm on one side. The sensing coil is 1.5 mm on top of the conductive object. It has 1mm thickness in Z direction. There is 1mm distance between two coils and exciting coil also has 1mm thickness. In figure 1b the geometry of coil is been shown. The volume which is used in calculation is 90 mm x 45 mm x 24 mm, and each of the two conductor blocks are 45 mm x 45 mm x 24 mm. The main regions are vacuum and the conducting tissue. A mesh with equidistant mesh points and a separation of 1 mm is used. The frequency of the alternating current is 1 MHz. The simulations are made with the sensor in the middle to avoid the edge effects. 5.1 Simulation with no Tissue 5 RESULTS Table 1 shows the induced voltage results of both simulations; they are both in the same range. Using an analytical calculation as a crosscheck, we calculate the magnetic field of a round current loop with I = 4.4 A, and R = m by using Biot-Savart s law. From this, at a distance of z = m we have a 277
4 Figure 1: Arrangement of the coil system and the tissue block scale in (m). The excitation and the measurement coils are simulated as rectangular coils. The tissue block is divided into two halves. For the simulation, the tissue was moved under the sensor. dimensions of coil from above. calculated magnetic field B = (T) in the centreline. We use this field to subsequently calculate both the flux Φ and the induced voltage, so that the induced voltage U ind of one winding is -2 given approximately by U = iω Φ = i 2 π f B a = i V, which is 2.22 V for 22 turns coil. ind Induced voltage in the secondary coil Finite Integration Technique i 1.97 V Edge Finite Element Technique i 2.21 V Table 1: Induced voltage in the secondary coil with no tissue under the sensor. The results are calculated with the Finite Integration Technique and the Edge Finite Element Technique. 5.2 Simulation with the Sensor in the Middle Over the Tissue For the following results, the sensor is placed in the middle of the tissue block, as shown in figure 2. So the sensor is right over the area where the conducting parts come together. Induced voltage in the secondary coil Finite Integration Technique Edge Finite Element Technique i i 2.08 Table 2: Induced voltage in the secondary coil with tissue. The results are calculated with the Finite Integration Technique and the Edge Finite Element Technique Figure 2: Edge FEM software Imaginary part of the Magnetic flux density in x-z plane, for a current density 1 (A/m 2 ) from exciting coil and magnetic vector potential for free space in x-y plane. 278
5 Figure 3: Eddy current density FIT shows the x-y plane with the imagery part of the current density in the top layer of the tissue. The maximum arrow presents 8.6 A/m 2 shows the y component of the current density cut in the x z plane in A/m 2. Figure 4: Real part of the magnetic flux density FIT: shows the x z plane. The maximum arrow represents Vs/m 2. shows the magnitude of the real part of the magnetic flux density in 3 Vs/m 2. 6 CONCLUSION The main objective of this study is to compare the numerical results from two forward solutions programs and also better understanding of the numerical behaviour of the eddy current problem. One is the commercial software MAFIA, the other is a home made edge based finite element program. The numerical comparisons show a good agreement. One of the main differences is the simulation of the coil and we suggest that the differences in results are mainly due to this. In numerical simulation of the forward problem of MIT, the modelling of the coil is important in order to have accurate prediction of the induced voltage. Modelling of the coil also has effect on the speed of the forward solver. If we mesh the coil region we must use a dense enough mesh. To compute an accurate forward model it is also important to have a dense enough (typically at least 3 elements) in skin depth region. Too coarse a mesh in skin depth region produces inaccurate sensed voltages. 7 REFERENCES BIRO O., (1999), Edge Element Formulations of Eddy Current Problems, Computer methods in applied mechanics and engineering, 169, pp GENCER N.G. and TEK. N., (1998), Imaging Tissue Conductivity via Contactless Measurements: A Feasibility Study, Turkish Journal of Electrical Engineering ELEKTRIK, 6(3) pp GENCER N.G. and TEK. M.N., (1999), Forward Problem Solution for Electrical Conductivity Imaging via Contactless Measurements, Physiological Measurement, 44, pp
6 GRIFFITHS H., STEWART W.R. and GOUGH. W., (1998), Magnetic Induction Tomography: Measurements with a Single Channel, In Proc: X-ICEBI, Barcelona, pp GRIFFITHS H., STEWART W.R. and GOUGH W., (1999), Magnetic Induction Tomography Electrical Bioimpedance Methods: Applications to Medicine and Biotechnology, 873, pp GRIFFITHS H., (2001), Magnetic Induction Tomography, Meas. Sci. Technol., 12(8) pp HAHNE P. and WEILAND T., (1992a), 3D Eddy Current Computation in the Frequency Domain regarding the Displacement Current, IEEE Transactions on Magnetics, 28(2) pp HAHNE P., (1992b), Zur numerischen Berechnung zeitharmonischer elektromagnetischer Felder, Ph.D. Thises, Techn. Univ. Darmstadt HOLLHAUS K., MAGELE C., BRANDSTÄTTER B., MERWA R., and SCHARFETTER H., (2002), Numerical Simulation of the Forward Problem in Magnetic Induction Tomography of Biological Tissue, In Proc: 10th IGTE Symposium, pp MERWA R., HOLLAUS K., BRANDSTATTER B. and SCHARFETTER H., (2003), Numerical Solution of the General 3D Eddy Current Problem for Magnetic Induction Tomography (Spectroscopy), Physiological. Measurement, 24(2), pp NETZ J., FORNER E. and HAAGEMANN S., (1993), Contactless Impedance Measurement by Magnetic Induction, a Possible Method for Investigation of Brain Impedance, Physiological Measurement, 14(4) pp REN Z., (1996), Influence of the R.H.S on the Convergence Behaviour of the Curl-Curl Equation, IEEE Trans. MAG, 32(3) May. SCHARFETTER H., LACKNER H.K. and ROSELL J., (2000), High Resolution Hardware for Magnetic Induction Tomography (MIT) at Low Frequencies, In Proc: EIT Conference, London, UK. TARJAN P.P. and MCFEE R., (1968), Electrodeless Measurements of the Effective Resistivity of the Human Torso and Head by Magnetic Induction, IEEE Trans. BME, 15(4), pp WEILAND T., (1996), MAFIA version 4, In Proc: Computational Accelerator Physics Conference (CAP), pp
Sensitivity Analysis of 3D Magnetic Induction Tomography (MIT)
Sensitivity Analysis of 3D Magnetic Induction Tomography (MIT) W R B Lionheart 1, M Soleimani 1, A J Peyton 2 1 Department of Mathematics UMIST, Manchester, UK, Email: bill.lionheart@umist.ac.uk, 2 Department
More informationPRELIMINARY RESULTS ON BRAIN MONITORING OF MENINGITIS USING 16 CHANNELS MAGNETIC IN- DUCTION TOMOGRAPHY MEASUREMENT SYSTEM
Progress In Electromagnetics Research M, Vol. 24, 57 68, 212 PRELIMINARY RESULTS ON BRAIN MONITORING OF MENINGITIS USING 16 CHANNELS MAGNETIC IN- DUCTION TOMOGRAPHY MEASUREMENT SYSTEM H. J. Luo 1, *, W.
More informationElectromagnetic Tomography: Real-Time Imaging using Linear Approaches
Electromagnetic Tomography: Real-Time Imaging using Linear Approaches Caeiros, J. M. S. and Martins, R. C. October 17, 2010 Abstract The objective of the present work is to apply linear image reconstruction
More informationForward problem solution for electrical conductivity imaging via contactless measurements
Phys. Med. Biol. 44 (1999) 927 940. Printed in the UK PII: S0031-9155(99)96674-5 Forward problem solution for electrical conductivity imaging via contactless measurements Nevzat G Gençer and M Nejat Tek
More informationAnalysis and Imaging in Magnetic Induction Tomography using the Impedance Method
Journal of Physics: Conference Series Analysis and Imaging in Magnetic Induction Tomography using the Impedance Method To cite this article: Julia Grasiela B Wolff et al 2012 J. Phys.: Conf. Ser. 407 012022
More informationUnit-1 Electrostatics-1
1. Describe about Co-ordinate Systems. Co-ordinate Systems Unit-1 Electrostatics-1 In order to describe the spatial variations of the quantities, we require using appropriate coordinate system. A point
More informationfiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES
Content-ELECTRICITY AND MAGNETISM 1. Electrostatics (1-58) 1.1 Coulomb s Law and Superposition Principle 1.1.1 Electric field 1.2 Gauss s law 1.2.1 Field lines and Electric flux 1.2.2 Applications 1.3
More informationDHANALAKSHMI SRINIVASAN INSTITUTE OF RESEARCH AND TECHNOLOGY
DHANALAKSHMI SRINIVASAN INSTITUTE OF RESEARCH AND TECHNOLOGY SIRUVACHUR-621113 ELECTRICAL AND ELECTRONICS DEPARTMENT 2 MARK QUESTIONS AND ANSWERS SUBJECT CODE: EE 6302 SUBJECT NAME: ELECTROMAGNETIC THEORY
More informationImage and Shape Reconstruction Methods in. Magnetic Induction and Electrical Impedance. Tomography
Image and Shape Reconstruction Methods in Magnetic Induction and Electrical Impedance Tomography A thesis submitted to the University of Manchester for the degree of Doctor of Philosophy in Faculty of
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationNon-Iterative Object Detection Methods in Electrical Tomography for Robotic Applications
Non-Iterative Object Detection Methods in Electrical Tomography for Robotic Applications Stephan Mühlbacher-Karrer, Juliana P. Leitzke, Lisa-Marie Faller and Hubert Zangl Institute of Smart Systems Technologies,
More informationCHAPTER 2. COULOMB S LAW AND ELECTRONIC FIELD INTENSITY. 2.3 Field Due to a Continuous Volume Charge Distribution
CONTENTS CHAPTER 1. VECTOR ANALYSIS 1. Scalars and Vectors 2. Vector Algebra 3. The Cartesian Coordinate System 4. Vector Cartesian Coordinate System 5. The Vector Field 6. The Dot Product 7. The Cross
More informationMedical Physics & Science Applications
Power Conversion & Electromechanical Devices Medical Physics & Science Applications Transportation Power Systems 1-5: Introduction to the Finite Element Method Introduction Finite Element Method is used
More informationExperimental Evaluation of Conductive Flow Imaging Using Magnetic Induction Tomography
Experimental Evaluation of Conductive Flow Imaging Using Magnetic Induction Tomography Lu Ma 1,AndyHunt 2 and Manuchehr Soleimani 1 1 Engineering Tomography Laboratory(ETL), Department of Electronic and
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More informationELECTROMAGNETISM. Second Edition. I. S. Grant W. R. Phillips. John Wiley & Sons. Department of Physics University of Manchester
ELECTROMAGNETISM Second Edition I. S. Grant W. R. Phillips Department of Physics University of Manchester John Wiley & Sons CHICHESTER NEW YORK BRISBANE TORONTO SINGAPORE Flow diagram inside front cover
More informationUNIT-I INTRODUCTION TO COORDINATE SYSTEMS AND VECTOR ALGEBRA
SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) Subject with Code : EMF(16EE214) Sem: II-B.Tech & II-Sem Course & Branch: B.Tech - EEE Year
More informationA Simulation Study for Electrical Impedance Tomography
A Simulation Study for Electrical Impedance Tomography Md Ahaduzzaman Khan, Dr Ahsan-Ul-Ambia, Md Ali Hossain, and Md Robiul Hoque Abstract Electrical Impedance Tomography (EIT) is the technique to visualize
More informationEvaluation of Material Plate Proprieties Using Inverse Problem NDT Techniques
Evaluation of Material Plate Proprieties Using Inverse Problem NDT Techniques M. CHEBOUT, A.SADOU, L. AOMAR, M.R.MEKIDECHE E-mail chebout_med@yahoo.fr Laboratoire d Etudes et de Modélisation en Electrotechnique,
More informationOptimization of Skin Impedance Sensor Design with Finite Element Simulations
Excerpt from the Proceedings of the COMSOL Conference 28 Hannover Optimization of Skin Impedance Sensor Design with Finite Element Simulations F. Dewarrat, L. Falco, A. Caduff *, and M. Talary Solianis
More informationSensibility Analysis of Inductance Involving an E-core Magnetic Circuit for Non Homogeneous Material
Sensibility Analysis of Inductance Involving an E-core Magnetic Circuit for Non Homogeneous Material K. Z. Gomes *1, T. A. G. Tolosa 1, E. V. S. Pouzada 1 1 Mauá Institute of Technology, São Caetano do
More informationElectromagnetism. 1 ENGN6521 / ENGN4521: Embedded Wireless
Electromagnetism 1 ENGN6521 / ENGN4521: Embedded Wireless Radio Spectrum use for Communications 2 ENGN6521 / ENGN4521: Embedded Wireless 3 ENGN6521 / ENGN4521: Embedded Wireless Electromagnetism I Gauss
More informationUpdate On The Electromagnetism Module In LS-DYNA
12 th International LS-DYNA Users Conference Electromagnetic(1) Update On The Electromagnetism Module In LS-DYNA Pierre L'Eplattenier Iñaki Çaldichoury Livermore Software Technology Corporation 7374 Las
More informationSampling methods for low-frequency electromagnetic imaging
Sampling methods for low-frequency electromagnetic imaging Bastian Gebauer gebauer@math.uni-mainz.de Institut für Mathematik, Joh. Gutenberg-Universität Mainz, Germany Joint work with Jin Keun Seo, Yonsei
More informationMIDSUMMER EXAMINATIONS 2001
No. of Pages: 7 No. of Questions: 10 MIDSUMMER EXAMINATIONS 2001 Subject PHYSICS, PHYSICS WITH ASTROPHYSICS, PHYSICS WITH SPACE SCIENCE & TECHNOLOGY, PHYSICS WITH MEDICAL PHYSICS Title of Paper MODULE
More informationFinite Element Modeling of Electromagnetic Systems
Finite Element Modeling of Electromagnetic Systems Mathematical and numerical tools Unit of Applied and Computational Electromagnetics (ACE) Dept. of Electrical Engineering - University of Liège - Belgium
More informationUNIT I ELECTROSTATIC FIELDS
UNIT I ELECTROSTATIC FIELDS 1) Define electric potential and potential difference. 2) Name few applications of gauss law in electrostatics. 3) State point form of Ohm s Law. 4) State Divergence Theorem.
More informationCoupling Impedance of Ferrite Devices Description of Simulation Approach
Coupling Impedance of Ferrite Devices Description of Simulation Approach 09 May 2012 TU Darmstadt Fachbereich 18 Institut Theorie Elektromagnetischer Felder Uwe Niedermayer 1 Content Coupling Impedance
More informationUniversity of Huddersfield Repository
University of Huddersfield Repository Leeungculsatien, Teerachai, Lucas, Gary and Zhao, Xin A numerical approach to determine the magnetic field distribution of an electromagnetic flow meter. Original
More informationScientific Computing
Lecture on Scientific Computing Dr. Kersten Schmidt Lecture 4 Technische Universität Berlin Institut für Mathematik Wintersemester 2014/2015 Syllabus Linear Regression Fast Fourier transform Modelling
More informationUse of the magnetic field generated by the internal distribution of injected currents for Electrical Impedance Tomography (MR-EIT)
ELEKTRİK, VOL.6, NO.3, 1998, c TÜBİTAK Use of the magnetic field generated by the internal distribution of injected currents for Electrical Impedance Tomography (MR-EIT) Y. Ziya İder & Özlem Birgül Middle
More informationElectromagnetic Analysis Applied to the Prediction of Stray Losses in Power Transformer
Electromagnetic Analysis Applied to the Prediction of Stray Losses in Power Transformer L. Susnjic 1), Z. Haznadar 2) and Z. Valkovic 3) 1) Faculty of Engineering Vukovarska 58, 5 Rijeka, Croatia, e-mail:
More informationIntroduction to Electromagnetic Theory
Introduction to Electromagnetic Theory Lecture topics Laws of magnetism and electricity Meaning of Maxwell s equations Solution of Maxwell s equations Electromagnetic radiation: wave model James Clerk
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationGeneral review: - a) Dot Product
General review: - a) Dot Product If θ is the angle between the vectors a and b, then a b = a b cos θ NOTE: Two vectors a and b are orthogonal, if and only if a b = 0. Properties of the Dot Product If a,
More informationComparison between Analytic Calculation and Finite Element Modeling in the Study of Winding Geometry Effect on Copper Losses
Comparison between Analytic Calculation Finite Element Modeling in the Study of Winding Geometry Effect on Copper Losses M. Al Eit, F. Bouillault, C. March, G. Krebs GeePs Group of electrical engineering
More informationPart IB Electromagnetism
Part IB Electromagnetism Theorems Based on lectures by D. Tong Notes taken by Dexter Chua Lent 2015 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after lectures.
More informationAssessing COMSOL Performances on Typical Electromagnetic Problems Faced by Turbo-Generator Manufacturers
1 Assessing COMSOL Performances on Typical Electromagnetic Problems Faced by Turbo-Generator Manufacturers E. A. Badea *1, J. Oesterheld 1, S. Friedel 2 and M. Olsson 2 1 Alstom (Switzerland) LTD, 2 COMSOL
More informationINSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad Electronics and Communicaton Engineering
INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad - 00 04 Electronics and Communicaton Engineering Question Bank Course Name : Electromagnetic Theory and Transmission Lines (EMTL) Course Code :
More informationEDDY-CURRENT nondestructive testing is commonly
IEEE TRANSACTIONS ON MAGNETICS, VOL. 34, NO. 2, MARCH 1998 515 Evaluation of Probe Impedance Due to Thin-Skin Eddy-Current Interaction with Surface Cracks J. R. Bowler and N. Harfield Abstract Crack detection
More informationMutual Resistance in Spicelink
. Introduction Mutual Resistance in Spicelink J. Eric Bracken, Ph.D. Ansoft Corporation September 8, 000 In this note, we discuss the mutual resistance phenomenon and investigate why it occurs. In order
More informationA Brief Revision of Vector Calculus and Maxwell s Equations
A Brief Revision of Vector Calculus and Maxwell s Equations Debapratim Ghosh Electronic Systems Group Department of Electrical Engineering Indian Institute of Technology Bombay e-mail: dghosh@ee.iitb.ac.in
More informationCHETTINAD COLLEGE OF ENGINEERING & TECHNOLOGY NH-67, TRICHY MAIN ROAD, PULIYUR, C.F , KARUR DT.
CHETTINAD COLLEGE OF ENGINEERING & TECHNOLOGY NH-67, TRICHY MAIN ROAD, PULIYUR, C.F. 639 114, KARUR DT. DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING COURSE MATERIAL Subject Name: Electromagnetic
More informationDevelopment and analysis of radial force waves in electrical rotating machines
DOI: 10.24352/UB.OVGU-2017-098 TECHNISCHE MECHANIK, 37, 2-5, (2017), 218 225 submitted: June 20, 2017 Development and analysis of radial force waves in electrical rotating machines S. Haas, K. Ellermann
More informationVirtual Prototyping of Electrodynamic Loudspeakers by Utilizing a Finite Element Method
Virtual Prototyping of Electrodynamic Loudspeakers by Utilizing a Finite Element Method R. Lerch a, M. Kaltenbacher a and M. Meiler b a Univ. Erlangen-Nuremberg, Dept. of Sensor Technology, Paul-Gordan-Str.
More informationElectromagnetics in Medical Physics
Electromagnetics in Medical Physics Part 4. Biomagnetism Tong In Oh Department of Biomedical Engineering Impedance Imaging Research Center (IIRC) Kyung Hee University Korea tioh@khu.ac.kr Dot Product (Scalar
More informationSimulation study of the sensing field in electromagnetic tomography for two-phase flow measurement
Flow Measurement and Instrumentation 16 (2005) 199 204 www.elsevier.com/locate/flowmeasinst Simulation study of the sensing field in electromagnetic tomography for two-phase flow measurement Ze Liu a,,minhe
More informationElectromagnetic Field Theory Chapter 9: Time-varying EM Fields
Electromagnetic Field Theory Chapter 9: Time-varying EM Fields Faraday s law of induction We have learned that a constant current induces magnetic field and a constant charge (or a voltage) makes an electric
More informationTECHNO INDIA BATANAGAR
TECHNO INDIA BATANAGAR ( DEPARTMENT OF ELECTRONICS & COMMUNICATION ENGINEERING) QUESTION BANK- 2018 1.Vector Calculus Assistant Professor 9432183958.mukherjee@tib.edu.in 1. When the operator operates on
More informationHaus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, ISBN:
MIT OpenCourseWare http://ocw.mit.edu Haus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, 1989. ISBN: 9780132490207. Please use the following
More informationCONSIDER a simply connected magnetic body of permeability
IEEE TRANSACTIONS ON MAGNETICS, VOL. 50, NO. 1, JANUARY 2014 7000306 Scalar Potential Formulations for Magnetic Fields Produced by Arbitrary Electric Current Distributions in the Presence of Ferromagnetic
More informationElectrical Impedance Tomography Based on Vibration Excitation in Magnetic Resonance System
16 International Conference on Electrical Engineering and Automation (ICEEA 16) ISBN: 978-1-6595-47-3 Electrical Impedance Tomography Based on Vibration Excitation in Magnetic Resonance System Shi-qiang
More informationModelling of Specimen Interaction with Ferrite Cored Coils by Coupling Semi-Analytical and Numerical Techniques
Modelling of Specimen Interaction with Ferrite Cored Coils by Coupling Semi-Analytical and Numerical Techniques A. Skarlatos, E. Demaldent, A. Vigneron and C. Reboud 25 th 28 th of June, Bratislava, Slovak
More informationAN INDEPENDENT LOOPS SEARCH ALGORITHM FOR SOLVING INDUCTIVE PEEC LARGE PROBLEMS
Progress In Electromagnetics Research M, Vol. 23, 53 63, 2012 AN INDEPENDENT LOOPS SEARCH ALGORITHM FOR SOLVING INDUCTIVE PEEC LARGE PROBLEMS T.-S. Nguyen *, J.-M. Guichon, O. Chadebec, G. Meunier, and
More informationApproximate Analytical Calculation of the Skin Effect in Rectangular Wires
Approimate Analytical Calculation of the Skin Effect in Rectangular Wires Dieter Gerling University of Federal Defense Munich, 85579 Neubiberg, Germany Dieter.Gerling@unibw.de Abstract- It is well-known
More informationCHAPTER 7 ELECTRODYNAMICS
CHAPTER 7 ELECTRODYNAMICS Outlines 1. Electromotive Force 2. Electromagnetic Induction 3. Maxwell s Equations Michael Faraday James C. Maxwell 2 Summary of Electrostatics and Magnetostatics ρ/ε This semester,
More informationKeywords: Electric Machines, Rotating Machinery, Stator faults, Fault tolerant control, Field Weakening, Anisotropy, Dual rotor, 3D modeling
Analysis of Electromagnetic Behavior of Permanent Magnetized Electrical Machines in Fault Modes M. U. Hassan 1, R. Nilssen 1, A. Røkke 2 1. Department of Electrical Power Engineering, Norwegian University
More informationarxiv: v3 [physics.med-ph] 20 Mar 2012
Hybrid tomography for conductivity imaging Thomas Widlak 1, Otmar Scherzer 1,2 1 Computational Science Center, University of Vienna, Nordbergstr. 15, A-1090 Vienna, Austria 2 RICAM, Austrian Academy of
More informationA 2-Dimensional Finite-Element Method for Transient Magnetic Field Computation Taking Into Account Parasitic Capacitive Effects W. N. Fu and S. L.
This article has been accepted for inclusion in a future issue of this journal Content is final as presented, with the exception of pagination IEEE TRANSACTIONS ON APPLIED SUPERCONDUCTIVITY 1 A 2-Dimensional
More informationCOLLEGE PHYSICS Chapter 23 ELECTROMAGNETIC INDUCTION, AC CIRCUITS, AND ELECTRICAL TECHNOLOGIES
COLLEGE PHYSICS Chapter 23 ELECTROMAGNETIC INDUCTION, AC CIRCUITS, AND ELECTRICAL TECHNOLOGIES Induced emf: Faraday s Law and Lenz s Law We observe that, when a magnet is moved near a conducting loop,
More informationHandling Nonlinearity by the Polarization Method and the Newton-Raphson Technique
Handling Nonlinearity by the Polarization Method and the Newton-Raphson Technique M. Kuczmann Laboratory of Electromagnetic Fields, Department of Telecommunications, Széchenyi István University, Egyetem
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK
KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK SUB.NAME : ELECTROMAGNETIC FIELDS SUBJECT CODE : EC 2253 YEAR / SEMESTER : II / IV UNIT- I - STATIC ELECTRIC
More informationThe Linear Induction Motor, a Useful Model for examining Finite Element Methods on General Induction Machines
The Linear Induction Motor, a Useful Model for examining Finite Element Methods on General Induction Machines Herbert De Gersem, Bruno Renier, Kay Hameyer and Ronnie Belmans Katholieke Universiteit Leuven
More informationNumerical investigation of the electric field distribution induced in the brain by transcranial magnetic stimulation
Numerical investigation of the electric field distribution induced in the brain by transcranial magnetic stimulation D.-H. Kim, N. Loukaides, J.K. Sykulski and G.E. Georghiou Abstract: Results are presented
More informationPulse eddy currents using an integral-fem formulation for cracks detection
International Journal of Applied Electromagnetics and Mechanics 33 (2010) 1225 1229 1225 DOI 10.3233/JAE-2010-1242 IOS Press Pulse eddy currents using an integral-fem formulation for cracks detection Gabriel
More informationAntenna Theory (Engineering 9816) Course Notes. Winter 2016
Antenna Theory (Engineering 9816) Course Notes Winter 2016 by E.W. Gill, Ph.D., P.Eng. Unit 1 Electromagnetics Review (Mostly) 1.1 Introduction Antennas act as transducers associated with the region of
More informationElectromagnetism. Christopher R Prior. ASTeC Intense Beams Group Rutherford Appleton Laboratory
lectromagnetism Christopher R Prior Fellow and Tutor in Mathematics Trinity College, Oxford ASTeC Intense Beams Group Rutherford Appleton Laboratory Contents Review of Maxwell s equations and Lorentz Force
More informationMagnetic Leakage Fields as Indicators of Eddy Current Testing
ECNDT 006 - We.4.3.4 Magnetic Leakage Fields as Indicators of Eddy Current Testing Božidar BRUDAR, International Center for Sustainable Development, Ljubljana, Slovenia Abstract: With eddy current testing
More informationCoupling of eddy-current and circuit problems
Coupling of eddy-current and circuit problems Ana Alonso Rodríguez*, ALBERTO VALLI*, Rafael Vázquez Hernández** * Department of Mathematics, University of Trento ** Department of Applied Mathematics, University
More information1 Fundamentals of laser energy absorption
1 Fundamentals of laser energy absorption 1.1 Classical electromagnetic-theory concepts 1.1.1 Electric and magnetic properties of materials Electric and magnetic fields can exert forces directly on atoms
More informationIndustrial Heating System Creating Given Temperature Distribution
SERBIAN JOURNAL OF ELECTRICAL ENGINEERING Vol. 5, No. 1, May 2008, 57-66 Industrial Heating System Creating Given Temperature Distribution Ilona Iatcheva 1, Ilonka Lilianova 2, Hristophor Tahrilov 2, Rumena
More informationHIGH VOLTAGE TECHNIQUES REVİEW: Electrostatics & Magnetostatics
HIGH VOLTAGE TECHNIQUES REVİEW: Electrostatics & Magnetostatics Zap You walk across the rug, reach for the doorknob and...zap!!! In the winter, when you change your pullover you hear and/or see sparks...
More informationCURRENT-CARRYING CONDUCTORS / MOVING CHARGES / CHARGED PARTICLES IN CIRCULAR ORBITS
PHYSICS A2 UNIT 4 SECTION 4: MAGNETIC FIELDS CURRENT-CARRYING CONDUCTORS / MOVING CHARGES / CHARGED PARTICLES IN CIRCULAR ORBITS # Questions MAGNETIC FLUX DENSITY 1 What is a magnetic field? A region in
More informationAn impedance method for spatial sensing of 3D cell constructs towards applications in tissue engineering
Electronic Supplementary Material (ESI) for Analyst. This journal is The Royal Society of Chemistry 2015 Supporting Information: An impedance method for spatial sensing of 3D cell constructs towards applications
More informationEngineering Electromagnetic Fields and Waves
CARL T. A. JOHNK Professor of Electrical Engineering University of Colorado, Boulder Engineering Electromagnetic Fields and Waves JOHN WILEY & SONS New York Chichester Brisbane Toronto Singapore CHAPTER
More informationIEEE TRANSACTIONS ON POWER DELIVERY, VOL. 22, NO. 1, JANUARY /$ IEEE
IEEE TRANSACTIONS ON POWER DELIVERY, VOL. 22, NO. 1, JANUARY 2007 195 Analysis of Half-Turn Effect in Power Transformers Using Nonlinear-Transient FE Formulation G. B. Kumbhar, S. V. Kulkarni, Member,
More informationHomogenization of the Eddy Current Problem in 2D
Homogenization of the Eddy Current Problem in 2D K. Hollaus, and J. Schöberl Institute for Analysis and Scientific Computing, Wiedner Hauptstr. 8, A-14 Vienna, Austria E-mail: karl.hollaus@tuwien.ac.at
More informationChapter 5. Magnetostatics
Chapter 5. Magnetostatics 5.1 The Lorentz Force Law 5.1.1 Magnetic Fields Consider the forces between charges in motion Attraction of parallel currents and Repulsion of antiparallel ones: How do you explain
More informationElectromagnetics in COMSOL Multiphysics is extended by add-on Modules
AC/DC Module Electromagnetics in COMSOL Multiphysics is extended by add-on Modules 1) Start Here 2) Add Modules based upon your needs 3) Additional Modules extend the physics you can address 4) Interface
More informationMagnetic Field Analysis
NISA - EMAG EMAG is the electromagnetic module of the family of general purpose finite element based program NISA. It can determine electric and magnetic field distributions in a wide class of electromagnetic
More informationELECTRICITY AND MAGNETISM
ELECTRICITY AND MAGNETISM Chapter 1. Electric Fields 1.1 Introduction 1.2 Triboelectric Effect 1.3 Experiments with Pith Balls 1.4 Experiments with a Gold-leaf Electroscope 1.5 Coulomb s Law 1.6 Electric
More information444 Index Boundary condition at transmission line short circuit, 234 for normal component of B, 170, 180 for normal component of D, 169, 180 for tange
Index A. see Magnetic vector potential. Acceptor, 193 Addition of complex numbers, 19 of vectors, 3, 4 Admittance characteristic, 251 input, 211 line, 251 Ampere, definition of, 427 Ampere s circuital
More informationST.JOSEPH COLLEGE OF ENGINEERING,DEPARTMENT OF ECE
EC6403 -ELECTROMAGNETIC FIELDS CLASS/SEM: II ECE/IV SEM UNIT I - STATIC ELECTRIC FIELD Part A - Two Marks 1. Define scalar field? A field is a system in which a particular physical function has a value
More informationMagnetic Induction Tomography Modeling in Biological Tissue Imaging Using Two-Port Network Technique
Sensors & Transducers, Vol. 50, Issue 3, March 03, pp. -9 Sensors & Transducers 03 by IFSA http://www.sensorsportal.com Magnetic Induction Tomography Modeling in Biological Tissue Imaging Using Two-Port
More informationExperiment and Simulation of the Eddy Current NDT on an Aluminium Plate Using a Uniform Field Probe
Experiment and Simulation of the Eddy Current NDT on an Aluminium Plate Using a Uniform Field Probe Luka Kufrin 1,2, A. Lopes Ribeiro 1,2, H. Geirinhas Ramos 1,2, O. Postolache 1 1 Instituto de Telecomunicações,
More informationarxiv: v2 [physics.acc-ph] 27 Oct 2014
Maxwell s equations for magnets A. Wolski University of Liverpool, Liverpool, UK and the Cockcroft Institute, Daresbury, UK arxiv:1103.0713v2 [physics.acc-ph] 27 Oct 2014 Abstract Magnetostatic fields
More information3D Magnetic Scalar Potential Finite Element Formulation for Conducting Shells Coupled with an External Circuit
3D Magnetic Scalar Potential Finite Element Formulation for Conducting Shells Coupled with an External Circuit Christophe Guérin CEDRAT, Meylan, France Gérard Meunier Grenoble Electrical Engineering Laboratory
More informationPhysics GRE: Electromagnetism. G. J. Loges 1. University of Rochester Dept. of Physics & Astronomy. xkcd.com/567/
Physics GRE: Electromagnetism G. J. Loges University of Rochester Dept. of Physics & stronomy xkcd.com/567/ c Gregory Loges, 206 Contents Electrostatics 2 Magnetostatics 2 3 Method of Images 3 4 Lorentz
More informationHomogeneous Bianisotropic Medium, Dissipation and the Non-constancy of Speed of Light in Vacuum for Different Galilean Reference Systems
Progress In Electromagnetics Research Symposium Proceedings, Cambridge, USA, July 5 8, 2010 489 Homogeneous Bianisotropic Medium, Dissipation and the Non-constancy of Speed of Light in Vacuum for Different
More informationr r 1 r r 1 2 = q 1 p = qd and it points from the negative charge to the positive charge.
MP204, Important Equations page 1 Below is a list of important equations that we meet in our study of Electromagnetism in the MP204 module. For your exam, you are expected to understand all of these, and
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationEddy Current Losses in the Tank Wall of Power Transformers
Eddy Current Losses in the Tank Wall of Power Transformers Erich Schmidt Institute of Electrical Drives and Machines, Vienna University of Technology A 14 Vienna, Austria, Gusshausstrasse 25 29 Phone:
More informationElectromagnetic Theory: PHAS3201, Winter 2008 Preliminaries D. R. Bowler drb/teaching.
Electromagnetic Theory: PHA3201, Winter 2008 Preliminaries D. R. Bowler david.bowler@ucl.ac.uk http://www.cmmp.ucl.ac.uk/ drb/teaching.html 1 yllabus The course can be split into three main areas: electric
More informationThe initial magnetization curve shows the magnetic flux density that would result when an increasing magnetic field is applied to an initially
MAGNETIC CIRCUITS The study of magnetic circuits is important in the study of energy systems since the operation of key components such as transformers and rotating machines (DC machines, induction machines,
More informationDHANALAKSHMI SRINIVASAN COLLEGE OF ENGINEERING AND TECHNOLOGY
DHANALAKSHMI SRINIVASAN COLLEGE OF ENGINEERING AND TECHNOLOGY DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING QUESTION BANK III SEMESTER EE 8391 ELECTROMAGNETIC THEORY Regulation 2017 Academic Year
More informationSub-domain finite element method for efficiently considering strong skin and proximity effects
Sub-domain finite element method for efficiently considering strong skin and proximity effects Patrick Dular, Ruth Sabariego, Johan Gyselinck, Laurent Krähenbühl To cite this version: Patrick Dular, Ruth
More informationElectrodynamics Qualifier Examination
Electrodynamics Qualifier Examination January 10, 2007 1. This problem deals with magnetostatics, described by a time-independent magnetic field, produced by a current density which is divergenceless,
More informationINTRODUCTION TO ELECTRODYNAMICS
INTRODUCTION TO ELECTRODYNAMICS Second Edition DAVID J. GRIFFITHS Department of Physics Reed College PRENTICE HALL, Englewood Cliffs, New Jersey 07632 CONTENTS Preface xi Advertisement 1 1 Vector Analysis
More informationNONLINEAR FINITE ELEMENT METHOD IN MAGNETISM
POLLACK PERIODICA An International Journal for Engineering and Information Sciences DOI: 10.1556/Pollack.4.2009.2.2 Vol. 4, No. 2, pp. 13 24 (2009) www.akademiai.com NONLINEAR FINITE ELEMENT METHOD IN
More informationInduction_P1. 1. [1 mark]
Induction_P1 1. [1 mark] Two identical circular coils are placed one below the other so that their planes are both horizontal. The top coil is connected to a cell and a switch. The switch is closed and
More information