并联多通道内核热耦合流动不稳定性的研究毕业论文.doc

并联多通道内核热耦合流动不稳定性的研究毕业论文 目录 1 绪论 1 1.1 流动不稳定性及其危害 1 1.2 研究现状 2 1.3 研究内容 4 2 点堆模型和核热耦合模型 6 3 并联多通道理论模型 9 3.1 单相水热工水力微分方程 10 3.2 两相水热工水力微分方程 11 3.3 压降模型 14 3.3.1 入口段 14 3.3.2 加热段 14 3.3.3 上升段 15 3.4 无量纲化 16 3.5 并联通道模型 18 4 数值方法 20 4.1 系统的刚性 20 4.2 Gear 方法 24 4.2.1 简单迭代法 24 4.2.2 Newton迭代方法 25 4.2.3 Gear方法 26 4.2.4 Gear算法的阶与步长的控制 28 5 程序计算结果及分析 33 5.1 不考虑核热耦合时的计算结果 34 5.2 核热耦合时的计算结果 37 5.3 计算结果对比及分析 40 6 结论与展望 42 参考文献 43 附录 45 致谢 78 绪论 流动不稳定性及其危害 在加热的流动系统中，如果流体发生相变即出现两相流时，流体以非均匀形态所出现的大的体积变化可能导致流动的不稳定性。这里所说的流动不稳定性，是指在一个质量流密度、压力降和空泡之间存在着耦合的两相系统中，流体受到一个微小的扰动后所产生的流量漂移或者以某一频率的恒定振幅或变振幅进行的流量振荡[]。这种现象与机械系统中的振动很相似：质量流密度、压降和空泡可以看作是机械系统中的质量、激发力和弹簧，在这中间，质量流密度和压降之间的关系起着重要的作用。流动不稳定性不仅在热源有变动的情况下会发生，而且在热源保持恒定的情况下也会发生。 流动不稳定性在反应堆领域广泛存在，流动不稳定性会严重影响反应堆的安全运行，这是由于： 部件可能遭受有害的强迫机械振动，机械振动和局部热应力周期性变化会导致疲劳破坏； 会引起控制问题。对于流体冷却的反应堆，当冷却剂兼作慢化剂时尤为重要，还会引起核反应性变化耦合反馈效应； 会影响局部传热特性，或者引起沸腾危机提前出现，及临界热流密度降低。 两相流不稳定性大致可以分为两大类：静力学不稳定性和动力学不稳定性。 静力学不稳定性是非周期性的改变系统的稳态工作运行点，它的基本特征是系统在经受一个微小扰动后，会从原来的稳态工作点转变到另一个不相同的稳态工作点运行。这类不稳定性是由于系统的流量和压降之间的变化、流型转换或传热机理的变化所引起的。 动力学不稳定性是周期性地改变系统的稳态工作状况，这里惯性和反馈效应是制约流动过程的主要因素。它的基本特征是当系统经受某一瞬间的扰动时，在以声速传播的压力扰动和以流动速度传播的流量扰动之间的滞后和反馈作用下，流动发生周期性振荡。这类不稳定性的产生主要是由于系统的流量、密度、压降之间的延迟与反馈效应，热力学不平衡性以及流型的转换等原因引起的。本文研究的流动不稳定性为动力学不稳定性。 研究现状 从上可见两相流动不稳定性的研究对于反应堆的安全运行以及事故安全分析都具有十分重要的意义。关于这方面的研究已经产生了很多成果。然而，大部分这方面的研究，都只是集中于单一的流体动力学不稳定性，而没有考虑到核热耦合的相互作用。这种作用通过中子(即反应性)，把热流密度与空泡份额联系在一起，结果也就与流量联系起来了。1988年发生在美国的LaSalle核电站紧急停堆事故即是由于核热耦合产生不断的振荡产生的。LaSalle事故被认为是由耦合的核热水力不稳定性产生的，但它的机制还未被完全了解[]。 目前关于核热耦合流动不稳定性的研究主要有时域法跟频域法两种。这两种方法各有利弊。时域法得出的结果直接随时间而变，可以很直观地看出是否发生了流动不稳定。但是时域法需要占用大量的计算时间，而且还容易产生由于数值计算方法引入的计算振荡问题。频域法不产生各种各样的变量在振荡中表现的瞬时，结果因此不会产生大量的中间计算结果，缺点是不直观。 在实际的动力系统中，很多装置具有上下两个联箱间并联一组管道的结构，可以通过平均管分析得到关于整体的不稳定性，但是有一种情况却是系统总流量保持不变，而管道之间已经发生了脉动，这是一种管道之间的相互作用，非常隐蔽，更加危险，在近年来成为不稳定性研究中的一个热点和难点问题。一般认为其求解的边界条件应该是总流量不变和两端压降不变，但是尚未见到实验证明这一边界条件的合理性。在很多的模型求解中，一般仅仅将两通道之间的压降保持相等，而这一压降是变化的。在文献[]中提到的方法尽管保证了压降不变，但在同一个时间步长内两个通道的压降是不等的。这一矛盾在解方程组时无法避免，在这一点上需要更多的实验验证和理论研究。本课题的研究基于前一种方法，即是求解时保持各个通道之间的压降相等。 阮养强[]基于多通量频域法，针对多个通道相互并联的情况