Double Magnetic Loop and Methods for Calculating Its Inductance

1Traffic Control Systems Group, ITACA Institute, Universitat Politècnica de València, Valencia 46022, Spain 2Department of Electronic Engineering, ITACA Institute, Universitat Politècnica de València, Valencia 46022, Spain 3Chemical and Nuclear Engineering Department, Institute of Industrial, Radiological and Environmental Safety, Universitat Politècnica de València, Valencia 46022, Spain 4Electronic and Telecommunications Engineering Department, Universidad Politécnica de Tulancingo, Tulancingo 43629, Mexico


Introduction
Magnetic loops are the most common sensors on roads around the world since they are an affordable and highly developed technology with a simple operation that is not affected by environmental conditions [1][2][3][4][5][6][7].Although these ones imply to drill and work on the road for their installation and possible future repairs like the rest of intrusive sensors [8], in practice, magnetic loops still have a long future ahead.Even though they might seem outdated, these are actually a widely extended and well-known reliable technology that offers good performance at a low price.Proof of this is that today they continue to be installed on the roads and they are even fundamental elements in the new algorithms for traffic management [9][10][11][12][13][14][15].
Their operation is straightforward, since it is based on the impedance variation that is recorded in the magnetic loops during the passage of vehicles over them, and as shown in Figure 1, an entire system usually consists of three parts [16]: (I) A magnetic loop formed by a wire with one or more turns superficially buried in the pavement.(II) A cable that links the magnetic loop with the control booth, which is also buried in the pavement.(III) An electronic unit located in the control booth that contains an oscillator and amplifiers to excite the inductive loop.
In order to have a better understanding of how they work, there are many publications and bibliography [1] since they are one of the most widespread sensors.However, a brief physical explanation is provided in the following points: (i) The electronic unit together with the magnetic loop forms an oscillator circuit.(ii) The current which passes through the loop produces a magnetic field  →  around the cable as shown in (1), where  is the number of turns of the loop,  is the current expressed in Amperes, and  is the length of the loop expressed in meters.
(iii) This magnetic field  →  produces a magnetic flux 0 through the loop as shown in (2), where  →  is the magnetic flux density,  →  is the surface enclosed by the loop,   is the relative magnetic permeability of the medium, and   is a constant value (4 ⋅ 10 −7 N / A 2 ).
(iv) The result is that the inductance of a common single loop  expressed in Henrys is obtained as follows: In this way, when a vehicle or any object built with ferromagnetic materials passes through the magnetic field generated by a magnetic loop buried on the road with a surface area  →  , a number of turns , and a current intensity  as shown in Figure 2, there is a decrease in the global magnetic field because of the currents that are induced in the vehicle.As seen in (3), the loop inductance is proportional to the magnetic flux, which causes that when passing a vehicle over it, the inductance also decreases.Moreover, like any oscillator circuit, the oscillation frequency of the whole system will be given by where  is a constant that depends on the characteristics of the electronic components used in the construction of the oscillator circuit.Thus, when a vehicle passes over a loop, we can obtain what is commonly known as "the vehicle magnetic profile" or "the vehicle inductance signature" by analyzing the inductance or frequency variation recorded.This magnetic profile is different for each type of vehicle as seen in Figure 3, which allows classifying the different vehicles as motorcycles, cars, trucks and buses.However, while the vehicle magnetic profiles for the single loops are widely known, the magnetic profiles left by the passage of vehicles over the double loops have not yet been studied, although we can anticipate that these new magnetic profiles will offer much more information than the previous ones.
To estimate the vehicle speed and classify them, nowadays it is necessary to use two single loops, since a single one is not able to get all the necessary parameters to do it.After analyzing the magnetic profile, there are two unknown data (vehicle length and vehicle speed) with the variation of a single parameter (inductance or oscillation frequency).
For that reason, there are usually two loops per lane separated by a certain distance.In this manner, the passage of a vehicle over the first loop is recorded in the detector, and after a short interval of time, the vehicle passes again over the second loop where it is also recorded [8,17].Then, as the distance between both loops is known by design, the vehicle speed, the direction of traffic and the vehicle length can be finally estimated, as well as the vehicle axle detection [14].Nevertheless, with the use of the double loops this problem would be solved and it would only be necessary to use one loop to find out all the previous data, since they have a simpler, more compact and more economical electronics.Moreover, having a single signal instead of two would facilitate the implementation of the measurement system.
Therefore, our work will aim to present and describe the characteristics of the double loop, to offer different methods to calculate its inductance, to verify which one provides better results and to improve the functional characteristics of the popular single loops, which despite being the most installed sensor on the roads, they are actually only dedicated to count vehicles.The presentation of the new vehicle magnetic profiles, the parameters that can be extracted with them and the advantages offered over the conventional loop will be the subject of the following paper.

