FLUENT海洋条件下自然循环流动分析（二）

1年前浏览4218

正文共： 1673字 15图预计阅读时间： 5分钟

1 前言

随着船舶核动力技术的发展，人们开始系统地研究海洋条件对反应堆热工水力特性的影响。海洋浮动核动力装置不同于陆上固定核动力装置，洋面由于时刻处于运动中，核动力装置实质上也在不停的运动，这就使其热工水力特性有别于陆上。在海洋条件下，船舶的基本运动形式为：平动、升降（或称升潜）及摆（横摇和艏摇），实际的运动状态可能耦合了其中的几种。在这些运动下，流体除了受到重力加速度影响外，还可能受到平动加速度、角加速度、向心加速度以及科氏加速度（相关知识需要脑补运动学和动力学原理），这些额外的加速度引起了附加的惯性力从而影响流体的流动，特别是自然循环流体。其基本原理如下：

建立非惯性坐标系Oxyz（船舶）和惯性坐标系Ox’y’z’（地球），假设非惯性坐标系相对于惯性坐标系同时作平动和转动（通常将旋转中心设在非惯性坐标系的原点上），流体质点在非惯性坐标系Oxyz的径向量为r=xi+yj+zk。

在非惯性坐标系下的动量方程还应该考虑附加的加速度项a。

之前我们做了一个海洋条件下自然循环回路的流动特性分析案例，该案例尚未考虑温度影响，今天我们继续做一个完整的海洋条件下自然循环回路分析。

2 建模与网格

建立如下的二维平面自然回路模型（注：坐标原点位于右下角），流道直径100mm，回路外轮廓1000mm×1000mm，左下角和右上角分别有一段热源壁面和冷源壁面，长度均为100mm。划分四边形结构化网格，节点数19404，最小正交质量0.999。为了方便后续监测回路的流量变化，在回路的适当位置创建一个内部面。

3 边界与求解设置

设置介质为液态水，按如下模型设置水的密度，UDF代码如下

DEFINE_PROPERTY(rhowater,c,t){real rho0=998.2;real T0=293.15;real beta=0.0002;real temp=C_T(c,t);real rho;rho=rho0/(1+beta*(temp-T0));return rho;}

设置热源温度200℃，冷源温度30℃。

为了对比分析，我们先计算陆地条件下的自然循环回路，此时设置重力加速度为-Y方向，对于这种闭式自然循环回路计算，可先采用稳态方法计算，然后转为瞬态计算，通常情况下瞬态求解才能获得收敛解。可通过监测净功率来判断收敛情况，通过默认的残差标准可能无法达到收敛判定条件。本案例采用稳态计算也可以获得收敛解，只不过能量残差需要设置到1e-8。

对于海洋条件下的自然循环计算，可以用以上的陆地条件计算结果作为初始值，引入海洋条件运动参数（可点击文末“阅读原文”查看相应的内容）。

创建几个监视器，分别监测回路的质量流量，监测对象为建模时创建的内部面；监测热源和冷源的功率值，达到收敛时，这两个值之和应接近于零。

设置压力的离散格式为PRESTO!，对于封闭空间的自然循环问题，采用该格式可以更为准确地计算近壁面的速度矢量。

以零速度、零压力初始化流场（静止状态）。

4 计算结果

首先，看一下陆地条件下，残差曲线、回路的质量流量的迭代曲线以及热源和冷源的功率监测曲线，可以看出本案例稳态计算收敛性还是很好的，达到稳定状态时，回路的流量为1.2kg/s，功率约为9.2kW。

稳态收敛速度分布如下，最大值约为0.0338m/s。

海洋条件下，我们定义一个参数叫流量比率，也就是海洋条件瞬态的流量和陆地条件稳态流量的比值，回路的流量比率随时间的变化曲线如下。可以看出，对于本案例的自然循环回路配置，海洋条件引起的流动甚至显著大于自然循环流量。注意，本案例的回路尺寸均为假想，没有实际工程意义。

参考文献

[1]宫厚军,等. 海洋运动对自然循环流动影响的理论分析[J].

[2]钱立波,等. 一维冷却剂通道海洋条件附加力模型研究[J].

[3]Sipeng Wang, et al. CFD analysis on flow and heat transfer characteristic in natural circulation system with subchannels under rolling conditions[J].

[4]杜思佳,等. 海洋条件对单相强迫流动影响的理论研究[J].

来源：仿真与工程

UDF 船舶海洋理论

著作权归作者所有，欢迎分享，未经许可，不得转载

首次发布时间：2023-07-05

FLUENT海洋条件下自然循环流动分析（二）

西京学院王凯讲师：岩石热破裂分析的时空高精度数值流形法

大学本科课程 现代通信系统概论（上）

河南理工大学董卓讲师：岩石水力裂缝起裂与扩展特性的理论及数值模拟研究

ALIAS汽车设计职业培训班-汽车A面从入门到精通（提高班-持续加餐更新）

各种超单元法在产品开发仿真应用及实战技巧27讲（附学习模型）

第七届CMHL船舶与海洋工程计算水动力学国际研讨会（回放）

ANSYSWokrbench显示动力学案例6子弹冲击侵彻油箱流固耦合仿真及理论分析油箱破坏燃烧机理

多学科跨软件联合仿真28讲：Abaqus与Star-CCM 耦合海上浮动弹性板流固耦合仿真

大学本科课程现代通信系统概论（上）