关于反应堆棒束通道内流体的湍流特性的研讨.docVIP

关于反应堆棒束通道内流体的湍流特性的研讨 　　摘要：在当今资源不断缺乏的大环境下，各个国家都在大力发展核电，但由于核电站事故的严重性，就会发现核电安全显得尤为重要。要想保证核电站的安全就需要对反应堆棒束通道内的流体特性进行分析和研究，了解流体的速度场和压力场，从而使反应堆能够安全可靠的运行。 　　关键词：CFD 湍流特性 数值模拟 棒束通道 　　一、ANSYS CFD软件简介 　　CFD（Computational Fluid Dynamics计算流体力学）是计算机技术、数值计算技术、流体实验、计算机模拟实验为一体的模拟软件。基本原理是数值求解控制流动的微分方程，得出流场在连续区域上的离散分布，从而近似模拟出流动情况。 　　CFD模拟是利用计算机求解各种守恒控制偏微分方程组的技术。涉及流体力学（湍流力学）、数值方法乃至计算机图形学等多学科。且会因所要解决问题的不同，建立的模型方程与数值方法也会有所不同，如可压缩气体的亚音速流动、不可压缩气体的低速流动等。 　　二、棒束通道内的湍流流动理论基础 　　粘性流体的运动存在两种完全不同的流动状态：层流状态和湍流（紊流）状态。1883年雷诺通过实验方法观察到了这两种不同的流动状态。 　　层流流线之间没有相互交叉，呈现出一种十分有规律的分布，不存在横向的交混流动。湍流组是一种十分复杂的三维非稳态、带旋转的不规则流动[1]，不仅有轴向流动，还会伴有很强横向流动和漩涡流。 　　众所周知，流体的流动状态和雷诺数有很大关系，这是由雷诺通过进行大量的试验总结出由层流到湍流的转变是由Re决定的，当Re大于一定数值时，流动就由层流过渡到湍流。Rec一般取为2000。本文所研究的棒束通道内的流动流体的Re在2.85×104到1.64×105之间，同时存在横向交混流动，所以为强烈的湍流流动。 　　在湍流流动中流体的各种物理参数，都会随着时间和空间的变化不断的随机的改变。一般把湍流流动当成是由各种不同的涡旋叠合而成的流动状态，然而这些涡旋的尺寸及旋转轴的方向都是不确定的有很强的随机性。大尺度涡主要由流动的边界条件决定，其尺寸可以与流场的大小相比拟，主要是受惯性影响而存在，是引起低频脉动的原因；小尺度的涡主要由粘性力决定，其尺寸可能只有流场尺度的千分之一的量级，是引起高频脉动的原因[2]。 　　三、涡粘性模型 　　（1）标准k-ε模型 　　标准的k-ε模型是由湍动能k方程和ε的方程组成，这样便形成了标准k-ε模型，标准的k-ε模型在1972年是由Launder以及Spalding提出的[3]。 　　各向均同性的小漩涡的机械能转变为热能的速率，也就是紊流中单位质量的流体湍动能的耗散率，定义式由下式给出： 　　经分析得： 　　式中，Cμ为经验系数。 　　对于不可压缩的流体，约束k和ε的方程为： 　　经验系数在模型中的取值为： 　　由于其使用广泛并且很适用，标准k-ε模型在已经工程上成为使用最广泛的湍流模型。但对于其适用性仍需作如下说明： 　　1.在模型中的所有的相关系数一般是在特定条件下进行试验获得，一些系数可以说具有一定的代表性，但是需要仔细的对计算时所设置的特定条件查阅有关的资料进行合理的设置。 　　2.所建立的标准k-ε模型，就是为了计算流体在充分发展阶段的湍流特性，所以它就是为了计算大雷诺数湍流流动而产生的湍流计算模型，尤其是在流体远离壁面特定距离的湍流流动区域。 　　3.标准k-ε模型是在零方程模型和一方程模型的基础上做了总结与创新，通过了大量的科学研究和在工程实际中的应用，对标准k-ε模型进行了检验和应用的目的，但在强旋流动和弯曲壁面流动会产生一些误差。那是由于在标准k-ε模型中μt为各向同性。而实际中 μt为各向异性。 　　（2）SST k-ω湍流模型 　　SST k-ω模型是建立在BSL k-ω模型的基础之上[4]。Wilcox的k-ω湍流模型有缺点，由于对进口的ω值很敏感， ω值不同，所输出的结果也会不一致，且变化很大。我们并不希望这样，于是Menter选择了BSL k-ω模型，在壁表面采用k-ω模型，而在紊流发展可观的中心部位采用k-ω模型。这就存在将k-ε模型向k-ω模型转变的需求。k-ω模型乘以调和函数F1，k-ε模型乘以（1-F1），于是变成BSL k-ω模型。 　　在BSL k-ω模型和SST k-ω模型中，在壁表面采用k-ω模型，湍流充分发展区采用k-ε模型，这就需要确定壁面距离（Wall Distance），在这个距离处开始转变湍流模型。在CFX- Solve：窗口中也有一个监视“Wall Scale”的视图。 　　四、不同湍流模型的比较 　　为了分析出四种湍流模型的适用性，做出每个模型C截面处的速度流线图来分析轴向交混程度，如下