斜入射分层线性各向异性等离子体电磁散射时域有限差分方法分析

提出了斜入射分层线性各向异性等离子体电磁散射的时域有限差分(FDTD)方法,通过将二维麦克斯韦方程等价地转换为一维麦克斯韦方程,避免了用二维时域有限差分方法分析该散射问题,极大地提高了计算效率.分析推导了TEz和TMz波斜入射线性分层各向异性等离子体电磁散射的FDTD方法,然后通过该方法计算不同入射角的各向异性等离子板的电磁波反射系数,并与其解析解进行比较,结果表明该方法的准确性和有效性.最后,将该算法应用于计算涂覆分层各向异性等离子体金属板在不同入射角下的反射系数,分析了不同入射角对反射系数的影响.     

用MEI方法研究脉冲线源辐射问题

 脉冲线电流源的辐射问题是研究高功率微波辐射、传播和散射问题的基础。用时域有限差分方法结合时域MEI吸收边界条件来对该问题进行求解，通过与不同种类边界条件配合得出的数值解的比较可以证明，利用时域MEI方法的一阶吸收边界条件对线源辐射问题进行处理时，可以使截断边界离源更近，从而减少了计算量，并达到提高计算精度的目的。     

Linear Full Decoupling,Velocity Correction Method for Unsteady Thermally Coupled Incompressible Magneto-Hydrodynamic Equations

We propose and analyze an effective decoupling algorithm for unsteady thermally coupled magneto-hydrodynamic equations in this paper. The proposed method is a first-order velocity correction projection algorithms in time marching, including standard velocity correction and rotation velocity correction, which can completely decouple all variables in the model. Meanwhile, the schemes are not only linear and only need to solve a series of linear partial differential equations with constant coefficients at each time step, but also the standard velocity correction algorithm can produce the Neumann boundary condition for the pressure field, but the rotational velocity correction algorithm can produce the consistent boundary which improve the accuracy of the pressure field. Thus, improving our computational efficiency. Then, we give the energy stability of the algorithms and give a detailed proofs. The key idea to establish the stability results of the rotation velocity correction algorithm is to transform the rotation term into a telescopic symmetric form by means of the Gauge–Uzawa formula. Finally, numerical experiments show that the rotation velocity correction projection algorithm is efficient to solve the thermally coupled magneto-hydrodynamic equations.     

三维周期结构弱无条件稳定时域有限差分算法

从麦克斯韦方程出发,导出了适用于周期结构的三维弱无条件稳定时域有限差分（FDTD）算法的迭代式,在理论上证明了其稳定性条件,其稳定性条件比普通FDTD的稳定性条件要宽松,并将周期边界条件应用到迭代式中,得到了可直接编程迭代的方程,给出了对其特有的一类非三对角阵的解决办法。最后通过算例将计算结果与常规FDTD及ADI-FDTD计算结果比较,验证了算法的精确性和有效性。     

Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves

This paper is devoted to time domain numerical solutions of two-dimensional (2D) material interface problems governed by the transverse magnetic (TM) and transverse electric (TE) Maxwell's equations with discontinuous electromagnetic solutions. Due to the discontinuity in wave solutions across the interface, the usual numerical methods will converge slowly or even fail to converge. This calls for the development of advanced interface treatments for popular Maxwell solvers. We will investigate such interface treatments by considering two typical Maxwell solvers – one based on collocation formulation and the other based on Galerkin formulation. To restore the accuracy reduction of the collocation finite-difference time-domain (FDTD) algorithm near an interface, the physical jump conditions relating discontinuous wave solutions on both sides of the interface must be rigorously enforced. For this purpose, a novel matched interface and boundary (MIB) scheme is proposed in this work, in which new jump conditions are derived so that the discontinuous and staggered features of electric and magnetic field components can be accommodated. The resulting MIB time-domain (MIBTD) scheme satisfies the jump conditions locally and suppresses the staircase approximation errors completely over the Yee lattices. In the discontinuous Galerkin time-domain (DGTD) algorithm – a popular Galerkin Maxwell solver, a proper numerical flux can be designed to accurately capture the jumps in the electromagnetic waves across the interface and automatically preserves the discontinuity in the explicit time integration. The DGTD solution to Maxwell interface problems is explored in this work, by considering a nodal based high order discontinuous Galerkin method. In benchmark TM and TE tests with analytical solutions, both MIBTD and DGTD schemes achieve the second order of accuracy in solving circular interfaces. In comparison, the numerical convergence of the MIBTD method is slightly more uniform, while the DGTD method is more flexible and robust.     

