Help with differential equation

In summary: I am not sure how you got that formula, but I will take your word for it. So, I am thinking I should go back and try to solve it by hand, but I am not sure I have the time or desire to do so.
  • #1
PhyPsy
39
0

Homework Statement


[itex]\frac{dr}{\sqrt{r^2- 2mr}}= \frac{d\rho}{\rho}[/itex]
Solve for r.

Homework Equations





The Attempt at a Solution


For those familiar with GR, this is an attempt to come up with the isotropic form of the Schwarzschild solution to the field equations, but I'm having trouble with the math, so I just need someone to check my math. The solution I am supposed to get is [itex]r=\rho(1+ \frac{1}{2}m/\rho)^2[/itex], but I am coming up with [itex]r=\rho(\frac{1}{2}+ m/\rho)^2[/itex].

After integration, the equation I came up with is [itex]2\ln(\sqrt{r-2m}+ \sqrt{r})= \ln\rho[/itex]
[itex]e^{2\ln(\sqrt{r-2m}+ \sqrt{r})}= e^{\ln\rho}[/itex]
[itex]e^{\ln(\sqrt{r-2m}+ \sqrt{r})+ \ln(\sqrt{r-2m}+ \sqrt{r})}= e^{\ln\rho}[/itex]
[itex]e^{\ln(\sqrt{r-2m}+ \sqrt{r})}e^{\ln(\sqrt{r-2m}+ \sqrt{r})}= e^{\ln\rho}[/itex]
[itex](\sqrt{r-2m}+ \sqrt{r})^2= \rho[/itex]
[itex]r- 2m+ 2\sqrt{r^2- 2mr}+ r= \rho[/itex]

I moved all terms on the left side of the equal sign to the right except for the square root term and squared both sides.
[itex]4(r^2- 2mr)= (\rho+ 2m- 2r)^2= \rho^2+ 4m\rho- 4r\rho+ 4m^2- 8mr+ 4r^2[/itex]

The terms [itex]4r^2[/itex] and [itex]-8mr[/itex] cancelled, and then I moved the term [itex]-4r\rho[/itex] to the left side of the equal sign.
[itex]4r\rho= \rho^2+ 4m\rho+ 4m^2[/itex]
[itex]r\rho= (\frac{1}{2}\rho+ m)^2[/itex]
[itex]r= \rho^2(\frac{1}{2}+ m/\rho)^2/\rho= \rho(\frac{1}{2}+ m/\rho)^2[/itex]

Does anyone see a problem?
 
Physics news on Phys.org
  • #2
PhyPsy said:

Homework Statement


[itex]\frac{dr}{\sqrt{r^2- 2mr}}= \frac{d\rho}{\rho}[/itex]
Solve for r.

Homework Equations





The Attempt at a Solution


For those familiar with GR, this is an attempt to come up with the isotropic form of the Schwarzschild solution to the field equations, but I'm having trouble with the math, so I just need someone to check my math. The solution I am supposed to get is [itex]r=\rho(1+ \frac{1}{2}m/\rho)^2[/itex], but I am coming up with [itex]r=\rho(\frac{1}{2}+ m/\rho)^2[/itex].

After integration, the equation I came up with is [itex]2\ln(\sqrt{r-2m}+ \sqrt{r})= \ln\rho[/itex]
[itex]e^{2\ln(\sqrt{r-2m}+ \sqrt{r})}= e^{\ln\rho}[/itex]
[itex]e^{\ln(\sqrt{r-2m}+ \sqrt{r})+ \ln(\sqrt{r-2m}+ \sqrt{r})}= e^{\ln\rho}[/itex]
[itex]e^{\ln(\sqrt{r-2m}+ \sqrt{r})}e^{\ln(\sqrt{r-2m}+ \sqrt{r})}= e^{\ln\rho}[/itex]
[itex](\sqrt{r-2m}+ \sqrt{r})^2= \rho[/itex]
[itex]r- 2m+ 2\sqrt{r^2- 2mr}+ r= \rho[/itex]

I moved all terms on the left side of the equal sign to the right except for the square root term and squared both sides.
[itex]4(r^2- 2mr)= (\rho+ 2m- 2r)^2= \rho^2+ 4m\rho- 4r\rho+ 4m^2- 8mr+ 4r^2[/itex]

The terms [itex]4r^2[/itex] and [itex]-8mr[/itex] cancelled, and then I moved the term [itex]-4r\rho[/itex] to the left side of the equal sign.
[itex]4r\rho= \rho^2+ 4m\rho+ 4m^2[/itex]
[itex]r\rho= (\frac{1}{2}\rho+ m)^2[/itex]
[itex]r= \rho^2(\frac{1}{2}+ m/\rho)^2/\rho= \rho(\frac{1}{2}+ m/\rho)^2[/itex]

Does anyone see a problem?

I see one. The integration will give you [itex]2\ln(\sqrt{r-2m}+ \sqrt{r})= \ln\rho +C[/itex]. There's an undetermined constant, C. I don't think you can assume it's zero. I think a choice of that constant will give you the books answer. Don't you have some way to determine the constant?
 
  • #3
Doing the math in my head, the constant would have to be [itex]\ln(\rho-m)[/itex]. The constant from integration of the left side of the equation allows for the rho term. I risk leaving the topic of this forum by asking this, but every book and article I have seen uses the specific transformation [itex]r=\rho(1+\frac{1}{2}m/\rho)^2[/itex] when deriving the isotropic form of the Schwarzschild solution. Why could I not use 0 as the constant of integration and have the transformation be [itex]r=\rho(\frac{1}{2}+m/\rho)^2[/itex]? It seems to me that also leads to a legitimate isotropic construct of the Schwarzschild solution.
 
  • #4
