Is the Fourier Transform Correctly Applied in Solving This Laplace Equation?

lriuui0x0
Messages
101
Reaction score
25
Homework Statement
Solve the laplace equation ##u_{xx} + u_{yy} = 0##, with ##x \ge 0, y \ge 0##, using Fourier transform.

Subject to the boundary conditions:

$$
\begin{aligned}
u(x, 0) &= \begin{cases}1 & 0 < x < 1 \\ 0 & \text{otherwise}\end{cases} \\
u(0, y) &= \lim_{x\to\infty} u(x,y) = \lim_{y\to\infty} u(x,y) = 0
\end{aligned}
$$
Relevant Equations
Fourier transform and inverse transform
I have tried to Fourier transform in ##x## and get the result in the transformed coordinates, please check my result:

$$
\tilde{u}(k, y) = \frac{1-e^{-ik}}{ik}e^{-ky}
$$

However, I'm having some problems with the inverse transform:

$$
\frac{1}{2\pi}\int_{-\infty}^\infty \frac{1-e^{-ik}}{ik}e^{-yk}dk
$$

Not sure how to do this integral. The solution says it's

$$
\frac{4xy}{\pi}\int_0^1 \frac{vdv}{[(x-v)^2+y^2][(x+v)^2+y^2]}
$$
 
  • Like
Likes   Reactions: Delta2
Physics news on Phys.org
Sorry the inverse FT integral should be:

$$
\frac{1}{2\pi}\int_{-\infty}^\infty \frac{1-e^{-ik}}{ik}e^{-yk} e^{ixk}dk
$$
 
I think you want to say that [tex] \hat u(k,y) = \hat u_0(k)e^{-|k|y}.[/tex] (We need [itex]\hat u \to 0[/itex] as [itex]y \to \infty[/itex] for both positive and negative [itex]k[/itex]). This is a product of transforms, so when you invert it you obtain a convolution: [tex] u(x,y) = \int_{-\infty}^\infty u_0(\xi) F(x - \xi,y)\,d\xi[/tex] where [itex]F[/itex] is the inverse transform of [itex]e^{-|k|y}[/itex]. But you also need [itex]u(0,y) = 0[/itex], which you can ensure by assuming the solution to be odd in [itex]x[/itex] on [itex](-\infty, \infty)[/itex]. So you need to take [itex] u_0(\xi) = -u_0(-\xi)[/itex] for [itex]\xi < 0[/itex].
 
Last edited:
pasmith said:
I think you want to say that [tex] \hat u(k,y) = \hat u_0(k)e^{-|k|y}.[/tex] (We need [itex]\hat u \to 0[/itex] as [itex]y \to \infty[/itex] for both positive and negative [itex]k[/itex]). This is a product of transforms, so when you invert it you obtain a convolution: [tex] u(x,y) = \int_{-\infty}^\infty u_0(\xi) F(x - \xi,y)\,d\xi[/tex] where [itex]F[/itex] is the inverse transform of [itex]e^{-|k|y}[/itex]. But you also need [itex]u(0,y) = 0[/itex], which you can ensure by assuming the solution to be odd in [itex]x[/itex] on [itex](-\infty, \infty)[/itex]. So you need to take [itex] u_0(\xi) = -u_0(-\xi)[/itex] for [itex]\xi < 0[/itex].
How do we get the ##e^{-|k|y}## term? The way I derived my result is:

$$
\begin{aligned}
\mathcal{F}(\frac{\partial^2u}{\partial x^2} + \frac{\partial^2u}{\partial y^2}) &= 0 \\
\frac{\partial^2\tilde{u}}{\partial y^2} -k^2 \tilde{u} &= 0 \\
\tilde{u}(k, y) &= A(k)e^{ky} + B(k)e^{-ky} \\
\end{aligned}
$$

And I was thinking ##e^{ky}## term must vanish since we require ##\lim_{y\to\infty} \tilde{u}(k,y) = 0##. However thinking about it again, this isn't necessarily true, as we can allow ##k \to -\infty##. But how did you derive the ##e^{-|k|y}## term?
 
Last edited:
lriuui0x0 said:
But how do you deal derive the ##e^{-|k|y}## term?

We need [itex]\hat{u}(k,y) \to 0[/itex] as [itex]y \to \infty[/itex] for every [itex]k \in \mathbb{R}[/itex]. For [itex]k > 0[/itex] that means we can only use the [itex]e^{-ky}[/itex] solution, for [itex]k < 0[/itex] we can only use the [itex]e^{ky}[/itex] solution, and for [itex]k = 0[/itex] we can only use the constant solution. Hence [tex] \begin{align*}<br /> \hat u(k,y) &= \begin{cases} \hat{u}_0(k)e^{-ky} & k > 0 \\<br /> \hat{u}_0(0) & k = 0 \\<br /> \hat{u}_0(k)e^{ky} & k < 0 \end{cases} \\<br /> &= \hat{u}_0(k)e^{-|k|y}.\end{align*}[/tex]
 
pasmith said:
I think you want to say that [tex] \hat u(k,y) = \hat u_0(k)e^{-|k|y}.[/tex] (We need [itex]\hat u \to 0[/itex] as [itex]y \to \infty[/itex] for both positive and negative [itex]k[/itex]). This is a product of transforms, so when you invert it you obtain a convolution: [tex] u(x,y) = \int_{-\infty}^\infty u_0(\xi) F(x - \xi,y)\,d\xi[/tex] where [itex]F[/itex] is the inverse transform of [itex]e^{-|k|y}[/itex]. But you also need [itex]u(0,y) = 0[/itex], which you can ensure by assuming the solution to be odd in [itex]x[/itex] on [itex](-\infty, \infty)[/itex]. So you need to take [itex] u_0(\xi) = -u_0(-\xi)[/itex] for [itex]\xi < 0[/itex].
How is convolution operation defined for functions only on ##\mathbf{R}^+##? Why do I need to extend the domain of ##u## to include the whole real line?
 

Similar threads

Replies
1
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 6 ·
Replies
6
Views
4K
Replies
8
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
5
Views
2K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 10 ·
Replies
10
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K