控制棒落棒时间分析软件CRAC的研发与验证

紧急停堆的落棒时间对反应堆安全至关重要。为适应华龙一号堆型的新型燃料组件设计,中国核动力研究设计院研制出一款落棒时间分析软件CRAC。采用一维流体力学公式结合经验机械阻力模型的方法,构建出CRAC软件理论框架,通过软件开发标准流程完成设计编码,并利用落棒试验数据开展了CRAC软件的验证。结果表明软件计算精度与保守性能满足华龙一号堆型安全停堆时间准则分析的需求。     

压水堆驱动线落棒历程计算

控制棒落棒性能验证是核电厂安全分析的重要部分,研制驱动线落棒历程计算程序有利于验证和改进控制棒驱动线设计。基于驱动线结构特点,分析运动组件的受力情况并进行分解,选择理论或数值方法逐一求取各分力的瞬态值,从而建立驱动线落棒历程的循环步进计算程序。利用秦山核电二期工程驱动线落棒性能试验数据对理论模型和程序计算结果进行对比验证。结果证明：所建立的驱动线落棒历程计算程序适用于压水堆驱动线系统,能正确地对运动组件落棒受力与运动历程进行模拟。     

Transient and stability analysis of a BWR core with thorium–uranium fuel

The kinetic response of a boiling water reactor (BWR) equilibrium core using thorium as a nuclear material, in an integrated blanket–seed assembly, is presented in this work. Additionally an in-house code was developed to evaluate this core under steady state and transient conditions including a stability analysis. The code has two modules: (a) the time domain module for transient analysis and (b) the frequency domain module for stability analysis. The thermal–hydraulic process is modeled by a set of five equations, considering no homogeneous flow with drift-flux approximation and non-equilibrium thermodynamic. The neutronic process is calculated with a point kinetics model. Typical BWR reactivity effects are considered: void fraction, fuel temperature, moderator temperature and control rod density. Collapsed parameters were included in the code to represent the core using an average fuel channel. For the stability analysis, in the frequency domain, the transfer function is determined by applying Laplace-transforming to the calculated pressure drop perturbations in each of the considered regions where a constant total pressure drop was considered. The transfer function was used to study the system response in the frequency domain when an inlet flow perturbation is applied. The results show that the neutronic behavior of the core with thorium uranium fuel is similar to a UO₂ core, even during transient conditions. The stability and transient analysis show that the thorium–uranium fuel can be operated safely in current BWRs.     

窄流道中柔性单板流固耦合数值模拟

板状燃料组件在先进核反应堆中得到了广泛应用。流体以一定流速轴向掠过平行板组件可能导致板的流致振动（FIV）,而板的振动又会影响流场的重新分布,两者之间构成强烈的流固耦合（FSI）关系。针对板状燃料组件的FIV现象开发了计算程序。程序基于物理组成贴体坐标系（PCBFC）,结合任意拉格朗日欧拉坐标法（ALE）实现网格的移动。本工作详细模拟了在窄通道中移动边界条件下流场的分布;数值求解板在流体压力下的梁式振动方程,从而实现窄流道中柔性单板流固耦合的数值模拟。     

基于双区域模型的钠冷快堆组件子通道分析程序的开发与验证

子通道分析程序是钠冷快堆堆芯热工水力设计和安全分析的重要工具。本文为计算和分析钠冷快堆组件在径向均匀与倾斜功率分布工况下的热工水力特性,利用双区域绕丝交混模型开发了一款适用于钠冷快堆组件分析的子通道程序SPLICA,并与FFM2A 19棒束实验数据与WARD 61棒束实验数据进行了对比验证。由于本文开发的子通道分析程序SPLICA使用了详细的绕丝交混模型,与经过二次开发后的COBRA程序的计算结果相比,对于FFM2A实验SPLICA程序计算得到的结果与实验结果符合得更好。这两个实验数据的验证结果证明了本文开发的子通道分析程序的准确性以及对高流量工况和低流量工况均具有良好的适用性。本程序能为钠冷快堆组件热工水力分析提供有效的设计和研究手段。     

数值反应堆堆芯通道级三维热工水力程序CorTAF开发及初步验证

