Spectral Analysis of the Derivation of Green’s Function of Helmholtz Integral Equation via Dirac-Delta Function and Cauchy Residual Approach

Md. Shuzon Ali¹, Sherajum Monira Bristy¹, Md. Abdullah Al Asad^2*, Nur Hasan Mahmud Shahen³

¹Department of Mathematics, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj, Bangladesh.
²Department of Electrical and Electronic Enginering, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj, Bangladesh.
³Department of Arts and Sciences, Bangladesh Army University of Science and Technology, Saidpur, Bangladesh.

DOI: 10.4236/oalib.1110245                

Abstract

In this study, we have determined Green’s functions for Helmholtz integral equations in a spherical polar coordinate system in the whole plane domain with the aid of spectral Fourier transform technique. Our intended Green’s function solution has a dominant role to represent wave propagation with a high quantum wave number. The Diracdelta function also plays an important role here to represent the scattering region for wave propagation. The evaluation of the improper double integrals in the complex plane furnishes our desired Green’s functions. The applied technique allows us to obtain all of the possible Green’s functions by using Somerfield radiation condition. With the help of computational software package MATLAB, we have drawn the solution plot that can express the analogue of wave propagation features. By using the MATLAB software package, we have drawn the solution plot that can express the analogue of wave propagation features.

Keywords

Helmholtz Integral Equation, Summerfield Radiation Condition, Cauchy Residual Approach, Classical Schrödinger Wave Equation, Green’s Function

1. Introduction

In the present world, an integral equation is observed to be a well-built branch of science and technology providing outcomes which advantageously employed by experts working in various fields such as wave acoustics, electromagnetics, plasma physics, and fluid dynamics and so on. One of the most important integral equations of special status is Helmholtz integral equation which displays such kind of wave phenomena of the fundamental particle with Somerfield radiation condition. Besides, scattering of buildings always governs contrary intention in the domain of noise fences.

The Green’s function solution is one of the most significant analytical tools for solving partial differential equations in the field of electromagnetism, wave mechanics, optics, quantum mechanics, fluid dynamics, relativistic particle dynamics, general relativity and so on [1] [2] . In order to exclude the integration over an infinite plane, one can use the formulations [3] [4] [5] utilizing the half-space Green’s function. In physical or mathematical, acoustics can be utilized for the study of solutions of the wave equation [6] . The wave phenomenon has been protracted to include even more particles and quasi-particles such as photons denoted by optical waves [7] [8] . On the other hand, phonons have represented by acoustic or elastic waves [9] . The acoustic or elastic wave phenomena of wave equation are the science of sound that deals with the generation of sound, propagation, and collaboration with the matter. Sources of quantum particles are known to be an indispensable element in all scattering experiments.

In previous, researchers have utilized different schemes to clarify the wave propagation of an acoustic wave. A phase-shift method is an accurate scheme for laterally invariant velocity structure but breaks down for arbitrary velocity variation, due to approximations made to the solutions of the wave equation [10] . In underwater acoustics, propagation, radiation, and scattering of sound are usually modeled by obtaining analytical and/or numerical solutions of a wave equation for the acoustic pressure, complemented by appropriate initial and boundary conditions [11] .

When the source moves even closer to the observer, then analysis of a wave field produce a non-asymptotic form of Green’s function that help us to elucidate near field analysis. It is vital to attain a correct causal form for the 3D transform, as otherwise the outline of forward and backward propagating waves does not follow logically [12] [13] [14] . This ability to treat modeling of Green’s function solution is achieved by making use of some early ideas on the non-iterative solution of the Lippmann-Schrödinger equation in quantum scattering [15] [16] as well as some very simple standard ideas from differential equations [17] . However, many studies have been already carried out on wave propagation throughout the world. In this work, we have considered the classical time-independent Schrödinger equation of finding the Green’s function solution for Helmholtz integral equation subject to Sommerfeld radiation conditions.

2. Green’s Function for Helmholtz Integral Wave Equation

Helmholtz boundary value problems can be put into one of three types: direct, eigenvalue or inverse. Among these methods, inverse is appropriate for analytical series methods that are established for a direct problem, with Dirichlet boundary conditions [18] . Helmholtz equation can be explained in many methods, but in our study, we introduce complex valued Green’s Function to solve non-homogeneous acoustic Helmholtz integral equation. Moreover, the general way to solve non-homogeneous acoustic Helmholtz integral equation for wave propagation is to form a Green’s function [19] . We used the one-dimensional (1D) Schrödinger classical wave equation as the following form

$\frac{{\partial }^{2}u}{\partial x^2}=\frac{1}{v^2}\frac{{\partial }^{2}u}{\partial t^2}$ (1)

Using separation of variables $u\left(x,t\right)=\psi \left(x\right)f\left(t\right)$ (where $\psi \left(x\right)$ is the Schrödinger wave function), the suitable wave equation solutions for $f\left(t\right)$ such as ${\text{e}}^{i\omega t}$ , take the following form

$\frac{{\text{d}}^2\psi \left(x\right)}{\text{d}{x}^2}=-\frac{{\omega }^2}{v^2}\psi \left(x\right)$ (2)

The above equation signifies an ordinary differential equation describing the spatial amplitude of the matter wave as a function of position. The entire energy of a particle is the sum of kinetic and potential energy and mathematically

$E=\frac{p^2}{2m}+V\left(x\right)$ (3)

But we know that $\omega =2\pi \upsilon$ and $v=\upsilon \lambda$ . Using these values in Equation (2), we get $\frac{{\text{d}}^2\psi \left(x\right)}{\text{d}{x}^2}+\frac{2m}{{\hslash }^2}\left[E-V\left(x\right)\right]\psi \left(x\right)=0$

