工程热力学课程设计(论文)模板_V2.0

1、工程热力学课程设计（论文）学号 姓名 院系 专业 完成日期 授课教师得分要求：1. 论文内容（任选其中一类） ：a）新能源/可再生能源开发与利用现状综述（论文类）；b） 国防技术（航空、航天、航海）领域中的新兴动力技术（论文类）；c）热力学研究进展（论文类）；d）节能减排装置/系统设计（设计类）。2. 不能从网上下载，抄袭；3. 13 周周五（本学期最后一次课）交论文，需同时提交打印版（需装订）和电 子版（Word文档，Word2000以上版本）各1份。提交电子版时，文件题目 格式：学号_姓名，如“00 1张_ 三”，由各班负责人收齐后提交。不按要求提交 论文按无效提交处理。设计类参考格式1总

2、体要求全文控制在8页A4纸以内，设计作品类应包括下列内容：作品背景、设计 制作中解决的关键技术问题的描述、 作品或模型的图片 、创新特色、 预计应用前 景等、参考文献等。采用 word 2000及以上版本编排。2页面要求A4页面。页边距：上25mm，下25mm，左、右各20mm。正文采用小四号 字体，标准字间距，单倍行间距。不要设置页眉，页码位于页面底部居中。3图表要求插图按序编号，并加图名（位于图下方） ，采用嵌入型版式。图中文字用小 五号宋体，符号用小五号 Times New Roman （矢量、矩阵用黑斜体）；坐标图的 横纵坐标应标注对应量的名称和符号 /单位。表格按序编号，并加表题（位

3、于表上方） 。采用三线表，必要时可加辅助线。 4字号、字体要求（仅作参考）家用电器节能系统设计说明书学号：xxx；班号：xx>设计者：XXX（空一行） 作品内容简介通过实验设计了一套家用电器节能系统 （400600字以内）（空一行）1 研制背景及意义2 设计方案2.1 电器控制电器部分采用自动控制比较容易实现，考虑到电器元件易发热等问题 2.2 机械部分机械部分设计如图 1 所示， 设计时考虑的主要问题：3 理论设计计算4 工作原理及性能分析完成制作后，作品实物外形照片见图 9。5 创新点及应用1) 适用于不同类型家庭电器。2) 操作和控制简便，容易地使用它。3) 。在全国大中城市， 家

4、用电器普及数量很多， 所有电器都有待在节能措施实施 改进，因此应用前景很广。图 1 家用节能机构原理图正文中表示物理量的符号，表示点、线、面的字母均用 Times New Roman 斜体；表示法定计量单位、词头的符号、函数等，化学元素符号均用 Times NewRoma n正体。(空一行)参考文献1 xxx ，xxx 家用电器节能现状和发展节能机械，2001，23(3)：275-2792 xxx 节能技术基础 xxxxxx 出版社， 1996： 15-473 xxxxx ， xxxx xxx ，xxx 译机器人操作的数学导论 xxxx 出版社， 1998： 11-674 Lee H Y, R

5、einholtz C F. Inverse kinematics of serial-chain manipulatorsJ. ASME Journal ofMechanical Design. 1996, 118(3): 396-404论文类参考格式1总体要求不少于3000字，全文控制在8页A4纸以内。论文应包括下列内容：摘要、关键词、引言、综述、结论、参考文献等。采用 word 2000 及以上版本编排。2页面要求A4页面。页边距：上25mm，下25mm，左、右各20mm。正文采用小四号 字体，标准字间距，单倍行间距。不要设置页眉，页码位于页面底部居中。3图表要求插图按序编号，并加图名（位

6、于图下方） ，采用嵌入型版式。图中文字用小 五号宋体，符号用小五号 Times New Roman （矢量、矩阵用黑斜体）；坐标图的 横纵坐标应标注对应量的名称和符号 /单位。表格按序编号，并加表题（位于表上方） 。采用三线表，必要时可加辅助线。 4论文结构（仅作参考）考虑温度反馈的点堆中子动力学求解学号：XXX；班号：XXX；作者：XXX摘要： 本文分别用采用绝热模型的幂级数法和 Simulink 动态仿真的方法，对考虑简单温度 反馈的点堆中子动力学方程的进行求解。 得到了考虑简单温度反馈时， 反应性、 反应堆功率、 缓发中子先驱核浓度， 反应堆温度随时间的变化情况， 并通过对引入小阶跃反应

