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

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

2.5维模拟中时域有限差分近-远场变换方法

为进行从高功率微波源到天线辐射整个系统的2.5维模拟,在近场模拟的基础上对2.5维模拟中时域有限差分(FDTD)的近-远场变换方法进行了研究。先将2.5维近场模拟区的输出边界绕对称轴旋转一周,得到封闭的3维圆柱面;然后根据等效原理由该封闭面上的FDTD近场求得等效电磁流;再根据自由空间中电流及磁流的时域辐射场公式计算出电磁流的远区辐射场。给出了2.5维模拟中TM01模式及TE01模式激励下圆锥喇叭远区辐射算例,所得计算结果与3维全电磁模拟软件CST的计算结果符合较好。     

2.5维模拟中时域有限差分近-远场变换方法

 为进行从高功率微波源到天线辐射整个系统的2.5维模拟,在近场模拟的基础上对2.5维模拟中时域有限差分(FDTD)的近-远场变换方法进行了研究。先将2.5维近场模拟区的输出边界绕对称轴旋转一周,得到封闭的3维圆柱面;然后根据等效原理由该封闭面上的FDTD近场求得等效电磁流;再根据自由空间中电流及磁流的时域辐射场公式计算出电磁流的远区辐射场。给出了2.5维模拟中TM01模式及TE01模式激励下圆锥喇叭远区辐射算例,所得计算结果与3维全电磁模拟软件CST的计算结果符合较好。     

截断各向异性等离子体的修正各向异性完全匹配层吸收边界

基于磁化等离子体本构方程,提出了一种截断各向异性色散介质的修正的各向异性完全匹配层(M-UPML)吸收边界算法。通过等效相对介电常数,将UPML推广到截断各向异性等离子体介质的情形,并推导了其时域有限差分方法(FDTD)迭代式。用该方法计算了半空间磁化与非磁化等离子体的反射系数,计算结果与解析解相一致,表明该吸收边界具有良好的吸收效果。     

利用时域有限差分法模拟金纳米球及球壳的光学特性

利用时域有限差分方法研究了金纳米球、金纳米球壳及多层球壳的消光特性及电场分布.结果表明:金纳米颗粒的几何参量对消光峰的位置有显著影响.随着SiO2核心半径的增大,金纳米壳的消光峰显著红移.随着金核心半径的增大,gold-silica-gold多层球壳消光谱的低能峰显著红移,而高能峰微弱蓝移.     

截断各向异性等离子体的修正各向异性完全匹配层吸收边界

 基于磁化等离子体本构方程,提出了一种截断各向异性色散介质的修正的各向异性完全匹配层(M-UPML)吸收边界算法。通过等效相对介电常数,将UPML推广到截断各向异性等离子体介质的情形,并推导了其时域有限差分方法(FDTD)迭代式。用该方法计算了半空间磁化与非磁化等离子体的反射系数,计算结果与解析解相一致,表明该吸收边界具有良好的吸收效果。     

色散介质电磁特性时域有限差分分析的Newmark方法

根据Debye模型、Drude模型和Lorentz模型3种常见色散介质模型频域极化率的特点,利用频域到时域的转换关系jω→?/?t,将极化矢量P与电场强度E的频域关系转换成时域内关于P的二阶微分方程,其对3种色散介质模型皆适用,具有统一的形式.然后采用相比于中心差分具有更高精度的Newmark两步算法(Newmark-β-γ法)求解该方程,进而得到E→P的递推公式,再结合本构关系得到D→E的时域递推式.实现了色散介质电磁场量的时域有限差分迭代计算.数值计算结果表明该方法是适用于3种色散介质模型的通用算法,并且相比于移位算子时域有限差分方法等以中心差分为基础的离散方案具有更高的计算精度.     

Investigation of three-pulse photon echo in thick crystal using finite-difference time-domain method

This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo's amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm~(3+):YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp(α'L),which is previously employed to describe the relationship between echo's efficiency and thickness,should be modified as 1.3 · 0.09 exp(2.4 ·α'L) with the propagation of echo considered.     

Exact boundary condition for time-dependent wave equation based on boundary integral