Electromagnetic Analysis
The design, shape, and construction of a rectangular or circular single loop are well-known worldwide [3].However, although every day there are more algorithms and new systems to regulate traffic based on magnetic loops [18,19], the double loop capable of providing better performance than the current loops at lower cost is still unknown.
A double loop is no more than the union of two rectangular loops, which can have different dimensions and turns (not to be confused with two single loops spaced at a certain distance as described above).How to implement this type of loop can be quite varied, but a classic way to proceed would be to build the outer loop and then a smaller inner loop located at one extreme.However, another way to build it could be to construct a small loop and then put another small one next to the first one.In all cases, the direction of the current in each loop can be chosen with the aim of generating different types of configurations.At any rate, in our study we will present a general theoretical analysis capable of simulating any type of design.Therefore, in order to analyze this type of loop, the space will be divided into three sections as shown in Figure 4. (i) The first section,  1 , corresponds to the three red segments located in the plane of the negative values of .Two segments parallel to the -axis with a length of  and one segment parallel to the -axis with a length of 2.
(ii) The second section,  2 , corresponds to the three turquoise segments located in the plane of the positive values of .Two segments parallel to the -axis with a length of  and one segment parallel to the -axis with a length of 2.
(iii) The third section,  3 , corresponds to the blue segment located on the -axis at  = 0 which has a length of 2.
In this way, point  where the magnetic field will be calculated will be determined by its coordinates (, , ).Consequently and as shown in Figure 4, the distances from any point of the double loop (  ,   ,   ) to the point of analysis of the magnetic field (, , ) will be defined as follows: In order to perform the electromagnetic analysis, the starting point will be the physical phenomenon of magnetic field generation due to the electrical currents flowing through a conductor.Maxwell's equations [20] revealed that the divergence of  →  is zero: This indicates that  →  has the solenoidal property (with no divergence), which means that the magnetic field  →  can be represented using an auxiliary vector function  →  as follows: And thus The auxiliary vector  → A is called magnetic potential vector [21] and it is related to the sources of the stable current density J which are responsible for generating the magnetic field.Then, in the case of a linear conductor,  → A is given by where  is the current in the linear conductor and  is the distance from the conductor to the analysis point.However, (9) represents a solution for the case of thin conductors, but the general solution must include a volumetric integral.In this way, the magnetic potential vector for each segment of our inductive double loop is obtained as follows: In the previous expressions (  ) it must be noted that, as shown in Figure 4, the subindex  refers to the considered segment of the section , which in turn corresponds to the -axis.For example,  →  32 would mean the third segment of section  2 located on the -axis, and the values ,  and  would correspond to the dimensions of the loop described above.However, before proceeding with the electromagnetic analysis of the double loop, we will define a series of terms to simplify the previous equations.These terms will be In this way, combining (10a), (10b), and (10c) with (11), we can obtain the magnetics potentials vectors in a simpler way as seen in the following: For the section  1 , located in the plane of the negative values of , the magnetic potential vector would have two components, x and ŷ: For the section  2 , located in the plane of the positive values of , the magnetic potential vector would also have two components, x and ŷ: Nevertheless, in the case of the section  3 , located on the axis  at  = 0, it would only have component ŷ: Thus, once the magnetic potential vector has been calculated, the magnetic field could already be obtained by applying the curl to the magnetic potential vector in the same way as (7).However, with the aim of obtaining the total magnetic field, the calculation must be performed for each loop section.In this manner, the equations of the three sections are presented in ( 16), (17), and ( 18) respectively. ) Hence, all this double loop electromagnetic analysis have led us to the conclusion that the total magnetic field produced by a double loop of dimensions ,  and 2 at a certain point  (, , ) as shown in Figure 4 will be the sum of the components obtained in ( 16), (17), and (18).The resulting expression is shown in the following:

