一维周期结构散射近-远场外推的两种方法

首先讨论了自由空间中一维周期结构近-远场外推的Floquet模方法,即对FDTD计算所得散射近场中输出边界上散射场进行级数展开,求出各阶Floquet模的复数幅值,再求得远区场。接着介绍了求解一维周期结构远区散射场的周期Green函数方法,即根据FDTD计算所得散射近场中输出边界上散射场,求得等效面电磁流后,再借助周期Green函数进行外推得到远区场。两种方法均仅用一个周期单元内的散射近场进行外推。计算结果验证了上述两种方法外推的有效性。     

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的计算结果符合较好。     

基于半空间FDTD的近远场外推方法: TE情形

研究基于半空间时域有限差分(FDTD)计算中TE情形的远区散射场计算方法.先将半空间Green函数代入势函数,推导远区纵向磁场的计算公式,对各空间内电磁流和对应的相位进行详细分析.再计算两个半空间模型的远区散射场,结果表明:该方法是正确和有效的.针对日本96式装甲车的计算发现,单站雷达探测系统位于战车正上方时具有最好的探测效果.     

基于时域有限差分法的液晶光学特性模拟

研究了模拟光在液晶中传播的时域有限差分(FDTD)法。该方法采用完全匹配层吸收边界条件来截断计算区域,通过独立计算总场区域与散射场区域的左连接边界和右连接边界上的入射场场值来实现非周期结构的平面波波源引入,并提出了分裂场一维辅助FDTD方法来计算含有非均匀各向异性介质的边界上的入射场场值。实现了平面波源入射下具有非周期结构的液晶器件光学特性计算。对开态下扭曲向列相液晶盒的仿真结果表明,30°入射时散射区透射率的最大泄漏误差小于3%,可知该方法对非周期、非均匀各向异性介质光学特性分析是有效的。     

垂直极化平行板有界波电磁脉冲模拟器辐射近场的快速估算方法

为快速估算出垂直极化平行板有界波电磁脉冲(EMP)模拟器的时域近场,将散射传递函数法应用于该类型模拟器近场的时域计算中,即对于给定的脉冲源,先寻找有效频谱范围能覆盖该源的高斯脉冲源,并应用时域有限差分(FDTD)方法计算该高斯脉冲源激励时模拟器中测试点场的时域响应,再利用傅里叶变换、系统的传递函数及傅里叶逆变换计算得到给定脉冲源激励时各测试点场的瞬态响应。所得计算结果与直接使用给定脉冲源激励时FDTD方法的计算结果符合较好。所述方法可用于同一模拟器在不同脉冲源激励时辐射近场的快速估算,能大大减少FDTD模拟计算的次数,尤其对于中大型模拟器能有效减少计算时间和内存。     

色散周期结构的辅助场时域有限差分法分析

利用辅助场时域有限差分(FDTD)方法分析色散周期结构的斜入射问题. 主要基于Floquet定理, 在FDTD迭代式中引入辅助场量, 结合变换的周期及吸收边界条件,解决了周期结构斜入射问题的电磁特性仿真. 进一步将辅助场FDTD方法引入至色散周期结构的模拟, 结合Z变换方法,给出了基于Drude色散模型的具体迭代公式. 数值计算结果表明所构建方法的有效性及普适应.     

分析水下目标的大比例变换总场散射场源FDTD方法

为解决探测水下目标的电磁散射问题,提出大比例变换的总场-散射场源时域有限差分（FDTD）方法.该方法包含两次FDTD计算：第一次计算采用细网格得到激励源周围的近场值;第二次计算采用粗网格得到远距离的电磁场值.两次FDTD计算通过总场-散射场边界建立联系.实现细粗网格的大比例变换,例如变换比例N=10,大大节省了计算时间,降低了计算内存的消耗,提高了计算效率.通过算例验证该方法的正确性和有效性.最后,计算水下岩层中存在异常体时的电磁响应,指出当岩层中异常体电导率不同时,接收点处电磁场的幅值和相位均不相同.     

