8 mm回旋行波管和返波管的三折螺旋波纹波导色散分析

 三折螺旋波纹波导使TE11模和TE21模相互耦合，其色散曲线能够在宽频带内与电子回旋共振，因此需要对色散方程进行研究。提出了适合工程实用的行波模式和返波模式的色散方程，并对方程中的耦合系数进行了简化，误差在1%左右。利用CST软件的VBA语言对螺旋波纹波导进行建模和计算，根据模拟得到的传输特性曲线的特点，提出一种模拟方法，将模拟与理论计算得到的色散曲线作对比，误差在5%左右。     

CFD在高速飞行器热气弹问题中的应用

采用气动力/热/结构耦合的方法对高速细长体飞行器结构热静气动弹性问题进行了研究.为保证耦合计算精度，达到准确预测热气动弹性特性的能力，气动力和气动热计算采用CFD数值模拟方法，热应力和热变形计算采用有限元方法并通过热考核试验验证.以该简单细长体飞行器模型为研究对象，对其热静气动弹性特性进行了计算与分析，计算结果表明:CFD/CSD耦合可准确模拟热气弹问题，且气动加热造成结构温升不均衡是结构变形的主导因素，力热耦合静气弹变形与单纯受力分析变形形式不同，对飞行器气动特性影响规律不同.准确预测飞行器热气动弹性特性对飞行器结构设计十分必要.      

高频腔本征值问题的三维数值模拟

 运用有限积分法软件MAFIA和有限元软件ANSYS对两类高频腔的本征值问题进行数值模拟，将模拟结果与部分实验数据进行验证，得到了比较一致的结果，其平均计算精度在5%以内，证实了两种电磁软件都可用于高频腔本征值问题的数值模拟，在计算精度上为满足工程设计需要做了一些研究。另外将两种模拟软件在计算时间上作了对比。     

耦合腔慢波特性及其注-波互作用的三维软件模拟

 利用MAFIA软件分别对休斯型耦合腔慢波结构进行了由冷腔到热腔的计算机模拟，并将冷腔模拟结果分别与理论计算结果和实验测试结果进行比较，发现模拟结果更加接近实验值，色散特性的平均误差在0.02%左右，耦合阻抗的平均误差在5%左右，基本满足设计要求；由冷腔模拟结果导入MAFIA三维粒子模拟模块中进行注 波互作用模拟，对其工作机理进行了定性分析。     

水陆两栖飞机静水面高速滑行性能数值计算与试验分析

针对水陆两栖飞机静水面高速滑行过程的运动响应大、流场强非线性等问题，提出了一种基于传统动网格技术的"状态预估——精确计算"的数值模拟方法:通过求解Reynolds平均N-S方程结合运动方程来模拟飞机静水面滑行时的流场特征和运动特性，数值模拟方法为隐式有限体积法，湍流模型采用k-ω（SST Menter）结合壁函数进行处理，自由液面捕捉采用VOF方法；数值计算时，首先采用粗网格对简化后的飞机在不同航速下的姿态和升沉进行快速预估，再将飞机置于预估状态下进行精确网格划分，最后进行精确数值计算分析.为了验证数值模拟结果的正确性，在物理水池中进行了静水拖曳试验，将数值计算结果与试验结果进行对比分析可得:数值计算与水池试验的流场特征吻合，且阻力、姿态和升沉的计算精度达到90%，验证了数值模拟方法的可行性.      

基于时域多分辨算法的非球形气溶胶散射特性仿真模拟

