Confirming Green's function for homogeneous Helmholtz equation (3D)

Click For Summary
SUMMARY

The discussion focuses on confirming Green's function for the homogeneous Helmholtz equation in three dimensions. The equation analyzed is ##\nabla^2 G + k^2 G = 0##, where the proposed solution leads to confusion regarding the delta function behavior at r=0. The participants clarify that the calculations are valid only for r>0 due to the singularity at r=0 and suggest integrating over a small sphere to properly evaluate the delta function. Additionally, they reference the special case of Poisson’s equation, indicating that it has a solution proportional to 1/r, which exhibits a delta inhomogeneity at the origin.

PREREQUISITES
  • Understanding of the Helmholtz equation in three dimensions
  • Familiarity with Green's functions and their applications
  • Knowledge of delta functions and their properties
  • Basic concepts of singularities in coordinate systems
NEXT STEPS
  • Study the derivation and properties of Green's functions for the Helmholtz equation
  • Learn about the integration techniques for singularities in differential equations
  • Explore the implications of Poisson’s equation and its solutions
  • Investigate the behavior of solutions near singular points in various coordinate systems
USEFUL FOR

Mathematicians, physicists, and engineers working with differential equations, particularly those focused on wave propagation and potential theory.

PhDeezNutz
Messages
851
Reaction score
561
Homework Statement
The 3D Helmholtz equation is

##\left(\nabla^2 + k^2 \right) \Psi \left( r \right)= 0##

Supposedly the Green's function for this equation is

##G\left(r \right) = - \frac{1}{4 \pi} \frac{e^{ikr}}{r}##
Relevant Equations
A green's function is defined as the solution to the following

##\left( \nabla^2 + k^2 \right) G = \delta \left( r \right)##

The laplacian in spherical coordinates (for purely radial dependence) is

##\nabla^2 t = \frac{1}{r^2} \frac{\partial }{\partial r} \left(\frac{1}{r^2}\frac{\partial t}{\partial r} \right)##
Plugging in the supposed ##G## into the delta function equation

##\nabla^2 G = -\frac{1}{4 \pi} \frac{1}{r^2} \frac{\partial}{\partial r} \left(\frac{r^2 \left(ikr e^{ikr} - e^{ikr} \right)}{r^2} \right)##

##= -\frac{1}{4 \pi} \frac{1}{r^2} \left[ike^{ikr} - rk^2 e^{ikr} - ike^{ikr} \right]##

##= \frac{k^2 e^{ikr}}{4 \pi r}##

##k^2 G## is simply

##k^2 G = - \frac{k^2 e^{ikr}}{4 \pi r}##

So we get

##\left( \nabla^2 + k^2 \right) G = 0##

I know the delta function is zero everywhere else besides r = 0 where it is infinity, but I'm getting 0 across the board instead of a delta function.

Thanks for any help in advanced.
 
Physics news on Phys.org
PhDeezNutz said:
but I'm getting 0 across the board instead of a delta function
No you are not. What you are doing is only valid for r>0 as your coordinate system is singular at r=0.

Try integrating the differential equation over a small sphere.

Edit: Also consider that the special case of Poisson’s equation (ie, Helmholtz with k=0) has a solution proportional to 1/r. This is known to have a delta inhomogeneity at r=0.
 
  • Like
Likes   Reactions: troglodyte and PhDeezNutz

Similar threads

Replies
2
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
11
Views
1K
  • · Replies 2 ·
Replies
2
Views
3K
Replies
4
Views
3K
Replies
3
Views
1K
  • · Replies 17 ·
Replies
17
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 7 ·
Replies
7
Views
806