Hello,(adsbygoogle = window.adsbygoogle || []).push({});

I have two different discrepancies to this system:

a) How and when is possible to have more solution of differential eq. or their system for same initial problem? For example this is happening in following system. It is written about this system:

"Different value of constant [itex]\dot{M}=4\pi r^2u\rho=const.[/itex] lead to physically distinct classes of solution for the same boundary condition at infinity". Can anybody explain it to me?

[itex]\frac{\rho'}{\rho}+\frac{u'}{u}+\frac{2}{r}=0[/itex]

[itex]uu'+a^2(r)\frac{\rho'}{\rho}+\frac{\alpha}{r^2}=0[/itex], where [itex]a(r)=a_0\big(\frac{\rho(r)}{\rho_0}\big)^{(\Gamma-1)/2}[/itex], [itex]\rho(r)[/itex] and [itex]u(r)[/itex] are function ofrand [itex]\Gamma, a_0, \rho_0, \alpha[/itex] are constant.[itex][/itex]

b) How to solve this system? I am quit sure that is possible just numerically and I tried. I expressed this equation in this way:

[itex]u'=\frac{D_1}{D}[/itex]

[itex]u'=\frac{D_2}{D}[/itex], where [itex]D_1=\frac{2a^2/r-\alpha/r^2}{\rho}[/itex], [itex]D_2=\frac{2u^2/r-\alpha/r^2}{u}[/itex] and [itex]D=\frac{u^2-a^2}{\rho u}[/itex].

No I see if I want smooth solution I need [itex]D_1=D_2=D=0[/itex] in same critical point [itex]r_c[/itex]. These are two boundary condition (Im not sure what difference is between boundary and initial condition?). I find out [itex]r_c[/itex] and [itex]u_c[/itex], value of [itex]u(r_c)=u_c[/itex]. So I have to condition [itex]a(r_c)=u(r_c)[/itex] and [itex]u(r_c)=u_c[/itex] Am I right till this time?

Now is coming my question:

Im not much familiar with numerical methods or numerical solving so I used Mathematica soft, concretely function NDSolve with conditions mentioned above. And Mathematica said "no". By "no" I mean it diverges in [itex]r_c[/itex]. What could be wrong?

**Physics Forums | Science Articles, Homework Help, Discussion**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Different solution for same initial (boundary) condition and numerical solution with critical point

Loading...

**Physics Forums | Science Articles, Homework Help, Discussion**