Differential equation system, got stuck in a physics problem

[ramses728]

hi guys,
I'm solving a pretty complex problem: calculating a trajectory of a charged particle in a custom magnetic field. I arrive to the point where this very nice equation system blocks my way :P

http://ramses728.altervista.org/img/phys.jpg

hope you can help me somehow

 

[dextercioby]

It's not a valid link, there's no picture in it. I'm using Firefox as a webbrowser.

 

[ramses728]

Strange... But if you copy and paste it, it works 100%

 

[CompuChip]

As I understand it, the equation is
m x'''(t) = \frac{q z''(t) y(t)}{r(t)}, m y'''(t) = \frac{q z''(t) x(t)}{r(t)}, m z'''(t) = \frac{ q\left( x'(t) x(t) + y'(t) y(t) \right) }{ r(t) }
(the last y'''(t) in the original should be a z'''(t) I presume), where
r(t) = \sqrt{ x(t)^2 + y(t)^2.

This looks very complicated. Are you sure it is correct? I would either try polar or spherical coordinates (especially if there is some kind of symmetry in the system), or otherwise to first solve x(r(t)) or at least the functions in a more convenient variable than t. But actually, I have no idea

 

[ramses728]

Yeah it's complicated indeed... and you got that right with z'''(t)

I actually know what it should look like when it's solved (graphically i mean) and there is some pretty serious symmetry in the system... I'll give a shot with the cilindrical coordinates, but there will be much work at it especially because I have to reconstruct the field function with polar coordinates (mabye esier... who knows...)

anyway if someone has any ideas how to solve this (also with maple or similar), I'd relly apreciate it

as soon as i get the cilindrical set of equations, I'll post them.

thanks, ramses

Edit: ah yeah forgot to say it's an equation system, so three equations together

 

[HallsofIvy]

How did you get third derivatives? Kinematic equations always involve acceleration, the second derivative.

 

[ramses728]

... you're right... i made some confusion added a derivation grade for acceleration and velocity XD ... so, here are the correct ones:

<br /> \left\{<br /> \begin{array}[l]{l}<br /> m x(t)&#039;&#039; = q z(t)&#039; \frac{x(t)}{r(t)}\\<br /> m y(t)&#039;&#039; = q z(t)&#039; \frac{y(t)}{r(t)}\\<br /> m z(t)&#039;&#039; = -\frac{q}{r(t)}(x(t)&#039; x(t) + y(t)&#039; y(t))\\<br /> \end{array}<br />

where

r(t)=\sqrt{x(t)^2 + y(t)^2}

still working on those polar ones

thanks for your interest

(P.S. I finally realized that fourmulaes can be written with latex! XD)

 

[Matthew Rodman]

So, as

r^{\prime}(t) = \frac{1}{2} \frac{(2 x(t) x^{\prime}(t) + 2 y(t) y^{\prime}(t))}{\sqrt{x^2(t) + y^2(t)}} = \frac{x(t)x^{\prime}(t) + y(t) y^{\prime}(t)}{\sqrt{x^2(t) + y^2(t)}}

You can write

m z^{\prime \prime}(t) = -q r^{\prime}(t)

Hence, for z you get

z^{\prime}(t) = \kappa - \frac{q}{m} r(t)

where \kappa is a constant.

Now, for \kappa = 0, the solution is straightforward. You will then get

x^{\prime \prime} = -(\frac{q}{m})^2 x(t)

and

y^{\prime \prime} = -(\frac{q}{m})^2 y(t)

But in the more general case, \kappa may not be taken as 0, so you'll get

x^{\prime \prime}(t) = \frac{q}{m} x(t) (\frac{\kappa}{r(t)} - \frac{q}{m})

and

y^{\prime \prime}(t) = \frac{q}{m} y(t) (\frac{\kappa}{r(t)} - \frac{q}{m})

 

[ramses728]

you got it! Well thanks very much!

just for curiosity are you a math student?

 

[Matthew Rodman]

I've only solved it for \kappa =0. For \kappa \ne 0, it is still a nasty mess.

And, no, I'm not a student.

 

[ramses728]

It's still a great advance, because I can set the initial conditions for z velocity (so z(t)&#039;) to be 0 forcing the constant to be zero. So I have the problem solved for this specifical case. But as you say it still needs work for the more general case...

Hope you didn't take me thinking you were a student as offensive, I just thought that you were closely involved on math so supposed that. If any offense was taken I apologise. Thanks again for this partial solution

 

[Matthew Rodman]

No worries, mate.