An efficient and stable spectral method for electromagnetic scattering from a layered periodic structure

The scattering of acoustic and electromagnetic waves by periodic structures plays an important role in a wide range of problems of scientific and technological interest. This contribution focuses upon the stable and high-order numerical simulation of the interaction of time-harmonic electromagnetic waves incident upon a periodic doubly layered dielectric media with sharp, irregular interface. We describe a boundary perturbation method for this problem which avoids not only the need for specialized quadrature rules but also the dense linear systems characteristic of boundary integral/element methods. Additionally, it is a provably stable algorithm as opposed to other boundary perturbation approaches such as Bruno and Reitich’s “method of field expansions” or Milder’s “method of operator expansions”. Our spectrally accurate approach is a natural extension of the “method of transformed field expansions” originally described by Nicholls and Reitich (and later refined to other geometries by the authors) in the single-layer case.     

A high accuracy algorithm for 3D periodic electromagnetic scattering

We develop a highly accurate numerical method for scattering of 3D electromagnetic waves by doubly periodic structures. We approximate scattered fields using the Müller boundary integral formulation of Maxwell’s equations. The accuracy is achieved as singularities are isolated through the use of partitions of unity, leaving smooth, periodic integrands that can be evaluated with high accuracy using trapezoid sums. The removed singularities are resolved through a transformation to polar coordinates. The method relies on the ideas used in the free space scattering algorithm of Bruno and Kunyansky.     

Experimental observation of strong photon localization in disordered photonic crystal waveguides

We demonstrate experimentally that structural perturbations imposed on highly dispersive photonic crystal-based waveguides give rise to spectral features that bear signatures of Anderson localization. Sharp resonances with effective Q's of over 30 000 are found in scattering spectra of disordered waveguides. The resonances are observed in a approximately 20-nm bandwidth centered at the cutoff of slowly guided Bloch modes. The origin of the spectral features can be explained by the interference of coherently scattered electromagnetic waves which results in the formation of a narrow impurity (or localization) band populated with spectrally distinct quasistates. Standard photon localization criteria are fulfilled in the localization band.     

Improving the Accuracy of the Fast Inverse Square Root by Modifying Newton–Raphson Corrections

Direct computation of functions using low-complexity algorithms can be applied both for hardware constraints and in systems where storage capacity is a challenge for processing a large volume of data. We present improved algorithms for fast calculation of the inverse square root function for single-precision and double-precision floating-point numbers. Higher precision is also discussed. Our approach consists in minimizing maximal errors by finding optimal magic constants and modifying the Newton–Raphson coefficients. The obtained algorithms are much more accurate than the original fast inverse square root algorithm and have similar very low computational costs.     

应用改进的特征基函数法和自适应交叉近似算法快速分析导体目标电磁散射特性

特征基函数的构造是特征基函数法的关键步骤之一, 传统方法在构造特征基函数时, 需要在每个子域设置足够多的入射波激励, 生成的特征基函数个数较多, 奇异值分解时间较长. 为了加快特征基函数的构造, 本文提出了一种改进的特征基函数法. 该方法充分考虑每个子域之间的耦合作用, 求出每个子域的次要特征基函数, 从而降低入射波激励的个数, 大大减少了特征基函数的个数; 并且结合自适应交叉近似算法加速阻抗矩阵元素的计算, 提高了次要特征基函数求解和缩减矩阵构建过程中的矩阵矢量相乘的速度. 数值结果证明了本文方法的精确性和高效性. 关键词： 电磁散射 矩量法 特征基函数法 自适应交叉近似     

A new multichannel time reversal focusing method for circumferential Lamb waves and its applications for defect detection in thick-walled pipe with large-diameter

This paper proposes a new multichannel time reversal focusing (MTRF) method for circumferential Lamb waves which is based on modified time reversal algorithm and applies this method for detecting different kinds of defects in thick-walled pipe with large-diameter. The principle of time reversal of circumferential Lamb waves in pipe is presented along with the influence from multiple guided wave modes and propagation paths. Experimental study is carried out in a thick-walled and large-diameter pipe with three artificial defects, namely two axial notches on its inner and outer surface respectively, and a corrosion-like defect on its outer surface. By using the proposed MTRF method, the multichannel signals focus at the defects, leading to the amplitude improvement of the defect scattered signal. Besides, another energy focus arises in the direct signal due to the partial compensation of dispersion and multimode of circumferential Lamb waves, alongside the multichannel focusing, during MTRF process. By taking the direct focus as a time base, accurate defect localization is implemented. Secondly, a new phenomenon is exhibited in this paper that defect scattered wave packet appears just before the right boundary of truncation window after time reversal, and to which two feasible explanations are given. Moreover, this phenomenon can be used as the theoretical basis in the determination of defect scattered waves in time reversal response signal. At last, in order to detect defects without prior knowing their exact position, a large-range truncation window is used in the proposed method. As a result, the experimental operation of MTRF method is simplified and defect detection and localization are well accomplished.     

