Roy Lewallen wrote:

Define:

eta=120 Pi
k=2/lambda

reactance = (eta/(4*Pi)) (2 SinIntegral[k l] +
Cos[k l]*(2 SinIntegral[k l] - SinIntegral[2 k l]) -
Sin[k l]*(2 CosIntegral[k l] - CosIntegral[2 k l] -
CosIntegral[(2 k a^2)/l]));

where 'a' is the radius.


. . .

This is the formulation by S.A. Schelkunoff, which is one of many
approximations to the general problem of finding the reactance of a
simple cylindrical dipole of arbitrary length and diameter. The general
problem was attacked for decades by some very skilled mathematicians and
engineers including R.W.P. King, David Middleton, Charles Harrison, G.H.
Brown, D. D. King, F. G. Blake, M.C. Gray, and others. You'll find their
works scattered about the IRE (now IEEE), British IEE, and various
physics journals. The problem can't be solved in closed form, so all
these people proposed various approximations, some of which work better
in some situations and others in others. A good overview can be found in
"The Thin Cylindrical Antenna: A Comparison of Theories, by David
Middleton and Rolond King, in _J. of Applied Physics_, Vol. 17, April 1946.


Thank you for that. If by chance you have that as a PDF, perhaps you can
mail it to me. But if not, I'll try to get it for interest sake. I
needed this for a piece of work, but the work will have finished by the
time I get much more done. But at least I have a better understanding now.


. . .

Does anyone have any comments? Any idea if Balanis's work is more
accurate? It is more up to date, but perhaps its an approximation and
the amateur radio book is more accurate. (The ham book seems quite
well researched, and is not full of the voodoo that appears in a lot
of ham books).

As I mentioned above, some approximations are better in some
circumstances (e.g., dipoles of moderate diameter near a half wave in
length) and some in others (e.g. fat dipoles or ones near multiples of a
half wave in length). I don't know which is better for your particular
question. The easy way to find out is to get one of the readily
available antenna modeling programs, any of which is capable of
calculating the answer to very high accuracy, and compare this correct
answer with the various approximations you find published.

OK. I'm just a bit suspicious of computer programs some times, as
someone will have to choose an algorithm of some sort. But I assume you
are talking of something like NEC which breaks antennas into segments.

BTW, I'm also looking for an exact formula for input resistance of a
dipole of arbitrary length. I know is 73.13 Ohms when 0.5 wavelengths
long, but I'm not sure exactly how much it varies when the length
changes (I believe it is not a lot).

There is no exact formula for that, either. Calculating an exact answer
requires knowledge of the current distribution, which is a function of
wire diameter. Assuming a sinusoidal distribution gets you very close
for thin dipoles, but it's not exact. You'll find calculations based on
this assumption in just about any antenna text such as Balanis or Kraus.

Balanis has it, but leaves it as an integral, without simplifying like
he does for the real part. Yet the formuals for hte real and imaginary
parts look very similar. I might be able to attack it with a computer
algebra system - maths never was my strongest subject.

I thought I'd looked in Krauss and not found it, but perhaps it is
there. I think there is a relatively new version of Kraus, but my copy
is quite old.

But again, you can get extremely accurate results from readily available
antenna modeling programs.

OK, thank you for that.

Roy Lewallen, W7EL

Dave wrote:
. . .

OK. I'm just a bit suspicious of computer programs some times, as
someone will have to choose an algorithm of some sort.

.. . .

I might be able to attack it with a computer
algebra system - maths never was my strongest subject.

I suppose it's natural to be more suspicious of others' work than your
own. I've personally found the opposite to often be more appropriate.

. . .

But I assume you
are talking of something like NEC which breaks antennas into segments.

Yes. You can find a good description of the method of moments in the
second and later editions of Kraus. The fundamental equation can only be
solved numerically, and the method of moments, used by NEC and MININEC,
is an efficient way to do it.

BTW, I'm also looking for an exact formula for input resistance of a
dipole of arbitrary length. I know is 73.13 Ohms when 0.5 wavelengths
long, but I'm not sure exactly how much it varies when the length
changes (I believe it is not a lot).

There is no exact formula for that, either. Calculating an exact
answer requires knowledge of the current distribution, which is a
function of wire diameter. Assuming a sinusoidal distribution gets you
very close for thin dipoles, but it's not exact. You'll find
calculations based on this assumption in just about any antenna text
such as Balanis or Kraus.

Balanis has it, but leaves it as an integral, without simplifying like
he does for the real part. Yet the formuals for hte real and imaginary
parts look very similar. I might be able to attack it with a computer
algebra system - maths never was my strongest subject.

Hallen's integral equation is exact, but it's not a formula, since you
can't plug numbers into one side and get a result on the other. Nor can
it be solved in closed form at all. That's why so much work was done on
approximate solutions and on developing numerical solution methods. Feel
free to write your own program to solve it, but such programs have
existed for decades and have been verified countless times as well as
being highly optimized.