非球形气溶胶的散射特性是影响辐射传输模拟准确性的重要因素.为实现非球形、非均质气溶胶散射特性的模拟,基于MRTD(multi-resolution time-domain)方法建立了一个新的气溶胶散射模型.采用MRTD技术实现了近场电磁场的计算;考虑气溶胶的特殊性,推导了基于体积积分方法的近远场外推方法,实现了粒子散射振幅矩阵和穆勒矩阵的仿真;构建了粒子吸收和消光截面的计算模型,实现了粒子积分散射特性的高精度模拟.将MRTD散射模型的结果与Mie理论、T矩阵法进行了对比,验证了模型的准确性;讨论了空间网格粗细对模拟精度的影响,并定量分析了模型的运行效率.结果表明,MRTD散射模型的相函数模拟误差小于8%,其中前向散射方向小于4%;当粒径与入射光波长相当时,消光和散射效率因子的相对误差小于0.1%;空间网格粗细对模拟精度影响显著,当粒子尺度参数小于20时,在相同模拟精度要求下,所需网格尺寸随尺度参数呈先增大后减小的特征.     

燃烧过程的多尺度模拟与动态自适应机理方法

介绍了一种多时间尺度模拟方法和动态自适应机理方法。在模拟过程中,根据当时的热化学状态动态获得即时简化机理,通过计算各个组分的特性时间,依据其特征时间将基元组分和反应划分为不同的计算组,各个组采用不同的时间步长进行计算。基于该方法,采用详细化学机理对癸烷及庚烷自燃着火过程进行模拟,并将模拟结果与全过程采用详细机理的隐式积分算法进行比较。结果证明该方法具有很高的精度,同时,该方法能大大提高计算的速度。该方法的应用将有助于实现燃烧问题的高效精确模拟。     

表面散射光对激光三角法位移测量的影响

三角测量法是一种位移测量方法, 其测量精度不但受到传感器本身因素如光源、探测器、机械结构特性等因素的影响, 还受到被测表面特性如色泽、材料、粗糙度、倾斜度以及外界环境的影响。国内外的一些学者提出了许多提高传感器测量精度的方法, 但这些研究都是针对某一特性表面进行的, 没有涉及到对不同表面的测量时存在的问题。利用随机行走理论对物体表面的散射场进行了分析, 给出了不同粗糙度表面下弱散射光强与散射角及入射角的关系。由于设计制造的传感器量程范围内接收散射光的角度在15°至25°内变化, 因此取20°作为散射角角度, 理论计算模拟的三种不同粗糙面散射光强变化最大能达到300%左右。该结果对于优化激光三角传感器的设计和提高测量精度有一定的意义。     

用5种不同的电子光学软件对微波管电子枪进行模拟计算

 对TWTCAD微机版、工作站版本、EGUN，CAMEO以及Orprogr电子光学软件背景和特点作了一个简单介绍，并利用其模拟计算了11支不同结构的电子枪。将计算得到的导流系数、射程和注腰半径与实验数据列表对比，结果表明：不同电子光学软件计算不同类型电子枪其计算精度也不一样。在导流系数方面，EGUN和Orprogr的计算精度较高，绝对平均误差分别为7.18%和7.4%；在射程方面，TWTCAD的计算精度最高，绝对平均误差为5.69%；在注腰半径方面，Orprogr与实验值偏差较大，绝对平均误差为47.43%，而工作站版本的计算精度最高，绝对误差仅为2.29%。     

基于墨量限制样本的打印机光谱特性化模型修正方法

打印机色彩特性化是打印色彩管理技术的关键。打印机特性化模型修正，是指针对因打印介质及墨盒更换等原因而造成的特性化模型精度下降问题，借助特定修正样本实现原始模型精度校正的过程。为进一步提高此类修正方法的建模精度与效率，提出了一种基于墨量限制样本的打印机光谱特性化修正方法，该方法充分利用新介质墨量限制过程所制备的墨量限制样本，在无需额外修正样本的条件下，实现了原始特性化模型的有效修正。研究以三类不同类型打印介质为例，以现有典型特性化模型修正方法为对照，通过实验对本文方法的有效性进行了验证。结果显示，该方法相较现有方法展现出更为优异的修正效果，在无需额外修正样本的条件下，其光谱修正精度可提升15%～20%，色度修正精度可提升10%～20%，实现了修正精度与效率的同步优化。     

天文导航中星体高精度细分定位方法研究

