Can a basic knowledge of perturbation theory solve this?

[ozone]

Hello all,

I have boiled a very long physics problem down to the point that I need to solve the coupled equations

\frac{\partial^2 x}{\partial u^2} + xf(u) + yg(u) = 0

\frac{\partial^2 y}{\partial u^2} + yf(u) - xg(u) = 0

We may assume that|f| ,|g| << 1. and that both f and g are periodic on the same interval T, i.e. f(u) = f(u+T),g(u) = g(u+T)
I was wondering if this is something that could be solved with little knowledge of perturbation theory or if this was a very advanced problem, and I was wondering if anyone had any useful resources. So far I have found that this would be relatively easy to solve if they weren't coupled, but this is a level of complexity which is beyond me

Thank you.

 

[mfb]

Set a=x+y, b=x-y
This should decouple your system.

 

[ozone]

I tried this substitution but I'm very certain it will not actually uncouple these ode's, unless I am missing a very interesting algebra trick

 

[mfb]

Oh sorry, I missed the minus sign in the second equation.
Hmm... it is impossible to decouple them with linear transformations.

Edit: Oh, complex numbers are great.

 

[pasmith]

Hello all,

I have boiled a very long physics problem down to the point that I need to solve the coupled equations

\frac{\partial^2 x}{\partial u^2} + xf(u) + yg(u) = 0

\frac{\partial^2 y}{\partial u^2} + yf(u) - xg(u) = 0

We may assume that|f| ,|g| << 1. and that both f and g are periodic on the same interval T, i.e. f(u) = f(u+T),g(u) = g(u+T)
I was wondering if this is something that could be solved with little knowledge of perturbation theory or if this was a very advanced problem, and I was wondering if anyone had any useful resources. So far I have found that this would be relatively easy to solve if they weren't coupled, but this is a level of complexity which is beyond me

Thank you.

If you set z = x + iy and (EDIT:) w = f - ig then your system becomes
<br /> \frac{d^2 z}{d u^2} + w(u)z = 0<br />
which may be easier to analyze.

 

[bigfooted]

If you set z = x + iy and w = f + ig then your system becomes
<br /> \frac{d^2 z}{d u^2} + w(u)z = 0<br />
which may be easier to analyze.

Depending on w(u), it may still be difficult to find the general solution.
If this second order ODE admits Liouvillian solutions, they can be found using Kovacic' algorithm, you can do this with e.g. Maple.

 

[the_wolfman]

Hello all,

I have boiled a very long physics problem down to the point that I need to solve the coupled equations

\frac{\partial^2 x}{\partial u^2} + xf(u) + yg(u) = 0

\frac{\partial^2 y}{\partial u^2} + yf(u) - xg(u) = 0

We may assume that|f| ,|g| &lt;&lt; 1. and that both f and g are periodic on the same interval T, i.e. f(u) = f(u+T),g(u) = g(u+T)
I was wondering if this is something that could be solved with little knowledge of perturbation theory or if this was a very advanced problem, and I was wondering if anyone had any useful resources. So far I have found that this would be relatively easy to solve if they weren't coupled, but this is a level of complexity which is beyond me

Thank you.

Yes you can probably use perturbation theory, but you have to first define your ordering. Its not enough that f and g are small. They have to be small compared to something.

As an example let's assume xg /\frac{\partial^2 x}{\partial u^2} is order \epsilon small. We also assume that x and y are of comparable magnitude, as are f and g. The validity of these assumptions depend on the problem, and you might have to alter the following approach depending on the problem. And its important to check that your final answer is consistent with this ordering.

With these assumptions we next expand x (and y) in terms of order \epsilon: x=\epsilon^0 x_0 + \epsilon^1 x_1 + \epsilon^2 x_2 + \dots. Here \epsilon is just an marker of smallness. We also note that \frac{\partial^2 }{\partial u^2} is order 1/\epsilon large and and it is helpful to write it as \frac{1}{\epsilon}\frac{\partial^2 }{\partial u^2}

Now we plug the expansions for x and y into the equation, and equate terms of the same order in \epsilon. This creates a hierarchy of equation that we can solve.

To lowest order we have:
\frac{\partial^2 x_0}{\partial u^2} = 0

\frac{\partial^2 y_0}{\partial u^2} = 0

The next order equation is
\frac{\partial^2 x_1}{\partial u^2} + x_0f(u) + y_0 g(u) = 0

\frac{\partial^2 y_1}{\partial u^2} + y_0f(u) - x_0 g(u) = 0

And all higher order equations are of this form:
\frac{\partial^2 x_j}{\partial u^2} + x_{j-1}f(u) + y_{j-1} g(u) = 0

\frac{\partial^2 y_j}{\partial u^2} + y_{j-1}f(u) - x_{j-1} g(u) = 0

Note that in each equation for x_j,y_j the terms x_{j-1},y_{j-1} are known from the previous equation.