A robust negative obstacle detection method using seed-growing and dynamic programming for visually-impaired/blind persons

Though a significant amount of work has been done on detecting obstacles, not much attention has been given to the detection of drop offs, e.g., sidewalk curbs, downward stairs, and other hazards. In this paper, we propose algorithms for detecting negative obstacles in an urban setting using stereo vision and two-stage dynamic programming (TSDP) technique. We are developing computer vision algorithms for sensing important terrain features as an aid to blind navigation, which interpret visual information obtained from images collected by cameras mounted on camera legs nearly as high as young person. This paper focuses specifically on a novel computer vision algorithm for detecting negative obstacles (i.e. anything below the level of the ground, such as holes and drop-offs), which are important and ubiquitous features on and near sidewalks and other walkways. The proposed algorithm is compared to other algorithms such as belief propagation and random growing correspondence seeds (GCS). According to the results, the proposed method achieves higher speed, more accurate disparity map and lower RMS errors. The speed of the proposed algorithm is about 28% higher than the random GCS algorithm. We demonstrate experimental results on typical sidewalk scenes to show the effectiveness of the proposed method.     

基于电磁拓扑的多腔体屏蔽效能快速算法

提出基于电磁拓扑理论计算开孔多腔体屏蔽效能的快速方法.首先给出双腔体等效电路和电磁拓扑信号流图,并推导孔缝节点处的散射矩阵,给出拓扑网络的散射矩阵方程和传输矩阵方程,获得双腔体的广义Baum-Liu-Tesche(BLT)方程.在此基础上研究了开孔三腔体,包括串型级联三腔体和串并型混合级联三腔体的广义BLT方程.对于串型级联三腔体,其电磁拓扑网络和广义BLT方程在双腔体基础上直接扩展即可获得.而对于串并型混合级联三腔体,通过将位于三腔体公共面上的孔缝等效为三端口网络节点,并根据三端口网络散射参数定义推导获得该节点的散射矩阵,最终得到串并型混合级联三腔体的广义BLT方程.本文方法对双腔体的计算结果与文献结果和实验结果相符合,对3组不同类型和尺寸开孔腔的屏蔽效能的计算结果与时域有限差分法计算结果符合较好.该算法不仅效率高,通过对所有计算结果和实验结果的误差统计分析,表明该算法具有较高的计算准确度.     

左手材料薄板波导中模式之间的正交关系

基于电磁场理论,采用直接计算的方法证明了左手材料薄板波导中各个模式是相互正交的,即不同传播模之间、传播模与辐射模以及不同辐射模之间是相互正交的.这为研究左手材料薄板波导中的模式耦合奠定了基础. 关键词： 左手材料 薄板波导 模式 正交关系     

海洋可控源三维电磁响应显式灵敏度矩阵的快速算法

借助电场耦合势三维有限体积法与直接求解技术,研究建立了一套海洋可控源三维电磁响应显式灵敏度矩阵(或称为Fréchet导数)高效算法.首先,利用Yee氏交错网格和有限体积法对电场混合势Helmholtz方程进行离散处理,建立与移动源电磁场正演模拟相对应的大型代数方程组,再应用直接法得到的逆矩阵和三维线性插值技术事先确定插值算子和投影算子,并利用投影算子与各个发射源离散向量的乘积计算多发射源电磁响应,极大地提高了多发射源电磁场正演模拟效率.在此基础上,根据块状模型和像素模型中异常体电导率分片常数分布特征,将电导率摄动产生的一次散射电流场表示成Yee氏剖分网格上散射电流元的叠加,由投影算子与散射电流元的离散向量的乘积快速计算出电场强度与磁场强度的显式灵敏度矩阵.最后,通过数值计算检验算法的有效性,并通过块状模型与像素模型分别研究海洋可控源电磁响应特征.     

