Solving an Object Falling with Non-Linear Power Series ODEs

[Euler1707]

Homework Statement

I have not had luck in finding a solution that describes an object falling. Forces include gravitational force which is constant and a vicous force directly proportional to the cube of the velocity. I am supposed to find v as a function of time.

Homework Equations

v' + a(v)^3 = b where a and b are constants.

The Attempt at a Solution

I have a lot of paper crumbled on the floor. I might add that I am just doing this for fun. I am practicing and learning some math. It is my hobby. :)

 

[Steely Dan]

If you are attempting a power series solution, why don't you tell us what your guess for the series is, and what the ODE looks like after substituting in the series?

 

[Euler1707]

If you are attempting a power series solution, why don't you tell us what your guess for the series is, and what the ODE looks like after substituting in the series?

Well,

I have attempted to guessv = \Sigmac_{n}t ^{n}

But my problem is the cubed term. Series solutions seem to only be easy when dealing with polynomial coefficients, but this time the dependent variable is causing the non-linearity of the ODE.

I have ofcourse substituted the derivative terms of the power series but I do not know how to cube the series in order to then find the coefficients. Cubing the series would result in a mess.

Thanks for any input.By the way this problem is from an OLD book in mechanics. Very cool read if anyone is interested

Slater and Frank Mechanics 1947.

 

[Dickfore]

Your equation should be:
<br /> m \, \dot{v} = m \, g - k \, v^3<br />

It is a 1st order ODE in velocity! The one you have is 2nd order. If you introduce dimensionless variables for time and velocity:
<br /> t = a \, x, v = b \, y, \ y = y(x)<br />
then, the equation reduces to:
<br /> m \, \frac{b}{a} \, y&#039;(x) = m \, g - k \, b^3 \, y^3<br />
Let us choose a, and b so that:
<br /> m \, \frac{b}{a} = m \, g = k \, b^3<br />
<br /> a = \left( \frac{m}{k \, g^2} \right)^{\frac{1}{3}}, \ b = \left( \frac{m \, g}{k} \right)^{\frac{1}{3}}<br />
then the ODE simplifies in form:
<br /> y&#039; = 1 - y^3<br />
This equation is with separable variables:
<br /> \frac{dy}{1 - y^3} = dx<br />
Integrating:
<br /> \int_{y_0}^{y}{\frac{dy&#039;}{1 - (y&#039;)^3}} = x, \ y(x = 0) = y_0<br />
Do the partial fraction decomposition:
<br /> \frac{1}{1 - (y&#039;)^3} = \frac{1}{(1 - y&#039;) (1 + y&#039; + (y&#039;)^2)} = \frac{A}{1 - y&#039;} + \frac{B \, y&#039; + C}{1 + y&#039; + (y&#039;)^2}<br />
<br /> 1 = A (1 + y&#039; + (y&#039;)^2) + (B \, y&#039; + C)(1 - y&#039;)<br />
<br /> 1 = (A + C) + (A + B - C) \, y&#039; + (A - B) \, (y&#039;)^2<br />
<br /> \left\lbrace\begin{array}{lcl}<br /> A + C &amp; = &amp; 1 \\<br /> A + B - C &amp; = &amp; 0 \\<br /> A - B &amp; = &amp; 0<br /> \end{array} \right. \Leftrightarrow A = \frac{1}{3} \, B = \frac{1}{3} \, C = \frac{2}{3}<br />
Thus, you have the integrals:
<br /> \frac{1}{3} \, \int_{y_0}^{y}{\frac{dy&#039;}{1 - y&#039;}} + \frac{1}{3} \, \int_{y_0}^{y}{dy&#039; \, \frac{y&#039; + 2}{1 + y&#039; + (y&#039;)^2}} = x<br />
Do the integrals (assuming 0 \le y, \ y_0 &lt; 1

 

[Euler1707]

Thank you

Your solution is very elegant. Although the book calls for a power series solution ( at least first 4 terms). I was reviewing from my ODE books but like always math books turn out not to have the thing you are looking for. hehe

Nevertheless thanks for taking the time to solve it.

Also yeah, I did realized that I got carried away with the primes. Ofcourse the physical description was enough to figure it out.

 

[Euler1707]

Oh, I think I might have found a better method that involves a polynomial solution. Maybe just by iteration using the Taylor series. The problem does not say anything about I.C s but I cannot think of anything else. xD