or, $\left({\nabla }^2+k^2\right)\psi \left(r\right)=\lambda U\left(r\right)\psi \left(r\right)$ (4)

where, $U\left(r\right)=\frac{2m}{{\hslash }^2}V(\; r\; )$

With Green’s function $G\left(r,r^{\prime }\right)$ , we can write

$\left({\nabla }^2+k^2\right)G\left(r,r^{\prime }\right)=-\delta \left(r-r^{\prime }\right)$

or, $\left({\nabla }^2+k^2\right)G\left(R\right)=-{\delta }^3\left(R\right)$ (5)

This is Helmholtz integral equation. Where, $\delta$ means Dirac delta function and $R=r-r^{\prime }$ .

The Fourier integral transform is defined as

$\hat{G}\left({\xi }_1,{\xi }_2,{\xi }_3\right)={\displaystyle \underset{-\infty }{\iiint }G\left(x_1,x_2,x_3\right){\mathrm{e}}^{-i\left({\xi }_1{x}_1+{\xi }_2{x}_2+{\xi }_3{x}_3\right)}\text{d}{x}_1\text{d}{x}_2\text{d}{x}_3}$

$\Rightarrow \hat{G}\left(\xi \right)={\displaystyle {\int }_{R^3}G\left(x\right)\exp \left(-i\xi \cdot x\right)\text{d}x}$ (6)

By utilizing Fourier integral formula (6), the solution of the above Equation (5) is

$G\left(x\right)=\frac{1}{8{\pi }^3}{\displaystyle {\int }_{R^3}\hat{G}\left(\xi \right)\exp \left(i\xi \cdot x\right)\text{d}\xi }$ (7)

Assume that

$\underset{x_1\to \pm \infty }{\lim }\frac{{\partial }^n}{\partial x_1^n}G\left(x\right)=0,\left(n=0,1,2,\cdots \right)$ Sommerfeld radiation condition (8)

The Fourier Transform of $\frac{\partial }{\partial x_1}G\left(x\right)$ is defined by

$\begin{array}{l}={\displaystyle {\int }_{R^3}\frac{\partial }{\partial x_1}G\left(x\right){\text{e}}^{-i\xi \cdot x}\text{d}x}\\ ={\displaystyle {\int }_{-\infty }^{\infty }{\displaystyle {\int }_{-\infty }^{\infty }{\text{e}}^{-i\left({\xi }_2{x}_2+{\xi }_3{x}_3\right)}\text{d}{x}_2\text{d}{x}_3{\displaystyle {\int }_{-\infty }^{\infty }\frac{\partial }{\partial x_1}G\left(x\right){\text{e}}^{-i{\xi }_1{x}_1}\text{d}{x}_1}}}\\ ={\displaystyle {\int }_{-\infty }^{\infty }{\displaystyle {\int }_{-\infty }^{\infty }{\text{e}}^{-i\left({\xi }_2{x}_2+{\xi }_3{x}_3\right)}\text{d}{x}_2\text{d}{x}_3\times \left[\underset{M\to \infty }{\lim }{\left[G\left(x\right){\text{e}}^{-i{\xi }_1{x}_1}\right]}_{x_1=-M}^{x_1=M}-{\displaystyle {\int }_{-\infty }^{\infty }G\left(x\right)\left(-i{\xi }_1\right){\text{e}}^{-i{\xi }_1{x}_1}\text{d}{x}_1}\right]}}\end{array}$

(Integral by Parts)

$\begin{array}{l}=\left(i{\xi }_1\right){\displaystyle {\int }_{-\infty }^{\infty }{\displaystyle {\int }_{-\infty }^{\infty }{\text{e}}^{-i\left({\xi }_2{x}_2+{\xi }_3{x}_3\right)}\text{d}{x}_2\text{d}{x}_3}{\displaystyle {\int }_{-\infty }^{\infty }G\left(x\right){\text{e}}^{-i{\xi }_1{x}_1}\text{d}{x}_1}}\\ =\left(i{\xi }_1\right){\displaystyle {\int }_{R^3}G\left(x\right){\text{e}}^{-i\xi \cdot x}\text{d}x}\\ =\left(i{\xi }_1\right)\hat{G}(\; \xi \; )\end{array}$

Now the Fourier Transform of $\frac{{\partial }^2}{\partial x_1^2}G\left(x\right)$ is

$\begin{array}{l}={\displaystyle {\int }_{R^3}\frac{{\partial }^2}{\partial x_1^2}G\left(x\right)\exp \left(-i\xi \cdot x\right)\text{d}x}\\ =i{\xi }_1\frac{\partial }{\partial x_1}{\displaystyle {\int }_{R^3}\frac{\partial }{\partial x_1}G\left(x\right)\exp \left(-i\xi \cdot x\right)\text{d}x}\\ =i{\xi }_1{\displaystyle {\int }_{-\infty }^{\infty }{\displaystyle {\int }_{-\infty }^{\infty }{\text{e}}^{-i\left({\xi }_2{x}_2+{\xi }_3{x}_3\right)}\text{d}{x}_2\text{d}{x}_3}{\displaystyle {\int }_{-\infty }^{\infty }\frac{\partial }{\partial x_1}G\left(x\right){\text{e}}^{-i{\xi }_1{x}_1}\text{d}{x}_1}}\\ =\left(i{\xi }_1\right)\cdot \left(i{\xi }_1\right){\displaystyle {\int }_{-\infty }^{\infty }{\displaystyle {\int }_{-\infty }^{\infty }{\text{e}}^{-\left({\xi }_2{x}_2+{\xi }_3{x}_3\right)}\text{d}{x}_2\text{d}{x}_3}{\displaystyle {\int }_{-\infty }^{\infty }G\left(x\right){\text{e}}^{-i{\xi }_1{x}_1}\text{d}{x}_1}}\\ =-{\xi }_1^2{\displaystyle {\int }_{R^3}G\left(x\right)\exp \left(-i\xi \cdot x\right)\text{d}x}\\ =-{\xi }_1^2\hat{G}\left(\xi \right)\end{array}$ (9a)

