基于JMCT-JBURN的燃料棒径向功率分布计算北大核心CSCD

研究给出基于蒙特卡罗粒子输运软件JMCT耦合燃耗分析软件JBURN计算燃料棒径向功率分布的方法。介绍了计算模型的建立、输运及燃耗计算的相关设置以及裂变功率、俘获功率的统计方法。UO_(2)燃料棒径向功率分布的计算结果表明,采用JMCT-JBURN软件的结果与工程常用软件符合较好,最大偏差不超过3%,证明了计算方法和计算软件的适用性。该方法可用于工程设计。     

液态水中重离子径向剂量分布的计算

根据建立的重离子径迹结构Monte Carlo计算模型,计算了重离子在液态水介质中的径向剂量分布和最大作用范围,计算结果同实验或其它δ射线理论计算值进行了比较,其目的在于建立细胞损伤强度和损伤范围与重离子参量之间的关系,从理论上探索重离子生物效应的机理.     

中国实验快堆燃料组件精细功率分布计算

在热工中,需要中子学计算给出燃料组件内各元件棒功率的相对分布。利用蒙特卡罗程序对中国实验快堆（CEFR）燃料组件内元件棒功率分布进行了理论计算分析,并保证计算结果相对统计误差小于0.8%。使用另一个基于六角形节块扩散理论的钠冷快堆中子学设计软件NAS计算得到的结果对蒙特卡罗程序计算结果进行了对比计算。结果表明,蒙特卡罗程序与NAS计算得到的元件棒相对功率分布结果的最大相对偏差小于3%。使用蒙特卡罗程序对CEFR燃料组件内精细功率分布的计算是可靠的,可用于设计计算当中。     

中国实验快堆燃料组件精细功率分布计算

在热工中,需要中子学计算给出燃料组件内各元件棒功率的相对分布。利用蒙特卡罗程序对中国实验快堆（CEFR）燃料组件内元件棒功率分布进行了理论计算分析,并保证计算结果相对统计误差小于0.8%。使用另一个基于六角形节块扩散理论的钠冷快堆中子学设计软件NAS计算得到的结果对蒙特卡罗程序计算结果进行了对比计算。结果表明,蒙特卡罗程序与NAS计算得到的元件棒相对功率分布结果的最大相对偏差小于3%。使用蒙特卡罗程序对CEFR燃料组件内精细功率分布的计算是可靠的,可用于设计计算当中。     

Terfenol棒径向开缝数目的计算

本文给出了一种简便的计算方法，用以计算满足极限工作频率要求的切割份数，并附以算例，与实际相符合。     

粒子的径迹结构与径向剂量分布

概述了重离子径向剂量分布的实验和用连续慢化方法计算径向剂量方面的进展。     

MCNP5在固态燃料熔盐堆功率分布计算的应用

采用蒙特卡罗输运程序MCNP5对固态燃料熔盐实验堆(TMSR-SF1)能量沉积比例及功率分布进行了计算分析。针对MCNP5不能处理缓发β及缓发γ的能量沉积问题进行了类比等效处理。对固态燃料熔盐实验堆在寿期初、寿期中、寿期末相应的能量沉积比例及功率分布进行了研究。通过计算发现,固态燃料熔盐实验堆内燃料球相比于压水堆棒状燃料元件(95%~97%左右)而言,能量沉积比例有所偏小,约为93%。同时,由于堆芯功率分布均匀,功率峰因子较小(约1.5),堆芯安全性较好。     

聚氨酯结构与径向分布的分子动力学模拟

采用分子动力学方法计算得到DHI-乙烯醇聚合体系统的结构和径向分布函数.讨论了系统结构和径向分布函数与温度和压力之间的关系.结果表明粘合系统的空间分布一般地随着温度和压力的增加而收窄,对增加聚氨酯系统的粘合性具有积极的意义.     

聚氨酯结构与径向分布的分子动力学模拟

采用分子动力学方法计算得到DHI-乙烯醇聚合体系统的结构和径向分布函数。讨论了系统结构和径向分布函数与温度和压力之间的关系。结果表明粘合系统的空间分布一般地随着温度和压力的增加而收窄,对增加聚氨酯系统的粘合性具有积极的意义。     

基于Dakota的燃料棒设计验证不确定度研究

