Numerical modeling and evaluation of the small magnetometer in low-mass experiment (SMILE)

Full text

(1)2007:034. MASTER'S THESIS. Numerical Modeling and Evaluation of the Small Magnetometer in Low-Mass Experiment (SMILE). Israel Alejandro Arriaga Trejo. Luleå University of Technology Master Thesis, Continuation Courses Space Science and Technology Department of Space Science, Kiruna. 2007:034 - ISSN: 1653-0187 - ISRN: LTU-PB-EX--07/034--SE.

(2) Numerical Modeling and Evaluation of the Small Magnetometer in Low-Mass Experiment (SMILE). Israel Alejandro Arriaga Trejo. Department of Space Science Lule˚ a University of Technology. Julius Maximilians Universität W¨ urzburg. Thesis submitted to the Joint European Master in Space Science and Technology to obtain the degree of Master in Science. September 2007.

(3) Supervisor: Examiner: Examiner:. Dr. Nickolay Ivchenko, Royal Institute of Technology Dr. Johnny Ejemalm, Lule˚ a University of Technology Prof. Dr. Wolfgang Dröge, Julius Maximilians Universität W¨ urzburg.

(4) ACKNOWLEDGEMENTS. The present document reflects the result of 2 years of graduate education abroad my country, Mexico, in the field of Space Science within the Joint European Master in Space Science and Technology coordinated by Lule˚ a University of Technology (LTU) in Sweden. I would not have been able to conclude it without the valuable guidance and supervision of Dr. Nickolay Ivchenko from the Royal Institute of Technology in Stockholm, Sweden and the support of my family and friends. I would like to thank my parents, Leticia and Mario, for the patience they had with me during the long time I prepared in Mexico to pursue a master degree. I always counted with their support and help during the difficult moments experienced in this adventure. Also I would like to thank the motivation I received from the rest of my family in specific from Elizabeth, Socorro, Carmen, Concepción, Silvia, Julieta, Obdulia, Cecilia (both of them), Concepción, Mario and close cousins. I do not want to miss the opportunity to thank also close friends in Mexico that all the time were in contact with me during my stay in Europe: Marisol, Izlian, Joel, Rodrigo, Narciso, Luis Alberto and Ricardo as well colleagues from IPN (UPIITA) and UNAM (Faculty of Science). Regarding my teachers who instructed me in the path of science and engineering areas I should thank Fernado Tellez who made me appreciate the beauty of Fourier/Wavelet analysis and Eugenio Ley Koo from whom I learnt to appreciate symmetries and elegance in physics. Here in Europe I would like to thank the coordinators of the Joint European Master in Space Science and Technology for giving me the opportunity to participate in it, especially Dr. Asta Pellinen-Wannberg, Dr. Victoria Barabash and Mr. Sven Molin. After having concluded studies in W¨ urzburg, Germany I decided to stay in Sweden for my specialisation and I must confess that Kiruna had a profound impact on my life. I will not forget the darkness of the winter night being illuminated by the auroras and the mist which appears to be like a rain of tiny crystals. The year I spent in Kiruna was like living a dream from which I did not want to wake up. At the Institutet för rymdvetenskap, IRV, I met nice lecturers: Martin Bohm, Johnny Ejemalm and Anita Enmark. During these 2 years of the programme I had the opportunity to make good friends: Madhuri Nukala, Daga, Lin Chen, Prut Srissupalak, Hubert, Xao Long. I would like to thank again Dr. Nickolay Ivchenko for all the help I received from the first time I contacted him to ask for the possibility to work under his supervision. I do not have words to thank all his help and guidance I received during my work at the Royal Institute of Technology. It was at KTH where I.

(5) ACKNOWLEDGEMENTS. ii. spent the last 6 months of the programme working my thesis project. It took me less than five minutes to meet there people like Nickolay, Göran Olsson, Ola, Lars, Elizabeth, Anita, Ingeborg, Maria Wedlund (ma amie sophistiquée), Amalia Victorin and Marta Baselga (from whom I learnt the correct usage of the words cutre and patoso, not common in the Spanish spoken in Mexico) but certainly it would take me a whole life to forget them..

(6) ABSTRACT. Fluxgate magnetometers have played a major role in space missions due to their stability, range of operation and low energy consumption. Their principle of operation is relatively simple and easy to implement, a nonlinear magnetic material is driven into saturation by an alternating excitation current inducing a voltage that is modulated by the external field intended to be measured. With the increasing use of nanosatellites the instruments and payload on board have been reduced considerably in size and weight. The Small Magnetometer in Low-Mass Experiment, SMILE, is a miniaturised triaxial fluxgate magnetometer with volume compensation incorporating efficient signal processing algorithms within a field programmable gate array (FPGA). SMILE was designed in collaboration between the Lviv Centre of Institute of Space Research in Ukraine where the sensor was developed and the Royal Institute of Technology (KTH) in Stockholm, Sweden where the electronics used to operate the instrument were designed and programmed. The characteristic dimensions of the SMILE magnetometer and geometry of its parts make impractical the task to find an analytical expression for the voltages induced in the pick-up coils to evaluate its performance. In this thesis, the results of numerical simulations of the SMILE magnetometer using a commercial finite element method (FEM) based software are presented. The results obtained are compared with the experimental data available and will serve as a first step to understand the behaviour of the nonlinear components that could lead to improvements of its design in a future..

(7) CONTENTS. 1. Theoretical Background . . . . . . . . . . . . 1.1 Maxwell’s Equations . . . . . . . . . . . 1.1.1 Magnetostatics . . . . . . . . . . 1.1.2 Induced currents . . . . . . . . . 1.1.3 Inductance . . . . . . . . . . . . 1.1.4 Magnetism in matter . . . . . . . 1.1.5 Boundary conditions . . . . . . . 1.2 A glimpse into the history of magnetism. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. . . . . . . . .. 7 7 8 10 11 11 14 15. 2. Fluxgate Magnetometers . . 2.1 Principle of operation 2.2 Single core sensor . . . 2.2.1 Induced voltage 2.3 Double core sensor . . 2.3.1 Induced voltage 2.4 Volume Compensation. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. 18 18 19 21 23 23 24. 3. Simulation software . . . . . . . . . . . . 3.1 The COMSOL Multiphysics software 3.1.1 3D Electromagnetics Module 3.1.2 Mesh generation . . . . . . . 3.1.3 Solvers . . . . . . . . . . . . . 3.2 Finite Element Method . . . . . . . 3.2.1 Distinctive Features . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. 27 27 27 28 28 29 33. . . . . . . . . . . . . . . in a uniform magnetic . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 34. 5. The SMILE Magnetometer . . . . . . . . . . . . . . . . . . . . . . . . 5.1 Sensor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5.2 Electronics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 41 41 41. 4. Accuracy estimation . . . 4.1 Sphere of permeable flux density . . . . . 4.2 Theoretical solution 4.3 Numerical Solution .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . material . . . . . . . . . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . .. . . . . . . inmersed . . . . . . . . . . . . . . . . . .. 34 34 36.

(8) Contents. 2. . . . . .. 45 45 52 55 63. 7. Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 67. Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 68. Appendices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 69. 8. Code for sphere of permeable material . . . . . . . . . . . . . . . . . .. 70. 9. Code for the simulation of the double core . . . . . . . . . . . . . . . .. 76. 6. Results . . . . . . . . . 6.1 Fluxgate element . 6.2 Double core . . . . 6.3 Excitation circuit . 6.4 Compensation coils. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . ..

