- #1

- 5

- 0

[tex]u_t=u_{xx}+2u_x[/tex]

0<=x<=L, t>=0, u(x,0)=f(x), u_x(0,t)=u_x(L,t)=0

How to do this?

0<=x<=L, t>=0, u(x,0)=f(x), u_x(0,t)=u_x(L,t)=0

How to do this?

You are using an out of date browser. It may not display this or other websites correctly.

You should upgrade or use an alternative browser.

You should upgrade or use an alternative browser.

- Thread starter useruseruser
- Start date

- #1

- 5

- 0

[tex]u_t=u_{xx}+2u_x[/tex]

0<=x<=L, t>=0, u(x,0)=f(x), u_x(0,t)=u_x(L,t)=0

How to do this?

0<=x<=L, t>=0, u(x,0)=f(x), u_x(0,t)=u_x(L,t)=0

How to do this?

- #2

HallsofIvy

Science Advisor

Homework Helper

- 41,847

- 966

Assume u is of the form u(x,t)= A(x)B(t). Then u

Divide both sides by AB:

[tex]\frac{B'}{B}= \frac{A"+2A'}{A}[/tex]

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 "[itex]\lambda[/itex]".

Then we have [itex]B'/B= \lambda[/itex] or [itex]B'= \lambda B[/itex] and [itex](A"+ 2A')/A= \lambda[/itex] or [itex]A"+ 2A'= \lambda A[/itex].

Now look at the boundary conditions: [itex]u_x(0,t)= A(0)B'(t)= 0[/itex] for all t so we must have A(0)= 0. Also, [itex]u_x(L,t)= A(0)B'(t)= 0[/itex] for all t so A(L)= 0.

Solve the boundary problem [itex]A"+ 2A'- \lamba A= 0[/itex], A(0)= 0, A'(0)= 0 to determine what [itex]\lambda[/itex] can be. You should see that this is an

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 [itex]Lv= \lambda v[/itex] 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

[tex]\sum P_n(t)sin(\frac{n\pi}{L}x)+ Q_n(t)cos(\frac{n\pi}{L}x)[/tex]

And since the derivative of sin is never 0 at multiples of [itex]\pi[/itex], you could see immediately that the coefficients of [itex]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

[tex]u(x,t)= \sum_{n=0}^\infty A_n(t)sin(\frac{n\pi}{L}x)[/tex]

Share: