Solutions to 2 similar ODE's. Why does they differ that much?

  • Context: Graduate 
  • Thread starter Thread starter fluidistic
  • Start date Start date
Click For Summary
SUMMARY

The discussion focuses on the differences between the solutions of two similar ordinary differential equations (ODEs): R''+(1/r)R'-KR=0 and R''+(2/r)R'-KR=0. The first ODE's solution involves Bessel functions of order 0, specifically R(r)=c_1J_0(λr), while the second yields R(r)=c_2( sin(√C r) / r), which does not involve Bessel functions. The key takeaway is that despite superficial similarities, the underlying mathematical structures differ significantly, with spherical Bessel functions related to sin(x)/x appearing in the context of the Helmholtz equation.

PREREQUISITES
  • Understanding of ordinary differential equations (ODEs)
  • Familiarity with Bessel functions and their properties
  • Knowledge of spherical Bessel functions and their applications
  • Basic grasp of the Helmholtz equation in spherical coordinates
NEXT STEPS
  • Study the properties and applications of Bessel functions of the first kind, particularly J_0(λr)
  • Explore the derivation and solutions of the Helmholtz equation in spherical coordinates
  • Investigate the relationship between spherical Bessel functions and trigonometric functions like sin(x)/x
  • Examine the physical interpretations of solutions to different forms of ODEs in mathematical physics
USEFUL FOR

Mathematicians, physicists, and engineering students interested in differential equations, particularly those studying wave equations and vibrations in circular domains.

fluidistic
Gold Member
Messages
3,934
Reaction score
286
I'm reading on wikipedia (http://en.wikipedia.org/wiki/Vibrations_of_a_circular_drum#The_radially_symmetric_case) that the solution to the ODE [itex]R''+(1/r)R'-KR=0[/itex], by letting [itex]K=-\lambda ^2[/itex] is a linear combination of the Bessel functions of order 0, namely [itex]R(r)=c_1 J_0 (\lambda r ) + c_2 Y_0 (\lambda r)[/itex]. Then they say that c_2 must equal 0 for the solution to make physical sense (because Y_0 diverges at r=0) and thus the solution to the ODE is of the form [itex]R(r)=c_1J_0 (\lambda r )[/itex]. So far so good.
However in the problem of my thread https://www.physicsforums.com/showthread.php?t=652277, I had the equation [itex]R''+\frac{2R'}{r}+CR=0[/itex] which is very similar to the ODE in wikipedia and I've determined that the physically acceptable solution is of the form [itex]R(r)=c_2\frac{\sin (\sqrt C r )}{r}[/itex] (I've checked that it indeed satisfies the ODE). I find this unbelievable that the solution to the ODE of my problem differs so much from the one in wikipedia. What's going on here? I'm totally clueless.

Edit: Here's a graph between J_0(x) and sin (x) /x : http://www.wolframalpha.com/input/?i=graph+J_0+(x)+and+sin+(x)+/x+from+0+to+100. They are very similar! Is it possible to write the Bessel function of the first kind of order 0 as a constant times sin (kx) /x?
 
Last edited:
Physics news on Phys.org
According to http://en.wikipedia.org/wiki/Bessel_function, "The graphs of Bessel functions look roughly like oscillating sine or cosine functions that decay proportionally to 1/√x, although their roots are not generally periodic, except asymptotically for large x."
 
haruspex said:
According to http://en.wikipedia.org/wiki/Bessel_function, "The graphs of Bessel functions look roughly like oscillating sine or cosine functions that decay proportionally to 1/√x, although their roots are not generally periodic, except asymptotically for large x."

Ah thank you! I missed this part!
So basically although the solutions to [itex]R''+(1/r)R'-KR=0[/itex] and [itex]R''+(2/r)R'-KR=0[/itex] can look similar, there's still some "huge" difference in that Bessel functions are involved in the first ODE but not in the second. And that we cannot write the Bessel function of the first kind of order 0 as a constant times sin (kx) /x.
Wow, impressive.

P.S.:I'm just seeing http://en.wikipedia.org/wiki/Spherical_Bessel_function#Spherical_Bessel_functions:_jn.2C_yn. Apparently [itex]\sin x /x[/itex] is a spherical Bessel function (of order 0 I think)! That appears when solving the Helmholtz equation in spherical coordinates. Wow.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 4 ·
Replies
4
Views
4K
Replies
7
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K