One Dimensional Diffusion Equation

Click For Summary
SUMMARY

The discussion focuses on solving the one-dimensional diffusion equation given by \(\frac{\partial u}{\partial t} = k\frac{\partial^2 u}{\partial x^2}\) with specific initial and non-homogeneous boundary conditions. The initial condition is defined as \(u(x,0)=f(x)\) where \(f(x)\) is piecewise defined. The solution approach involves using Fourier series and separation of variables, but challenges arise in satisfying both the boundary and initial conditions. A suggested method to address the non-homogeneous boundary conditions is to "homogenize" them by introducing a new function \(v(x,t)\) that transforms the problem into one with homogeneous boundary conditions.

PREREQUISITES
  • Understanding of partial differential equations (PDEs)
  • Familiarity with Fourier series and their applications
  • Knowledge of boundary value problems
  • Concept of homogenization in mathematical analysis
NEXT STEPS
  • Study the method of separation of variables for PDEs
  • Learn how to express functions as Fourier series
  • Explore techniques for solving non-homogeneous boundary value problems
  • Investigate the concept of homogenization in differential equations
USEFUL FOR

Mathematics students, educators, and researchers focusing on partial differential equations, particularly those dealing with diffusion processes and boundary value problems.

NeoDevin
Messages
334
Reaction score
2

Homework Statement


Solve:

\frac{\partial u}{\partial t} = k\frac{\partial^2 u}{\partial x^2}, 0<x<\pi, t>0

with initial condition

u(x,0)=f(x)=\left\{\begin{array}{cc} 1,& 0\leq x< \pi/2 \\ 0, &\pi/2 \leq x < \pi \end{array}\right

and with non-homogeneous boundary conditions

u(0,t) = 1, u(\pi,t)=0

The Attempt at a Solution



I've tried this one a couple different ways, I tried separation of variables, and Fourier series. I can get a solution (or infinitely many) to the equation, but I can't seem to make them fit the boundary/initial conditions both.

My solution for the Fourier series method is:

u(x,0)=\frac{A_0}{2} + \sum^{\infty}_{n=1} e^{-kn^2t}(A_n \cos(nx) + B_n \sin(nx)

Any ideas how to make this (or any other solution) match both the boundary and initial conditions? Or is the problem incosistent?
 
Physics news on Phys.org
it would be easier to see this if you express the function f(x), the initial condition, as a Fourier series.
 
There is a standard way of handling problems with non-homogeneous boundary conditions: "homogenize" them!

The simple linear function 1- \frac{1}{\pi}x satisfies the boundary conditions. Now let v(x,t)= u(x,t)- 1+ \frac{1}{\pi}x
It is easy to see that v(x,t) also satisfies
\frac{\partial v}{\partial t} = k\frac{\partial^2 v}{\partial x^2}
but the boundary conditions for v are v(0, t)= 0, v(\pi,t)= 0.
Of course, you have to change the initial condition to
v(x,0)=\left\{\begin{array}{cc} \frac{1}{\pi}x,& 0\leq x< \pi/2 \\ \frac{1}{\pi}x-1, &\pi/2 \leq x < \pi \end{array}\right

Now you can write v in a purely sine series:
v(x,t)=\sum^{\infty}_{n=1} e^{-kn^2t}B_n \sin(nx)

Of course, as mjsd said, you will eventually have to write the initial condition as a Fourier series in order to find the Bn
 

Similar threads

  • · Replies 11 ·
Replies
11
Views
3K
Replies
0
Views
1K
  • · Replies 5 ·
Replies
5
Views
1K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 105 ·
4
Replies
105
Views
11K
  • · Replies 14 ·
Replies
14
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K
Replies
6
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K