周期对称结构R-FDTD方法及其减元优化

论述减缩时域有限差分方法(R-FDTD)中暂存场分量边值补充计算的必要性,提出周期对称结构R-FDTD方法,基于对称关系和周期边界条件(PBC),给出需要补充计算的暂存分量表达式.利用对称性,将计算空间缩减为原来的1/4,对称面外侧场分量由对称关系得到,1/4空间周期对称结构R-FDTD可更进一步将计算区域的内存使用量降为FDTD算法的1/6,且不影响计算精度.计算无限大钢筋网和钢筋混凝土墙壁的电磁脉冲响应,结果与FDTD的计算结果吻合.改进的算法在内存使用和计算时间上具有明显的优势.     

光子晶体传输特性的时域精细积分法分析

将精细积分法应用于时域有限差分法中,提出了一种求解光子晶体传输特性的时域精细积分法,并对其计算精度及稳定性进行了分析.从一阶麦克斯韦方程出发,在空间上采用Yee元胞进行差分离散,结合吸收边界条件及激励源表达式将方程整理为标准的一阶常微分方程组形式.通过时间步长的精细划分和指数矩阵的加法定理,在时间上利用精细积分法对齐次微分方程进行积分求解,并结合激励向量的特解得到空间离散的场分量,最终通过傅里叶变换求得方程的解.利用时域精细积分法对光子晶体进行了实例计算,并将其结果分别与时域有限差分法和四阶龙格库塔法在精度、稳定性等方面进行了比较,结果表明时域精细积分法具有更高的计算精度,并且克服了时域有限差分法以及四阶龙格库塔法在计算稳定性上对时间步长的限制.提出的方法具有精确、稳定的特点,为光子晶体传输特性的研究提供了一种新的有效的分析方法.     

Reducing the calculation workload of the Green function for electromagnetic scattering in a Schwarzschild gravitational field

When the finite difference time domain(FDTD) method is used to solve electromagnetic scattering problems in Schwarzschild space-time, the Green functions linking source/observer to every surface element on connection/output boundary must be calculated.When the scatterer is electrically extended, a huge amount of calculation is required due to a large number of surface elements on the connection/output boundary.In this paper, a method for reducing the calculation workload of Green function is proposed.The Taylor approximation is applied for the calculation of Green function.New transport equations are deduced.The numerical results verify the effectiveness of this method.     

A Hybrid Technique Combining FDTD and Series Solution for Near-to-Far-Field Transformation

A new near-to-far-field transformation that combines the FDTD and the series solution is presented. In this method, near-field data is calculated by FDTD and far-field data is also expressed by series solution. The boundary condition is applied to a virtual surface, which encloses arbitrarily shaped scatterers. In order to verify this method, scattering patterns from this transformation are compared with scattering patterns from exact solution. Numerical results show that errors in the far-field data from this method is less than surface equivalence theorem. Thus, this method is valid for obtaining the far-field data.     

On multidimensional ultrasonic scattering in an inhomogeneous elastic background.

