Magnetic Coupling Visualization Using Sensitivity to Magnetic Permeability

This letter presents a method for visualizing magnetic coupling based on its sensitivity to magnetic permeability. Because a magnetic flux causes magnetic coupling that is proportional to permeability, an increase in permeability along the propagation path strengthens the magnetic coupling, whereas an increase in permeability outside the path has little effect on the coupling. The propagation path for magnetic coupling can thus be visualized based on the sensitivity to permeability. Calculating this sensitivity is faster than performing an electromagnetic simulation using the adjoint variable method. Two examples, namely, a loop pair and a simple electromagnetic interference filter, are used to show that the propagation path due to magnetic coupling can be visualized using the proposed method.


Letters
Magnetic Coupling Visualization Using Sensitivity to Magnetic Permeability Katsuya Nomura , Member, IEEE Abstract-This letter presents a method for visualizing magnetic coupling based on its sensitivity to magnetic permeability.Because a magnetic flux causes magnetic coupling that is proportional to permeability, an increase in permeability along the propagation path strengthens the magnetic coupling, whereas an increase in permeability outside the path has little effect on the coupling.The propagation path for magnetic coupling can thus be visualized based on the sensitivity to permeability.Calculating this sensitivity is faster than performing an electromagnetic simulation using the adjoint variable method.Two examples, namely, a loop pair and a simple electromagnetic interference filter, are used to show that the propagation path due to magnetic coupling can be visualized using the proposed method.

I. INTRODUCTION
T HE identification of dominant parasitic coupling is impor- tant for electromagnetic compatibility (EMC) design [1], [2].In magnetic field coupling [3], which is a major type of coupling, the magnetic flux from a magnetic field source flows along a propagation path to the receiving point and generates an induced current by penetrating the parasitic receiving loop.Unlike conducted noise, noise due to magnetic coupling is not transmitted along conductors, making it difficult to optimize the design of a system in terms of EMC.
An electromagnetic interference (EMI) scanning system [4], [5] or an electromagnetic field analysis [6] is commonly used to identify the locations of strong magnetic fields.However, a strong magnetic field is not always the source of the radiation that causes EMC problems.In addition, these methods cannot identify the propagation path and the receiving point.
Methods for identifying and visualizing coupling paths using energy parcels [7], [8] or impedance density [9], [10], [11] have been developed.Although these methods can directly visualize energy propagation or mutual coupling, they use a physical quantity defined by the product of the electric field and the magnetic field, making it difficult to distinguish between magnetic field  [12].Thus, it is important for EMC design to distinguish magnetic field coupling from electric field coupling.This letter presents a method for visualizing magnetic coupling based on the sensitivity to magnetic permeability.The rest of this letter is organized as follows.The principle of the visualization method and a fast sensitivity analysis method are explained in Section II.Two examples used for verifying the proposed method are given in Section III.Finally, Section IV concludes this article.

II. METHOD
Noise propagation characteristics are evaluated using the S-parameter S 21 , whose value depends on various factors, including magnetic field coupling.If magnetic field coupling is dominant, the magnet flux, which is proportional to permeability, along the propagation path greatly affects S 21 ; an increase in permeability along the propagation path strengthens the coupling and increases |S 21 |, whereas an increase in permeability outside the path has little effect on |S 21 |.Therefore, magnetic field coupling can be visualized using the partial derivative of |S 21 | to magnetic permeability, which is called the sensitivity.
The sensitivity can be written as where μ r_i is the relative permeability of the ith element.This expression implies that the sensitivity can be obtained numerically and that the analysis domain can be divided into many elements by discretization.The finite element method (FEM) is used here as the numerical method.When calculating the sensitivity using the difference method, the number of required forward analyses almost equals the number of elements.In contrast, when calculating the sensitivity using the adjoint variable method, only one forward analysis is required.In the adjoint variable method, the extended objective function J * is defined as where J is the objective function, which corresponds to |S 21 | in this case, K is the stiffness matrix derived from the FEM, e is the electric field vector, f is the nodal force vector, and λ is the adjoint variable vector.The sensitivity can be written as By solving the adjoint equation the calculation of ∂e/∂μ r_i can be omitted and the sensitivity can be evaluated using only the second term, which greatly reduces the computational cost.Sensitivity analysis using the adjoint variable method is a very well known technique in topology optimization [13], which is a structural optimization method that iteratively updates the structure based on sensitivity information.Topology optimization, originally developed for structural design, has been applied to electromagnetic design problems and optimization using the FEM [14], the finite-difference time-domain (FDTD) method [15], and the method of moments [16].Although the sensitivity can be obtained using these three methods, only the FEM and the FDTD method are suitable for magnetic path visualization in the proposed method because sensitivity information regarding the air domain around the target structure is important.
In practice, the value of sensitivity in (1) depends on the volume of the element; a larger element leads to a larger value.Therefore, the following normalized sensitivity is evaluated for visualization: where V i is the volume of the ith element.