Relationship between the amplitude and phase of a signal scattered by a point-like acoustic inhomogeneity

Wave scattering by a point-like inhomogeneity, i.e., a strong inhomogeneity with infinitesimal dimensions, is described. This type of inhomogeneity model is used in investigating the point-spread functions of different algorithms and systems. Two approaches are used to derive the rigorous relationship between the amplitude and phase of a signal scattered by a point-like acoustic inhomogeneity. The first approach is based on a Marchenko-type equation. The second approach uses the scattering by a scatterer whose size decreases simultaneously with an increase in its contrast. It is shown that the retarded and advanced waves are scattered differently despite the relationship between the phases of the corresponding scattered waves.     

A reduced basis method for electromagnetic scattering by multiple particles in three dimensions

We consider the development of efficient and fast computational methods for parametrized electromagnetic scattering problems involving many scattering three dimensional bodies. The parametrization may describe the location, orientation, size, shape and number of scattering bodies as well as properties of the source field such as frequency, polarization and incident direction. The emphasis is on problems that need to be solved rapidly to accurately simulate the interaction of scattered fields under parametric variation, e.g., for design, detection, or uncertainty quantification. For such problems, the use of a brute force approach is often ruled out due to the computational cost associated with solving the problem for each parameter value.In this work, we propose an iterative reduced basis method based on a boundary element discretization of few reference scatterers to resolve the computationally challenging large scale problem. The approach includes (i) a computationally intensive offline procedure to create a selection of a set of snapshot parameters and the construction of an associated reduced basis for each reference scatterer and (ii) an inexpensive online algorithm to generate the surface current and scattered field of the parametrized configuration, for any choice of parameters within the parameter domains used in the offline procedure. Comparison of our numerical results with directly measured results for some benchmark configurations demonstrate the power of our method to rapidly simulate the interacting electromagnetic fields under parametric variation of the overall multiple particle configuration.     

Accelerated Stochastic Simulation of Large Chemical Systems

For efficient simulation of chemical systems with large number of reactions, we report a fast and exact algorithm for direct simulation of chemical discrete Markov processes. The approach adopts the scheme of organizing the reactions into hierarchical groups. By generating a random number, the selection of the next reaction that actually occurs is accomplished by a few successive selections in the hierarchical groups. The algorithm which is suited for simulating systems with large number of reactions is much faster than the direct method or the optimized direct method. For a demonstration of its efficiency, the accelerated algorithm is applied to simulate the reaction-diffusion Brusselator model on a discretized space.     

Fast volumetric integral-equation solver for acoustic wave propagation through inhomogeneous media

Elements are described of a volumetric integral-equation-based algorithm applicable to accurate large-scale simulations of scattering and propagation of sound waves through inhomogeneous media. The considered algorithm makes possible simulations involving realistic geometries characterized by highly subwavelength details, large density contrasts, and described in terms of several million unknowns. The algorithm achieves its competitive performance, characterized by O(N log N) solution complexity and O(N) memory requirements, where N is the number of unknowns, through a fast and nonlossy fast Fourier transform based matrix compression technique, the adaptive integral method, previously developed for solving large-scale electromagnetic problems. Because of its ability of handling large problems with complex geometries, the developed solver may constitute an efficient and high fidelity numerical simulation tool for calculating acoustic field distributions in anatomically realistic models, e.g., in investigating acoustic energy transfer to the inner ear via nonairborne pathways in the human head. Examples of calculations of acoustic field distribution in a human head, which require solving linear systems of equations involving several million unknowns, are presented.     

Optical near-field microscopy of light focusing through a photonic crystal flat lens

We report here the direct observation by using a scanning near-field microscopy technique of the light focusing through a photonic crystal flat lens designed and fabricated to operate at optical frequencies. The lens is fabricated using a III-V semiconductor slab, and we directly visualize the propagation of the electromagnetic waves by using a scanning near-field optical microscope. We directly evidence spatially, as well as spectrally, the focusing operating regime of the lens. At last, in light of the experimental scanning near-field optical microscope pictures, we discuss the lens ability to focus light at a subwavelength scale.     

A novel CN–ICCG–FDTD algorithm research of plasma reflection and transmission characteristics of electromagnetic wave

