# Time-domain real-valued TM-modal waves in lossy waveguides.

1. INTRODUCTION

A layout of any analytical study of the waveguide time-domain problems consists of two principal parts. The first one applies to derivation of the waveguide modal fields, the amplitudes of which depend on time, t, and an axial waveguide coordinate, z. At this step, the one-dimensional Klein-Gordon Equation (KGE) appears, eventually, which specifies the modal amplitudes. The second part relates to solving the KGE and has directly to do with the physical analysis of the dynamic processes, in which the modal amplitudes participate.

In the first part of our layout, we apply the simple mathematical technique of the vector analysis for a straightforward derivation ([dagger]) of a complete set of the TM-time-domain modal waves from the system of Maxwell's equations with the time derivative, [[partial derivative].sub.t]. Klein-Gordon equation appears in this way naturally enough. Our main goal in the time-domain studies is a fresh look at the dynamic physical properties of the waveguide waves.

All the electromagnetic quantities in Maxwell's equations are measurable. Therefore, the solutions to Maxwell's equations, with their pertinent physical content, should be found in a class of the real-valued functions.

The classical time-harmonic field approach operates in a space of complex-valued solutions. The real-valued fields, which are needed for the physical analysis, can be obtained then via superposition of two solutions as

E (R,t) = E (R, [omega]) exp (-i[omega]t) + [E.sup.*] (R, [omega]) exp (i[omega]t) (1)

where E(R, t) and E(R, [omega]) are the electric field strength and its phasor, respectively; R is a position vector at a point of observation, t an observation time, -[infinity] < t < [infinity], and u a frequency parameter, -[infinity] < [omega] < [infinity].

Just the superposition in (1) presupposes linearity of Maxwell's equations.

We find the solutions to this problem within the framework of a waveguide version of the evolutionary approach to electromagnetics ([double dagger]), see [2-5]. The TM-wave solutions are obtained directly in the Hilbert space, L2, of the real-valued functions in a form, which can be exhibited symbolically as follows:

H(R,t) = I(z, t)H(r) and [epsilon](R,t) = V(z, t)E(r) + e(z, t)Z(r), (2)

where H(r) and E(r) are the two-component basis vectors in the waveguide cross section; Z(r) is a one-component basis vector with the unit vector z, z the axial variable, and r a projection of R onto the waveguide cross section. The scalar factors, I(z, t), V(z, t) and e(z, t), are the amplitudes, physically.

We had an opportunity to discuss a background history of the time-domain electromagnetics in our previous publications, e.g., [7-10]. A reader can find additional information on this topic in other publications, e.g., [11-17].

The article is organized as follows:

In Section 2, the time-domain problem is stated for the transverse-longitudinal decompositions of Maxwell's equations with time derivative, [[partial derivative].sub.t]. Hilbert space, [L.sub.2], of the real-valued functions is defined as the space of solutions.

In Section 3, this problem is solved, and main results are listed. A complete set of the basis elements, [{[H.sub.m], [E.sub.m], [Z.sub.m]}.sup.[infinity].sub.m=1], is derived with needed their physical dimensions. The "magnetic" basis elements, [H.sub.m](r), are obtained with dimension ampere per meter, whereas the "electric" elements, [E.sub.m](r) and [Z.sub.m](r), are derived with dimension volt per meter, both. The modal amplitudes, [I.sub.m](z, t), [V.sub.m](z, t) and [e.sub.m](z, t), which are attached to the basis elements (like in Eq. (2)), have dimensionless quantities. The modal amplitude problem is obtained in a general form. The modal basis problem and modal amplitude problem are obtained as the autonomous ones.

In Section 4, the modal amplitude problem is analyzed. Eventually, this problem comes to solving the Klein-Gordon Equation for the modal amplitude [e.sub.m](z, t). As soon as [e.sub.m](z, t) is obtained, the other amplitudes, [I.sub.m](z, t) and [V.sub.m](z, t), can be found as t- and z-derivatives of [e.sub.m](z, t), respectively. Implementation of this scheme is illustrated by two examples where the variables (z, t) are separated and exact explicit solutions are obtained.

In Section 5, a new version of separation of the variables in the KGE is considered. This version is based on Miller's concept about existence of ten hidden "orbits of symmetry" in the KGE [20]. One of the orbits is taken for analysis herein. A final physical result can be shortly announced as follows. The symmetry of KGE on the chosen orbit discloses existence of a new countable set of the modal amplitudes oscillating with the same cut-off frequency. Besides, these elements have some polynomial factors, which involve z and t variables. The lowest element coincides with that one what yields separation of variables on the orbit (z, t). All other elements belong to another orbit of symmetry.

In Section 6, the time-domain energetic quantities are listed. The conservation of energy law is presented as a continuity equation with z- and t-derivatives. The modal energetic waves, which propagate accompany with the field waves, are discovered and illustrated graphically.

In Section 7, a short summary of new research findings is given.

2. THE NOTATIONS AND THE GOVERNING EQUATIONS

2.1. A Description of the Waveguide and the Notations

The waveguide is geometrically homogeneous along its axis, Oz. Its invariable cross-section domain, S, is bounded by a closed singly-connected contour, L. The shape of L may be rather arbitrary provided that none of the possible its inner angles (i.e., measured within S) exceed [pi]. The standard waveguides with the rectangular cross section satisfy this requirement. A right-handed triplet {z, l, n} (where z x l = n and so on) of the mutually orthogonal unit vectors is introduced. The vector z is oriented along the waveguide axis, Oz; the vector l is tangential to the contour, L; and n is the outer normal to the domain S. The waveguide surface has the properties of the perfect electric conductor. The waveguide is filled with a lossy medium specified by its conductivity, [sigma]. The relative permittivity and permeability of the medium within the waveguide are taken as [epsilon] = [mu] = 1.

2.2. The Governing Equations

The electromagnetic field quantities are the functions of four independent variables, namely: two transverse coordinates, which are accumulated in the vector r, and the variables (z, t). In the form (2) of the expected solution, every field component is presented as a product of two functions. In every product, one function depends on (z, t), only, and the other one depends on r, solely. This trick is known in the partial differential equation theory as an incomplete separation of the variables [18]. In order to apply this method effectively, we should first rearrange Maxwells equations. Projecting these equations onto the waveguide cross section and the waveguide axis results in

[[partial derivative].sub.z]E + [[mu].sub.0][[partial derivative].sub.t] [H x z] = [[nabla].sub.[perpendicular to][].sub.Ez] (3a)

[[nabla].sub.[perpendicular to]] x [z x E] = 0 (3b)

