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

Click For Summary
SUMMARY

The discussion focuses on the application of the Fourier Transform in solving a Laplace equation, specifically addressing the inverse transform of the function $$\tilde{u}(k, y) = \frac{1-e^{-ik}}{ik}e^{-ky}$$. Participants clarify that the correct inverse transform involves the convolution of the original function with the inverse transform of $$e^{-|k|y}$$, leading to the expression $$u(x,y) = \int_{-\infty}^\infty u_0(\xi) F(x - \xi,y)\,d\xi$$. Key insights include the necessity for the solution to approach zero as $$y$$ approaches infinity and the requirement for the function to be odd in $$x$$ to satisfy boundary conditions.

PREREQUISITES
  • Fourier Transform theory
  • Laplace equation fundamentals
  • Convolution operations in functional analysis
  • Boundary conditions in partial differential equations
NEXT STEPS
  • Study the properties of the Fourier Transform in solving partial differential equations
  • Learn about convolution integrals and their applications in signal processing
  • Explore boundary value problems and their solutions in mathematical physics
  • Investigate the implications of odd and even functions in Fourier analysis
USEFUL FOR

Mathematicians, physicists, and engineers working on differential equations, particularly those interested in Fourier analysis and boundary value problems.

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 <br /> \hat u(k,y) = \hat u_0(k)e^{-|k|y}. (We need \hat u \to 0 as y \to \infty for both positive and negative k). This is a product of transforms, so when you invert it you obtain a convolution: <br /> u(x,y) = \int_{-\infty}^\infty u_0(\xi) F(x - \xi,y)\,d\xi where F is the inverse transform of e^{-|k|y}. But you also need u(0,y) = 0, which you can ensure by assuming the solution to be odd in x on (-\infty, \infty). So you need to take <br /> u_0(\xi) = -u_0(-\xi) for \xi &lt; 0.
 
Last edited:
pasmith said:
I think you want to say that <br /> \hat u(k,y) = \hat u_0(k)e^{-|k|y}. (We need \hat u \to 0 as y \to \infty for both positive and negative k). This is a product of transforms, so when you invert it you obtain a convolution: <br /> u(x,y) = \int_{-\infty}^\infty u_0(\xi) F(x - \xi,y)\,d\xi where F is the inverse transform of e^{-|k|y}. But you also need u(0,y) = 0, which you can ensure by assuming the solution to be odd in x on (-\infty, \infty). So you need to take <br /> u_0(\xi) = -u_0(-\xi) for \xi &lt; 0.
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 \hat{u}(k,y) \to 0 as y \to \infty for every k \in \mathbb{R}. For k &gt; 0 that means we can only use the e^{-ky} solution, for k &lt; 0 we can only use the e^{ky} solution, and for k = 0 we can only use the constant solution. Hence <br /> \begin{align*}<br /> \hat u(k,y) &amp;= \begin{cases} \hat{u}_0(k)e^{-ky} &amp; k &gt; 0 \\<br /> \hat{u}_0(0) &amp; k = 0 \\<br /> \hat{u}_0(k)e^{ky} &amp; k &lt; 0 \end{cases} \\<br /> &amp;= \hat{u}_0(k)e^{-|k|y}.\end{align*}
 
pasmith said:
I think you want to say that <br /> \hat u(k,y) = \hat u_0(k)e^{-|k|y}. (We need \hat u \to 0 as y \to \infty for both positive and negative k). This is a product of transforms, so when you invert it you obtain a convolution: <br /> u(x,y) = \int_{-\infty}^\infty u_0(\xi) F(x - \xi,y)\,d\xi where F is the inverse transform of e^{-|k|y}. But you also need u(0,y) = 0, which you can ensure by assuming the solution to be odd in x on (-\infty, \infty). So you need to take <br /> u_0(\xi) = -u_0(-\xi) for \xi &lt; 0.
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
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 6 ·
Replies
6
Views
3K
Replies
8
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
5
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 10 ·
Replies
10
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K