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

Click For Summary
SUMMARY

The discussion focuses on solving the second-order partial differential equation (PDE) given by the heat equation, \( u_t = u_{xx} \), with boundary condition \( u(0,t) = \sin t \) and initial condition \( u(x,0) = 0 \). The Laplace transform is applied, leading to the transformed equation \( su(x,s) - u(x,0) = u_{xx}(x,s) \) and the boundary condition \( u(0,s) = \frac{1}{1+s^2} \). The solution is derived as \( u(x,s) = \frac{1}{1+s^2} e^{-x\sqrt{s}} \), followed by the application of the inverse Laplace transform, which involves convolution, resulting in the final solution \( u(x,t) = (\sin(t) * \frac{x e^{-\frac{x^2}{4t}}}{2\sqrt{\pi} t^{3/2}}) \).

PREREQUISITES
  • Understanding of Laplace transforms and their properties
  • Knowledge of solving second-order partial differential equations
  • Familiarity with convolution theorem in the context of Laplace transforms
  • Basic concepts of boundary and initial value problems
NEXT STEPS
  • Study the properties of Laplace transforms, particularly the inverse transform
  • Learn about convolution in the context of Laplace transforms
  • Explore methods for solving boundary value problems for PDEs
  • Investigate the application of Fourier transforms in solving similar PDEs
USEFUL FOR

Mathematicians, physicists, and engineers dealing with heat equations, boundary value problems, and those looking to deepen their understanding of Laplace transforms and their applications 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 3 ·
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
10K
  • · Replies 7 ·
Replies
7
Views
6K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
10K
  • · Replies 1 ·
Replies
1
Views
10K