Experimental Measurements
After obtaining the expression that describes the magnetic field generated by a double loop, the next step was to verify that the theoretical results coincided with the experimental ones.For this reason, a double loop was constructed in the laboratory by our research team (Group of Traffic Control System, ITACA Institute, Universitat Politècnica de València, Spain).It was implemented with an external loop of  1 turns and, inside and over it, a smaller one of  2 turns located in the negative half-plane, both with the same direction of circulation.The scheme is shown in Figure 5.
With these conditions, we only were interested in the component of the magnetic field perpendicular to the surface of the loop (  →   ), as the surface vector of the loop only has component along the -axis.Then, taking into account the turns of each of the loop sections, this value was easily obtained from (19) as follows: Signal type is rectangular.Current intensity through the loop () is 51.80 .

Region
Where the Readings of the Magnetic Field Were Carried Out.Height above the plane of the loop is 0.0825 .
Position is along the -axis in the center of the loop ( = 0).
With the above-mentioned characteristics, the theoretical calculation in the region of the measurements was carried out by applying and programming the expression obtained in  (20) in Matlab ©, while the experimental measurements were performed in our laboratory in order to verify the goodness of the theoretical model developed for the calculation of the double loop magnetic field.The instrument used to measure this magnetic field was the Exposure Level Tester ELT-400 shown in Figure 6.
This device contains a series of turns with a diameter of 0.125  and measures the magnetic field by means of this spherical sensor.It is able to detect frequencies from 1 to 400 , although this can be set as wished.To perform the measurements, the ELT-400 was configured as shown in Figure 7 and as follows: Selected frequency range: 30  − 400 .
After collecting and processing all the information, the comparison between the calculated and measured values of the magnetic field as well as the tolerance of the measuring instrument (±35 ) was made, which is shown in Figure 8.
It can be observed that the differences between the measured and calculated magnetic field values are, except in specific points, within the tolerance range of the instrument.The difference between the theoretical and measured values within the contour of the loop is below 20% of the reading and the mean value is below 8%.Therefore, it can be concluded that the theoretical model for double loops developed in this paper predicts with a good precision the behavior of the magnetic field.In addition, different types of tests were also carried out with other types of loops, both single and double, varying the type and amplitude of the applied current, and very similar results were obtained.

