光学观测中太阳夹角的分析与计算

对影响光学观测效果的太阳夹角计算方法进行了研究。首先给出了太阳夹角的定义及计算公式。然后重点研究了太阳视位置的计算方法,比较了几种赤纬角计算方法,给出了一种改进的时角计算方法,实现了高精度的太阳视位置计算,并对其计算精度进行了分析。最后提出了在实际光学观测应用中的太阳夹角计算算法。实验表明,太阳视位置计算精度比传统数值近似方法高一个数量级,能够满足光学观测中太阳光照分析的需要。     

快速富里叶变换算法的计算优化

对快速富里叶变换算法的计算进行了优化,介绍了计算优化的方法和计算效率,在各种计算机上进行计算对比,计算速度平均提高了7到10倍。     

船用增压锅炉的热平衡计算方法研究

研究了船用增压锅炉的热平衡计算方法。文章简要叙述了增压锅炉烟气侧的工作过程,重点分析了涡轮增压机组在不同锅炉负荷下功率、热力计算及热平衡计算方法。在实例计算部分,验证了增压机组的热损失数学表达式,并分析其计算结果,最终证明计算方法的正确性。     

蒙卡中子输运程序JMCT和子通道热工水力程序COBRA-EN耦合计算

介绍了中子输运蒙特卡罗方法与热工水力耦合计算的流程。开发了一套蒙卡中子输运程序JMCT和子通道分析程序COBRA-EN耦合接口。通过33棒束模型的计算展示了考虑耦合计算和不考虑耦合计算的差异,论证了耦合计算在反应堆分析中的重要性。通过对反应堆组件的模拟计算,测试了耦合计算的正确性。最后分析了蒙卡计算的统计涨落和迭代计算过程中收敛标准的关系,讨论了蒙卡中子输运和热工水力耦合过程中收敛标准设置的方案和可行性。     

计算物理前沿及其与计算技术的交叉

计算物理学是随着计算技术发展而形成的一门新兴的交叉学科。它是用现代计算技术武装起来的“实验的”理论物理学和“虚拟的”实验物理学。文章回顾了计算物理的发展历史，介绍了计算物理研究状况和前沿问题，并展望了新世纪计算物理面临的挑战。     

LNG球形贮罐绝热计算

分析了贮罐绝热体的传热规律,并建立了模型。通过近似理论计算得到了通过绝热体的传热量;利用计算软件对绝热体内传热过程进行了有限元计算,获得了通过贮罐绝热体的传热量以及绝热体内温度场的分布。以传热量为基准,将近似理论计算与有限元计算结果进行比较,表明有限元计算结果更为准确。     

基于脉管制冷机的二维可压缩流动数值计算程序的开发与应用

本文在SIMPLEC计算程序的基础上,开发了适用于各种流动Mach的计算程序,既能用于计算不可压缩流体,又可以计算可压缩流动。利用此程序计算了一些典型算例,和有关文献结果符合较好。最后利用可压缩程序计算了小孔型脉管制冷机当小孔处于最佳开度时小孔处的气流分布。     

非共轴光学系统传递函数计算

用坐标变换计算非共轴光学系统是十分有效的。本文讨论了非共轴光学系统传递函数(OTF)的计算要点,并给出了平行光束传递函数的计算方法。计算了一个无焦系统的波动光学传递函数。     

超高速发射实验模型的数值计算

 采用三阶精度PPM（Parabolic Piecewise Method）方法和VOF（Volume of Fluid）相结合,运用Lagrange-Remapping算法,编制了多介质流体高精度欧拉计算程序MFPPM,可以对多介质复杂流场进行数值计算。对程序计算精度进行了测试,并应用于超高速发射实验模型数值计算,对Sandia实验室不带会聚作用的实验模型进行了一维计算以及带会聚作用的实验模型进行二维计算并与其实验结果和CTH计算结果进行对比,其中一维计算最大相对误差1%,二维计算相对误差4.7%,在此基础上对超高速发射设计模型进行了初步计算。     

涡轮叶栅气热耦合数值模拟

本文建立了透平机械叶栅流场与叶片温度场的气热耦合计算模型,编制了三维非正交曲线坐标系下的温度场求解程序,计算了实心及空心叶片的温度场分布。介绍了区域分解方法在耦合计算及并行计算中的应用,以矩形区域温度场求解程序为例研究了程序的并行化,最后进仃了叶栅流场计算与叶片导热计算程序的分区耦合计算。     

反应堆二次中子源源强计算及验证

研究建立了反应堆二次中子源源强的机理性计算方法,并计算了某核电厂反应堆第2循环装料时的二次中子源源强。为验证二次中子源源强的计算结果,计算了反应堆装料第2步完成后堆外源量程探测器的计数率,并与实测结果进行对比。源量程探测器计数率计算结果与实测结果符合得非常好,验证了二次中子源源强计算结果的正确性以及所建立的计算方法的合理性。     

任意倾斜角度下的永磁导轨感应磁场解析计算