(9) LIST OF FIGURES. 1.1. 2.1. 2.2 2.3. 2.4. 2.5. 2.6. 2.7 2.8 3.1 3.2. The upper figure shows the dependance of the relative permeability µr with the magnetic field for a nonlinear material. In the bottom figure, a characteristic B − H curve is shown. . . . . . . In the single core configuration a nonlinear magnetic material is surrounded by excitation coils where the current that saturates the core periodically flows. . . . . . . . . . . . . . . . . . . . . . . B − H curve of the core, the saturation is attained for Hmax = 10A/m with Bsat = 0.5T. . . . . . . . . . . . . . . . . . . . . . . (Top) Excitation field with amplitude H0 = 15A/m and period τ = 125µs. (Centre) Magnetic flux density inside the magnetic material. No external magnetic field is present. (Bottom) Induced voltage in the sensing coil for N = 1 and Asns = 1 cm2 . . (Top) Excitation field with amplitude H0 = 15A/m and period τ = 125µs. (Centre) Magnetic flux density inside the magnetic material. An external magnetic field Hext = 3A/m is present. (Bottom) Induced voltage in the sensing coil for N = 1 and Asns = 1 cm2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Double core sensor used in the SMILE magnetometer. The cores are two tapes of an amorphous alloy with dimension 16x1x0.02mm. The excitation coils are wound around each core and are connected in series (800turns). The pick-up coils surround the whole configuration. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Excitation field driving the cores. Both cores are wound by the same coil in opposite direction and are consdiered to be identical, the excitation fields satisfy the relation Hexc1 = −Hexc2 . The period of the excitation field is τ = 125µs. . . . . . . . . . . . . . Magnetic flux density inside the cores when an external magnetic field Hext = 3A/m is present. . . . . . . . . . . . . . . . . . . . . Induced voltage in the sensing coil for N = 1, Asns = 2 cm2 and τ = 125µs for the double core sensor. . . . . . . . . . . . . . . . Mesh of the SMILE model in COMSOL. The geometry has been decomposed in a total of 68873 elements. . . . . . . . . . . . . . Base function and its derivative for k = 2, when the geometry has been meshed with n = 5 elements. . . . . . . . . . . . . . . .. 15. 19 22. 22. 23. 24. 25 25 26 29 30.

(10) List of Figures. 3.3. 4.1. 4.2. 4.3. 4.4. 4.5 5.1 5.2. 5.3 6.1 6.2 6.3. 6.4 6.5. 6.6 6.7. Solution of Poisson equation in the unit interval [0, 1] using the FEM for ρ = 1 and ǫ = 1. The approximate solutions were calculated using 5 and 10 mesh elements of the unit interval. . . Sphere of permeable material inmersed in an uniform magnetic flux density of 1mT. The radius of the sphere is 1cm and µr = 4 × 104 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Comparison of the numerical and theoretical solution at z = 45mm. The length of the simulation box is 30cm and the SPOOLES solver together with a fine mesh were used to obtain the numerical solution. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Comparison of the numerical and theoretical solution at z = 20mm. The length of the simulation box is 30cm and the SPOOLES solver together with a fine mesh were used to obtain the numerical solution. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Comparison of the numerical and theoretical solution at z = 0mm. The length of the simulation box is 30cm and the SPOOLES solver together with a fine mesh were used to obtain the numerical solution. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Estimated error for the numerical solutions obtained using different meshing modes. . . . . . . . . . . . . . . . . . . . . . . . . . . SMILE magnetometer. The sensor has a mass of 21g and dimensions of 2cm per side. . . . . . . . . . . . . . . . . . . . . . . . . Electronic board containing the FPGA, microcontroller and additional components used to operate the instrument. The board was designed at the deprtament of Space and Plasma Physics in The Royal Institute of Technology (KTH). . . . . . . . . . . . . Block diagram of the different modules implemented in the electronic board. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . COMSOL model of the fluxgate element with corresponding pickup coils used in the SMILE magnetometer. . . . . . . . . . . . . Measured inductance of the three pick-up coils of Sensor 6. . . . Magnetic flux density produced when a current of 10µA flows through the pick-up coils. The number of elements used to mesh the geometry is 43088. . . . . . . . . . . . . . . . . . . . . . . . . Inductance of the pick-up coils as a function of the relative permeability of the cores. . . . . . . . . . . . . . . . . . . . . . . . . Magnetic flux through a circular region with radius b = 8cm situated at z = 0cm. The radius of the sphere is a = 1cm and the magnitude of the induction field B0 = 1mT. . . . . . . . . . . B-H curves for different values of the α parameter. The maximum induction field is 0.36T. . . . . . . . . . . . . . . . . . . . . . . . Inductance of the pick-up coils for different values of the α parameter in model (6.6). . . . . . . . . . . . . . . . . . . . . . . .. 4. 32. 35. 37. 38. 39 40 42. 43 43 46 46. 47 48. 49 50 51.

(11) List of Figures. 6.8 6.9 6.10 6.11 6.12 6.13 6.14 6.15 6.16 6.17 6.18 6.19. 6.20 6.21 6.22. 6.23 6.24 6.25. 6.26 6.27. Measured and simulated inductance of the pick-up coils for Sensor 6. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Measured and simulated inductance for Sensor 6 for different turns in the pick-up coils. . . . . . . . . . . . . . . . . . . . . . . Measured and simulated inductance for Sensor 6 for different saturation values. . . . . . . . . . . . . . . . . . . . . . . . . . . . . Nonlinear cores with respective excitation coils. Each excitation coils consists of 800 turns. . . . . . . . . . . . . . . . . . . . . . Magnetic flux through the pick-up coils for different values of the excitation current and external induction field. . . . . . . . . . . Distribution of the magnetic flux density [T] along the axis of one core when no external magnetic field is present. . . . . . . . . . . Magnetic flux density (norm) [T] inside the cores for different excitation currents and external fields. . . . . . . . . . . . . . . . Magnetic field (norm) [A/m] inside the cores for different excitation currents and external fields. . . . . . . . . . . . . . . . . . . Relative permeability of the cores for different excitation currents and external fields. . . . . . . . . . . . . . . . . . . . . . . . . . . Excitation current used to saturate the cores. . . . . . . . . . . . Magnetic flux and induced voltage in the pick-up coils when different external induction fields are present. . . . . . . . . . . . . Resonant circuit used to generate the excitation circuit. The circuit elements r and Lc represent the input resistance and inductance of the excitation coils respectively. . . . . . . . . . . . . Inductance of the excitation coils. . . . . . . . . . . . . . . . . . . Block diagram model of the resonant circuit shown in Figure 6.19 used in Simulink to determine the excitation current. . . . . . . . (Top). A periodic rectangular pulse modulated with a ramp is used to drive the resonant circuit. (Bottom) Excitation current generated. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Excitation current generated for different amplitudes of the voltage source. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Compensation coils. . . . . . . . . . . . . . . . . . . . . . . . . . Compensation field [T] produced when a 1mA current is used to bias the coils. The geometry was partitioned with a total of 29257 elements. . . . . . . . . . . . . . . . . . . . . . . . . . . . . Compensation field [T] along the axis of axis of symmetry of the fluxgate element. . . . . . . . . . . . . . . . . . . . . . . . . . . . Linear model used to fit the data in Table 6.2. . . . . . . . . . .. 5. 53 53 54 54 55 56 57 58 59 60 60. 61 61 62. 62 63 64. 64 65 66.

(12) LIST OF TABLES. 6.1 6.2. Inductance of the pick-up coils. . . . . . . . . . . . . . . . . . . . Compensating field for different currents driving the compensation coils. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 47 66.

(13) 1. THEORETICAL BACKGROUND. The present chapter is intended to settle the theory upon which the SMILE magnetometer bases its principle of operation. An overview of magnetism of steady currents is adresssed at the beginning followed by a description of magnetic materials. The scope used to discuss the properties of magnetism in matter is a macroscopic one (phenomenological) without going in details of quantum properties of atoms that constitute magnetic materials. The material introduced in this chapter is based on references [15], [8], [23], [14] and [17]. For the historical approach references [26] and [27] were consulted.. 1.1 Maxwell’s Equations It is an experimental fact that electromagnetic interactions in a physical system can be described by Maxwell’s equations ρ ǫ0. (1.1). ∇·B= 0. (1.2). ∇·E=. ∇×E+. ∂B =0 ∂t. (1.3). 1 ∂E = µ0 J (1.4) c2 ∂t with appropiate boundary conditions according to the geometry analysed and the so called constitutive relations which provide information about the medium. ∇×B−. D = ǫE. (1.5). B = µH. (1.6). J = σE. (1.7). The E and B fields completely charaterise the electric and magnetic fields in the system..