三维随机粗糙面与目标复合电磁散射的FDTD方法

用时域有限差分方法(FDTD)研究三维周期性延拓的随机粗糙面与上方目标复合电磁散射.用周期性延拓消除数值计算中截取有限大小粗糙面产生的边缘效应,讨论一个周期单元粗糙面的边长与其相关长度之间的关系.给出在FDTD方法中向粗糙面加载入射波的方法,建立了粗糙面上单个三维目标的复合散射FDTD计算模型.数值结果给出粗糙面与目标散射的近场分布,应用近远场变换得到全方位散射角的双站散射系数.比较了三维与二维散射模型的区别.结果显示当粗糙面上放置目标时,其后向散射显著增强.     

大耦合孔同轴输出腔的3维解析分析

 在高功率微波器件中，通常采用大耦合孔的同轴微波输出腔，该腔为3维结构。采用场等效原理将腔体、耦合孔及其输出负载(行波边界)进行分区,每个区域中的场由界面上的磁流密度决定。采用格林函数积分法可得每个区域中的3维场,再由各个区域的场匹配方程求得腔体的谐振频率、特性阻抗、有载品质因数、模式分布等参数。为腔体的计算和设计提供一种计算模型。     

双间隙谐振腔的高频分析

从麦克斯韦方程组出发,运用电磁场的匹配条件,导出了双间隙谐振腔等效周期结构横磁模的色散关系及场分布表达式。然后通过数值计算,求出了双间隙谐振腔的谐振频率和场分布,并分析了谐振腔的结构尺寸对频率的影响,为C波段双间隙输出腔的设计提供了理论依据。     

晶间界面附近区域中的弹性场——各向异性理论

根据现代关于晶间界面的研究可知,晶界和相界多具有二维周期结构。由此出发,采用Fourier级数展开的表示,在给定了界面结构,在数学上即是形式地给定了边界条件的情况下,求出了两个各向异性物相间相界附近的弹性场。并以Ni中扭转晶界为例,作了应力场的数值计算。本文的方法和方程对于研究各种类型的晶界和相界均普遍可用,为处理这类问题提供了一种有效的简捷手段。 关键词：     

New boundary conformal field theories indexed by the simply laced Lie algebras

We consider the field theory of N massless bosons which are free except for an interaction localized on the boundary of their (1+1)-dimensional world. The boundary action is the sum of two pieces: a periodic potential and a coupling to a uniform abelian gauge field. Such models arise in open-string theory and dissipative quantum mechanics, and possibly in edge state tunneling in the fractional quantized Hall effect. We explicitly show that conformal invariance is unbroken for certain special choices of the gauge field and the periodic potential. These special cases are naturally indexed by semi-simple, simply laced Lie algebras. For each such algebra, we have a discrete series of conformally invariant theories where the potential and gauge field are conveniently given in terms of the weight lattice of the algebra. We compute the exact boundary state for these theories, which explicitly shows the group structure. The partition function and correlation functions are easily computed using the boundary state result.     

Continuous‐wave vertically emitting photonic crystal terahertz laser

Compact semiconductor light sources with high performance continuous‐wave (CW) and single mode operation are highly demanded for many applications in the terahertz (THz) frequency range. Distributed feedback (DFB) and photonic crystal (PhC) quantum cascade (QC) lasers are amongst the leading candidates in this field. Absorbing boundary condition is a commonly used method to control the optical performance of a laser in double‐metal confinement. However, this approach increases the total loss in the device and results in a large threshold current density, limiting the CW maximum output power and operating temperature. In this letter, a robust surface emitting continuous‐wave terahertz QC laser is realized in a two‐dimensional PhC structure by a second order Bragg grating extractor that simultaneously provides the boundary condition necessary for mode selection. This results in a 3.12 THz single mode CW operation with a 3 mW output power and a maximum operation temperature (T_max) of 100 K. Also, a highly collimated far‐field pattern is demonstrated, which is an important step towards real world applications.     

