Ground-penetrating radar for soil and underground pipelines using the forward modeling simulation method

Ground-penetrating radar was simulated using the FDTD method with uniaxial anisotropic perfectly matched layers as absorbing boundary conditions. To study underground pipeline detection, the electromagnetic response properties of underground pipelines, under known soil conditions, were determined. Pipelines at different depths and with different diameters and different compositions were simulated and researched for their response. The simulation results were consistent with research and analysis, i.e., the deviation was in an acceptable range. This research is of great value to guide the detection of pipelines in cities and to explain the detection of materials underground.     

A hybrid modeling for the theoretical analysis of reflections in a multimode-interference coupler based on silicon-on-insulator nanowires

A hybrid method combining the three-dimensional (3D) beam propagation method (BPM) and a 3D finite-difference time-domain (FDTD) method is presented for the simulation of multimode-interference (MMI) couplers based on silicon-on-insulator (SOI) nanowires. This hybrid method makes it possible to have a robust analysis for the reflection properties while avoiding a time-consumed computation with a simplex 3D-FDTD simulation. The reflected power is sensitive to both the length and the width of the MMI coupler and thus one should choose the optimal values for good self-imaging. Furthermore, the reflection could be reduced by introducing tapers between the MMI section and the access waveguides.     

Analysis of femtosecond optical pulse propagation in one-dimensional nonlinear photonic crystals using finite-difference time-domain method

In this paper, we have used the finite-difference time-domain (FDTD) method to analyze the optical pulse propagation in a nonlinear, one-dimensional photonic crystal (1DPC). Hyperbolic secant pulses with various carrier wavelengths are utilized in this study. In a nonlinear regime, a 1DPC introduces a photonic band-gap whose central wavelength and width depend on the input pulse intensity. In the present work, three different cases are considered. These correspond to the carrier wavelengths of the incident pulses being out of, near to, and partially in the band-gap. For each case, the effect of nonlinearity on pulse propagation is investigated. Also, we have analyzed the two-frequency regime, in which each of the two pulses has a different carrier frequency (wavelength). This kind of study can be done directly with FDTD without any further computational burden but it is somewhat complicated using nonlinear coupled-mode equations (NLCME) and nonlinear Schrödinger equation (NLSE), which require separate treatments for each carrier wavelength.     

Ultrasound propagation in trabecular bone: A numerical study of the influence of microcracks

The accumulation of microdamage in trabecular bone tissue is suspected of being a predictive indicator of osteoporosis diagnosis. To quantify this microdamage, the Dynamic AcoustoElastic Testing (DAET) method measures the Time Of Flight (TOF) and amplitude variations of transmitted ultrasound (US) pulses, while the bone sample is submitted to a low frequency sinusoidal hydrostatic pressure (opening/closing of microcracks). However, DAET is both sensitive to viscoelastic properties changes and microcracks density. To verify the microcracks density contribution on DAET results, a numerical approach is proposed. Multliple configurations of microdamaged trabecular bone-tissue-like mesh have been simulated. A 2D pseudo-spectral time domain numerical model was then developed to simulate linear wave propagation in heterogeneous solids. The influence of the microcracks number and orientation on the US TOF was particularly investigated. Results are discussed and compared with experimental data obtained from DAET measurements in trabecular bone samples.     

时域有限差分法在建筑声学中的应用及前景

介绍了声波方程的基本差分格式及稳定条件、数值色散、吸收边界条件等数值计算理论，例举了前人用时域有限差分法对噪声传播过程的模拟和室内声学中座椅吸声低谷效应模拟的模拟结果。本文指出，由于时域有限差分法的特点使其具有在模拟脉冲响应方而的特别优势，因而，应用这一技术研究厅堂的声学特性，尤其是低频特性，将有重要意义。     

Design of ultrasonic wedge transducer

Cones and wedges inserted between an ultrasonic transducer and the specimen provide the transducer (circular or rectangular shape) with enhanced capability for point or line contact with the specimen. Such an arrangement is useful in that the transducer can be used for transmitting to and receiving from a point (or line) source, and that it can eliminate the undesirable aperture effect that makes the transducer blind to waves traveling in certain directions and those of certain frequencies. In this paper, a comprehensive numerical analysis based on a wave propagation model is carried out for the study of characteristics and parameters of cones and wedges influencing their performance. We study the effect of the dimensions, shape and aperture on the frequency response and the angle of incidence of the wave. For computational accuracy and efficiency, the boundary element method is used in the analysis.     

Improving the performance of a photonic crystal ring-resonator-based channel drop filter using particle swarm optimization method

The optical channel drop filters (CDFs) based on photonic crystals (PCs) are believed to be the essential components for compact photonic integrated circuits and WDM systems. One of the promising designs for a PC-CDF is based on the photonic crystal ring resonators (PCRRs). In this study, different parameters of a PCRR based CDF is optimized by using the particle swarm optimization (PSO) method in conjunction with the numerical method of two-dimensional (2D) finite-difference time-domain (FDTD) in a square lattice dielectric-rod 2D-PC. Hence a PC-CDF with ideal drop efficiency, large free spectral range (FSR) and acceptable quality factor will be proposed.     

