Accessibility / Report Error

Analysis of Electromagnetic Problems in the Presence of Non-uniform Movements

Abstract

Electromagnetic problems with moving structures are considered. For objects moving at a uniform velocity, the most popular technique used in the literature consists of a change of the reference frame and the utilization of Lorentz transformation of Maxwell’s electric and magnetic fields and Voigt-Lorentz transformations of space and time variables. In order to account for non-uniform motions, a numerical approach based on the finitedifference time-domain (FDTD) method is used. In this approach, Maxwell’s equations are applied without modification and the motion of objects is implemented by changing their positions in the FDTD time loop. Three types of movements are considered: vibration, rotation, and acceleration. Full-wave numerical results are reported for different objects in motion: observer, source, and reflecting electromagnetic surface. In all the problems analyzed, interesting results are obtained and these results are discussed in detail.

Index Terms
Acceleration; FDTD method; Rotation; Vibration.

I. INTRODUCTION

The problems of electromagnetic wave propagation and radiation in the presence of moving objects have attracted a lot of attention for a long time, due to their wide application in many areas such as RF Doppler radars, astrophysics, global positioning system (GPS), and optical gyroscopes. Numerous investigations have been carried out in these domains, which are interesting from a theoretical and practical point of view [1[1] M. Marvasti and H. Boutayeb, “FDTD analysis of the Sagnac effect employed in the global positioning system,” IEEE Transactions on Antennas and Propagation, vol. 71, pp. 9119-9123, 2023. DOI: 10.1109/TAP.2023.3296903.
https://doi.org/10.1109/TAP.2023.3296903...
][2] J. Cooper, “Scattering of electromagnetic fields by a moving boundary: the one-dimensional case,” IEEE Transactions on Antennas and Propagation, vol. 28, no. 6, pp. 791-795, 1980. [3] B. Michielsen, G. Herman, A. De Hoop, and D. De Zutter, “Three-dimensional relativistic scattering of electromagnetic waves by an object in uniform translational motion,” Journal of Mathematical Physics, vol. 22, no. 11, pp. 2716-2722, 1981. [4] N. Engheta, M. W. Kowarz, and D. L. Jaggard, “Effect of chirality on the Doppler shift and aberration of light waves,” Journal of Applied Physics, vol. 66, no. 6, pp. 2274-2277, 1989. [5] C. Neipp, A. Hernández, J. Rodes, A. Márquez, T. Beléndez, and A. Beléndez, “An analysis of the classical Doppler effect,” European Journal of Physics, vol. 24, no. 5, p. 497, 2003. [6] S. Sahrani and M. Kuroda, “Numerical analysis of the electromagnetic wave scattering from a moving dielectric body by overset grid generation method,” in 2012 IEEE Asia-Pacific Conference on Applied Electromagnetics (APACE), pp. 6-10, 2012. [7] J. Cooper, “Long-time behavior and energy growth for electromagnetic waves reflected by a moving boundary,” IEEE Transactions on Antennas and Propagation, vol. 41, no. 10, pp. 1365-1370, 1993. [8] S. Borkar and R. Yang, “Scattering of electromagnetic waves from rough oscillating surfaces using spectral Fourier method,” IEEE Transactions on Antennas and Propagation, vol. 21, no. 5, pp. 734-736, 1973. [9] R. Kleinman and R. Mack, “Scattering by linearly vibrating objects,” IEEE Transactions on Antennas and Propagation, vol. 27, no. 3, pp. 344-352, 1979. [10] D. De Zutter, “Reflections from linearly vibrating objects: plane mirror at oblique incidence,” IEEE Transactions on Antennas and Propagation, vol. 30, no. 5, pp. 898-903, 1982. [11] M. Ho, “Numerical simulation of scattering of electromagnetic waves from traveling and/or vibrating perfect conducting planes,” IEEE Transactions on Antennas and Propagation, vol. 54, no. 1, pp. 152-156, 2006. [12] V. C. Chen, F. Li, S.-S. Ho, and H. Wechsler, “Micro-Doppler effect in radar: phenomenon, model, and simulation study,” IEEE Transactions on Aerospace and Electronic Systems, vol. 42, no. 1, pp. 2-21, 2006. [13] H. Gao, L. Xie, S. Wen, and Y. Kuang, “Micro-Doppler signature extraction from ballistic target with micro-motions,” IEEE Transactions on Aerospace and Electronic Systems, vol. 46, no. 4, pp. 1969-1982, 2010. [14] N. Engheta, A. Mickelson, and C. Papas, “On the near-zone inverse Doppler effect,” IEEE Transactions on Antennas and Propagation, vol. 28, no. 4, pp. 519-522, 1980. [15] M. Ho, “Simulation of scattered EM fields from rotating cylinder using passing center swing back grids technique in two dimensions,” Progress In Electromagnetics Research, vol. 92, pp. 79-90, 2009.-[16[16] M. Wang, M. Yu, Z. Xu, G. Li, B. Jiang, and J. Xu, “Propagation properties of terahertz waves in a time-varying dusty plasma slab using FDTD,” IEEE Transactions on Plasma Science, vol. 43, no. 12, pp. 4182-4186, 2015.].