An exact non-reflecting boundary conditions based on a boundary integral equation or a modified Kirchhoff-type formula is derived for exterior three-dimensional wave equations. The Kirchhoff-type non-reflecting boundary condition is originally proposed by L. Ting and M.J. Miksis [J. Acoust. Soc. Am. 80 (1986) 1825] and numerically tested by D. Givoli and D. Cohen [J. Comput. Phys. 117 (1995) 102] for a spherically symmetric problem. The computational advantage of Ting–Miksis boundary condition is that its temporal non-locality is limited to a fixed amount of past information. However, a long-time instability is exhibited in testing numerical solutions by using a standard non-dissipative finite-difference scheme. The main purpose of this work is to present a new exact boundary condition and to eliminate the long-time instability. The proposed exact boundary condition can be considered as a limit case of Ting–Miksis boundary condition when the two artificial boundaries used in their method approach each other. Our boundary condition is actually a boundary integral equation on a single artificial boundary for wave equations, which is to be solved in conjunction with the interior wave equation. The new boundary condition needs only one artificial boundary, which can be of any shape, i.e., sphere, cubic surface, etc. It keeps all merits of the original Kirchhoff boundary condition such as restricting the temporal non-locality, free of numerical evaluation of any special functions and so on. Numerical approximation to the artificial boundary condition on cubic surface is derived and three-dimensional numerical tests are carried out on the cubic computational domain.     

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

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

A curved no-slip boundary condition for the lattice Boltzmann method

We present a generalization of the no-slip boundary condition by Lätt et al. [J. Lätt, B. Chopard, O. Malaspinas, M. Deville, A. Michler, Straight velocity boundaries in the lattice Boltzmann method, Physical Review E 77 (5) (2008) 056703] from straight to curved geometries for the lattice Boltzmann Bhatnager–Gross–Krook method (LBGK). The boundary condition is based on a reconstruction of the populations from the density, velocity and rate of strain. For curved boundaries, the reconstruction reduces the question of accuracy to a technical issue of interpolation. We present a method of interpolation allowing a very accurate representation of the curved boundary. The resulting boundary condition is verified for three different test cases: Taylor–Couette flow in-between rotating cylinders, laminar flow around a cylinder and flow past an impulsively started cylinder, demonstrating its second order accuracy and low error constant. The present boundary is stable for relaxation frequencies close to two.     

Optical simulation of in-plane-switching blue phase liquid crystal display using the finite-difference time-domain method

The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell's equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.     

Modeling of elastic wave propagation on a curved free surface using an improved finite-difference algorithm

Modelling and understanding the effect of elastic wave propagation along a curved free surface has been one of the important issues in seismic exploration[1—3], earth seis-mology[4], and non-destructive ultrasonic detection[5]. Several approaches have beenproposed for simulating wave propagation in heterogeneous media with a topographic stress-release boundary. These include finite-element methods (FEM), boundary element methods (BEM), finite-difference methods (FDM), pseudo-spectral metho…     

Application of the thermodynamics of curved boundary layers to the scaled particle theory of hard-sphere fluids

The thermodynamics of curved boundary layers, with the assumption that the distance between the surface of a fluid cavity and its surface of tension is a quadratic function of the cavity radius, is applied to the exact thermo-dynamic expression forG, the central function of scaled particle theory. The coefficients in the quadratic representation are determined so as to satisfyall five of the known exact conditions onG valid for cavity radii between one-half and one molecular diameter. The results of the calculation are displayed as the hard-sphere equation of state, the boundary tension associated with the surface of tension, and the distance between the cavity surface and the surface of tension. Although the hard-sphere equation of state obtained by this method using all five conditions onG is more accurate than in the case where only two or three conditions are used, the original scaled particle theory, in whichG itself was represented simply by a quadratic function of inverse powers of cavity radius, still yields the more accurate equation of state. Nevertheless, the present approach limits approximations to the distance between the cavity surface and the surface of tension, a small quantity in itself. The path to a still more improved theory remains well defined, contingent upon the discovery of additional exact conditions, and does not depend, as in the original scaled particle theory, on a form forG arrived at in a somewhat ad hoc manner.Research supported under NSF Grant #GP-12408.     

磁化铁氧体材料电磁散射递推卷积-时域有限差分方法分析

