【精选】点堆中子动力学方程的蒙特卡罗方法.pdfVIP

第 36卷 第 6期 核 动 力 工 程 VbI_36．NO．6 2015年 l2 月 NuclearPowerEngineering D ec． 20l5 文章编号 ：0258-0926(2015)06-0004-06；doi：10．13832／j．jnpe．2015．06．0004 点堆中子动力学方程的蒙特卡罗方法 杨俊云，肖 刚冰，应阳君 北京应用物理与计算数学研究所，北京，100094 摘要：基于广义半马尔科夫过程 (GSMP)模拟方法实现点堆中子动力学方程的蒙特卡罗求解。该方法 模拟裂变系统内中子数和缓发中子先驱核数 目的瞬态演化过程 ，并计算出任意时刻裂变功率和缓发中子源强 等物理量。利用本文提出的方法研究快中子增殖堆 (简称 “快堆”)和热堆参数下的点堆动力学方程，对反 应性的阶跃输入、斜坡输入和振荡输入的点堆中子场瞬态过程进行模拟，并与传统数值算法的计算结果进行 比较。该方法不存在数值计算的刚性问题，能方便地对复杂反应性输入过程进行计算，并能充分考虑瞬态过 程中反应性变化对中子代时间的影响。 关键词：点堆中子动力学；GSMP模拟方法；反应堆启动 中图分类号：TL41l．4；TL327 文献标志码：A AMonteCarloMethodforSolvingPointReactorKineticsEquation YangJunyun，XiaoGang，YingYnagjun ]nstituteofAppliedPhysicsandComputationalMathematics，Beijing，100094，China Abstract：A MonteCralomethodbasedonthesimulationofthegeneralizedsemi—Mrakov processmethodiSdescribedandinvestigatedinthisPaDeLThetransientsofneutronpopulationand precursorpopulationinthefissionableassemblyraesimulatedbvthismethod．Thefissionpower andthedelayedneutronsourceintensityatanytimeraeobtained．ThepointkineticsequationsOf fastreactornadthermalreactorraestudiedbvthismethod．Thetrna sientsofreactorswith step reactivity,ramp reactivity nad oscillatory reactivity rae simulated，which rae systematically compraedwith thetraditionalnumericalresults．ThereiSno stiffnessproblem forthismethod． whichCna convenientlysimulatethecomplexreactivityinputprocess，nadfullyconsidertheeffect ofthereactivitychangeonthekineticsparameters． KeyWOrds：Pointreactorkineticsequations，SimulationofGSMPmethod，Reactorstart—up 0 引 言 点堆中子动力学方程 (PKEs)的准确求解还 中子场瞬态演化过程贯穿于反应堆启动、运 存在困难。在多组缓发中子情况下，很难给出 行、停闭以及事故或应急等工况中，自反应堆研 PKEs的解析解，数值方法成为PKEs求解的重要