7、性情况下计 算结果的分析得出了绝热模型的适用范围。关键词： 温度反馈、幂级数法、 Simulink 仿真、绝热模型1 引言在核反应堆实际的运行过程中， 反应性是不断变化的， 影响其变化的因素有很多， 在中 子动力学分析中注重的是反应堆在启动、 停堆和功率调节过程中短期内中子密度随着时间的 变化关系， 因此可以不考虑燃料的燃耗和裂变产物的积累对反应性的影响， 只需要考虑温度 对反应性的影响。 考虑温度反馈可以对反应堆功率的瞬态行为和其它反应堆堆芯系统变量进 行估计，这些变量都是紧密耦合的，考虑温度反馈的点堆中子动力学方程为非线性方程组， 求解比较困难。 以前尝试的求解方法多数只适用于瞬发临界。

8、目前对于考虑温度反馈且引入 大阶跃反应性（反应堆处于超瞬发临界）的情况，主要采用诺德黑姆-福赫斯模型 1 进行处理，可得到较好的结果， 但其只针对考虑单组缓发中子效应的问题。比较常见的反应性反馈的数学模型还有准稳态模型和绝热模型 【2，其中绝热模型因其形式简单， 能大大降低方程求 解的难度，得到了广泛地应用。但由于其忽略了反应堆的热量损失， 在一些情况下并不适用。本文在求解考虑温度反馈的点堆中子动力学方程时就采用了绝热模型，并应用求解点堆中子动力学方程综合性能较好的幕级数法3进行编程求解，并对所得结果进行分析。此外，还在不应用简化模型的情况下，对考虑温度反馈的点堆中子动力学方程中的各变化量进行

9、 Simulink动态仿真，得到了反应性、反应堆功率、缓发中子先驱核浓度，反应堆温度随时间 的变化情况。2采用绝热模型的幕级数法考虑温度对反应性的反馈作用时，反应性二不但是时间的函数，还是中子通量密度或中子密度n(t)的函数，因此考虑 6组缓发中子效应的点堆中子动力学方程变为非线性方程组。 不考虑外加中子源时，方程组具有如下形式：dn tdt6dG (t) dt dT(t) dt訓)iCi(t)Kn (t) T(t) Tc(t)o T(t) T(t)i 1,丄，6Ct式中，T(t)为核反应堆温度，K ; T(0)为T(t)的初始值，K ; Tc为堆芯平均冷却剂温度(不随时间变化)，K ; 1/

10、 丫可被理解为热能由燃料传递给冷却剂的平均时间，S；二为初始时引入的阶跃反应性；二为反应性温度系数飞'，K-1。采用绝热模型对方程组进行简化，认为热量传递的时间常数 泳浙刘与功率偏离的时间相比很大，可忽略热量损失，即-弋。则有呼 Kn(t)由和式可得:丁 Kn(t)将方程(1)、(2)、(6)改写为如下矩阵形式：警 M (t)y(t)dt式中，y(t)为8维列向量，M(t)为(8) (8)矩阵：y(t)n(t)C(t),C2(t)丄 3), (t)T(t)12L60A1A10L00M(t)2A02L00MMMOMM6A00L60K00000利用幕级数作为基函数对式在一个时间步长内展开为

11、:(8)y(t) AA1 LAntnAntnn 0式中，叫-、I：、八：均为8维列向量。 由初始条件，根据（8）式可得式中,y(0)A0y'(0)A1Ly(n)(0)(n!)Ann 1y(n)(t)cn1M (i)(t) y(ni 1(t)i 01 d (t)OM' (t)A dtOO另由（39）式可得(9)(10)M1 d(t)1Ot=0时，由方程（9）和（10）可得(ii)M(i 1)(0)An因此，要想求得丨.，只需求解由（11）式分别取n =1,2, ,k组成一个线性方程组即可（n 为整数，是幕级数的阶数），再将I代入（4-22）式即可求得该步长内点堆中子动力学方程组

12、的解，然后将所得末点值作为下一个时间步长的初始值重新计算，如此迭代循环可以得到整个时间段内的解。当迭代计算进行第一步时，M' (0),L , M (i)(0)均为零矩阵，则(11)式可化简为1AnM (0)An 1(12)n1 i 1L oA hOO进行后续迭代计算时，则有4:M' (0)且M(0)(亡2)均为零矩阵，则(11)式可化简为：1An -(M (0)An1 M'(0)An2)(n 1)(13)n其中，A 10,0,0,0,0,0,0,0 T。计算时采用参考文献5中给出的反应堆特性参数及初始条件。平均中子代时间圧5 X10-5s,缓发中子总份额3=0.0065

13、 , 6组缓发中子份额及先驱核衰变常数如表1所示。=-0.00306 3 K=7.6886 氷0-5K/J，初始条件： 处=0.043 3 no=O.O1W。表1六组缓发中子份额及先驱核衰变常数组数i123456份额3(沁0-3)0.221.421.272.570.750.27衰变常数Ms-1)0.01240.03050.1110.3011.143.01将计算得到的100min内方程组的解随时间变化的部分结果列于表2中。参考文献5中对考虑简单温度反馈的点堆中子动力学方程求解时没有采用简化模型，其计算结果还是比较准确的，将其作为标准参考值与自己用幕级数法采用绝热模型编程求得的结果进行比较。 参考

