一种受限空间三维流场的数值模拟方法、存储介质及系统与流程

本发明属于核设施放射性气溶胶输运基础研究，具体涉及一种受限空间三维流场的数值模拟方法、存储介质及系统。

背景技术：

1、在核设施的典型事故中，如泄漏事故、燃烧事故、爆炸事故等均可能涉及放射性气溶胶的输运。在核设施的设计、安全分析和事故应急中，放射性气溶胶在厂房内部输运至厂房外部之前，其中一个非常重要的环节是在受限空间内的输运，不同的事故情景、事故原因将导致受限空间内的建筑结构和流场状态参数随之发生变化，放射性气溶胶在受限空间内的输运过程也将受到严重影响，只有在准确得到流场状态参数之后才能准确计算放射性气溶胶的浓度分布。

2、目前，navier-stokes方程的求解思路有两种：一是求解一个含有速度和压力的动量方程和一个关于压力的泊松方程，即压力-速度方程；二是通过将方程两边取旋度，得到涡量输运方程，从动量方程中消去求解压力项带来的困难，即涡量-速度方程。涡量速度的形式的边界条件容易找到，并且解得的涡量能直接表示出流动的特性。

3、控制方程的数值求解方法有多种，其中以有限差分方法、有限体积方法和有限元方法应用最为广泛。有限差分方法思路简捷明了、编程容易和格式构造方便灵活、精度高；系数矩阵为稀疏矩阵，便于求解。有限体积方法对复杂物理模型的高度适应性；积分型守恒方程；控制方程无法做到维度解耦，构造高精度格式较为困难。有限元方法计算精度可提高，可求解复杂的物理模型；但其数值间断处存在虚假震荡，没有局部守恒，只能保证全局守恒。

技术实现思路

1、针对现有技术中存在的缺陷，本发明的目的是提供一种受限空间三维流场的数值模拟方法、存储介质及系统，采用直接数值模拟模型，得到受限空间三维流场的数值模拟程序算法。

2、为达到以上目的，本发明采用的技术方案是：一种受限空间三维流场的数值模拟方法，包括步骤：进行控制方程的推导及求解；进行参数无量纲化；建立数值模型，并确定初始条件及边界条件；利用数值模型进行方程组求解。

3、进一步，所述控制方程的推导及求解时，采用涡量-速度方法求解三维非定常不可压缩burgers模型方程，以涡量和速度作为求解变量，通过超松弛迭代方法反复迭代求解涡量输运方程和速度泊松方程，得到涡量和速度的数值解。

4、进一步，所述进行参数无量纲化时，引入特征长度、特征速度及特征时间，将所述控制方程中坐标及速度表示为无量纲形式。

5、进一步，所述边界条件包括涡量的边界条件及速度的边界条件，所述速度的边界条件采用第一类边界，所述涡量的边界条件须采用第二类边界。

6、进一步，求解任意一组方程时，采用升级后的sor进行求解，形式为：

7、

8、其中，表示方程左端未知项；、分别表示方程的数量和未知数的数量；为松弛因子；为方程右端已知项；为方程左端系数矩阵的元素，i，j表示系数矩阵a的行和列，imax＝jmax＝xmax×ymax×zmax，xmax、ymax、zmax分别表示笛卡尔坐标系下三维空间x、y、z方向网格数量的最大值。

9、进一步，所述升级后的sor求解时，未知数的编号至多为7个。

10、进一步，所述数值模型包括1个主程序及7个子程序，所述子程序包括3个涡量子程序、3个速度子程序及1个超松弛迭代子程序。

11、进一步，所述主程序用于参数定义、网格划分、参数初始化、迭代循环及结果输出。

12、本发明还提供一种存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现一种受限空间三维流场的数值模拟方法。

13、本发明还提供一种受限空间三维流场的数值模拟系统，包括：推导单元，用于进行控制方程的推导及求解；无量纲化单元，用于进行参数无量纲化；模型单元，用于建立数值模型，并确定初始条件及边界条件；求解单元，用于利用数值模型进行方程组求解。

14、本发明的效果在于：可针对核设施事故工况下放射性气溶胶在受限空间内的输运特点，采用直接数值模拟模型，以有限差分方法为基础，基于涡量-速度方法推导了三维非定常不可压缩burgers模型方程的求解方法，得到受限空间三维流场的数值模拟程序算法。

技术特征：

1.一种受限空间三维流场的数值模拟方法，其特征在于，包括：

2.如权利要求1所述的一种受限空间三维流场的数值模拟方法，其特征在于：

3.如权利要求1所述的一种受限空间三维流场的数值模拟方法，其特征在：

4.如权利要求1所述的一种受限空间三维流场的数值模拟方法，其特征在于：

5.如权利要求1所述的一种受限空间三维流场的数值模拟方法，其特征在于：

6.如权利要求4所述的一种受限空间三维流场的数值模拟方法，其特征在于：

7.如权利要求1所述的一种受限空间三维流场的数值模拟方法，其特征在于：

8.如权利要求7所述的一种受限空间三维流场的数值模拟方法，其特征在于：

9.一种存储介质，其上存储有计算机程序，其特征在于：所述计算机程序被处理器执行时实现权利要求1至8任一项所述的一种受限空间三维流场的数值模拟方法。

10.一种受限空间三维流场的数值模拟系统，其特征在于，包括：

技术总结
本发明涉及一种受限空间三维流场的数值模拟方法，包括步骤：进行控制方程的推导及求解；进行参数无量纲化；建立数值模型，并确定初始条件及边界条件；利用数值模型进行方程组求解。本发明还提供一种存储介质及受限空间三维流场的数值模拟系统，采用一种受限空间三维流场的数值模拟方法、存储介质及系统可采用直接数值模拟模型，以有限差分方法为基础，基于涡量‑速度方法推导了三维非定常不可压缩Burgers模型方程的求解方法，得到受限空间三维流场的数值模拟程序算法。

技术研发人员：朱业明,孟东原,杨亚鹏,梁博宁,孙洪超,冯宗洋,李国强,孙树堂,张建岗
受保护的技术使用者：中国辐射防护研究院
技术研发日：
技术公布日：2024/1/11