Analysis on the efficiency of parallel FDTD method and its application in two-dimensional photonic crystal

A new kind of parallel finite-difference time-domain (FDTD) method is used to calculate and analyze the transmission properties of photonic crystal. The effectiveness and the accuracy of parallel FDTD method were proved by the typical two-dimensional sine wave. The result shows that parallel FDTD method, which can save time effectively in electromagnetic problems with large size and long time, performs as accurate as the series one. This method, whose results shows that the photonic band gap normalized frequency moves to the low-frequency direction with the enlargement of the section area of the square cylinder and that the band gap gets thinner, is applied to simulate and analyze the two-dimensional photonic crystal and gives calculating method and analyzing accordance to calculate and analyze large-sized photonic crystal structure.     

Efficient inverse radiation analysis in a cylindrical geometry using a combined method of hybrid genetic algorithm and finite-difference Newton method

An inverse radiation problem was considered to estimate boundary conditions such as temperature distribution and emissivity in axisymmetric absorbing, emitting, and scattering medium, given the measured incident radiative heat fluxes. The finite-volume method was employed to solve a direct radiative transfer equation for a two-dimensional axisymmetric geometry. Various parameter estimators, such as conjugate-gradient method, hybrid genetic algorithm, and finite-difference Newton method, were employed to solve the inverse problems, while discussing their performances in terms of estimation accuracy and computational efficiency. Based on this, we proposed, as a best inverse analysis tool, a new combined method that adopted the hybrid genetic algorithm as an initial value selector and used the finite-difference Newton method as a parameter estimator.     

Three-dimensional FDTD method for optical pulse propagation analysis in microstructured optical fibers

We propose a three-dimensional (3D) finite-difference time-domain (FDTD) method to analyze the pulse propagation characteristics in microstructured optical fibers (MOFs). The computation domain size is greatly reduced by adopting the technique of moving problem space. The propagating pulse is virtually held in the buffer cell of the problem space as simulation continues. This method is capable to investigate the temporal evolution of the propagating pulse. Spectral information can be obtained by Fourier analysis. As an example, the influence of the kerr nonlinearity on the optical pulse propagation in a Lorentz dispersive MOF is demonstrated. The model is also used to simulate the nonlinear interactions between the pump spectral broadening and third harmonic generations in a highly nonlinear fused silica nanowire with good agreement with the generalized nonlinear envelop equation (GNEE) model.     

Dual-beam-reflection phenomenon due to leaky modes in a photonic band gap

We show a dual-beam-reflection phenomenon for a Gaussian beam illuminating a hybrid structure of a dielectric waveguide and photonic crystal (WG-PC) inside the photonic band gap by numerical simulations. One reflection beam has a giant negative lateral shift, but the other has a positive lateral shift. The finite-difference time-domain (FDTD) simulations show that this phenomenon has a time delay effect and comes from the leaky surface mode of the hybrid structure. Field profile of the leaky mode demonstrates a strong localized stationary field in the higher dielectric medium. Furthermore, the maximum lateral shift is almost two times of the waist of the incident beam.     

Extension of the distributed point source method for ultrasonic field modeling

The distributed point source method (DPSM) was recently proposed for ultrasonic field modeling and other applications. This method uses distributed point sources, placed slightly behind transducer surface, to model the ultrasound field. The acoustic strength of each point source is obtained through matrix inversion that requires the number of target points on the transducer surface to be equal to the number of point sources. In this work, DPSM was extended and further developed to overcome the limitations of the original method and provide a solid mathematical explanation of the physical principle behind the method. With the extension, the acoustic strength of the point sources was calculated as the solution to the least squares minimization problem instead of using direct matrix inversion. As numerical examples, the ultrasound fields of circular and rectangular transducers were calculated using the extended and original DPSMs which were then systematically compared with the results calculated using the theoretical solution and the exact spatial impulse response method. The numerical results showed the extended method can model ultrasonic fields accurately without the scaling step required by the original method. The extended method has potential applications in ultrasonic field modeling, tissue characterization, nondestructive testing, and ultrasound system optimization.     

Design and Simulation of a Wide-Beam High-Power Parabolic Antenna with Dual-Source

For resolving the essential problem that the beam is too narrow in the application of the high-power UWB (Ultra-Wideband) antenna with single-source, a novel wide-beam high-power parabolic antenna with dual-source has been designed, and the edge of the paraboloid was cut in order to reduce the size of the antenna. Radiation properties of this proposed antenna are studied with finite-difference time-domain (FDTD) numerical method. The radiation properties of the proposed antenna can be improved by changing the structure of the parabolic antenna and adjusting the angle between the source and the axis of the antenna. The result shows that the beam-width has been improved significantly under the precondition of maintaining the peak power, and also the size of the antenna was minished by cutting the edge of the paraboloid.     

