2nd order Differential Equation

[xago]

Homework Statement

Ugh I feel really stupid for posting this but for some reason I can't remember how to solve it.
I am trying to solve the diff eqn: D*\frac{d^{2}\phi}{dx^{2}} - \Sigma_{a}*\phi(x) = -q_{t}
(thermal diffusion equation for neutrons slowing down to thermal energy)

The Attempt at a Solution

Anyways I've taken the roots of the LHS as if it was homogeneous:
(\phi(x)^{(2)}-\frac{\Sigma_{a}}{D}) = 0

which gives the 2 roots ±\sqrt{\frac{\Sigma_{a}}{D}}

and therefore \phi(x) = C1*e^{\sqrt{\frac{\Sigma_{a}}{D}}x} + C2*e^{-\sqrt{\frac{\Sigma_{a}}{D}}x}

so the only thing that's bothering me is the -q_{t} on the RHS which I don't know what to do with. I've used Maple to get the solution and it gives me
\phi(x) = C1*e^{\sqrt{\frac{\Sigma_{a}}{D}}x} + C2*e^{-\sqrt{\frac{\Sigma_{a}}{D}}x} + \frac{q_{t}}{\Sigma_{a}}

I would just like to know where that last term comes from and how to get it.

 

[vivek.iitd]

you can get the solution using Method of variation of parameters. Using Wronskian and all.

 

[vivek.iitd]

After you calculate the solution to the homogeneous differential equations, the two arbitrary constants are replaced by two functions, A(x) and B(x) of 'x', which calculated by finding Wronskian, and when you will put those functions into the solution obtained for the homogeneous differential equation you will get the last term.

 

[xago]

Ok I've done the problem with method of variation of parameters but I get the particular solution of \frac{D*q_{t}}{\Sigma_{a}} and not \frac{q_{t}}{\Sigma_{a}} like the solution that Maple gives me..

 

[HallsofIvy]

What, exactly, is q_t? If it is not a function of x, even if it is a function of another variable, such as t, in an equation in which the differentiation is with respect to x, it can be treated like a constant. If that is true then look for a constant "particular solution": \phi= A so that d^2\phi/dt^2= 0. The equation reduces to -\Sigma_\alpha A= -q_t which immediately gives
A= \frac{q_t}{\Sigma}

That, I presume, was how Maple interpreted this problem. If q_t is a function of x, then the solution will depend strongly on exactly what that solution is.

 

[xago]

Yes qt is a constant and your solution is much more practical, but now I'm wondering why I couldn't get the same solution with variation of parameters. It should be possible right?

Edit: I think I see now that variation of parameters cannot be used unless the non-homogenous term is a function of x, unlike here where it is just a constant.

 

[Ray Vickson]

Yes qt is a constant and your solution is much more practical, but now I'm wondering why I couldn't get the same solution with variation of parameters. It should be possible right?

Edit: I think I see now that variation of parameters cannot be used unless the non-homogenous term is a function of x, unlike here where it is just a constant.

If q_t and \Sigma_a are constants, let y = \phi - q_t/ \Sigma_a, and look at the DE for y. That is one of the standard tricks used in dealing with DEs.

RGV

 

[HallsofIvy]

Yes qt is a constant and your solution is much more practical, but now I'm wondering why I couldn't get the same solution with variation of parameters. It should be possible right?

Edit: I think I see now that variation of parameters cannot be used unless the non-homogenous term is a function of x, unlike here where it is just a constant.

No, variation of parameters works in that case- although it is much harder! We can't say why you did not get that solution with variation of parameters because you did not show what you did. When I try using variation of parameters, I get precisely q_t/\Sigma_\alpha as the "particular solution.

 

[vivek.iitd]

Hey Xago, I got the mistake which you are doing,

Look carefully, that when you are calculating Wronskian, you may be taking the g(x) to be qt, which is wrong because you will have to take it qt/D, as you have earlier divided the whole equation with D to get the characteristic solution (solution to the homogeneous part).
now if you take the g(x) to be qt/D, you will see that the D will go away automatically.