How Do You Solve a PDE with Boundary Conditions Using Separation of Variables?

[useruseruser]

u_t=u_{xx}+2u_x
0<=x<=L, t>=0, u(x,0)=f(x), u_x(0,t)=u_x(L,t)=0

How to do this?

 

[HallsofIvy]

Just a standard "separation of variables" will work.

Assume u is of the form u(x,t)= A(x)B(t). Then u_x= A'B, u_t= AB', and u_xx= A"B. Putting those into the equation, AB'= A"B+ 2A'B.

Divide both sides by AB:
\frac{B&#039;}{B}= \frac{A&quot;+2A&#039;}{A}

Since the left side depends upon t only while the right side depends upon x only, in oder to be equal they must both be equal to the same constant. Call that constant "\lambda".

Then we have B&#039;/B= \lambda or B&#039;= \lambda B and (A&quot;+ 2A&#039;)/A= \lambda or A&quot;+ 2A&#039;= \lambda A.

Now look at the boundary conditions: u_x(0,t)= A(0)B&#039;(t)= 0 for all t so we must have A(0)= 0. Also, u_x(L,t)= A(0)B&#039;(t)= 0 for all t so A(L)= 0.

Solve the boundary problem A&quot;+ 2A&#039;- \lamba A= 0, A(0)= 0, A'(0)= 0 to determine what \lambda can be. You should see that this is an eigenvalue problem. A&quot;+ 2A&#039;= \lambda A is trivially satisfied for A(x) identically equal to 0 but then we can't satisfy "u(x,0)= f(x)". The values of \lambda for which there are not trivial solutions are the "eigenvalues" for this problem and the non-trivial solutions are the "eigenvectors" or, in this case, "eigenfunctions". Once you have determined what the eigenvalues are, you can put them into B"= \lamba B to solve that.

As you may have learned in Linear Algebra (which should be a prerequisite for both differential equations and partial differential equations) the set of all vectors satisfying Lv= \lambda v for a subspace so there will be an infinite number of eigenfunctions. Generally you will need to write the entire solution as a sum of eigenfunctions.

Once you have all of that theory, you can start talking about Fourier series. Solutions to linear differential equations with constant coefficients can be written as Fourier series
\sum P_n(t)sin(\frac{n\pi}{L}x)+ Q_n(t)cos(\frac{n\pi}{L}x)

And since the derivative of sin is never 0 at multiples of \pi, you could see immediately that the coefficients of sin((n\pi/L)x) must be 0. The simplest way to do this problem, jumping over a lot of preliminary theory, is to start by assuming a function of the form <br /> u(x,t)= \sum_{n=0}^\infty A_n(t)sin(\frac{n\pi}{L}x)