基于声学有限元反应堆失水事故流固耦合模拟方法及系统技术方案

技术编号:39491468 阅读:12 留言:0更新日期:2023-11-24 11:13
本发明专利技术公开了基于声学有限元反应堆失水事故流固耦合模拟方法及系统,涉及反应堆流固耦合力学领域,其技术方案要点是:建立固体区域与流体区域的几何模型,并对几何模型进行网格划分,得到固体区域与流体区域的有限元模型;对失水事故工况下压力容器破口的压力变化进行求解,并获取压力边界条件;对有限元模型中固体区域进行二次建模;对有限元模型中流体区域进行二次建模;对完成固体区域和固体区域二次建模后的流固耦合模型进行流固耦合计算,得到流固耦合模拟分析结果

【技术实现步骤摘要】
基于声学有限元反应堆失水事故流固耦合模拟方法及系统


[0001]本专利技术涉及反应堆流固耦合力学领域,更具体地说,它涉及基于声学有限元反应堆失水事故流固耦合模拟方法及系统


技术介绍

[0002]反应堆是反应堆系统的核心,其性能对反应堆安全起至关重要的作用

冷却剂失水事故
(Loss of Coolant Accident

LOCA)
是一种管道破裂事故,发生失水事故后,冷却剂向外喷射或流失,产生喷射推力和引起压力容器内压力变化,卸压波穿过一回路冷却系统并向反应堆压力容器传播,使反应堆冷却剂系统受到很大的外载荷作用,引起堆内结构如吊篮产生振动,对系统安全运行带来不利影响,甚至导致部件和设备损坏的严重后果
。LOCA
动力响应分析的目的是通过对发生失水事故以后反应堆冷却系统

堆内构件可能发生的破坏情况的研究,获得堆内结构在卸压波传播过程中的水力载荷,计算出堆内结构的动态响应

动应力以及非线性接触力,从而评估反应堆的结构完整性和安全可靠性

该物理问题中涉及了流体动力学,结构动力学与非线性力学,使得求解过程复杂,计算资源消耗量大

工程设计需要一种对于该极端工况下反应堆动响应进行快速且准确的数值模拟方法

[0003]现有技术中公开了一种失水事故下含初始压力的压力容器流固耦合模拟方法,实现了对于极端工况下反应堆流固耦合动响应高精度仿真,其主要包括三个阶段:首先使用
Star

CD CFD
代码和
Abaqus/Standard FEM

code
针对
HDR v32
试验进行建模计算,随后使用
APROS
系统代码获得破口的压力边界条件作为两相计算,最后使用第三方
MPCCI
耦合软件处理
CFD

FEM
代码之间的耦合,实现流固耦合数值模拟

此外,现有技术还公开有采用
ADINA
软件反应堆失水事故
(LOCA)
流固耦合模拟分析方法,分别利用
ADINA
中结合
structures、CFD

ADINA

FSI
模块与
ADINA

Potential based formulation
模块,实现了基于
CFD

CSD
流固耦合与声固耦合的
LOCA
数值模拟,对比分析了两种方法的计算精度与效率,验证了声学方法在保证精度前提下的高效性

[0004]然而,采用失水事故下含初始压力的压力容器流固耦合模拟方法需要两个软件进行同时计算,并且需要两者之间数据传递,从而导致计算成本很大,计算效率极低;并且真实反应堆内存在大量非线性间隙,这对
CFD

CSD
流固耦合计算的收敛问题和计算资源提出巨大挑战

而采用
ADINA
软件反应堆失水事故
(LOCA)
流固耦合模拟分析方法虽然采用了一个
ADINA
软件,但
ADINA
软件工程领域应用较少,使用范围有所限制,因此上述方法不便于在工程单位中进行推广

此外,现有的失水事故数值模拟技术均未考虑实际工程中堆内结构件之间的非线性力学行为,而含非线性碰撞的流固耦合与声固耦合是工程中关注的重点以及数值模拟的难点

因此,如何研究设计一种能够克服上述缺陷的基于声学有限元反应堆失水事故流固耦合模拟方法及系统是我们目前急需解决的问题


技术实现思路

[0005]为解决现有技术中的不足,本专利技术的目的是提供基于声学有限元反应堆失水事故
流固耦合模拟方法及系统,采用对反应堆压力容器固体与流体区域分别进行模型的相对精确建立,并分别对固体与流体区域进行实时数据传输的仿真模拟,确保了对于压力容器在极端工况下
(
失水事故
)
堆内动响应的数值模拟具备较高的精度,并大大降低时间成本,实现对于
LOCA
的精确高效的分析

