We will solve the problem in next full description
1) Differential equation
u t = c 2 u x x u_t=c^2u_{xx} u t = c 2 u xx 0 < x < l . t > 0 0<x<l.\space t>0 0 < x < l . t > 0
2) Boundary conditions (zero conditions for example)
u(0,t)=0, u(l,t)=0;
3) Initial condition
u(x,0)=ϕ ( x ) , 0 < x < l \phi(x),0<x<l ϕ ( x ) , 0 < x < l
Solution:
Let u(x,y)=X(x)T(t) or have separeted variables.
Inserting this form in differential equation we have:
X ( x ) T ′ ( t ) = c 2 T ( t ) X ′ ′ ( x ) X(x)T'(t)=c^2T(t)X''(x) X ( x ) T ′ ( t ) = c 2 T ( t ) X ′′ ( x )
Dividing both paths by c 2 X ( x ) T ( t ) c^2X(x)T(t) c 2 X ( x ) T ( t ) we have
1 c 2 ⋅ T ′ ( t ) T ( t ) = X ′ ′ ( x ) X ( x ) \frac{1}{c^2}\cdot \frac{T'(t)}{T(t)}=\frac{X''(x)}{X(x)} c 2 1 ⋅ T ( t ) T ′ ( t ) = X ( x ) X ′′ ( x )
In this equation left part depend on t bur right part depend on x therefore both part are equal some constant which may be denoted as -λ \lambda λ
Thus we have the system:
x = { T ′ + λ c 2 T ( t ) = 0 X ′ ′ ( x ) + λ ⋅ X ( x ) = 0 x = \begin{cases}
T'+\lambda c^2T(t)=0\\
X''(x)+\lambda\cdot X(x)=0
\end{cases} x = { T ′ + λ c 2 T ( t ) = 0 X ′′ ( x ) + λ ⋅ X ( x ) = 0
Let we consider second eqution
X ′ ′ ( x ) + λ ⋅ X ( x ) = 0 X''(x)+\lambda\cdot X(x)=0 X ′′ ( x ) + λ ⋅ X ( x ) = 0
Boundary conditions at x=0 and x=l take the form:
X(0)=X(l)=0.
1) λ < 0 \lambda<0 λ < 0
General solution of ode is X ( x ) = C 1 e − − λ x + C 2 e λ x X(x)=C_1e^{-\sqrt{-\lambda}x}+C_2e^{\sqrt{\lambda }x} X ( x ) = C 1 e − − λ x + C 2 e λ x
Then X ( 0 ) = C 1 + C 2 = 0 X(0)=C_1+C_2=0 X ( 0 ) = C 1 + C 2 = 0 ,C 2 = − C 1 C_2=-C_1 C 2 = − C 1
Therefore X ( x ) = C ⋅ ( e − − λ x − e − λ x ) X(x)=C\cdot (e^{-\sqrt{-\lambda}x}-e^{\sqrt{-\lambda}x}) X ( x ) = C ⋅ ( e − − λ x − e − λ x )
Condition X(l)=0 implies
C ⋅ ( e − − λ l − e − λ l ) = 0 C\cdot (e^{-\sqrt{-\lambda}l}-e^{\sqrt{-\lambda}l})=0 C ⋅ ( e − − λ l − e − λ l ) = 0 and hence C=0, therefore X(x)≡ 0 \equiv0 ≡ 0
but we need nontrival solution hence case λ < 0 \lambda<0 λ < 0 is impossible, case λ = 0 \lambda=0 λ = 0 is impossible too with similar argues, typical case is λ > 0 \lambda>0 λ > 0
In this case X ( x ) = C 1 s i n ( λ x ) + C 2 c o s ( λ x ) X(x)=C_1sin(\sqrt \lambda x)+C_2cos(\sqrt \lambda x) X ( x ) = C 1 s in ( λ x ) + C 2 cos ( λ x )
X(0)= C 2 = 0 = > C_2=0=> C 2 = 0 => X(x)=C s i n ( λ x ) Csin(\sqrt \lambda x) C s in ( λ x )
X(l)=0=>s i n ( λ l ) = 0 = > λ l = π ⋅ n . n > 0 sin(\sqrt\lambda l)=0=>\sqrt \lambda l=\pi\cdot n.n>0 s in ( λ l ) = 0 => λ l = π ⋅ n . n > 0
λ n = ( π ⋅ n l ) 2 \lambda_n=\left(\frac{\pi\cdot n}{l} \right)^2 λ n = ( l π ⋅ n ) 2 and X n ( x ) = s i n ( π ⋅ n ⋅ x l ) X_n(x)=sin\left(\frac{\pi\cdot n\cdot x}{l}\right) X n ( x ) = s in ( l π ⋅ n ⋅ x )
For this λ n \lambda_n λ n an equation T ′ + λ c 2 T ( t ) = 0 T'+\lambda c^2T(t)=0 T ′ + λ c 2 T ( t ) = 0 has solution
T n ( t ) = e − λ n ⋅ c 2 ⋅ t T_n(t)=e^{-\lambda_n\cdot c^2\cdot t} T n ( t ) = e − λ n ⋅ c 2 ⋅ t and u + n ( x , t ) = T n ( t ) X n ( x ) = e − λ n c 2 t s i n ( π n x l ) u+n(x,t)=T_n(t)X_n(x)=e^{-\lambda_nc^2t}sin\left( \frac{\pi nx}{l}\right) u + n ( x , t ) = T n ( t ) X n ( x ) = e − λ n c 2 t s in ( l πn x )
After that we search solution in the form
∑ n = 1 ∞ C n e − λ n c 2 t s i n ( π n x l ) \sum_{n=1}^{\infty}C_ne^{-\lambda_nc^2t}sin\left( \frac{\pi nx}{l}\right) ∑ n = 1 ∞ C n e − λ n c 2 t s in ( l πn x ) for initial condition to satisfy.
Insertinf t=0 we have
∑ n = 1 ∞ C n s i n ( π n x l ) = ϕ ( x ) , 0 < x < l \sum_{n=1}^{\infty}C_nsin\left( \frac{\pi nx}{l}\right)=\phi(x),0<x<l ∑ n = 1 ∞ C n s in ( l πn x ) = ϕ ( x ) , 0 < x < l
Then C n = 2 l ∫ 0 l ϕ ( x ) s i n π n z l d z C_n=\frac{2}{l}\int_0^l\phi(x)sin\frac{\pi n z}{l}dz C n = l 2 ∫ 0 l ϕ ( x ) s in l πn z d z -Fourier coefficients. With such coefficients solution of the problem is:
∑ n = 1 ∞ C n e − λ n c 2 t s i n ( π n x l ) \sum_{n=1}^{\infty}C_ne^{-\lambda_nc^2t}sin\left( \frac{\pi nx}{l}\right) ∑ n = 1 ∞ C n e − λ n c 2 t s in ( l πn x )
Comments