[[partial derivative].sub.z] [z x H] = [[epsilon].sub.0][[partial derivative].sub.t]E + [sigma]E (3c)

[[nabla].sub.[perpendicular to]] x [H x z] = [[epsilon].sub.0][[partial derivative].sub.t][E.sub.z] + [sigma][E.sub.z] (3d)

[[nabla].sub.[perpendicular to]] x H = 0 (3e)

[[nabla].sub.[perpendicular to]] x E + [[partial derivative].sub.z][E.sub.z] = [[sigma]/[[epsilon].sub.0]] [[integral].sup.t.sub.0] ([[nabla].sub.[perpendicular to]] x E + [[partial derivative].sub.z][E.sub.z]) dt' (3f)

where [[nabla].sub.[perpendicular to]] is the transverse part of the operator nabla, [nabla]. This operator, [[nabla].sub.[perpendicular to]], acts on the transverse variables, (r), only. Eqs. (3a) and (3b) are the projections of the vector equation, [nabla] x E = - [[mu].sub.0][[partial derivative].sub.t]H, onto the waveguide cross-section, S, and axis Oz, respectively. Eqs. (3c) and (3d) are the projections of equation [nabla] x H = [[epsilon].sub.0][[partial derivative].sub.t]E + [sigma]E onto the domain S and the axis Oz, respectively. Eq. (3e) is the divergent equation, [[mu].sub.0][nabla] x H = 0. Combination of the equation [[epsilon].sub.0][nabla] x [epsilon] = [rho] and the continuity equation, [nabla] x J = -[[partial derivative].sub.t][rho], for the induced current and charge densities, J = [sigma] [epsilon] and [rho], results in Eq. (3f). A complete set of the Maxwell's differential Eqs. (3a)-(3f) should be supplemented with the boundary conditions as

[E.sub.z] [|.sub.L] = 0, l x E[|.sub.L] = 0, n x H[|.sub.L] = 0. (4)

2.3. The Space of Solutions

A Hilbert space, [L.sub.2], of the real-valued functions is chosen as a space of solutions and defined by an inner product as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

where [X.sub.1] = col([[epsilon].sub.1], [H.sub.1]) and [X.sub.2] = col([[epsilon].sub.2], [H.sub.2]) is a pair of the column-vectors from the space of solutions, [z.sub.1] [less than or equal to] z [less than or equal to] [z.sub.2] and [t.sub.1] [less than or equal to] t [less than or equal to] [t.sub.2]. The free-space constants, [[epsilon].sub.0] and [[mu].sub.0], play role of the weighting coefficients herein. Notice that the procedure of complex conjugation is absent in the integrand in Eq. (5).

3. SOLVING THE PROBLEM (3)-(4)

Equations (3b) and (3e) suggest to search out the vectors E and H as

E (r,z,t) = V (z,t) [[epsilon].sup.1/2.sub.0] [[nabla].sub.[perpendicular to]][phi] (r)

H (r,z,t) = I (z,t) [0 2 [[mu].sup.-1/2.sub.0] [z x [[nabla].sub.[perpendicular to]] [phi] (r)] (6)

where the scalar functions, V(z, t), I(z, t) and [phi] (r), should be found afterwards. The free-space constants, [[epsilon].sub.0] and [[mu].sub.0], are installed in (6) specially in order to provide later the fields E and H with their physical dimensions [Vm.sup.-1] (volt per meter) and [Am.sup.-1] (ampere per meter), respectively. The longitudinal field component, [E.sub.z], can be presented similarly as

[E.sub.z] (r,z,t) = E (z,t) [[epsilon].sup.-1/2/.sub.0] [phi] (r) (7)

where E(z, t) is one more unknown scalar function, as yet.

Substitution of H and [E.sub.z] from Eqs. (6) and (7) to Eq. (3d) results in

I (z, t) [[nabla].sup.2.[perpendicular to]] [phi] (r) = [[c.sup.-1][[partial derivative].sub.t]E (z, t) + 2[gamma][c.sup.-1]E (z, t)] [phi] (r) (8)

where [c.sup.-1] = [square root of ([[epsilon].sub.0][[mu].sub.0])] is the light speed, and 2[gamma] = [sigma]/[[epsilon].sub.0] is a lossy parameter.

Mention here that the boundary condition [E.sub.z][|.sub.L] = 0 results in [phi](r)[|.sub.L] = 0. This fact, jointly with presence of the factor [[nabla].sup.2.sub.[perpendicular to]] [phi] (r) in Eq. (8), suggests to introduce into consideration the well-studied Dirichlet boundary-eigenvalue problem for the transverse Laplacian, [[nabla].sup.2.sub.[perpendicular to]], as

[[nabla].sup.2.sub.[perpendicular to]] [[phi].sub.m] (r) + [[kappa].sup.2.sub.m] [[phi].sub.m] (r) = 0 and [[phi].sub.m] (r)[|.sub.L] = 0 (9)

where [{[[kappa].sup.2.sub.m]}.sup.[infinity].sub.m=1] is a set of the real-valued eigenvalues; [[kappa].sup.2.sub.0] = 0 is also the eigenvalue, and an eigensolution [[phi].sub.0](r), corresponding to [[kappa].sup.2.sub.0], is a solution to the problem as {[[nabla].sup.2.sub.[perpendicular to]][[phi].sub.0] = 0, [[phi].sub.0] [|.sub.L] = 0}. This solution is equal to zero in accordance with the minimum-maximum principle for the harmonic functions ([section]). The subscript m = 1, 2, ... regulates positions of the other eigenvalues, [[kappa].sup.2.sub.m] > 0, on a real axis in order of increasing their numerical values. The set of eigenfunctions, [{[[phi].sub.m]}.sup.[infinity].sub.=1], corresponding to all the eigenvalues, [[kappa].sup.2.sub.m], is complete. This set originates a basis in Hilbert space [L.sub.2](S) provided that the problem (9) is supplemented with the appropriate normalization conditions for the functions [[phi].sub.m](r). Hence, any twice-differentiable function, [phi](r), satisfying the same boundary condition, [phi](r)[|.sub.L] = 0, can be expanded in terms of the basis elements.

Let us take the eigenfunctions [[phi].sub.m](r) as the potential [phi](r) in formulas (6) and (7), and scale by [[kappa].sub.m] the factor E(z,t) as follows: E(z,t) = [[kappa].sub.m][e.sup.m](z,t). Then Eqs. (7) and 8) take the form as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)

where m = 1, 2, ... Eq. (3f) yields a relationship between [V.sub.m] and [e.sub.m] as

