分层有耗手征介质中斜入射电磁波的传播矩阵

根据相位匹配条件,推导了平面电磁波斜入射情形下分层有耗手征介质中本征波复数波矢量的实部和虚部,二者方向不同使得在介质中传播的本征波为非均匀平面波.然后通过波矢量实部确定了介质中本征波的折射角.最后根据边界条件和本征波场方程给出了分层有耗手征介质中斜入射电磁波的传播矩阵,该传播矩阵可用于解析分析任意入射角情形下分层手征介质的反射透射和电波传播特性.     

无限长非均匀圆柱对任意入射平面波的散射特性

基于无限长非均匀介质圆柱对斜入射平面波散射场公式推导,研究发现当入射角为90°时,斜入射散射公式可以简化为垂直入射的散射公式,所用的算法也可以统一.并用已有算法模拟了无限长非均匀柱粒子对任意入射平面波散射场的强度分布,并且与已有文献结果进行比较,验证了上述结论的正确性.     

б�����Ų��ڷǾ�����ײ���������е�����ͷ�������

提出碰撞吸收型等离子体的等效折射指数,重新计算了折射角,并与利用常规方法得到的折射角进行了比较。分析了等离子体密度和碰撞频率对折射角的影响。考察了斜入射电磁波在吸收型等离子体中的折射特性,并以一维非均匀碰撞吸收型等离子体为例,对斜入射电磁波的反射特性进行了考察。结果表明,等效折射指数的概念涵盖了入射角度的影响,在考察斜入射时更准确、便捷。     

1维介质光栅近场及其衍射的FDTD分析

 应用FDTD方法计算了单色平面波斜入射时1维介质光栅的近场，进而求出介质光栅的衍射效率。借助于周期边界条件，整个计算区域为周期结构的一个单元。考虑入射波在两种介质界面上会产生反射和透射，故在总场边界上引入入射波、反射波和透射波。给出了光栅单元截面为矩形和梯形的算例。该方法可用于斜入射情形下具有复杂结构和任意单元截面的介质栅近场及其衍射特性分析。     

柱状固体复合结构对斜入射超声波的散射

较详细阐述了三层柱状固体复合结构对斜入射超声波的散射理论．厚度远小于声波波长和内层半径的中间圆环可模拟为内外层媒质之间的界面薄层．将声波斜入射时中间国环层力学量的传递矩阵作线性近似展开，我们首次建立了描述界面薄层的三维弹簧物理模型，由此进一步分析了固体界面处的声学边界条件．利用弹簧模型，我们最后从数值上研究了玻璃纤维／铝基复合结构中界面薄层力学参量的变化对斜入射超声散射截面积的影响．     

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

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

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

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

层状介质时域有限差分方法斜入射平面波引入新方式

应用二维时域有限差分方法分析层状介质中的目标散射时,在总场-散射场边界斜入射平面波源用常规方法难以引入,因为在总场-散射场边界处设置的入射波实际上包含了入射脉冲以及各分层界面的反射和多次反射.为解决这个问题,提出了斜入射平面波的混合引入方式,即对总场-散射场的四个边界面采取不同的处理方式.对于总场-散射场的纵向侧边界,用含有斜入射角度的修正一维时域有限差分方法,只要在自由空间位置加入入射脉冲就会自行产生由各分层界面形成的反射波,包括多次反射.同时,把纵向总场-散射场侧边界向下延伸,使得总场-散射场下边界位于完全匹配层内,这样透射波和散射波均为外向行波而被吸收.对于总场-散射场的上边界,由于完全位于自由空间中,边界上各点的入射波将是总场-散射场纵向边界角点处入射波的带有时间延迟的复制.数值模拟结果表明了本文所提出方法的正确性和有效性. 关键词： 时域有限差分 层状介质 斜入射平面波 修正一维麦克斯韦方程     

脉冲平面波斜入射旋转对称体电磁散射的FDTD算法

提出一种应用时域有限差分算法求解脉冲平面波斜入射时旋转对称体散射问题的方法.对斜入射脉冲平面波在圆柱坐标系下进行分解、应用总场/散射场连接边界、PML吸收边界条件.通过与矩量法分析结果的比较,验证了该方法的有效性.     

氧化膜对一维随机粗糙铝表面双向反射特性影响

本文采用时域有限差分方法模拟了一维随机粗糙铝表面双向反射分布函数,分析了表面氧化膜厚度对双向反射分布函数的影响.研究结果表明,氧化膜对一维随机粗糙铝表面的双向反射分布特性有很大的影响。双向反射分布函数受入射波的偏振态影响,横电(TE)波产生的影响比横磁(TM)波更大。在TE波入射时会产生镜反射高峰值现象,但在TM波斜入射时会产生较强的倒反射现象.随机粗糙铝表面在不同的自相关长度(T)条件下的BRDF的峰值均随着T的增大而略微有所增大.     

2(1/2)维等离子体粒子模拟分布式并行程序设计

设计了两种粒子模拟的并行算法,并对其进行了比较,基于消息传递环境开发2(1/2)维粒子模拟并行程序,测试并分析了并行性能.粒子模拟算法中,给出了一个初始化粒子三维Maxwell速度分布的算法,并对常用的电磁场的Lindman吸收边界作了推广.最后对激光钻孔问题作了粒子模拟计算,验证了该并行程序.     

