公开(公告)号：CN116644628A

公开(公告)日：2023-08-25

申请号：CN202310566449.3

申请日：2023-05-19

Applicant: 西安交通大学

Inventor： 张魁 ,  余颂娇 ,  刘毓民 ,  祁添翼 ,  田文喜 ,  苏光辉 ,  秋穗正

IPC: G06F30/23 ,  G06F30/28 ,  G06F111/10 ,  G06F113/08 ,  G06F119/08 ,  G06F119/14

Abstract: 本发明公开了一种判定弥散型板燃料元件起泡与沸腾临界先后顺序的数值模拟方法。包括以下步骤：1、填写并读入输入卡；2、设定时间循环次数与计算网格结构；3、计算燃料元件径向及轴向功率分布；4、进行稳态热工参数计算，更新燃耗和芯块最高温度；5、更新包壳和燃料板厚度；6、更新流道尺寸，计算燃料组件温度与热流密度；7、计算流道出口含汽率，根据沸腾临界类型计算临界热流密度；8、进行热流密度收敛判断，计算沸腾临界阈值温度；9、进行燃料颗粒开裂判断，计算裂变气体释放量；10、进行起泡判断，计算起泡阈值温度；11、将步骤8与步骤10温度进行对比，判定弥散型板燃料元件起泡与沸腾临界发生先后顺序。

公开(公告)号：CN116644628B

公开(公告)日：2023-11-07

申请号：CN202310566449.3

申请日：2023-05-19

Applicant: 西安交通大学

Inventor： 张魁 ,  余颂娇 ,  刘毓民 ,  祁添翼 ,  田文喜 ,  苏光辉 ,  秋穗正

IPC: G06F30/23 ,  G06F30/28 ,  G06F111/10 ,  G06F113/08 ,  G06F119/08 ,  G06F119/14

Abstract: 本发明公开了一种判定弥散型板燃料元件起泡与沸腾临界先后顺序的数值模拟方法。包括以下步骤：1、填写并读入输入卡；2、设定时间循环次数与计算网格结构；3、计算燃料元件径向及轴向功率分布；4、进行稳态热工参数计算，更新燃耗和芯块最高温度；5、更新包壳和燃料板厚度；6、更新流道尺寸，计算燃料组件温度与热流密度；7、计算流道出口含汽率，根据沸腾临界类型计算临界热流密度；8、进行热流密度收敛判断，计算沸腾临界阈值温度；9、进行燃料颗粒开裂判断，计算裂变气体释放量；10、进行起泡判断，计算起泡阈值温度；11、将步骤8与步骤10温度进行对比，判定弥散型板燃料元件起泡与沸腾临界发生先后顺序。

公开(公告)号：CN119227579A

公开(公告)日：2024-12-31

申请号：CN202411352088.3

申请日：2024-09-26

Applicant: 西安交通大学

Inventor： 张魁 ,  刘毓民 ,  余颂娇 ,  田文喜

IPC: G06F30/28 ,  G06F30/23 ,  G06F111/10 ,  G06F113/08 ,  G06F119/08 ,  G06F119/14 ,  G06F113/24

Abstract: 本发明公开了一种高燃耗起泡板型核燃料临界热流密度数值模拟方法，1、将高燃耗起泡板型核燃料几何模型及网格模型导入至商用CFD软件的ANSYS fluent中。2、设置计算模型并设置相应参数。3、根据需要计算的工况设置物性参数。4、设置模型边界条件。5、开启计算参数监测。6、在燃料芯块设置初始体积热源加热功率值进行加热，通过不断提高燃料芯块体积热源以提升流体域壁面热流密度，待加热壁面的最高温度和最高空泡份额读数收敛后，再提高燃料芯块功率密度。7、重复步骤6不断提高热源功率密度，当监测到加热壁面最高温度发生温度飞升或者加热壁面最大空泡份额大于0.8时，即认为此时发生临界，获取所得临界热流密度。本发明数值模拟方法对现有试验数据模拟结果良好，误差在‑25％到+10％。