高等流体力学课程论文_第1页
高等流体力学课程论文_第2页
已阅读5页,还剩5页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、科目高等流体力学教师:叶建姓名:杨伟,陆祺,张鹏青4应殳事研究生课程考核试卷(适用于课程论文、提交报告)学号:20131002051,20131002021t,20131002026t专业:动力工程及工程热物理,核科学与技术,动力工程及工程热物理类别:学术,学术,学术上课时间:2013年9月至2013年11月考生成绩:卷面成绩平时成绩课程综合成绩阅卷评语:阅卷教师(签名)气冷堆燃料元件表面流场及温度场模拟1.本文主要研究内容国内外已经有许多学者利用计算流体力学(CFD)对模块式高温气冷堆球床内的流场、温度场做过大量研究,并取得了不少成果。这方面的研究涉及了燃料元件与壁面或燃料元件之间接触和间隙

2、情况下的传热和流动的特征,孔隙率对温度分布均匀性的影响等等。本文的主要研究内容是数值计算获得堆芯元件表面冷却剂氦气的速度及温度分布规律,并研究不同堆积方式对分布规律的影响。2.球床内流动与换热数学模型实际上,在球床对内,燃料元件的分布是随机不规则分布的,有些部分可能是正四棱锥堆积,有些部分是立方体形式堆积;当然,还有可能是其它的堆积方式。本文主要研究的是正四棱锥堆积于立方体堆积时对温度、速度分布的影响。因此,建立的物理模型由立方体和正四棱锥两种,采用的计算软件为Fluent。2.1物理模型在球床对内,球形燃料元件的不同堆积方式对氦气的流动及其温度分布会产生的影响不同。因此,我们需要建立模型模拟

3、以下情况:(1)球床内温度、压力和速度分布规律;(2)堆积方式变化时温度分布和流动的改变。基于这两方面,我们需要建立不同堆积方式的模型。图2.1正四棱锥堆积图2.2和图2.3分别为立方体堆积模型和正四棱锥堆积模型,其中可以明显看到,左边为燃料元件的分布,右边为流体区域。在两个模型当中,燃料球的尺寸均为60mm。图2.1中为正四棱锥堆积方式的平面尺寸图。第一层燃料球球心之间的距离为7.002cm,第一层与第二层燃料球球心的距离为3.3912cm。通过对以上两个简单的模型进行数值计算,可以获得球床内冷却剂氦气的速度及温度分布规律。通过对以上两个模型的计算结果比较和分析不同堆积方式对球床内温度分布及

4、氦气流动的影响。图2.2立方体堆积模型图2.3正四棱锥堆积模型图2.4非结构网格2.2-£模型为了使方程封闭,特别是在高雷诺数下的计算,应采用久-三标准模型,该模型是在湍动能k方程基础之上引入湍流脉动动能的耗散率三方程形成的,于1972年提出。其中k为湍动能,三为湍流耗散率,分别来自湍流动能方程式和耗散方程式。它们的计算公式如下:久=匚:.、;-;二(2.4)式中u为流体流速,I为湍流强度。(2.5)由量纲分析可以得粘性系数与湍动能k和耗散率三的关系式如下:=兀、(2.6)其中有:三=茫-二(2.7)式中二为无量纲常数,通常取0.09,k为湍动能,1为湍流尺度。2.3物理模型及求解器

5、设置在求解控制方程时采用Fluent软件,选用基于压力的求解器,物理模型采用非稳态模型,由于本文三维模型网格单元基本为多面体结构,因此选least-quarescellbased这种压力梯度来计算控制方程中的导数项。对于近壁面的处理方式采用标准壁面函数,其中压力-速度解耦一般使用SIMPLE和SIMPLEC,对非定常流动问题或者包含比平均网倾斜度更高的网格则使用PISO算法,基于本模型网格不规则且是非定常流动模型,因此选用PISO算法,其计算结果会更精确一些;但是这种算法相比较起SIMPLEC算法,可能收敛较困难,因此需要根据实际需要调整松弛因子,该因子减小有利于收敛,但是会增加收敛时间,所以

6、需谨慎使用。由于球床内的燃料元件堆积形成的空隙通道不规则,流体压力变化剧烈,PRESTO!主要用于高旋流,压力急剧变化流。因此压力基分离求解器的压力插值算法建议使用PRESTO!。这里可以先采用一阶方程进行初始迭代,再进行二阶迭代计算,若MUSCL算法能够满足收敛条件则亦用MUSCL算法进行更精确的计2.4边界条件设定及其实现在本模型当中,由于只取局部进行分析,设定氦气为在高压下的局部模型当中为不可压缩流体。根据HTR-10的运行参数知,运行压力为3MP,由于在堆芯取的局部模型压降相对于运行压力来说比较小,相对压力取0Pa时有利于分析。进出口温度分别为250°C和700°C