基于网格片的氧碘化学激光器多块并行数值模拟

为实现氧碘化学激光器(COIL)喷管流场大规模数值模拟,采用VICON程序中的基本方程与数值算法,应用JASMIN框架中的基本数据结构——网格片,以及多块结构网格拼接并行算法,发展了三维多块COIL并行模拟程序。数值实验结果展示了该并行模拟程序的正确性及可扩展性,并在2048个处理器核上模拟450万网格单元算例,加速比超过420。     

短波自由电子激光放大器数值模拟

介绍在理论上研究自由电子激光(FEL)的MICFEL程序。MICEFL是以单粒子理论为基础的二维放大器程序。利用MICFEL计算的Paladin实验结果与FRED程序计算结果以及25米Paladin实验结果相比均基本符合。理论计算的结果显示了FEL增益光导的特性。     

激光聚变中受激布里渊散射的混合模拟研究

在激光等离子体相互作用中, 受激布里渊散射的模拟耗时较长. 为缩短模拟时间, 对电子采用Boltzmann描述, 开发了一维混合Vlasov程序, 利用离子束双流不稳定性理论结果对程序进行了考核.在此基础上对激光聚变中的受激布里渊散射进行了研究, 结果表明程序正确, 模拟结果可信.     

北京谱仪触发判选系统的模拟

介绍了北京谱仪(BES)触发判选系统的模拟程序,并用实验数据对其正确性作了验证.运用该程序,首次用BES的蒙特卡罗(Monte Carlo)数据从系统角度直接给出了触发判选效率,所得结果与实验方法测定的触发效率相一致,说明触发判选系统符合BES物理取数的要求.该模拟软件可以用来完备BES的Monte Carlo程序库,并可为今后触发系统的研究提供参考数据.     

Simulation of BES Trigger

A simulation program of BES trigger is described and a comparison of its results with experimental data is presented. With this program, the efficiencies of BES trigger system for different kinds of events generated by Monte-Carlo simulation are given for the first time and they are consistent with those measured by the experimental method. It proves that this trigger simulation code is suitable for studying the performance of BES trigger system and is a Complement of BES Monte-Carlo simulation code library.     

Ka波段TE01模回旋速调管输出腔

 利用模式展开和场匹配理论,建立了具有突变结构谐振腔的散射矩阵。在此基础上编制了FORTRAN计算程序,并对Ka波段TE₀₁模回旋速调管的具有突变结构的开放式输出腔进行了数值计算。分析了回旋速调管输出腔的腔体长度和半径变化,以及输出耦合孔的半径和厚度变化对输出腔的品质因数和频率的影响。结果表明：由程序模拟输出腔所得到的计算结果与HFSS软件的结果一致;与冷测实验结果相比,品质因数的相对偏差为1.65%,谐振频率的偏差为10 MHz。     

基于低密度奇偶校验码的超强前向纠错码设计

光纤中的传输效应(如色散、偏振模色散(PMD)和非线性效应等)会严重影响传输速率和传输距离的进一步提高。因而有必要研究性能更好的新超强前向纠错(Super-FEC) 码型,使其在光通信系统中获得更高的编码增益和更好的纠错性能。简单介绍了低密度奇偶校验(LDPC)码, 然后提出光通信系统中一种基于LDPC码的超强前向纠错(Super-FEC)码型,构造了冗余度为6.68%的新颖LDPC(3969,3720)码, 并给出超强FEC码型的仿真结果。通过与RS(255,239)及BCH(3860,3824)+BCH(2040,1930)进行对比分析, 还分析与探讨了新LDPC码型的编译码电路的设计实现问题。通过理论分析与仿真结果表明这种超强 FEC 码型具有良好的性能, 可以节省硬件开销, 比较适用于光通信系统中, 可作为超强 FEC 码型的候选码。     

中子时空动力学计算程序CORCA-K的进展

主要介绍了自主化中子时空动力学程序CORCA-K的相关研发进展。CORCA-K程序是中核集团自主研发的NESTOR软件包中的三维瞬态中子学计算软件。其核心功能是采用数值计算方法,求解三维瞬态中子扩散方程。其空间离散采用第二类边界条件节块格林函数方法,与三维稳态扩散计算程序CORCA-3D保持自洽性;时间离散采用对角线隐式龙格库塔格式,可在向后欧拉格式的基础上叠加实现,且具备时步自适应功能。以三维LMW基准问题、三维LRA基准问题和秦山核电厂3号机组第5循环为计算对象,开展了瞬态计算验证。结果表明,CORCA-K计算结果与参考程序符合较好,且在保证计算精度的条件下,具有更高的计算效率。     

SSS程序在电磁内爆中的数值模拟

 对SSS程序的边界条件进行修改,将修改后的SSS程序用于固体套筒内爆实验的数值模拟,并与现有的一维磁流体力学程序的计算结果进行比较,对比结果表明,对SSS程序的改造是成功的,也可以预见将SSS程序改造成一维磁流体力学计算程序是可行的。