Magnetic dipole radiation through a multilayered structure

Calculation of the radiation of a magnetic dipole

Description

The function implements an horizontal magnetic dipole (along the x-axis) and calculates the x and z component of the magnetic field for y=0. The script allows to set a wanted value of multiple layers, with their own relative electric permittivity, relative magnetic permeability, conductivity (see the LAYERS PARAMETERS section in the script). The script finally produces an animation of the radiated field.

• MATLAB release: MATLAB 7.4 (R2007a)

[wei]

How could it work?
Can Lorenzo explain to us ?

[Lorenzo Luini]

@wei: the radiation is calculated as follows:
- first the exact solution of the magnetic dipole radiation is decomposed into several more simple wave functions (plane waves or cylindrical waves)
- each simple wave interacts with the multilayered structure according to the customary EM boundary conditions
- the final radiation pattern is obtained by summing up all the contributions of the simple waves.

[Alexander Patrushev]

Very interesting!

Lorenzo, could you tell is that possible to compute full EM vector in the similar way?

[Alexander Patrushev]

Also, what do primary and secondary fields mean?

[Lorenzo Luini]

@Alexander: The primary field is the radiation emitted by the EM source, in this case, the magnetic dipole; the secondary field originates from the interation of the primary field with the multilayered structure, i.e. consists in the reflected and transmitted components of the field. Coming to your first question (if I have understood it correctly), the method can be applied to most of the radiation patterns that are known in a closed-form solution (both concerning E and H, therefore it is possible to calculate the full EM verctor) . The tough part may be to mathematically decompose the primary field into several elementary waves. Afterwards, the propagation of simpler waves is relatively easy.

[Alexander Patrushev]

@Lorenzo: thanks for explanation. Is there any doc or article to understand how you get used math?

[Alexander Patrushev]

@Lorenzo: could you specify which boundary conditions used and where?

[Lorenzo Luini]

@Alexander: unfortunately I have no documentation that could be useful to explain the whole procedure.
Concerning the boundary conditions: they are the usual ones for EM fields, i.e the conservation of the total tangential H component and of the total normal B component at each interface between 2 layers.

[Liping]

The result seem very interesting, can you please gave any literature i.e. books or papers on which you based to do the modelling?

[Lorenzo Luini]

@Liping: unfortunately the only source I have is in Italian and it is published by my university (Politecnico di Milano). However, the rationale is described above as an aswer to the question of "wei".