(14) 1. Theoretical Background. 8. In the systems which are the object of our study the fields do not vary rapidly in time and the effects of the finite velocity of propagation of electric and magnetic fields can be considered to be instantaneous at two distinct points, known as the quasistatic approximation. With these assumptions is possible to simplify Maxwell’s equations: ρ ǫ0. (1.8). ∇·B= 0. (1.9). ∇·E=. ∇×E+. ∂B =0 ∂t. ∇ × B = µ0 J. (1.10) (1.11). The set of equations (1.8)-(1.11) together with (1.5) - (1.7) will be used to describe quantitatively the principle of operation of fluxgate magnetometers. 1.1.1 Magnetostatics Although it is completely valid to assert that a magnetic field is just a relativistic effect of an electric field described in a moving reference frame with respect to an inertial frame where no magnetic interaction is detected, a more conventional treatment is followed here. During the year of 1820 Andre-Marie Ampère, after having attended the presentation of Fran¸cois Arago at the Paris Academy of Sciences, completed the work initiated by Hans Christian Ørsted by conducting experiments with current carrying wires[27]. In one of his experiments Ampère discovered that two rectilinear wires attract or repell each other depending on the direction of the currents that flow through them. This effect is not restricted to straight wires only, indeed any two closed circuits carrying steady currents interact with each other as Ampère realised. The force excerted by a circuit C1 carrying a current i1 on a circuit C2 with current i2 is given by the Biot-Savart law Z Z µ0 dl2 × (dl1 × r) FC1 →C2 = i1 i2 (1.12) 4π r3 C1 C2 the constant µ0 is the permeability of free space and is defined as µ0 = 4π × 10−7 newton/ampere2 It is possible to rewrite (1.12) as Z Z µ0 dl1 × r FC1 →C2 = i2 dl2 × i1 4π r3 C2 C1. (1.13). (1.14). The integral on the right hand side of equation (1.14) defines a vector valued function B : S ⊂ R3 → R3 , the magnetic flux density,.

(15) 1. Theoretical Background. µ0 B= i1 4π. Z. C1. dl1 × r r3. 9. (1.15). which is specified by the geometry and current flowing through C1 . By symmetry, the circuit C2 also produces a magnetic flux density as a result of the current circulating in it. The magnetic flux density is measured in Teslas [T] in the International System of Units. In general, for a distribution of steady currents in a region V in space, the magnetic flux density is given by Z µ0 J×r 3 B= d r (1.16) 4π V r3. where J is the current density distribution. As can be shown by direct substitution, the magnetic flux density defined by (1.16) is consistent with (1.9). The fact that the divergence of B is always zero has important physical consequences. First of all, the non existence of free magnetic monopoles1 and second, the field lines of the magnetic flux density always form closed trajectories, i.e., they can not originate from a point and diverge. Since every vector valued function with zero divergence can be expressed as the curl of a vector field it is possible to write B = ∇ × (A + ∇ψ). (1.17). for some vector field A and scalar field ψ. The vector valued function A : V ⊂ R3 → R3 proposed in (1.17) is known in the literature as the magnetic vector potential and it has been shown that its usage does not change the physics of the problem[8]. The equation that should satisfy the magnetic vector potential A is obtained by substituting (1.17) in (1.11): ∇2 A = ∇(∇ · A) − µ0 J. (1.18). If the Coulomb gauge is used, ∇ · A = 0, Poisson’s equation results and the magnetic potential can be calculated as Z µ0 J 3 A= d r (1.19) 4π V r the unit of A is the Ampere [A]. When the vector field B of a distribution of currents in space is calculated at points situated far from the sources, the denominator in (1.16) or equivalently in (1.19) can be expanded in series leading to an expansion of B (or A). For the magnetic vector potential we have, 1 Several experiments have been conducted to detect magnetic monoples. On February 14, 1982 Blass Cabrera at Stanford University recorded an event that would have indicated the pass of a monopole through his laboratory with his instruments [3]..

(16) 1. Theoretical Background. A = AM + AD + AO + . . .. 10. (1.20). The first term of the expansion in A always vanishes meanwhile the second term depends on the distribution of the currents in space. This second term is referred as the magnetic dipole moment of the configuration µ0 m × r 4π r3 where the vector m is the dipole moment and is defined as Z 1 m= r × Jd3 r 2 V AM =. (1.21). (1.22). The magnetic dipole moment plays an important role in physics as the dominant term in expression (1.20). It is the term that has the major contribution to the magnetic field at distances far from the sources since the remaining ones decrease considerably fast. Earth’s magnetic field is well described with a dipole moment m = 8 × 1022Am2 , tilted about 11 degrees with respect to the rotation axis. The magnetic field at the surface of the Earth is about 30µT at the equator and up to 60µT near the poles. It is on the top of this geomagnetic background that small disturbances due to auroral currents (space physics) or crustal mineral deposits (in exploration geophysics) need to be measured. Accurate vector measurements including stable and precise orientation of the vector magnetometer axes are required since the disturbances created by auroral currents are in the order of 50 − 500nT. The Earth is not the only planet with a characteristic dipole field, in other planets the dipole field dictates the dynamics of charged particles in vast region of space (that can be much larger than the planet itself), called magnetosphere. For instance, the magnetosphere of Mercury has been also approximated with a dipole displaced from its centre[1]. 1.1.2 Induced currents When experiments showed that charge in motion produced magnetic fields that influenced the orientation of permanent magnets and magnetised needles, physicists started to wonder if the inverse effect could be possible, this is, that a magnetic field originated the motion of charges. During the year 1820, the galvanometer was already invented and measure the current through a wire was possible. Experiments where a conductor was placed near a permanent magnet were performed expecting that the field produced by the magnet induced a current in the conductor. If the configuration remained static no induced current was detected. Michael Faraday realised that is possible to induce a current but in order to achieve it the configuration should not remain static. Indeed he noticed that when the magnetic flux through the circuit composed by the conductor and galvanometer changed with time, a deflection in the needle of the instrument was observed. Let Γ represent a closed path that delimits a surface S, the flux through S is defined as.

(17) 1. Theoretical Background. Φ=. Z. B · da. 11. (1.23). S. Faraday induction law states that the line integral of E around Γ equals the rate of change of the flux Φ through S Z dΦ E · dl = − (1.24) dt Γ The integral in the left side of (1.24) is called the induced electromotance and is measured in volts [V]. 1.1.3 Inductance When two electrical circuits C1 and C2 are placed near one to the other, the magnetic flux density produced by the current flowing through C1 will contribute to the magnetic flux generated in C2 by the current flowing through it. The symmetric effect is also observed, the magnetic flux desity due to the current in C2 affects the flux produced by the induction field in C1 . It is possible to characterise quantitatively this magnetic coupling of circuits with the mutual inductance, Z Z dl1 · dl2 µ0 (1.25) M12 = 4π C1 C2 r and the magnetic flux produced by circuit C1 on C2 can be determined by Φ12 = M12 i1. (1.26). The integral in (1.25) is known in literature as Neumann integral and it is symmetric with respect to the subscripts as one would expect from the physics of the problem. Since every circuit is coupled to itself by its own magnetic flux it is natural to introduce the concept of self-inductance, Φ = Li. (1.27). where i is the current flowing through the circuit. The mutual and selfinductance depend only on the geometry of the circuits and are measured in henrys [H]. 1.1.4 Magnetism in matter Experiments have shown that three kinds of magnetic materials exist[23]. In two of these groups the interaction with magnetic fields is feeble even if they are in presence of strong magnetic fields, on the other hand in the remaining group the interaction is quite strong and can be perceived macroscopically without any special instrument. Paramagnetic and diamagnetic materials constitute the first two groups mentioned and the other one is constituted by ferromagnetic.

(18) 1. Theoretical Background. 12. materials. Magnetism is a quantum effect and henceforth quantum mechanics gives complete account of magnetism in matter, nonetheless classical models can be used to model macroscopically magnetised bodies. Ferromagnetic materials can be modeled as a distribution of magnetic dipole moments. The vector valued function M : S ⊂ R3 → R3 which gives the density of these dipoles at every point within the magnetised body is referred in literature as magnetisation. Using (1.21) it is possible to calculate the magnetic vector potential at any point outside the magnetised body. The magnetic potential of a body with magnetisation M is given by Z Z µ0 M×n 2 µ0 ∇×M 3 A= d r+ d r (1.28) 4π S r 4π V r where V denotes the volume occupied by the body and S the surface delimiting this volume. The potential given by (1.28) is composed of two terms and it is equivalent to the potential that a surface current density distributed on S with value λm = M × n. (1.29). and a volume current density in V given by Jm = ∇ × M. (1.30). would produce. A letter m was added as a subscript in expressions (1.29) and (1.30) to emphasise that these are not conduction currents as the ones present in conductors. The current given by (1.30) is the result of the magnetisation of the body. The magnetic flux density produced by this current (and by the surface current density) is obtained by taking the curl of (1.28). For a physical system where conduction currents and currents due to magnetisation are present, the magnetic flux density is calculated as the superposition of the respective fields produced by each current. If the nature of the different kind of currents is distinguished, (1.11) can be written as ∇ × B = µ0 (Je + Jm ) (1.31) where Je denotes the conduction currents. Using (1.30) results in ∇×(. B − M) = Je µ0. (1.32). The term in parenthesis in (1.32) is the vector valued function H : R3 → R3 defined as the magnetic field H=. B −M µ0. (1.33). The units of the magnetic field vector are ampere/meter. By rewriting (1.32) in terms of the magnetic field H one arrives to.

