*5.2. A Pratical Device Example*

The second test problem is the modelling of a practical device proposed by the EDF (Electricité de France) [26]. The power station "Folies" is equipped with a three-phase reactance to limit short-circuit currents. In this case, the current in each reactance is 1000 A phase-shifted with 120 degrees, the frequency is 50 Hz (Figure 5).

**Figure 5.** Three phase reactance in the power station "Folies".

The parameters are presented as below (Figure 6):


**Figure 6.** Parameters of three reactances.

In nominal conditions, a three-phase reactance generate a leakage induction in the neighborhoods. In order to minimize this electromagnetic disturbance, an electromagnetic shielding and a passive loop are added between the magnetic field source and the protected area (Figure 7). The device parameters are presented as below:

	- -Permeability μ<sup>r</sup> = 1,
	- -Electrical conductivity <sup>σ</sup> <sup>=</sup> 3.03 <sup>×</sup> 107 <sup>S</sup>/<sup>m</sup>
	- -Section radius rs = 9.25 mm
	- -Zloop = 3.85 m
	- -Permeability μr = 20,000,
	- -Electrical conductivity σ = 2.2 × 106 S/m
	- -Thickness = 3.5 mm
	- -Zshielding = 5.15 m
	- -The skin depth δ = 0.339mm is thinner than the shielding's thickness.

The last case is tested by our integral method and the FEM 3D with shell elements. Let us note that this example is in 3-dimensional space and cannot be modelled by the 2D FEM. The current values in the passive loop and the current distribution in the shells are also compared.

The obtained results from the coupling method converge quite close to the current values presented in Table 2. Figure 8 also shows that the surface distribution of the current in the shell is quite similar to the two methods. The results achieved by the coupling integral method are very encouraging. The convergence is reached with few elements (about 1000 elements).

**Figure 7.** Geometry of the test case "Folies".

**Table 2.** Current values in the passive loop, where j is an imaginary unit.


**Figure 8.** Surface distribution of the current (A/m) in the thin shell for the test case "Folies": (**a**) FEM 3D with shell elements; (**b**) Integral method.

However, some matrices of Equation (38) are fully populated and compression algorithms must be applied if there is a large number of elements. The coupling of the model order reduction techniques or the matrix compression algorithms with some integral methods clearly demonstrates its efficiency. For example, the coupling of a matrix compression algorithm like the FMM and the MoM or the PEEC method has reduced the computation time and the memory requirements down to more than 10 times and the compressed ratio is more than 80 percent [18,27]. Moreover, the acquired model can be reused to build a real circuit, which is easy to employ in all conventional SPICEs-like circuit solvers [28,29].