In 1887, for the analysis of the Doppler effect in elastic incompressible media, Voigt imposed the convective wave equation (with moving observer) to have the same form as the wave equation with the observer at rest [17[17] W. Voigt, “Ueber das doppler’sche princip,” Nachr Ges Wiss, vol. 41, 1887.]. This made him use auxiliary variables for the three variables of space (in Cartesian coordinates) and the variable of time. Lorentz used Voigt’s auxiliary variables for the investigation of different electromagnetic phenomena with moving objects, based on Maxwell’s electrodynamics [18[18] H. A. Lorentz, “Attempt of a theory of electrical and optical phenomena in moving bodies,” Versuch einer Theorie der Electrischen und Optischen Erscheinungen in Bewegten Körpern. EJ Brill, 1895.]. To account for his hypothesis of length contraction in the Michelson-Morley experiment and for the apparent mass increase of particles in particle accelerators, Lorentz multiplied these variables with the Lorentz factor [19[19] H. A. Lorentz, “Electromagnetic phenomena in a system moving with any velocity smaller than that of light,” Proceedings of the Royal Netherlands Academy of Arts and Sciences, vol. 6, pp. 809-831, 1904.]. Today, the investigations of electromagnetic problems with moving objects are usually based on the utilization of Voigt-Lorentz transformations.

Voigt-Lorentz transformations are applied to calculate Maxwell’s electric and magnetic fields in a frame moving with uniform velocity, from the fields which are known in a reference frame. However, the problem becomes more complex if multiple objects moving at different speeds need to be considered or if non-unform motions need to be studied.

The analysis of electromagnetic problems with objects having non-uniform motion, such as vibration, can be useful in many recent applications, for example, in time-varying waveguides [20[20] S. Taravati and A. A. Kishk, “Space-time modulation: Principles and applications,” IEEE Microwave Magazine, vol. 21, no. 4, pp. 30-56, 2020.] and Doppler radars for the detection of vital signs [21[21] L. Chioukh, H. Boutayeb, D. Deslandes, and K. Wu, “Noise and sensitivity of harmonic radar architecture for remote sensing and detection of vital signs,” IEEE Transactions on Microwave Theory and Techniques, vol. 62, no. 9, pp. 1847-1855, 2014.], [22[22] M. S. Rabbani, J. Churm, and A. P. Feresidis, “Fabry-Pérot beam scanning antenna for remote vital sign detection at 60 GHz,” IEEE Transactions on Antennas and Propagation, vol. 69, no. 6, pp. 3115-3124, 2021.].

The finite-difference time-domain (FDTD) algorithm was introduced in 1966 by Yee, based on finite differences and Yee’s cell, for the numerical resolution of Maxwell’s equations [23[23] K. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Transactions on Antennas and Propagation, vol. 14, no. 3, pp. 302-307, 1966.]. In 1975, the FDTD method was applied to analyze the effect of electromagnetic radiation on human eyes [24[24] A. Taflove and M. E. Brodwin, “Numerical solution of steady-state electromagnetic scattering problems using the timedependent Maxwell’s equations,” IEEE Transactions on Microwave Theory and Techniques, vol. 23, no. 8, pp. 623-630, 1975.]. Nowadays, thanks to its development by numerous pioneers [25[25] A. Reineix and B. Jecko, “Analysis of microstrip patch antennas using finite difference time domain method,” IEEE Transactions on Antennas and Propagation, vol. 37, no. 11, pp. 1361-1369, 1989.][26] D. M. Sheen, S. M. Ali, M. D. Abouzahra, and J.-A. Kong, “Application of the three-dimensional finite-difference timedomain method to the analysis of planar microstrip circuits,” IEEE Transactions on Microwave Theory and Techniques, vol. 38, no. 7, pp. 849-857, 1990.-[27[27] K. S. Yee, D. Ingham, and K. Shlager, “Time-domain extrapolation to the far field based on FDTD calculations,” IEEE Transactions on Antennas and Propagation, vol. 39, no. 3, pp. 410-413, 1991.], the FDTD method is used in a wide range of applications from DC to optics.

