差分收敛法对双原子分子高J值转动谱线的预言

为了更精确地预言转动量子数J ≥ 100时双原子分子R支和Q支的振转跃迁谱线, 本文在考虑了转动能级展开式中高阶小项H_v的前提下重新推导出能更好地预言R支和Q支跃迁谱线的物理解析公式. 另一方面, 通过对差分收敛法计算过程的细致分析, 从物理误差的角度提出了一个在没有实验数据作为参照时仍然能有效收敛的重要物理判据. 应用这些新公式和新判据对TiF和CO分子R支振转跃迁谱线和TiF分子Q支振转跃迁谱线进行了研究. 结果表明: 包含高阶小项H_v的新公式的预言结果精度比原有不含H_v的公式的结果提高了一个数量级; 新判据的使用能更有效地减小预言谱线的可能误差, 提高预言结果精度. 最后通过与最小二乘法计算结果的对比, 进一步说明新公式和新判据在预言R支和Q支振转跃迁谱线数据方面的有效性和准确性.     

激光熔蚀金属平板的三维有限体积法数值分析

为了研究气流条件下强激光对金属靶的熔蚀效应,采用有限体积方法建立了数值模型,并开发了三维Fortran计算程序。综合考虑强激光与材料耦合规律、光束能量空间分布、材料高温热物理性能以及熔蚀界面移动等关键影响因素,模拟了激光辐照下金属靶板升温、熔化和剥蚀的复杂物理过程。最后,将计算结果与试验数据进行了比较,验证了计算模型和程序的有效性。结果表明,计算模型能够反映强激光熔蚀金属平板的基本规律,熔蚀深度和后表面温度计算值与试验吻合较好,并且自编计算程序简单高效。     

双原子正方晶格光子晶体边界模式的光传输特性

利用平面波展开法,发现双原子正方晶格光子晶体中ΓM方向边界面存在着快慢两类边界模式,并且通过计算色散关系和电场分布研究了边界参量对这两类边界模式传输特性的影响.依据两种模式的色散关系,计算了群指数和群速度色散参量,结果表明边界参量的变化对第一类边界模式传输特性的影响较小,该模式的平均群指数始终维持在5.0左右;第二类边界模式与第一类模式明显不同,边界参量的变化能够有效地影响到这种模式的传输特性,该模式的最大平均群指数可达178左右.利用时域有限差分法记录了不同时刻电场强度在边界附近的分布及监测点处的电场幅度变化情况,结果表明,两类模式都能够被限制在边界附近并向前传播,时域有限差分法得到的群速度与平面波展开法的结果完全吻合.     

高功率微波照射下半空间上方天线罩耦合特性

利用两次傅里叶变换得到半空间界面的反射波后,将直接入射波和界面反射波两个时域脉冲源通过等效原理引入到包含目标的时域有限差分(FDTD)的计算域中。计算了高功率微波入射情形下,改变入射面、极化方式时天线罩内部观察点的能流密度的变化,并总结其规律。数值算例表明,对于贴近地面的目标来讲,半空间的反射波会给目标内部带来新的能量冲击。该方法可以用于计算目标内部耦合近场及功率密度。     

平行板传输线特性阻抗仿真计算及解析修正

提出了一种基于时域反射原理的平行板传输线特性阻抗仿真计算方法,利用该方法借助CST软件构建模型对平行板传输线的特性阻抗值进行了仿真计算,同时重点研究了下极板展宽结构对特性阻抗值的影响。以仿真计算结果为基础,对原有特性阻抗公式进行了修正,提出了包含下极板展宽结构在内的不对称结构的平行板传输线阻抗计算公式,该公式可作为计算有界波模拟器等电磁脉冲模拟装置特性阻抗的参考公式。     

基于JASMIN框架的抛物方程有限差分解法并行计算及其应用

抛物方程有限差分解法的网格步长严格受波长限制,在求解城市小区电波传播问题时,计算速度明显变慢,为此,基于JASMIN框架研究了抛物方程有限差分解法的并行方法,通过将同一步进面划分成多个网格片,并分配到不同的处理器进行运算,实现了抛物方程有限差分解法的并行计算。与解析解的对比验证了并行程序的正确性,同时通过实例分析了并行程序的高效性,算例表明,抛物方程有限差分解法的求解效率得到了有效的提高。最后,模拟和分析了某一电信基站天线在包含9栋规则建筑物的城市小区环境中的电磁特性,结果表明,该方法能够得到基站在空间各处的信号覆盖强弱,可以为基站选址提供参考。     

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

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

基于涡量矩理论的绕振荡水翼涡动力学分析

文章采用标准k-ω SST湍流模型和动网格技术, 实现了绕俯仰振荡NACA66水翼非定常流动结构与水动力特性的数值模拟, 并基于有限域涡量矩理论定量表征了局部旋涡结构对水翼动力特性的影响. 研究结果表明: 在水翼升程阶段, 当攻角较小时, 层流向湍流的转捩点由水翼尾缘向前缘移动; 在较大攻角时, 顺时针尾缘涡?TEV在水翼吸力面上生成并向前缘发展, 同时与吸力面上的顺时针前缘涡?LEV融合发展为附着在整个吸力面上的新前缘涡?LEV, 新的?LEV与逆时针尾缘涡+TEV相互作用直至完全脱落, 直接导致了水翼的动力失速, 在回程阶段, 绕振荡水翼的流场结构逐渐由湍流转变为层流. 基于有限域涡量矩理论的定量分析发现, 有限域内附着的?LEV和?TEV提供正升力, 当?LEV发展覆盖整个吸力面时对升力的贡献最大, 占总升力近50%, 而+TEV提供负升力. 同时发现, 有限域内各旋涡内部的不同区域提供的升力有正有负; 而逸出有限域的旋涡内部不同区域提供的升力方向均保持一致, 其中顺时针涡提供正升力, 而逆时针涡提供负升力. 在失速阶段, 域外旋涡整体对升力贡献较小且存在小幅波动, 体现了流动的非定常性.      

基于FDM-DEM耦合的冲击损伤大理岩静态断裂力学特征研究

为探究循环冲击损伤后大理岩的静态断裂力学特征,基于有限差分（finite difference method,FDM）-离散元（discrete element method,DEM）耦合的建模技术构建了三维分离式霍普金森压杆（split Hopkinson pressure bar,SHPB）数值模型,其中杆件系统和岩石试件分别采用FLAC3D和PFC3D程序建模。利用该模型对中心直切槽半圆盘（NSCB）试样进行了恒定子弹速度下的循环冲击,随后对受损试样进行静态三点弯曲断裂实验。通过编写Fish程序,提取试样断裂面数据,对断裂面进行重构并定量计算表面粗糙度。通过与相关室内实验结果的对比分析,验证了本文数值分析的合理性与可靠性。模拟结果表明,随着循环冲击次数的增加,试样内部微裂纹、破碎颗粒均增加。连接力场分布混乱,部分力链发生断裂。力链的变化是试样力学性能劣化的根本原因。在静态三点弯曲断裂实验中,冲击5次后试样的静态断裂韧度较天然试样产生一定程度的降低。试样在静载过程中产生的微裂纹和碎块的数量随循环冲击次数的增加而增加,断裂面粗糙度随循环冲击次数的增加而增加。