核反应堆动力学时空动力学.ppt

点堆动力学 比事先计算形状函数的绝热近似方法更简化的办法就是采用点堆动力学了。 从头到尾只用一套动态参数（根据初始形状函数算出，或从实验获得） 从时空动力学到点堆模型 点堆方法虽然位居末尾，有诸多缺点，但是优点也不少。 可以模拟堆的许多动态行为．反应性，温度系数，功率系数，空泡系数等等物理量，都能方便得用于点堆方程。不能用于时空动力学方程。 Why？因为:在时空动力学方程中,任何量的变化都是以截面的改变来表示的,没有反应性这个量. 时空动力学方程的直接数值解法 转word文件  格林函数节块法等先进粗网节块法也可以用来解时空动力学方程。 与空间有关的中子学现象 主要有三个方面 通量分布形状函数随时间的变化。 缓发中子对通量倾斜（或曰畸变，英语为 tilt）的滞阻作用。 静态反应性与动态反应性的区别 动态过程中通量的空间分布是变化的 在动态过程中，通量的空间分布，与在该堆内人为地减小裂变率而达到临界而得出的稳态分布（稳态计算就是这样干的）显著不同。 缓发中子的滞阻效应 缓发中子先驱核的空间分布情况取决于早先时刻的裂变率，因而取决于早先时刻的通量。因此缓发中子源具有的空间分布更反映以前的而不是当前的裂变率分布，因此它对通量的畸变具有滞阻作用。 静态反应性与动态反应性 这个提法似乎有些奇怪。反应性总是和动态联系在一起的，静态时还有什么反应性？就是零嘛！ 不奇怪。过去我们进行静态堆芯物理计算时，算的就是通量分布和有效增值系数 k，而反应性就是用 k 来定义的。 两种反应性皆用微扰公式表示 　 罗璋琳先生在“实验反应堆物理”中从另一个角度来区分动态反应性和静态反应性 静态反应性 在反应堆物理实验中，反应性的测量是把中子探测器放在堆中某一确定的位置，测量中子数随时间的变化，从而求出反应性。 凡是这种以反应堆动态方程为出发点，根据堆中子的时间行为来确定的反应性称为动态反应性。这里反应性被当作一个动态参数。 动态反应性 总结 1。 时空动力学方程中用不到反应性、缓发中子有效份额等物理量。WHY？ （点堆方程是集总参数模型，时空动力学方程是分布参数模型。正如我们研究电路时用电阻、电容、电感等集总参数，研究电场则用麦克斯韦方程） 2.时空动力学的准静态处理，需要用到微扰理论、价值理论。推导出的精确点堆方程中的各项动态参数（反应性等）是与实际动态过程的形状函数有关的。  3.时空动力学问题的直接数值解法还有许多，但那都是数学技巧问题，由于时间有限，不一一列举。 要求大家懂得 1. 静态反应性与与动态过程毫无关系，是一个以微分方程（扩散或输运方程）特征值来定义的量。 2. 动态反应性是与动态过程有关的量。因为其定义与时变的形状函数有关。 3。 虽然直接求解时空动力学方程是纯数学问题，而且能得到最精确的结果。但通过时空动力学的准静态处理，使我们对于实际动态过程有更深入了解。 4.你可以不用准静态方法来解时空动力学方程，但是你必须了解与空间有关的各种动力学现象。 用蒙特卡罗方法求解时空动力学方程 我系有一位博士研究生（沈华韵）在做此项研究． 点堆方程的数值解法 由于点堆方程的重要性，而在一般情况下不能得到其解析解，故在此较详细讨论一下其数值解法。 在此之前，先对与常微分方程数值解有关的的理论问题做一番介绍。 转word文件 厄米特分段多项式方法 三阶厄米特分段多项式方法是求解点堆方程的有效方法。精度高，计算速度快，适应性强。 转word文件 编程作业 时空动力学比较复杂，编程比较困难、耗时。 编一个解点堆方程的程序。算几个例题。 对求解点堆方程的算法进行调研，写一个简单报告。 如果能编一个一维时空动力学方程，可以免期末考试。 参考书 Karl O.Ott 的书 STACEY的书 黄祖洽的书 谢仲生的书 核反应堆动力学 反应堆动力学研究堆的动态行为。 在初等反应堆物理中，只是很简略的介绍了点堆动力学。 国外大学核工系大多单独开设反应堆动力学课程，因为其内容很丰富。国内黄祖洽院士写了一本核反应堆动力学的专著 ，有数百页的篇幅。 我们限于学时，只讲3小时。 但是基本内容都有了。 点堆模型局限性 不能描述与空间分布有关的动态现象 当堆内发生非均匀扰动时不宜用点堆模型分析。如果采用点堆模型，则 误差很大； 计算结果不一定是偏于保守的。 对上页的说明 左图： 右图： 何谓中子时空动力学？ 所谓中子时空动力学，不过是考虑了缓发中子的中子扩散方程或中子输运方程而已。 我们已经知道如何求解扩散方程和输运方程，原则上也知道如何解时空动力学方程。 时空动力学的直接解法 对时间导数项进行