The FDTD method has been successfully used for analyzing electromagnetic problems involving moving objects with uniform speeds [1[1] M. Marvasti and H. Boutayeb, “FDTD analysis of the Sagnac effect employed in the global positioning system,” IEEE Transactions on Antennas and Propagation, vol. 71, pp. 9119-9123, 2023. DOI: 10.1109/TAP.2023.3296903.
https://doi.org/10.1109/TAP.2023.3296903...
], [28[28] F. Harfoush, A. Taflove, and G. A. Kriegsmann, “A numerical technique for analyzing electromagnetic wave scattering from moving surfaces in one and two dimensions,” IEEE Transactions on Antennas and Propagation, vol. 37, no. 1, pp. 55-63, 1989.][29] M. J. Inman, A. Z. Elsherbeni, and C. Smith, “Finite difference time domain simulation of moving objects,” in Proceedings of the 2003 IEEE Radar Conference (Cat. No. 03CH37474), pp. 439-445, 2003. [30] K. Zheng, X. Liu, Z. Mu, and G. Wei, “Analysis of scattering fields from moving multilayered dielectric slab illuminated by an impulse source,” IEEE Antennas and Wireless Propagation Letters, vol. 16, pp. 2130-2133, 2017. [31] K.-S. Zheng, J.-Z. Li, G. Wei, and J.-D. Xu, “Analysis of Doppler effect of moving conducting surfaces with LorentzFDTD method,” Journal of Electromagnetic Waves and Applications, vol. 27, no. 2, pp. 149-159, 2013. [32] Y. Liu, K. Zheng, Z. Mu, and X. Liu, “Reflection and transmission coefficients of moving dielectric in half space,” in 2016 11th International Symposium on Antennas, Propagation and EM Theory (ISAPE), pp. 485-487, 2016. [33] K. Zheng, Z. Mu, H. Luo, and G. Wei, “Electromagnetic properties from moving dielectric in high speed with lorentzfdtd,” IEEE Antennas and Wireless Propagation Letters, vol. 15, pp. 934-937, 2015. [34] Y. Li, K. Zheng, Y. Liu, and L. Xu, “Radiated fields of a high-speed moving dipole at oblique incidence,” in 2017 International Applied Computational Electromagnetics Society Symposium (ACES), pp. 1-2, 2017. [35] K. Zheng, Y. Li, X. Tu, and G. Wei, “Scattered fields from a three-dimensional complex target moving at high speed,” in 2018 International Applied Computational Electromagnetics Society Symposium-China (ACES), pp. 1-2, 2018. [36] K. Zheng, Y. Li, L. Xu, J. Li, and G. Wei, “Electromagnetic properties of a complex pyramid-shaped target moving at high speed,” IEEE Transactions on Antennas and Propagation, vol. 66, no. 12, pp. 7472-7476, 2018. [37] K. Zheng, Y. Li, S. Qin, K. An, and G. Wei, “Analysis of micromotion characteristics from moving conical-shaped targets using the Lorentz-FDTD method,” IEEE Transactions on Antennas and Propagation, vol. 67, no. 11, pp. 7174-7179, 2019. [38] M. Marvasti and H. Boutayeb, “Analysis of moving dielectric half-space with oblique plane wave incidence using the Finite Difference Time Domain method,” Progress In Electromagnetics Research M, vol. 115, pp. 119-128, 2023. [39] M. Marvasti and H. Boutayeb, “Analysis of moving bodies with the FDTD method,” in 2023 17th European Conference on Antennas and Propagation (EuCAP), pp. 1-5, 2023. [40] M. Marvasti and H. Boutayeb, “Non-relativistic finite difference time domain method for electromagnetic problems with moving bodies,” in 2023 IEEE MTT-S International Conference on Numerical Electromagnetic and Multiphysics Modeling and Optimization (NEMO), pp. 180-183, 2023. [41] M. Marvasti and H. Boutayeb, “Numerical study of electromagnetic waves with sources, observer, and scattering objects in motion,” IEEE Transactions on Microwave Theory and Techniques, 2023.-[42[42] M. Marvasti and H. Boutayeb, “Electromagnetic analysis of moving structures in a moving reference frame,” The Journal of Engineering, vol. 2023, no. 11, p. e12302, 2023.]. In some of these works, such as in [37[37] K. Zheng, Y. Li, S. Qin, K. An, and G. Wei, “Analysis of micromotion characteristics from moving conical-shaped targets using the Lorentz-FDTD method,” IEEE Transactions on Antennas and Propagation, vol. 67, no. 11, pp. 7174-7179, 2019.], techniques were proposed for the implementation of Voigt-Lorentz transformations in FDTD. In other papers, such as [1[1] M. Marvasti and H. Boutayeb, “FDTD analysis of the Sagnac effect employed in the global positioning system,” IEEE Transactions on Antennas and Propagation, vol. 71, pp. 9119-9123, 2023. DOI: 10.1109/TAP.2023.3296903.
https://doi.org/10.1109/TAP.2023.3296903...
], [38[38] M. Marvasti and H. Boutayeb, “Analysis of moving dielectric half-space with oblique plane wave incidence using the Finite Difference Time Domain method,” Progress In Electromagnetics Research M, vol. 115, pp. 119-128, 2023.]-[42[42] M. Marvasti and H. Boutayeb, “Electromagnetic analysis of moving structures in a moving reference frame,” The Journal of Engineering, vol. 2023, no. 11, p. e12302, 2023.], interesting and meaningful results have been obtained by using a direct FDTD approach that did not include Voigt-Lorentz transformations.

The objective of this paper is to use a direct FDTD method to analyze electromagnetic problems with non-uniform movements including vibration, rotation, and acceleration. The theoretical and numerical approaches are described in detail. A series of full-wave simulations are carried out, for the following problems with non-uniform motions: vibrating metallic slab, rotating observation point, rotating line source, accelerating observation point, accelerating plane wave source, and accelerating electromagnetic reflecting surface. For all these problems, the numerical results provide insight into the physical processes involved.

