有温度反馈阶跃引入负反应性瞬变的解

考虑有温度反馈反应堆阶跃引入负反应性过程的变化规律,对点堆中子动力学方程的解进行研究,给出不同的分析模型,分析发现,（1）用于瞬发超临界的分析模型及其结果不能用于引入大的负反应性过程的分析.（2）在相同的初始功率时,先驱核瞬跳近似模型的结果要优于目前比较普遍的中子瞬跳近似模型的结果.（3）无论初始功率与引入负阶跃反应性的大小,温度瞬跳近似模型的结果总是优于先驱核瞬跳近似与中子瞬跳近似模型的结果.故而温度瞬跳近似模型是目前最为准确的解析模型. 关键词： 中子动力学 温度反馈 反应性 点堆     

有温度反馈阶跃引入负反应性瞬变的解

考虑有温度反馈反应堆阶跃引入负反应性过程的变化规律,对点堆中子动力学方程的解进行研究,给出不同的分析模型,分析发现,（1）用于瞬发超临界的分析模型及其结果不能用于引入大的负反应性过程的分析.（2）在相同的初始功率时,先驱核瞬跳近似模型的结果要优于目前比较普遍的中子瞬跳近似模型的结果.（3）无论初始功率与引入负阶跃反应性的大小,温度瞬跳近似模型的结果总是优于先驱核瞬跳近似与中子瞬跳近似模型的结果.故而温度瞬跳近似模型是目前最为准确的解析模型.     

考虑热传递时有温度反馈的缓发超临界过程分析

对输入小阶跃反应性（ρ0〈β）、有热传递和温度反馈时反应堆缓发超临界过程进行了研究。提出了一个新的物理模型,由数值计算求出任意初始功率条件下反应堆反应性、功率随时间的变化规律,并进行分析讨论,给出了一些有益的新结果。The delayed supercritical process of nuclear reactor with temperature feedback and heat transfer while inserting small step reactivity（ρ0〈β）is analyzed. A new model is proposed. For an initial power, the variations of output power and reactivity with time are obtained by numerical method. The results are analyzed and discussed. Some useful new conclusions are drawn.     

随机中子动力学及其数值模拟研究进展

随机中子动力学是核动力设计和核反应堆安全中的重要课题,本文从随机中子动力学的基础概念和研究方法出发,介绍随机中子动力学研究的历史发展和研究现状.裂变中子与光子的多重性是反应堆零功率中子噪声主要来源,对中子涨落的方程描述及其求解,演化出零功率中子噪声与功率反应堆噪声的随机理论.随机中子动力学的重要应用包括反应性微观测量、功率反应堆噪声测量和分析、核临界漂移分析和核材料识别与检测等.在半个多世纪的研究中,以脉冲堆点火过程的脉冲爆发等待时间分布为代表的随机性,一直缺乏定量分析方法和工具.直到近几年,模拟随机中子动力学过程的广义半马尔科夫过程模拟方法取得了重要进展,很好地揭示了脉冲堆实验中子点火规律.最后讨论随机中子动力学研究中有待解决的研究课题.     

基于FLUENT的核热耦合程序反应性反馈参数敏感性

将计算流体力学模型与中子动力学模型耦合来进行反应堆瞬态安全分析的方法,由于可以开展复杂几何结构的三维流动传热分析,因此受到很大的关注。基于FLUENT用户自定义功能（UDF）开发了一套可用于池式铅堆瞬态安全分析的核热耦合程序,程序耦合了临界/次临界点堆中子动力学模型和燃料棒模型。由于反应堆处于不同寿期时,随着燃料燃耗、可燃毒物积累等因素导致反应性反馈系数有较大变化,因此使用开发的核热耦合程序对中国科学技术大学提出的小型自然循环铅冷快堆进行不同关键反馈系数下无保护的瞬态超功率事故安全分析。调整点堆模块考虑到的四个反应性反馈系数,可以发现燃料多普勒系数对堆安全的影响最大,同时定量的分析结果表明超功率事故引入时间长短对事故演化有重要影响。     

高温气冷堆堆芯中子时空动力学模拟计算