基于时空滤波理论的低频声场干涉结构研究

基于时空滤波理论, 采用相干函数、功率响应函数、扩展分布函数、双频函数四个系统函数较波导不变量更全面地描述了低频声场的相干结构. 理论分析、仿真研究和海试数据处理都验证了低频声场存在稳定的干涉结构, 当目标由远及近、又由近及远做匀速直线运动时, 相干函数(LOFAR图)的干涉结构为一簇类双曲线; 功率响应函数的干涉图表征了简正波群时延的差分, 体现了其频散特性; 扩展分布函数能反映LOFAR图上干涉条纹的斜率大小; 而双频函数能体现相慢度的差分. 各个系统函数均能特别突出干涉结构的某方面的特征, 各有特色.     

Computational and experimental investigation of the fields generated by a 1-3 piezocomposite transducer

Large-scale three-dimensional numerical simulations using the finite-difference time domain technique are used to compute the continuous wave fields associated with a composite transducer. The interior of the transducer is made of a periodic array of square rods. This lattice causes elastic wave Bragg diffraction similar to electrons in a periodic lattice. A low frequency mode shape is assumed for the rods. This prescribed motion includes longitudinal and transverse components. It is shown that the transverse motion in the rod gives rise to shear waves causing standing waves (lateral resonances) in the polymer regions. This is also confirmed by experimental results presented here and other independent analytical and experimental work. The full-scale numerical simulation is performed on a large parallel supercomputer and permits modeling of not only the composite transducer but the radiated pressure from near to far field. In addition, cover plates and edge effects are included, unlike analytical treatments. Although only mechanical effects are included, the wave propagation approach captures many essential features.     

Causal impulse response for circular sources in viscous media

The causal impulse response of the velocity potential for the Stokes wave equation is derived for calculations of transient velocity potential fields generated by circular pistons in viscous media. The causal Green's function is numerically verified using the material impulse response function approach. The causal, lossy impulse response for a baffled circular piston is then calculated within the near field and the far field regions using expressions previously derived for the fast near field method. Transient velocity potential fields in viscous media are computed with the causal, lossy impulse response and compared to results obtained with the lossless impulse response. The numerical error in the computed velocity potential field is quantitatively analyzed for a range of viscous relaxation times and piston radii. Results show that the largest errors are generated in locations near the piston face and for large relaxation times, and errors are relatively small otherwise. Unlike previous frequency-domain methods that require numerical inverse Fourier transforms for the evaluation of the lossy impulse response, the present approach calculates the lossy impulse response directly in the time domain. The results indicate that this causal impulse response is ideal for time-domain calculations that simultaneously account for diffraction and quadratic frequency-dependent attenuation in viscous media.     

Computational simulation of wave propagation problems in infinite domains

This paper deals with the computational simulation of both scalar wave and vector wave propagation problems in infinite domains. Due to its advantages in simulating complicated geometry and complex material properties, the finite element method is used to simulate the near field of a wave propagation problem involving an infinite domain. To avoid wave reflection and refraction at the common boundary between the near field and the far field of an infinite domain, we have to use some special treatments to this boundary. For a wave radiation problem, a wave absorbing boundary can be applied to the common boundary between the near field and the far field of an infinite domain, while for a wave scattering problem, the dynamic infinite element can be used to propagate the incident wave from the near field to the far field of the infinite domain. For the sake of illustrating how these two different approaches are used to simulate the effect of the far field, a mathematical expression for a wave absorbing boundary of high-order accuracy is derived from a two-dimensional scalar wave radiation problem in an infinite domain, while the detailed mathematical formulation of the dynamic infinite element is derived from a two-dimensional vector wave scattering problem in an infinite domain. Finally, the coupled method of finite elements and dynamic infinite elements is used to investigate the effects of topographical conditions on the free field motion along the surface of a canyon.