CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

hankel negative order

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 16, 2006, 06:06
Default hankel negative order
  #1
sankar
Guest
 
Posts: n/a
HI all, I am using hankel functions, and come across a probelm in Matlab while applying for negative order. eg: besselh(-7,0.05)

ans =

-2.2976e-003 +3.7553e+013i

But i found some error when checked using the formula given in Abramowitz & Stegun (class 9.16) which is the source for this H(-n,z) = exp(n*pi*i)*H(n,z) where n is order and z is argument

>> exp(7*pi*i)*besselh(7,0.05)

ans =

2.9895e-002 +3.7553e+013i

even though the imaginary part is same..there is some error in real part.For certain orders this is working good in matlab. can any one comment on this which one to use.. Is there any fortran code to handle negative order

sankar.
  Reply With Quote

Old   June 16, 2006, 07:12
Default Re: hankel negative order
  #2
Tom
Guest
 
Posts: n/a
A play with a calculator should make you suspicious since

H_n(x) = J_n(x) + iY_n(x)

Now for n=7 and x=0.05 we expect J_n(x)~1.21e-15 and Y_n(x)~-3.755e+13 (lead order term in asympotic series as x->0 - which is a good approximation because of the relative largeness of n).

Now exp(i.n.pi) = -1 so H_n(x)~-1.21e-15 + i.3.755e+13.

Tom.
  Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
1st order vs 2nd order ken FLUENT 8 March 14, 2013 04:43
channelFoam for a 3D pipe AlmostSurelyRob OpenFOAM 3 June 24, 2011 14:06
blockMesh error ... balkrishna OpenFOAM Pre-Processing 0 August 17, 2010 03:39
[blockMesh] BlockMesh FOAM warning gaottino OpenFOAM Meshing & Mesh Conversion 7 July 19, 2010 15:11
[blockMesh] Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Meshing & Mesh Conversion 10 April 2, 2007 15:00


All times are GMT -4. The time now is 07:11.