堆芯是核动力系统的核心部件,其完整性是反应堆安全运行的重要前提。传统核反应堆堆芯热工水力分析方法无法满足未来先进核动力系统的高精度模拟需求。本文依托开源CFD平台OpenFOAM,针对压水堆堆芯棒束结构特点建立了冷却剂流动换热模型、燃料棒导热模型和耦合换热模型,开发了一套基于有限体积法的压水堆全堆芯通道级热工水力特性分析程序CorTAF。选取GE3×3、Weiss和PNL2×6燃料组件流动换热实验开展模型验证,计算结果与实验数据基本符合,表明该程序适用于棒束燃料组件内冷却剂流动换热特性预测。本工作对压水堆堆芯安全分析工具开发具有参考和借鉴意义。     

Thermal-hydraulic and neutron-physical characteristics of a new SCWR fuel assembly

A new fuel assembly design for a thermal supercritical water cooled reactor (SCWR) core is proposed. Compared to the existing fuel assemblies, the present fuel assembly has two-rows of fuel rods between the moderator channels, to achieve a more uniform moderation for all fuel rod cells, and subsequently, a more uniform radial power distribution. In addition, a neutron-kinetics/thermal-hydraulics coupling method is developed, to analyze the neutron-physical and thermal-hydraulic behavior of the fuel assembly designs. This coupling method is based on the sub-channel analysis code COBRA-IV for thermal-hydraulics and the neutron-kinetics code SKETCH-N for neutron-physics. Both the COBRA-IV code and the SKETCH-N code are accordingly modified. An interface is established for the data transfer between these two codes. This coupling method is applied to both the one-row fuel assemblies (previous design) and the two-row fuel assemblies (present design). The performance of the two types of fuel assemblies is compared. The results show clearly that the two-row fuel assembly has more favorable neutron-physical and thermal-hydraulic characteristics than the one-row fuel assembly. The effect of various parameters on the fuel assembly performance is discussed. The coupling method is proven to be well suitable for further applications to SCWR fuel assembly design analysis.     

超临界水冷堆双排燃料组件子通道分析

研究基于Cobra-IV程序,开发了适用于超临界水冷堆燃料组件分析的子通道程序.针对超临界水冷堆慢谱双排组件,进行了稳态计算,获取了相关组件热工水力参数.在此基础上,针对单一通道进行了瞬态计算,分析了燃料棒线功率变化和冷却剂流量变化条件下,超临界水冷堆燃料组件的流动和传热的动态响应,为超临界水冷堆组件的优化设计提供了参考.     

燃料组件5×5格架多跨模型CFD模拟方法研究

本文详细描述了某典型燃料组件5×5格架模型CFD分析的几何模型简化、网格划分、求解及后处理等过程。在5×5结构单跨模型上研究了弹簧刚突对搅混特性及压降的影响,并采用简化弹簧刚突的5×5格架模型实现了包含11层格架的多跨模型计算。单跨模型计算结果表明,弹簧刚突结构强化了横向流动,利于换热,Nu提高了8%,但弹簧刚突格架模型较简化弹簧刚突模型压降损失增加了40%。多跨模型计算得到了多层格架全程流动换热特性,为燃料组件自主研发中定位格架数量及布置的设计优化以及DNB预测计算提供了有效的CFD分析方法。     

International Solvent Extraction Conference ISEC ′90

The supercritical-water-cooled power reactor (SCPR) is expected to reduce power costs compared with those of current LWRs because of its high thermal efficiency and simple reactor system. The high thermal efficiency is obtained by supercritical pressure water cooling. The fuel cladding surface temperature increases locally due to a synergistic effect from the increased coolant temperature, the expanded flow deflection due to coolant density change and the decreased heat transfer coefficient, if the coolant flow distribution is non-uniform in the fuel assembly. Therefore, the SCPR fuel assembly is designed using a subchannel analysis code based on the SILFEED code for BWRs.

The SCPR fuel assembly has many square-shaped water rods. The fuel rods are arranged around these water rods. The fuel rod pitch and diameter are 11.2 mm and 10.2 mm, respectively. Since coolant flow distribution in the fuel assembly strongly depends on the gap width between the fuel rod and the water rod, the proper gap width is examined. Subchannel analysis shows that the coolant flow distribution becomes uniform when the gap width is 1.0 mm. The maximum fuel cladding surface temperature is lower than 600°C and the temperature margin of the fuel cladding is increased in the design.     

数值堆热工流体程序CVR-PACA验证及典型应用

