密立根油滴实验数据的处理方法

依据欧几里德算法，由计算机处理密立根油滴实验数据，得到基本电荷电量，并且给出了用Fortran编写的程序．     

类钠CrⅩⅣ-FmXC离子3s^2S1／2-3p^2P1／2，3／2的跃迁计算

用多组态HFR方法系统地计算了高剥离粪钠CrⅩⅣ-FmXC离子（Z=24-100）单电子组态的能级，相应谱线的跃迁波长，跃迁概率和振子强度，在已有实验的基础上，用自编的Fortran程序对HF计算的3p^2P1／2,3／2能级结果进行了最小二乘拟合计算，从而预测计算了CrⅩⅣ-FmXC离子3s^2S1／2-3p^2P1／2,3／2跃迁的谱线波长。计算结果与已有的实验值均符舍得很好，尤其是与最近测量的实验值相当一致：而与其他理论预测值也基本吻合。     

离心叶轮及无叶扩器内湍流的非结构化网格数值解法

利用非结构化网格上的有限体积法对三维湍流进行了数值研究,提出了不用求解单元顶点处变量值的二阶混合差分格式,既避免了复杂耗时的单元顶点处变量值的计算,又避免了普通差分格式在非结构化网格上易于引起网格界面方向相关性问题。最后利用非结构化网格计算程序,数值求解了一台离心风机内叶轮及无叶扩压器通道内三维相对定常湍流,计算结果与实验值的比较表明速度的计算值与实验值吻合较好。     

一维磁化套筒惯性聚变模拟程序的设计与校验

在调研跟踪国外磁化套筒惯性聚变(MagLIF)研究领域最新进展的基础上,以一维磁流体力学方程为基础,结合氘氚燃料能量交换方程,建立一维磁化套筒惯性聚变物理模型,并使用Fortran语言编写完成一维数值模拟程序MagLIF-1D;通过与圣地亚实验室LASNEX,HYDRA等程序计算结果对比展示,完成程序校验工作,讨论后认为程序计算结果存在差异的主要原因可能来自不同程序对于材料状态方程库的选择;通过计算,MagLIF-1D程序可以直接获得内爆速度、燃料压力、燃料密度、聚变产额等关键物理量,这为后续更好地开展磁化套筒惯性聚变实验设计提供了有力工具。     

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

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

工质热物理性质的计算方法及程序设计

文中就工质物性计算程序中常用的公式法和数据库法进行了讨论 ,并着重介绍了公式法的应用原理和编程方法。以低温工质为例 ,介绍了程序计算模块和界面模块各自的特点和相互关系 ,并分别采用 Fortran和 VC+ +两种语言混合编程的方法 ,实现了程序模块间的连接。此方法为DOS和 Windows不同平台间程序的移植提供了借鉴思路     

Li-Al靶件产量计算方法

ORIGEN2程序经过改造后，已经适合了反应堆的运行工况，并在现有乏燃料的计算中取得了比较满意的结果。因此选用ORIGEN2程序作为Li-Al件产量的计算软件。ORIGEN2程序的结果数据文件，包含了大量原始和结果数据。由于数据繁杂，不利于分析，而且多数数据并非所需，因此需要编制结果数据提取程序，提取需要的数据。     

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

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

精确简洁解析的原子组态相互作用波函数的计算

使用广义Laguerre类型轨道,在组态相互作用框架内,发展并完善了一套精确、简洁、解析的原子(离子)波函数从头计算的Fortran程序,计算原子中电子的结构.通过对指定能态能量求其最小值实现对波函数节点位置和空间范围的优化,从而求出最优波函数.作为例子,计算了⁵G态下的Ti的能量最小值和优化波函数.     

推求j—j耦合谱项的一种新方法—因式化法

推求ｊ－ｊ耦合谱项对于研究重核原子的结构和光谱十分重要，但用人工方法推算其谱项困难很大，其难点和关键在于推算组合组态（ｊ）＾ｎ）ｊ的Ｊ值，我们采用一种新的方法“因式化法”解决了这一难题，该法利用结构化的ＴＲＵＥＢＡＳＩＣ语言编制了计算程序，并在ＰＣ－３８６ＤＸ微机上通过。利用本程序可快速，准确，清晰地推算各种等效电子组态ｊ－ｊ耦合谱项的Ｊ值。同时在程序中对计算结果进行了验证。     

强流脉冲束在螺旋管透镜中的非线性传输计算

为计算强流脉冲束在螺旋管透镜6维相空间中的非线性传输,用Visual Fortran 6.5语言设计了一个计算程序,计算由漂浮空间、螺旋管透镜等元件组成的束流光学系统。程序在计算非强流脉冲束流的线性传输时,粒子的轨迹通过矩阵的直接相乘计算得出;程序在计算强流脉冲束流的非线性传输时,需要考虑束流中的空间电荷效应对束流传输的影响,在束流运动过程中,空间电荷场也在不断地变化,而且粒子运动的轨迹与空间电荷势又是相互依赖的,因此需要求得一个自洽的解,先把元件分成若干均等的区间,把电流分成若干等份,后采用束流电流迭代与元件区间迭代的计算方法。程序运行结束时,横向和纵向相图以及束流光学系统的束流包络线可以在微机屏幕上直观地显示出来。     