The remainder of the paper is organized as follows. Section II presents and discusses the main points of the theoretical aspects: Maxwell’s equations, the FDTD method, the approach based on a change of the reference frame for moving bodies, and the proposed theoretical/numerical approach. Numerical results are shown and discussed in Section III. Different electromagnetic problems are analyzed: the illumination by a plane wave of a vibrating metallic plate with different vibration frequencies, frequency spectrum for a rotating observer and for a rotating line source, electric field observed by an accelerating observer, the waves radiated by an accelerating ideal or resistive plane wave source, and reflected waves from an accelerating electromagnetic surface illuminated by a plane wave. Concluding remarks are given in Section IV.

II. THEORY

A. Maxwell’s equations

In an isotropic medium, Maxwell’s equations can be written as:

(1) × E = - μ H t
(2) × H = σ E + ε E t

where E is the electric field and H is the magnetic field. µ (permeability), ɛ (permittivity), and σ (conductivity) are the constitutive parameters of the medium. In Cartesian coordinates, six partial differential equations are obtained:

(3) H x t = 1 μ ( E y z - E z y )
(4) H y t = 1 μ ( E z x - E x z )
(5) H z t = 1 μ ( E x y - E y x )
(6) E x t = 1 ε ( H z y - H y z - σ E x )
(7) E y t = 1 ε ( H x z - H z x - σ E y )
(8) E z t = 1 ε ( H y x - H x y - σ E z )

B. FDTD method

FDTD is based on a discretization of Maxwell’s equation in time and space:

(9) F n ( i , j , k ) = ( i δ x , j δ y , k δ z , n δ t )

where F is an electric or magnetic field component. The application of second-order error central finite difference for a variable of space can be written:

(10) F n ( i , j , k ) x = F n ( i + 1 / 2 , j , k ) - F n ( i - 1 / 2 , j , k ) δ x + O ( δ x 2 )

For time variable, the following can be obtained:

(11) F n ( i , j , k ) t = F n + 1 / 2 ( i , j , k ) - F n - 1 / 2 ( i , j , k ) δ t + O ( δ t 2 )

The application of finite differences in (3)-(8) gives equations for the numerical updates of electric and magnetic field components. For example, the update equation for Hx in Yee’s algorithm is:

(12) H x n + 1 / 2 ( i , j + 1 / 2 , k + 1 / 2 ) = H x n - 1 / 2 ( i , j + 1 / 2 , k + 1 / 2 ) + δ t μ ( i , j + 1 / 2 , k + 1 / 2 ) δ × { E y n ( i , j + 1 / 2 , k + 1 ) - E y n ( i , j + 1 / 2 , k ) + E z n ( i , j , k + 1 / 2 ) - E z n ( i , j + 1 , k + 1 / 2 ) }

The Yee’s cell [23[23] K. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Transactions on Antennas and Propagation, vol. 14, no. 3, pp. 302-307, 1966.] is used for the positions of the different field components. Fig. 1 shows Yee’s algorithm with black curves. The algorithm presents a time loop. At each time step n, the electric field En can be updated using magnetic field Hn-1/2, and electric field En-1, from previous time step. Then, Hn+1/2 is calculated by using Hn-1/2 and En. The next iteration is obtained by the incrementation n = n + 1. The algorithm ends when n > Nbiteration.

C. Method based on a change of the reference frame and Voigt-Lorentz transformations for the analysis of moving bodies

The general method considered in the literature for analyzing electromagnetic problems involving motion is based on Voigt-Lorentz transformations between a rest frame and a moving frame (with uniform velocity v=vx^). The time variable and the space variables are transformed as:

(13) t = γ ( t - v x c 2 )
(14) x = γ ( x - v t )
(15) y = y
(16) z = z

where the γ=1/1-v2c2 is the Lorentz factor. The transformations of the components of the electromagnetic fields can be written [43[43] A. Einstein, “On the electrodynamics of moving bodies,” Annalen der physik, vol. 17, no. 10, pp. 891-921, 1905.], [44[44] W. Engelhardt, “On the relativistic transformation of electromagnetic fields,” Apeiron, vol. 11, no. 2, p. 309, 2004.]:

(17) E x = E x
(18) E y = γ ( E y - v c H z )
(19) E z = γ ( E z + v c H y )
(20) H x = H x
(21) H y = γ ( H y + v c E z )
(22) H z = γ ( H z - v c E y )

Finally, the transformed Maxwell’s equations can be expressed as follows:

(23) × E = - μ H t
(24) × H = σ E + ε E t

where ∇ denotes the derivatives with respect to transformed space variables.

D. Proposed approach

Lorentz aether theory [19[19] H. A. Lorentz, “Electromagnetic phenomena in a system moving with any velocity smaller than that of light,” Proceedings of the Royal Netherlands Academy of Arts and Sciences, vol. 6, pp. 809-831, 1904.] and Einstein’s special theory of relativity [43[43] A. Einstein, “On the electrodynamics of moving bodies,” Annalen der physik, vol. 17, no. 10, pp. 891-921, 1905.] use both the same approach described in II-C. However, these two theories differ in the physical interpretation of Voigt-Lorentz transformations of space and time variables. In Lorentz aether theory, except for the γ factor which adds a physical effect that is not present in Maxwell’s equations, the transformed space and time variables are not physical. In the special theory of relativity, Voigt-Lorentz transformations are physical.