(19) 1. Theoretical Background. ∇ × H = Je. 13. (1.34). When there are no free currents and only currents resulting from magnetisation are present (1.34) reduces to ∇×H=0. (1.35). and the scalar magnetic potential φm : R3 → R proposed by Pierre Simon de Laplace can be introduced. The scalar magnetic potential is not uniquely determined just as the potential used in electrostatics. The magnetic field can be obtained by taking the gradient of the scalar potential. H = −∇φm. (1.36). Equation (1.36) is mathematically identical to the electrostatic situation where the electric field can be derived from a scalar potential by taking its gradient. Indeed if the analogy is continued the magnetic charge density ρm , can be determined by taking the divergence of H. It should be stressed that this is only a mathematical artifice that does not have physical reality. ρm = ∇ · H. (1.37). Even though the magnetic charge density is a mathematical artifice, it can serve to introduce and visualise a physical field of major importance in the design of fluxgate sensors, the demagnetising field HD : R3 → R3 . It is well known that a uniformly magnetised sphere can be modeled as a distribution of magnetic charge on its surface with one hemisphere containing charge of one sign and the remaing of opposite sign. The magnetic field in the interior of the sphere as a result of this charge distribution is uniform and it is directed in opposite direction to the magnetisation field. The field is said to try to oppose the magnetisation field and henceforth its name. In general the demagnetising field is not antiparallel to the magnetisation field, only for a uniformly ellipsoid (the sphere being a particular case) this is true. The relation between HD and M is given by HD = −DM. (1.38). where D is the demagnetising tensor. Linear Isotropic materials In general there is no physical law deduced from basic principles that relates the B and H fields within a magnetised sample, nonetheless for linear isotropic materials it is possible to introduce the magnetic susceptibility as M = χm H. (1.39).

(20) 1. Theoretical Background. 14. From (1.33) one can write B in terms of the magnetisation and the magnetic field B = µ0 (H + M). (1.40). by substituting (1.39) in (1.40) it is possible to relate the magnetic flux density and the magnetic field within an isotropic magnetic material B = µ0 µr H. (1.41). µr = 1 + χm. (1.42). where. is the relative permeability. It should be noticed from (1.39) and (1.42) that χm and µr are dimensionless constants. Ferromagnetic materials For ferromagnetic substances like iron, the constitutive relation among B and H is nonlinear and non single valued[14], B = B(H). (1.43). A typical curve showing the relationship given by (1.43) is depicted in Figure 1.1 and it is called hysteresis loop. The most important feature of the curve is the induction saturation Bsat , which plays a major role in the principle of operation of fluxgate magnetometers. The width of the hysteresis loop determines the magnetic hardness of the material (its sensitivity for past magnetisation history). It is common to use soft magnetic materials as cores in the fluxgate magnetometer for a robust operation and low losses. In specialised literature it is common to find an additional definition for the permeability of a material[19]. The differential permeability is defined as the slope at every point in the B − H curve, µd =. 1 dB µ0 dH. (1.44). 1.1.5 Boundary conditions In order to find a unique solution to Maxwell’s equations that describe the phenomena observed in a given region of study where two different media are present, a set of restrictions that the magnetic field (or the magnetic flux density) should satisfy in the boundary that delimits the interface must be imposed. These set of restrictions are referred as boundary conditions, (B2 − B1 ) · n = 0. (1.45).

(21) 1. Theoretical Background. 15. Relative permeability µ. r. 4. 4. x 10. 3 2 1 0 −1. −0.5. 0 Magnetic Field H[A/m]. 0.5. −1. −0.5 0 0.5 Magnetic Field H[A/m]. 1. 1. Magnetic Flux Density B[T]. −3. 1. x 10. 0.5 0 −0.5 −1 −2. −1.5. 1.5. 2. Fig. 1.1: The upper figure shows the dependance of the relative permeability µr with the magnetic field for a nonlinear material. In the bottom figure, a characteristic B − H curve is shown.. B2 × n =. µ2 B1 × n µ1. (1.46). where the subscripts denote the value of the magnetic flux density and permeability for each medium. Equation (1.45) states that the normal component of the magnetic flux density at the interface is continuous (in the mathematical sense) meanwhile the tangential component is discontinuous.. 1.2 A glimpse into the history of magnetism Those who cannot remember the past are condemned to repeat it George Santayana The following section is intended to serve as a brief reference to the history of magnetism. The history of magnetism is vast and interesting and is the result of the contribution of many illustrated as well as affortunated persons that were eager to unravel the misteries that Nature hides. Nowadays, most of the technological achievements that we enjoy in society are the result of the discoveries of those men that were able to understand the basic principles of an unified aspect of nature, electromagnetism. Magnets and their nature have been known for long time to the humankind. The latin poet Titus Lucretius Carus, author of De rerum natura poem, wrote about the inherent attraction of iron by magnets whose name derived from the.

(22) 1. Theoretical Background. 16. northern part of Greece, Magnesia[7]. The repulsion and attraction of magnets was an empirical fact also observed as well as the intrinsic aligment of lodestones with the geographical north-south direction. Besides their physical properties, magnets and lodestones were endowed with additional properties and powers with the pass of time. It was said that their usage restored husbands to wives, healed poissoned women, expelled demons and detected gold among other things. During the year of 1269, in the city of Lucera in Italy, Petrus Peregrinus de Maharncuria, wrote the Epistola Petri Peregrini de Maricourt ad Sygerum de Foucaucort, miltem, de magnete (Letter on the Magnet of Peter Peregrinus of Maricourt to Sygerius of Foucaucourt, Soldier) where he summarised all the knowldege about lodestones and instruments using them till then. It is important to mention that apparentely Petrus Peregrinus was the firts to shape a piece of lodestone into a sphere, since in the letter he explained how an iron needle interacted with it. Influenced by the dogmas of his epoch, he attributed the orientation of hanging lodestones to heaven forces in the celestial sphere. Almost 300 years had to pass for the human kind until another step was made toward the deciphering of the real nature of magnetism. In the year 1600, William Gilbert (1544-1603) the royal physician of Queen Elizabeth I published his De Magnete, Magneticisque Corporibus, et de Magno Magnete Tellure (On the Magnet: Magnetic Bodies Also, and On the Great Magnet the Earth). The importance of Gilbert’s treatise is that he confronted the superstitions in which magnetism was held captive. Gilbert correctly proposed in his work that the Earth itself behaves as a big magnet. Besides that, Gilbert carried on a series of experiments to study in a objective manner the magnetic field. Nonetheless Gilbert also believed that inside magnets there was a kind of essence that perturbated the dormant soul of iron with their presence in order to explain the attraction observed. During the eighteenth century admirable progress was done with electricity (phenomenon considered to be completely independent of magnetism). In this century Peter van Musschenbrok (1692-1761) discovered by accident that electricity can be stored, constructing the first capacitor in the Netherlands (Leyden jars). In France, Charles Augustus Coulumb (1736-1806) quantitatively proved the inverse square nature of electric attraction with his torsion balance. But the major achievement reached during this epoch is without hesitation the voltaic cell developed in Italy by Alessandro Volta (1745-1827) after having heard of the experiments conducted by Luigi Galvani (1737-1798) professor of anatomy at the University of Bologna and his assistant Giovanni Aldini. Advances in magnetism would have to wait until the nineteenth century when experiments with electricty were done in laboratories in Europe and in part of America. In Denmark, a man influenced by the famous german philosopher Immanuel Kant (1724-1804) started to wonder about the possibility of a relation between electricity and magnetism. His name was Hans Christian Ørsted (1777-1851) who during a public lecture placed a compass near a thin platinum wire that conducted electricity and noticed together with the audience the deflection of the compass. The results of this public demonstration.