In this similar way we can write

$\frac{{\partial }^2}{\partial x_2^2}G\left(x\right)=-{\xi }_2^2\hat{G}\left(\xi \right)$ (9b)

$\frac{{\partial }^2}{\partial x_3^2}G\left(x\right)=-{\xi }_3^2\hat{G}\left(\xi \right)$ (9c)

Plugging Fourier transform on both sides of equation (5), we can write

$F\left\{\left(\frac{{\partial }^2}{\partial x_1^2}+\frac{{\partial }^2}{\partial x_2^2}+\frac{{\partial }^2}{\partial x_3^2}+k^2\right)G\left(x_1,x_2,x_3\right)\right\}=-F\left\{\delta \left(x_1\right)\delta \left(x_2\right)\delta \left(x_3\right)\right\}$

$\Rightarrow {\displaystyle {\int }_{R^3}\left(\frac{{\partial }^2}{\partial x_1^2}+\frac{{\partial }^2}{\partial x_2^2}+\frac{{\partial }^2}{\partial x_3^2}\right)G\left(x\right)\exp \left(-i\xi \cdot x\right)\text{d}x}+k^2{\displaystyle {\int }_{R^3}G\left(x\right)\exp \left(-i\xi \cdot x\right)\text{d}x}=-1$

$\Rightarrow -\left({\xi }_1^2+{\xi }_2^2+{\xi }_3^2\right)\hat{G}\left(\xi \right)+k^2\hat{G}\left(\xi \right)=-1$

$\Rightarrow -{\xi }^2\hat{G}\left(\xi \right)+k^2\hat{G}\left(\xi \right)=-1$

$\Rightarrow \hat{G}\left(\xi \right)\left\{k^2-{\xi }^2\right\}=-1$ $\left[i^2=-1\right]$

Therefore,

$\hat{G}\left(\xi \right)=\frac{1}{{\xi }^2-k^2}$ (10)

The above equation satisfies the spherical symmetric potential and even function criteria respectively (Details in appendix A and B).

Now, the Equation (7) can be written as

$G\left(x\right)=\frac{1}{8{\pi }^3}{\displaystyle {\int }_{R^3}\frac{1}{{\xi }^2-k^2}\exp \left(i\xi \cdot x\right)\text{d}\xi }$

$\Rightarrow G\left(x\right)=\frac{1}{4{\pi }^2}{\displaystyle {\int }_0^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}\text{d}\xi }{\displaystyle {\int }_0^{\pi }{\mathrm{e}}^{i\xi r\cos \varphi }\sin \varphi \text{d}\varphi }$

$\Rightarrow G\left(x\right)=\frac{1}{4{\pi }^2}{\displaystyle {\int }_0^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}G\left(\xi \right)\text{d}\xi }$ ;

where, $G\left(\xi \right)=\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]$

$\Rightarrow G\left(x\right)=\frac{1}{4{\pi }^2}{\displaystyle {\int }_0^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}}\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]\text{d}\xi$

$\Rightarrow G\left(x\right)=\frac{1}{4{\pi }^2}\frac{1}{2}{\displaystyle {\int }_{-\infty }^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}}\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]\text{d}\xi$

Therefore

