Use the laplace transform to solve the one dimensional boundary value problem
∂u/∂t=∂2u/∂x2, 0<x<1, t>0
subject to u(x,0)=sin πx, u(0,t)=0, u(1,t)=0
Solution;
Take the Laplace transform of both the left hand side as;
"L[\\frac{\\partial u}{\\partial t}]=sU(x,s)-u(x,0)"
And of the right hand side as;
"L[\\frac{\\partial^2u}{\\partial x^2}]=\\frac{\\partial^2U}{{\\partial x^2}}(x,s)"
We have;
"\\frac{\\partial ^2U}{\\partial x^2}(x,s)=sU(x,s)-u(x,0)"
From the problem;
"u(x,o)=sin(\u03c0x)"
Substitute;
"\\frac{\\partial^2U}{\\partial x^2}(x,s)-sU(x,s)=-sin(\u03c0x)"
This is a ODE whose general solution can be written as;
"U(x,s)=U_h(x,s)+U_p(x,s)"
where "U_h(x,s)" is the general solution of the homogeneous problem obtained as;
"m^2-s=0"
"m=\\displaystyle_-^+\\sqrt s"
Hence;
"U_h(x,s)=c_1e^{\\sqrt s x}+c_2e^{-\\sqrt s x}"
and "U_p(x,s)" is any particular solution of the non-homogeneous problem,which is;
"U_p(x,s)=Acos(\u03c0x)+Bsin(\u03c0x)"
We first use the method of undetermined coefficients to find A and B. To this end we have;
"\\frac{d}{dx}U_p(x,s)=-\u03c0Asin(\u03c0x)+\u03c0Bcos(\u03c0x)"
The second derivative;
"\\frac{d^2}{dx^2}U_p(x,s)=-\u03c0^2Acos(\u03c0x)-\u03c0^2Bcos(\u03c0x)"
Therefore;
"\\frac{d^2}{dx^2}U_p(x,s)-U_p(x,s)=-sin(\u03c0x)"
By substitution;
"(-\u03c0^2-s)(Acos(\u03c0x)+Bsin(\u03c0x))=-sin(\u03c0x)"
From which we conclude that;
"-(\u03c0^2+s)A=0" ; A=0
Also;
"-(\u03c0^2+s)B=-1" ;"B=\\frac{1}{s+\u03c0^2}"
Now we have the general solution as;
"U(x,s)=c_1e^{\\sqrt sx}+c_2e^{-\\sqrt sx}+\\frac{1}{s+\u03c0^2}sin(\u03c0x)"
We note the the Laplace transforms of the boundary conditions give;
"u(0,t)=0 \\implies U(0,s)=0"
"u(1,t)=0\\implies U(1,s)=0"
So we have;
"U(0,s)=0=c_1+c_2"
"U(1,s)=0=c_1e^{\\sqrt s}+c_2e^{-\\sqrt s}"
Which gives "c_1=c_2=0" and we have;
"U(x,s)=\\frac{1}{s+\u03c0^2}sin(\u03c0x)"
To find our solution we apply the inverse Laplace transform;
"u(x,t)=sin(\u03c0x)L^{-1}[\\frac{1}{s+\u03c0^2}]"
"u(x,t)=e^{-\u03c0^2t}sin(\u03c0x)"
Comments
Leave a comment