Analytical Solution: Infinite Magnetic Line-source

A comparison of the 2D “analytic” solution of an Infinite Magnetic Line-source sine wave (with a carrier frequency of 5 GHz) to the output of the 2D FDTD Yee simulation.


Analytical Solution

The analytical solution to the Infinite Magnetic Line Source problem is from Balanis, Advanced Engineering Electromagnetics, pp.570-574 and is:

860343133088ed552a404fd68d5ea0c7.png Eqn:6.1


df5e12fd03cdee0c14946806ab633179.pngmagnetic line source constant current
ee53ed02e07929d204850b677c1d6929.pngangular frequency
a1f21f33e4cc224c1aa90a478a5339d5.pngelectric permittivity (farads/meter)
2b3a144c31b05e0a773b0b60dc13e991.pngmagnetic permeability (Henrys/meter)
4e5a86924aeba8376b7eaaafbdde58ab.pngdistance from source
666ca07012f4b30d0b419f35349c5042.pngHankel function of the second kind of order 0.

Scilab Code

The algorithm was run using scilab (an open source numerical computation software package, similar to Matlab).

Scilab code:

//  Code for infinite magnetic line source
//  usage: scilab
//    getf hello.sci
//    hello
function hello

  mu = 1.256637061e-06;
  epsilon = 8.854187818e-12;
  omega = 3.141592654e+10;
  dx = 0.003;
  current = 7.0;

  bta = omega * sqrt( mu * epsilon );
  temp = bta * bta / (4.0 * omega * mu);
  phaseadj = 2.0 * dx * bta;

  fd = mopen( 'foo1.txt' , 'w');
  for istep = 1.0 : 1.0 : 500
    x = istep * dx;
    j = besselj(0,x);
    y = bessely(0,x);
    h = besselh(0,2, bta * x + phaseadj);  // hankel function 2nd kind 0th order
    hz = -current * temp * h;
//    mfprintf(fd, 'i: %f x: %f  j0: %f   y0: %f   h0: %f hz:%fn', istep,x,j,y,h,hz);
    mfprintf(fd, '%d %fn', istep,hz);



A 2d FDTD (TE) simulation was run with a grid area of 1000 x 1000.

The source excitation was placed int the center (x=500,y=500) and was an Infinite Magnetic Line-source sine wave with a carrier frequency of 5 GHz.

The simulation was run for 5000 time steps (to allow the simulation to reach steady state) and then a “snapshot” of the values from x = 500 to 1000 was recorded.

The results of the FDTD simulation are compared with the analytical solution:


As can be seen, the FDTD simulation matches the analytical solution.

See Also


External Links