$G\left(x\right)=\{\begin{array}{l}\frac{1}{4\pi }\left(\frac{\cos kr+i\sin kr}{r}\right);I_m\left(k\right)>0\\ \frac{1}{4\pi }\left(\frac{\cos kr-i\sin kr}{r}\right);I_m\left(k\right)<0\end{array}$ (11)

(Details in appendix C) Equation (11) is the desired Green’s function solution for Helmholtz integral equation.

3. Results & Discussion

Our derived Green’s function solution for acoustic Helmholtz integral wave equation that yield good outcomes at adequately high wave numbers. This wave number k is often taking some real or complex constant, and Green’s function gets more accuracy for increasing the wave number, as depicted in Figures 1(a)-(d). The wave number can be measured complex if the medium of propagation is energy absorbing or a function of space. We have seen that, at very low wave numbers, it behaves very much like the Laplace equation and solutions at large wave number are highly oscillatory. This causes a great increase in complexity of analytical and numerical method. We also imposed Sommerfeld radiation condition or scattered condition for the solution of Helmholtz equation where confirms that the scattered wave is outgoing i.e. propagates away from the obstacle, as shown in the Figure 2(a) & Figure 2(b). Figure 3(a) & Figure 3(b)

Figure 1. Graphical Green’s function (a) 3D Plot of Green’s function for $k=100$ , $\mathrm{Im}\left(k\right)>0$ ; (b) 3D Plot of Green’s function for $k=100$ , $\mathrm{Im}\left(k\right)<0$ ; (c) Contour two dimensional; (d) curtain mesh plot of Green’s function when wave number $k=10000$ .

Figure 2. Scattering region of Helmholtz integral wave equation for wave number k = 50 and k = 500 respectively.

shows, the Helmholtz integral wave equation is an equation of the elliptic type, for which it is usual to consider boundary value problems. The wave front emitted from the source can no longer be regarded as a sphere. If the light emitted from a source exhibits sharp directivity then the wave front near the source point will act like a superimposed form of a plane wave and a spherical wave, but, it is clear that we need more and more time period at very high wave number to get accurate result for the solution of Helmholtz equation, as shown in Figure 4. Overall, Green’s function of Helmholtz integral equation for wave propagation are the solution of monochromatic mechanical and longitudinal plane waves (sink) when combined with ${\text{e}}^{-ikr}$ and diverging wave (source) when combined with ${\text{e}}^{ikr}$ respectively.

Figure 3. (a) Green’s function solution for Helmholtz integral equation when k = 10 and number of iteration n = 1000 and (b) when k = 25 and number of iteration n = 1000.

Figure 4. Green’s function solution for Helmholtz integral equation when k =200 and number of iteration n = 1000.

4. Conclusion

We have applied here an arguably more direct and general scheme to determine the generalized form of Green’s functions. A potential advantage of this approach is the fact that it may allow one to find a general form for Green’s functions of a given partial differential equation with a high quantum wave number k. The imposition of suitable Sommerfeld radiation condition allows us to find Green’s function solution that helps to represent wave propagation of Helmholtz integral equation. Our observation of this study allows us that monochromatic longitudinal plane wave is an idealization and represents a complex valued physical situation. The skills for proving the required facts about the involved generalized functions are likely to be beneficial or solution of many similar problems of acoustic, electrodynamics and mathematical physics.

Acknowledgements

This work is partially funded by Research and Development R & D Project with Grand No. 39.00.0000.012.20.011.22.56 (GO’s serial No. 15, Registration No. 23).

Appendix

A. Spherically Symmetric Potential Test

The Fourier transform of a function is spherically symmetric if and only if the function is spherically symmetric itself. This is one of the notorious properties of the Fourier transform.

Let k be defined as $k=\left\{\left({\xi }_1,{\xi }_2,{\xi }_3\right)\in {\Re }^3|{\xi }_1^2+{\xi }_2^2+{\xi }_3^2\le R,\left(R>0\right)\right\}$ . Also, let us consider

$V={\displaystyle \int {\displaystyle \int {\displaystyle {\int }_K\text{d}{\xi }_1\text{d}{\xi }_2\text{d}{\xi }_3}}}$ (1)

In order to attain a solution of acoustic Helmholtz integral equation subject to the Sommerfeld radiation condition, we need a spherically symmetric ansatz.

In spherical coordinate system, we have

$\begin{array}{l}x={\xi }_1=\rho \sin \varphi \cos \theta \\ y={\xi }_2=\rho \sin \varphi \sin \theta \\ z={\xi }_3=\rho \cos \varphi \end{array}$ (2)

Squaring and adding each of expressions in Equation (2), we get

$\begin{array}{c}{\xi }_1^2+{\xi }_2^2+{\xi }_3^2={\rho }^2{\sin }^2\varphi \left({\cos }^2\theta +{\sin }^2\theta \right)+{\rho }^2{\cos }^2\varphi \\ ={\rho }^2{\sin }^2\varphi +{\rho }^2{\cos }^2\varphi \\ ={\rho }^2\left({\sin }^2\varphi +{\cos }^2\varphi \right)\end{array}$

$\therefore {\xi }_1^2+{\xi }_2^2+{\xi }_3^2={\rho }^2$

Because of $\left({\cos }^2\theta +{\sin }^2\theta \right)=1$ and $\left({\sin }^2\varphi +{\cos }^2\varphi \right)=1$ , $\tan \theta =\frac{{\xi }_2}{{\xi }_1}$ and $\tan \varphi =\frac{\sqrt{{\xi }_1^2+{\xi }_2^2}}{{\xi }_3}$ .

Now, the determinant of the Jacobian matrix of $\left({\xi }_1,{\xi }_2,{\xi }_3\right)$ with respect to $\left(\rho ,\theta ,\varphi \right)$ is defined as

$J=\frac{\partial \left({\xi }_1,{\xi }_2,{\xi }_3\right)}{\partial \left(\rho ,\theta ,\varphi \right)}=\left|\begin{array}{ccc}\frac{\partial {\xi }_1}{\partial \rho }& \frac{\partial {\xi }_1}{\partial \theta }& \frac{\partial {\xi }_1}{\partial \varphi }\\ \frac{\partial {\xi }_2}{\partial \rho }& \frac{\partial {\xi }_2}{\partial \theta }& \frac{\partial {\xi }_2}{\partial \varphi }\\ \frac{\partial {\xi }_3}{\partial \rho }& \frac{\partial {\xi }_3}{\partial \theta }& \frac{\partial {\xi }_3}{\partial \varphi }\end{array}\right|$ (3)

Now, from Equation (21) we can write

$V={\displaystyle \int {\displaystyle \int {\displaystyle \int \text{d}{\xi }_1\text{d}{\xi }_2\text{d}{\xi }_3}}}={\displaystyle \int {\displaystyle \int {\displaystyle \int J\frac{\partial \left({\xi }_1,{\xi }_2,{\xi }_3\right)}{\partial \left(\rho ,\theta ,\varphi \right)}\text{d}\xi }}}$ (4)

Using Equations (3) in (4) and the range of spherical coordinates $\left(\rho ,\theta ,\varphi \right)$ are as follows $\rho \ge 0$ , $0\le \theta \le 2\pi$ , $0\le \varphi \le \pi$ . We get

$V={\displaystyle {\int }_{\rho =0}^R{\displaystyle {\int }_{\theta =0}^{2\pi }{\displaystyle {\int }_{\varphi =0}^{\pi }\left|\begin{array}{ccc}\sin \varphi \cos \theta & -\rho \sin \varphi \sin \theta & \rho \cos \varphi \cos \theta \\ \sin \varphi \sin \theta & \rho \sin \varphi \cos \theta & \rho \cos \varphi \sin \theta \\ \cos \varphi & 0& -\rho \sin \varphi \end{array}\right|\text{d}\rho \text{d}\theta \text{d}\varphi }}}$

$\begin{array}{l}\Rightarrow V={\displaystyle {\int }_{\rho =0}^R{\displaystyle {\int }_{\theta =0}^{2\pi }{\displaystyle {\int }_{\varphi =0}^{\pi }[\sin \varphi \cos \theta \left(-{\rho }^2{\sin }^2\varphi \cos \theta \right)}}}\\ +\rho \sin \varphi \sin \theta \left(-\rho {\sin }^2\varphi \sin \theta -\rho {\cos }^2\varphi \sin \theta \right)\\ +\rho \cos \varphi \cos \theta \left(-\rho \sin \varphi \cos \varphi \cos \theta \right)]\text{d}\rho \text{d}\theta \text{d}\varphi \end{array}$

$\begin{array}{l}\Rightarrow V={\displaystyle {\int }_{\rho =0}^R{\displaystyle {\int }_{\theta =0}^{2\pi }{\displaystyle {\int }_{\varphi =0}^{\pi }[-{\rho }^2{\sin }^3\varphi \left({\sin }^2\theta +{\cos }^2\theta \right)}}}\\ -{\rho }^2{\cos }^2\varphi \sin \varphi \left({\sin }^2\theta +{\cos }^2\theta \right)]\text{d}\rho \text{d}\theta \text{d}\varphi \end{array}$

$\Rightarrow V={\displaystyle {\int }_{\rho =0}^R{\displaystyle {\int }_{\theta =0}^{2\pi }{\displaystyle {\int }_{\varphi =0}^{\pi }\left[-{\rho }^2{\sin }^3\varphi -{\rho }^2{\cos }^2\varphi \sin \varphi \right]\text{d}\rho \text{d}\theta \text{d}\varphi }}}$ $\left[\because {\sin }^2\theta +{\cos }^2\theta =1\right]$

$\Rightarrow V={\displaystyle {\int }_{\rho =0}^R{\displaystyle {\int }_{\theta =0}^{2\pi }{\displaystyle {\int }_{\varphi =0}^{\pi }\left[-{\rho }^2\sin \varphi \right]\text{d}\rho \text{d}\theta \text{d}\varphi }}}$

$\Rightarrow V={\displaystyle {\int }_{\rho =0}^R{\displaystyle {\int }_{\theta =0}^{2\pi }{\displaystyle {\int }_{\varphi =0}^{\pi }{\rho }^2\sin \varphi \text{d}\rho \text{d}\theta \text{d}\varphi }}}$ $\left[\because \left|J\right|=\left|-{\rho }^2\sin \varphi \right|={\rho }^2\sin \varphi \right]$

$\Rightarrow V={\displaystyle {\int }_{\rho =0}^R{\displaystyle {\int }_{\theta =0}^{2\pi }{\rho }^2{\left[-\cos \varphi \right]}_0^{\pi }\text{d}\rho \text{d}\theta }}$

$\Rightarrow V=2{\displaystyle {\int }_{\rho =0}^R{\rho }^2{\left[\theta \right]}_0^{2\pi }\text{d}\rho }$

Therefore, $V=4\pi {\left[\frac{{\rho }^3}{3}\right]}_0^R=\frac{4}{3}\pi R^3$

Since V does not depends on $\theta$ and $\varphi$ but it is a function of R only. Hence, we can summarize that V is spherically symmetric. We can write

$G\left(x\right)=\frac{1}{4{\pi }^2}{\displaystyle {\int }_0^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}\exp \left(i\xi \cdot x\right)\text{d}\xi }{\displaystyle {\int }_0^{\pi }\sin \varphi \text{d}\varphi }$