(23) 1. Theoretical Background. 17. were presented in Paris on September 4, 1820 by Fran¸cois Arago (1786-1853) to the Academie des Sciences. They were really astonishing to the scientific community of that time since electricty and magnetism were considered two unrelated phenomena. Among the attendants to the meeting in Paris was Andre Marie Ampère (1775-1836). It is said that within a few weeks Ampère had reproduced Ørsted’s experiments and discovered that a short solenoid behaved in the same way as a permanent magnet near a current carrying wire; this encouraged him to affirm that magnetism was the result of circular currents flowing inside magnets. A major breakthrough in the progress of magnetism was acheived in England by Michael Faraday (1791-1867). Faraday’s life can be considered a Cinderella science tale, where the courage of a man to discover the basic principles behind nature made him succeed in life. Faraday lacked of a formal education in science, indeed it was by reading newspapers that he started to get in touch with it. He worked in a binding shop and one of his clients invited him to attend a series of lectures given by Sir Humphry Davy (1778-1829). Faraday took notes of the lectures and bind them. Davy was impressed by the work done by Faraday and this opened him the doors to work as an assistant in his laboratory. The lack of a mathematical training in Faraday’s formation was more than compensated with his ability to perform experiments and explain in a neat way the basic principles behind them. In the beginning Faraday was not attracted by magnetism; he started to work on it by an invitation of his friend Richard Phillips, editor of the Annals of Philosophy, who persuaded him to write an article about the knowldege by then of electricty and magnetism. It is well known the results of this work, Faraday performed a considerable number of experiments to undertsand the action at distance observed in interactions with magnets. One of the most remarkable contributions was the concept of field, a corner stone concept in the study of physics. Faraday also discovered the law of induction, giving the basis for a enormous stride in technological achievements (generation of electricity). The man who integrated all the knowledge in a single elegant frame was the Scotish scientist James Clerk Maxwell(1831-1879). Maxwell was the opposite to Farday, he was an excellent theoretical physicist. From his youth Maxwell started to show his mathematical skills when he correctly stated in an essay contest that the rings of Saturn could not be solid bodies. He gave solid mathematical fundation to Faraday’s ideas. When he formulated the equations of electromagnetism in mathematical terms he was able to predict theoretically a missing term in Ampère’s law, the displacemet current. This was one of the milestones in human history, when the basic laws of electromagnetism were started to be understood. The understanding of the electromagnetic field has lead us to enjoy the benefits of the generation of electicity, mobile communications, medical applications, exploration of space, among others achievements..

(24) 2. FLUXGATE MAGNETOMETERS. Fluxgate magnetometers are sensors designed to measure magnetic fields and date back to the beginning of the 1900’s. They were used during the Second World War to detect submarines and since the 1950’s, with the start of the space era, became an essential payload in space missions. Even though different configurations have been tested and reported in literature their principle of operation is basically the same and relatively simple. A nonlinear magnetic material is driven periodically into saturation by an alternating current inducing a voltage rich in even harmonics of the excitation current frequency when an external induction field is present. Using different signal processing techniques in the time or frequency domain, the information contained in the induced voltage can be used to determine the magnitude of the external field. The performance of the fluxgate magnetometers is considerably improved when volume compensation is used. With volume compensation a set of coils nullify the external magnetic field in the region where the fluxgate element is situated, eliminating nonlinearities in the sensor and allowing it to work as a zero level indicator. In this chapter the theory of operation of the single and double core geometries are described. The analysis is restricted to these geometries since the Small Magnetometer in Low-Mass Experiment (SMILE) sensor uses these configurations to measure magnetic fields.. 2.1 Principle of operation The basic components that constitute an elementary fluxgate magnetometer are a nonlinear magnetic material surrounded by a primary winding and a secondary set of coils used to measure the induced voltage. Through the winding flows a current that saturates periodically the magnetic core in both directions, inducing a voltage in the sensing coils as a result of the change in magnetic flux with time. The external magnetic flux through the sensing coils is gated with changes of the relative permeability due to the saturation of the core by the excitation current. In the absence of an external magnetic field the voltage induced will contain the odd harmonics of the excitation current frequency. If the sensor is placed in a region where a magnetic field is present the symmetry found previously will be broken and this will be manifested in the voltage induced which will contain harmonics not found in the driving signal. All the information about the external magnetic field can be deduced by analysing the harmonics added to the sensing coils signal. Different configurations have.

(25) 2. Fluxgate Magnetometers. 19. Fig. 2.1: In the single core configuration a nonlinear magnetic material is surrounded by excitation coils where the current that saturates the core periodically flows.. been proposed to measure magnetic fields. In specialised literature two groups are distinguished[20]: parallel and orthogonal sensors. With parallel sensors, the magnetic field produced by the excitation current is parallel to the external field component to be measured meanwhile with orthogonal geometries the direction of the aforementioned fields are perpendicular to each other. In this report the analysis is focused on a specific pair of parallel sensors: the single and double core sensor (also known as Vacquier sensor) since the SMILE magnetometer bases its principle of operation on them.. 2.2 Single core sensor The most basic fluxgate parallel detector is the single core sensor which consists of a nonlinear magnetic core (a soft magnetic material) surrounded by excitation and sensing coils. The geometry of the configuration is shown in Figure 2.1. The excitation coils are driven by an alternating current that produces a magnetic field which varies the permeability of the core an henceby saturates it. The core is saturated equally in both directions by the positive and negative excursions of the excitation current in the absence of an external magnetic field. With an external field the core remains saturated for a longer time in one direction than the other; unbalance shown in the harmonic content of the induced voltage. It is possible to determine the waveform of the induced voltage at the pickup coils for the single core geometry using Maxwell’s equations and further.

(26) 2. Fluxgate Magnetometers. 20. assumptions as it is shown in refereces [20] and [21]. The magnetic flux density inside the magnetic core is given by (1.40) which is writen again for convenience B = µ0 (H + M). (2.1). When the sensor is operating in the linear region (1.39) can be used to relate the magnetizing field and the magnetic field in the interior of the core. A demagnetising field will be present due to the geometry of the core (open ends), M = χm H. (2.2). HD = −DM. (2.3). H = Hext − DM. (2.4). In general the demagnetising and the magnetisation fields are not collinear vectors as was pointed out in Section 1.1.4. For the simplicity of the analysis it is found in the literature that the demagnetising factor is considered to be a scalar. From (2.1), (2.2), (2.3) and (2.4) it is possible to obtain M=. µr − 1 Hext 1 + D(µr − 1). (2.5). Hext 1 + D(µr − 1). (2.6). H= B = µ0. µr Hext 1 + D(µr − 1). (2.7). where the expression relating the magnetic susceptibility and the relative permeability has been used. From (2.6) it can be seen that the magnetic field inside the core differs from the external magnetic field, the demagnetising factor being responsible of this departure. Equation (2.7) is commonly rewritten as, B=. µr Bext 1 + D(µr − 1). (2.8). with the apparent permeability µa : R → R defined by[20], µa =. µr 1 + D(µr − 1). (2.9). The voltage induced in the sensing coils is obtained using Faraday’s law and can be approximated by dB (2.10) dt where the subscript sns was added to indicate parameters associated with the geometry of the pick-up (or sensing) coils; N denotes the number of turns Vsns = −Nsns Asns.