[V.sub.m](z,t) - [[kappa].sup.-1.sub.m] [[partial derivative].sub.z] [e.sub.m] (z,t) + 2[gamma] [[integral].sup.t.sub.0] [[V.sub.m] {z,t') - [[kappa].sup.-1.sub.m] [[partial derivative].sub.z] [e.sub.m] {z,t')] dt' = 0. (11)

Factually, this relationship is much simpler. If we introduce a function as

f (z, t) = [[integral].sup.t.sub.0] [[V.sub.m] {z, t') - [[kappa].sup.-1.sub.m] [[partial derivative].sub.z][e.sub.m] (z,t')]dt' (12)

and notice that the first pair in Eq. (11) is [[partial derivative].sub.t]f (z, t), then Eq. (11) yields an initial-value problem for the function f(z, t) as

[[partial derivative].sub.t]f (z,t) + 2[gamma]f (z,t) = 0 and f (z,t) [|.sub.t=0] = 0 (13)

where the initial condition follows from the mean-value theorem applied to the integral in Eq. (12). Solving this problem results in f (z, t) = 0. Hence,

[V.sub.m] (z,t) = [[kappa].sup.-1.sub.m] [[partial derivative].sub.z][e.sub.m] (z,t). (14)

We shall operate henceforward with the field components taken as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

Substitution of the formulas (15) to Eq. (3a) yields

[[partial derivative].sub.z] [V.sub.m] (z, t) + [c.sup.-1] [[partial derivative].sub.t][I.sub.m] (z, t) - [[kappa].sub.m][e.sub.m] (z, t) = 0. (16)

Finally, the formulas (10) and (14), being substituted to Eq. (16), results in the well known Telegraph equation [19], i.e.,

[([[kappa].sub.m]c).sup.-2] [[partial derivative].sup.2.sub.t][e.sub.m] + 2[gamma][([[kappa].sub.m]c).sup.-2] [[partial derivative].sub.t][e.sub.m] - [[kappa].sup.-2.sub.m] [[partial derivative].sup.2.sub.z][e.sub.m] + [e.sub.m] = 0. (17)

Substitution of the fields (15) to Eq. (3c) yields identity, 0 = 0.

3.1. Main Results

3.1.1. Normalization of the Solutions to Dirichlet Problem

The TM-modal wave problem starts with solving the Dirichlet boundary-value problem. Supplement that with the normalization condition, which we choose as

[[nabla].sup.2.sub.[perpendicular to]] [[phi].sub.m] (r) + [[kappa].sup.2.sub.m][[phi].sub.m] (r) = 0, [[phi].sub.m] (r) [|.sub.r[member of]L] = 0, [[N.sup.2.sub.m]/S] [[integral].sub.S] [[phi].sup.2.sub.m][[phi].sup.2.sub.m] (r) ds = N, (18)

where [N.sub.m] is the normalization constant, and N is the "nominated" number 1. That is, we take this number, 1, and assign the force unit N = kg[ms.sup.-2] (newton). Evidently, that the constant, [N.sub.m], has the physical dimension [N.sup.1/2] provided that the potential [[phi].sub.m](r) is a dimensionless quantity, see Example 1.

Example 1 Solve the problem (18) for a standard waveguide with the rectangular cross section specified as 0 [less than or equal to] x [less than or equal to] a and 0 [less than or equal to] y [less than or equal to] b. Separation of the variables in Helmholtz Eq. (18) yields the potential, [[phi].sub.m](r) [equivalent to] [[phi].sub.m](x,y), as

[[phi].sub.m] (x,y) [equivalent to] [[phi].sub.p,q] (x,y) = sin ([pi]px/a) sin ([pi]qy/b), (19)

which corresponds to the eigenvalues distinct from zero,

[[kappa].sup.2.sub.m] [equivalent to] [[kappa].sup.2.sub.p,q] = [[pi].sup.2] [[(p/a).sup.2] + [(q/b).sup.2]]. (20)

The subscript parameter, m, is a doublet, (p,q), composed of the integers, p = 1, 2, ... and q = 1, 2, ... The normalization condition yields [N.sub.m] = 2[N.sup.1/2].

3.1.2. The Waveguide Modal Basis

Complete set of the normalized solutions to the problem (18) generates a complete set of elements of a modal basis as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21)

where [N.sub.m] is proportional to [N.sup.1/2]. Products [[[mu].sup.-1/2.sub.0] [N.sup.1/2]] and [[[epsilon].sup.1/2.sub.0] [N.sup.1/2]] have physical dimensions A (ampere) and V = kg[m.sup.2] [A.sup.-1] [s.sup.-3] (volt), respectively. Square root of the eigenvalues, [square root of [[kappa].sub.2.sub.m]] = [[kappa].sub.m], and the operator [[nabla].sub.[perpendicular to]], have dimension [m.sup.-1] (inverse meter). Thus, the elements [H.sub.m] have dimension [Am.sup.-1] (ampere per meter), and [E.sub.m] and [Z.sub.m], both, have dimension [Vm.sup.-1] (volt per meter). Hence, the modal amplitudes in Eq. (15) are the dimensionless quantities.

3.1.3. The Waveguide Evolutionary Equations

It is convenient to operate with the dimensionless variables, which we introduce by scaling z and t as

[xi] = [[kappa].sub.m]z and [tau] = [[omega].sub.m]t (22)

where [[kappa].sub.m] and [[omega].sub.m] = [[kappa].sub.m]c are the cut-off wave number and the cut-off frequency, physically. Then the modal TM-waves are presentable as

[H.sub.m] = [I.sub.m]([xi],[tau])[H.sup.m](r) and [E.sub.m] = [V.sub.m]([xi],[tau])[E.sub.m](r) + [e.sup.m] ([xi],[tau]) [Z.sub.m](r). (23)

The modal amplitudes, dependent on ([xi], [tau]), are the dimensionless quantities. The amplitude [e.sub.m] ([xi], [tau]) should be found by solving equation

[[partial derivative].sup.2.sub.[tau]] [e.sub.m] ([xi], [tau]) + 2[??][[partial derivative].sub.t][e.sub.m] ([xi], [tau]) - [[partial derivative].sup.2.sub.[xi]] [e.sub.m] ([xi], [tau]) + [e.sup.m] ([xi], [tau]) = 0 (24)

where [??] = [gamma]/[[omega].sub.m] is a dimensionless lossy parameter. The other amplitudes are

[V.sub.m]([xi],[tau]) = [[partial derivative].sub.[xi]][e.sub.m]([xi], [tau]) and [I.sub.m]([xi],[tau]) = -[[partial derivative].sub.[tau]] [e.sup.m]([xi],T) - 2[??][e.sub.m]([xi], [tau]). (25)

