2 Homogeneous Heat PDE in 1D

 2.1 Boundary conditions both Dirichlet and Homogeneous
 2.2 Boundary conditions both Dirichlet, one homogeneous other end nonHomogeneous
 2.3 Boundary conditions both Neumann and homogeneous
 2.4 Boundary conditions one end Neumann, the other nonHomogeneous Dirichlet
 2.5 Boundary conditions at both ends homogeneous Neumann (insulated)
 2.6 Boundary conditions at one end non-homogeneous Neumannm other end homogeneous Dirichlet

2.1 Boundary conditions both Dirichlet and Homogeneous

2.1.1 Example 1

Solve

1-∂u-  ∂2u-
k ∂t = ∂x2    0 < x < L,t > 0

Boundary conditions

pict

Initial conditions

         {                 L-
u (x,0) =     x      0L < x < 2
           L − x    2 < x < L

Solution.      nπ
λn = L

             ∞  (           )
u (x,t) = 4L ∑     1-sin(nπ-)  exp(− kλ2t)sin(λ x)
        π2 n=odd  n2     2           n       n

  

2.1.2 Example 2

A bar length L = 10  cm with insulated sides is initially at 100  degrees. starting at t = 0  , the ends are held at zero degree. Find the temperature distribution in the bar at time t

1 ∂u    ∂2u
-2---=  --2-   0 < x < L,t > 0
α  ∂t   ∂x

Boundary conditions

pict

Initial conditions

u(x,0) = T0 = 1000

Solution

Let

u = X (x)T (t)

Substituting this into the PDE gives

1
-2T′X = X ′′T
α

Dividing by XT

1-T-′= X-′′
α2 T    X

Therefore

 1 T′   X′′
α2-T-=  X--= − λ

Where λ  is some constant. Hence we obtain two ODE’s to solve. They are

pict

And

T′ + λα2T = 0

case λ = 0

pict

From first boundary conditions, X (0) = 0  , the above becomes

0 = B

Hence X  = Ax  . From second boundary conditions 0 = AL  , or A = 0  , therefore trivial solution and λ = 0  is not an eigenvalue.

case λ > 0

  ′′
X  + λX  = 0

The solution is

         (  --)       ( -- )
X = A cos √ λx + B sin  √λx

From first B.C.  0 = A  and the solution becomes          (    )
X = B sin  √λx . From second B.C.

        (    )
0 = B sin √λL

For non-trivial solution, we want    (√ --)
sin    λL  = 0  or

pict

Therefore

             ( nπ )
Xn (x) = bn sin --x
               L

Now we solve the time ODE  T ′ +λn α2T = 0  . This has solution

Tn = ane−λnα2t

Therefore the solution becomes

pict

To find Bn  , we use orthoginality. At t = 0  , we are given u (x,0)  . Hence at t = 0  the solution becomes

         ∞       (    )
u (x,0) = ∑  Bn sin  nπx
         n=1       L

Multiplying both sides by    (    )
sin mπL-x and integrating

pict

Therefore

pict

Therefore, the solution is

pict

Here is an animation.

  

2.2 Boundary conditions both Dirichlet, one homogeneous other end nonHomogeneous

A bar length L

1 ∂u    ∂2u
α2-∂t = ∂x2-   0 < x < L,t > 0

Boundary conditions

pict

Initial conditions

u (x,0) = x

Solution,      nπ
λn = L-

                       ∞∑
u(x,t) = T0 − T0 x− 4T0    1-exp(− α2λ2nt) sin(λnx)
             L     π  n=evenn

  

2.3 Boundary conditions both Neumann and homogeneous

A bar length L

1 ∂u    ∂2u
α2-∂t = ∂x2-   0 < x < L,t > 0

Boundary conditions, insulated

pict

Initial conditions

u (x,0) = x

Solution, λn = nπ
     L

        L-  4L- ∑∞  -1    (  2 2 )
u(x,t) = 2 − π2     n2 exp − α λnt cos(λnx)
               n=odd

  

2.4 Boundary conditions one end Neumann, the other nonHomogeneous Dirichlet

A bar length L = 2

1 ∂u    ∂2u
α2-∂t = ∂x2-   0 < x < L,t > 0

Boundary conditions, left end only insulated

pict

Initial conditions

u(x,0) = 00

Solution,      2n+1π
λn = --2-L-

                ∞∑      n
u(x,t) = T0 − 4T0   (− 1)-exp(− α2λ2nt) cos (λnx)
             π  n=12n +1

  

2.5 Boundary conditions at both ends homogeneous Neumann (insulated)

See my solution for exam 1, Math 322. UW, Fall 2016, problem 3.

Solve, 0 < x < π

Solve

        2
∂u-= α ∂-u− βu    0 < x < L, t > 0
∂t     ∂x2

over 0 < x < π  with, a > 0,β > 0

pict

And initial conditions

u (x,0) = x

Solution

         π    (       )   2∑∞  (− 1n − 1)         (2   )
u (x,t) = 2 + c0 e−βt − 1 + π   ---n2----cos(nx)e− n α+β t
                           n=1

Solution is not unique, since there is unknown c0  . To find c0  , the solvability condition for   2
∇ u = 0  with Neumann B.C. is used, which says that total flux must be zero at steady state, which is the case here, since flux is zero as given (insulated). Since solvability condition is satisfied, then energy is conserved. Solution at t = ∞ from above is

         1
u(x,∞ ) = 2 π− c0

Since energy is conserved then comparing it to energy at initial conditions

∫ π             ∫ π
   ρcu(x,0)dx =    ρcu(x,∞ )dx
 0               0

But u(x,0) = x  , hence

pict

Hence the final solution is

              ∑∞    n              (    )
u(x,t) = π-+ 2   (− 1-−-1)cos(nx)e− n2α+ βt
        2   π n=1   n2

Using α = 2,β = 1  , here is animation for 3.5 seconds

  

2.6 Boundary conditions at one end non-homogeneous Neumannm other end homogeneous Dirichlet

This 1D heat PDE has one end with boundary condition that is time dependent.

pict

Now let

pict

Be some reference temperature distribution that only needs to satisfy the boundary conditions given. i.e. ∂∂rx (0,t) = et,r(L,t) = 0  . Then the difference temperature distribution is

v (x,t) = u (x,t)− r(x,t)

Since u,r  satisfy the nonhomogeneous B.C’s, then v  satisfies the homogeneous B.C which is ∂∂vx (0,t) = 0,v(L,t) = 0  . With the initial conditions

pict

Since v(x,t)  has homogeneous B.C. it can be easily solved. The solution can be found by separation of variables to be

          ∑∞         (∘ ---) − kλnt
v(x,t) =        An cos   λnx  e
        n=1,3,5,⋅⋅⋅

Where eigenvalue λn = (nπ)2
      2L  . An  is now found from initial conditions.

           ∞∑         ( nπ )
v (x,0) =        An cos  2Lx
        n=1,3,5,⋅⋅⋅

Multiplying both sides by    (m-π )
cos  2L x ,  integrating and changing the order of integration and summation on RHS gives

∫ L           (m π )        ∞∑       ∫ L   (m π )    (nπ  )
   (− x + L)cos---x  dx =        An    cos ---x  cos ---x dx
 0             2L        n=1,3,5,⋅⋅⋅    0     2L        2L

By orthogonality of cos  the above reduces to one term

pict

Hence the solution is

          ∑∞         (    )    (nπ)2
v (x,t) =        Ancos  nπx  e−k 2L- t
        n=1,3,5,⋅⋅⋅       2L

But u = v+ r  , therefore the final solution is

                     ∞∑         (    )    (  )
u(x,t) = et(x− L) +       An cos  nπx  e−k n2Lπ2t
                  n=1,3,5,⋅⋅⋅       2L

Here is animation using Mathematica, for L = 1,k = 0.01  for 2 seconds.