(27) 2. Fluxgate Magnetometers. 21. and A the area. The vector notation can be abandoned since all the vector fields considered are colinear. Using (2.8), (2.9) and (2.10) results in,. Vsns = −Nsns Asns. ". r (1 − D) dµ µr dBext dt Bext + 2 [1 + D(µr − 1)] 1 + D(µr − 1) dt. #. (2.11). Since fluxgate magnetometers measure slowly varying magnetic fields, the second term in (2.11) can be disregarded leading to r (1 − D) dµ dt Bext (2.12) [1 + D(µr − 1)]2 known in literature as the fluxgate equation [24]. The dynamics of the systems under consideration is completely specified by this equation.. Vsns = −Nsns Asns. 2.2.1 Induced voltage In order to determine the induced voltage in the pick-up coils it is necessary to solve (2.12) nonetheless, it is possible to estimate the shape of the induced signal by using simplified waveforms for the excitation field and further assumptions such as zero demagnetising factors. In Figure 2.2 an idealised B − H curve for a nonlinear magnetic core is shown; the behaviour of the magnetic material is described by a piecewise linear model and the demagnetising field ignored to simplify even more the analysis,   −Bsat Bsat H B (H) =  Hmax Bsat. if H ≤ −Hmax if −Hmax < H ≤ Hmax if H > Hmax. (2.13). Let us assume for the moment that no external field Hext is present and that the excitation current driving the core produces a field of the form,. Hexc (t) =.   . 4H0 τ t τ 0 − 4H τ (t − 4 ) + H0 4H0 3τ τ (t − 4 ) − H0. if t ≤ τ4 if τ4 < t ≤ if t > τ. 3τ 4. (2.14). where H0 and τ represent the maximum amplitude of the excitation field and the period of the driving signal respectively. Due to the nonlinear characteristics of the core, the external flux is gated out in the pick-up coils in a symmetric way since saturation is achieved in both directions during equal time intervals (see Figure 2.3). When an external magnetic field Hext is present, the field produced by the excitation current will have an offset which will unbalance the time intervals during which the core is saturated not being equal any more (see Figure 2.4). The induced voltage can be calculated from Farday’s law and using signal processing techniques in the time domain (correlation filter) or in the frequency domain (harmonic content analysis) the external field can determined..

(28) 2. Fluxgate Magnetometers. 22. 0.6. 0.4. B[T]. 0.2. 0. −0.2. −0.4. −0.6 −10. −5. 0 H[A/m]. 5. 10. Fig. 2.2: B − H curve of the core, the saturation is attained for Hmax = 10A/m with Bsat = 0.5T.. Excitation field[A/m]. 15 10 5 0 −5 −10 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 1 0.5 0 −0.5. Induced voltage(V). Magnetic flux density[T]. −15. −1. 2 1 0 −1 −2. Fig. 2.3: (Top) Excitation field with amplitude H0 = 15A/m and period τ = 125µs. (Centre) Magnetic flux density inside the magnetic material. No external magnetic field is present. (Bottom) Induced voltage in the sensing coil for N = 1 and Asns = 1 cm2 ..

(29) 2. Fluxgate Magnetometers. 23. Excitation field[A/m]. 15 10 5 0 −5 −10 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 1 0.5 0 −0.5. Induced voltage(V). Magnetic flux density[T]. −15. −1. 2 1 0 −1 −2. Fig. 2.4: (Top) Excitation field with amplitude H0 = 15A/m and period τ = 125µs. (Centre) Magnetic flux density inside the magnetic material. An external magnetic field Hext = 3A/m is present. (Bottom) Induced voltage in the sensing coil for N = 1 and Asns = 1 cm2 .. 2.3 Double core sensor The double core configuration or Vacquier sensor is composed of two cores wound by the same wire in opposite directions to ensure that the current flowing through the excitation coils will drive the cores into saturation in opposite directions in a symmetric way when no external field is present. The pick-up coils may surround each core independently or both of them. In the first case the sensing coils are connected in series. If the double core sensor is situated in a region where no external induction field is present the voltage measured in the pick-up coils will be zero, since the magnetic flux on each core cancels mutually. When an external magnetic field is present the flux will not cancel any more and a signal with twice the frequency of the excitation current is induced in the pick-up coils. Figure 2.5 shows one of the double core sensor that forms part of the SMILE sensor. 2.3.1 Induced voltage Like in the previous section, it is possible to illustrate the waveform induced in the pick-up coils by using piecewise linear models for the excitation fields; for an analytical approach of the Vacquier sensor refer to [19]. In the analysis of this configuration, identical cores are considered and are saturated in opposite directions by the same current. Assuming the same B −H characteristics for the single core sensor, the magnetic flux density inside the cores when no external field is present is plotted in Figure 2.6. It can be seen that it is necessary to break the symmetry with an external field to measure a non zero voltage in the.

(30) 2. Fluxgate Magnetometers. 24. Fig. 2.5: Double core sensor used in the SMILE magnetometer. The cores are two tapes of an amorphous alloy with dimension 16x1x0.02mm. The excitation coils are wound around each core and are connected in series (800turns). The pick-up coils surround the whole configuration.. pick-up coils (Figure 2.7).. 2.4 Volume Compensation The compensation technique is used to nullify the external magnetic field in the region where the fluxgate sensor is located. Two approaches have been proposed in literature to accomplish this[20]: component compensation and volume compensation. By using feedback coils, the external magnetic field component parallel to the axis of the fluxgate sensor can be cancelled. With this arrangement, each fluxgate element is equipped with its own set of feedback coils and the current used to drive them is extracted from the pick-up coils by means of a feedback loop. The main disadvantage of component compensation is that the fluxgate elements are sensitive to transverse magnetic fields, both present externally and those generated by compensating coils of other elements. This places strict requirements on mutual mechanical alignment and stability of the fluxgate elements. The effect of the external transverse fields cannot be removed, however. It is possible to avoid these drawbacks if the fluxgate elements are located in a common null field region with vector feedback. With this approach a set of compensation coils in three orthogonal directions create a field that will cancel completely an exterior field in the region where the fluxgate sensors are located. The advantage of this configuration is that the axes of of the sensor are defined by the orientation of the compensating coils only. For tri-axial sensors it.

(31) 2. Fluxgate Magnetometers. 25. H1[A/m]. 10 0 −10 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. B1[T]. 1. 0. −1. H2[A/m]. 10 0 −10. B2[T]. 1. 0. −1. Fig. 2.6: Excitation field driving the cores. Both cores are wound by the same coil in opposite direction and are consdiered to be identical, the excitation fields satisfy the relation Hexc1 = −Hexc2 . The period of the excitation field is τ = 125µs.. H1[A/m]. 10 0 −10 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. B1[T]. 1. 0. −1. H2[A/m]. 10 0 −10. B2[T]. 1. 0. −1. Fig. 2.7: Magnetic flux density inside the cores when an external magnetic field Hext = 3A/m is present..

(32) 2. Fluxgate Magnetometers. 26. 1. B [T]. 1 0 −1. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 0. 20. 40. 60 Time[µs]. 80. 100. 120. 2. B [T]. 1 0. Induced voltage(V). B1 + B2. −1 0.5 0 −0.5. 10 0 −10. Fig. 2.8: Induced voltage in the sensing coil for N = 1, Asns = 2 cm2 and τ = 125µs for the double core sensor.. has been found experimentally that a volume compensated sensor offers major advantages than independent coils for each axis [22]..

(33) 3. SIMULATION SOFTWARE. The software used to simulate the sensor elements that form part of the SMILE magnetometer was COMSOL Multiphysics. COMSOL Multiphysics is a commercial software for simulation of scientific and engineering models described by a system of partial differential equations using the Finite Element Method (FEM). In this chapter a brief introduction to COMSOL Multiphysics and the libraries used to model the components are given. For a detailed description of the software refer to [4].. 3.1 The COMSOL Multiphysics software COMSOL Multiphysics is an interactive software for modeling and simulation of scientific and engineering problems based on a partial differential equation formulation. The software contains a Graphic User Interface (GUI) that allows the user to create the geometry of the problem under study and additional tools to mesh the model and solve it. A set of predefined libraries are available and cover wide areas of physics such as: fluid mechanics, chemistry, photonics, quantum mechanics, electromagnetism among others. Once a model has been solved it is possible to visualise the solution obtained and if necessary use the post processing tools for further analysis of the data. It is also possible to create the geometry of the model with COMSOL Script, a programming language where the models can be saved in .M files and can be accesed with the interface to Matlab. 3.1.1 3D Electromagnetics Module All the models presented in this document were simulated with the Electromagnetics module of COMSOL v3.2 (AC/DC module in COMSOL v3.3). This is a specialised module containing a diversity of libraries designed to solve electromagnetic related problems. Since the purpose of the project is to model the components of the SMILE magnetometer for a static regime, the 3D Magnetostatic class was selected. This library solves simultaneoulsy the system of equations given by[5], −∇ · (−σv × (∇ × A) + σ∇V − Je ) = 0. (3.1). −1 e ∇ × (µ−1 0 µr ∇ × A − M) − σv × (∇ × A) + σ∇V = J. (3.2).