$\Rightarrow G\left(x\right)=\frac{1}{4{\pi }^2}{\displaystyle {\int }_0^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}\text{d}\xi }{\displaystyle {\int }_0^{\pi }\exp \left(i\xi \cdot x\right)\sin \varphi \text{d}\varphi }$

$\Rightarrow G\left(x\right)=\frac{1}{4{\pi }^2}{\displaystyle {\int }_0^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}\text{d}\xi }{\displaystyle {\int }_0^{\pi }{\mathrm{e}}^{i\xi r\cos \varphi }\sin \varphi \text{d}\varphi }$ (5)

[since ${\text{e}}^{i\xi x}={\text{e}}^{i\xi r\cos \varphi }$ ]

To solve ${\displaystyle {\int }_0^{\pi }{\mathrm{e}}^{i\xi r\cos \varphi }\sin \varphi \text{d}\varphi }$ let us consider $z=i\xi r\cos \varphi$ , then $\text{d}z=-i\xi r\sin \varphi \text{d}\varphi$ and the limit of z takes the following form.

When, $\varphi =0$ then $z=i\xi r\cos 0=i\xi r$ and

when, $\varphi =\pi$ then $z=i\xi r\cos \pi =-i\xi r$

$\begin{array}{c}{\displaystyle {\int }_0^{\pi }{\mathrm{e}}^{i\xi r\cos \varphi }\sin \varphi \text{d}\varphi }={\displaystyle {\int }_{i\xi r}^{-i\xi r}{\text{e}}^z\frac{\text{d}z}{-i\xi r}}\\ =-\frac{1}{i\xi r}{\displaystyle {\int }_{i\xi r}^{-i\xi r}{\text{e}}^z\text{d}z}\\ =-\frac{1}{i\xi r}{\left[{\text{e}}^z\right]}_{i\xi r}^{-i\xi r}\\ =-\frac{1}{i\xi r}\left[{\text{e}}^{-i\xi r}-{\text{e}}^{i\xi r}\right]\end{array}$

Therefore,

${\displaystyle {\int }_0^{\pi }{\mathrm{e}}^{i\xi r\cos \varphi }\sin \varphi \text{d}\varphi }=\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]$ (6)

Using this in Equation (25), we can write

$G\left(x\right)=\frac{1}{4{\pi }^2}{\displaystyle {\int }_0^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}}\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]\text{d}\xi$ (7)