生物细胞亚显微结构对光散射特性的影响

构建细胞模型、"细胞基质"模型以及"细胞基质-细胞核"模型等3个模型,分别表征不同层次细胞结构,采用时域有限差分方法进行仿真计算,分别比较了3个模型的远场雷达散射截面在5个散射区的数值差异,分析了细胞器、细胞核和细胞基质等细胞亚显微结构对细胞电磁散射特性的影响.研究结果表明,生物细胞整体电磁散射特性主要由其细胞基质决定,细胞核的存在显著增强了细胞在20°≤θ≤40°和150°≤θ≤180°两个散射区的散射能力,随机分布的细胞器增强了细胞的侧向散射,削弱了细胞核在上述两个散射区的散射效果,使细胞散射波分布趋于均匀.     

屏蔽电缆对脉冲X射线响应的数值计算

 结合蒙特卡洛方法和时域有限差分(FDTD)方法,计算了电缆受脉冲X射线辐照时介质层内的运流电流密度,并以此为麦克斯韦方程的源,计算得到了电缆两端接匹配负载时的芯线响应电流。该方法综合考虑了电缆芯线、介质层和屏蔽层的沉积电荷对芯线响应电流的影响。计算结果表明：芯线响应电流大小与电缆受辐照长度成正比,电流由辐照中心向两边流走;源区越靠近中心位置,电流幅度越小,源的中心位置处,电流为零,源区存在静电场;源区外,电流大小相等,方向相反。最后,利用有限差分法计算得到的电场强度反推出了芯线电荷数,与蒙特卡洛方法计算的结果相比,FDTD方法计算的要低20%,该误差可能由将3维问题近似为1维问题所引起。     

一维光子晶体激光器中模场的空间分布及其对阈值的影响

采用有限时域差分（FDTD）法和传输矩阵法（TMM）计算了一维光子晶体微腔中模式的阈值和模场的空间分布,并用局域化长度和模面积描述了其空间分布特性.着重研究了模场的空间分布对阈值特性的影响.模拟结果显示,激光器的阈值依赖于模式的空间分布,局域化长度和模面积最小的模式具有最低阈值.通过调节包括缺陷宽度、光学厚度比、折射率比等介质参量,可以有效地调节模式的局域化长度和模面积,从而降低激光器的阈值.     

宽波束双馈源超宽带抛物面天线设计与仿真

 针对单馈源超宽带抛物面天线波束较窄的问题，设计了一种宽波束双馈源超宽带抛物面天线。该抛物面天线是由一个普通的抛物面从中间对称分式，并在中间以平板连接而成的。馈源采用改进型平面TEM喇叭天线，同轴垂直馈电方式。两个馈源分别置于抛物面的两个焦点上。并利用时域有限差分数值方法分析了该天线的辐射性能。通过改变抛物面的结构，调节双馈源与天线主轴的夹角等方法可以有效改善天线的辐射性能。结果表明，所设计的宽波束双馈源抛物面天线可以在远区更宽的空间上形成峰值功率比较均匀的电磁脉冲，在不降低峰值功率的前提下有效地增加了波束的宽度；在抛物面天线两边进行少量切割不会对辐射场产生明显影响，节省了空间。     

各向异性海底地层海洋可控源电磁响应三维积分方程法数值模拟

采用积分方程法对各向异性海底地层海洋可控源电磁（MCSEM）响应进行三维数值模拟.首先利用算子理论给出各向异性海底地层中的压缩积分方程,由于满足压缩映射条件,该积分方程在任意参数下总是迭代收敛的.然后为了提高计算效率,引入分区域多重网格准线性近似技术.通过具体算例验证了所述算法在计算精度与效率方面的有效性.最后利用该算法考察并分析了海底地层的各向异性对MCSEM三维响应特征的影响. 关键词： 海洋可控源电磁法 各向异性 三维模拟 积分方程法     

A hybrid spatial-spectral denoising method for infrared hyperspectral images using 2DPCA

The traditional noise reduction methods for 3-D infrared hyperspectral images typically operate independently in either the spatial or spectral domain, and such methods overlook the relationship between the two domains. To address this issue, we propose a hybrid spatial-spectral method in this paper to link both domains. First, principal component analysis and bivariate wavelet shrinkage are performed in the 2-D spatial domain. Second, 2-D principal component analysis transformation is conducted in the 1-D spectral domain to separate the basic components from detail ones. The energy distribution of noise is unaffected by orthogonal transformation; therefore, the signal-to-noise ratio of each component is used as a criterion to determine whether a component should be protected from over-denoising or denoised with certain 1-D denoising methods. This study implements the 1-D wavelet shrinking threshold method based on Stein’s unbiased risk estimator, and the quantitative results on publicly available datasets demonstrate that our method can improve denoising performance more effectively than other state-of-the-art methods can.