Methods for Inductance Calculation
Throughout time, various methods to calculate the inductance of magnetic loops according to different geometric configurations have been proposed [22][23][24][25][26][27][28], such as the Mills' method based on Grover's equations [22].Most of the mathematical expressions for these calculations appear in classical texts [23] and based on these expressions, the inductance values presented by other different models of loops can be easily deduced.
However, the development of computer systems has allowed to implement numerical methods that make use of the intrinsic definition of the physical process of magnetic induction.In this regard, it would be convenient to include our previously presented work [25], in which the methods to calculate the inductance of a rectangular loop were analyzed and compared with other calculation techniques and real measures.
In this way, this time we will present three methods to calculate the inductance of the innovative double loops, since there are no studies about them.Therefore, these methods will be deeply analyzed and compared in order to know which method is more effective and if all of them are equally good.These three methods are as follows: (A) Electromagnetic Analysis Method (B) Numerical Integration Method (C) Mills and Grover's Method (A) Electromagnetic Analysis Method.The first method to calculate the inductance of a double loop is based on the electromagnetic analysis, which considers the use of a numerical method that employs the expression of the magnetic field generated by a double loop.This means to carry out a study applied to the case of a double loop with the same structure as Figure 4.
The procedure starts by making use of the flux calculation, since all the loops work in the same way.However, in this case, the magnetic field and the flux will present some peculiarities.
The magnetic flux will be obtained as the integral of the product of the magnetic field by the differential surface all along the entire surface of the loop: The expression of  →   used in (22), which only taken into account the field in the  direction, will be the same as in (20), and the differential surface can be easily replaced by the product of the differential length according to the -axis and -axis as shown in the following: To solve this expression, the integrals will be replaced by summations.For this, it is necessary to obtain the magnetic field in a series of points in space and consider a differential surface  around.
In this way, if we defined   as the number of points in which the magnetic field will be measured along the -axis and   as the number of points in which the magnetic field will be measured along the -axis, the  and  differentials of ( 23) would be given by the following: Consequently, the total flux through the loop would be as shown in the following: Summaries' limits have been designed to prevent the magnetic field measurements on the conductors, since at these points the value of   presents a singularity.The summation is extended to the rectangle defined by the points (−+, − ) in the -axis and (− + ,  − ) in the -axis, but this causes an error in the measurement that has been tried to be solved by increasing the differential surface by 50% within the limits of the summation points.For this reason,   and   factors have been introduced [23].The   factor takes a value 1 in all points except in those where  = 1 or  =   − 1 that has a value 1.5, and the   factor also takes a value 1 in all points except in those where  = 1 or  =   − 1 that has a value 1.5.
On the other hand, due to the abrupt change that appears in the   component of the magnetic field in the vicinity of the conductors, it is evident the importance of choosing the number of points   and   .In our previous article [25], it was shown that the optimal values of   and   to minimize the error in the calculation of the inductance for a single rectangular loop were given by where   is the radius of the conductor used to build the loop and ,  and 2, the dimensions of the loop shown in Figures 4 and 5.However, in this case, the space has been divided into three sections (two areas) as explained previously.For the positive region of the -axis the flux is concatenated by  1 turns, and for the negative region of the -axis, the flux is concatenated by  1 ±  2 turns as shown in Figure 5.The "+" sign would be used when the current goes through both loops in the same direction and the sign "−" when the current goes through the loops in opposite directions.Therefore, we could separate   into two sections as shown in Equation ( 27), which would represent the measured points, respectively, the negative and the positive ones at the -axis.
Finally, after all this mathematical and electromagnetic analysis, we could obtain the inductance of a double loop with dimensions ,  and 2 and with  1 and  1 ±  2 turns respectively by the following expression: In this expression,  represents the intensity used for the calculation of the magnetic field and the rest of values are known and have been described above.Nevertheless, it should be pointed out that the result of this first method will give the value of the inductance without considering the thickness of the loop, since this method is based only on the electromagnetic analysis of the magnetic field created by a double loop.Therefore, in this first method the spacing between the turns of the loop is considered null.
(B) Numerical Integration Method.After obtaining a first way to calculate the inductance of a double loop that does not take into account the spacing between turns, it seems reasonable that the second presented method does.These method, which apparently should provide more accurate values, is the utilization of numerical integration techniques but considering the size and spacing between the turns of the loop ( V ).However, on this occasion the calculation will be more complex than the last one, but these values are supposed to be much more real, since in reality the cables have a thickness that although it is minimal, it should be contemplated.
For the purpose of calculation mentioned, it is assumed an assembly in which the loops are stacked vertically.First, the larger  1 turns would be installed and after, over them, the smaller  2 , assuming that are equi-spaced a distance  V .In this way, for each turn of the loop, the flux through it would be produced by the current that flows through the turn itself plus the flux generated by each of the other turns of the loop and which is concatenated by the one that is being analyzed.
The different flux components mentioned above will be represented by the terms 0  , where the subindex  indicates the type of loop which generates the field (1 Big, 2 Small), the subindex  indicates the type of loop which is crossed by the magnetic field (1 Big, 2 Small) and the subindex  indicates the distance between the loop that generates de field and the one which detects it.Hence, we can define the following terms: (i) 0 110 is the flux generated by one of the turns of the largest loops and which goes through them.(ii) 0 220 is the flux generated by one of the turns of the smallest loops and which goes through them.(iii) 0 11 is the flux generated by a big turn which goes through the other big ones separated a distance  V .(iv) 0 12 is the flux generated by a big loop which goes through the other small ones separated a distance of  V from the turn that has generated the flux.
(v) 0 22 is the flux generated by a small turn which goes through the other small ones separated a distance of  V .
(vi) 0 21 is the flux generated by a small loop which goes through the other big ones separated a distance  V from the turn that has generated the flux.(vii)   will be the minimum value between  1 and  2 and   will be the maximum value between  1 and  2 .
With this nomenclature, the inductance of the double loop would be given by the following: The "+" sign would be used when the current goes through both loops in the same direction while the sign "−" would be used when the current goes through the loops in opposite directions.In any case, the fluxes described and shown in (29) are given by the following: In 0 11 ,  1 (, , ) is the magnetic field component along the -axis generated by a large turn, extending from − to + along the -axis and from − to + along the -axis at the point of space (, , ).In 0 12 ,  1 (, , ) is also the magnetic field component along the -axis generated by a large turn, extending from − to + along the -axis and from − to + along the -axis at the point of space (, , ).
In 0 22 ,  2 (, , ) is the magnetic field component along the -axis generated by a small turn, extending from − to 0 along the -axis and from − to + along the -axis at the point of space (, , ).
In 0 21 ,  2 (, , ) is also the magnetic field component along the -axis generated by a small turn, extending from − to 0 along the -axis and from − to + along the -axis at the point of space (, , ).
The values used in (29) and (30) ( 1 ,  1 , , ,   ,   ) are the same ones used to calculate the previous inductance and   1 takes the following value: As it can be seen, this method seems much more accurate than the previous one, but it is true that there is a high amount of operations and summations to perform, which will make it more difficult to implement and with higher computational cost.Therefore, taking into account the fact that the first method did not consider spacing between the turns of the loop, which is not physically correct, and that the second method, although it does, is computationally complicated because of the large number of operations and summations, it leads us to think of a third method that also considers spacing but is much simpler to implement.
(C) Mills and Grover's Method.In order to calculate the inductance of a double loop, Mills and Grover's method will be proposed.As said above, in this third and last method it will be taken into account again that the size of the conductors prevents all turns from being in the same coordinate .To consider this phenomenon, studies such as Mills' emerged [22,28], from which new expressions of inductance were deduced.
When working with double loops, the self-inductance of the  1 turns of the loop and, over them, the  2 turns, equally spaced a distance  V , is represented by the following: where (i)  10 is the self-inductance of a rectangular loop with only one large-size turn ( + )  2; (ii)  20 is the self-inductance of a rectangular loop with also only one small-size turn (  2); (iii)  11 ( V ) is the mutual inductance between two turns of a rectangular large loop separated between them at a distance of  V ; (iv)  22 ( V ) is the mutual inductance between two turns of a rectangular small loop separated between them at a distance of  V ; (v)  12 ( V ( 1 +  − )) are the mutual inductances between a big-size turn and a small-size turn separated at a distance of  V ( 1 +  − ).
It must be noted that the above-mentioned parameters  10 and  20 are obtained as the sum of the internal and external inductance of the conductors that constitute the loop as follows: where   is the inductance per unit length, which must consider the relationship between the inductance at a certain frequency   and the inductance at low frequency  0 .In order to carry out these operations, we will use Johnson's studies [23].In this manner, where (), (),   (), and   () are the imaginary and real parts of the Bessel function of first order and their derivatives, which are necessary for the calculations when this method is applied.On the other hand, the inductance per unit length at a low frequency is given by ( 35), but for copper conductors it takes a value of 0.5 ⋅ 10 −7 /: To obtain   [22][23][24][25][26][27], we focus on the mutual inductance of a pair of parallel conductors, whose expression is: As seen in Figure 9,  is the length of the filaments and  is the separation between them (both quantities expressed in meters), which results in an inductance expressed in Henrys.This expression takes a positive sign when the current in both cables has the same direction and takes a negative sign when the current has opposite directions.In addition, it is know that the external inductance of a pair of parallel conductors with the dimensions shown in Figure 10 is given by 1 and  2 are the self-inductances of the simple conductors and  12 and  21 are the mutual inductance measures between the centers of the conductors, assuming a uniform current distribution throughout the cross section of the conductor.Mutual inductances have been considered to have a negative sign as it has been assumed that the current in the two conductors have opposite directions.Therefore, if both conductors have the same dimensions: Thus: The external self-inductance of a conductor is obtained thanks to a method in which the conductor is replaced with two conductors with a null straight section separated by a distance equal to the radius of the conductor.In this waym In addition, the inductance of a rectangular loop with a single turn is given by the sum of the inductance of two pairs of parallel conductors as follows: where  1 =  +  and  2 = 2 for the big loop and  1 =  and  2 = 2 for the small loop.
As it can be deduced from this expression, the external inductance of a rectangular loop with one turn is equal to the mutual inductance of two identical coaxial rectangular loops separated by a distance equal to the radius of the conductor.In this way, the mutual inductance of two parallel rectangular loops as shown in Figure 11 can be obtained from mutual inductances between parallel conductors.
Therefore, the mutual inductance between the two rectangular loops with the same dimensions as shown in Figure 11 may be expressed as follows: In the above expression, the terms   represent the mutual inductance between the segment  of the bottom loop and the segment  of the top loop.The reason to multiply the expression by two is because the mutual inductances are all symmetric, that is, for all  and all ,   =   .
On the other hand, to calculate the mutual inductance between two parallel loops with different dimensions, Grover's equations must be used again, since they provide the mutual inductance between two parallel straight conductors as shown in Figure 12.According to Grover's formula, the mutual inductance between two parallel conductors with sizes  and  spaced a distance  and displaced at a distance  is given by where It should be pointed out that if the two conductors overlap partially or totally, the parameter  will have negative values.
In any case, from all these expressions it is finally possible to obtain the inductance between two parallel rectangular loops with different dimensions as shown in Figure 13, which is simply another way of seeing Figures 4 and 5.
Therefore, the mutual inductance between the two parallel loops is finally obtained as a sum of the mutual inductances of the parallel conductors where the terms   represent mutual inductances between two parallel rectilinear conductors with the same dimensions, and the terms   represent the mutual inductances between two parallel rectilinear conductors with different dimensions.This is shown in the following: Although this last equation may seem complicated and with a large number of operations, it should not be forgotten that in none of the previous expressions there are integrals or summations.The maximum operations that appear are square roots and trigonometric functions.Therefore, it apparently seems to be the method that best behaves computationally.Moreover, this also considered spacing between the turns of the loop as the second one.
Thus, once all methods have been presented, the inductance values obtained are going to be compared to each other to verify the similarity between them and know which the best methods are and which are not.