研究开发三维圆柱几何堆芯多群中子时空动力学改进准静态方法模拟计算程序.对给定的模块式高温气冷堆模型进行模拟计算.初始状态下,计算结果与中子扩散程序CITATION吻合很好.动态情况下,模拟堆芯反应性、堆内各能群中子平均注量率和堆芯相对功率等物理量随时间变化,计算结果与理论分析一致.     

一个新的求解点堆中子动力学方程组的数值方法

使用中子密度一阶泰劳多项式分段近似技术,给出一个新的求解点堆中子动力学方程组的数值方法并采用全隐格式以克服方程组的刚性,同时确保解的必要精度。数值结果表明:在隐式一阶多项式近似下,对合适的反应性输入能够取得足够精确的结果。当反应性给定时,对于求解反应堆动力学问题,能给出一个简法的计算过程。     

求解中子动力学方程的加权蒙特卡罗方法

为了实现基于蒙特卡罗方法的中子动力学计算,在传统的直接蒙特卡罗动力学方法的基础上,提出了一种加权蒙特卡罗动力学方法。该方法通过引入粒子权重的概念,隐式考虑中子俘获反应和裂变反应过程中中子数目的变化,避免了模拟粒子的数目随时间的变化,降低了统计偏差,消除了程序计算过程中粒子的存库操作,提高了计算精度。基于单能点堆模型,开发了中子动力学计算程序NECP-Dandi,进行了大量数值验证与分析,包括无缓发中子、单组缓发中子、六组缓发中子、正阶跃反应性引入、负阶跃反应性引入、正脉冲反应性、负脉冲反应性和正线性反应性引入等情况。数值结果表明,相比于直接蒙特卡罗动力学方法,加权蒙特卡罗动力学方法在计算结果的精度和计算效率上有较为明显的改进,程序结构更为简洁。     

散射中子对超瞬发临界系统动态行为的影响北大核心CSCD

在点反应堆理论的基础上,建立了研究散射中子影响超瞬发临界系统动态行为的模拟计算方法,利用CFBR-Ⅱ堆的脉冲实验验证了该方法的正确性。通过分析计算表明,在超瞬发临界状态,散射中子的影响包括增加反应性、使脉冲波形后沿展宽和增加脉冲裂变产额。建立的模拟计算方法为预测脉冲堆外各种物体产生的散射中子的动力学效应提供了途径,对提高脉冲堆的核安全有应用价值。     

中子辐照诱导钨再结晶的模拟研究

钨作为未来聚变堆偏滤器靶板的候选材料,需要长期服役在高温且受到高能中子辐照的严峻环境,这将导致钨发生中子辐照诱导再结晶,从而提高钨发生沿晶脆断的可能性,威胁偏滤器的运行安全,因此研究中子辐照诱导钨再结晶的物理机制具有重要意义.然而,与最近高通量同位素反应(HFIR)堆高温下中子辐照实验观察到的结果相比,目前考虑辐照增强再结晶驱动力效应的模型低估了中子辐照对再结晶的影响,结果表明仍有其他效应影响再结晶过程.基于此,本文在假设晶界迁移率与自体扩散系数成正比的前提下,引入辐照增强晶界迁移因子(R),建立了新的辐照诱导再结晶动力学模型.模拟结果显示,在综合考虑辐照增强再结晶驱动力和晶界迁移效应后,模型计算出的850℃下达到一半再结晶分数所需要的时间(t_X=0.5)和HFIR堆中子辐照实验结果相符,这表明辐照增强晶界迁移效应是影响辐照诱导再结晶现象的重要因素之一.另外,模型研究了不同辐照温度下钨的t_X=0.5.结果表明辐照与未辐照的t_X=0.5差别随温度升高而逐渐下降.这是因为随着温度的升高,辐照缺陷复合加剧,辐照缺陷对再结晶驱动...     

ErP5O14非晶玻璃的红外量子剪裁