燃料棒设计验证是评价燃料棒在反应堆内运行时安全性能的过程,其中输入参数的不确定度对评价结果有非常重要的影响。为了系统研究燃料棒设计验证的不确定度,使用Dakota中蒙特卡罗与拉丁超立方的非参数抽样方法,结合燃料棒性能分析软件开展了燃料棒设计验证计算,并与传统的不确定度计算方法进行了比较。结果表明,传统方法未充分考虑输入参数的不确定度,导致内压准则在正常运行条件下容易受到挑战,统计性的抽样方法弥补了这一缺陷,获得了较大的安全裕量,为燃料棒安全性以及经济性的提升提供了理论依据;同时,2种抽样方法所获得的燃料温度计算结果较传统方法更加具有参考意义;对于包壳腐蚀准则以及包壳应变准则,由于不确定度输入参数选取得当,抽样方法与传统方法的计算结果无明显区别。因此,基于非参数抽样的统计法对于评价燃料棒在反应堆内的安全性能更加具有实用性。     

角点检测不确定度评定方法

 提出一种评定角点检测不确定度的新方法。建立了一个像素强度不确定度数学模型,该数学模型在像素强度不确定度与像素强度和图像梯度之间建立了联系。给出了实现角点检测不确定度评定的具体方法：首先利用像素强度不确定度数学模型估计每个参与角点检测的像素的强度不确定度,然后再用蒙特卡罗方法估计出角点检测的不确定度。以经典SUSAN角点检测算子为例,对文中提出的角点检测不确定度评定方法进行了验证。实验结果表明,使用该方法能准确地估计出角点检测结果的不确定度。     

同一燃耗区域核燃料包壳破损建模分析

核反应堆燃料包壳是构成反应堆安全的重要屏障, 一旦发生破损,放射性裂变产物就会释放到一回路的冷却剂中。本文通过一回路134Cs/137Cs比值确定破损燃料棒所在燃耗区域,对该燃耗区域裂变产物的产生、裂变产物进入芯块间隙、一回路核素平衡分别进行建模,然后,利用建立的数学模型对破口处裂变产物向一回路释放过程进行分析,得出破口的大小和破损根数及所在燃耗区域,并以压水堆核电站燃料包壳破损的数据进行验证,证实了模型的合理性。     

高温堆乏燃料贮罐中子屏蔽性能计算

球床模块式高温气冷堆采用包覆颗粒球形燃料元件,在反应堆运行过程中,不断排出的乏燃料球将被装入乏燃料贮罐。乏燃料贮罐应选取适当的材料和厚度,对光子和中子进行有效屏蔽,使罐外的剂量率满足相应的限值要求。为此,使用张弛长度法和蒙特卡罗模拟法研究乏燃料贮罐的中子屏蔽性能。屏蔽材料为铁和含硼聚乙烯,计算了铁和不同B4C含量聚乙烯的屏蔽性能,并给出了乏燃料贮罐装满乏燃料球后,乏燃料球自吸收对贮罐外剂量率的影响。两种方法计算结果吻合很好,可以为实际工程中的屏蔽设计提供参考意见。     

高温堆乏燃料贮罐中子屏蔽性能计算

球床模块式高温气冷堆采用包覆颗粒球形燃料元件,在反应堆运行过程中,不断排出的乏燃料球将被装入乏燃料贮罐。乏燃料贮罐应选取适当的材料和厚度,对光子和中子进行有效屏蔽,使罐外的剂量率满足相应的限值要求。为此,使用张弛长度法和蒙特卡罗模拟法研究乏燃料贮罐的中子屏蔽性能。屏蔽材料为铁和含硼聚乙烯,计算了铁和不同B4C含量聚乙烯的屏蔽性能,并给出了乏燃料贮罐装满乏燃料球后,乏燃料球自吸收对贮罐外剂量率的影响。两种方法计算结果吻合很好,可以为实际工程中的屏蔽设计提供参考意见。     

Shadow hybrid Monte Carlo: an efficient propagator in phase space of macromolecules