or, $G\left(x\right)=\frac{1}{4{\pi }^2}{\displaystyle {\int }_0^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}G\left(\xi \right)\text{d}\xi }$ (7a)

where, $G\left(\xi \right)=\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]$ .

B. Even Function Test

The function $G\left(x\right)$ of the above can be considered as a correct form if $G\left(\xi \right)$ satisfy even function property. i.e. $G\left(\xi \right)=G\left(-\xi \right)$ .

Now, let us consider,

$G\left(\xi \right)=\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]$ (8)

If we replace $\xi$ by $-\xi$ , then we can write

$G\left(-\xi \right)=\frac{1}{-i\xi r}\left[{\text{e}}^{-i\xi r}-e^{i\xi r}\right]$

$\Rightarrow G\left(-\xi \right)=-\frac{1}{-i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]$

$\Rightarrow G\left(-\xi \right)=\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]$

Since $G\left(-\xi \right)=G\left(\xi \right)$ . So, we can say that $G\left(\xi \right)=\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]$ is an even function.

C. Contour Integration (Jordan Lemma)

We have, $G\left(x\right)=\frac{1}{4{\pi }^2}{\displaystyle {\int }_0^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}}\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]\text{d}\xi$

$\Rightarrow G\left(x\right)=\frac{1}{4{\pi }^2}\frac{1}{2}{\displaystyle {\int }_{-\infty }^{\infty }\frac{{\xi }^2}{{\xi }^2-k^2}}\frac{1}{i\xi r}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]\text{d}\xi$

Therefore,

$G\left(x\right)=\frac{1}{8{\pi }^{2}ir}{\displaystyle {\int }_{-\infty }^{\infty }\frac{\xi }{{\xi }^2-k^2}}\left[{\text{e}}^{i\xi r}-{\text{e}}^{-i\xi r}\right]\text{d}\xi$ (9)

From the above equation, ${\xi }^2-k^2=0$ . Therefore, $\xi =\pm k$ . So, here we observe two pole points for which $G\left(\xi \right)$ are not analytic. To remove this, we need contour integration.

According to Cauchy’s residue theorem, we know

${\displaystyle {\int }_{\Gamma }G\left(\xi \right)\text{d}\xi }=2\pi i\times \text{sumoftheresidues}$ (10)

In equation (29), there has two parts. One is ${\displaystyle {\int }_{-\infty }^{\infty }\frac{\xi {\text{e}}^{i\xi r}}{{\xi }^2-k^2}\text{d}\xi }$ and another part of the form ${\displaystyle {\int }_{-\infty }^{\infty }\frac{-\xi {\text{e}}^{-i\xi r}}{{\xi }^2-k^2}\text{d}\xi }$ .

Case I: First, Consider ${\displaystyle {\int }_{-\infty }^{\infty }\frac{\xi {\text{e}}^{i\xi r}}{{\xi }^2-k^2}\text{d}\xi }$

Now, ${\displaystyle {\int }_{-\infty }^{\infty }\frac{\xi {\text{e}}^{i\xi r}}{{\xi }^2-k^2}\text{d}\xi }=\underset{R\to \infty }{\lim }{\displaystyle {\int }_{-R}^R+}{\displaystyle {\int }_{\Gamma }^{}\frac{\xi {\text{e}}^{i\xi r}}{{\xi }^2-k^2}\text{d}\xi }$ (11)

In the above, $\Gamma$ is the semi-circular arc of radius R. The poles have obtained by using the equation. The integrand has simple poles $\xi =\pm k$ . If $I_m\left(\xi \right)>0$ , then we can only see the location of the pole $\xi =+k$ . Because $\xi =+k$ only lies inside the contour C. On the other hand, $\xi =-k$ lies outside the contour C. As a result, we cannot see the location of the pole $\xi =-k$ .

Now, according to Cauchy residue theorem, we have

$\begin{array}{c}{I}_1=2\pi i\times \left(\text{Residueofupper-halfplaneat}\xi =+k\right)\\ =2\pi i\underset{\xi \to k}{\lim }\left\{\left(\xi -k\right)\frac{\xi {\text{e}}^{i\xi r}}{{\xi }^2-k^2}\right\}\\ =2\pi i\underset{\xi \to k}{\lim }\left\{\frac{1}{2}\left(\xi -k\right)\left(\frac{1}{\xi +k}+\frac{1}{\xi -k}\right){\text{e}}^{i\xi r}\right\}\\ =2\pi i\underset{\xi \to k}{\lim }\left\{\frac{1}{2}\left(\xi -k\right)\left(\frac{1}{\xi -k}\right){\text{e}}^{i\xi r}\right\}\\ =\pi i\underset{\xi \to k}{\lim }{\text{e}}^{i\xi r}=\pi i{\text{e}}^{ikr}\end{array}$

Taking $\underset{R\to \infty }{\lim }$ and using the property of contour integral $\underset{R\to \infty }{\lim }{\displaystyle {\int }_{\Gamma }^{}{\text{e}}^{i\xi r}F\left(\xi \right)\text{d}\xi }=0$ .

i.e. $\underset{R\to \infty }{\lim }{\displaystyle {\int }_{\Gamma }^{}\frac{\xi {\text{e}}^{i\xi r}}{{\xi }^2-k^2}\text{d}\xi }=0$

Therefore, $I_1=\pi i{\text{e}}^{ikr}$ (12)

Second, Consider ${\displaystyle {\int }_{-\infty }^{\infty }\frac{-\xi {\text{e}}^{-i\xi r}}{{\xi }^2-k^2}\text{d}\xi }$ (13)