介绍了一种从星体图像中高精度提取恒星位置的方法.这种方法把星光成像看成是高斯点扩散函数模型,利用线性内插和最小二乘法拟合得到高斯曲面参数,从高斯曲面模型中得到亚像素级的恒星位置.由于直接进行高斯曲面拟合计算非常复杂,将二维曲面拟合转化为两个方向上的一维曲线拟合,分别得到不同的曲面系数.仿真结果表明当信噪比小于0.05时,定位精度能达到1/20像素,与质心法比较,高斯曲面拟合法精度高,抗噪声能力强.     

局部间断Petrov-Galerkin方法在大气污染模型中的应用

构造数值模拟两类大气污染模型的局部间断Petrov-Galerkin方法。首先通过变量代换将大气污染模型方程转化为与之等价的一阶微分方程组,再利用间断Petrov-Galerkin方法求解微分方程组。该方法既可以选取不同的检验函数和试探函数空间,又可以保持间断Petrov-Galerkin方法的优势。同局部间断有限元方法相比,该方法的计算公式较为简便,数值算例表明该方法具有三阶精度。与有限体积元方法相比,该方法具有较小的误差。本算法可为大气污染模型的数值模拟提供实用工具。     

核磁共振陀螺中原子气室温度场的研究

为研究核磁共振陀螺中加热机构对原子气室性能的影响,设计了5种典型加热方式。利用有限元分析软件ANSYS建立了原子气室的温度场模型,给出了原子气室表面的稳态温度场分布情况。同时设计了探测精度为0.01℃的测温电路,对原子气室表面不同位置的温度进行监控,获得了不同加热方式下原子气室表面的温度变化情况。将仿真和实验结果进行比较,发现误差在5%之内,验证了仿真模型的正确性。综合仿真和实验结果比较了不同加热方式下原子气室表面温度分布情况,获得了能够使原子气室表面温度分布最均匀的加热方式。     

基于JSNT程序的压水堆屏蔽计算

反应堆屏蔽计算是评估核电站安全性能的基础，是指导电站设计、运行的重要手段之一。JSNT程序是中物院高性能数值模拟软件中心研制的大规模并行离散纵标输运程序，具有较高的计算精度和计算效率。利用JSNT程序对某压水堆进行了建模计算，给出了中子通量密度的分布结果，并与实验测量值进行了对比。对比结果显示:无论是采用S8计算还是S16计算，计算结果都能满足工程要求；相比S8而言，采用S16计算可以显著提高计算精度，能够将某些测点处的相对误差降至1%以内。     

Computational comparison of the conventional multislice method and the real space multislice method for simulating exit wavefunctions

The simulation of exit wavefunctions of a TEM sample is carried out by both the conventional multislice (C-MS) method and the real space multislice (RS–MS) method, with a series of different slice thicknesses. Surprisingly, the results calculated by the RS–MS method remain constant to seven significant figures when the slice thicknesses are changed within a wide range, while those calculated by the C-MS method change significantly. From the mathematical point of view, the C-MS method is only a first order approximation of the RS–MS method. Thus, the calculation precision of the RS–MS method is much higher than that of the C-MS method. Or to achieve the same precision, the slice thickness for the RS–MS method can be taken as two orders of magnitude larger than for the C-MS method in the simulation. As a result, the calculation speed for the RS–MS method may be even faster than that for the C-MS method in some cases. For the two advantages stated above, the RS–MS method may be more suitable for precise calculations.     

考虑铺层角的复合材料圆柱壳自由振动三维弹性准确解

为了获得任意角度铺层的多层复合材料圆柱壳的自由振动准确解,在三维弹性理论的基础上,结合分层理论和状态空间法,建立横向位移和应力的传递矩阵,轴向和环向位移采用双螺旋模式的位移函数,对任意角度铺层复合材料圆柱壳简支边界条件下的自由振动进行了理论推导,得到了自由振动方程的精确形式。与文献理论解和有限元计算结果对比,结果表明,关注频率在2倍的环频率以下时,薄壳的固有频率计算精度能控制在1%以内,厚壳的固有频率计算精度能控制在2%以内。对于厚壳的计算可将壳体沿厚度方向划分为多层来处理,这样能有效提高计算精度。计算分析了铺层角对壳体固有频率的影响,环向模态数较低时,固有频率随着铺层角的增加呈抛物线变化趋势;环向模态数较高时,固有频率随着铺层角的增大单调递增。该理论方法同样适用于均质各向同性壳和正交各向异性圆柱壳。      