Results and Discussion
After making the description and performing the analysis of the three methods, the first thing that stands out is that there are two methods that consider the spacing between the turns ( V ) and one which does not.For that reason, what would be expected at first sight is that the accuracy and complexity of the two methods which do consider this are greater than the method which does not.Therefore, the inductance values obtained by using the electromagnetic analysis method should provide different values than the other ones, but a simple explanation for this phenomenon would be the fact of not having considered the flux lines that are lost between the loops because of the fact that they have a real thickness which is not zero.
To check all of the above, various test related to the inductance values of the different methods and types of loops were performed and are presented below.In fact, it was studied the effect of increasing the number of turns in a single and in a double loop and the effect of changing the dimensions of a double loop.

How the Number of Turns Affects the Value of the
Inductance of a Single Loop.For the purpose of this study, a single 22 loop located and centered in the  plane with a cable radius of 0.75 and a turn spacing of 1.9 was used.In this first test, the number of turns of the loop was gradually increased to see the effect that it has on the value of its inductance.With the nomenclature described in the paper, the characteristics of this loop were as follows: As mentioned, in the ideal case it would be expected for the three methods to provide similar results, but it is clearly observed in Figure 14(a) that there are two methods that are practically identical, and one method whose results are far from the other ones.Therefore, the first point to emphasize is that the Mills and Grover's method and the numerical integration method, those that consider the separation between turns, provide good and identical results, while the method based on the electromagnetic analysis begins to fail when the turns increase.
This means that the more turns the loops have, the more different the results are, because the error of not considering separation between turns increases for each turn.