7、,取其之间的任意温度均可,为方便计算,取模型入口温度为600K。整个堆芯入口质量流量为4.3kg/s,根据模型入口截面的大小可以算出正四棱锥模型入口流量为0.03312kg/s,立方体模型的入口流量为0.02454kg/s。对于进出口的湍动能k和湍流耗散率三,可由上述公式计算得出。在流体流动过程当中,球体壁面温度分布不均匀,不能设定球体表面温度,第一类边界条件不适用,这里采用第二类边界条件。同时选用无滑移光滑表面。从HTR-10的运行参数可知,其热功率为10MW,燃料球总数量为27000个,尺寸为60mm,有以下公式可得出:反应堆总热功率燃料球总表面和=32747.93W/m2所以,设定壁面的

8、热流量为32747.93W/m2.3计算结果与分析本章根据第三章所建立的数学模型进行计算,得到固定床反应堆堆芯内冷却剂氦气的速度及其温度分布规律,并讨论了燃料球尺寸、孔隙率等对上述规律的影响。温度分布卫压力分布+J扰动分布;图3.1速度、温度、压力及扰动分布3.1正四棱锥堆积方式从整个模型的分析可以看出,我们已得到温度、速度甚至压力变化的大概规律,但是要想更好更具体地分析冷却剂氦气的温度、速度以及压力分布,我们还需要从单一燃料元件着手,研究其表面压力以及温度等的分布情况。在整个球床堆芯内,除球床周围边界的燃料元件外,几乎所有燃料元件的温度、速度和压力的分布规律基本相同。因此,取模型当中的单一燃

9、料元件来分析即可,在此模型当中,分别取两种情况下其中心的单一燃料元件对比分析。如图3.2所示,在元件间接触点附近出现温度的极大值点。在燃料元件末端也出现了较高的温度,这是由于在这个地方冷却剂流速很低,且有涡流形成,使该区域内冷却剂出现滞留甚至倒流,导致对流换热很差,从而使元件的局部温度很高。图3.2元件表面温度分布如图3.3所示,在Z轴-1.6cm和1.6cm的地方都出现了压力的极小值点,这些地方出现这种现象的原因是由于这些位置间隙通道较小,且形成类似扩压管结构,导致压力迅速降低而流速迅速升高。图3.3元件表面压力分布如图3.4所示,该图为元件表面沿冷却剂流动方向的速度分布图,可发现速度最高点

10、与压力的极小值点相对应。原因如上所述,均是由于元件间狭小缝隙通道处形成扩压管结构,在此处压力迅速降低,同时速度升高。VelacfyMagjitixie(mTs5-lMi-eiIDO-ClRgoofisSte-OO(QQM0图3.4元件表面速度分布因此可以得出结论,在燃料元件堆积形式不变,改变尺寸或者增大流速、热流密度等的情况下,冷却剂在球床内的流动的规律基本一致,分布形式大致相同,只是在形式或者程度上发生了变化。3.2立方体堆积方式当燃料元件堆积形式发生改变时,必然对对球床燃料元件温度、速度和压力分布和变化的影响。在不改变燃料元件尺寸的前提下,对其堆积形式做一些改变,分析这一单一变量对上述分布

11、的影响。把上述正四棱锥堆积方式变为立方体堆积形式,其分布如下:如图3.5所示,该图是模型中心的单一元件温度分布图。在坐标-3cm和3cm处出现了温度的最高值,而在0点附近也出现了温度分布的较高点。其中坐标-3cm和3cm处出现温度极值的原因有两个:一是此处位于上下两层燃料元件接触点;二是此处流速较低,与正四棱锥模型一样,有涡流形成,该区域内冷却剂出现滞留甚至倒流,对流换热很差。0点附近出现温度较高点的原因单一,就是燃料元件之间在此处相互接触,接触点附近流体流动较慢,对流换热较差。图3.5元件表面温度分布如图3.6所示,沿流动方向,压力在元件表面的变化与正四棱锥堆积形式的变化有很大区别。在立方体

12、排布形式当中,元件上半部分压降变化很慢,但到下半部分时,压力会迅速降低,这与上面所述的原理一样,在燃料元件相互形成的较小空隙处形成了类似扩压管的通道,导致这个位置压力突然降低。两者有一个共同点就是这个压力迅速降低的位置都是在与周围元件形成的较小空隙处。但是元件堆积方式不同,这个较小间隙位置就处于元件不同位置。图3.6元件表面压力分布如图3.7所示,速度在0点附近出现了最大值点,在-3cm和3cm处最小,几乎为零。与图3.13相对应分析可知,在-3cm和3cm处为燃料元件接触点,也是燃料元件在流动方向的最低和最高位置,在这些位置会出现涡流甚至倒流,速度最低。而在0点附近周围元件形成了类似于扩压管的结构,压降最快,速度也增加最快。Velacty(mft)2Wt*«1PpHflfl(cm)图3.7元件表面速度分布4.小结本文不仅计算得到了燃料元件以正四棱锥方式堆积的球床内温度、速度和压力的分布规律,而且得到不同尺寸和不同元件堆积方式模型下的分布规律;讨论了孔隙率对

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论