The advantage of the ICCG method for solving large sparse matrix is taken in the CN–FDTD equation solving. The CN–ICCG–FDTD can accelerate iteration for numerical calculation, and reduce memory overhead. Maxwell's equation of the electromagnetic wave in dispersive medium plasma is deduced from the CN format after being differentiated in the time domain, the FDTD method become unconditionally stable. In this process, a large sparse matrix is generated. As for solving such matrix, ICCG method has sufficient advantages. In ICCG method, convergence is so fast and stable that it is quite easy for computer programming. In addition, data required by ICCG method are memorized in a computer in the format of the one-dimensional compression, which helps the computer memory save a large amount of capacity, especially when the issues are rather complex to cope with, this algorithm will display such strength. With these advantages, ICCG method is able to calculate the reflection coefficient and the transmission coefficient of electromagnetic waves in plasma as well as the phase angle of them. The calculation indicates that the time step required by CN–ICCG–FDTD method has eliminated the constraints of the CFL, thus shortening the required time, making the calculated result stable and accurate, and improving the efficiency of programming.     

Mathematical simulation of nonlinear waveguiding systems based on 2D photonic crystals

A mathematical model of nonlinear two-dimensional waveguiding systems based on photonic crystals was constructed using the finite-difference time-domain method (FDTD), the total field/scattered field (TF/SF) method, and the perfectly matched layer method (PML). The constructed numerical algorithm can also be used for investigation of other two-dimensional waveguiding systems. The results of simulating the elements of particular waveguiding systems are presented.     

Surface electromagnetic waves supported by nano conducting layers with inhomogeneities in the conductivity profile

Conducting interfaces and nano conducting layers can support surface electromagnetic waves. Uniform charge layers of non-zero thickness and their asymptotic behavior toward conducting interfaces of infinitely small thicknesses, where the thin charge layer is modeled via a surface conductivity σ _s , are already studied. Here, the possible effects of inhomogeneity in the conductivity profile of the thin conducting layers are investigated for the first time and a new approximate yet accurate enough analytical formulation for mode extraction in such structures is given. In order to rigorously analyze the structure and justify the proposed approximate formulation, the Galerkin’s method with Legendre polynomial basis functions is applied, i.e. the transverse electric field for the TE polarized surface waves and the transverse magnetic field for the TM polarized surface waves are each expanded in terms of Legendre polynomials and then each eigenmode; subjected to appropriate boundary conditions, is sought in the complete space spanned by Legendre basis functions. The proposed approximate solution is then proved to be accurate. In particular, sinusoidal fluctuations are introduced into formerly uniform conductivity profiles and it is numerically demonstrated that surface electromagnetic waves supported by nano conducting layers are not much sensitive to the very shape of conductivity profiles.     

Error control in the perfectly matched layer based multilevel fast multipole algorithm

The development of efficient algorithms to analyze complex electromagnetic structures is of topical interest. Application of these algorithms in commercial solvers requires rigorous error controllability. In this paper we focus on the perfectly matched layer based multilevel fast multipole algorithm (PML-MLFMA), a dedicated technique constructed to efficiently analyze large planar structures. More specifically the crux of the algorithm, viz. the pertinent layered medium Green functions, is under investigation. Therefore, particular attention is paid to the plane wave decomposition for 2-D homogeneous space Green functions in very lossy media, as needed in the PML-MLFMA. The result of the investigations is twofold. First, upper bounds expressing the required number of samples in the plane wave decomposition as a function of a preset accuracy are rigorously derived. These formulas can be used in 2-D homogeneous (lossy) media MLFMAs. Second, a more heuristic approach to control the error of the PML-MLFMA’s Green functions is presented. The theory is verified by means of several numerical experiments.     

水平层状横向同性地层中频率测深资料全参数快速迭代反演

提出水平层状横向同性地层中频率测深资料的全参数快速迭代反演算法,以便从测量资料中同时确定各个地层的横向、 纵向电导率及层界面深度.首先,利用水平层状介质中电磁场并矢Green函数在频率波数域中的解析解和Sommerfeld积分的 快速计算技术确定正演响应.然后,利用摄动理论和Fourier逆变换公式,研究建立一套快速求解全参数Fréchet导数的 有效算法,并利用规范化处理和奇异值分解技术给出迭代反演过程.最后,给出数值试验证明反演理论的有效性以及反演算法 的抗噪声能力.