发明名称 一种自然循环蒸汽发生器的热传导建模与计算方法
摘要 本发明公开了一种自然循环蒸汽发生器的热传导建模与计算方法。基于圆筒壁无内热源一维导热微分方程通过将自然循环蒸汽发生器传热管和污垢中的温度分布表示为传热管一次侧壁温、二次侧壁温和污垢外表面温度的函数,并假定动态过程中传热管和污垢处于准稳态过程,将求解偏微分方程的问题转换为求解常微分方程;另外,热传导和流体热工水力耦合时采用了松耦合计算方法,在保证模型精度的基础上,大大减少了仿真计算的工作量,实现了自然循环蒸汽发生器热传导的快速计算,有利于自然循环蒸汽发生器乃至核动力装置高精度实时仿真计算的实现,对于指导核电站机组的安全、可靠运行具有重要意义。
申请公布号 CN104091036A 申请公布日期 2014.10.08
申请号 CN201410373080.5 申请日期 2014.07.31
申请人 中国船舶重工集团公司第七一九研究所 发明人 刘现星;陈保同;刘英;刘志宏;王长杰;郭研;范迎九
分类号 G06F17/50(2006.01)I;G06F19/00(2011.01)I 主分类号 G06F17/50(2006.01)I
代理机构 武汉东喻专利代理事务所(普通合伙) 42224 代理人 李佑宏
主权项 一种自然循环蒸汽发生器的热传导建模与计算方法,其特征在于,包括如下步骤:(1)n=0,初始化当前时刻n的传热管一次侧壁温<img file="FDA0000546797450000011.GIF" wi="92" he="74" />传热管二次侧壁温<img file="FDA0000546797450000012.GIF" wi="66" he="76" />和污垢外表面温度<img file="FDA0000546797450000013.GIF" wi="91" he="71" />(2)根据当前时刻n的传热管一次侧壁温<img file="FDA0000546797450000014.GIF" wi="60" he="70" />和传热管二次侧壁温<img file="FDA0000546797450000015.GIF" wi="89" he="71" />计算得到当前时刻n的传热管定性温度<img file="FDA0000546797450000016.GIF" wi="902" he="173" />由当前时刻n的传热管定性温度<img file="FDA0000546797450000017.GIF" wi="64" he="79" />计算得到当前时刻n的传热管导热系数<img file="FDA0000546797450000018.GIF" wi="97" he="71" />传热管密度<img file="FDA0000546797450000019.GIF" wi="70" he="71" />和传热管比热<img file="FDA00005467974500000110.GIF" wi="104" he="77" />根据当前时刻n的传热管二次侧壁温<img file="FDA00005467974500000111.GIF" wi="56" he="77" />和污垢外表面温度<img file="FDA00005467974500000112.GIF" wi="88" he="71" />计算得到当前时刻n的污垢定性温度<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mover><mi>T</mi><mo>&OverBar;</mo></mover><mi>cr</mi><mi>n</mi></msubsup><mo>=</mo><msubsup><mi>T</mi><mn>2</mn><mi>n</mi></msubsup><mo>+</mo><mrow><mo>(</mo><msubsup><mi>T</mi><mn>3</mn><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>T</mi><mn>2</mn><mi>n</mi></msubsup><mo>)</mo></mrow><mrow><mo>(</mo><mfrac><msubsup><mi>r</mi><mn>3</mn><mn>2</mn></msubsup><mrow><msubsup><mi>r</mi><mn>3</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>r</mi><mn>2</mn><mn>2</mn></msubsup></mrow></mfrac><mo>-</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mo>&CenterDot;</mo><mi>ln</mi><mrow><mo>(</mo><msub><mi>r</mi><mn>3</mn></msub><mo>/</mo><msub><mi>r</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA00005467974500000113.GIF" wi="912" he="184" /></maths>由当前时刻n的污垢定性温度<img file="FDA00005467974500000114.GIF" wi="70" he="80" />计算得到当前时刻n的污垢导热系数<img file="FDA00005467974500000115.GIF" wi="80" he="72" />其中,r<sub>1</sub>为传热管一次侧半径,r<sub>2</sub>为传热管二次侧半径,r<sub>3</sub>为污垢外半径;(3)利用当前时刻n的传热管一次侧壁温<img file="FDA00005467974500000116.GIF" wi="99" he="64" />传热管二次侧壁温<img file="FDA00005467974500000117.GIF" wi="90" he="62" />传热管导热系数<img file="FDA00005467974500000118.GIF" wi="91" he="75" />传热管密度<img file="FDA00005467974500000119.GIF" wi="98" he="71" />传热管比热<img file="FDA00005467974500000120.GIF" wi="72" he="79" />和污垢导热系数<img file="FDA00005467974500000121.GIF" wi="93" he="70" />计算得到下一时刻n+1的传热管一次侧壁温<img file="FDA00005467974500000122.GIF" wi="123" he="79" />传热管二次侧壁温<img file="FDA00005467974500000123.GIF" wi="88" he="77" />和污垢外表面温度<img file="FDA00005467974500000124.GIF" wi="122" he="78" />(4)n=n+1;(5)重复执行步骤(2)至步骤(4),得到传热管一次侧壁温、传热管二次侧壁温和污垢外表面温度随时间的变化。
地址 430064 湖北省武汉市武昌区中山路450号