发明名称 模拟启堆时堆芯流道内冷却剂温度变化的新型计算方法
摘要 本发明属于解对流扩散方程的领域,特别涉及一种具有快速计算特点的模拟启堆时堆芯流道内冷却剂温度变化的新型计算方法。该数值计算方法的主要思想是:1、利用一个新的未知数将微分方程降阶;2、通过积分推导出相邻两节点处温度的关系;3、推导出任意节点与边界处温度的关系;4、利用边界条件解出各节点处温度。所述方法能够进行全局求解,在保证计算精度的同时,具有计算量小、占用计算机内存小的特点,能够快速得到流道内冷却剂温度分布情况。
申请公布号 CN106610917A 申请公布日期 2017.05.03
申请号 CN201610001360.2 申请日期 2016.01.05
申请人 华北电力大学 发明人 黄美;赵媛媛;欧阳晓平;汤建楠;张志俭
分类号 G06F17/13(2006.01)I;G06F17/12(2006.01)I;G06Q50/06(2012.01)I 主分类号 G06F17/13(2006.01)I
代理机构 代理人
主权项 模拟启堆时堆芯流道内冷却剂温度变化的新型计算方法是利用新引进的未知数对一维对流扩散方程的对流项进行降阶,其特征在于,将一个含有二阶导数一个变量两个自变量的偏微分方程变成一个只含有一阶导数两个变量两个自变量的偏微分方程组;模拟启堆时堆芯流道内冷却剂温度变化的新型计算方法的步骤如下:(1)简化压水堆堆芯冷却剂流道中冷却剂的对流传热为一维问题,认为燃料棒对冷却剂的加热为对流扩散方程中的源项;(2)利用新引进的未知数对一维对流扩散方程的对流项进行降阶,将一个含有二阶导数一个变量两个自变量的偏微分方程变成一个只含有一阶导数两个变量两个自变量的偏微分方程组;一维非线性对流扩散方程的表达式如下:<img file="845209dest_path_image001.GIF" wi="236" he="45" />引入变量v:<img file="174559dest_path_image002.GIF" wi="61" he="41" />得到方程组:<img file="373459dest_path_image003.GIF" wi="202" he="88" />(3) 将两方程联立得到方程组,对两个等式分别关于时间和空间进行积分以得到两变量在时间及空间上的关系;确定时间步长<img file="776759dest_path_image004.GIF" wi="17" he="9" />,以k,k+1,k+2……表示各个节点 ,利用下标k,k+1,k+2……表示节点处温度和温度梯度的空间位置,利用上标0,1,2……t表示各节点处温度和温度梯度的时间维度;这样就通过积分得到了关系式:<img file="353146dest_path_image005.GIF" wi="124" he="45" />其中,<img file="107475dest_path_image006.GIF" wi="44" he="21" />;<img file="906804dest_path_image007.GIF" wi="20" he="24" />表示<img file="430189dest_path_image008.GIF" wi="25" he="25" />,<img file="860033dest_path_image009.GIF" wi="22" he="24" />表示<img file="649129dest_path_image010.GIF" wi="25" he="25" />;<img file="189832dest_path_image011.GIF" wi="286" he="88" /><img file="833303dest_path_image012.GIF" wi="284" he="88" /><img file="230786dest_path_image013.GIF" wi="436" he="50" /><img file="959708dest_path_image014.GIF" wi="313" he="66" />通过上述关系式即可得到空间任意一点与<img file="350106dest_path_image015.GIF" wi="72" he="22" />之间的关系。
地址 102206 北京市昌平区回龙观北农路2号