(34) 3. Simulation software. 28. where σ, v and V represent the electrical conductivity, the velocity of the medium and electric potential respectively. Since only magnetostatics simulations were performed, the system of equations is simplified to −1 e ∇ × (µ−1 0 µr ∇ × A − M) = J. (3.3). Boundary conditions Different boudary conditions can be specified in the surfaces that delimit the geometry of the models specified to be consistent with the physical situations simulated. The boundary conditions available in the quasistatic library are: magnetic field, surface current, electric insulation, magnetic potential, magnetic insulation and continuity among others. For a detailed description of the different boundary conditions in COMSOL, refer to [5]. The magnetic insulation, n×A=0. (3.4). n × (H1 − H2 ) = 0. (3.5). and the continuity condition,. were the most used ones in the simulations performed. 3.1.2 Mesh generation In order to find a solution for a given model in COMSOL it is necessary to mesh the geometry before applying the finite element method. Different parameters can be specified independently by the user such as the maximum element size of the elements, their growth rate, and resolution of curved regions. A set of mesh modes are included in COMSOL v3.2 and v3.3 with predefined values for these parameters which in most of the cases sucessfully mesh the given geometry. The majority of the models used to simulate parts of the SMILE magnetometer were meshed with the extremely coarse, specially those where the cores are included. With this mode the number of elements produced was below 50000, for which the memory of the computer was sufficient to obtain a solution. Figure 3.1 shows a meshed model of the SMILE magnetometer including all its components. Due to the relative dimensions of the cores compared with other parts of the sensor, the number of elements increases considerably even though the coarsest mesh mode is used. 3.1.3 Solvers COMSOL Multiphysics includes a variety of solvers appropiate to the nature of the problem analysed. The results reported in this document were obtained with time independent linear and nonlinear solvers. The solver used was SPOOLES1 . COMSOL v3.3 includes a solver that uses efficiently double processor architectures and saves considerable time of computation, the PARDISO solver. 1. Sparse Object Oriented Linear Equations Solver..

(35) 3. Simulation software. 29. Fig. 3.1: Mesh of the SMILE model in COMSOL. The geometry has been decomposed in a total of 68873 elements.. 3.2 Finite Element Method The finite element method is a numerical method used to solve partial differential equations. It started to be applied in electromagnetic problems during the 1960s and it is well suited to analyse complex geometries[25]. The basic ideas behind courtaines are[25]: (i) discretise the analysed geometry in elements, (ii) define a set of basis functions over each element, (iii) construct an approximate solution using the basis functions, (iv) apply Garlekin’s method to obtain a system of equations and (v) solve the resultant system to determine the contribution of each basis function to the approximation[18]. The essence of the method is exemplified by applying it to a 1D problem. For this purpose, the Poisson equation will be solved in the unit interval [0, 1], −∇2 Φ =. ρ ǫ. (3.6). with boundary conditions Φ(0) = Φ(1) = 0. (3.7). Mesh It is necessary to discretise first the geometry where the solution is computed. For the case considered the simplest discretization is a uniform partition of the unit interval in n segments, being the norm of the partition h,.

(36) 3. Simulation software. 30. 1. φ2(x). 0.8 0.6 0.4 0.2 0. 0. 0.2. 0.4. 0.6. 0.8. 1. 0.6. 0.8. 1. x 6. 2. dφ /dx. 4 2 0 −2 −4 −6. 0. 0.2. 0.4 x. Fig. 3.2: Base function and its derivative for k = 2, when the geometry has been meshed with n = 5 elements.. 1 n and each node in the mesh labeled by xk , k = 0, 1, . . . , n. h=. (3.8). Basis Functions The set of basis functions used in the FEM are well localised2 . The most common are polynomial approximations of the form,. ϕk (x) =.  0    x−xk−1   . h xk+1 −x h. 0. if 0 ≤ x < xk−1 if xk−1 ≤ x < xk if xk ≤ x < xk+1 if xk+1 ≤ x ≤ 1. (3.9). Equation (3.9) defines a set of well localised functions in the unit interval. Figure 3.2 shows the second base function and its derivative when the geometry is divided in 5 subintervals (elements). Approximate solution The approximate solution to (3.6) is constructed with the basis functions {ϕk (x)} as 2. Compact support.

(37) 3. Simulation software. ˜ Φ(x) =. X. 31. n−1 k=1 ak ϕk (x). (3.10). where the coefficients ak are determined using the Garlekin method. Garlekin method The system of equations that the ak coefficients should satisfy is obtained by substituting (3.10) in (3.11) Z. 1. 0. ǫ. ˜ dϕk dΦ dx = dx dx. Z. 1. (3.11). ρϕk dx. 0. for k = 1, 2, . . . , n − 1. Equation (3.11) is a particular case of the most general integral used in Garlekin’s method. The Garlekin method is a technique developed by the russian mathematican Boris Grigoryevich Garlekin to determine the weighting coefficients of an expansion of the form ˜ Φ(x) = a1 ϕ1 (x) + a2 ϕ2 (x) + . . . + an ϕn (x). (3.12). for an approximate solution to the differential equation LΦ = f. (3.13). such that the error of the approximation is orthogonal to the basis functions, where L is a differential operator and f is the inhomogeneous term. The system of equations resulting from (3.11) is of the form Ma = f. (3.14). with .   1  M=  h  . 2 −1 0 . . . 0 −1 2 −1 . . . 0 0 −1 2 . . . 0 0 0 −1 . . . 0 .. .. .. .. .. . . . . . 0 0 0 . . . −1 . and.   a= . a1 a2 .. . an−1.     .         . (3.15). (3.16).

(38) 3. Simulation software. 32. 0.14 N=5 N = 10 Theoretical. 0.12. 0.1. Φ(x). 0.08. 0.06. 0.04. 0.02. 0. 0. 0.2. 0.4. 0.6. 0.8. 1. x. Fig. 3.3: Solution of Poisson equation in the unit interval [0, 1] using the FEM for ρ = 1 and ǫ = 1. The approximate solutions were calculated using 5 and 10 mesh elements of the unit interval.. . ρ  f =h  ǫ. 1 2 .. . n−1.     . (3.17). It can be seen that M is a sparse matrix and the elements in the diagonal are non vanishing. This is a characteristic feature of the FEM, the matrices obtained are tridiagonal. The advantage of sparse matrices is that efficient algorithms such as Gauss elimination can be implemented in computer programs to obtain an inverse. In Figure 3.3 the approximate and analytical solutions for (3.6) are shown. For the numerical approximation 5 and 10 elements were used to mesh the unit interval. As the number of elements increases the discrepancy between the real and the approximated solution is reduced. In general, the method previously outlined is extended to two and three dimensions. For 2D problems the mesh elements consist of triangular or quadrilateral regions and for 3D cases the geometry is discretised with tetrahedral elements. It is important to mention that for the 1D situation described above, the FEM reduced to the method of Finite Differences (FDM). In general these methods are not equal, the FEM is used to find approximations to solutions of differential equations meanwhile the FDM is an approximation to the differential equation itself [9]..

(39) 3. Simulation software. 33. 3.2.1 Distinctive Features The FEM has proven to be a powerful technique to solve not only EM problems. Among its most distinctive features are that it can be used to analyse complex geometries taking into account inhomogeneties of the medium. Besides that, the resultant system of equations obtained is characterised by sparse matrices with most of their entries equal to zero allowing the use efficient computer methods to solve them. Like every numerical method the FEM has also disadvantages. The most critical probably being the mesh of the geometry even though considerable advances have been achieved nowadays with computer programs in this area. A direct consequence of the discretization of the geometry is the dimensions of the matrices used to calculate the coeffiecients of the basis functions which can surpass the memory required to find a solution..