14、文献5中给出的100min内方程组的解随时间变化的部分结果见表3。由表2和表3可以看出，在取四位有效数字的情况下，前30min内自己用幕级数法编程求得的结果与参考文献5中给出的结果基本完全相同，而30min后的计算结果与文献5中的结果开始产生明显的差距。可见绝热模型在阶跃反应性引入的初始一段时间还是非常适用的。为了更清楚地了解采用绝热模型的幕级数法的计算精度，求取其计算结果与参考文献5中给出结果的相对误差，得到相对误差随时间变化的曲线如图1至图3所示。表2幕级数法计算结果t(min)p($)n(W)C1C2C3C4C5C600.0430.013.5489.3112.2881.7080.1316

15、0.01794100.042990.103328.2385.7822.8817.431.3550.1851200.042940.9514260.1790.1210.7160.512.481.705300.042458.6522371719619171460113.515.50400.0382070.421.975 X045.926 X041.566 X041.190 :W4924.0126.2500.01426278.88.930 X042.503 X056.322 X044.745 1043665500.060-0.024194.37.876 X041.919 X054.521 X043.3

16、39 1042561348.970-0.039146.242.084 X044.723 X041.085 X047968609.883.0280-0.04238.5053925874419981466112.215.2790-0.04291.493692.01537350.8257.419.702.681100-0.0430.2600120.6267.661.0944.833.4290.4668t(min)p($)n(W)C1C2C3C4C5C600.0430.013.5489.3112.2881.7080.13160.01794100.042990.103328.2385.7822.8817

17、.431.3550.1851200.042940.9515260.1790.2210.7160.512.481.705300.042458.6542372719719171460113.515.51400.0382570.511.976 X045.933 X041.568 1041.191 104925.2126.4500.01468282.29.018 X042.531 X056.396 1044.802 1043709506.060-0.0224205.28.243 X042.019 X054.768 1043.524 1042704368.370-0.035553.992.382 X04