This work is concerned with the modeling of elastic wave scattering by solid or fluid-filled objects embedded in an inhomogeneous elastic background. The medium is probed by a monochromatic force and the scattered field is computed (forward problem) or observed (inverse problem) at some known receiver locations. Based on vector integral equations for elastic scattering, a general framework is developed, independent of both the problem geometry and the transmitter-receiver characteristics. This framework encompasses both forward and inverse modeling. In the forward model, a Born approximation for an inhomogeneous background is applied to obtain a closed form expression for the scattered field. In the inverse model, this approximation is also invoked to linearize for the multiparameter characteristic of the object. Finally, an iterative inversion scheme alternating forward and inverse modeling is proposed to improve the resolution and accuracy of the reconstruction algorithm.     

New method for the theory of hadron scattering

A new method is suggested for analyzing the scattering of particles whose interaction is mediated by the exchange of field quanta. The conventional approach uses the Lehmann, Symanzik, and Zimmerman (LSZ) method for formulating the asymptotic boundary condition constraints. Then after making the appropriate renormalization subtractions one finds Bethe-Salpeter (BS) equations for the transition amplitudes. The covariant form of the BS equations makes the construction of solutions arduous. Approximate reductions to noncovariant form are generally employed. Our method is based on an alternative to the LSZ technique and leads directly to a set of coupled noncovariant integral equations for the transition operators. No approximation is required to achieve this result. Our equations are similar to the Lippmann-Schwinger equations for many-body scattering. The complexity of these equations is governed by the maximum number of field quanta permitted to be present at any instant.     

三维燃烧介质和壁面温度的非接触联合重建研究

温度分布在线实时测量对于燃烧过程优化和污染物控制具有重要意义, 针对以往非接触三维温度分布重建过程的耗时性问题和忽略壁面辐射的不足, 本文提出了一种新的离散重建模型, 用于三维吸收、 发射和散射性高温燃烧介质以及壁面温度的快速联合非接触测量. 该模型以四个CCD(Charge Coupled Device) 为测量传感器, 通过构建辐射逆问题求解方程, 从CCD输出的辐射投影图像重建温度分布. 介质中不同投影方向内的辐射传递过程通过离散传递法来描述, 介质的散射和壁面反射则通过离散坐标法来近似. 离散后计算局部辐射强度的病态方程通过最小二乘余量法来求解, 论文对其计算速度进行了优化. 通过非对称温度分布测量算例分析了该模型的有效性, 讨论了测量噪音、 介质和壁面辐射特性对重建精度的影响, 并与其他方法对比分析了模型的重建速度. 计算结果表明本文提出的离散模型可以有效地用于大型高温燃烧介质和壁面温度分布的联合非接触测量. 即使在有噪声的情况下, 该模型也能获得准确的测量结果, 与其他计算方法相比, 采用改进的最小二乘余量法, 能有效地提高温度分布的重建计算速度.     

一种新时域交替隐式差分算法在散射问题中的应用

提出一种新时域交替隐式有限差分(ADI-FDTD)算法格式. 传统ADI-FDTD算法的 场量步进方程涉及周围若干网格的较多场量，导致两个区域的步进方程处理较困难：一个是 邻近完全匹配层(PML)和散射场交界区，另一个是邻近连接边界区. 特别是后者，考虑入射 波影响需对场量所在区域判断，根据不同情况对原有方程进行修正，一维和二维散射问题 相对简单，可三维问题修正式有数十种之多而几乎无法完成. 本方法基于分裂场形式的ADI- FDTD技术，使得散射场区和PML吸收层区的表达形式完全一致，从而忽略两者差别.另 关键词： 时域交替隐式有限差分算法 电磁散射     

A shifted Jacobi collocation algorithm for wave type equations with non-local conservation conditions

In this paper, we propose an efficient spectral collocation algorithm to solve numerically wave type equations subject to initial, boundary and non-local conservation conditions. The shifted Jacobi pseudospectral approximation is investigated for the discretization of the spatial variable of such equations. It possesses spectral accuracy in the spatial variable. The shifted Jacobi-Gauss-Lobatto (SJ-GL) quadrature rule is established for treating the non-local conservation conditions, and then the problem with its initial and non-local boundary conditions are reduced to a system of second-order ordinary differential equations in temporal variable. This system is solved by two-stage forth-order A-stable implicit RK scheme. Five numerical examples with comparisons are given. The computational results demonstrate that the proposed algorithm is more accurate than finite difference method, method of lines and spline collocation approach     