III. NUMERICAL EXAMPLES
The proposed method is verified using two examples, namely a loop pair and an EMI filter.The loop pair is an example of pure magnetic field coupling, while the EMI filter is an example of inductive noise and conductive noise mixed together.COMSOL multiphysics and its RF and optimization modules were used for the electromagnetic and sensitivity analyses.All computation was conducted by COMSOL by following three steps.
1) Discretize the weak form of the wave equation derived from Maxwell's equations, and solve the matrix equation, that is, Ke = f .This step determines the electric field vector e. 2) Solve the adjoint (4).This step determines the adjoint variable vector λ. 3) Calculate the sensitivity ∂J ∂μ r_i as λ T ∂K ∂μ r_i e − ∂f ∂μ r_i .This step determines the normalized sensitivity 1 The sensitivity calculation was conducted using symbolic partial differentiation and the chain rule as stated in the COMSOL user's guide [17].

A. Loop Pair
The loop pair structure is shown in Fig. 1.In an air with dimensions of 2000, 2000, and 500 mm, two torus-shaped metal loops with a major radius of 100 mm and a minor radius of 1 mm are set with a separation distance of 1000 mm.Ports 1 and 2 are set as the left and right loops, respectively.The left loop is excited at 10 MHz and the normalized sensitivity in ( 5) is calculated.At this frequency, the wavelength is almost 30 m and the loops are thus regarded as small loop antennas.
Electromagnetic and sensitivity analyses were conducted on a computer with a dual-core CPU (Xeon 4214R) and 256 GB of RAM.The computational times were 16 and 10 s, respectively.As shown in Fig. 2, although the flux density amplitude is high only around the left loop, which is the magnetic source, the normalized sensitivity is high around both loops and between the loops, as shown in Fig. 3.This confirms that propagation due to magnetic coupling can be visualized using the sensitivity.

B. EMI Filter
The proposed method is applied to a simple EMI filter.The circuit diagram and structure of the filter are shown in Figs. 4  and 5, respectively.In this C-L-C-L-C filter circuit, the dominant noise depends on the circuit parameters [18].When the inductances of L1 and L2 are small, the conduction noise directly flows from port 1 to port 2 via L1 and L2; magnetic coupling is negligible in this case.On the other hand, when the inductances of L1 and L2 are sufficiently large, the conduction noise is suppressed and instead the induction current due to magnetic coupling between the input and output loops becomes dominant.Fig. 6(a) and (b) shows the noise propagation path for these two cases, respectively.
The situations corresponding to these two cases were realized by a simulation at 10 MHz with the circuit parameters given in Table I.First, to confirm that the dominant paths are different Authorized licensed use limited to the terms of the applicable license agreement with IEEE.Restrictions apply.A sensitivity analysis was conducted at 10 MHz under the two parameter settings given in Table I.For both settings, the computational times for the electromagnetic and sensitivity analyses were 37 and 21 s, respectively, on the same computer used in the loop pair example.Figs.7 and 8 show the sensitivity for settings 1 and 2, respectively.As shown in Fig. 7(a) and (b), when conductive noise is dominant, the sensitivity is greater above the conductors connected to the capacitors.This implies that increasing the permeability in this region increases the parasitic inductance, which degrades the bypass performance and increases |S 21 |.On the other hand, as shown in Fig. 8(a) and (b), when inductive noise is dominant, a region with high sensitivity appears between the input and output loops; the magnetic coupling can thus be visualized.Figs.7(c) and 8(c) support these arguments.In Fig. 7(c), the sensitivity is high above the three patterns connected to the capacitor.In Fig. 8(c), the sensitivity is high inside the input and output loops but low outside the loops; this indicates that higher permeability inside the loops strengthens the magnetic field coupling, whereas higher permeability outside the loops causes the magnetic flux to deviate from the loop.A negative sensitivity means that S 21 decreases as permeability increases, as can be explained from the sensitivity definition.
Thus, the sensitivity distribution strongly depends on the noise propagation path, which can help identify the magnetic coupling.

IV. CONCLUSION AND FUTURE WORK
In this letter, a method for visualizing magnetic coupling was proposed and verified.The method uses the sensitivity distribution of |S 21 | with respect to magnetic permeability to evaluate the effect of magnetic flux on the transmission characteristics.Two simple examples, namely, a loop pair and a C-L-C-L-C EMI filter, confirmed the effectiveness of the proposed method.In the future, this method will be applied to more realistic structures.It may be possible to visualize the electric field coupling by using the sensitivity of the dielectric permittivity.This is a topic for future work.Reflecting the results of the visualization for product design is also a subject for future work.

Fig. 3 .
Fig. 3. Normalized sensitivity (a) in analytic domain and (b) on plane.In (a), only domains with values greater than 0.001 are plotted.

Fig. 5 .
Fig. 5. Structure of simple EMI filter.TABLE I CIRCUIT PARAMETERS FOR SIMPLE EMI FILTER

Fig. 7 .
Fig. 7. Normalized sensitivity (a) in analytic domain, (b) on x-y plane, and (c) on y-z plane for setting 1.In (a) and (c), only values and absolute values greater than 0.001 are plotted, respectively.

Fig. 8 .
Fig. 8. Normalized sensitivity (a) in analytic domain, (b) on x-y plane, and (c) on y-z plane for setting 2. In (a) and (c), only values and absolute values greater than 0.001 are plotted, respectively.

TABLE II |S
21 | VALUES UNDER ORIGINAL AND DOUBLED INDUCTANCESAuthorized licensed use limited to the terms of the applicable with IEEE.Restrictions apply.