What is the correct method for solving a homogeneous ODE with substitution?

  • Thread starter Beer-monster
  • Start date
  • Tags
    Ode
In summary, Beer thinks you might be missing an x in the first expression. The last step is not appropriate. You can solve it using the method described.
  • #1
Beer-monster
296
0
Okay, so I'm trying to work out a streamline of a 2D fluid, I apply the usual relation ship between the diffential and the flow velocities and get the following differential equation.
[tex] \frac {dy } {dy} = \frac {(y^2-x^2)} {2xy} [/tex]
I know this is a homogeneous equation (the lecturer included that handy tip), so I try the substitustion y = xv. Play with the equation until I get end up with
[tex] \frac {dv} {dx} = -\frac {v} {2} - \frac {1} {v} [/tex]
Which my best attempts to separate gets me (note I should have thought about this tex stuff first:cry: )
[tex] \ln y = -\frac {y^2} {2x^2} + C [/tex]
Which is not the given answer. Any idea where I'm going wrong. Is it my method of solving the equation? Or did I make the wrong substitution? Is the degree of the equation different from what I though (and what do they mean by degree anyway can anyone explain...I seem to have forgottten:blushing: )
Thanks for any help
~Beery
 
Last edited:
Physics news on Phys.org
  • #2
I think you might get a different result if you tried it again. I would guess that the problem is 'careless' errors in your algebra.
 
  • #3
Beer-monster said:
Okay, so I'm trying to work out a streamline of a 2D fluid, I apply the usual relation ship between the diffential and the flow velocities and get the following differential equation.
[tex] \frac {dy } {dy} = \frac {(y^2-x^2)} {2xy} [/tex]
I know this is a homogeneous equation (the lecturer included that handy tip), so I try the substitustion y = xv. Play with the equation until I get end up with
[tex] \frac {dv} {dx} = -\frac {v} {2} - \frac {1} {v} [/tex]
Which my best attempts to separate gets me (note I should have thought about this tex stuff first:cry: )
[tex] \ln y = -\frac {y^2} {2x^2} + C [/tex]
Which is not the given answer. Any idea where I'm going wrong. Is it my method of solving the equation? Or did I make the wrong substitution? Is the degree of the equation different from what I though (and what do they mean by degree anyway can anyone explain...I seem to have forgottten:blushing: )
Thanks for any help
~Beery

So y=xv and when I substitute that into the DE, I get:

[tex](x^2v^2-x^2)dx-2x(xv)(xdv+vdx)=0[/tex]

Can you simplify that down to:

[tex](v^2-1-2v^2)dx-2xvdv=0[/tex]

and then re-arrange it to:

[tex]\frac{dx}{x}+\frac{2v}{v^2+1}dv=0[/tex]

?
 
  • #4
Um sorta. The second part was a bit fiddly but I think I got there.
Your method is a little unfamilar took some sorting.
Subbing for y=xv and working through using the method I'm more familar with I get something similiar
[tex] \frac{dv}{dx} +v = \frac{(x^2v^2-x^2)}{2x^2v}[/tex]
Cancelling for x^2

[tex] \frac{dv}{dx} +v = \frac{v^2-1}{2v} [/tex]

This is the bit I'm not sure about. I deduct v from both side then invert the whole thing to get.

[tex] \frac{1}{x} \frac{dx}{dv} = \frac{2v}{v^2 -1} - \frac {1}{v} [/tex]
Which I'm not sure is correct (and I also -despite many attempts- never got the hang of partial fractions)
Does this make sense?
 
Last edited:
  • #5
Beer-monster said:
Um sorta. The second part was a bit fiddly but I think I got there.
Your method is a little unfamilar took some sorting.
Subbing for y=xv and working through using the method I'm more familar with I get something similiar
[tex] \frac{dv}{dx} +v = \frac{(x^2v^2-x^2)}{2x^2v}[/tex]
Cancelling for x^2
[tex] \frac{dv}{dx} +v = \frac{v^2-1}{2v} [/tex]
This is the bit I'm not sure about. I deduct v from both side then invert the whole thing to get.
[tex] \frac{1}{x} \frac{dx}{dv} = \frac{2v}{v^2 -1} - \frac {1}{v} [/tex]
Which I'm not sure is correct (and I also -despite many attempts- never got the hang of partial fractions)
Does this make sense?

Beer, not sure what you're doing. Seem to be missing an x in the first expression. That last step is not appropriate. This is the correct way:

[tex]\frac{dy}{dx}=\frac{y^2-x^2}{2xy}[/tex]

so multiplying by 2xydx leads to:

[tex]2xydy=(y^2-x^2)dx[/tex]

or:

[tex](y^2-x^2)dx-(2xy)dy=0\quad\tag{1}[/tex]

Now let y=xv so that:

[tex]dy=xdv+vdx[/tex]

Now just substitute this and the expression y=xv into (1) to get:

[tex](x^2v^2-x^2)dx-2x(xv)(xdv+vdx)=0[/tex]

You can simplify this right? Then just divide by [itex]x^2[/itex] to obtain:

[tex]v^2dx-dx-2xvdv-2v^2dx=0[/tex]

rearrage:

[tex](v^2-1-2v^2)dx-2xvdv=0[/tex]

divide by:

[tex]x(v^2-1-2v^2)[/tex]