In the above integral, a negative sign can be measured because with a positive parameter $\xi$ , then Jordan’s lemma states the following upper bound for the contour integral: where equal sign is when $G\left(\xi \right)$ vanishes everywhere. An analogous proclamation for a semicircular contour in the lower half-plane holds when $\xi <0$ .

Now, ${\displaystyle {\int }_{-\infty }^{\infty }\frac{-\xi {\text{e}}^{-i\xi r}}{{\xi }^2-k^2}\text{d}\xi }=\underset{R\to \infty }{\lim }{\displaystyle {\int }_{-R}^R+}{\displaystyle {\int }_{{\Gamma }^{\prime }}^{}\frac{-\xi {\text{e}}^{-i\xi r}}{{\xi }^2-k^2}\text{d}\xi }$ (14)

where ${\Gamma }^{\prime }$ is the semi-circular arc of radius R. In the similar process of the above, the integrand has simple poles at $\xi =\pm k$ . Here, $\xi =-k$ only lies inside the contour C. On the other hand, $\xi =k$ lies outside the contour C. As a result, we cannot see the location of the pole $\xi =k$ .

According to Cauchy residue theorem, we know

$\begin{array}{c}{I}_2=-2\pi i\times \left(\text{Residueoflowerhalf-planeat}\xi =-k\right)\\ =-2\pi i\underset{\xi \to -k}{\lim }\left\{\left(\xi +k\right)\frac{-\xi {\text{e}}^{-i\xi r}}{{\xi }^2-k^2}\right\}\\ =-2\pi i\underset{\xi \to -k}{\lim }\left\{\left(-\frac{1}{2}\right)\left(\xi +k\right)\left(\frac{1}{\xi +k}+\frac{1}{\xi -k}\right){\text{e}}^{-i\xi r}\right\}\\ =-2\pi i\underset{\xi \to -k}{\lim }\left\{\left(-\frac{1}{2}\right)\left(\xi +k\right)\left(\frac{1}{\xi +k}\right){\text{e}}^{-i\xi r}\right\}\\ =\pi i\underset{\xi \to -k}{\lim }{\text{e}}^{-i\xi r}=\pi i{\text{e}}^{-ikr}\end{array}$

In the similar process of the above, we can write

$I_2=\pi i{\text{e}}^{ikr}=\pi i{\text{e}}^{ikr}$ (15)

We have

$I=I_1+I_2=\pi i{\text{e}}^{ikr}+\pi i{\text{e}}^{ikr}=2\pi i{\text{e}}^{ikr}$

Using this value in Equation (9), we get

$G\left(x\right)=\frac{1}{8{\pi }^{2}ir}2\pi i{\text{e}}^{ikr}$

$\Rightarrow G\left(x\right)=\frac{1}{4\pi }\frac{{\text{e}}^{ikr}}{r}$

Therefore,

$G\left(x\right)=\frac{1}{4\pi }\left(\frac{\cos kr+i\sin kr}{r}\right)$ (16)

Case II:

First, Consider ${\displaystyle {\int }_{-\infty }^{\infty }\frac{\xi {\text{e}}^{i\xi r}}{{\xi }^2-k^2}\text{d}\xi }$

If $I_m\left(\xi \right)<0$ , then we get only the location of the pole $\xi =-k$ . In case of upper half-plane, $\xi =-k$ only lies inside the contour C. On the other hand, $\xi =+k$ lies outside the contour C. As a result, we cannot see the location of the pole $\xi =+k$ .

Now, according to Cauchy residue theorem, we have

$\begin{array}{c}{I}_1=2\pi i\times \left(\text{Residueofupper-halfplaneat}\xi =-k\right)\\ =2\pi i\underset{\xi \to -k}{\lim }\left\{\left(\xi +k\right)\frac{\xi {\text{e}}^{i\xi r}}{{\xi }^2-k^2}\right\}\\ =2\pi i\underset{\xi \to -k}{\lim }\left\{\frac{1}{2}\left(\xi +k\right)\left(\frac{1}{\xi +k}+\frac{1}{\xi -k}\right){\text{e}}^{i\xi r}\right\}\\ =2\pi i\underset{\xi \to -k}{\lim }\left\{\frac{1}{2}\left(\xi +k\right)\left(\frac{1}{\xi +k}\right){\text{e}}^{i\xi r}\right\}\\ =\pi i\underset{\xi \to -k}{\lim }{\text{e}}^{i\xi r}=\pi i{\text{e}}^{-ikr}\end{array}$

Therefore, we get

$I_1=\pi i{\text{e}}^{-ikr}$ (17)

In this alike process, according to Cauchy residue theorem, we know

$\begin{array}{c}{I}_2=-2\pi i\times \left(\text{Residueoflower-halfplaneat}\xi =k\right)\\ =2\pi i\underset{\xi \to k}{\lim }\left\{\left(\xi -k\right)\frac{-\xi {\text{e}}^{-i\xi r}}{{\xi }^2-k^2}\right\}\\ =2\pi i\underset{\xi \to k}{\lim }\left\{\left(-\frac{1}{2}\right)\left(\xi -k\right)\left(\frac{1}{\xi +k}+\frac{1}{\xi -k}\right){\text{e}}^{-i\xi r}\right\}\\ =\pi i\underset{\xi \to k}{\lim }\left\{\left(\xi -k\right)\left(\frac{1}{\xi -k}\right){\text{e}}^{-i\xi r}\right\}\\ =\pi i\underset{\xi \to k}{\lim }{\text{e}}^{-i\xi r}=\pi i{\text{e}}^{-ikr}\end{array}$

Therefore, we have

$I_2=\pi i{\text{e}}^{-ikr}$ (18)

Now $I=I_1+I_2=\pi i{\text{e}}^{-ikr}+\pi i{\text{e}}^{-ikr}=2\pi i{\text{e}}^{-ikr}$

Using this value, we have