双数态输入型M-Z干涉仪中量子相位的最大似然估计和贝叶斯分析

本文采用最大似然估计和贝叶斯分析两种方法，对双数态输入型马赫-曾德尔干涉仪中的量子相位估计进行研究.通过对干涉仪输出端粒子数差和宇称进行理论计算和数值模拟，发现采用贝叶斯分析和粒子数差测量的组合，可以在全相位空间实现最优测量，即待估相位精度达到由量子克拉美罗下界所给定的测量极限，同时，利用贝叶斯分析进行相位估计所需样本数较最大似然估计更少.在通过宇称测量进行相位估计的研究中，发现无法利用贝叶斯分析实现相位估计，但可采用最大似然估计进行研究.借助蒙特卡洛数值模拟分析，我们验证了理论计算结果，即相位估计精度会随待估相位θ₀的变化而改变，将其与粒子数差测量结果对比，进一步确认了粒子数差测量方案的优越性.     

点堆中子动力学方程组曲率权重法的解

针对核反应堆安全工程对某些数值计算结果要求较高的精度和正的误差, 以及舰船核反应堆机动性对计算速度的要求, 需要从数学上寻找一种新的数值计算方法, 以满足实际曲线向上凸或向下凹时计算值总是略高于真实值, 且误差不大于欧拉法和改进的欧拉法所得值. 本文研究曲率权重法求解点堆中子动力学方程组, 该方法是在曲率圆法的基础上引入权重的思想来衡量间隔步长上两个曲率对该步长曲率平均值的贡献. 与欧拉法和改进的欧拉法比较, 曲率权重法的计算结果总是能够高于真实值或有正的误差, 且精度和计算速度得到明显提升. 将该方法用于次临界堆阶跃和线性引入反应性时中子密度的求解, 能够快速得到满足计算要求和高精度的数值结果.     

Precision Calculations of Atomic Polarizabilities: A Relevant Physical Quantity in Modern Atomic Frequency Standard

Electric dipole polarizabilities of atoms are very important in many different physical applications, such as the precision atomic frequency standard. Calculations of these properties are very important and challenging. We propose a calculation strategy to calculate the frequency dependent dipole polarizabilities with high precision variationally by using a set of high quality orbital bases where the electron correlations can be taken into account adequately. The static polarizabilitiez of the ground state of Na are calculated accurately by such a method and can be compared with precision experiment measurement directly. The calculation result is in excellent agreement with the available experimental measurements within about 0.1~, which demonstrates the validity of our strategy. Our calculation strategy has a wide usage, not only in polarizibilies, but also in other fields such as theoretical treatment of electron-atom scattering processes. Using the same orbital bases, we carry out precision calculation of Na- affinities. Our calculated affinity is in excellent agreement with precision laser spectroscopy measurements within 0. 1%.     

水下爆炸冲击波数值仿真研究

由于在水下爆炸冲击波的数值仿真研究中,水的状态方程、人工黏性系数和网格尺寸对数值计算结果影响很大,采用常规TNT炸药的水下爆炸为例,以冲击波的峰值压力和比冲量为衡量指标,研究了这3个主要影响因素对数值仿真结果的影响。首先,通过采用常用的5种水的状态方程进行系列仿真,给出了各种状态方程的适用范围;其次,讨论了人工黏性系数对计算结果的影响,并给出了一次与二次人工黏性系数的建议取值范围;最后,通过对不同炸药当量及不同网格尺寸开展系列运算,从而得到不同炸药当量在满足工程计算精度要求下所对应的建议网格尺寸,并得到了不同炸药当量所对应的建议网格尺寸的表达式。