PhyPsy said:
Doing the math in my head, the constant would have to be [itex]\ln(\rho-m)[/itex]. The constant from integration of the left side of the equation allows for the rho term. I risk leaving the topic of this forum by asking this, but every book and article I have seen uses the specific transformation [itex]r=\rho(1+\frac{1}{2}m/\rho)^2[/itex] when deriving the isotropic form of the Schwarzschild solution. Why could I not use 0 as the constant of integration and have the transformation be [itex]r=\rho(\frac{1}{2}+m/\rho)^2[/itex]? It seems to me that also leads to a legitimate isotropic construct of the Schwarzschild solution.

[itex]\ln(\rho-m)[/itex] isn't a constant. You can't add that. You generally have some boundary conditions you want your coordinates to satisfy. In the books solution for large [itex]\rho, r[/itex] you have [itex]\rho \approx r[/itex]. Your version wouldn't have that.
 
  • #5
you could try differentiating your answer to see what factors got changed then go back thru your steps to see if you properly accounted for them.

in undergrad we once did an eperiment to measure the speed of light via intereference and a rotating magnetic and got 1/2 the seepd of light. We never figured out why and the prof couldn't find anything wrong with what we did. After many years, I realized we may not have factered in that the mirror was double sided which meant it effectively spun at twice the speed. Oh well, that's why I went into programming that and the money even though physics was way more interesting.
 
  • #6
PhyPsy said:
After integration, the equation I came up with is [itex]2\ln(\sqrt{r-2m}+ \sqrt{r})= \ln\rho[/itex]
how did you get this? Seems you should have had
[tex]
\ln(\sqrt{(r-m)^2-m^2} + r-m) = \ln(C \rho)
[/tex]
choosing [itex]C=2[/itex] you get your answer
 
  • #7
Dick said:
[itex]\ln(\rho-m)[/itex] isn't a constant. You can't add that. You generally have some boundary conditions you want your coordinates to satisfy. In the books solution for large [itex]\rho, r[/itex] you have [itex]\rho \approx r[/itex]. Your version wouldn't have that.
Bah, you're right, since [itex]\rho[/itex] is a function of r, I can't use it as a constant. My bad.

jedishrfu said:
you could try differentiating your answer to see what factors got changed then go back thru your steps to see if you properly accounted for them.

in undergrad we once did an eperiment to measure the speed of light via intereference and a rotating magnetic and got 1/2 the seepd of light. We never figured out why and the prof couldn't find anything wrong with what we did. After many years, I realized we may not have factered in that the mirror was double sided which meant it effectively spun at twice the speed. Oh well, that's why I went into programming that and the money even though physics was way more interesting.
Ha ha, I'm actually in IT myself, though I am not a hardcore programmer; I do some here and there. Physics is a very interesting subject to me, but a career in it, on the other hand, doesn't appeal to me too much.

quZz said:
how did you get this? Seems you should have had
[tex]
\ln(\sqrt{(r-m)^2-m^2} + r-m) = \ln(C \rho)
[/tex]
choosing [itex]C=2[/itex] you get your answer
Well, to be honest, I cheated and used Wolfram to compute the integral for the left side. I thought it might be a bad idea, but I looked into using substitution and integration by parts and couldn't get anywhere using those techniques. Yes, differentiating your solution does indeed return [itex]\frac{dr}{\sqrt{r^2-2mr}}[/itex]. I am curious to know how you came up with that.

Also, I was not aware you would get a constant multiplied by [itex]\rho[/itex] on the right side. I don't remember learning about that in calculus, but it has been 7 years since I took it, so maybe I just completely forgot.
 
  • #8
Ok, first you notice that [itex]r^2-2mr = (r-m)^2-m^2[/itex] and then you substitute [itex]x\equiv r-m[/itex], and [itex]dr=dx[/itex]. The integral then becomes
[tex]
\int\frac{dx}{\sqrt{x^2-m^2}} = \ln(\sqrt{x^2-m^2}+x)
[/tex]
Second, on the right hand side you have
[tex]
\int\frac{d\rho}{\rho}=\ln\rho + C_1
[/tex]
Let's write [itex]C_1=\ln C[/itex] so that we have pure [itex]\ln[/itex]s on both sides
[tex]
\ln\rho + C_1 = \ln(C\rho)
[/tex]
 

FAQ: Help with differential equation

What is a differential equation?

A differential equation is a mathematical equation that describes the relationship between a function and its derivatives. It involves one or more independent variables and the derivatives of the dependent variable with respect to those independent variables.

Why are differential equations important?

Differential equations are important because they are used to model and analyze a wide range of phenomena in various fields such as physics, chemistry, biology, economics, and engineering. They provide a mathematical framework for understanding and predicting the behavior of complex systems.

How do I solve a differential equation?

The method for solving a differential equation depends on the type of equation. Some common techniques include separation of variables, substitution, and integration. It is important to first identify the type of equation and then choose an appropriate method to solve it.

Can differential equations be solved analytically?

Some differential equations can be solved analytically, meaning an explicit formula or expression can be found for the solution. However, there are many differential equations that do not have analytical solutions and require numerical or approximate methods for solving.

How are differential equations used in real-life applications?

Differential equations are used extensively in real-life applications to model and understand various phenomena, such as population growth, heat transfer, chemical reactions, and electrical circuits. They are also used in engineering fields to design and optimize systems and in economics to study economic trends and behavior.

Back
Top