I thought I'd looked in Krauss and not found it, but perhaps it is
there. I think there is a relatively new version of Kraus, but my copy
is quite old.

Getting the resistance is pretty straightforward once you assume the
shape of the current distribution. Assume some arbitrary current at the
feedpoint which, along with the assumed current distribution, gives you
the field strength in any direction. With the impedance of free space,
this directly gives the power density. Integrate the power density over
all space to get the total radiated power. Then you know how much power
is radiated per ampere of current at the feedpoint, from which you can
calculate the feedpoint resistance.

This calculation is done in all editions of Kraus, I'm sure; I have only
the first and second, but I can't imagine it was deleted in later ones.
Be careful when reading Kraus, however. Unlike many authors, he uses a
uniform, rather than triangular, current distribution for his short
elemental dipole examples. This is equivalent to a very short dipole
with huge end hats, not just a plain short dipole. The half wavelength
and other dipoles in his text are conventional.

. . .

Roy Lewallen, W7EL

Roy Lewallen wrote:

But I assume you
are talking of something like NEC which breaks antennas into segments.

Yes. You can find a good description of the method of moments in the
second and later editions of Kraus. The fundamental equation can only be
solved numerically, and the method of moments, used by NEC and MININEC,
is an efficient way to do it.

Since you clearly know more about this stuff than me, do you know of the
best freely available software for this which works under Unix? (I use
Sun's Solaris for 99% of the things I do, including sending this
message. I use Solaris on my laptop too, rather than Windows).

Hence I'm almost certainly looking for source code in either C, C++ or
Fortran. Anything that works under Linux would almost certainly be able
to be compiled for Solaris without too much effort.
I found this page:

http://www.si-list.net/swindex.html

which has some source. I downloaded one

http://www.si-list.net/NEC_Archives/necpp-1.1.1.tar.gz

It would not compile immediately on my Sun. gcc 4.3.1 complained about
some ambiguous code. gcc 3.4.1 did not, so I got past that bit.

It then tries to link with the 'blas', 'atlas' and 'lapack_atlas'
libraries, none of which my Sun has.

I then swapped to the Sun C/C++ and Fortran compilers, removed
references to 'blas', 'atlas' and 'lapack_atlas' , and replaced them
with 'sublibperf' which is the optimised library on Solaris. That worked
ok, and I had an executable:

$ ./nec2++
usage: nec2++ [-iinput-file-name] [-ooutput-file-name]
-g: print maximum gain to stdout.
-b: Perform NEC++ Benchmark.
-h: print this usage information and exit.
-v: print nec2++ version number and exit.


I've not done any more than that at this point, but proved it will
compile on Solaris with little effort.

Anyway, if you have any recommendations for the best freely available
Unix/Linux code, I would be interested.


Hallen's integral equation is exact, but it's not a formula, since you
can't plug numbers into one side and get a result on the other. Nor can
it be solved in closed form at all. That's why so much work was done on
approximate solutions and on developing numerical solution methods. Feel
free to write your own program to solve it, but such programs have
existed for decades and have been verified countless times as well as
being highly optimized.

OK, I understand that.

Getting the resistance is pretty straightforward once you assume the
shape of the current distribution. Assume some arbitrary current at the
feedpoint which, along with the assumed current distribution, gives you
the field strength in any direction. With the impedance of free space,
this directly gives the power density. Integrate the power density over
all space to get the total radiated power. Then you know how much power
is radiated per ampere of current at the feedpoint, from which you can
calculate the feedpoint resistance.


This calculation is done in all editions of Kraus, I'm sure; I have only
the first and second, but I can't imagine it was deleted in later ones.
Be careful when reading Kraus, however. Unlike many authors, he uses a
uniform, rather than triangular, current distribution for his short
elemental dipole examples. This is equivalent to a very short dipole
with huge end hats, not just a plain short dipole. The half wavelength
and other dipoles in his text are conventional.

I think I found what I was looking for in either Kraus or Balanis last
night. The book is beside the bed, and as my wife is still asleep I'm
not going to look for it.

Dave wrote:

Since you clearly know more about this stuff than me, do you know of the
best freely available software for this which works under Unix? (I use
Sun's Solaris for 99% of the things I do, including sending this
message. I use Solaris on my laptop too, rather than Windows).

Hence I'm almost certainly looking for source code in either C, C++ or
Fortran. Anything that works under Linux would almost certainly be able
to be compiled for Solaris without too much effort.
I found this page:

http://www.si-list.net/swindex.html

which has some source. I downloaded one

http://www.si-list.net/NEC_Archives/necpp-1.1.1.tar.gz
. . .

Sorry, my knowledge doesn't extend to that of programs suitable for Unix
or Linux. The "Unofficial NEC archives" site is the best I know of for
various compilations. Hopefully some of the other readers can help you out.