谱元方法是一种高精度的数值计算方法,采用该方法开发了数值堆高精度热工水力并行CFD计算程序CVR-PACA。应用CVR-PACA对单棒光棒通道湍流流场、3×3光棒棒束湍流流场、Matis-H压水堆棒束通道基准题、19棒带绕丝组件通道湍流流场进行了仿真计算。通过与实验测量值对比,研究定量验证了大涡模拟(LES)模型及非稳态雷诺时均(URANS)模型对各类棒束通道流场预测的准确性。算例在建模过程中采用网格分裂技术实现了复杂几何的纯六面体网格划分,用于支撑谱元方法计算。研究较为全面地积累了高精度谱元方法模拟流场流动及换热的建模经验,获取了各类棒束通道内丰富的流动和换热细节,获得的建模经验能更加精准有力地指导相关设计的优化改进。     

三维耦合分析软件的落棒事故分析

作为数值反应堆中必不可少的物理和热工部分,中广核研究院有限公司开发了三维物理热工耦合分析软件,通过动态链接库技术实现了自主研发的核反应堆系统瞬态分析软件和三维核设计软件的耦合,并已与国际基准题结果对比验证。本文为耦合软件的应用,围绕华龙一号的落棒分析问题,开展不同落棒组合的耦合计算分析,并研究停堆棒组落棒和温度调节棒（R）棒组两组落棒对堆芯功率的影响。分析结果表明,非中心对称的棒组落棒事故会导致堆芯径向功率出现不对称,并使得堆芯出口回路温度不同。落棒反应性价值越大,R棒调节后的稳态功率回升相比初始稳态差异越大,DNBR公式计算值的变化趋势与功率呈现相反规律。     

快堆燃料组件热工流体力学计算研究

对于钠冷快堆,在燃料和包壳最高温度等设计限值下,为获得较高的堆芯出口温度,需深入分析燃料组件内的热工流体力学问题,准确预测组件内的冷却剂温度分布。本文在CRT模型和F.C.Engel等人工作的基础上,提出了ICRT压降关系式,用以计算冷却剂在湍流区、过渡流区和层流区的棒束压降;引入CRT模型和WEST对流传热模型,改进了SUPERENERGY子通道分析程序,并将改进程序与原程序计算结果进行了对比,结果表明：最热子通道出口温度略有降低,液膜温压略有增加;并用计算流体力学软件CFX对中国实验快堆单盒燃料组件活性段进行了三维数值模拟,将计算结果用CRT模型、ICRT压降关系式及改进后的SUPERENERGY子通道分析程序进行了验证,相互符合较好。     

Methodology for a numerical simulation of an insertion or a drop of the rod cluster control assembly in a PWR

In a pressurised water reactor, the rod cluster control assembly is a system which controls the neutronic activity of the core. It consists of long rods, connected by a spider fixture and a cylindrical system for the control drive mechanism. At its withdrawn position, the activity of the core is maximum, and at its completely inserted position, the activity of the core vanishes. In case of emergency, an effective way to shutdown the reactor is to let it drop under its own weight. An other way to verify the efficiency of the rod cluster control assembly is the insertion test. It consists in inserting the rod into its guides and in checking if the reaction friction force is not high enough to block the movement of the rod cluster control assembly.We present in this paper a methodology for a numerical simulation of an insertion or a drop of the rod cluster control assembly into its guides (discontinuous and continuous guides, guide thimble). A numerical model is elaborated in which many loads are taken into account: fluid load, gravity and friction force between the rod and the guide. The numerical results are compared to experimental measurements obtained from a full-scale structure. A good agreement between the calculated and the measured data is observed.The numerical model takes into account the possible deflection of the guide. It shows clearly that the friction force cannot be neglected when the guide is bowed. So one can locate a faulty guiding system by examining the reaction force during the insertion test. Then, the numerical model can help the decider to make his choice among different rod cluster/fuel assembly components.     

CSR1000堆芯流量分配的研究

以中国百万千瓦级超临界水冷堆(CSR1000)堆芯为研究对象,建立热工水力计算模型,计算出冷却剂和慢化剂温度分布、堆芯功率分布、燃料组件出口压力及流量分配等参数。计算结果表明,适当增加堆芯内部燃料组件流量比例,可以有利于径向功率展平,内外燃料组件通道出口压降,呈现"N"型变化,增大内部燃料组件的堆芯入口功率,内部组件内的流量分配也将减少,而外部燃料组件通道中的流量将增加,适当调整堆芯入口流量初始分配比例,可以使各通道功率分布展平。     

超临界水冷堆环形燃料组件核热耦合分析

