利用移动粒子有限容积法研究高瑞利数熔融池换热特性的方法技术

一种利用移动粒子有限容积法研究高瑞利数熔融池换热特性的方法，主要步骤如下：1、针对具体问题进行粒子建模，设定粒子初始布置；2、根据各粒子位置划分网格；3、在时间步长Δt之后，使用有限容积法求解能量方程；4、粒子法处理显式计算动量方程中的粘性项和源项，估算粒子速度及位置；5、计算粒子数密度，求解压力Poisson方程；6、计算速度的修正值，并计算各粒子在下一时层的速度及位置；7、输出计算结果，并根据求解的具体问题输出壁面热流密度和熔融池内物质的相态。本方法具有粒子法精确捕捉自由液面及模拟相变问题的优点，又具有较高的计算精度，能够通过此方法来研究高瑞利数熔融池换热特性，为核电厂反应堆严重事故安全特性研究提供重要依据。

Study on Heat Transfer Characteristics of High Rayleigh Number Melting Pool by Moving Particle Finite Volume Method

【技术实现步骤摘要】
利用移动粒子有限容积法研究高瑞利数熔融池换热特性的方法
本专利技术涉及核电厂严重事故下封头熔融池换热特性研究
，具体涉及一种利用移动粒子有限容积法研究高瑞利数熔融池换热特性的方法。
技术介绍
当核电厂压水堆发生严重事故时，堆芯如果不能得到有效冷却，会快速升温导致大量燃料组件和结构材料发生熔化并迁移至压力容器底部，进入到下封头内的高温熔融物逐渐积累形成碎片床或熔融池，不断向下封头壁面传递衰变热，当衰变热不能充分导出时，熔融池会将下封头壁面加热到很高温度，严重威胁压力容器的完整性。反应堆压力容器下封头内熔融池的流动还热特性对下封头壁面热负荷分布以及外部冷却能力有着重要的影响，熔融物堆内持留技术是缓解堆芯熔化事故后果的一项关键措施，运用不同的管理策略，可将熔融物滞留在压力容器下封头内，壁面后续事故的发生和大量放射性物质释放到外界环境中，从而可以有效终止反应堆严重事故。关于熔融池换热特性的研究，国内外已经开展了一些研究，包括实验研究和数值模拟研究。对于实验研究，由于不同的研究侧重点和实验条件的限制，各实验采用不同比例的实验装置、不同的熔融物模拟物、不同的加热方式以及不同的边界条件，因此也具有不同的瑞利数范围和实验结果。瑞利数作为对流换热问题的关键参数，对熔融池内的对流换热特性具有重要的影响。对比不同实验得到的熔融池换热关系式可以发现，在瑞利数相对较小的范围内不同实验获得的熔融池换热特性参数差别较小。但是在高瑞利数范围内，实验结果的差别开始增大，特别是在反应堆量级下(1016～1017)，不确定性更大。所以国内外也开展了一些大型实验，但相应的实验花费也更大。除了实验研究，国内外还开展了许多理论研究及数值模拟研究，以求在低消耗低成本的基础上对高瑞利数熔融池换热特性进行比较深入的研究。目前使用的方法大多限于商业软件如Fluent等先对下封头熔融池进行建模，然后进行网格划分，其中使用较多的方法是有限体积法，因其具有很好的守恒性，对网格的适应性很好，可以很好的解决复杂的工程问题。但对于反应堆下封头熔融池换热特性而言，其内可能存在氧化层和金属层，还有因为冷却变化和再次受衰变热加热而导致的硬壳生成及再熔化的相变问题，网格法在处理这些问题上不及粒子法的优势大，而且粒子法建模简单，而且能够处理复杂形状的流体流动问题，在处理湍流问题上具有更好的鲁棒性。尽管粒子法在处理流体流动及存在相变的问题上的上有其独特的优势，但是粒子法在计算传热问题时，在计算域的边界因为支持域被截断而存在较大的误差。综上所述，网格法中的有限容积法在计算传热问题上精确度较高，粒子法在求解复杂流体流动问题及相变问题上优势较大。本专利技术综合粒子法和有限容积法提出了一种研究高瑞利数熔融池传热特性的移动粒子有限容积法。
技术实现思路
为了克服网格法与粒子法的缺陷，本专利技术的目的是为了精确并高效地对下封头熔融池换热特性进行研究，提供一种利用移动粒子有限容积法研究高瑞利数熔融池换热特性的方法，该方法能够对存在多种流体、存在相变问题、湍流流动的熔融池换热进行研究，获得熔融池内的温度场分布和壁面热流密度分布情况、熔融池内的硬壳分布特性，为核电厂反应堆严重事故安全特性研究提供重要依据。为了实现上述目的，本专利技术采取了以下的技术方案予以实施：一种利用移动粒子有限容积法研究高瑞利数熔融池换热特性的方法，步骤如下：步骤1：对高瑞利数熔融池进行粒子建模，不同种类的粒子代表不同的物质，每种粒子根据所表示的物质具有其对应的质量、密度和熔化焓，使用粒子i来表示某一粒子，则其质量、密度、熔化焓对应标记分别为mi、ρi、设定粒子的初始布置，包括粒子的位置、速度、压力和温度，对应标记分别为ri、ui、pi和Ti，粒子初始布置时相邻两粒子间的距离记为l0；步骤2：根据各粒子的位置划分网格：在本步骤中，设置一个粒子搜索半径，对于每一个粒子i，其坐标位置为(xi，yi，zi)，与其间距lij≤1.2l0的所有粒子作为搜索的对象，式中lij为j粒子与i粒子的距离，j粒子的坐标标记为(xj，yj，zj)，lij使用公式(1)进行计算，并将搜索到的粒子使用j1，j2，……，jn符号进行标记，其中下标n表示搜索到的粒子总数；将j1，j2，……，jn粒子中的每三个粒子进行组合，将有种组合，由公式(2)进行计算，对每个组合中的粒子进行考察，在这些组合中，每个组合内的三个坐标点必须满足构成面的要求，取其中的一个组合进行解析，其中一种组合粒子标记为j1，j2，j3，粒子中心分别标记为O1，O2，O3，空间坐标对应为(x1，y1，z1)，(x2，y2，z2)，(x3，y3，z3)，连接O1、O2为直线O1O2，连接O1、O3为直线O1O3，连接O2、O3为直线O2O3，若要判断O1、O2、O3能否构成一个平面，只需要判断直线O1O2与直线O1O3之间存在一个不为0的夹角，则只需满足公式(3)，对所有的组合进行判断，将能够构成面的组合保留，不能构成面的组合删去，随后考虑i粒子与每个组合中三个粒子的关系，i粒子与组合中三个粒子必须能够构成四面体，因为组合中的三个粒子能够构成平面，所以要使这四个粒子能够构成四面体，只需满足i粒子不在组合中三个粒子所构成的平面上即可，添加i粒子，粒子中心标记为O，空间坐标定义为(xi，yi，zi)，若要使O、O1、O2、O3四个点构成四面体，因为O1、O2、O3三个点已经满足构成面的要求，则只需要满足i粒子到平面O1O2O3的距离不为0，使用公式(4)进行判定，式中：——i粒子到平面O1O2O3的距离m；a、b、c、d——平面O1O2O3的平面方程的四个系数，通过公式(5)获得，i粒子与j1、j2、j3粒子构成的平面O1O2O3的距离若满足公式(4)，则说明i粒子能够与三角形O1O2O3构成四面体OO1O2O3，四面体OO1O2O3的四个面分别为三角形OO1O2、三角形OO1O3、三角形OO2O3和三角形O1O2O3；使用海伦公式计算四个三角形的面积，四个面的面积对应标记为海伦公式如公式(6)所示，式中，S——三角形的面积；L1、L2、L3——三角形的三边长；p——三角形的半周长，即p＝(L1+L2+L3)/2；四面体OO1O2O3的体积通过公式(7)进行计算，式中，V——四面体体积，在四面体OO1O2O3中，标记为S——三角形O1O2O3的面积，标记为使用公式(6)计算；h——O点与三角形所在平面O1O2O3的距离，记为使用公式(4)计算；在使用粒子位置划分网格时，还有几个参量需要求解，其一是四面体OO1O2O3中四个顶点O、O1、O2、O3分别对四面体总体积的分属控制体积，对应标记为VO、这个分属控制体积由公式(8)进行计算，实际上，四面体OO1O2O3的四个点的分属控制体积与其相对的三角形面积成反比，即有公式(9)所示的关系，其二是对四面体OO1O2O3中四个面的三角形进行面积划分，在三角形O1O2O3中，三角形O1O2O3的面积通过公式(6)求出；点M12为线段O1O2的中点，点M13为线段O1O3的中点，点M23为线段O2O3的中点。此时需要在三角形中寻找一点P，连接P、O1，连接P、O2，连接P、O3，连接P、M12，连接P、M13，连接P、M23，此时三角形O1O2O3被分成本文档来自技高网...