How the Number of Turns Affects the Value of the
Inductance of a Double Loop.The next test was the same as the previous one but instead of working with a single loop, we will work with a double loop.This time, it was a double 22 loop formed by an external loop of  1 turns located and centered in the  plane and a smaller one of  2 turns located in the half-plane of the negative values of , both with a cable radius of 0.75 and a turn spacing of 1.9.In the previous analysis, only the number of turns  1 was varied since it was a single loop and  2 = 0.However, as we were working with a double loop, to check the effect of increasing the number of turns, the value of  1 was kept fixed and  2 was the value that we were gradually increasing.In this way, with the nomenclature described in the paper, the characteristics of this double loop were as follows: The result obtained was really similar to the previous one.As it can be seen in Figure 14(b), the two methods that consider spacing provided practically identical values, while the electromagnetic analysis method, as it does not, provided different values, this time greater.The fact that the results were even greater is because as we had concluded before, the more turns the loops have, the more different the results are, and there are more turns in a double loop ( 1 in the first loop and  1 +  2 in the second one) than in a single one (only  1 ).Thus, at this point we could already begin to draw conclusions, but in order to make a study as rigorous as possible, we will show more results from our analysis.

How the Dimensions of the Loop Affect the Value of the
Inductance of a Double Loop.Once we had analyzed how the number of turns affects the inductance value, in the two remaining experiments we analyzed what happens when varying the dimensions of the loops.For this purpose, we worked with a double loop formed by  1 large turns, centered with respect to the -axis, extending from − to + along the -axis, and  2 small turns, also centered with respect to the -axis, extending from − to 0 along the -axis.Both had a cable radius of 0.75 and a turn spacing of 1.9.In this way, the number of turns was kept fixed and the value that we varied was the length , which means increasing the length of the smallest loop, which has  2 turns and is located in the negative half-plane.In addition, it should be noted that this fact also increases the length of the biggest loop, since its dimensions are ( + )  2.Then, with the nomenclature described in the paper, the characteristics of this loop were as follows: In Figure 14(c) it can be observed very clearly that although the trend of the three methods was the same, the results were different.Mills and Grover's method as well as the numerical integration one had a practically linear relationship between the increase of  and the increase of the inductance.However, it is noted that the first method, which had extracted its equations from the analysis of the electromagnetic field without considering spacing, gave results distant from the real ones from the beginning, being this error greater as the length of  is bigger.
On the other hand, after analyzing what happens when the length of  is increased, it was analyzed what happens precisely if we increase the length of , with whose variation we only be modify the length of the largest loop.The double loop was again formed by  1 large turns, centered with respect to the -axis, extending from − to + along the axis, and  2 small turns, also centered with respect to the -axis and extending from − to 0 along the -axis.Both had a cable radius of 0.75 and a turn spacing of 1.9.The number of turns was kept fixed again and the value that we varied was .In this manner, with the nomenclature described in the paper, the characteristics of this loop were as follows: (i)  = 1  (ii)  = 1  (iii)  =  1  10  (iv)  1 = 3 In this test, it was seen how two methods followed the same trend again, providing almost the same results, but the other one did not behave properly.However, in Figure 14(d) it is clear that Mills and Grover's method and the method of numerical integration provide exactly the same values and the electromagnetic analysis method does not.
Therefore, after analyzing the previous results, we could conclude that of the three methods proposed, the electromagnetic method could be useful for very thin conductors with little separation between them and for complex geometries in which Mills and Grover can not be used, as it only serves for parallel and perpendicular conductors.When it was tried to increase the turns of the loop, it was the only one that differed from the rest.On the other hand, when the length of the loop was increased, it was also evident that this method did not behave correctly.For this reason, we can affirm that although it is a valid method and can give us an approximation if the separation between turns is big, it should be only used as a reliable source if that separation is minimal.
Regarding the two remaining methods, it must be noted that both offer good and similar results as it can be seen in Figures 14(a), 14(b), 14(c) and 14(d), but in the case of having to opt for one of them, we must emphasize that the numerical integration method carries a much higher computational cost than the Mills and Grover's one.In fact, it is clear by observing Figure 14 that this method had the best behavior in every of the cases.Then, we could conclude that because of its low computational cost and extreme accuracy, Mills and Grover's method would be the best way of calculating the inductance of a double magnetic loop, which is why it will be the one that we will use in our simulation programs.