在超临界水冷堆预概念设计中,组件设计是十分重要的,将影响堆芯性能。超临界水冷堆中水密度变化剧烈的特性要求必须进行核热耦合分析。从中子学及热工性能角度,使用三维核热耦合程序对环形燃料组件进行了优化设计。应用中子学计算程序FENNEL-N对环形燃料组件进行三维扩散计算,可得到组件内单棒功率分布,应用热工计算程序SUBSC对组件进行子通道分析。在计算过程中,分析了燃料棒间距及燃料棒与组件壁盒之间的间隙对组件性能的影响。计算结果显示,增大棒间距和棒壁间隙能提高组件k_inf,但会增大组件内功率峰因子;子通道受热不均匀性对组件热工性能影响较大,通过加入定位格架的方式能展平冷却剂出口温度,降低最大包壳温度。对环形燃料组件的安全分析表明,从中子学角度该组件是安全的。     

双面慢化环形燃料反应堆物理实验

环形燃料零功率反应堆是首个双面慢化环形燃料作为核燃料的反应堆。本文采用周期法、落棒法获取环形燃料零功率反应堆的临界参数、控制棒价值、元件价值、含Gd元件的反应性效应等关键参数,对环形燃料零功率反应堆的物理性能进行实验研究,验证环形燃料反应堆堆芯物理设计计算程序。结果表明：根据外推过程确定堆芯临界装载环形燃料元件96根,实心燃料元件172根,此时k_eff为1.000 40,堆芯调节棒价值为-247.5 pcm,安全棒价值为-1 358.4 pcm;元件价值与理论值平均偏差为1.3 pcm,含Gd元件反应性效应与理论值平均相对偏差为8.8%。本文结果为环形燃料的工程化设计程序提供关键数据支撑。     

Thermal-fluid flow analysis and demonstration test of a spent fuel storage system

Thermal-fluid flow analysis and demonstration test were performed for a spent fuel storage system. The commercial computational fluid dynamics (CFD) code, FLUENT was used for the numerical analysis. Effective thermal conductivities of a spent fuel assembly and a fuel basket were derived to optimize a thermal analysis model. Also, a porous model, which can simplify a complex configuration of a fuel assembly, was used in the thermal analysis. Demonstration test were performed to verify the thermal analysis method and procedure using a half scaled-down model and an electrically heated dummy fuel. The numerical analysis results were compared with the experimental data. Thermal analyses of the storage system were carried out for normal and off-normal conditions by using the verified analysis method.     

脉动流下棒束通道内相位差及瞬态流场研究

事故条件及海洋条件下反应堆处于非稳态工况,堆芯燃料组件内热工水力行为具有瞬变及多因素耦合特性,对反应堆的安全提出更高挑战,因此有必要对燃料组件内瞬态特性进行研究。本文通过测量棒状燃料组件内压降和流量之间延迟时间开展棒束通道脉动流条件下相位差研究,对比了相位差在不同振幅、不同流动状态下的变化特性,并分析了定位格架对脉动流相位差的作用特点。另外,基于粒子图像测速（PIV）技术开展了脉动流条件下棒束通道内流场分布特性研究,对比了相同流量条件下稳态工况与瞬态工况下流场分布差异,分析了主流具备不同加速度时棒束通道内流场分布特征。实验结果表明：定位格架可减小脉动流下棒束通道内相位差;棒束通道内流场演化滞后于主流量变化。实验结果有助于揭示燃料组件在非稳态条件下瞬态特性,并为燃料组件的设计和优化奠定基础。     

基于刚柔耦合的改进型控制棒组件变形通道落棒行为分析

燃料组件在堆芯内经历长期辐照后易产生弯曲形变，影响控制棒的安全落棒，因此亟需研究变形通道下控制棒落棒行为影响机制。通过数值模拟手段对导向管发生弯曲变形后的落棒行为规律进行分析研究，利用刚柔耦合方法分别计算直型、C型、S型导向管内的落棒行为，分析整个落棒行程、速度、加速度、沿程碰撞力随时间的变化情况，对比直型和2种不同变形通道对落棒行为的影响。研究结果表明，刚柔耦合方法可以较好地模拟变形通道下的落棒行为，C型落棒未发生卡滞，S型落棒卡滞于第2道弯折处。本研究将有助于为弯曲变形导致落棒卡涩问题的极限弯曲阈值提供判断依据，为工程设计提供参考。