results in:

[tex]\frac{dx}{x}+\frac{2v}{v^2+1}dv=0[/tex]

You can finish it. Please go through each step above by hand to make sure you understand it then since I did yours, you do this one:

[tex]\frac{dy}{dx}=\frac{xy}{x^2+3y^2}[/tex]
 
Last edited:
  • #6
Well I tried your problem and my maths must be much suckier than I thought. I followed your method as best I could and got something stupid:cry: i.e.

[tex] -\frac{1}{3v} + v = \ln{x} + C [/tex]

Which is probably miles off.:confused:

I probably chose the wrong substitution or divided by the wrong factor somewhere.

Times like these (when I spend a few hours on a problem and still can't solve it) I get quite stressed with myself.

*sigh*

Thanks for your help
 
  • #7
Beer-monster said:
Well I tried your problem and my maths must be much suckier than I thought. I followed your method as best I could and got something stupid:cry: i.e.
[tex] -\frac{1}{3v} + v = \ln{x} + C [/tex]
Which is probably miles off.:confused:
I probably chose the wrong substitution or divided by the wrong factor somewhere.
Times like these (when I spend a few hours on a problem and still can't solve it) I get quite stressed with myself.
*sigh*
Thanks for your help

Beer, sorry I take so long to get back with you. I seem to be on a different schedule as you. This is important:

Did you follow my example above in detail and understand every part of it? If so, can you then not look at my solution and solve the problem (up to the first order in v) by yourself? You can do that right? If not, where are you having problems? First get that one straight ok? Now for the second one:

Using y=xv, I get:

[tex](x^2+3v^2x^2)(vdx+xdv)=x^2vdx[/tex]

right?

Group them together:

[tex]x^2vdx+x^3dv+3v^3x^2dx+3v^2x^3dv=x^2vdx[/tex]

divide by x^2:

[tex]vdx+xdv+3v^3dx+3v^2xdv=vdx[/tex]

group the dv's and dx's:

[tex](v+3v^3-v)dx+(x+3v^3x)dv=0[/tex]

or:

[tex]3v^3dx+x(1+3v^2)dv=0[/tex]

divide by:

[tex]x3v^3[/tex]

[tex]\frac{dx}{x}+\frac{1+3v^2}{3v^3}dv=0[/tex]

Can you finish it?
 
Last edited:
  • #8
Yeah I tried it now and got the method, at first I wasn't quite sure how you got the divisions you did. But after trying it myself I got the same thing, but in more steps (slower but more intuitive for me).


My main problem then was the integration. Luckily it was tutorial session and my lecturer pointed out that the numerator was the differential of the denominator thus the integral was just the log of the denominator. That had honestly never occurred to me :blushing:

Don't you just hate it when the answer is really simple but your brain just never saw it *headdesks* ?

It was the same with your problem, I followed the method but couldn't solve the integral. Does the same apply, since the numerator is related to the differential of the denomonator, does that make the answer just the integral of the donminator.

Should I check for this whenever I get a similar function (i.e. an algebraic fraction I can't factorise)?
 
  • #9
If the numerator is the differntial of the denominator that is correct. But that's not the case with the v's on the RHS of this:

[tex]\frac{dx}{x}+\frac{1+3v^2}{3v^3}dv=0[/tex]

right?

the differential of the denominator is [itex]9v^2dv[/itex]

To integrate the first expression, I would just split the RHS up into two parts:

[tex]\frac{1+3v^2}{3v^3}=1/3v^{-3}+v^{-1}[/tex]

thus we have:

[tex]\frac{dx}{x}+(1/3v^{-3}+v^{-1})dv=0[/tex]

Integrating:

[tex]ln(x)+1/3(-1/2 v^{-2})+ln(v)=c[/tex]


letting v=y/x:

[tex]ln(y/x)-\frac{x^2}{6y^2}=c-ln(x)[/tex]

simplifying:

[tex]6y^2ln(y/c)=x^2[/tex]
 
  • #10
That had also occurred to me but it didn't work with the previous on so I didn't think it would now. I'm such a dope.

Thanks so much for all your help Saltydog, now that I got that method, I hopefully won't be as helpless with these sorts of equations next time
 

FAQ: What is the correct method for solving a homogeneous ODE with substitution?

What is an ODE?

An ODE stands for Ordinary Differential Equation. It is an equation that involves a function and its derivatives, and is used to model various physical phenomena in science and engineering.

Why are ODEs important?

ODEs are important because they allow us to mathematically describe and understand complex systems and processes in science and engineering. They are used in a wide range of fields, from physics and chemistry to economics and biology.

What are some real-world applications of ODEs?

ODEs have a wide range of real-world applications, including modeling population growth, predicting the spread of diseases, understanding chemical reactions, and designing control systems for engineering processes. They can also be used to analyze and optimize complex systems, such as in economics and finance.

What are the different types of ODEs?

There are several types of ODEs, including linear and nonlinear, first-order and higher-order, and autonomous and non-autonomous. Each type has its own unique characteristics and methods for solving.

How do you solve an ODE?

There are various methods for solving ODEs, including analytical methods such as separation of variables and power series, and numerical methods such as Euler's method and Runge-Kutta methods. The specific method used depends on the type and complexity of the ODE.

Back
Top