Conclusions
This article is aimed to be a presentation of the double loop, where geometry, construction, operating mode and three possible ways to calculate its inductance have been explained.After presenting these three above-mentioned methods, an analysis has clarified that if precision is required, Mills and Grover's method or the numerical integration method must be used, as they both take into account the separation between turns, although we recommend to choose the first one because of its low computational cost.
In future papers, we will focus more closely on the advantages offered by using this type of loops, what will help to understand the need and importance of this paper.It will focus on the new vehicle magnetic profiles, the parameters that can be extracted from them and the benefits of using them in comparison with the conventional loops.
The reality is that magnetic loops, despite being from the eighties, are still the most used technology to capture data from traffic.For that reason, we must improve the existing infrastructure and provide this sensor with greater potential and reliability.

Figure 4 :
Figure 4: Double loop presented in three sections.

Figure 5 :
Figure 5: Outline of the double loop constructed.

) 3 . 1 . 2 .
Characteristics of the Double Loop Constructed.Conductor type is tinned copper wire conductors individually insulated with polyvinyl chloride with a cross section of 0.28  2 .Number of turns is 4 exterior turns ( 1 ) of 1.20  0.46  and 5 interior turns ( 2 ) of 0.40  0.46 .Dimensions is 1.20  0.46 .With these values, the parameters of the loop according to the nomenclature used in Figures 4 and 5 would be as follows: Characteristics of the Signal Used to Energize the Loop.Frequency of the signal applied to the loop is 139.2 .

Figure 8 :
Figure 8: Calculated and measured values of the magnetic field |B|.The instrument tolerance is also considered.

Journal of Advanced Transportation l d Figure 9 :Figure 10 :
Figure 9: Features of two ideal parallel conductors of no straight section for measuring mutual induction.

Figure 11 :Figure 12 :
Figure 11: Geometry for calculating the mutual inductance between two parallel and coaxial rectangular loops with the same dimensions.

Figure 13 :
Figure 13: Disposition of two parallel loops of different dimensions.

Figure 14 :
Figure 14: Inductance values.(a) When increasing the number of turns in a single loop.(b) When increasing the number of turns in a double loop.(c) When increasing the length of "a" in a double loop.(d) When increasing the length of "d" in a double loop.