平行移动式自稳定高温超导块材磁悬浮数值模拟计算包括永磁导轨感应磁场的计算。无限长的永磁轨道结构具有多样性。通过对称分布永磁轨道磁场的解析计算导出了任意倾斜角度下永磁轨道磁感应强度解析计算方法。采用此方法,成功计算出了多块交错排列的Halbach结构永磁轨道的感应磁场计算。采用该方法和商用Ansoft软件的有限元法分别对Halbach永磁轨道磁场进行了计算,对比计算结果表明,此方法是可行的,并且计算效率也远优于有限元法。     

声学时域边界元法的时间卷积计算方法研究

核函数中保留Dirac函数的原型,形成关于时间的卷积积分,是声学时域边界元法中一种稳定、有效的时间数值积分计算方法 (CQ-BEM)。然而,传统CQ-BEM中卷积积分系数的获取有计算量大、耗时长,且对不同单元需要重新计算的问题,极大地降低了CQ-BEM法计算时域声场的效率。针对传统CQ-BEM积分系数计算效率低的问题,本文利用多项式展开定理给出了待求函数泰勒系数的解析表达与数值计算方法,建立了不同单元间待求系数的转换理论,可以在一次循环迭代内完成不同单元的积分系数的计算,大幅降低了计算量,提高了CQ-BEM方法的声场计算效率。脉动球源数值算例结果表明,在相同要求下,本文方法计算时间较传统方法减少50%以上,相对误差小5个数量级以上,且计算时间随单元数的增长率仅为传统方法的2.34%。因此,本文提出的系数计算方法能够有效提高CQ-BEM方法的时域声场计算效率,拓展了CQ-BEM在大型机电设备时域声场模拟的计算规模。     

多孔和复合媒质直流电导率的计算方法(Ⅱ)——应用于各向同性和各向异性复合材料

本文系统地讨论了复合媒质直流电导率的第一性原理计算方法,导出了计算定域各向异性的复合材料的有效电导率的积分方程,提出了计算有效电导率张量的简便方法。我们计算了定域各向异性材料的有效电导率张量。还发展了计算无周期性的更一般的复合材料的理论方法。另外,本文提出积分方程近似计算的一种新方法(双基展开法)。该方法可用于对孔隙度大的复杂结构作精确的理论计算,并具有节省计算时间的优点。 关键词：     

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

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

铀高压状态方程的第一原理研究北大核心CSCD

利用基于密度泛函理论的赝势平面波方法,研究了铀的状态方程。基态计算结果表明:在100GPa的压强范围内,α-U结构是最稳定的相。利用量子分子动力学方法(QMD),计算了有限温度下的状态方程,在此基础上计算了冲击Hugoniot点,并与实验数据及其他理论计算结果进行了比较,结果表明:在100GPa以下的压力区域内,QMD计算结果与实验结果符合得很好,而用快速计算状态方程(QEOS)方法计算的结果偏软。     

计算机控制光学抛光驻留时间求解中两类优化算法的分析

建立了基于矩阵计算的驻留时间计算模型,根据实际加工要求建立了最小二乘和最佳一致逼近最优化求解数学模型,总结了两类优化问题的求解方法。根据自研数学解法器,利用数值计算分析了这两类算法的计算特点。仿真结果显示,两种自研算法具有较高的计算精度,最小二乘逼近算法计算效率有待提高,对外界扰动和计算模型等误差不敏感,最佳一致逼近算法计算效率较高,但对误差比较敏感。实际加工时,如果面形精度已经比较高时,建议多采用最小二乘逼近算法。     

铀高压状态方程的第一原理研究

利用基于密度泛函理论的赝势平面波方法,研究了铀的状态方程。基态计算结果表明:在100GPa的压强范围内,α-U结构是最稳定的相。利用量子分子动力学方法(QMD),计算了有限温度下的状态方程,在此基础上计算了冲击Hugoniot点,并与实验数据及其他理论计算结果进行了比较,结果表明:在100GPa以下的压力区域内,QMD计算结果与实验结果符合得很好,而用快速计算状态方程(QEOS)方法计算的结果偏软。     

轴流式压气机级的特性计算

根据简化的二维流动方程,编写了正问题的计算机程序.对亚声速轴流式压气机级进行了计算,把计算结果与试验结果进行了比较.对这种计算方法作了评价并提出了特性计算中需要解决的问题。     

JMCT在CPR1000工程设计中的验证与确认

CPR1000系列反应堆是目前国内广泛应用的第二代压水堆型号之一，蒙特卡罗程序在CPR1000系列反应堆的验证与确认是该程序实现反应堆工程设计应用的关键环节。基于某CPR1000机组实际参数，使用由国内单位研发的蒙特卡罗程序JMCT在该机组开展了粒子输运建模计算，分别进行了临界计算和固定源计算，并进行了验证与确认。对于临界计算，采用JMCT建立了全堆芯pin-by-pin模型，计算了堆芯有效增殖因子和功率分布。对于固定源计算，建立适用于屏蔽分析的反应堆模型和辐照监督管精细结构模型，计算了两个核电机组多个循环的辐照监督管探测器位置累积快中子注量。通过将JMCT的计算结果与参考程序的计算结果、反应堆实际测量值进行了对比，验证了JMCT程序在CPR1000反应堆工程设计中的实际使用效果，证明了JMCT程序具备工程级的计算精度。