My program, EZNEC, has been reported to run under Linux using some
versions of the wine Windows emulator, but not with others. You'll have
a lot more to choose from if you can emulate Windows.

Roy Lewallen, W7EL

Roy Lewallen wrote:
Dave wrote:

Since you clearly know more about this stuff than me, do you know of
the best freely available software for this which works under Unix? (I
use Sun's Solaris for 99% of the things I do, including sending this
message. I use Solaris on my laptop too, rather than Windows).

Hence I'm almost certainly looking for source code in either C, C++ or
Fortran. Anything that works under Linux would almost certainly be
able to be compiled for Solaris without too much effort.
I found this page:

http://www.si-list.net/swindex.html

which has some source. I downloaded one

http://www.si-list.net/NEC_Archives/necpp-1.1.1.tar.gz
. . .

Sorry, my knowledge doesn't extend to that of programs suitable for Unix
or Linux. The "Unofficial NEC archives" site is the best I know of for
various compilations. Hopefully some of the other readers can help you out.

Thank you.

My program, EZNEC, has been reported to run under Linux using some
versions of the wine Windows emulator, but not with others. You'll have
a lot more to choose from if you can emulate Windows.

Roy Lewallen, W7EL


Thanks a lot. I'm just a bit anti-windows myself. Fed up with all the
hassle of viruses etc.

The CPUs in this machine are not even capable or running Windows, as
they are not AMD/Intel compatible.

I do of course have access to Windows machines, and might well check out
your code later. I assume analysing a simple dipole will be trivual and
not require a huge amount of work.

But you have helped me understand the differences in impedance between
the different sets of data I have seen, which was my main concern. For
the immediate future at least, I am not going to do any numerical
analysis, but I might well look at that for a later date.


Dave

Dave wrote:
Roy Lewallen wrote:


Since you clearly know more about this stuff than me, do you know of the
best freely available software for this which works under Unix? (I use
Sun's Solaris for 99% of the things I do, including sending this
message. I use Solaris on my laptop too, rather than Windows).

Hence I'm almost certainly looking for source code in either C, C++ or
Fortran. Anything that works under Linux would almost certainly be able
to be compiled for Solaris without too much effort.

FORTRAN would be the language of choice (since that's what NEC was
written *and validated* in.. one would be concerned about a C
translation, although I'm sure there are C versions out there which have
been validated)

It then tries to link with the 'blas', 'atlas' and 'lapack_atlas'
libraries, none of which my Sun has.

There should be versions out there that don't link with the matrix math
packages.


Anyway, if you have any recommendations for the best freely available
Unix/Linux code, I would be interested.

What you've got is probably as good as anything else, especially if
you're just looking for a table of Z vs length and diameter.

Jim Lux wrote:
Dave wrote:
Roy Lewallen wrote:


Since you clearly know more about this stuff than me, do you know of
the best freely available software for this which works under Unix? (I
use Sun's Solaris for 99% of the things I do, including sending this
message. I use Solaris on my laptop too, rather than Windows).

Hence I'm almost certainly looking for source code in either C, C++ or
Fortran. Anything that works under Linux would almost certainly be
able to be compiled for Solaris without too much effort.

FORTRAN would be the language of choice (since that's what NEC was
written *and validated* in.. one would be concerned about a C
translation, although I'm sure there are C versions out there which have
been validated)

Good point.

It then tries to link with the 'blas', 'atlas' and 'lapack_atlas'
libraries, none of which my Sun has.

There should be versions out there that don't link with the matrix math
packages.


The Sun library 'libsunperf' has all the functions of blas, atlas and
lapack_atlas (well at least alls those used by NEC). I simply needed to
link against that one library, rather than the other 3, and I soon had a
n executable. I've not used it yet, as I have more pressing things to do.

That sun library should be highly optimised for the UltraSPARC
processors in my workstation.


Anyway, if you have any recommendations for the best freely available
Unix/Linux code, I would be interested.

What you've got is probably as good as anything else, especially if
you're just looking for a table of Z vs length and diameter.



Thank you for that.

Dave writes:

Since you clearly know more about this stuff than me, do you know of
the best freely available software for this which works under Unix? (I
use Sun's Solaris for 99% of the things I do, including sending this
message. I use Solaris on my laptop too, rather than Windows).

Hence I'm almost certainly looking for source code in either C, C++ or
Fortran. Anything that works under Linux would almost certainly be
able to be compiled for Solaris without too much effort.
I found this page:

NEC-2 can be built for Linux, and binaries are available for some
distributions. I am sure NEC-4 also can, if you get a license.

However, I have seen strange results come out of NEC-2 on Debian
Linux. Looked like numerical instability, and could be caused by a
problem in the toolchain. I don't think these binaries are used much,
so they don't receive the amount of tender loving care that they
deserve.

73
LA4RT Jon, Trondheim, Norway