spacetime triangle diagram technique

1. (Electromagnetics) The system of Maxwell's equations is reduced to a second-order PDE for the field components, or potentials, or their derivatives.
2. The spatial variables are separated using convenient expansions into series and/or integral transforms—except one that remains bounded with the time variable, resulting in a PDE of hyperbolic type.
3. The resulting hyperbolic PDE and the simultaneously transformed initial conditions compose a problem, which is solved using the [http://mathworld.wolfram.com/Riemann-VolterraMethod.html Riemann–Volterra integral formula]. This yields the generic solution expressed via a double integral over a triangle domain in the bounded-coordinate—time space. Then this domain is replaced by a more complicated but smaller one, in which the integrant is essentially nonzero, found using a strictly formalized procedure involving specific spacetime triangle diagrams (see, e.g., Refs.

A.B. Utkin,

Localized Waves Emanated by Pulsed Sources: The Riemann–Volterra Approach.

In: Hugo E. Hernández-Figueroa, Erasmo Recami, and Michel Zamboni-Rached (eds.)

[https://books.google.com/books?isbn=3527671536 Non-diffracting Waves.]

Wiley-VCH: Berlin, {{ISBN|978-3-527-41195-5}},

pp. 287–306 (2013)

A.B. Utkin,

[http://www.sciencedirect.com/science/article/pii/S0165212511001430 The Riemann–Volterra time-domain technique for waveguides: A case study for elliptic geometry.] Wave Motion 49(2), 347–363 (2012), doi: 10.1016/j.wavemoti.2011.12.001

V.V. Borisov, A.V. Manankova, A.B. Utkin,

[http://iopscience.iop.org/0305-4470/29/15/020 Spherical harmonic representation of the electromagnetic field produced by a moving pulse of current density],

Journal of Physics A: Mathematical and General 29(15), 4493–4514 (1996), doi: 10.1088/0305-4470/29/15/020

).

1. In the majority of cases the obtained solutions, being multiplied by known functions of the previously separated variables, result in the expressions of a clear physical meaning (nonsteady-state modes). In many cases, however, more explicit solutions can be found summing up the expansions or doing the inverse integral transform.

The STTD technique belongs to the second among the two principal ansätze for theoretical treatment of waves — the frequency domain and the direct spacetime domain.

The most well-established method for the inhomogeneous (source-related) descriptive equations of wave motion is one based on the Green's function technique.J. D. Jackson,

[https://www.amazon.com/Classical-Electrodynamics-Third-Edition-Jackson/dp/047130932X Classical Electrodynamics],

3rd ed., Wiley, New York (1999) For the circumstances described in Section 6.4 and Chapter 14 of Jackson's Classical Electrodynamics, it can be reduced to calculation of the wave field via retarded potentials (in particular, the Liénard–Wiechert potentials).

Despite certain similarity between Green's and Riemann–Volterra methods (in some literature the Riemann function is called the Riemann–Green function

see, e.g., G. A. Korn and T. M. Korn, Mathematical Handbook for Scientists and Engineers, Courier Dover Publications, New York (2000)

), their application to the problems of wave motion results in distinct situations:

• The definitions of both Green's function and corresponding Green's solution are not unique as they leave room for addition of arbitrary solution of the homogeneous equation; in some circumstances the particular choice of Green's function and the final solution are defined by boundary condition(s) or plausibility and physical admissibility of the constructed wavefunctions.

A comprehensive discussion of this subject can found in H. Kleinert,

[http://www.worldscientific.com/worldscibooks/10.1142/7305 Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets], 5th ed., World Scientific, Singapore (2009)

The Riemann function is a solution of the homogeneous equation that additionally must take a certain value at the characteristics and thus is defined in a unique way.

• In contrast to Green's method that provides a particular solution of the inhomogeneous equation, the Riemann–Volterra method is related to the corresponding problem, comprising the PDE and initial conditions,

R. Courant and D. Hilbert, [https://books.google.com/books?id=SAdRAAAAMAAJ Methods of Mathematical Physics], Vol. 2, Wiley, New York (1989)

V.I. Smirnov, [https://archive.org/details/courseofhigherma0005smir A Course of Higher Mathematics, Vol. 4: Integral Equations and Partial Differential Equations], Pergamon Press, Oxford (1964)

and it was the Riemann–Volterra representation that Smirnov used in his Course of Higher Mathematics to prove the uniqueness of the solution to the above problem (see, item 143).

• In the general case, Green's formula implies integration over the entire domain of variation of coordinates and time, while integration in the Riemann–Volterra solution is carried out within a limited triangle region, assuring the boundness of the solution support.
• Causality of the (unique) Riemann–Volterra solution is provided automatically, without need to recur to additional considerations, such as the retarded nature of the argument, wave propagation in certain direction, specific choice of the integration path, etc. (Usually the descriptive equations, such as the classical scalar wave equation, possess the T-symmetry. It is the time-asymmetric initial conditions that define the arrow of time through the limitation of the integration domain in the Riemann formula to $t\; >\; 0$, see more in and a particular example given below.)
• Green's function can be readily derived from the Liénard–Wiechert potential of a moving point source, but concrete calculation of the wavefunction, inevitably involving the analysis of the retarded argument, may develop in a rather complicated task unless some special techniques, like the parametric method,

C.J. Chapman,

[http://rspa.royalsocietypublishing.org/content/431/1881/157 The spiral Green function in acoustics and electromagnetism],

Proc. Roy. Soc. A 431(1881), 157–167 (1990), doi: 10.1098/rspa.1990.0124

are invoked. The Riemann-Volterra approach presents the same or even more serious difficulties, especially when one deals with the bounded-support sources: here the actual limits of integration must be defined from the system of inequalities involving the space-time variables and parameters of the source term. However, this definition can be strictly formalized using the spacetime triangle diagrams. Playing the same role as the Feynman diagrams in particle physics, STTDs provide a strict and illustrative procedure for definition of areas with the same analytic representation of the integration domain in the 2D space spanned by the non-separated spatial variable and time.

• The method can only be applied to problems possessing known Riemann function.
• Application of the method and analysis of the results obtained require more profound knowledge of the special functions of mathematical physics (e.g., operating with the generalized functions, Mathieu functions of different kinds and [http://mathworld.wolfram.com/LommelFunction.html Lommel's functions of two variables]) than Green's function method.
• In some cases the final integrals require special consideration in the domains of rapid oscillation of the Riemann function.

Several efficient methods for scalarizing electromagnetic problems in the orthogonal coordinates $x\_1,\; x\_2,\; x\_3$ were discussed by Borisov in Ref.

V.V. Borisov, Electromagnetic Fields of Transient Currents. Leningrad State University Press: Leningrad (1996, in Russian)

The most important conditions of their applicability are $h\_3\; =\; 1$ and

$\backslash partial\_3(h\_1/h\_2)\; =\; 0$, where $h\_i,\; i=\; 1,2,3$ are the

metric (Lamé) coefficients (so that the squared length element is $ds^2\; =\; h\_1^2\; dx\_1^2+h\_2^2\; dx\_2^2+h\_3^2\; dx\_3^2$). Remarkably, this condition is met for the majority of practically important coordinate systems, including the Cartesian, general-type cylindrical and spherical ones.

For the problems of wave motion is free space, the basic method of separating spatial variables is the application of integral transforms, while for the problems of wave generation and propagation in the guiding systems the variables are usually separated using expansions in terms of the basic functions (modes) meeting the required boundary conditions at the surface of the guiding system.

In the Cartesian $\backslash left\backslash \{\; x\_1=x,\; \backslash ;\; x\_2=y,\; \backslash ;\; x\_3=z\; \backslash right\backslash \}$ and general-type cylindrical coordinates

$\backslash left\backslash \{\; x\_1=x\_1(x,y),\; \backslash ;\; x\_2=x\_2(x,y),\; \backslash ;\; x\_3=z\; \backslash right\backslash \}$ separation of the spatial variables result in the initial value problem for a hyperbolic PDE known as the 1D Klein–Gordon equation (KGE)

: $$

\begin{align}

& \left( \partial_\tau ^2 - \partial_z^2 + k^2 \right)\psi(\tau,z) = f(\tau,z) \\

& \psi (\tau ,z) = 0 \text{ for } \tau < 0

\end{align}

Here $\backslash tau$ is the time variable expressed in units of length using some characteristic velocity (e.g., speed of light or sound),

$k$ is a constant originated from the separation of variables, and $f(\backslash tau\; ,z)$ represents a part of the source

term in the initial wave equation that remains after application of the variable-separation procedures (a series coefficient or a result of

an integral transform).

The above problem possesses known Riemann function

: $$

R_k(\tau ,z;\tau ',z') = J_0\left( k\sqrt {(\tau - \tau')^2 - (z - z')^2} \right),

where $J\_0(\backslash cdot)$ is the Bessel function of the first kind of order zero.

{{multiple image

| align = right

| direction = horizontal

| footer = The simplest STTD representing a triangle integration domain resulted from the [http://mathworld.wolfram.com/Riemann-VolterraMethod.html Riemann–Volterra integral formula].

| footer_align = center

| image1 = F140910a_STTDmethod_F_RotatedDomain.gif

| alt1 = Canonical variables

| caption1 = Canonical variables ξ, η.

| width1 = 300

| image2 = F140910a_STTDmethod_F_zTauDomain.gif

| alt2 = Initial variables

| caption2 = Initial variables z, τ.

| width2 = 300

}}

Passing to the canonical variables $z,\; \backslash tau\; \backslash to\; \backslash xi,\; \backslash eta$ one gets the simplest STTD diagram reflecting straightforward application of the Riemann–Volterra method, with the fundamental integration domain represented by spacetime triangle MPQ (in dark grey).

Rotation of the STTD 45° counter clockwise yields more common form of the STTD in

the conventional spacetime $z,\; \backslash tau$.

For the homogeneous initial conditions the (unique) solution of the problem is given by the Riemann formula

: $$

\psi(\tau ,z) = \frac{1}{2} \iint\limits_{\triangle MPQ} \, d\tau' \, dz'R(\tau ,z;\tau',z')f(\tau',z').

Evolution of the wave process can be traced using a fixed observation point ($z\; =\; \backslash text\{const\}$) successively increasing the triangle height ($\backslash tau$) or, alternatively, taking "momentary picture" of the wavefunction $\backslash psi$ by shifting the spacetime triangle along the $z$ axis ($\backslash tau\; =\; \backslash text\{const\}$).

More useful and sophisticated STTDs correspond to pulsed sources whose support is limited in spacetime. Each limitation produce specific modifications in the STTD, resulting to smaller and more complicated integration domains in which the

integrand is essentially non-zero. Examples of most common modifications and their combined actions are illustrated below.

{{multiple image

| align = left

| direction = horizontal

| width = 300

| image1 = F140910a_STTDmethod_F_zTau_z0.gif

| width1 =

| alt1 =

| caption1 = STTD for a source limited from left by plane $z\; =\; 0$, i.e. , which is the case, e.g., for a travelling source propagating along a semi-infinite radiator $l$.

| image2 = F140910a_STTDmethod_F_zTau_l.gif

| width2 =

| alt2 =

| caption2 = STTD for a source limited from right by plane $z\; =\; l$, i.e. $f(\backslash tau,z)\; =\; 0\; \backslash text\{\; for\; \}\; z\; >\; l.$

| image3 = F140911a_STTDmethod_F_zTau_z0l.gif

| width3 =

| alt3 =

| caption3 = STTD for a source limited from both sides, i.e. $f(\backslash tau\; ,z)\; =\; 0\; \backslash text\{\; for\; \}\; z\; \backslash notin\; [0,l]$, which is the case, e.g., for a travelling source propagating along a radiator of finite length $l$.

| header = Static limitations to the source area

| header_align =

| header_background =

| footer =

| footer_align =

| footer_background =

| background color =

}}

{{multiple image

| align = left

| direction = horizontal

| width = 300

| image1 = F140911a_STTDmethod_F_zTau_front.gif

| width1 =

| alt1 =

| caption1 = STTD for a semi-infinite travelling source pulse.

| image2 = F140911a_STTDmethod_F_zTau_frntBack.gif

| width2 =

| alt2 =

| caption2 = STTD for a finite travelling source pulse.

| image3 = F140911a_STTDmethod_F_zTau_frntBack_z0.gif

| width3 =

| alt3 =

| caption3 = STTD for a finite travelling source pulse propagating along a semi-infinite radiator $z\; \backslash isin\; \backslash left[\; 0,\; \backslash infty\; \backslash right)$.

| header = Combined action of limitations of different type, see Refs.

V.V. Borisov and A.B. Utkin,

[http://iopscience.iop.org/0022-3727/28/4/003 The transient electromagnetic field produced by a moving pulse of line current], Journal of Physics D: Applied Physics 28(4), 614-622 (1995), doi: 10.1088/0022-3727/28/4/003

A.B. Utkin,

[http://www.opticsinfobase.org/josaa/abstract.cfm?uri=josaa-29-4-457 Droplet-shaped waves: causal finite-support analogs of X-shaped waves],

J. Opt. Soc. Am. A 29(4), 457-462 (2012), doi: 10.1364/JOSAA.29.000457

A.B. Utkin,

[https://ieeexplore.ieee.org/xpl/articleDetails.jsp?arnumber=6712820 Droplet-shape wave produced by line macroscopic current pulse of finite length],

IEEE Xplore DD-2013, {{ISBN|978-1-4799-1037-3}}, 145–150 (2013), doi: 10.1109/DD.2013.6712820

for details and more complicated examples

| header_align =

| header_background =

| footer =

| footer_align =

| footer_background =

| background color =

}}

{{multiple image

| align = left

| direction = vertical

| width = 910

| image1 = F140913a_STTDmethod_F_zTau_z0lT_S_all.gif

| width1 =

| alt1 =

| caption1 = A sequence of generic STTDs for a "short", finite source pulse of duration $T$ propagating along a finite radiator $z\; \backslash isin\; [0,l]$ with a constant velocity $\backslash beta$.{{Citation needed|date=December 2019|reason=removed citation to predatory publisher content}} In this case the source can be expressed in the form

: $f(\backslash tau,z)\; =\; f\_0(\backslash tau,z)\; \backslash times$

h(z) h(l - z)

h \left( \tau - \frac{z}{\beta} \right)

h \left( T - \tau + \frac{z}{\beta} \right),

    where $h\; (\; \backslash cdot\; )$ is the Heaviside step function.

| image2 = F140913a_STTDmethod_F_zTau_z0lT_L_all.gif

| width2 =

| alt2 =

| caption2 = The same STTD sequence for a "long" pulse.{{Citation needed|date=December 2019|reason=removed citation to predatory publisher content}}

| image3 =

| width3 =

| alt3 =

| caption3 =

| header =

| header_align =

| header_background =

| footer =

| footer_align =

| footer_background =

| background color =

}}

{{clear}}

In the spherical coordinate system — which in view of the General considerations must be represented in the sequence

$\backslash left\backslash \{x\_1=\backslash theta,\; \backslash ;\; x\_2=\backslash varphi,\; \backslash ;\; x\_3=r\; \backslash right\backslash \}$, assuring $h\_3\; =\; 1$ — one can scalarize problems for the transverse electric (TE) or transverse magnetic (TM) waves using the Borgnis functions, Debye potentials or Hertz vectors. Subsequent separation of the angular variables $\backslash theta,\; \backslash varphi$ via expansion of the initial wavefunction $\backslash psi\; \backslash left(\; \backslash tau,\; \backslash theta,\; \backslash varphi,\; r\; \backslash right)$ and the source

: $\backslash ;\; f\backslash left(\; \backslash tau,\; \backslash theta,\; \backslash varphi,\; r\; \backslash right)$ in terms of

: $\backslash left($

\begin{array}{c}

\sin m\varphi \\

\cos m\varphi

\end{array}

\right) P_n^{(m)}(\cos \theta),

where $P\_n^\{(m)\}\; \backslash !\; \backslash left(\; \backslash cdot\; \backslash right)$ is the associated Legendre polynomial of degree $n$ and order

$m$, results in the initial value problem for the hyperbolic

Euler–Poisson–Darboux equation

: $$

\begin{align}

& \left( \partial _\tau ^2 - \partial _r^2 + \frac{n(n + 1)}{r^2} \right)\psi_{nm}(\tau,r)

= f_{nm}(\tau, r) \\

& \psi _{nm}(\tau,r) = 0 \text{ for } \tau < 0

\end{align}

known to have the Riemann function

: $$

R(\tau,r;\tau',r') = {P_n}\left( \frac{r^2 + {r'}^2 - (\tau - \tau')^2}{2rr'} \right),

where $P\_n\; \backslash !\; \backslash left(\; \backslash cdot\; \backslash right)$ is the (ordinary) Legendre polynomial of degree $n$.

The STTD technique represents an alternative to the classical Green's function method. Due to uniqueness of the solution to the initial value problem in question, in the particular case of zero initial conditions the Riemann solution provided by the STTD technique must coincide with the convolution of the causal Green's function and the source term.

The two methods provide apparently different descriptions of the wavefunction: e.g., the Riemann function to the Klein–Gordon problem is a Bessel function (which must be integrated, together with the source term, over the restricted area represented by the fundamental triangle MPQ) while the retarded Green's function to the Klein–Gordon equation is a Fourier transform of the imaginary exponential term (to be integrated over the entire plane $z\text{'},\backslash tau\text{'}$, see, for example, Sec. 3.1. of Ref.W. Geyi,

[http://www.jpier.org/PIER/pier.php?paper=0510212 A time-domain theory of waveguide],

Progress in Electromagnetics Research 59, 267–297 (2006), doi: 10.2528/PIER05102102

) reducible to

: $$

G_k(\tau,z;\tau',z') =

- \frac{1}{(2\pi)^2}

\int\limits_{-\infty}^\infty dp \, {\rm e}^{ip(z - z')}

\int\limits_{-\infty}^\infty d\Omega \, \frac{{\rm e}^{i\Omega(\tau -\tau')}}{\Omega^2 - p^2 - k^2}.

Extending integration with respect to $\backslash Omega$ to the complex domain, using the residue theorem

(with the poles $\backslash Omega\_\{1,2\}$ chosen as $\backslash lim\_\{\backslash varepsilon\; \backslash to\; 0+\}\; \backslash left(\backslash pm\; \backslash sqrt\; \{p^2\; +\; k^2\; +\; i\backslash varepsilon\}\; \backslash right)$ to satisfy the causality conditions) one gets

: $$

G_k(\tau,z;\tau',z') =

\frac{1}{\pi}

\int\limits_0^\infty dp \, \frac{{\sin \left((\tau - \tau')\sqrt {p^2 + k^2} \right)}}{\sqrt {p^2 + k^2}}

\cos(p(z - z')).

Using formula 3.876-1 of Gradshteyn and Ryzhik,{{cite book |author-first1=Izrail Solomonovich |author-last1=Gradshteyn |author-link1=Izrail Solomonovich Gradshteyn |author-first2=Iosif Moiseevich |author-last2=Ryzhik |author-link2=Iosif Moiseevich Ryzhik |author-first3=Yuri Veniaminovich |author-last3=Geronimus |author-link3=Yuri Veniaminovich Geronimus |author-first4=Michail Yulyevich |author-last4=Tseytlin |author-link4=Michail Yulyevich Tseytlin |author-first5=Alan |author-last5=Jeffrey |editor-first1=Daniel |editor-last1=Zwillinger |editor-first2=Victor Hugo |editor-last2=Moll |editor-link2=Victor Hugo Moll |translator=Scripta Technica, Inc. |title=Table of Integrals, Series, and Products |publisher=Academic Press, Inc. |date=2015 |orig-year=October 2014 |edition=8th |language=English |isbn=978-0-12-384933-5 |lccn=2014010276 |title-link=Gradshteyn and Ryzhik |chapter=3.876. |page=486}}

: $$

\int\limits_0^\infty dx \, \frac{\sin(p\sqrt {x^2 + a^2})}{\sqrt {x^2 + a^2}} \cos(bx) =

\left\{ \begin{align}

& \frac{\pi}{2} J_0(a\sqrt {p^2 - b^2}) & \text{ for } & 0 < b < p \\

& 0 & \text{ for } & b > p > 0

\end{align}

\right.

\qquad (a>0),

the last Green's function representation reduces to the expressionApparently this result was first published by Geyi (2006: 275), merely as a way to simplify the Green's solution and reduce the domain of integration.

: $$

\frac{1}{2}

J_0\left( k\sqrt {(\tau - \tau ')^2 - (z - z')^2} \right)

h(\tau - \tau ' - |z - z'|),

in which 1/2 is the scaling factor of the Riemann formula and $J\_0(\backslash cdot)$ the Riemann function, while the Heaviside step function $h(\backslash cdot)$ reduces, for $\backslash tau\; >\; 0$, the area of integration to the fundamental triangle MPQ, making the Green's function solution equal to that provided by the STTD technique.

{{clear}}

• V.V. Borisov, N.M. Reutova, A.B. Utkin, [http://iopscience.iop.org/0305-4470/38/10/012 Electromagnetic waves produced by a travelling current pulse with high-frequency filling.] Journal of Physics A: Mathematical and General, 38(10), 2225–2240 (2005), doi: 10.1088/0305-4470/38/10/012
• V.V. Borisov, Nonsteady-State Electromagnetic Waves. Leningrad: Leningrad State University Press: Leningrad (1987, in Russian)

Category:Wave mechanics

Category:Electrodynamics