[0006]本专利技术的上述技术目的是通过以下技术方案得以实现的:
[0007]第一方面,提供了基于声学有限元反应堆失水事故流固耦合模拟方法,包括以下步骤:
[0008]依据反应堆的结构参数建立固体区域与流体区域的几何模型,并对几何模型进行网格划分,得到固体区域与流体区域的有限元模型;
[0009]基于有限元模型,对失水事故工况下压力容器破口的压力变化进行求解,并获取相应的水力载荷信息作为流固耦合模拟的压力边界条件;
[0010]对有限元模型中固体区域的固体边界条件

施加的非线性弹簧以及赋予固体的材料属性进行二次建模;
[0011]对有限元模型中流体区域的单元属性

流体介质参数

流体边界条件和流固耦合交界面进行二次建模;
[0012]对完成固体区域和固体区域二次建模后的流固耦合模型进行流固耦合计算,得到反应堆失水事故的流固耦合模拟分析结果

[0013]进一步的,所述结构参数包括几何尺寸

材料特性

接触刚度

支撑形式以及堆内构件的各阶振型与频率

[0014]进一步的,所述有限元模型中固体区域所对应的结构网格与流体区域所对应的流体网格的坐标位置相互匹配,且全为四面体单元

[0015]进一步的,所述有限元模型在构建过程中依据失水事故问题的对称性建立实际模型的一半,并通过对称处理得到最终的有限元模型

[0016]进一步的,所述有限元模型采用商用软件
Hypermesh
或商用软件
CATIA
进行建模

[0017]进一步的,所述水力载荷信息采用商用软件
Flowmaster
进行求解

[0018]进一步的,所述固体区域和
/
或流体区域采用商用软件
ANSYS
进行二次建模与耦合求解,二次建模时固体区域采用
SOLID185
单元,流体区域采用
FLUID30
单元

[0019]第二方面,提供了基于声学有限元反应堆失水事故流固耦合模拟系统,包括:
[0020]模型构建模块,用于依据反应堆的结构参数建立固体区域与流体区域的几何模型,并对几何模型进行网格划分,得到固体区域与流体区域的有限元模型;
[0021本文档来自技高网
...

【技术保护点】

【技术特征摘要】
1.
基于声学有限元反应堆失水事故流固耦合模拟方法,其特征是,包括以下步骤:依据反应堆的结构参数建立固体区域与流体区域的几何模型,并对几何模型进行网格划分,得到固体区域与流体区域的有限元模型;基于有限元模型,对失水事故工况下压力容器破口的压力变化进行求解,并获取相应的水力载荷信息作为流固耦合模拟的压力边界条件;对有限元模型中固体区域的固体边界条件

施加的非线性弹簧以及赋予固体的材料属性进行二次建模;对有限元模型中流体区域的单元属性

流体介质参数

流体边界条件和流固耦合交界面进行二次建模;对完成固体区域和固体区域二次建模后的流固耦合模型进行流固耦合计算,得到反应堆失水事故的流固耦合模拟分析结果
。2.
根据权利要求1所述的基于声学有限元反应堆失水事故流固耦合模拟方法,其特征是,所述结构参数包括几何尺寸

材料特性

接触刚度

支撑形式以及堆内构件的各阶振型与频率
。3.
根据权利要求1所述的基于声学有限元反应堆失水事故流固耦合模拟方法,其特征是,所述有限元模型中固体区域所对应的结构网格与流体区域所对应的流体网格的坐标位置相互匹配,且全为四面体单元
。4.
根据权利要求1所述的基于声学有限元反应堆失水事故流固耦合模拟方法,其特征是,所述有限元模型在构建过程中依据失水事故问题的对称性建立实际模型的一半,并通过对称处理得到最终的有限元模型
。5.
根据权利要求1所述的基于声学有限元反应堆失水事故流固耦合模拟方法,其特征是,所述有限元模型采用商用软件
Hypermesh
或商用软件
CATIA
进行建模
。6.
根据权利...

【专利技术属性】
技术研发人员:叶献辉刘文进张毅雄曾忠秀王新军卢喜丰党高健张渝冯志鹏曹阳吕勇波李丽娟
申请(专利权)人:中国核动力研究设计院
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1