(40) 4. ACCURACY ESTIMATION. In order to have an estimate of the accuracy of the solutions obtained with the numerical approximation, different simulations were performed with a model whose theoretical solution is well known. The model used is a sphere of permeable material immersed in a region pervaded by a uniform induction field. The results obtained indicate that the accuracy of the numerical approximations depends on the number of elements used in the mesh as well as the extent of region where the solution is computed.. 4.1 Sphere of permeable material inmersed in a uniform magnetic flux density The model used as reference to compare the theoretical solution with the numerical approximations was a sphere of relative permeability µr = 4 × 104 , radius of 1cm situated in a region where a uniform magnetic flux density of 1mT is present. The geometry of the problem is depicted in Figure 4.1.. 4.2 Theoretical solution It is possible to find a closed expression that characterises the magnetic flux density in the interior and exterior of the sphere shown in Figure 4.1. As can be noticed from the description of the problem, no free currents are present, therefore it is possible to express the solution via the scalar potential introduced in section 1.1.4. The symmetry of the configuration allows to use spherical coordinates (r, θ, φ). The magnetic scalar potential should satisfy Laplace’s equation and be of the form φim =. ∞ X. αn rn Pn (cosθ). (4.1). n=0. for points situated in the interior of the sphere and φem = −H0 rcosθ +. ∞ X. βn r−(n+1) Pn (cosθ). (4.2). n=0. for those situated in the exterior; where a and H0 denote the radius of the sphere and the external magnetic field respectively. The coefficients αn and βn can be determined by applying boundary conditions at r = a,.

(41) 4. Accuracy estimation. 35. Fig. 4.1: Sphere of permeable material inmersed in an uniform magnetic flux density of 1mT. The radius of the sphere is 1cm and µr = 4 × 104 .. ∂φem ∂φim (r = a) = (r = a) ∂θ ∂θ. (4.3). ∂φem ∂φi (r = a) = µ0 µr m (r = a) (4.4) ∂r ∂r which lead to the system of equations for the only two non-vanishing coefficients α1 and β1 µ0. β1 − aα1 = aH0 a2. (4.5). 2β1 = −H0 (4.6) a3 By solving the linear system given by (4.5) and (4.6) the magnetic scalar potential can be expressed in a closed form as µr α1 +. φim = −. 3H0 rcosθ µr + 2. φem = −H0 rcosθ + a3 H0. µr − 1 cosθ µr + 2 r2. (4.7) (4.8).

(42) 4. Accuracy estimation. 36. The components of the magnetic flux density and the magnetic field can be derived from (4.7) and (4.8). For the magnetic flux density, its components in cartesian coordinates are, Bxi = 0. (4.9). Byi = 0. (4.10). Bzi = 3 Bxe = 3a3 B0. µr B0 µr + 2. µr − 1 xz 2 µr + 2 (x + y 2 + z 2 ) 52. µr − 1 yz µr + 2 (x2 + y 2 + z 2 ) 52 " # 2 2 2 e 3 µr − 1 x + y − 2z Bz = 1 − a B0 µr + 2 (x2 + y 2 + z 2 ) 52 Bye = 3a3 B0. (4.11) (4.12) (4.13). (4.14). 4.3 Numerical Solution The numerical simulations were implemented using COMSOL Multiphysics v.3.2. With the results obtained it was realised that the accuracy of the solution depends on the number of mesh elements and the dimensions of the region where the model is meshed which we shall refer henceforward as the simulation box. The results improved as the dimensions of the simulation box were considerably greater than the radius of the sphere and a finer mesh was used, however this increased the time for computation and memory used. Figure 4.2 to 4.4 show the solution obtained with the simulation software and the theoretical one in a 30x30x30cm region. The magnitude of the magnetic flux density is plotted at slices parallel to the xy plane at three different heights. The plots in the upper part of the figures show the magnitude of the magnetic flux density, the calculated and the one obtained numerically. In the lower part, the absolute and percentual error are plotted. As can be seen from the aforementioned figures, with these parameters the absolute error is at least two orders of magnitude less than that of the solutions. It can be noticed also, that discrepancy of the solutions tends to increase as the slice used to compare them approaches a height of z = 0cm and decreasess for planes situated far from the sphere. In order to quantify the deviation of the numerical solution with the theoretical one, an estimated error using the norm of the magnetic flux density was calculated in a fixed region τ for different mesh modes and simulation boxes. With |BC |, |BT | : τ ⊂ R3 → R denoting the norm of the magnetic flux density of the numerical and theoretical solutions respectively in the τ region, the estimated error ε : τ ⊂ R3 → R was calculated as.

(43) 4. Accuracy estimation. 37. Slice at z =−0.045m. Slice at z =−0.045m. −3. −3. x 10. x 10. 1.04. 1.03. |[T]. 1.02 Comsol. 1.02 1.01. 1.01. |B. |B. Theory. |[T]. 1.03. 1. 1. 0.05. 0.05 0.05 0 y[m]. 0.05 0. 0 −0.05. −0.05. y[m]. x[m]. Slice at z =−0.045m. 0 −0.05. −0.05. x[m]. Slice at z =−0.045m. −6. x 10. ||*100%. 0.4 0.3. Theory. 3. Comsol. |/|B. 2 1. |1−|B. ||B. Theory. |−|B. Comsol. ||[T]. 4. 0 0.05. 0.2 0.1. 0 0.05. 0.05 0 y[m]. 0 −0.05. −0.05. x[m]. 0.05 0 y[m]. 0 −0.05. −0.05. x[m]. Fig. 4.2: Comparison of the numerical and theoretical solution at z = -45mm. The length of the simulation box is 30cm and the SPOOLES solver together with a fine mesh were used to obtain the numerical solution..

(44) 4. Accuracy estimation. 38. Slice at z =−0.02m. Slice at z =−0.02m. −3. −3. x 10. x 10. 1.15. 1.05. Comsol. 1.2. 1.15. |[T]. 1.25. 1.2. |B. |B. Theory. |[T]. 1.25. 1.1 1.05. 1.1. 1. 1. 0.95 0.05. 0.95 0.05 0.05 0 y[m]. 0.05 0. 0 −0.05. −0.05. y[m]. x[m]. Slice at z =−0.02m. 0 −0.05. −0.05. x[m]. Slice at z =−0.02m. −5. x 10. ||*100%. 2. Theory. 1.5. Comsol. |/|B. 1 0.5. |1−|B. ||B. Theory. |−|B. Comsol. ||[T]. 2 1.5. 0 0.05. 1 0.5. 0 0.05. 0.05 0 y[m]. 0 −0.05. −0.05. x[m]. 0.05 0 y[m]. 0 −0.05. −0.05. x[m]. Fig. 4.3: Comparison of the numerical and theoretical solution at z = -20mm. The length of the simulation box is 30cm and the SPOOLES solver together with a fine mesh were used to obtain the numerical solution..

(45) 4. Accuracy estimation. 39. Slice at z =0m. Slice at z =0m. −3. −3. x 10. x 10. 3. 4. 2.5 3 |BComsol|[T]. Theory. |[T]. 2 1.5. |B. 1. 2 1. 0.5 0 0.05. 0 0.05 0.05 0 y[m]. 0.05 0. 0 −0.05. −0.05. y[m]. x[m]. 0 −0.05. −0.05. x[m]. Slice at z =0m. Slice at z =0m −5. x 10. ||*100%. 12. 8. 4. |/|B. 3 2 1. 0 0.05. 10. Theory. 4. Comsol. 5. |1−|B. ||B. Theory. |−|B. Comsol. ||[T]. 6. 6. 2. 0 0.05. 0.05 0 y[m]. 0 −0.05. −0.05. x[m]. 0.05 0 y[m]. 0 −0.05. −0.05. x[m]. Fig. 4.4: Comparison of the numerical and theoretical solution at z = 0mm. The length of the simulation box is 30cm and the SPOOLES solver together with a fine mesh were used to obtain the numerical solution..

(46) 4. Accuracy estimation. 40. 0.09 normal coarse coarser extra coarse extremely coarse. 0.08. Relative error. 0.07 0.06 0.05 0.04 0.03 0.02 0.01 10. 15. 20. 25 30 35 Length of simulation box [cm]. 40. 45. 50. Fig. 4.5: Estimated error for the numerical solutions obtained using different meshing modes.. R R R

(47)

(48) τ

(49) 1 − RRR ε=. τ.

(50). |BC |

(51) |BT |

(52) dxdydz. dxdydz. (4.15). Figure 4.5 summarises the results of the simulations performed with different lengths of simulation boxes and mesh modes. The dimensions of the τ region were kept constant for all the cases shown to 5x5x5cm. As can be seen from the figure the results improve as the simulation box increases its size and a finer mesh is used..

No results found