MHB Solve Laplace Transform ODE: $\sin t, 0<t<\infty$

Click For Summary
The discussion focuses on solving the heat equation using the Laplace transform with boundary conditions involving a sine function. The initial Laplace transform leads to a second-order ordinary differential equation, and the solution is expressed in terms of exponential functions. The key challenge is applying the inverse Laplace transform correctly, particularly using the convolution theorem for products of transforms. The final solution is presented as a convolution of the sine function with a specific exponential term. The participants confirm the correctness of the steps taken, emphasizing the importance of proper application of transforms in solving differential equations.
Markov2
Messages
149
Reaction score
0
Solve

$\begin{eqnarray*}
{{u}_{t}}&=&{{u}_{xx}},\text{ }0<x<\infty ,\text{ }0<t<\infty \\
u(0,t)&=&\sin t,\text{ }0<t<\infty \\
u(x,0)&=&0,\text{ }0\le x<\infty .
\end{eqnarray*}$

I need to apply the Laplace transform to solve this, so by applying it I get $su(x,s)-u(x,0)=u_{xx}(x,s)$ and $u(0,s)=\frac1{1+s^2}.$ I need to treat the first equation as a second order ODE right? Then the solution would be $u(x,s)=c_1e^{x\sqrt s}+c_2e^{-x\sqrt s},$ here's where I'm getting stuck, I've been told that I must found bounded solutions, so in order to get it bounded I should consider $u(x,s)=c_2e^{-x\sqrt s},$ by applying the initial conditions I get $u(0,s)=c_2=\frac1{1+s^2}$ so we have $u(x,s)=\frac1{1+s^2}e^{-x\sqrt s},$ now I must apply the inverse Laplace transform, so first $u(x,s)=\frac s{1+s^2}\frac{e^{-x\sqrt s}}s=\mathcal L(\cos t)\mathcal L^{-1}\left(\frac{e^{-x\sqrt s}}s\right)$
I can't find the inverse of the remaining function.

Is this correct?
Thanks!
 
Physics news on Phys.org
Markov said:
Solve

$\begin{eqnarray*}
{{u}_{t}}&=&{{u}_{xx}},\text{ }0<x<\infty ,\text{ }0<t<\infty \\
u(0,t)&=&\sin t,\text{ }0<t<\infty \\
u(x,0)&=&0,\text{ }0\le x<\infty .
\end{eqnarray*}$

I need to apply the Laplace transform to solve this, so by applying it I get $su(x,s)-u(x,0)=u_{xx}(x,s)$ and $u(0,s)=\frac1{1+s^2}.$ I need to treat the first equation as a second order ODE right? Then the solution would be $u(x,s)=c_1e^{x\sqrt s}+c_2e^{-x\sqrt s},$ here's where I'm getting stuck, I've been told that I must found bounded solutions, so in order to get it bounded I should consider $u(x,s)=c_2e^{-x\sqrt s},$ by applying the initial conditions I get $u(0,s)=c_2=\frac1{1+s^2}$ so we have $u(x,s)=\frac1{1+s^2}e^{-x\sqrt s},$ now I must apply the inverse Laplace transform, so first $u(x,s)=\frac s{1+s^2}\frac{e^{-x\sqrt s}}s=\mathcal L(\cos t)\mathcal L^{-1}\left(\frac{e^{-x\sqrt s}}s\right)$

There are two problems with this approach.

1. Why insert the additional s's? You already know the inverse LT of $1/(s^{2}+1)$.
2. Second of all, the inverse LT of a product is NOT the product of the inverse LT's. Instead, you must use the convolution theorem, which says that

$$\mathcal{L}^{-1}\{F(s)\cdot G(s)\}=\int_{0}^{t}f(\tau)g(t-\tau)\,d\tau\equiv(f* g)(t).$$

As it turns out, the inverse LT of the remaining piece is

$$\mathcal{L}^{-1}\{e^{-x\sqrt{s}}\}=\frac{xe^{-\frac{x^{2}}{4t}}}{2\sqrt{\pi}\,t^{3/2}}.$$

Can you finish from here?
 
Oh yes, so it's actually $\mathcal L(\sin t)\mathcal L\left( {\frac{{x{e^{ - \frac{{{x^2}}}{{4t}}}}}}{{2\sqrt \pi {t^{\frac{3}{2}}}}}} \right),$ so the solution is $u(x,t)=(\sin (t) *{\frac{{x{e^{ - \frac{{{x^2}}}{{4t}}}}}}{{2\sqrt \pi {t^{\frac{3}{2}}}}}}),$ is it correct now?
 
Markov said:
Oh yes, so it's actually $\mathcal L(\sin t)\mathcal L\left( {\frac{{x{e^{ - \frac{{{x^2}}}{{4t}}}}}}{{2\sqrt \pi {t^{\frac{3}{2}}}}}} \right),$ so the solution is $u(x,t)=(\sin (t) *{\frac{{x{e^{ - \frac{{{x^2}}}{{4t}}}}}}{{2\sqrt \pi {t^{\frac{3}{2}}}}}}),$ is it correct now?

If by
$$\sin (t) *{\frac{{x{e^{ - \frac{{{x^2}}}{{4t}}}}}}{{2\sqrt \pi {t^{\frac{3}{2}}}}}}$$
you mean $\sin(t)$ convolved with
$${\frac{{x{e^{ - \frac{{{x^2}}}{{4t}}}}}}{{2\sqrt \pi {t^{\frac{3}{2}}}}}},$$
I would agree. You might want to write out what that is in terms of the integral. If the integral is tractable, you might even compute it.
 
Yes, I meant the convolution, thanks!
 
You're welcome. Have a good one!
 
Markov said:
I get $u(0,s)=c_2=\frac1{1+s^2}$ so we have $u(x,s)=\frac1{1+s^2}e^{-x\sqrt s},$
Wait a sec, is this step correct? I remember I applied the Laplace transform to $u(0,t),$ and substitute in the solution which I applied the Laplace transform? Is it correct? I'm not sure, because when dealing with Fourier transform, in this post: http://www.mathhelpboards.com/showthread.php?85-Solving-PDF-by-using-Fourier-transform&p=446&viewfull=1#post446

It wasn't correct, can you verify?
 
Markov said:
Wait a sec, is this step correct?

It's correct. You took the LT correctly, applied the initial conditions correctly.
 
I did the same on http://www.mathhelpboards.com/showthread.php?85-Solving-PDF-by-using-Fourier-transform&p=446&viewfull=1#post446

So why it wasn't correct there?
 
  • #10
Markov said:
I did the same on http://www.mathhelpboards.com/showthread.php?85-Solving-PDF-by-using-Fourier-transform&p=446&viewfull=1#post446

So why it wasn't correct there?

Because you didn't do the same thing in the other thread. You didn't take the FT of the initial condition, you just took the initial condition straight up and applied that. ThePerfectHacker pointed out that you needed to take the FT of the initial condition.
 

Similar threads

  • · Replies 17 ·
Replies
17
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
6K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 5 ·
Replies
5
Views
5K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 1 ·
Replies
1
Views
2K