基于表面阻抗边界条件的等离子体薄涂层电磁散射的时域有限差分分析

基于表面阻抗边界条件时域有限差分(FDTD)方法研究了一维斜入射情况下非磁化等离子体薄涂层涂敷金属材料的电磁散射特性, 该方法忽略对薄层背景材料进行网格剖分, 大大减少了计算量. 首先推导了理想导体涂敷等离子体薄涂层的表面阻抗频域表达式, 然后代入边界条件并变换到时域, 再用分段线性递推卷积方法将时域表达式离散得到FDTD迭代式. 编程计算了垂直及斜入射情形下的平行极化和垂直极化反射系数, 通过验证算例与解析解对比, 结果表明该方法的准确性和有效性. 最后利用该方法分析了不同入射角对反射系数的影响.     

Investigation of composite electromagnetic scattering from ship-like target on the randomly rough sea surface using FDTD method

Composite electromagnetic scattering from a two-dimensional (2D) ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain (FDTD) method. A uniaxial perfectly matched layer is adopted for truncation of FDTD lattices. The FDTD updated equations can be used for the total computation domain by choosing the uniaxial parameters properly. To validate the proposed numerical technique, a 2D infinitely long cylinder over the sea surface is taken into account first. The variation of angular distribution of the scattering changing with incident angle is calculated. The results show good agreement with the conventional moment method. Finally, the influence of the incident angle, the polarization, and the size of the ship-like target on the composite scattering coefficient is discussed in detail.     

一种适用于任意阶空间差分时域有限差分方法的色散介质通用吸收边界条件算法

基于伸展坐标完全匹配层方程和辅助微分方程 方法, 给出了一种在时域有限差分(FDTD) 计算中适用于常见色散介质的通用边界条件算法. 该算法适用于任意阶的FDTD空间差分, 并且由于所采用的D-H方程独立于计算区域, 所以可以被用于截断任意电介质. 数值试验结果表明, 与卷积完全匹配层 算法相比较, 所提出的吸收边界条件算法不仅通用性强、 计算复杂度低、 计算时间短, 并且吸收效果有明显的提高.     

基于交替隐式有限差分法的快速早期乳腺癌时域微波断层成像

采用时域微波断层成像技术进行早期乳腺癌检测能够准确地获得乳房的电参数分布,具有明确的物理解释和医学诊断价值.临床应用讲究即时性,为了提高检测的速度,本文将交替隐式有限差分法应用到乳腺癌检测中,基于正反演时间步进成像算法进行成像分析,结果显示在保证精度的前提下,采用交替隐式有限差分法的成像时间可缩短为传统时域有限差分法的23%,提高了微波断层成像技术的临床可应用性.     

随机表面散射光场的格林函数法与基尔霍夫近似的比较

从简谐光波满足的亥姆霍兹方程出发，将由格林定理得到的介质分界面上的积分方程转化为以表面上的光波及其导数为未知量的线性方程组，并对其进行数值求解，实现了光场的数值计算. 同时，由透射光场的格林函数积分得出了基尔霍夫近似下光场的表达式. 通过类比推导夫琅和费面上散斑场自相关函数的方法，提出了产生随机表面及其导数的傅里叶变换方法. 在此基础上，对采用基尔霍夫近似进行自仿射分形随机表面的散射光场数值计算的精确程度进行了研究. 发现在随机表面粗糙度比较小时，基尔霍夫近似的精度比较高；在粗糙度相同的情况下，表面的分形 关键词： 格林函数积分 基尔霍夫近似 自仿射分形随机表面