18、5.469 X041.264 1049295711.996.9380-0.035612.0953621.225 X0428302082159.521.7190-0.03282.91712722936681.5501.938.466.237100-0.02950.7945339.2793.7185.2136.610.471.426表3参考文献5中给出的计算结杲0.5000.40.30.20.110 203040506070t(min)8090100差误对相图1反应性的相对误差随时间变化O/T10203040506070 8090 1(t(min)图2功率的相对误差随时间变化7 6 5 o o O

19、4 3 2 o o O差误对相OVyrC1C2C3C4C5C6,6 5 4 3 20 0 0 0 0 差误对相0.1102030406070809010050t(min)图3缓发中子先驱核浓度的相对误差随时间变化由图1至图3可以看出由采用绝热模型的幕级数法计算得到前40min之内方程的结果与参考文献5中的结果间的相对误差非常小，但之后随着时间的增长相对误差明显增大， 100min时反应性的相对误差已达到约0.46，反应堆功率和缓发中子先驱核浓度的相对误差已达到约0.68。分析误差产生的主要原因应该是计算时采用了绝热模型这一简化处理，由于初始时引入了一个较小的阶跃反应性，中子的增殖或功率的增长要

20、在较长时间内完成，在开始的较短时间内可以忽略反应堆的热量损失，这时绝热模型是适用的，因此在采用绝热模型计算的前一段时间的结果还是非常准确的，但随着时间的增长，反应堆有足够的时间将热量传递出去， 此时采用绝热模型就会使求得的温度随时间变化比实际情况下快，从而使反应性的变化也快于实际情况。3基于Simulink的动态响应仿真创建动态系统 Simulink模型的一般步骤为：首先建立理论数学模型，列出描写系统动力学的全部方程。然后，根据理论数学模型，在Simulink模块库中选择所需的模块并用鼠标勾画各模块间的信号连线， 建成仿真系统。最后，根据所研究的具体问题设置各模块的参 数和仿真时间并进行仿真和

21、调试得到满足需要的Simulink模型。由于参考文献5中给出的考虑温度反馈的方程并未采用任何化简模型，而是较充分地考虑了反应堆的热量传递等情况，因此，采用参考文献5中所列出的考虑温度反馈的点堆中子动力学方程组作为理论数学模型。忽略外加中子源项，并考虑引入阶跃反应性的情况， 其形式如下：dTdt1cpn(t)AHD0 (T(t) Tc) /T(t)0 (14)dn(t) dt-(t n(t)iCi(t)Ai 1(15)dCi(t)Ln(t) iCi (t)ii1,2 丄，6(16)dtA(t)0T (t)Tc(17)其中A 1.649Kc0 - 75( fg?p/ )0.25式中，H、D分别为反

22、应堆的高度、直径，m； T为反应堆的温度，K； Cp为反应堆的热容，J/K; £为裂变能转化为热能的分数；Tc为冷却剂的温度，K ;Kc为冷却剂传热系数，Jm-1s-1K-1 ； P为冷却剂密度，kg/m3;卩为冷却剂的质量流速，kg m-1 s-1;Cp为冷却剂的比热容，J kg-1 K-1。图4 Simulink系统仿真模型采用参考文献5中给出的各物理量的值对各模块的非结构参数进行设置。其中反应堆 特性参数及初始条件已在上一节中列出，其它各量的值为：f1 ,A=17.52 ,H=0.23m, D=0.20m ,Cp=1800 tD2/4 H 1000, Tc=20=-0.0030

23、6 伶在用Simulink对点堆中子动力学方程组进行求解时可能会涉及到刚性问题，为保持求 解的稳定性，选择专门用于处理刚性问题的变步长解法ode456,7。设置仿真时间为 6000s,仿真结果可随时由 Scope模块显示出清晰图形以供分析。100min内求得的各物理量随时间变化如图5至8所示。0.05350250200150100mi仿真结果I1K!董参考文献5中的结1LTj-斗参JLXR<果0.040.030.020.01-0.01-0.02-0.038090p 010203070405060t(min)图5反应性随时间的变化100仿真结果H参考文献5中的结果1气4/)/1ifJX!&

24、#39;Il1300102030708090405060t(mi n)100图6功率随时间的变化504540353025202030708090100405060t(mi n)X kf / / /10图8温度随时间的变化为了检验仿真结果的准确性，图5和图6中不但给出了仿真的结果， 同时还给出了参考文献给出的结果，从图中可以看出，仿真得到的反应性和功率的结果与参考文献5中采用Fortran语言编程求得的结果吻合的非常好，这说明基于Simulink的动态响应仿真具有很好的计算精度，而且可以避开繁琐的程序编写，具有快捷、高效的特点。由图5至图8可以看到，在引入小的正阶跃反应性后反应堆功率（中子密度）

25、和缓发中子先驱核浓度开始缓慢增长，因而反应堆堆芯温度也缓慢升高，由于存在温度反馈且温度系数为负数，温度升高导致了反应性减小。当反应性减小到零，中子密度和缓发中子先驱核 浓度也就不再增加达到了最大值。此时，由于裂变反应释放的热量仍大于堆芯传出的热量， 反应堆温度继续上升， 反应性继续减小，其值将小于零中子密度和缓发中子先驱核浓度开始 减小。直到由于裂变反应释放的热量仍小于堆芯传出的热量，反应堆温度开始下降，反应性开始增大。4 结论在实际反应堆运行过程中， 反应性一般是随时间变化的。 在中子动力学分析中注重的是 短时间内中子密度随着时间的变化关系， 因此主要考虑温度对反应性的影响。 考虑简单的温

26、度反馈之后， 点堆动力学方程变为非线性方程组， 直接编程求解比较困难。 本文采用比较简 单且常用的绝热模型对方程进行化简， 并应用幂级数法编程计算得到结果， 分析该结果得到 在引入较小的正阶跃反应性后， 中子的增殖或功率的增长需要较长的时间， 在开始的较短时 间内可以忽略反应堆的热量损失， 这时应用绝热模型可以得到非常准确的结果的。 但随着时 间的增长， 反应堆有足够的时间将热量传递出去， 此时采用绝热模型后求得各物理量随时间 的变化速度明显快于实际情况，绝热模型不再适用。基于 Simulink 的动态响应仿真不但可 以避开繁琐的程序编写，而且具有很好的计算精度。从仿真结果看出，当温度系数为负数， 反应堆出现反应性扰动时， 堆内各物理量能通过负反馈效应自行调整达到另一平衡态， 反应 堆具有内在的稳定性，这对反应堆的调节和运行安全都具有重大意义。参考文献1 蔡章生 . 核动力反