In this work, the framework of Lorentz aether theory is preferred and the length contraction effect (γ factor) is ignored. Indeed, the γ factor can be considered equal to one for small v/c.

Furthermore, in 2004, Engelhardt raised important concerns about the relativistic equations for the transformed fields and emphasized the necessity for developing an electromagnetic theory for moving matter [44[44] W. Engelhardt, “On the relativistic transformation of electromagnetic fields,” Apeiron, vol. 11, no. 2, p. 309, 2004.].

An alternative approach consists of using numerical solutions of Maxwell’s equations and changing the positions of objects in the desired directions with time. In such a method, Maxwell’s equations are applied without modification. Fig. 1 shows in red the proposed implementation of moving objects in the FDTD algorithm. Each object (observer, source, scatterer) can move at a different speed. For example, after every mfixi=[δxviδt] iteration, an object moves by one cell in the algorithm. The electromagnetic fields can also move, as demonstrated in [1[1] M. Marvasti and H. Boutayeb, “FDTD analysis of the Sagnac effect employed in the global positioning system,” IEEE Transactions on Antennas and Propagation, vol. 71, pp. 9119-9123, 2023. DOI: 10.1109/TAP.2023.3296903.
https://doi.org/10.1109/TAP.2023.3296903...
], for the implementation of the Fresnel drag effect or the Sagnac effect in a dielectric medium. The proposed method is not limited to uniform motion and it can be used for problems with multiple objects moving at different speeds.

Fig. 1
Flowchart of the proposed FDTD approach for moving objects. In black: Yee’s algorithm. In red: implementation of movements (source, observer, scatterer, or field) ([.%.] denotes the Modulo operator)

III. NUMERICAL RESULTS

A. Vibratory motion

A vibrating metallic plate, as shown in Fig. 2, is considered. The boundary conditions include Absorbing Boundary Conditions (ABCs), Perfect Magnetic Conductors (PMCs), and Perfect Electric Conductors (PECs). f0 is the frequency of the exciting source and fm is the frequency of vibration of the metallic slab. Fig. 3a and Fig. 3b show the spectrum of the reflected wave for fm > f0 and for fm < f0, respectively. The spectrum shows the presence of the fundamental source frequency and frequencies due to intermodulation with the vibrating frequency. This is in agreement with experimental results obtained with a vibrating metallic plate and a Doppler radar [21[21] L. Chioukh, H. Boutayeb, D. Deslandes, and K. Wu, “Noise and sensitivity of harmonic radar architecture for remote sensing and detection of vital signs,” IEEE Transactions on Microwave Theory and Techniques, vol. 62, no. 9, pp. 1847-1855, 2014.]. As an extension of this work, it would be possible to analyze the signal received in the context of the motion of the human chest which is due to both respiration and heartbeat [21[21] L. Chioukh, H. Boutayeb, D. Deslandes, and K. Wu, “Noise and sensitivity of harmonic radar architecture for remote sensing and detection of vital signs,” IEEE Transactions on Microwave Theory and Techniques, vol. 62, no. 9, pp. 1847-1855, 2014.], [22[22] M. S. Rabbani, J. Churm, and A. P. Feresidis, “Fabry-Pérot beam scanning antenna for remote vital sign detection at 60 GHz,” IEEE Transactions on Antennas and Propagation, vol. 69, no. 6, pp. 3115-3124, 2021.]. Instead of using a simple sinusoid for the movement of the metallic slab, one could use a more realistic motion of the human chest.

Fig. 2
Vibrating metallic plate in FDTD.

Fig. 3
Simulated spectrum of the reflected wave, for vibrating metallic plate, with frequency of vibration (a) fm > f0 (f0=10GHz; fm=25GHz, 35GHz, or 60GHz), and (b) fm < f0 (f0=10GHz; fm=4GHz or 5GHz).

B. Rotatory motion

1) Rotating observation point:Fig. 4 shows a rotating observer illuminated by a plane wave. f0 is the frequency of the exciting source and fr is the frequency of rotation of the observer. Fig. 5a and Fig. 5b show the spectrum of the observed signal for fr > f0 and for fr < f0, respectively. In these results, the frequencies due to intermodulation can be identified.

Fig. 4
Rotating observer.

Fig. 5
Spectrum of observed signal, for rotating observation point, with the frequency of rotation (a) fr > f0 (f0=10GHz; fr=25GHz, 35GHz, or 60GHz), and (b) fr < f0 (f0=10GHz; fr=4GHz or 5GHz).

2) Rotating source: A line source is now rotating as shown in Fig. 6 and the observer is at rest. f0 is the frequency of the exciting source and fr is the frequency of rotation of the line source. Fig. 7a and Fig. 7b show the spectrum of the observed signal for fr > f0 and fr < f0, respectively. The results are again meaningful. Fig. 8a and Fig. 8b show the electric field distribution, at a time instant, for the rotating source, with fr > f0 and fr < f0, respectively.

Fig. 6
Rotating line source.

Fig. 7
Spectrum of observed signal, for rotating ideal line source, with frequency of rotation (a) fr > f0 (f0=10GHz; fr=25GHz, 35GHz, or 60GHz), and (b) fr < f0 (f0=10GHz; fr=4GHz or 5GHz).