Shadow hybrid Monte Carlo (SHMC) is a new method for sampling the phase space of large molecules, particularly biological molecules. It improves sampling of hybrid Monte Carlo (HMC) by allowing larger time steps and system sizes in the molecular dynamics (MD) step. The acceptance rate of HMC decreases exponentially with increasing system size N or time step δt. This is due to discretization errors introduced by the numerical integrator. SHMC achieves an asymptotic O(N^1/4) speedup over HMC by sampling from all of phase space using high order approximations to a shadow or modified Hamiltonian exactly integrated by a symplectic MD integrator. SHMC satisfies microscopic reversibility and is a rigorous sampling method. SHMC requires extra storage, modest computational overhead, and a reweighting step to obtain averages from the canonical ensemble. This is validated by numerical experiments that compute observables for different molecules, ranging from a small n-alkane butane with four united atoms to a larger solvated protein with 14,281 atoms. In these experiments, SHMC achieves an order magnitude speedup in sampling efficiency for medium sized proteins. Sampling efficiency is measured by monitoring the rate at which different conformations of the molecules' dihedral angles are visited, and by computing ergodic measures of some observables.     

基于能量密度分布的辐射源粒子空间抽样方法研究

利用隐式蒙特卡罗方法模拟热辐射光子在物质中的输运过程时,物质辐射源粒子是需要细致处理的物理量.传统的物质辐射源粒子抽样方法是体平均抽样方法,对于大多数问题,这样处理不会带来大的偏差.但是对于一些辐射吸收截面大、单一网格内温差显著的问题,体平均抽样方法的计算结果偏差较大.分析了产生偏差原因,提出一种基于辐射能量密度分布的辐射源粒子空间位置抽样方法,并推导了相应的抽样公式以解决此类问题.数值实验表明,新方法计算结果明显优于原方法且与解析结果基本一致.     

基于蒙特卡罗方法的中子屏蔽材料设计

基于蒙特卡罗粒子输运程序MCNP,设计了一种强度高、密度低、具有优异中子屏蔽性能的新型玻璃纤维/B4C/环氧树脂复合材料,模拟计算了镅-铍（Am-Be）中子源产生中子对该材料的透射率;研究了该材料的中子屏蔽性能与传统屏蔽材料的差异以及不同B4C质量分数对该材料的屏蔽性能影响;根据模拟结果分析了该材料对不同能区中子（慢中子、中能中子、快中子）具有的不同屏蔽性能。研究发现：B4C质量分数为10%的该种新型玻璃纤维/B4C/环氧树脂复合材料的中子屏蔽性能,尤其是慢中子屏蔽性能较传统的含硼聚乙烯和Al-B4C合金材料更为优异;但随着B4C质量分数的增大,屏蔽性能提升不明显。结果验证了蒙特卡罗方法用于中子屏蔽材料优化设计的可行性。     

On dynamics of elastic rod based on exact Cosserat model

The analysis of kinematics and dynamics of an elastic rod with circular cross section is studied on the basis of exact Cosserat model under consideration of the tension and shear deformation of the rod. The dynamical equations of a rod with arbitrary initial shape are established in general form. The dynamics of a straight rod under axial tension and torsion is discussed as an example. In discussion of static stability in the space domain the Greenhill criteria of stability and the Euler load are corrected by the influence of tension and shear strain. In analysis of dynamical stability in the time domain it is shown that the Lyapunov and Euler stability conditions of the rod in space domain are the necessary conditions of Lyapunov's stability in the time domain. The longitudinal, torsional and lateral vibrations of a straight rod based on exact model are discussed, and an exact formula of free frequency of lateral vibration is obtained. The free frequency formulas of various simplified models, such as the Rayleigh beam, the Kirchhoff rod, and the Timoshenko beam, can be seen as special cases of the exact formula under different conditions of simplification.     

基于高斯原理的Cosserat弹性杆动力学模型

在动力学普遍原理中, 高斯最小拘束原理的特点是可通过寻求函数极值的变分方法直接得出运动规律, 而无须建立动力学微分方程. Kirchhoff动力学比拟方法以刚性截面的姿态表述弹性细杆的几何形态, 并发展为以弧坐标s和时间t为自变量的弹性杆分析力学. 由于截面姿态的局部微小改变沿弧坐标的积累不受限制, Kirchhoff模型适合描述弹性杆的超大变形. Cosserat弹性杆模型考虑了Kirchhoff模型忽略的截面剪切变形、中心线伸缩变形和分布力等因素, 是更符合实际弹性杆的动力学模型. 建立了基于高斯原理的Cosserat弹性杆的分析力学模型, 导出拘束函数的普遍形式, 以平面运动为例进行讨论. 关于弹性杆空间不可自相侵占的特殊问题, 给出相应的约束条件对可能运动施加限制, 以避免自相侵占情况发生.