【技术保护点】
1.一种利用移动粒子有限容积法研究高瑞利数熔融池换热特性的方法，其特征在于：步骤如下：步骤1：对高瑞利数熔融池进行粒子建模，不同种类的粒子代表不同的物质，每种粒子根据所表示的物质具有其对应的质量、密度和熔化焓，使用粒子i来表示某一粒子，则其质量、密度、熔化焓对应标记分别为mi、ρi、

【技术特征摘要】
1.一种利用移动粒子有限容积法研究高瑞利数熔融池换热特性的方法，其特征在于：步骤如下：步骤1：对高瑞利数熔融池进行粒子建模，不同种类的粒子代表不同的物质，每种粒子根据所表示的物质具有其对应的质量、密度和熔化焓，使用粒子i来表示某一粒子，则其质量、密度、熔化焓对应标记分别为mi、ρi、设定粒子的初始布置，包括粒子的位置、速度、压力和温度，对应标记分别为ri、ui、pi和Ti，粒子初始布置时相邻两粒子间的距离记为l0；步骤2：根据各粒子的位置划分网格：在本步骤中，设置一个粒子搜索半径，对于每一个粒子i，其坐标位置为(xi，yi，zi)，与其间距lij≤1.2l0的所有粒子作为搜索的对象，式中lij为j粒子与i粒子的距离，j粒子的坐标标记为(xj，yj，zj)，lij使用公式(1)进行计算，并将搜索到的粒子使用j1，j2，……，jn符号进行标记，其中下标n表示搜索到的粒子总数；将j1，j2，……，jn粒子中的每三个粒子进行组合，将有种组合，由公式(2)进行计算，对每个组合中的粒子进行考察，在这些组合中，每个组合内的三个坐标点必须满足构成面的要求，取其中的一个组合进行解析，其中一种组合粒子标记为j1，j2，j3，粒子中心分别标记为O1，O2，O3，空间坐标对应为(x1，y1，z1)，(x2，y2，z2)，(x3，y3，z3)，连接O1、O2为直线O1O2，连接O1、O3为直线O1O3，连接O2、O3为直线O2O3，若要判断O1、O2、O3能否构成一个平面，只需要判断直线O1O2与直线O1O3之间存在一个不为0的夹角，则只需满足公式(3)，对所有的组合进行判断，将能够构成面的组合保留，不能构成面的组合删去，随后考虑i粒子与每个组合中三个粒子的关系，i粒子与组合中三个粒子必须能够构成四面体，因为组合中的三个粒子能够构成平面，所以要使这四个粒子能够构成四面体，只需满足i粒子不在组合中三个粒子所构成的平面上即可，添加i粒子，粒子中心标记为O，空间坐标定义为(xi，yi，zi)，若要使O、O1、O2、O3四个点构成四面体，因为O1、O2、O3三个点已经满足构成面的要求，则只需要满足i粒子到平面O1O2O3的距离不为0，使用公式(4)进行判定，式中：——i粒子到平面O1O2O3的距离m；a、b、c、d——平面O1O2O3的平面方程的四个系数，通过公式(5)获得，i粒子与j1、j2、j3粒子构成的平面O1O2O3的距离若满足公式(4)，则说明i粒子能够与三角形O1O2O3构成四面体OO1O2O3，四面体OO1O2O3的四个面分别为三角形OO1O2、三角形OO1O3、三角形OO2O3和三角形O1O2O3；使用海伦公式计算四个三角形的面积，四个面的面积对应标记为海伦公式如公式(6)所示，式中，S——三角形的面积；L1、L2、L3——三角形的三边长；p——三角形的半周长，即p＝(L1+L2+L3)/2；四面体OO1O2O3的体积通过公式(7)进行计算，式中，V——四面体体积，在四面体OO1O2O3中，标记为S——三角形O1O2O3的面积，标记为使用公式(6)计算；h——O点与三角形所在平面O1O2O3的距离，记为使用公式(4)计算；在使用粒子位置划分网格时，还有几个参量需要求解，其一是四面体OO1O2O3中四个顶点O、O1、O2、O3分别对四面体总体积的分属控制体积，对应标记为VO、这个分属控制体积由公式(8)进行计算，实际上，四面体OO1O2O3的四个点的分属控制体积与其相对的三角形面积成反比，即有公式(9)所示的关系，其二是对四面体OO1O2O3中四个面的三角形进行面积划分，在三角形O1O2O3中，三角形O1O2O3的面积通过公式(6)求出；点M12为线段O1O2的中点，点M13为线段O1O3的中点，点M23为线段O2O3的中点。此时需要在三角形中寻找一点P，连接P、O1，连接P、O2，连接P、O3，连接P、M12，连接P、M13，连接P、M23，此时三角形O1O2O3被分成了6个小三角形，分别标记为三角形O1PM12，三角形O1PM13，三角形O2PM12，三角形O2PM23，三角形O3PM13，三角形O3PM23；因为点M13为线段O1O3的中点，所以三角形O1PM13与三角形O3PM13面积相等，标记为同样的，三角形O1PM12与三角形O2PM12面积相等，标记为三角形O2PM23与三角形O3PM23面积相等，标记为对于各个三角形的面积，有如公式(10)和公式(11)所示的几何关系，三角形的三个角的角度由公式(3)所示的余弦定理求出；由公式(10)和公式(11)求出三角形O1O2O3所分成的6个小三角形的面积；同理也能求得四面体OO1O2O3另外三个面所分成的共计18个小三角形的面积；在四面体OO1O2O3中寻找一点Q，对于顶点O的控制体积而言，连接QP1，QP2，QP3，QM1，QM2，QM3，QO，至此点O...

【专利技术属性】
技术研发人员：田文喜，李勇霖，陈荣华，苏光辉，秋穗正，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：陕西,61