3.1.4. A Conclusion

The modal amplitude problem (24), (25) and the modal basis problem (21) are autonomous. Tracing back analysis given in this Section, a reader can be certain that the problems (21) and (24), (25) were extracted from Maxwell's equations with [[partial derivative].sub.t] rigorously, in the mathematical sense.

4. EXAMPLES OF SOLVING KLEIN-GORDON EQUATION

Time derivative of the first order, [[partial derivative].sub.[tau]], can be eliminated from Eq. (24) by applying a substitution for the expected solution, [e.sub.m]([xi], [tau]), as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (26)

where [[??].sub.m] ([xi], [tau]) is a new unknown function. Simple manipulations with Eqs. (26) and (24) result in canonical Klein-Gordon Equation (KGE) as

[[partial derivative].sup.2.sub.[tau]] [[??].sub.m] ([xi], [tau]) - [[partial derivative].sup.2.sub.[xi]] [[??].sub.m] ([xi], [tau]) + [[eta].sup.2][[??].sub.m] = 0 (27)

and slightly changed formulas (25) to the form of

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (28)

where [eta] = [square root of (1 - [[??].sup.2])] [greater than or equal to] 0, [??] = [gamma]/[[omega].sub.m], and [gamma] = [sigma]/(2[[epsilon].sub.0]) is the lossy parameter.

The modal basis problem, (18) and (21), is common for the time-domain fields and also for the time-harmonic fields presentable in the form (23). The modal amplitude problem rests upon solving the KGE (27). This solution yields then the amplitude of the longitudinal field component by formula (26). The amplitudes of the transverse components can be calculated by formulas (28).

As an implementation of this scheme, let us look for the amplitudes of the time-harmonic fields. The method of separation of [xi]- and [tau]-variables can be applied for solving the KGE (27). Exact explicit solutions for the real-valued time-harmonic modal amplitudes are obtained in two examples below.

Example 2 Let an expected solution to Eq. (27) be in the form of

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (29)

where [??] = [omega]/[[omega].sub.m] is a dimensionless frequency, [omega] is a frequency parameter, [[omega].sub.m] = [[kappa].sub.m]c, [[kappa].sub.m] is square root of an eigenvalue from Dirichlet problem (9), c is the light speed. The constants, [A.sub.m] and [[phi].sub.m], are real-valued free numerical parameters ([parallel]). Substitution of the function [[??].sub.m]([xi], [tau]) to Eq. (27) yields [[GAMMA].sub.m] = [+ or -] [square root of ([[??].sup.2] - [[eta].sup.2])] provided that [[??].sup.2] - [[eta].sup.2] [greater than or equal to] 0. In the sign doublet, ([+ or -]), the upper sign corresponds to the wave propagation lengthwise Oz-axis, and the lower sign corresponds to the opposite direction. In order to simplify notations henceforth, we now introduce two "phase functions" as

[[PHI].sub.m] ([xi], [tau]s) = [??][tau] - [xi]][[GAMMA].sub.m] + [[phi].sub.m] [equivalent to] [[omega].sub.t] - (z/c) [square root of ([[omega].sup.2] - [[omega].sup.2.sub.m] + [[gamma].sup.2])] + [[phi].sub.m] (30a)

[[PHI].sub.m] ([xi],[tau]) = [[PHI].sub.m]([xi],[tau]) - [[??].sub.m] (30b)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is a "lossy" phase shift. If [??] = [gamma]/[[omega].sub.m] = 0, (i.e., [gamma] = 0, the waveguide is lossless), then [[??].sub.m] = 0 and [[PHI].sub.m] ([xi],[tau]) = [[PHI].sub.m] ([xi],[tau]).