根据递推卷积原理，将磁化铁氧体材料的频域磁导率过渡到时域，通过引入时域复数磁化率张量和时域复数磁感应强度矢量，提出了磁化铁氧体材料电磁散射的三维时域有限差分方法.为了验证该方法，用它计算了磁化铁氧体球的后向雷达散射截面，与文献结果符合很好.计算结果表明，该方法是分析磁化铁氧体材料电磁散射一种可行的方法. 关键词： 递推卷积 磁化铁氧体 电磁散射 时域有限差分方法     

一种滑移区气体流动的格子Boltzmann曲边界处理新格式

针对滑移区复杂气-固边界存在速度滑移现象,提出了一种基于格子Boltzmann方法的非平衡态外推与有限差分相结合的曲边界处理新格式.该格式具有可考虑实际物理边界与网格线偏移量的优势,较传统half-way DBB(diffusive bounce-back)格式更能准确反映实际边界情况,同时还可获取壁面处气体宏观量及其法向梯度等信息.采用本文所提曲边界处理格式模拟分析了滑移区气体平直/倾斜微通道Poiseuille流、微圆柱绕流和同心微圆柱面旋转Couette流问题.研究结果表明,采用曲边界处理新格式所得结果与理论值以及文献结果符合良好,适用于滑移区气体流动的复杂边界处理,且比half-way DBB格式具有更高的精度,较修正DBB格式具有更好的适应性.     

一种计算多导体传输线电磁脉冲响应的时域方法

基于波形松弛技术, 提出一种计算外界电磁脉冲激励下理想大地上无损多导体传输线瞬态响应的时域迭代方法。首先利用波形松弛技术对复频域内多导体传输线的电报方程进行解耦, 其中相邻导线的耦合作用等效为线上的分布源, 从而使电报方程转换为一系列关于独立导线的解耦方程组; 然后将复频域内传输线的解耦方程转换到时域, 根据时域方程建立相应的等效电路; 最后利用电路仿真软件PSCAD计算电磁脉冲激励下多导体传输线的瞬态响应。本文时域方法的计算结果与时域有限差分(FDTD) 法计算的结果进行对比, 证实了该时域方法的有效性和准确性, 这为工程和科研人员快速评估、分析电磁脉冲激励下多导体传输线的瞬态响应问题提供了一种可靠方法。     

基于直接离散方式的磁化铁氧体材料电磁散射的时域有限差分方法分析

根据磁化铁氧体材料的磁场强度和磁感应强度之间色散的本构关系式,利用时间微分算子/t和jω的时域和频域对应关系,将磁化铁氧体材料频域的本构关系转化为时域的本构关系,然后将时间微分算子/t在时域采用直接离散的方式,得到磁场强度和磁感应强度的时域有限差分迭代式.数值结果表明,该方法易于实现,简单可行,并节约内存. 关键词： 电磁散射 磁化铁氧体 时域有限差分方法 直接离散方法     

Bandgap characteristics of 2D plasma photonic crystal with oblique incidence:TM case

A novel periodic boundary condition (PBC), that is the constant transverse wavenumber (CTW) method, is introduced to solve the time delay in the transverse plane with oblique incidence. Based on the novel PBC, the FDTD/PBC algorithm is proposed to study periodic structure consisting of plasma and vacuum. Then the reflection coefficient for the plasma slab from the FDTD/PBC algorithm is compared with the analytic results to show the validity of our technique. Finally, the reflection coefficients for the plasma photonic crystals are calculated using the FDTD/PBC algorithm to study the variation of bandgap characteristics with the incident angle and the plasma parameters. Thus it has provided the guiding sense for the actual manufacturing plasma photonic crystal.     

时域混合算法在一维海面与舰船目标复合电磁散射中的应用

采用时域积分方程(TDIE)与时域基尔霍夫近似(TDKA)的混合算法研究粗糙海面与舰船目标的复合瞬态电磁散射.该方法将舰船目标及其近邻海面划分为TDIE区域,用TDIE方法精确求解;将剩余电大尺寸的粗糙海面划分为TDKA区域,采用高效的TDKA电流近似求解.通过混合算法和传统TDIE算法结果的对比,表明TDIE-TDKA混合算法能保证计算的精度,同时具有较高的计算效率.最后,讨论了海面上方有无目标、海面上方风速、电磁脉冲入射角、舰船目标尺寸、吃水深度对后向散射磁场的影响.