Abstract
The transmission of information by propagating or diffusive waves is common to many fields of engineering and physics. Such physical phenomena are governed by a Helmholtz (real wavenumber) or pseudo-Helmholtz (complex wavenumber) equation. Since these equations are linear, it would be useful to be able to use tools from signal theory in solving related problems. The aim of this paper is to derive multidimensional input/output transfer function relationships in the spatial domain for these equations in order to permit such a signal theoretic approach to problem solving. This paper presents such transfer function relationships for the spatial (not Fourier) domain within appropriate coordinate systems. It is shown that the relationships assume particularly simple and computationally useful forms once the appropriate curvilinear version of a multidimensional spatial Fourier transform is used. These results are shown for both real and complex wavenumbers. Fourier inversion of these formulas would have applications for tomographic problems in various modalities. In the case of real wavenumbers, these inversion formulas are presented in closed form, whereby an input can be calculated from a given or measured wavefield.
1. Introduction
The transmission of information over space and time is often governed by the theory of waves. Many important physical phenomena are described by a Helmholtz equation, for example, in the fields of electromagnetism, acoustics, and optics [1–4]. Other physical phenomena are well described by the use of damped or diffusion waves, such as the propagation of heat or photonic waves. For example, ultrasound has long been used for medical imaging [3] while the use of optical radiation via diffuse photon density waves for imaging inhomogeneities in turbid media is a newer development [5–8]. Similarly, photothermal tomographic imaging methods have also been used for nondestructive evaluation [9–12], and more recently for biomedical imaging [13, 14]. Furthermore, a combination of ultrasonic and optical techniques known as photoacoustics has also gained prominence in recent years [4, 15].
Modelled as a linear process, the propagation and scattering of waves obey the principles of superposition and homogeneity. One common method of studying linear processes is to view them as linear systems with an input and output. The input/output relationships of the system can then be characterized via the system transfer function, which implies interpreting the problem as a signal theory problem and thus enabling the application of powerful system theoretic results and concepts. This type of approach applied to propagating or diffusion wave problems has already yielded promising results [16–18] and the goal of this paper is to expand the scope of this approach.
The seeds of such an approach to the solution of tomographic-type problems can be seen in [1, 8, 10, 19, 20]; however, in those papers there is no clear reference to any such system theoretic concepts. As such, a unifying framework for a signal theoretic approach to wavefields does not exist. The goal of this work is to propose a clear and cohesive signal theoretic framework for working with acoustic, thermal, photonic, or other wavefield-related problems. This would be particularly powerful, enabling the application of a vast array of results from a signal processing point of view and in particular, recent novel results from algebraic signal processing [21–27] should find ground-breaking applications.
With these ideas in mind, the goal of this exposition is to derive multidimensional spatial transfer function relationships for the Helmholtz and pseudo-Helmholtz equation so that the tools of linear system theory can be used to solve related problems.
We consider the input quantity to be the time- and space-dependent inhomogeneity term (the “forcing” term in the relevant partial differential equation). The output quantity is then taken as the resulting field present at any point in space. The goal is to develop transfer function relationships in spatial multidimensions so that general results of the theory of signals and systems may be used. It is known that the relationship between the input (inhomogeneity term) and the output (resulting wavefield) is one of convolution with the impulse response (Green's function in classical terms). In the Fourier domain, the equivalent relationship is a transfer function relationship, that is to say one of multiplication. In this paper, key transfer function relationships are shown to hold for the spatial (not Fourier) domain with the use of appropriate coordinate systems. That is to say, the relationship between input and output is shown to be one of multiplication (not convolution) even in the spatial domain, hence yielding the signal theoretic construction that is the aim of the paper.
Rather than focusing on a specific application area (acoustics, medical imaging, seismic, etc.), the end goal of this paper is a unifying framework. The backbone of this framework is Fourier transforms in multidimensions which are also transcribed into curvilinear coordinate systems. The representations in this paper are embedded in a time-dependent, spatially one- to three-dimensional description.
The outline of this paper is as follows. Sections 2 to 7 establish the background preliminaries, definitions, and sign conventions that are important to the developments in the rest of the paper. The definition of multidimensional Fourier transforms is given in Section 2. Section 3 presents the mathematical forms of the types of signal to which this development applies, namely, any signal governed by a Helmholtz or pseudo-Helmholtz equation. Section 4 introduces the Green's and transfer functions for the (generalized) Helmholtz equation and points out that the spatial transfer function assumes a particularly simple form in curvilinear coordinates. Section 5 explains the sign conventions that will be employed in this paper to ensure consistency and simplicity of results. Sections 6 and 7 develop multidimensional Fourier transforms in curvilinear coordinates. Sections 8, 9, and 10 present the relevant transfer function relationships in curvilinear coordinates in 3D, 2D, and 1D, respectively, and form the core results of this paper. Section 11 summarizes the results of the previous three sections in a tabular format that is easy to refer to. Sections 12 and 13 present some applications while Section 14 concludes the paper.
2. Fourier Transforms in -Dimensions
The theory of Fourier transforms can be extended to in a way that is completely analogous to the treatment of one-dimensional Fourier transforms. Let elements of be denoted by or more generally , and elements of the corresponding Fourier dual space be denoted by . Under the suitability of integration of , then the Fourier transform in is defined for in as Under suitable conditions, the function can be recovered from the inverse transform through Various other conventions are possible regarding the location of the positive and negative signs and also the factors of in (2.1) and (2.2).
3. The Helmholtz and Pseudo-Helmholtz Equation
All wave fields governed by the wave equation (such as acoustic or electromagnetic waves) lead to the Helmholtz equation once a Fourier transform is used to transform the time domain to the frequency domain: where . Here, is the temporal Fourier transform of the inhomogeneous time- and space-dependent source term for the wave equation. From a signal theory perspective, this is considered to be the input to the system. The variable represents a physical variable that is governed by the wave equation, for example, acoustic pressure, and is considered to be the output from a signal theory point of view. Both and are functions of position, , and (temporal) frequency, . The variable represents the speed of the wave, which for an acoustic wave would be the speed of sound and for an electromagnetic wave would be the speed of light. For wavefields governed by the wave equation, is a real (and positive) quantity.
Other physical phenomena lead to a “pseudo-” Helmholtz equation upon Fourier transformation of the time variable to a (temporal) frequency variable. For example, the equation for a diffuse photon density wave (DPDW) which describes the photon density in a solid due to incident energy intensity (optical source function) is given in the time domain by [2]
In the above equation, is the optical absorption coefficient (m−1), is the speed of light in the turbid medium (m/s), and is the optical diffusion coefficient (m). Taking the temporal Fourier transform (FT) of (3.2) to transform time to frequency gives the following: where is the complex wavenumber. Equation (3.3) is a pseudo-Helmholtz equation. The term “pseudo” is used because although (3.3) has the form of a Helmholtz equation similar to (3.1), the wavenumber is a complex variable, with the imaginary part indicating a decaying or damped wave. In this development the terms “Helmholtz” and “pseudo-Helmholtz” will generally be used interchangeably unless otherwise noted. Similarly, the standard heat equation is given by where is the time- and space-dependent heat source, describes the temperature in the material as a function of time and space, and is the thermal diffusivity of a material. As done previously, we take the temporal FT of (3.4) to obtain where is the complex wavenumber. Once again (3.5) has the form of a Helmholtz equation, although it is a pseudo-Helmholtz equation to be exact since the wavenumber is complex and indicates a damped wave.
Thus, we see that the general Helmholtz form of (3.1) can be used to describe several different physical phenomena, from the propagation of light or acoustic waves to the heavily damped nature of photonic or thermal waves. The exact form of the wavenumber in each case is the best indicator as to the propagation characteristics of a wave with a real indicating a propagating wave and a complex indicating a damped wave.
4. Green's Function and Transfer Function for the Helmholtz Equation
Taking the full spatial Fourier transform of the Helmholtz (or pseudo-Helmholtz) equation and rearranging yields where we have used the shorthand notation to denote a point in 3D spatial Fourier frequency space. The wavenumber may be real or complex. A capital letter is used to denote the full 3D spatial and temporal Fourier transform of a function, although it should be clear from the arguments which is being indicated. In cases where it may not be clear, a tilde will be used to denote the function in spatial Fourier space. For shorthand notation, let , so that is the length of the spatial Fourier vector. The Fourier space has the same spatial frequency dimension as the spatial dimension of the problem. So in 2D space, a spatial frequency variable would not be required and so forth.
By inverse spatial Fourier transformation of (4.1), the wavefield in dimensions is given by Using the definition of the Fourier transform of , the above equation can be rewritten as Let us define the spatial Green’s function in -dimensional space as This is the Green's function for the (pseudo-) Helmholtz equation. By switching the order of integration, (4.3) can be rewritten such that it can be clearly interpreted as a spatial convolution of Green's function with the input source: The notation has been used to denote a (multidimensional) space-only convolution. The frequency domain equivalent to (4.5) is (4.1), which can be interpreted as a multiplication in terms of the spatial transfer function as where the Fourier transform of the spatial Green's function is the spatial transfer function: The dependence of the spatial Green's function on temporal frequency is via the wavenumber of the wavefield (which may be real or complex) and this is emphasized by writing as a function of spatial and temporal frequencies.
4.1. Spatial Transfer Function in Curvilinear Coordinates
Equations (4.4) and (4.7) hold in any -dimensional space with Green's function of (4.4) being the -dimensional spatial inverse Fourier transform of the spatial transfer function of (4.7). More specifically, in 1D, , and in 2D . Transforming to polar coordinates in the 2D spatial frequency domain via the transformation , gives where is the radial frequency variable in polar coordinates. Similarly, in 3D, becomes where is the spherical radial frequency variable in spherical polar coordinates in spatial frequency space via the coordinate transformation given by , and . Thus, the spatial transfer functions assume particularly simple, radially symmetric forms in polar coordinates for 2D space and spherical polar coordinates for 3D space. These particularly simple forms for the 2D and 3D transfer functions in polar coordinates motivates the development of the general 2D and 3D Fourier transforms in curvilinear coordinates [28, 29] so that this simple form of the transfer function may be exploited.
5. Notation and Sign Conventions
Wavenumbers are defined by the squares of their quantities and arise as a result of taking the Fourier transform of the corresponding propagation equation (be it acoustic, thermal, or otherwise), which in turn leads to a Helmholtz or pseudo-Helmholtz equation. These wavenumbers are defined so that the wave propagation equation transformed to the temporal frequency domain all have the form of a pseudo-Helmholz equation as given by (3.1) or (3.5) with a complex or real wavenumber. For the rest of the paper, we will consider that represents a generic complex wavenumber while represents a generic real wavenumber.
Rather than the squared wavenumber, the quantity of interest will prove to be the wavenumber itself, namely, and , which are the square roots of the given squared wavenumber in the Helmholtz equation. Each can be considered as the sum of a real and an imaginary part, so that with denoting the real part of and denoting the imaginary part. Since the square root of any can be , we will use the convention that for a given complex , the required square root of the corresponding is defined such that the imaginary part of is negative. Hence, is chosen as the square root of such that and so forth. If this sign convention is adopted, then it was shown in [36] that the many results for complex wavenumber use the same mathematical form of travelling wave solution as for the real wavenumbers.This makes the notation and book-keeping considerably simpler. It is noted that this sign convention is the opposite of what this author adopts in [30].
5.1. Sommerfeld Radiation Condition
To aid in the selection of a causal solution, the Sommerfeld radiation condition is required. The Sommerfeld radiation condition states that the sources in the field must be sources not sinks of energy. Therefore, energy radiated from sources must scatter to infinity and cannot radiate from infinity into the field. Mathematically, a solution , where is the spatial variable, to the Helmholtz equation is considered to be radiating if it satisfies where is the dimension of the space and is the wavenumber in the Helmholtz equation. This is the radiation equation based on an implied time variation of which is implicit in our chosen definition of the Fourier transform, (2.1). For example, with the standard definition of the Fourier transform, the transform of is , clearly implying the dependence. Had a different definition of the Fourier transform been used, the implied time variation would have been , in which case the sign of in (5.1) would be reversed.
6. 2D Fourier Transforms in Terms of Polar Coordinates
Given the desire to exploit the simple form of the spatial transfer function in curvilinear coordinates, we consider 2D and 3D Fourier transform in terms of curvilinear coordinates. Let us first consider the Fourier transforms in 2D in polar coordinates. In order to define this, some preliminary definitions of Hankel transforms are required first.
6.1. Hankel Transforms
The Hankel transform of order is defined by the integral [31] where is the th order Bessel function. If , the transform is self-reciprocating and the inversion formula is given by
6.2. Connection between the 2D Fourier Transform and Hankel Transform
To develop Fourier transforms in 2D in polar coordinates, both the function and its Fourier transform are expressed in polar coordinates. In general is not radially symmetric and is a function of both and so that the dependence can be expanded into a Fourier series due to the periodicity of the function in : where the Fourier coefficients can be found from Similarly, the 2D Fourier transform can also be expanded into its own Fourier series so that and where those Fourier coefficients are similarly found from
The relationship between the Fourier coordinates in normal space in (6.3) and the spatial frequency space Fourier coordinates in (6.5) is desired. The details of this are given in [28], and the results are summarized here. It is emphasized that the relationship between and is not a Fourier transform. The relationship between them is given by where is the th order Hankel transform. The reverse relationship is given by The th term in the Fourier series for the original function will Hankel transform into the th term of the Fourier series of the Fourier transform. However, it is an th order Hankel transform for the th term, namely, all the terms are not equivalently transformed. The mapping from to is one of th-order Hankel transform, which in general is not a 2D Fourier transform.
The operation of taking the 2D Fourier transform of a function is thus equivalent to (1) first finding its Fourier series expansion in the angular variable and (2) then finding the th-order Hankel transform (of the radial variable to the spatial radial variable) of the th coefficient in the Fourier series. Since each of these operations involves integration over one variable only with the others being considered parameters vis-à-vis the integration, the order in which these operations are performed is interchangeable.
7. Spherical Hankel Transform
We introduce the spherical Hankel transform, which will form part of the 3D Fourier transform in spherical coordinates.
7.1. Definition of the Spherical Hankel Transform
The spherical Hankel transform can then be defined as [32, 33] is used to specifically denote the spherical Hankel transform of order . The inverse transform is given by [29] The spherical Hankel transform is particularly useful for problems involving spherical symmetry.
7.2. Spherical Harmonics
The spherical harmonics are the solution to the angular portion of Laplace's equation in spherical polar coordinates and can be shown to be orthogonal. These spherical harmonics are given by [32] where is called a spherical harmonic function of degree and order , is an associated Legendre function, represents the colatitude and represents the longitude. With the normalization of the spherical harmonics as given in (7.3), the spherical harmonics are orthonormal so that Here is the kronecker delta and the overbar indicates the complex conjugate. It is important to note that there are several different normalizations of the spherical harmonics that are possible, that will differ from (7.3) and thus lead to a different version of (7.4). Several of these are nicely catalogued in [34], including which disciplines tend to use which normalization. This is important to note since any result that uses orthogonality will differ slightly depending on the choice of normalization.
The spherical harmonics form a complete set of orthonormal functions and thus form a vector space. When restricted to the surface of a sphere, functions may be expanded on the sphere into a series approximation much like a Fourier series. This is in fact a spherical harmonic series. Any square-integrable function may be expanded as where The coefficients are sometimes referred to as the spherical Fourier transform of [35].
7.3. 3D Fourier Transforms in Spherical Polar Coordinates
To find 3D Fourier transforms in spherical polar coordinates, both the function and its 3D Fourier transform are written in terms of polar coordinates in the spatial and spatial frequency domains. That is, a function in 3D space is expressed as as a function of spherical polar coordinates and its 3D Fourier transform is written as in frequency spherical polar coordinates. Certain relationships can be shown to hold. The relationship between the function and its transform are summarized here, with the relevant details omitted for brevity. The function itself can be expanded as a series in terms of the spherical harmonics as where are the spherical harmonics and the Fourier coefficients are given by The 3D Fourier transform of can also be written in Fourier space in spherical polar coordinates as where We emphasize again that the relationship between and is not that of a Fourier transform. In fact, this relationship is given by [29] where denotes a spherical Hankel transform of order . The inverse relationship is given by
8. Spatial Transfer Function Relationship in 3D
With this long set of preliminaries, definitions, and conventions established, we can now begin the derivation of the 3D spatial transfer function relationship between input and output . The foregoing section makes use of several key results in [36] which will be presented here without proof. The results in [36] address general cases, not all of which are relevant to the work here. The versions of the theorems required herein are stated based on the sign conventions presented in Section 5 with regards to the sign conventions for the complex wavenumbers and also for the chosen definition of the Fourier transform. These affect the implied time dependence, the Sommerfeld radiation condition, and thus the physically correct solution that is chosen from several mathematical possibilities. The versions of the theorems presented here ensure that the chosen waves are outwardly propagating (Sommerfeld radiation condition) and also decay to zero at infinity for a damped wavefield, thus ensuring bounded and physically meaningful solutions.
Theorem 8.1. It is shown in [36] that the following result holds true: Here, is any analytic function defined on the positive real line that remains bounded as goes to infinity, is a spherical Bessel function of order , is a spherical Hankel function of order , and is a wavenumber which may be real or complex, chosen such that the imaginary part of is negative. Given the definition of the Fourier transform that is being currently used, the presented result satisfies the Sommerfeld radiation condition, ensuring an outwardly propagating wave.
8.1. 3D Transfer Function Relationship in Spherical Polar Coordinates
From (4.7) and using the conversion to spherical polar coordinates in 3D, the transfer function for the Helmholtz equation in the Fourier domain can be written as which depends on the frequency spherical radial variable only.
In general, as the imaginary part of the wavenumber gets smaller, the transfer function becomes more frequency selective, in the sense of strongly passing a smaller bandwidth. As the imaginary part of the wavenumber gets larger, the transfer function becomes more low-pass in nature, passing lower frequencies more strongly than the higher frequencies, with a corresponding larger bandwidth of frequencies being passed. This is physically meaningful as the imaginary part of represents the damping inherent in this wave modality. Wavenumbers with larger imaginary parts are more heavily damped and that damping tends to affect the higher frequencies more—that is, the higher frequencies are attenuated and the lower ones are passed through, implying a low-pass nature to the physical system.
Recall that is the magnitude of the spatial frequency variables. The case for the imaginary part of equal to zero (in other words, a real wavenumber) is easy to visualize. For the case of a real , the function is discontinuous at , meaning that only those spatial frequencies where exactly are passed and the rest are attenuated. In other words, the resulting wave must have . While these statements are fairly mathematical in nature, they explain the nature of the resolution of various imaging modalities. For example, it is far more difficult to achieve the resolution of acoustic imaging with thermal imaging and one of the explanations for this is the “low-pass” blurring nature of the complex wavenumber of a thermal wave versus the highly selective band-pass nature of the acoustic wavenumber. This is discussed further in other papers, for example, in [18].
8.2. 3D Green Function Coefficients
Let us define a set of functions that will be referred to as the Green function coefficients. The actual Green function for the system is the full 3D inverse Fourier transform of the transfer function, which for a spherically symmetric function is equivalent to a spherical Hankel transform of order zero only [33]. With the help of Theorem 8.1, we define these Green function coefficients as where has been used to denote the th order Green function coefficient, with the subscript indicating the order of the coefficient and the wavenumber included as a parameter of the function.
8.3. Transfer Function Relationship in Spherical Polar Coordinates
The 3D Fourier transform of the input (source) function is written in polar coordinates and expanded in terms of a spherical harmonic series as
Hence, (4.1) for the output wavefield becomes The symbol is used to denote the wavenumber. The spatial inverse Fourier transform and thus the output in spatial coordinates is given by The temporal frequency is a constant as far as the spatial inverse Fourier transformation is concerned. We recognize that the quantity within the curly brackets can be evaluated with the help of Theorem 8.1 and using the definition of the Green function coefficients in (8.3) to obtain It now follows that the wavefield expression becomes If the measured wavefield itself, namely, the left-hand side of (8.8), is expanded as in a spherical harmonic series so that then (8.8) gives us the simple input-output transfer function relationship we seek: Note how this relationship is in the spatial domain (not the frequency spatial domain) and gives a multiplicative transfer function relationship between the input coefficients and the resulting wavefield .
Using the relationships proposed between the coefficients of the function in the spatial domain and the coefficients in the Fourier domain as given in (7.11), the general relationship between Hankel and Fourier transforms is given by
Hence (8.10) can be written as so that the value of the measured wavefield is related to the spherical Hankel transform of the input function, evaluated at the wavenumber of the wavefield in question. Using the definition of the spherical Hankel transforms, this becomes a very compact yet powerful expression: The value of defining the Green function coefficients now becomes apparent in that it permits the relationship between and to have a simple multiplicative transfer function relationship. In other words, the relationship between output wave coefficients and input is one of multiplication with the transfer function coefficients instead of a complicated convolution-type relationship.
If we also assume that the input function is separable in time and space (a fairly general assumption) so that . This would be the case for a spatial inhomogeneity and a temporal input. In this case, then and the formula for the forward problem is given by
8.4. Discussion and Relationship with the Fourier Diffraction Theorem
Equation (8.13) is the key input-output relationship that is sought. It gives the relationship between input and output and the relationship is a transfer function (multiplication) type of relationship, even in the spatial domain where the normal relationship to be expected is one of convolution. This relationship is for the spherical harmonic expansion coefficients, not between the full functions themselves. However, because it is a direct, proportional type of relationship between the th term of the input and the th term of the output, it is particularly useful and simple to apply.
In particular, (8.13) states that the wavefield (output) coefficients are directly proportional to the transfer function coefficients and the proportionality factor between them is the spherical Hankel transform of the input coefficients, evaluated on the sphere . Note that th order output coefficients are related to th order transfer function coefficients in proportion to the th order spherical Hankel transform of the th order input coefficient. This is in fact a generalization of the Fourier diffraction theorem of tomography [3] which loosely states that the output wavefield is proportional to the Fourier transform of the input inhomogeneity evaluated somewhere on the sphere in Fourier space. This is still true but we have expressed a more precise version of this theorem in the sense that we have removed any ambiguity regarding the “somewhere” on the sphere and replaced it with a relationship (8.14) that does not depend on angular location. This was enabled by the use of the Fourier transform in curvilinear coordinates so that a full Fourier transform requires an set of Fourier coefficients and spherical Hankel transforms. Furthermore, the idea of proportionality between Fourier transform of the inhomogeneity and output wavefield is also made more precise by stating that the proportionality between them is actually the coefficients of the Green function itself. In essence, the output wavefield can be seen as being the Green function coefficients, with each term weighted by the spherical Hankel transform of the input inhomogeneity.
9. 2D Transfer Function Relationship in 2D Polar Coordinates
We proceed to develop the transfer function expressions for the 2D case in polar coordinates. As for the 3D case, we begin with a necessary theorem which is stated without proof. As previously mentioned, the results in [36] address general cases, and the versions of these theorems required herein are stated here based on the sign conventions presented in Section 5 for the complex wavenumbers and also for the Fourier transform which affects the implied time dependence and thus the Sommerfeld radiation condition.
Theorem 9.1. It is shown in [36] that the following result holds true: Here, is an analytic function defined on the positive real line that remains bounded as goes to infinity, is a Bessel function of order , is a Hankel function of order , and is a wavenumber which may be real or complex, chosen such that the imaginary part of is negative. Given the definition of the Fourier transform that is being currently used, the presented result satisfies the Sommerfeld radiation condition, ensuring an outwardly propagating wave.
9.1. 2D Green Function Coefficients
The overall system transfer function in 2D is given by which depends on the frequency radial variable only.
As for the 3D case, we define the 2D Green function coefficients as the th-order inverse Hankel transform of the overall system spatial transfer function. These will be shown to be the required transfer functions for working in the 2D polar formulation. From the definition of the inverse Hankel transform, these are given by Recall that the actual (full) Green function for the system is the full 2D inverse Fourier transform of the transfer function, which is an inverse Hankel transform of order zero only. In (9.3), has been used to denote the nth order Green function coefficients with the subscript indicating the order.
9.2. Fourier Theorem for 2D Fourier Transforms in Polar Coordinates
The 2D Fourier transform of the input function is written in polar coordinates as where can be found from The output wavefield is similarly expanded as a series so that in the spatial domain we can write while in the spatial Fourier domain this is The output wavefield is given by and since is radially symmetric and does not require a full series, the output wavefield in the Fourier domain is given by or equivalently as The equivalent expression to (9.8) in the spatial domain is given by or The temporal frequency is a constant as far as the spatial inverse Fourier transformation is concerned. The quantity within the curly brackets can be evaluated with the help of Theorem 9.1 as Hence, the output wavefield expression is given by The definition of forward and inverse Hankel transforms, the general relationship between Hankel and Fourier transforms is given by and may be used to simplify (9.13) to yield or more compactly as It is in (9.15) that the interpretation of the evaluation of the 2D Fourier transform becomes apparent, namely through the evaluation of the Bessel function at the (real or complex) wavenumbers. Using the definition of the Hankel transform, (9.16) can be written in the compact and powerful formulation of If it is also further assumed that the inhomogeneity function is separable in time and space (a fairly general assumption) so that , and then (9.16) can be further reduced to It is noted that the value of the measured wavefield at any position is related to the Hankel transform of the heterogeneity function evaluated at the wavenumber of the wavefield in question, with the Green function coefficient acting as the proportionality term.
Equation (9.19) is the 2D equivalent of (8.14). The same comments can be made as those made in the discussion after the 3D version of this problem, as the relationship is identical, save for translating 3D tools such as the spherical harmonic expansions and spherical Hankel transforms into 2D tools involving Fourier series and Hankel transforms.
10. Transfer Function Relationship in 1D
We proceed to develop the equivalent relationship for the 1D case. We have seen from (4.7) that the 1D transfer function in the Fourier domain is given by which depends on only a single spatial frequency. Equation (4.2) in 1D then becomes The temporal frequency is a constant as far as the spatial inverse Fourier transformation is concerned. We are thus interested in calculating integrals of the form where is an analytic function defined on the positive real line that approaches zero as goes to infinity. The required theorem is given in [36] and the relevant result is presented below.
Theorem 10.1. It is shown in [36] that the following is true Here, is an analytic function defined on the positive real line that remains bounded as goes to infinity and is a wavenumber which may be real or complex, chosen such that the imaginary part of is negative. Given the definition of the Fourier transform that is being currently used, the presented result satisfies the Sommerfeld radiation condition, ensuring an outwardly propagating wave.
10.1. 1D Green's Functions via Inverse Fourier Transformation of the Spatial Transfer Functions
The Green's function for the Helmholtz equation is 1D and is given by inverse Fourier transform of the spatial transfer function as and can be evaluated with the help of Theorem 10.1 to give
10.2. Transfer Function Relationship in 1D
Equation (10.2) can now be evaluated with the help of Theorem 10.1 as well as the result in (10.6) as where . This is the 1D version of the transfer function relationship, with the wavefield being directly proportional to the Fourier transform of the object function evaluated at . This is similar to the results for the 2D and 3D cases. For a real wavenumber, a “sphere” in 1D becomes the two points on the real line at and the proportionality term is the Green's function for the space in question.
As before, we assume that the inhomogeneity function is separable in time and space (a fairly general assumption) so that , where is the 1D Fourier transform of . The relevant result now reads This result is directly applicable for computational purposes.
11. Summary of Results
The relationships in the previous sections are summarized in Table 1.
12. Applications of the Transfer Function Relationships to the Wave Equation
The transfer function relationships can be used to find time domain Green's functions where the temporal Fourier integral can be easily inverted. In particular, we consider the case where the input to the standard wave equation is a Dirac-delta function at the origin
This corresponds to the Helmholtz equation above with . The temporal and spatial Fourier transform of is 1, which is spherically symmetric so that the transfer function relationships given above only need the zeroth-order component and simplify considerably. In 3D, the transfer function relationship is given by In 2D, the relationship is given by In 1D, the relationship is These three relationships can now be inverse Fourier transformed in time.
For the 3D case, we use the fact that so that the inverse Fourier transform of (12.2) gives The 1D time response can be found from the inverse Fourier transform of (12.4) where is the Heaviside unit step function.
The 2D case is a little more complicated but can nevertheless be evaluated in closed form by an inverse Fourier transform of (12.3). To do this, we write the zeroth-order Hankel function in its integral form as [37] Hence from (12.7) and (12.3), the inverse Fourier transform of (12.3) gives Changing the order of integration gives But so that (12.9) becomes Changing variables so that gives which yields the desired result in closed form.
13. Application: Inversion Formulas for Tomographic Applications with Real Wavenumbers
In the case of a real wavenumber , an immediate application of the above formulas is for closed-form inversion formulas which are useful for tomographic applications. These inversion formulas are applicable to applications where the wavefield (e.g., acoustic field) is measured and the goal is to reconstruct the source (input) that led to that measured wavefield (output). As we have seen, formulas above give the output (wavefield) in terms of the transform of the input evaluated at . Thus, if the output wavefield is measured, this implies knowledge of the transform of the input. An inverse transform then leads to the input itself.
13.1. Inversion Formula in 3D
Equation (8.14) admits an inversion useful for tomographic applications. It can be written as where has been used as a dummy integration variable in order to avoid possible confusion and has been used as the radial variable and indicates the position where a measurement of the wavefield is made. Multiplying both sides by and integrating over gives
Using the orthogonality of the spherical Bessel functions, it follows that (13.2) becomes Recalling that , the inversion formula for the spatial source becomes If the measured wavefield is also separable in space and time so that , then (13.4) can be further simplified to Equation (13.5) gives a simple form for finding the spatial source function at any position once having measured the temporal response of the wavefield at a position . However, since the spherical harmonic transform of the wavefield is required, this implies that sufficient angular information about the wavefield must be obtained in order to permit this spherical harmonic calculation.
13.2. Inversion Formula in 2D
Equation (9.19) admits an inversion leading to the input function. We assume that the source function is separable in time and space (a fairly general assumption) so that . In this case, then and (8.12) becomes where has been used as a dummy integration variable in order to avoid possible confusion and indicates the radial position where the wavefield is measured. Multiplying both sides of (13.6) by and integrating over gives Using the orthogonality of the spherical Bessel functions, it follows that (13.7) becomes Recalling that , the inversion formula for the inhomogeneity becomes If the measured wavefield is also separable in space and time so that , then (13.4) can be further simplified to Equation (13.10) gives a simple form for finding the spatial source function at any position once having measured the wavefield response at a radial position . Comparing (13.10) and (13.5), it is noted that they have identical forms with the exception of the replacement of Hankel and Bessel functions for the 2D case with spherical Hankel and Bessel functions for the 3D case.
13.3. Inversion Formula in 1D
We can also further assume that the wavefield is also separable in time and space so that , which leads to to Using the orthogonality of the Fourier kernel, both sides are multiplied by and integrated over all to give and similarly Several points need to be made regarding (13.12) and (13.13). First, they both give the source function as a function of position along the real line. The interpretation of the variable is that of the position at which the measurement is made and is considered to be a fixed quantity. Both (13.12) and (13.13) are inverse Fourier transforms of , evaluated at or , depending on whether measurements are made in transmission or reflection. Clearly is the time taken for a wave to travel the distance .
The 2D and 3D equivalent to (13.12) and (13.13) are (13.5) and (13.10) which also similarly involve an inverse transform of but have a different kernel for the integration. The reason for the differences in the nature of the kernels from the 2D/3D cases to the 1D case is that in the 2D and 3D cases the kernels used for the Fourier transforms are the Bessel and spherical Bessel functions, respectively, which are the standing wave solutions in those dimensions. However, Green's functions in those dimensions are the Hankel and spherical Hankel functions which are the travelling wave solutions. We note that the kernels for the inversion formulas of (13.5) and (13.10) involve inverse Hankel and spherical Hankel functions, which represent the inverse of Green's function. In contrast, the travelling wave solutions in 1D are the complex exponential and those are also the kernel of 1D Fourier transform. The ratio of the two complex exponentials (one represents the Fourier kernel, the other represents the inverse of 1D Green's function) finally give another complex exponential which represents the shift .
14. Summary and Conclusions
This paper presents the derivation of spatial transfer function relationships for the Helmholtz equation. The focus has been on deriving forward transfer function relationships so that once given a particular input function, the resulting output wavefield can be calculated at any point in space. These are termed “transfer function” relationships because the relationship between the input and output quantities is one of multiplication and no convolution is involved, even in the spatial domain where normally a convolution would be required. Interestingly, instead of Green's function itself, Green's function coefficients for the space are required. These Green function coefficients form the kernel of the output response. Each element of the kernel is then weighted (multiplicatively) by the relevant Fourier/Hankel/Spherical Hankel transform of the input, evaluated on the -dimensional sphere of radius , where is the wavenumber and is the dimension of the space. Combined together, these form the final output response. This is true for dimensions , 2, 3. These simple but powerful results serve to cast the entire problem as an input-output problem with a transfer function relating input to output. In this view of the problem, the input (inhomogeneity) and output (resulting wavefield) in space are related by a simple multiplicative transfer function relationship and not via a convolution. Some applications of these results were shown in the manuscript.