A Pin by Pin transient analysis method for PWR steps are as follows: 1, the core is calculated as the initial state, the initial point of transient analysis; 2, on the next time step, considering the external disturbance information, at the same time using the neutron flux density on a time step update section; SP3 state the coefficients of the equation and the SP3 equation source term; step tn+1 3, using the exponential function expansion nodal method for space Pin by Pin numerical solution of the discrete time, the iterative convergence is obtained after the new time step gate element of zero order and two order neutron flux density; 4, update the delayed neutron precursor nuclear concentration; 5, determine the current time step is the last time step, if not, repeat the 2 4; if so, stop the calculation, the output transient parameters. Compared with the method of component homogenization, the invention does not need to consider the error of power reconstruction, and can only obtain zero order neutron flux density relative to the diffusion approximation, and more core information is reserved. At the same time, the solution time has not increased significantly.
【技术实现步骤摘要】
一种Pin-by-Pin分析压水堆堆芯瞬态的方法
本专利技术涉及核反应堆设计和反应堆物理计算领域,具体涉及一种Pin-by-Pin分析压水堆瞬态的方法。
技术介绍
在反应堆运行过程中,需要实时掌控堆内的真实情况,以便于操作员可以及时进行反应性调节,避免事故的发生,或在事故发生后控制事故的范围和等级,保障工作人员与周围群众的生命财产安全。这需要研究一种实时且准确分析压水堆堆芯瞬态过程的方法,即中子动态学计算。堆芯与事件相关的现象有很多种,根据其过程时间常数的不同大致可以分成三类:短时间瞬态、中等时间瞬态和长时间瞬态。短时间瞬态持续时间通常为毫秒到秒量级,主要是由于人为地或者事故地引入反应性而导致的中子通量密度迅速改变,比如瞬态实验,启动,停堆过程等。中等时间瞬态通常为持续数小时或者数天的过程,主要是热中子堆中裂变产物(Xe135,Sm149等)的积累、消耗以及β衰变所导致的瞬态过程。长时间瞬态的时间跨度达到数个月甚至数年,主要是由于易裂变产物的积累和燃耗,和可裂变核素的积累、燃耗以及衰变。用确定论方法直接求解中子输运方程是很耗时的,无法满足实时性要求,因此传统方法都是求解角度P1近似的中子扩散方程,空间上通常采用组件均匀化思想,先求得组件均匀化的中子通量密度,再通过功率重构得到组件内部的功率分布。但这样求得的中子通量密度和功率分布往往与实际的真实结果误差较大,特别是对于局部的非均匀效应。SP3近似方程是相对于中子扩散方程更加精确的中子输运方程形式,通过空间上Pin-by-Pin(以栅元为典型网格大小)求解SP3近似方程可以获得每个栅元内部的中子通量密度和功率 ...
【技术保护点】
一种Pin‑by‑Pin分析压水堆瞬态的方法,其特征在于:步骤如下:步骤1:计算堆芯初始状态,作为瞬态分析的初始点;得到每个栅元内部的零阶和二阶中子通量密度和修正后的多群微观总截面,微观散射截面和微观裂变截面;具体包括如下步骤:1)读取输入卡片中堆芯几何参数进行几何建模,读取核素的原始宏观总截面,宏观散射截面,缓发中子份额,衰变常数,能群信息,中子能谱,单次裂变产生中子数,读取瞬态时间步长信息和扰动信息;2)基于输入参数和堆芯几何计算初始的缓发中子先驱核浓度、SP3稳态方程的系数和源项;计算公式如下所示:
【技术特征摘要】
1.一种Pin-by-Pin分析压水堆瞬态的方法,其特征在于:步骤如下:步骤1:计算堆芯初始状态,作为瞬态分析的初始点;得到每个栅元内部的零阶和二阶中子通量密度和修正后的多群微观总截面,微观散射截面和微观裂变截面;具体包括如下步骤:1)读取输入卡片中堆芯几何参数进行几何建模,读取核素的原始宏观总截面,宏观散射截面,缓发中子份额,衰变常数,能群信息,中子能谱,单次裂变产生中子数,读取瞬态时间步长信息和扰动信息;2)基于输入参数和堆芯几何计算初始的缓发中子先驱核浓度、SP3稳态方程的系数和源项;计算公式如下所示:式中:g′——碰撞前中子能群;t0——初始时间;x——中子位置;Sf(t0)——初始时刻的裂变源;ν(x)——单次裂变产生中子数;Σf,g′(x,t0)——初始时刻在x处的第g′群中子宏观裂变截面;ν(x)Σf,g′(x,t0)——初始时刻在x处的第g′群中子宏观产生截面;——初始时刻的0阶第g′群中子通量密度矩;式中:i——缓发中子组数;Nd——缓发中子总组数;Ci(x,t0)——初始时刻在x处的第i组缓发中子先驱核浓度;λi——第i组缓发中子衰变常数;βi(x)——在x处的第i组缓发中子份额;式中:Sd——缓发中子源项;χd,g,i(x)——在x处的第i组缓发中子能谱;式中:D1——等效扩散系数;Σt,g(x)——在x处的第g群中子宏观总截面;Σr,g(x)=Σt,g(x)-Σs,g→g(x)公式(5)式中:Σr,g(x)——在x处的第g群中子宏观移出截面;Σs,g→g(x)——在x处的第g群中子宏观自散射截面;式中:Sg——等效中子源;μ0——初始角度;g——碰撞后中子能群;Σs,g′→g(x,μ0)——在x处角度为μ0的中子从第g′群散射到第g群的宏观截面;keff——有效增殖系数;χg(x)——第g群中子能谱;3)利用指数函数展开节块方法进行Pin-by-Pin稳态计算,得到初始时间步每个栅元内部的零阶和二阶中子通量密度和有效增值系数;SP3稳态方程如下所示:式中:——初始时刻的2阶第g′群中子通量密度矩;4)利用有效增值系数修正原始截面,得到更新后的截面参数,如公式所示:步骤2:对下一个时间步,考虑外部扰动信息,同时利用上一个时间步的...
【专利技术属性】
技术研发人员:曹良志,谢伟华,李云召,赵传奇,
申请(专利权)人:西安交通大学,
类型:发明
国别省市:陕西,61
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。