{"title":"非均匀介质中的矢量亥姆霍兹电磁波传播器","authors":"Laurence Keefe;Austin McDaniel;Max Cubillos;Ilya Zilberter;Timothy Madden","doi":"10.1109/TAP.2025.3562744","DOIUrl":null,"url":null,"abstract":"The vector electric field Helmholtz equation, containing cross-polarization terms, is factored to produce both pseudo-differential and exponential operator forms of a 3-D, one-way, vector, wave equation for propagation through inhomogeneous media. From this operator factorization, we develop a high-order approximate, vector Helmholtz propagator that correctly handles forward arc, high-angle scattering, and diffraction from inhomogeneities at all resolved length scales and seamlessly includes evanescent waves. Our implementation of the exponential operator form of the one-way propagator is discussed extensively. A rational approximation/partial fraction decomposition of the exponential operator converts the propagator into a moderate number of large, sparse, linear solves whose results are summed together at each step to advance the electric field in space. We use a new AAA-Lawson rational interpolant for this approximation, rather than the more common Padé expansions that have appeared in the seismic and ocean acoustics literature previously. GMRES is used to solve these large systems. A direct-solve, free-space propagation method proves to be an effective preconditioner for GMRES, but can also serve as a standalone propagator in homogeneous media. Scalar computational examples shown include plane-wave diffraction by a circular aperture and Gaussian beam propagation through sine-product and homogeneous refractive index fields. The sine-product example compares its results to that of paraxial propagation through the same media and demonstrates the substantial differences between these propagator paradigms when the scale of the inhomogeneities is of the order of the fundamental wavelength in the Helmholtz equation. We also examine the convergence of the homogeneous media beam results to fields generated by Clenshaw-Curtis (C-C) evaluation of the first Rayleigh-Sommerfeld integral for the same initial conditions.","PeriodicalId":13102,"journal":{"name":"IEEE Transactions on Antennas and Propagation","volume":"73 8","pages":"5823-5838"},"PeriodicalIF":5.8000,"publicationDate":"2025-04-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"A Vector Helmholtz Electromagnetic Wave Propagator for Inhomogeneous Media\",\"authors\":\"Laurence Keefe;Austin McDaniel;Max Cubillos;Ilya Zilberter;Timothy Madden\",\"doi\":\"10.1109/TAP.2025.3562744\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"The vector electric field Helmholtz equation, containing cross-polarization terms, is factored to produce both pseudo-differential and exponential operator forms of a 3-D, one-way, vector, wave equation for propagation through inhomogeneous media. From this operator factorization, we develop a high-order approximate, vector Helmholtz propagator that correctly handles forward arc, high-angle scattering, and diffraction from inhomogeneities at all resolved length scales and seamlessly includes evanescent waves. Our implementation of the exponential operator form of the one-way propagator is discussed extensively. A rational approximation/partial fraction decomposition of the exponential operator converts the propagator into a moderate number of large, sparse, linear solves whose results are summed together at each step to advance the electric field in space. We use a new AAA-Lawson rational interpolant for this approximation, rather than the more common Padé expansions that have appeared in the seismic and ocean acoustics literature previously. GMRES is used to solve these large systems. A direct-solve, free-space propagation method proves to be an effective preconditioner for GMRES, but can also serve as a standalone propagator in homogeneous media. Scalar computational examples shown include plane-wave diffraction by a circular aperture and Gaussian beam propagation through sine-product and homogeneous refractive index fields. The sine-product example compares its results to that of paraxial propagation through the same media and demonstrates the substantial differences between these propagator paradigms when the scale of the inhomogeneities is of the order of the fundamental wavelength in the Helmholtz equation. We also examine the convergence of the homogeneous media beam results to fields generated by Clenshaw-Curtis (C-C) evaluation of the first Rayleigh-Sommerfeld integral for the same initial conditions.\",\"PeriodicalId\":13102,\"journal\":{\"name\":\"IEEE Transactions on Antennas and Propagation\",\"volume\":\"73 8\",\"pages\":\"5823-5838\"},\"PeriodicalIF\":5.8000,\"publicationDate\":\"2025-04-25\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"IEEE Transactions on Antennas and Propagation\",\"FirstCategoryId\":\"94\",\"ListUrlMain\":\"https://ieeexplore.ieee.org/document/10977779/\",\"RegionNum\":1,\"RegionCategory\":\"计算机科学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"ENGINEERING, ELECTRICAL & ELECTRONIC\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"IEEE Transactions on Antennas and Propagation","FirstCategoryId":"94","ListUrlMain":"https://ieeexplore.ieee.org/document/10977779/","RegionNum":1,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"ENGINEERING, ELECTRICAL & ELECTRONIC","Score":null,"Total":0}
A Vector Helmholtz Electromagnetic Wave Propagator for Inhomogeneous Media
The vector electric field Helmholtz equation, containing cross-polarization terms, is factored to produce both pseudo-differential and exponential operator forms of a 3-D, one-way, vector, wave equation for propagation through inhomogeneous media. From this operator factorization, we develop a high-order approximate, vector Helmholtz propagator that correctly handles forward arc, high-angle scattering, and diffraction from inhomogeneities at all resolved length scales and seamlessly includes evanescent waves. Our implementation of the exponential operator form of the one-way propagator is discussed extensively. A rational approximation/partial fraction decomposition of the exponential operator converts the propagator into a moderate number of large, sparse, linear solves whose results are summed together at each step to advance the electric field in space. We use a new AAA-Lawson rational interpolant for this approximation, rather than the more common Padé expansions that have appeared in the seismic and ocean acoustics literature previously. GMRES is used to solve these large systems. A direct-solve, free-space propagation method proves to be an effective preconditioner for GMRES, but can also serve as a standalone propagator in homogeneous media. Scalar computational examples shown include plane-wave diffraction by a circular aperture and Gaussian beam propagation through sine-product and homogeneous refractive index fields. The sine-product example compares its results to that of paraxial propagation through the same media and demonstrates the substantial differences between these propagator paradigms when the scale of the inhomogeneities is of the order of the fundamental wavelength in the Helmholtz equation. We also examine the convergence of the homogeneous media beam results to fields generated by Clenshaw-Curtis (C-C) evaluation of the first Rayleigh-Sommerfeld integral for the same initial conditions.
期刊介绍:
IEEE Transactions on Antennas and Propagation includes theoretical and experimental advances in antennas, including design and development, and in the propagation of electromagnetic waves, including scattering, diffraction, and interaction with continuous media; and applications pertaining to antennas and propagation, such as remote sensing, applied optics, and millimeter and submillimeter wave techniques