Fig. 8
Simulated electric field distribution for a rotating ideal line source: (a) fr > f0 and (b) fr < f0.

C. Acceleration

Acceleration can be implemented in the proposed numerical approach by changing the speed of motion of an object in the FDTD time loop.

1) Accelerating observer: An observer is accelerating as illustrated in Fig. 9. As shown in Fig. 10, depending on the direction (±xˆ) of the observer motion (the plane wave is propagating toward +xˆ direction), the frequency of the signal decreases or increases with time, as one can expect. This results in clear chirp signals, which have applications in radars.

Fig. 9
Accelerating observation point.

Fig. 10
Observed signals in the time domain, for accelerating observer.

2) Accelerating plane wave source: In this subsection, a monochromatic plane wave source is accelerating and the observer is at rest, as shown in Fig. 11. A distinction is made between the ideal plane wave source (made of current sources) and the resistive plane wave source (made of current sources with resistors having low resistance). The wave observed in the time domain is shown in Fig. 12 for the two cases (the source moves toward the observer). With the ideal source, which is not realistic, the field can grow without bounds if the source accelerates toward the observer. The field does not increase in the case of a resistive plane wave source. Fig. 13a, Fig. 13b, and Fig. 13c show the numerical field distributions at different time instants for an accelerating ideal plane wave source. Two chirp waves are propagating in -xˆ and +xˆ directions, with the speed of light c.

Fig. 11
Accelerating source.

Fig. 12
Radiated waves in the time domain, for accelerating ideal or resistive plane wave sources.

3) Accelerating reflecting electromagnetic surface: A periodic reflecting electromagnetic surface, also called a partially reflecting surface (PRS), is now considered. The PRS is made of infinitely long metallic wires and it is accelerating toward the observer, as shown in Fig. 14. The observer and the monochromatic source are at rest. Fig. 15 shows the reflected wave in the time domain. One can note that the frequency of the signal increases logically with time (up-chirp).

Fig. 13
Electric field distribution for accelerating ideal plane wave source toward -x direction, at different time instants: (a) t = 0.2ns, (b) t = 0.6ns, and (c) t = 1ns. Two chirp waves are propagating in -x and +x directions, with the speed of light c.

Fig. 14
Accelerating PRS in FDTD.

Fig. 15
Reflected plane wave from accelerating PRS.

IV. CONCLUSION

The concept of the change of the reference frame and the utilization of Voigt-Lorentz transformations in Maxwell’s equations are usually applied in the literature for electromagnetic problems with uniform motion, where only one structure or system is moving with a single constant speed. In addition to these restrictions, important concerns have been raised about such an approach. In order to account for more non-uniform motions, a general electromagnetic approach for moving bodies has been proposed in this paper. The proposed method consists of using Maxwell’s equations without modification and moving the objects numerically. With the proposed FDTD algorithm, non-uniform motions (vibration, rotation, acceleration) have been considered. Different problems have been analyzed numerically: a vibrating metallic slab, a rotating observer, and a rotating source. Physical insight has been provided by the spectrums of the reflected wave from a vibrating metallic slab, as well as the observed signal for a rotating observer or a rotating line source. This paper has also presented the analysis of accelerating observer, accelerating plane wave source, and accelerating electromagnetic reflecting surface. The ideal and resistive plane wave sources have been distinguished. It is shown that the amplitude of the field increases with an ideal (non-realistic) approaching source but does not increase with a resistive plane wave source. Furthermore, accelerating plane reflectors could be used to generate chirp signals which are used in radars.

Future work will concentrate on the utilization of the proposed method for the development of Doppler radars used in the detection of vital signs. With the FDTD method, it could be possible to model a realistic motion of the human chest based on respiration and heartbeat.