研究了Er_1.0P₅O₁₄铒非晶玻璃的红外量子剪裁现象. 从吸收谱和激发光谱的计算比较中肯定了Er_1.0P₅O₁₄非晶 玻璃的1537.0 nm红外荧光为多光子量子剪裁荧光. 从Er_1.0P₅O₁₄非晶玻璃的可见和红外荧光发射光谱中发现激发²H_11/2, ⁴G_11/2和⁴G_9/2能级所导致的⁴I_13/2→⁴I_15/2量子剪裁红外荧光很强;基于自发辐射速率、无辐射弛豫速率和能量传递速率等参数的计算,对其量子剪裁机理进行了分析.发现起源于基态的强下转换能量传递{²H_11/2→⁴I_9/2,⁴I_15/2→⁴I_13/2},{⁴G_11/2→⁴I_13/2, ⁴I_15/2→²H_11/2},{⁴G_9/2→⁴F_7/2,⁴I_15/2→⁴I_13/2}和{⁴G_9/2→⁴I_13/2, ⁴I_15/2→²H_11/2}是导致Er_1.0P₅O₁₄非晶玻璃具有强的三光子和四光子量子剪裁红外荧光的原因.研究结果对改善太阳能电池效率有一定意义.     

核磁共振碳谱研究:烷烃化学位移和(CSS)与分子路径指数矢量(VPM)

系统研究了核磁共振碳谱和化学位移规律及其定量构谱关系(QSSR).本文研究了一组十元素分子路径指数矢量VPM,并发现它与烷烃化学位移和CCS有良好线性相关性.采用多元线性回归进行准确估计与预测,结果优良.     

The ΔNγ electromagnetic transition form factors GM(q) and GE(q)

A parameter-free, nonperturbative calculation of the ΔNγ electromagnetic transition amplitudes GM*(q2), GE*(q2), and the resonant multipole ratio REM(q2)≡E1+3/2(q2)/M1+3/2(q2) is performed in terms of the well-known nucleon isovector Sachs form factor GMV. Our methods are fully relativistic with conservation of the electromagnetic current guaranteed. We find that GM*(q2) decreases more rapidly than the nucleon dipole form factor when −q21 GeV2/c2 and that REM(q2) remains small even for very high four-momentum transfer implying that the perturbative QCD prediction REM(q2)→1 is purely asymptotic and is valid only for extremely high |q2|.     

Interrelations of discrete Painlevé equations through limiting procedures

We study the discrete Painlevé equations associated to the affine Weyl group which can be obtained by the implementation of a special limits of -associated equations. This study is motivated by the existence of two -associated discrete both having a double ternary dependence in their coefficients and which have not been related before. We show here that two equations correspond to two different limits of a -associated discrete Painlevé equation. Applying the same limiting procedures to other -associated equations we obtained several -related equations most of which have not been previously derived.     

Phase formation and magnetic properties of Nd2Fe14B-type Nd16Co76B8−xCx alloys and their hydrides

An attempt is made to synthesize Nd₂Co₁₄C compound by mechanical alloying Nd₁₆Co₇₆B_8−xC_x (0x8) alloys and subsequent annealing. Phase formation and magnetic properties of Nd₂Fe₁₄B-type Nd₁₆Co₇₆B_8−xC_x alloys and their hydrides are investigated. The Nd₂Co₁₄(B,C) phase with Nd₂Fe₁₄B-type structure is formed for Nd₁₆Co₇₆B_8−xC_x (0x7) alloys, while NdCo₇C_δ phase with TbCu₇-type structure is observed in Nd₁₆Co₇₆C₈ alloy. The lattice parameter c of the Nd₂Co₁₄(B,C) phase decreases with increasing the carbon content. A limit volume of the unit cell to form the Nd₂Fe₁₄B-type structure is estimated to be 0.870 nm³. The spin-reorientation temperature T_SR increases with increasing the carbon content, due to an enhancement of magnetocrystalline anisotropy caused by carbon substitution for boron. After hydrogenation, the lattice expansion is observed for Nd₁₆Co₇₆B_8−xC_x (0x7) alloys. The spin-reorientation temperature of Nd₁₆Co₇₆B_8−xC_xH_y (0x7) is much lower than that of the host alloys. Some structural and magnetic properties of hypothetic Nd₂Co₁₄C and Nd₂Co₁₄CH_y compounds are estimated by extrapolation.