$G\left(x\right)=\frac{1}{8{\pi }^{2}ir}2\pi i{\text{e}}^{-ikr}$

$\Rightarrow G\left(x\right)=\frac{1}{4\pi }\frac{{\text{e}}^{-ikr}}{r}$

Therefore,

$G\left(x\right)=\frac{1}{4\pi }\left(\frac{\cos kr-i\sin kr}{r}\right)$ (19)

So, we can write from Equations (36) and (39)

$G\left(x\right)=\{\begin{array}{l}\frac{1}{4\pi }\left(\frac{\cos kr+i\sin kr}{r}\right);I_m\left(k\right)>0\\ \frac{1}{4\pi }\left(\frac{\cos kr-i\sin kr}{r}\right);I_m\left(k\right)<0\end{array}$ (20)

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Delkhosh, M., Delkhosh, M. and Jamali, M. (2012) Greens Function and Its Applications. Journal of Basic and Applied Scientific Research, 2, 8865-8876. https://doi.org/10.1155/2012/180806
[2]	Van Gijzen, M.B., Erlangga, Y.A. and Vuik, C. (2007) Spectral Analysis of the Discrete Helmholtz Operator Preconditioned with a Shifted Laplacian. SIAM Journal on Scientific Computing, 29, 1942-1958. https://doi.org/10.1137/060661491
[3]	Seybert, A.F. and Soenarko, B. (1988) Radiation and Scattering of Acoustic Waves from Bodies of Arbitrary Shape in a Three-Dimensional Half Space. Journal of Vibration Acoustics Stress and Reliability in Design, 110, 112-117. https://doi.org/10.1115/1.3269465
[4]	Seybert, A.F. and Wu, T.W. (1989) Modified Helmholtz Integral Equation for Bodies Sitting on an Infinite Plane. The Journal of the Acoustical Society of America, 85, 19-23. https://doi.org/10.1121/1.397716
[5]	Shankar, R. (1994) Renormalization-Group Approach to Interacting Fermions. Reviews of Modern Physics, 66, 129-191. https://doi.org/10.1103/RevModPhys.66.129
[6]	Wang, Z., Gerstein, M. and Snyder, M. (2009) RNA-Seq: A Revolutionary Tool for Transcriptomics. Nature Reviews Genetics, 10, 57-63. https://doi.org/10.1038/nrg2484
[7]	Wang, Z. and Fan, S. (2005) Optical Circulators in Two-Dimensional Magneto-Optical Photonic Crystals. Optics Letters, 30, 1989-1991. https://doi.org/10.1364/OL.30.001989
[8]	Fleury, R., Sounas, D.L., Sieck, C.F., Haberman, M.R. and Alù, A. (2014) Sound Isolation and Giant Linear Nonreciprocity in a Compact Acoustic Circulator. Science, 343, 516-519. https://doi.org/10.1126/science.1246957
[9]	Gazdag, J. (1978) Wave Equation Migration with the Phase-Shift Method. Geophysics, 43, 1342-1351. https://doi.org/10.1190/1.1440899
[10]	Brekhovskikh, L.M. and Godin, O.A. (1999) The Lateral Wave. In: Brekhovskikh, L.M. and Godin, O.A., Eds., Acoustics of Layered Media II, Springer, Berlin, 81-119. https://doi.org/10.1007/978-3-662-03889-5_3
[11]	Jensen, F.B. (2003) Nitrite Disrupts Multiple Physiological Functions in Aquatic Animals. Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology, 135, 9-24. https://doi.org/10.1016/S1095-6433(02)00323-9
[12]	Lin, J., Rodríguez-Herrera, O.G., Kenny, F., Lara, D. and Dainty, J.C. (2012) Fast Vectorial Calculation of the Volumetric Focused Field Distribution by Using a Three-Dimensional Fourier Transform. Optics Express, 20, 1060-1069. https://doi.org/10.1364/OE.20.001060
[13]	Sheppard, C.J.R., Lin, J. and Kou, S.S. (2013) Rayleigh-Sommerfeld Diffraction Formula in k Space. JOSA A, 30, 1180-1183. https://doi.org/10.1364/JOSAA.30.001180
[14]	Kou, S.S., Sheppard, C.J.R. and Lin, J. (2013) Evaluation of the Rayleigh-Sommerfeld Diffraction Formula with 3D Convolution: The 3D Angular Spectrum (3D-AS) Method. Optics Letters, 38, 5296-5299. https://doi.org/10.1364/OL.38.005296
[15]	Sams, W.N. and Kouri, D.J. (1969) Noniterative Solutions of Integral Equations for Scattering. I. Single Channels. The Journal of Chemical Physics, 51, 4809-4814. https://doi.org/10.1063/1.1671871
[16]	Sams, W.N. and Kouri, D.J. (1970) Noniterative Solutions of Integral Equations for Scattering. IV. Preliminary Calculations for Coupled Open Channels and Coupled Eigenvalue Problems. The Journal of Chemical Physics, 53, 496-501. https://doi.org/10.1063/1.1674015
[17]	Kells, H.R. and Lear, S.A. (1960) Thermal Death Time Curve of Mycobacterium tuberculosis var. bovis in Artificially Infected Milk. Applied Microbiology, 8, 234-236. https://doi.org/10.1128/am.8.4.234-236.1960
[18]	Read, W.W. and Sneddon, G.E. (1995) Analytic Solutions for Solute Transport in Hillside Seepage. Proceedings 12th Australasian Fluid Mechanics Conference, Sydney, 10-15 December 1995, 191-194.
[19]	Khater, I. (2018) Analyse mathématique d’un modèle d’équation de réaction-diffusion à retard. Doctoral Dissertation, Depot institutionnel de l’Universite Abou Bekr Belkaid Tlemcen UABT.