REFERENCES

  • [1]
    M. Marvasti and H. Boutayeb, “FDTD analysis of the Sagnac effect employed in the global positioning system,” IEEE Transactions on Antennas and Propagation, vol. 71, pp. 9119-9123, 2023. DOI: 10.1109/TAP.2023.3296903.
    » https://doi.org/10.1109/TAP.2023.3296903.
  • [2]
    J. Cooper, “Scattering of electromagnetic fields by a moving boundary: the one-dimensional case,” IEEE Transactions on Antennas and Propagation, vol. 28, no. 6, pp. 791-795, 1980.
  • [3]
    B. Michielsen, G. Herman, A. De Hoop, and D. De Zutter, “Three-dimensional relativistic scattering of electromagnetic waves by an object in uniform translational motion,” Journal of Mathematical Physics, vol. 22, no. 11, pp. 2716-2722, 1981.
  • [4]
    N. Engheta, M. W. Kowarz, and D. L. Jaggard, “Effect of chirality on the Doppler shift and aberration of light waves,” Journal of Applied Physics, vol. 66, no. 6, pp. 2274-2277, 1989.
  • [5]
    C. Neipp, A. Hernández, J. Rodes, A. Márquez, T. Beléndez, and A. Beléndez, “An analysis of the classical Doppler effect,” European Journal of Physics, vol. 24, no. 5, p. 497, 2003.
  • [6]
    S. Sahrani and M. Kuroda, “Numerical analysis of the electromagnetic wave scattering from a moving dielectric body by overset grid generation method,” in 2012 IEEE Asia-Pacific Conference on Applied Electromagnetics (APACE), pp. 6-10, 2012.
  • [7]
    J. Cooper, “Long-time behavior and energy growth for electromagnetic waves reflected by a moving boundary,” IEEE Transactions on Antennas and Propagation, vol. 41, no. 10, pp. 1365-1370, 1993.
  • [8]
    S. Borkar and R. Yang, “Scattering of electromagnetic waves from rough oscillating surfaces using spectral Fourier method,” IEEE Transactions on Antennas and Propagation, vol. 21, no. 5, pp. 734-736, 1973.
  • [9]
    R. Kleinman and R. Mack, “Scattering by linearly vibrating objects,” IEEE Transactions on Antennas and Propagation, vol. 27, no. 3, pp. 344-352, 1979.
  • [10]
    D. De Zutter, “Reflections from linearly vibrating objects: plane mirror at oblique incidence,” IEEE Transactions on Antennas and Propagation, vol. 30, no. 5, pp. 898-903, 1982.
  • [11]
    M. Ho, “Numerical simulation of scattering of electromagnetic waves from traveling and/or vibrating perfect conducting planes,” IEEE Transactions on Antennas and Propagation, vol. 54, no. 1, pp. 152-156, 2006.
  • [12]
    V. C. Chen, F. Li, S.-S. Ho, and H. Wechsler, “Micro-Doppler effect in radar: phenomenon, model, and simulation study,” IEEE Transactions on Aerospace and Electronic Systems, vol. 42, no. 1, pp. 2-21, 2006.
  • [13]
    H. Gao, L. Xie, S. Wen, and Y. Kuang, “Micro-Doppler signature extraction from ballistic target with micro-motions,” IEEE Transactions on Aerospace and Electronic Systems, vol. 46, no. 4, pp. 1969-1982, 2010.
  • [14]
    N. Engheta, A. Mickelson, and C. Papas, “On the near-zone inverse Doppler effect,” IEEE Transactions on Antennas and Propagation, vol. 28, no. 4, pp. 519-522, 1980.
  • [15]
    M. Ho, “Simulation of scattered EM fields from rotating cylinder using passing center swing back grids technique in two dimensions,” Progress In Electromagnetics Research, vol. 92, pp. 79-90, 2009.
  • [16]
    M. Wang, M. Yu, Z. Xu, G. Li, B. Jiang, and J. Xu, “Propagation properties of terahertz waves in a time-varying dusty plasma slab using FDTD,” IEEE Transactions on Plasma Science, vol. 43, no. 12, pp. 4182-4186, 2015.
  • [17]
    W. Voigt, “Ueber das doppler’sche princip,” Nachr Ges Wiss, vol. 41, 1887.
  • [18]
    H. A. Lorentz, “Attempt of a theory of electrical and optical phenomena in moving bodies,” Versuch einer Theorie der Electrischen und Optischen Erscheinungen in Bewegten Körpern. EJ Brill, 1895.
  • [19]
    H. A. Lorentz, “Electromagnetic phenomena in a system moving with any velocity smaller than that of light,” Proceedings of the Royal Netherlands Academy of Arts and Sciences, vol. 6, pp. 809-831, 1904.
  • [20]
    S. Taravati and A. A. Kishk, “Space-time modulation: Principles and applications,” IEEE Microwave Magazine, vol. 21, no. 4, pp. 30-56, 2020.
  • [21]
    L. Chioukh, H. Boutayeb, D. Deslandes, and K. Wu, “Noise and sensitivity of harmonic radar architecture for remote sensing and detection of vital signs,” IEEE Transactions on Microwave Theory and Techniques, vol. 62, no. 9, pp. 1847-1855, 2014.
  • [22]
    M. S. Rabbani, J. Churm, and A. P. Feresidis, “Fabry-Pérot beam scanning antenna for remote vital sign detection at 60 GHz,” IEEE Transactions on Antennas and Propagation, vol. 69, no. 6, pp. 3115-3124, 2021.
  • [23]
    K. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Transactions on Antennas and Propagation, vol. 14, no. 3, pp. 302-307, 1966.
  • [24]
    A. Taflove and M. E. Brodwin, “Numerical solution of steady-state electromagnetic scattering problems using the timedependent Maxwell’s equations,” IEEE Transactions on Microwave Theory and Techniques, vol. 23, no. 8, pp. 623-630, 1975.
  • [25]
    A. Reineix and B. Jecko, “Analysis of microstrip patch antennas using finite difference time domain method,” IEEE Transactions on Antennas and Propagation, vol. 37, no. 11, pp. 1361-1369, 1989.
  • [26]
    D. M. Sheen, S. M. Ali, M. D. Abouzahra, and J.-A. Kong, “Application of the three-dimensional finite-difference timedomain method to the analysis of planar microstrip circuits,” IEEE Transactions on Microwave Theory and Techniques, vol. 38, no. 7, pp. 849-857, 1990.
  • [27]
    K. S. Yee, D. Ingham, and K. Shlager, “Time-domain extrapolation to the far field based on FDTD calculations,” IEEE Transactions on Antennas and Propagation, vol. 39, no. 3, pp. 410-413, 1991.
  • [28]
    F. Harfoush, A. Taflove, and G. A. Kriegsmann, “A numerical technique for analyzing electromagnetic wave scattering from moving surfaces in one and two dimensions,” IEEE Transactions on Antennas and Propagation, vol. 37, no. 1, pp. 55-63, 1989.
  • [29]
    M. J. Inman, A. Z. Elsherbeni, and C. Smith, “Finite difference time domain simulation of moving objects,” in Proceedings of the 2003 IEEE Radar Conference (Cat. No. 03CH37474), pp. 439-445, 2003.
  • [30]
    K. Zheng, X. Liu, Z. Mu, and G. Wei, “Analysis of scattering fields from moving multilayered dielectric slab illuminated by an impulse source,” IEEE Antennas and Wireless Propagation Letters, vol. 16, pp. 2130-2133, 2017.
  • [31]
    K.-S. Zheng, J.-Z. Li, G. Wei, and J.-D. Xu, “Analysis of Doppler effect of moving conducting surfaces with LorentzFDTD method,” Journal of Electromagnetic Waves and Applications, vol. 27, no. 2, pp. 149-159, 2013.
  • [32]
    Y. Liu, K. Zheng, Z. Mu, and X. Liu, “Reflection and transmission coefficients of moving dielectric in half space,” in 2016 11th International Symposium on Antennas, Propagation and EM Theory (ISAPE), pp. 485-487, 2016.
  • [33]
    K. Zheng, Z. Mu, H. Luo, and G. Wei, “Electromagnetic properties from moving dielectric in high speed with lorentzfdtd,” IEEE Antennas and Wireless Propagation Letters, vol. 15, pp. 934-937, 2015.
  • [34]
    Y. Li, K. Zheng, Y. Liu, and L. Xu, “Radiated fields of a high-speed moving dipole at oblique incidence,” in 2017 International Applied Computational Electromagnetics Society Symposium (ACES), pp. 1-2, 2017.
  • [35]
    K. Zheng, Y. Li, X. Tu, and G. Wei, “Scattered fields from a three-dimensional complex target moving at high speed,” in 2018 International Applied Computational Electromagnetics Society Symposium-China (ACES), pp. 1-2, 2018.
  • [36]
    K. Zheng, Y. Li, L. Xu, J. Li, and G. Wei, “Electromagnetic properties of a complex pyramid-shaped target moving at high speed,” IEEE Transactions on Antennas and Propagation, vol. 66, no. 12, pp. 7472-7476, 2018.
  • [37]
    K. Zheng, Y. Li, S. Qin, K. An, and G. Wei, “Analysis of micromotion characteristics from moving conical-shaped targets using the Lorentz-FDTD method,” IEEE Transactions on Antennas and Propagation, vol. 67, no. 11, pp. 7174-7179, 2019.
  • [38]
    M. Marvasti and H. Boutayeb, “Analysis of moving dielectric half-space with oblique plane wave incidence using the Finite Difference Time Domain method,” Progress In Electromagnetics Research M, vol. 115, pp. 119-128, 2023.
  • [39]
    M. Marvasti and H. Boutayeb, “Analysis of moving bodies with the FDTD method,” in 2023 17th European Conference on Antennas and Propagation (EuCAP), pp. 1-5, 2023.
  • [40]
    M. Marvasti and H. Boutayeb, “Non-relativistic finite difference time domain method for electromagnetic problems with moving bodies,” in 2023 IEEE MTT-S International Conference on Numerical Electromagnetic and Multiphysics Modeling and Optimization (NEMO), pp. 180-183, 2023.
  • [41]
    M. Marvasti and H. Boutayeb, “Numerical study of electromagnetic waves with sources, observer, and scattering objects in motion,” IEEE Transactions on Microwave Theory and Techniques, 2023.
  • [42]
    M. Marvasti and H. Boutayeb, “Electromagnetic analysis of moving structures in a moving reference frame,” The Journal of Engineering, vol. 2023, no. 11, p. e12302, 2023.
  • [43]
    A. Einstein, “On the electrodynamics of moving bodies,” Annalen der physik, vol. 17, no. 10, pp. 891-921, 1905.
  • [44]
    W. Engelhardt, “On the relativistic transformation of electromagnetic fields,” Apeiron, vol. 11, no. 2, p. 309, 2004.

Publication Dates

  • Publication in this collection
    30 Aug 2024
  • Date of issue
    2024

History

  • Received
    11 Nov 2023
  • Reviewed
    29 Nov 2023
  • Accepted
    21 June 2024
Sociedade Brasileira de Microondas e Optoeletrônica e Sociedade Brasileira de Eletromagnetismo Praça Mauá, n°1, 09580-900 São Caetano do Sul - S. Paulo/Brasil, Tel./Fax: (55 11) 4238 8988 - São Caetano do Sul - SP - Brazil
E-mail: editor_jmoe@sbmo.org.br