红外光谱计算机辅助解析程序

本文介绍一种人工智能红外光谱解析程序,它是我们提出的1R1S红外光谱解析系统的关键部份。此程序在NLC-1180E微计算机上运行时使用Fortran语言,在Perkin-Elmer 3600数据站上运行时使用Basic语言。建立程序的过程中,我们着重作如下三方面努力:(1)在大量实     

推求j－j耦合谱项的一种新方法──因式化法

推求ｊ－ｊ耦合谱项对于研究重核原子的结构和光谱十分重要，但用人工方法推算其谱项困难很大。其难点和关键在于推算组合组态［（ｊ）ｎ］Ｊ的Ｊ值。我们采用一种新的方法“因式化法”解决了这一难题。该法利用结构化的ＴＲＵＥＢＡＳＩＣ语言编制了计算程序，并在ＰＣ－３８６ＤＸ微机上通过。利用本程序可快速、准确、清晰地推算各种等效电子组态ｊ－ｊ耦合谱项的Ｊ值。同时在程序中对计算结果进行了验证。     

类钠CrⅩⅣ-FmXC离子3s2S1/2-3p2P1/2,3/2的跃迁计算

用多组态HFR方法系统地计算了高剥离类钠CrⅩⅣ-FmXC离子(Z=24～100)单电子组态的能级,相应谱线的跃迁波长,跃迁概率和振子强度.在已有实验的基础上,用自编的Fortran程序对HFR计算的3p2P1/2,3/2能级结果进行了最小二乘拟合计算.从而预测计算了CrⅩⅣ-FmXC离子3s2S1/2-3p2P1/2,3/2跃迁的谱线波长.计算结果与已有的实验值均符合得很好,尤其是与最近测量的实验值相当一致;而与其他理论预测值也基本吻合.     

压气机变形叶片的网格生成及区域分解方法

针对压气机优化设计和多学科分析对网格自动可靠生成方法的需求,建立了一套集叶片几何参数化建模、多块结构化网格生成、区域自动分解、叶片自由变形于一体的多功能网格生成系统。基于椭圆型微分方程法和O4H型拓扑结构,建立了压气机单叶道流域的三维分块结构化网格生成方法,改进了适合于变形叶片及网格的自由变形方法,自主开发了相应的计算程序,并用NASA Stage35对方法及程序进行了验证。结果表明网格生成质量较好,特别是近壁网格正交性和光滑性俱佳;改进的自由变形方法能够保持网格周期性和轮毂轮盖几何形状;流场计算结果与实验数据吻合良好。研究工作对发展多级高速轴流压气机优化设计和多学科分析系统具有重要参考价值。     

夸克胶子物质内的qq-↑q→qq-↑q反应

夸克胶子物质内的从夸克-夸克-反夸克到夸克-夸克-反夸克的散射过程会影响夸克物质和反夸克物质的热平衡化，其散射振幅是研究热平衡化的一个基础．依据微扰QCD理论，编制Fortran程序，推导了αs^4阶的各个qq-↑q→qq-↑q散射费曼图的振幅平方和不同费曼图间的干涉项．     

MCNP-4C多粒子输运蒙特卡罗程序的MPI并行化

MCNP程序是由美国Los Alamos国家实验室研制的一个大型、多功能的粒子输运蒙特卡罗程序，可计算任意三维复杂几何系统内的中子、光子、电子或中子-光子-电子耦合输运问题，还可计算临界系统的多种本征值问题。MCNP程序的用户遍及全球，国内用户保守估计也在百家以上，过去主要应用在核科学领域，如今已推广到包括医学在内的许多领域。由于蒙特卡罗计算具有数据独立、循环粒度大、负载均衡的特性。因此，很适合作并行计算。虽然从MCNP-4A程序开始，MCNP程序具有了PVM并行计算功能，但并行系统的开发一直存在这样那样的问题，以致无法正常运行。由于MCNP程序有巨大的计算需求和计算量，我们每年使用MCNP程序完成的计算量就超过万CPU小时。许多问题采用串行计算，时间周期太长。因此，迫切需要研究缩短计算周期的并行程序。     

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

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

非结构网格的生成及新型数据类型的应用

本文基于 Ｄｅｌａｕｎａｙ三角化方法及局部重组的阵面推进法，给出了任意形状平面区域的非结构化网格剖分方法。文中对加点方法及数据结构和存取模式进行了改进，同时应用链表结构、动态数组和派生类型的数据变量，对算法进行优化，节省了计算内存，提高了网格的生成效率。本文研究为特殊计算区域的三维非结构化网格的生成提供了基本方法。     

S=1/2体系的EPR波谱模拟

用一个以Fortran语言编制的S=1/2体系的EPR模拟谱程序.并引入了弛豫项,图谱线型得到了改善,使模拟谱与实验谱更为接近;对超超精细分裂的考虑,使该程序可适合于复杂图谱的模拟;求积分方法的近似处理,缩短了模拟谱的耗时,使复杂图谱的模拟在微机上成为现实.