Let us take the upper sign in ([+ or -] [[GAMMA].sub.m]) and substitute solution (29) to Eqs. (26) and (28). It results in the modal amplitudes of the TM-waves as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (31a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (31b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (31c)

Condition [[GAMMA].sub.m] = 0 yields an equation for the frequency parameter, [[??].sub.cut-off]. Solving this equation specifies so-called "cut-off" frequencies as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (32)

where [??] = [gamma]/[[omega].sub.m] [less than or equal to] 1, [[omega].sub.m] = [[kappa].sub.m]c. Nominally admissible value of the parameter [gamma] = [[kappa].sub.m]c yields [[omega].sub.cut-off] = 0. The values of [[kappa].sub.m], m = 1, 2 ..., are called as the "cut-off wave numbers" of the lossless waveguides.

Notice that all information about the shape and size of the contour, L, of the waveguide cross section, S, is accumulated just in the numbers [[kappa].sub.m]. For the cavity with a rectangular cross section, the values of [[kappa].sub.m] are obtained in Example 1. In a similar way, one can easily find (analytically or numerically) these parameters, [[kappa].sub.m], for the other singly-connected contours.

It should be noticed, as well, that the graphical presentations of the amplitude in terms of [xi]- and [tau]-variables are preferable. All specificity of the waveguide contour L is already put away in Km, see Eq. (22). Thus, the amplitudes for the rectangular and circular waveguides look graphically equally.

Example 3 Let a cavity be a piece of the short-circuited waveguide where 0 [less than or equal to] z [less than or equal to] l. Apply the boundary condition, zxE (r, z, t) = 0, at z = 0 and z = l; r [member of] S. This condition is equivalent to [[partial derivative].sub.[xi]][[??].sub.m][|.sub.z=0,l] = 0. Separation of the variables ([xi], [tau]) by factorization of an expected solution to (27) as [[??].sub.m] = [T.sub.m]([tau])[X.sub.m](l) yields

[T.sup.-1.sub.m] ([tau]) [[d.sup.2]/d[[tau].sup.2]] [T.sub.m] ([tau]) + [[eta].sup.2] = [X.sup.-1.sub.m] ([xi]) [[d.sup.2]/d[[xi].sup.2]] [X.sub.m] ([xi]) = -[[LAMBDA].sup.2.sub.m] (33)

where [[LAMBDA].sup.2.sub.m] > 0 is a real-valued constant of separation of the variables. The first problem, {[X.sub.m]([xi]) : [[d.sup.2]/d[[xi].sup.2]] [X.sub.m] + [[LAMBDA].sup.2.sub.m] [X.sub.m] = 0, [[partial derivative].sub.[xi]] [X.sub.m] [|.sub.z=0,l] = 0}, yields [X.sub.m] ([xi]) = cos([[LAMBDA].sub.m]l[xi]) where [[LAMBDA].sub.m] = [[pi].sub.s]/([[kappa].sub.m]l), s = 1, 2, ..., and [[kappa].sub.m] is the waveguide cut-off wave umber. The other problem, {[T.sub.m]([tau]) : [[d.sup.2]/d[[tau].sup.2]] [T.sub.m] + ([[LAMBDA].sup.2.sub.m] + [[eta].sup.2])[T.sub.m] = 0}, yields [T.sub.m]([tau]) = [B.sub.m] sin([tau] [square root of ([[LAMBDA].sup.2.sub.m] + [[eta].sup.2])] + [[??].sub.m]) provided that [square root of ([[LAMBDA].sup.2.sub.m] + [[eta].sup.2])] [greater than or equal to] 0. The parameters, [B.sub.m] and [[??].sub.m], can be specified by applying the appropriate initial conditions.

Introduce one more pair of the "phase functions" as

[[PSI].sub.ms] ([tau]) = [tau][[OMEGA].sub.ms] + [[??].sub.m] [equivalent to] t[[omega].sub.m][[OMEGA].sub.ms] + [[??].sub.m] (34a)

[[??].sub.ms] ([tau]) = [[PSI].sub.ms] ([xi], [tau]) - [[delta].sub.s] (34b)

where a dimensionless frequency, [[OMEGA].sub.ms], and a "lossy" phase shift, [[delta].sub.s], are

[[OMEGA].sub.ms] = [square root of (1 + [[[pi]s/([[kappa].sub.m]l)].sup.2] - [[??].sup.2])], [[delta].sub.s] = [sin.sup.-1] [[??]/[square root of (1 + [[pi]s/([[kappa].sub.m]l)].sup.2])]. (35)

Then the amplitudes of the TM-cavity modes can be shortly written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (36a)

[V.sup.(c).sub.ms] (z,[tau]) = -[B.sub.m] [[pi]s/([[kappa].sub.m]l)] [e.sup.-[??][tau]] sin [[[PSI].sub.ms] ([tau])] sin ([pi] sz/l) (36b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (36c)

where the superscript, [sup.(c)], implies "cavity". Nominally admissible values of the lossy parameter [??] =[gamma]/[[omega].sub.m] are within the interval 0 [less than or equal to] [??] [less than or equal to] [square root of (1 + [[[pi]s/([[kappa].sub.m]l)].sup.2])].

5. THE KGE AND PARABOLIC CYLINDER FUNCTIONS

Klein-Gordon Eq. (27) is relativistic. Hence, that must maintain its form under action of the relativistic Lorentz transformations in any inertial reference frame. Therefore, the solutions to KGE obey specific properties of symmetry, as well. These symmetries were studied by Miller within the framework of the Group Theory in [20]. Some Miller's ideas were used in our previous publications [5, 9]. In this article, we continue expansion of the group-theoretical results on development of the time-domain electromagnetics.

The Group Theory operates with a so-called point transformation. This implies, technically, that in addition to the independent variables in Eq. (27), i.e., [xi] and [tau], one should introduce a pair of dependent variables as u [equivalent to] u([xi], [tau]) and v [equivalent ti] v([xi], [tau]). Then solution to Eq. (27), [[??].sub.m]([xi], [tau]), can be interpreted as [[??].sub.m] [u([xi], [tau]), v([xi],[tau])]. Formal substitution of the latter to Eq. (27) yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (37)

where the dependent variables are undefined as yet. Right now, one can say only that (u, [u.sub.[xi]], [u.sub.[tau]], [u.sub.[xi][tau]], [u.sub.[xi][xi]], [u.sub.[tau][tau]]) and (v, [v.sub.[xi]], [v.sub.[tau]], [v.sub.[xi][tau]], [v.sub.[xi][xi]], [v.sub.[tau][tau]]) should exist.

Just the Group Theory proposes a way for definition of u([xi], [tau]) and v([xi], [tau]) proceeding from the symmetry of (27). In this article, we consider one of the possible other 10 cases, see Appendix in [5]. Specifically, that looks as

[eta][tau] = ([u.sup.2] + [v.sup.2])/2 and [eta][xi] = uv (38)

where -[infinity] < u < [infinity] and 0 [less than or equal to] v < [infinity]. Inversion of Eq. (38) yields

u ([xi],[tau]) = [square root of ([eta]/2)] ([square root of ([tau] + [xi])] + [square root of ([tau] - [xi])]) and v ([xi],[tau]) = [square root of ([eta]/2)] ([square root of ([tau] + [xi]) - [square root of ([tau] - [xi])]) (39)

where we suppose initially that [tau] [greater than or equal to] 0 and 0 [less than or equal to] [xi] [less than or equal to] [tau].

We can calculate now the variable coefficients, which are placed in the square brackets in Eq. (37). The first pair of the coefficients, which stands in front of the derivatives, [[partial derivative].sup.2.sub.u] and [[partial derivative].sup.2.sub.v], are distinct from zero and, respectively, are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (40)

All the other coefficients are equal to zero. Finally, Eq. (37) looks as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (41)

Virtually, the chosen symmetry transformation of Eq. (27) rearranges that to its equivalent form (41). In other words, (41) holds whenever (27) holds.

It is evident that Eq. (41) can be solved by separation of the variables, u and v. To this aim, the solution should be factorized as follows:

[[??].sub.m] (u,v) = [U.sub.m] (u) [V.sub.m] (v). (42)

Substitution of the product (42) to Eq. (41) yields

[[d.sup.2]/d[u.sup.2]] [U.sub.m] + ([u.sup.2] + [lambda]) [U.sub.m] = 0 and [[d.sup.2]/d[v.sup.2]] [V.sub.m] + ([v.sup.2] + [lambda]) [V.sub.m] = 0 (43)

where [lambda] is a constant (possibly, complex-valued) of separation of the variables.

Both Eq. (43) are pertained to the type of ordinary differential equations for the parabolic cylinder functions. However, it is necessary to rearrange their form to a canonical one (see [21]) which looks as

[[d.sup.2]/[dx.sup.2]] y (x) - ([1/4][x.sup.2] + a) y (x) = 0 (44)

In order to put in order Eq. (43), we slightly change the variables, u and v, and the parameter [lambda], appropriately, what yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (45)

where [alpha] is a new parameter, subscript (m) is omitted. The new variables are

[??] = u[square root of (i2)] = [square root of (i[eta])] ([square root of ([tau] + [xi])] + [square root of ([tau] - [xi])] (46a)

[??] = v[square root of (i2)] = [square root of (i[eta])] ([square root of ([tau] + [xi])] + [square root of ([tau] - [xi])] (46b)

where i = [square root of (-1)] is the imaginary unit.

Both equations in Eq. (45) coincide (with accuracy to notations). Therefore, we can take their solution as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (47)

where U([alpha], *) is one of the possible form of the [solutions.sup.+] to Eq. (44), parameter [alpha] has a half-integer value, namely: [alpha] = -(2n + 1)/2, n = 0, 1, 2, ... .

Finally, we are ready to write down two linearly independent solutions to the equations in (45). They are expressible via products of the parabolic cylinder functions (47). Therefore, they are parametrized by superscript [sup.([alpha])] as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (48a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (48b)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] complex conjugated and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] complex conjugated. These solutions satisfy Eq. (27). So long as the parabolic cylinder functions with

[alpha] = - (2n + 1)/2 where n = 0,1,2, ... (48c)

are convertible to the Hermite polynomials, calculations by formulas (48a)-(48c) result in the simple explicit expressions. In particular, if n = 0 then

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (49)

One can form a linear combination of these solutions as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (50)

where [a.sub.m] and [b.sub.m] are numerical parameters which specify [A.sub.m] and [[phi].sub.m] as

[A.sub.m] = [square root of ([a.sup.2.sub.m] + [b.sup.2.sub.m]) and [[phi].sub.m] = [sin.sup.-1] ([a.sub.m]/[A.sub.m]) = [cos.sup.-1] ([b.sub.m]/[A.sub.m]). (51)

This linear combination one-to-one corresponds to Eq. (31a) provided that the frequency parameter, [??], coincides with the cut-off frequency, n, see Eq. (32).

The modal amplitudes [V.sup.[eta].sub.m] and [I.sup.[eta].sub.m], which are generated by [e.sup.[eta].sub.m], can be found by formulas (25). The results of calculations of [V.sup.[eta].sub.m] and [I.sup.[eta].sub.m] one-to-one correspond to Eqs. (31b), (31c) provided that [??] = [eta]. The lossy phase shift, [[??].sub.m], from Eq. (30b) is found as [[??].sub.m] = [sin.sup.-1] [??] = [cos.sup.-1] [eta] for this case.

As it follows from Eqs. (48a)-(48c), a countable set of other solutions exists for the frequency [??] = [eta]. If n = 1 in (48c), the formulas (48a), (48b) result in

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (52a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (52b)

The case n = 2 yields one more pair of the solutions as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (53a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (53b)

And so on. Thus, in the large, we obtained a new countable set of exact solutions to the relativistic Maxwell's equations with time derivative.

Remark The orbit of symmetry of KGE, where u([xi], [tau]) = [xi] and v([xi], [tau]) = t hold, absorbs the time-harmonic solutions like (31a)-(31c) and (36a)-(36c). In these solutions, the frequency parameter does not suffer restrictions. In the orbit of symmetry, where Eqs. (46a), (46b) hold, the frequency parameter is restricted by condition (32). The solutions, oscillating with the cut-off frequency, appear to be arranged in a new countable set of the fields generated by Eqs. (48a)-(48c).

6. THE INSTANT ENERGETIC CHARACTERISTICS

6.1. The General Energetic Relationships

Introduce a control volume, [delta]V, limited by two consecutive waveguide cross sections located at coordinates z and z + [delta]z. Apply Poynting's theorem to the TM-modal field (23). Take the limit when [delta]z [right arrow] 0. On this way **, two energetic characteristics appear as

[P.sub.mz] ([xi],[tau]) = [[1/S] [[integral].sub.S] z x [E.sub.m] x [H.sub.m] ds] [V.sub.m][I.sub.m] = [c[[kappa].sup.2.sub.m] N] [P.sub.mz] ([xi],[tau]) (54a)

[W.sub.m] ([xi],[tau]) = [[[kappa].sup.2.sub.m]N] [W.sub.m] ([xi],[tau]) (54b)

where factor [c[[kappa].sup.2.sub.m]N] has physical dimension [Wm.sup.-2] (watt per [meter.sup.2]) and factor [[[kappa].sup.2.sub.m]N] has physical dimension [Jm.sup.-3] (joule per [meter.sup.3]). If [delta]V[|.sub.[sigma]z[right arrow]0] [right arrow] 0, dimensional quantity [P.sub.mz] specifies the modal power flow through a fixed waveguide cross section. The dimensional quantity [W.sub.m] specifies the modal field energy density stored at the same cross section. Dimensionless energetic quantities are presentable via the dimensionless modal amplitudes as

[P.sub.mz] ([xi],[tau]) = [I.sub.m] ([xi],[tau]) [V.sub.m] ([xi],[tau]) (55a)

[W.sub.m] ([xi],[tau]) = [W.sup.m.sub.m] ([xi],[tau]) + [W.sup.e.sub.m] ([xi],[tau]) (55b)

[W.sup.e.sub.m] ([xi],[tau]) = [I.sup.2.sub.m] ([xi],[tau])/2 (55c)

[W.sup.2.sub.m] ([xi],[tau]) = [[V.sup.2.sub.m] ([xi],[tau]) + [e.sup.2.sub.m] ([xi],[tau])]/2 (55d)

where [W.sup.m.sub.m] and [W.sup.e.sub.m] are the energy densities stored in the magnetic and electric parts of the modal field, respectively. Mathematically, characteristics (55a)-(55d) specify the global field properties in the space of solutions.

Poynting's theorem yields an energetic continuity equation as

[[partial derivative].sub.[xi]] [P.sub.mz] ([xi],[tau]) + [[partial derivative].sub.[tau]] [W.sub.m] ([xi],[tau]) + 2[??] ([xi],[tau]) + [e.sup.2.sub.m], ([xi],[tau])] = 0. (56)

This is the time-domain law of conservation of the modal field energy, physically. Mathematically, that specifies the local field properties in the space solutions.

Notice in passing that the solutions, which were obtained in Examples 2 and 3, satisfy the energetic continuity Eq. (56). Besides, all the solutions expressible via the parabolic cylinder functions satisfy the conservation of energy law (56), as well, for arbitrary [alpha] = -(2n + 1)/2, n [greater than or equal to] 0.

Combination of Umov theorem [22] and Poynting's theorem [23] yields the definition for instant velocity of transportation of the modal field energy as

[v.sub.m] ([xi],[tau]) = c [[P.sub.mz] ([xi],[tau])]/[[W.sub.m] ([xi],[tau])] (57)

where c is the light speed. One can find discussion on this topic in [10].

6.2. The Time-harmonic Energetic Waveguide Waves

Exact explicit solutions (31a)-(31c) make possible for a fresh look at the energetic field properties in the time-domain. For the time being, we ignore the phase shift, [[??].sub.m], in Eq. (30b). This results in [[??].sub.m] = [[PHI].sub.m] in Eq. (31a). Observation of Eqs. (31b), (31c) under this supposition suggests to introduce a new energetic quantity as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (58)

Observation of Eq. (31a) suggests to introduce one more energetic quantity as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (59)

Notice that the quantity [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] has the exact physical sense provided that [??] = 0 because dm [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. That quantity species difference of the energy density stored in the transverse components of the magnetic and electric fields, see Eqs. (55c), (55d). In the meanwhile, [w.sup.[??].sub.m]([xi],[tau]) is the energy density stored in the longitudinal component of the electric field. So, propagation of the TM-modal electromagnetic wave with its amplitudes (31a)-(31c) is accompanied with propagation of two antiphase energetic waves, (58) and (59). This means, physically, that a periodical exchange by energy occurs between [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Averaging these quantities over the period of oscillations, T = 2[pi]/[omega], yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (60)

So, the averaged values depend neither on time [t.sub.0] nor on coordinate [xi]

Drop that supposition, [[??].sub.m] = [[PHI].sub.m], and for the general case introduce a new energetic quantity as

[S.sup.[??].sub.m] ([xi],[tau]) = (1/2) ([I.sup.[??]2.sub.m] - [V.sup.[??]2.sub.m]) (61)

where [V.sup.[??].sub.m] and [I.sup.[??].sub,.m] are specified in Eqs. (31b), (31c). It seems natural to name this quantity, [S.sup.[??].sub.m] ([xi],[tau]), as a "surplus" of the energy density stored in the transverse field components of a modal wave.

Let us put [A.sub.m] = 1 in formula (59) (for the sake of simplicity) and do the same in Eqs. (31a)-(31c) and (61). One can then study the variations in time dependence of the energetic quantities [S.sup.[??].sub.m] ([xi],[tau]) and [w.sup.[??].sub.m] ([xi],[tau]) in any waveguide cross section by specifying a chosen coordinate [xi].

In Fig. 1, time dependence of the quantities [S.sup.[??].sub.m] (0, [tau]) and [w.sup.[??].sub.m] (0, t) are presented. It is evident that these quantities accomplish antiphased oscillations. Hence, propagation of the TM-modal wave with its amplitudes (31a)-(31c) is accompanied by the energetic wave process where exchange by energy between [S.sup.[??].sub.m] ([xi],[tau]) and [w.sup.[??].sub.m]([xi],[tau]) occurs.

In Fig. 2, time dependence of the power flow, [P.sup.[??].sub.mz] ([xi],[tau]), and the energy density, [W.sup.[??].sub.m] ([xi],[tau]), are presented for the same values [xi] = 0 and [A.sub.m] = 1. Case (a) corresponds to the lossless waveguide, and case (b) corresponds to the value [??] = 0.05 of the dimensionless lossy parameter [??] = [gamma]/[[omega].sub.m], where [gamma] = [sigma]/(2[[epsilon].sup.0]) and [[omega].sub.m] is the cut-off frequency.

Formula for the normalized by c velocity of transportation of the time-harmonic modal field energy was derived by definition (57) with making usage of formulas (55a)-(55d) and (31a)-(31c). The final result is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (62)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is the lossy phase shift.

In Fig. 3, variations of [v.sup.[??].sub.m] (0, [tau])/c are presented. The maximal values of [v.sup.[??].sub.m]/c never exceed 1. The minimal values of [v.sup.[??].sub.m]/c are actually slightly negative if [??] [not equal to] 0. This is caused by presence of the lossy phase shift, [[??].sub.m], in (62).

6.3. The Time-harmonic Energetic Cavity Oscillations

Calculations by formulas (55c), (55d) and (36a)-(36c) of the energy densities stored in the electric and magnetic parts, individually, for the cavity field yield

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (63a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (63b)

Perform averaging Eqs. (63a), (63b) as (1/l) [[integral].sup.l.sub.0] (*) dz. It is evident that

[1/l] [[integral].sup.l.sub.0] [cos.sup.2] ([pi]sz/l) dz = [1/l] [[integral].sup.l.sub.0] [sin.sup.2] ([pi]sz/l) (dz = 1/2, s = 1, 2, ... (64)

Thus, the averaged values, [[bar.W].sup.e.sub.m] and [[bar.W].sup.m.sub.m], depend on time as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (65a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (65b)

where [B.sub.ms] = [1/4] [B.sup.2.sub.m] {1 + [[[pi]s/([[kappa].sup.m]l)].sup.2]}.

If the cavity is lossless (i.e., if [??] = 0), then [[delta].sub.s] = 0, and hence, [[??].sub.ms] = [[PSI].sup.ms]. In this case, Eqs. (65a) and (65b) describe antiphased oscillations of the electric and magnetic field energy in the cavity. This process is analogous to the oscillations of the kinetic and potential energies of a pendulum.

Figure 4 exhibits just these oscillations of the averaged quantities (65a), (65b).

In Fig. 5, the instant velocity of transportation of energy by the cavity field is presented. Calculations were performed by formula (57). The needed amplitudes were taken from Eqs. (36a)-(36c). The waveguide cross section is fixed by specifying [xi] = 1. The normalized by c velocity, [v.sub.m]/c, oscillates with respect to its averaged in time value [[bar.v].sub.m]/c = 0. This implies, physically, that we observe a standing energetic wave in the cavity.

7. SUMMARY OF THE RESEARCH FINDINGS

The system of Maxwell's equations with [[partial derivative].sub.t] is solved explicitly by a simple method via the straightforward calculations. The solution is obtained in Hilbert space [L.sub.2] of the real-valued functions. The time-domain modal fields are found in a form of the transverse-longitudinal decompositions. Each of the field components is a product of two factors. One factor is a vectorial element of the modal basis, dependent only on transverse waveguide coordinates. The basis elements are obtained with the required physical dimensions. The other factors are the appropriate modal amplitudes (dimensionless), each of which dependent solely on time, t, and axial waveguide coordinate, z. The modal basis and modal amplitude problems are autonomous.

Besides the time-domain fields, their instant energetic characteristics are obtained explicitly. The energetic waves, which propagate accompanying the field waves, are established and analyzed. The velocity of transportation of the field energy is obtained explicitly as the function of the variables (t, z).

Existence of a new countable set of the waveguide modes is established. Every modal field from this set oscillates with the same cut-off frequency, but different modes have different amplitudes as the polynomials in t and z.

REFERENCES

[1.] Tretyakov, O. A. and F. Erden, "Evolutionary approach to electromagnetics as an alternative to the time-harmonic field method," 2012 IEEE International Symposium on Antennas and Propagation and USNC-URSI National Radio Science Meeting, Chicago, United States, Jul. 8-14, 2012.

[2.] Tretyakov, O. A., "Evolutionary equations for the theory of waveguides," IEEE AP-S Int. Symp. Dig., 2465-2471, Seattle, Jun. 1994.

[3.] Tretyakov, O. A., "Evolutionary waveguide equations," Soviet Journal on Communication Technology and Electronics (English Translation of Elektrosvyaz i Radiotekhnika), Vol. 35, No. 2, 7-17, 1990.

[4.] Tretyakov, O. A., "Essentials of nonstationary and nonlinear electromagnetic field theory," Analytical and Numerical Methods in Electromagnetic Wave Theory, M. Hashimoto, M. Idemen, O. A. Tretyakov, Eds., Chapter 3, Science House Co. Ltd., Tokyo, 1993.

[5.] Aksoy, S. and O. A. Tretyakov, "Evolution equations for analytical study of digital signals in waveguides," Journal of Electromagnetic Waves and Applications, Vol. 17, No. 12, 1665-1682, 2003.

[6.] http://www.springer.com/birkhauser/mathematics/journal/28.

[7.] Tretyakov, O. A. and F. Erden, "Separation of the instantaneous and dynamic polarizations in studies of dispersive dielectrics," MSMW'07 Symposium Proceedings, 42-48, Kharkov, Ukraine, Jun. 25-30, 2007.

[8.] Erden, F. and O. A. Tretyakov, "Excitation by a transient signal of the real-valued electromagnetic fields in a cavity," Phys. Rev., E, Vol. 77, 056605, May 2008.

[9.] Tretyakov, O. A. and O. Akgun, "Derivation of Klein-Gordon equation from Maxwell's Equations and study of relativistic time-domain waveguide modes," Progress In Electromagnetics Research, Vol. 105, 171-191, 2010.

[10.] Tretyakov, O. A. and M. Kaya, "The real-valued time-domain TE-modes in lossy waveguides," Progress In Electromagnetics Research, Vol. 127, 405-426, 2012.

[11.] Gabriel, G. J., "Theory of electromagnetic transmission structures, Part I: Relativistic foundation and network formalism," Proc. IEEE, Vol. 68, No. 3, 354-366, 1980.

[12.] Borisov, V. V., Transient Electromagnetic Waves, Leningrad Univ. Press, 1987 (in Russian).

[13.] Kristensson, G., "Transient electromagnetic wave propagation in waveguides," Journal of Electromagnetic Waves and Applications, Vol. 9, Nos. 5-6, 645-671, Sept. 1995.

[14.] Shvartsburg, A. B., "Single-cycle waveforms and non-periodic waves in dispersive media (exactly solvable models)," Phys. Usp., Vol. 41, No. 1, 77-94, Jan. 1998.

[15.] Slivinski, A. and E. Heyman, "Time-domain near-field analysis of short-pulse antennas--Part I: Spherical wave (multipole) expansion," IEEE Trans. on Antenn. and Propag., Vol. 47, 271-279, Feb. 1999.

[16.] Geyi, W., "A time-domain theory of waveguides," Progress In Electromagnetics Research, Vol. 59, 267-297, 2006.

[17.] Dusseaux, R., "Telegraphist's equations for rectangular waveguides and analysis in nonorthogonal coordinates," Progress In Electromagnetics Research, Vol. 88, 53-71, 2008.

[18.] Polyanin, A. D. and A. V. Manzhirov, Handbook of Mathematics for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, FL, 2006.

[19.] Polyanin, A. D., Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, FL, 2002.

[20.] Miller, Jr., W., Symmetry and Separation of Variables, Addison Wesley Publication Co., Boston, MA, 1977.

[21.] Abramowitz, M. and I. A. Stegun, Handbook of Mathematical Functions, Chapter 19, Dover Publications, Inc., 1965.

[22.] Umov, N. A., "Ein theorem iiber die wechselwirkungen in endlichen entfernungen," Zeitschrift Fur Mathematik Und Physik, 97, XIX, 1874.

[23.] Poynting, J. H., "On the transfer of energy in the electromagnetic field," Philos. Trans. of the Royal Society of London, Vol. 175, 343-361, 1884.

([double dagger]) Postulation of the harmonic in time varying the fields (like in Eq. (1)) "kills" the time derivative, [[partial derivative].sub.t], in Maxwell's equations at the very beginning. We keep [[partial derivative].sub.t] from the beginning and up to the end in our analysis. In this way, we obtain the KGE where [[partial derivative].sup.2.sub.t] participates. Mathematicians call all the differential equations with time derivatives as evolution (or evolutionary) equations, see, for example, Journal of Evolution Equations [6]. Apparently, this term was induced by Cauchy theorem for dynamic systems. Any solution, obtained via the Cauchy theorem, exhibits how a process progresses in time (i.e., evolves, shortly) starting from a given initial state and up to the state at a time of observation. That is why we name this method as Evolutionary Approach to Electromagnetics (EAE).

([section]) It is so because the cross section S is chosen as a singly-connected domain.

([parallel]) Involvement of two parameters, Am and [[phi].sub.m], in the solution (29) is equivalent to presentation that via two linearly independent sine and cosine functions as [[??].sub.m([xi], [tau]) = [a.sub.m] sin(*) + [b.sub.m cos(*) where (*) is ([tau][omega]/[[omega].sub.m] - [xi][[GAMMA].sub.m]), [a.sub.m] = [A.sub.m] cos [[phi].sub.m], [b.sub.m] = [A.sub.m] sin [[phi].sub.m].

([paragraph]) Miller has established ten of so-called "orbits of symmetry." These result in definition for the eleven pairs of the dependent variables, {u([xi], [tau]), v([xi], [tau])}. The first pair as u([xi], [tau]) = [xi] and v(c, [tau]) = [tau] is trivial. This case directs electromagnetics to the time-harmonic field approach. The other ten pairs are capable of elucidation of new aspects in the time-domain electromagnetics. One can find a complete list of these pairs in Appendix [5].

+ The solutions to Eq. (44) can be expressed via confluent hypergeometric series (Kummer's and Whittaker's functions). We use standard solutions denoted in [21] as U([alpha], x).

* One can find details of similar procedure in article [10].

Oleg A. Tretyakov * and Mehmet Kaya

Department of Electronics Engineering, Gebze Institute of Technology, Gebze, Kocaeli, Turkey

* Corresponding author: Oleg A. Tretyakov (o.tretyakov@gmail.com).

Received 2 March 2013, Accepted 11 April 2013, Scheduled 17 April 2013