CN109522624B - 建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法 - Google Patents

建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法 Download PDF

Info

Publication number
CN109522624B
CN109522624B CN201811291092.8A CN201811291092A CN109522624B CN 109522624 B CN109522624 B CN 109522624B CN 201811291092 A CN201811291092 A CN 201811291092A CN 109522624 B CN109522624 B CN 109522624B
Authority
CN
China
Prior art keywords
sealing plate
model
grid
dryer
oblique sealing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811291092.8A
Other languages
English (en)
Other versions
CN109522624A (zh
Inventor
王明军
房迪
田文喜
苏光辉
秋穗正
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Tongyuan Software & Control Technology Co ltd
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201811291092.8A priority Critical patent/CN109522624B/zh
Publication of CN109522624A publication Critical patent/CN109522624A/zh
Application granted granted Critical
Publication of CN109522624B publication Critical patent/CN109522624B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Drying Of Solid Materials (AREA)

Abstract

本发明公开了一种建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法,包括以下步骤:建立不包含密封板的干燥器组件几何结构;建立不包含斜密封板的干燥器组件的流体域模型;建立适于上述流体域模型的三维长方体拓扑块;在拓扑块上设置节点数目及分布,形成结构化六面体网格模型;检查每个网格单元中心的坐标是否满足斜密封板的投影直线方程,据此标记斜密封板的网格单元;将标记出的网格单元设置为固体,形成斜密封板的干燥器组件网格模型;本发明解决了对带斜密封板干燥器组件进行网格刻画时的几何模型建立困难、采用非结构网格时网格量较大、网格质量低等问题,同时为改变其倾斜角度而建立新的网格模型提供了便利。

Description

建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法
技术领域
本发明属于方法发明技术领域,具体涉及到一种建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法。
背景技术
蒸汽发生器是核电站系统重要大型换热设备,能够实现一回路与二回路之间的能量交换。在蒸汽发生器中,饱和蒸汽自传热管束区产生,但由于相间界面处液膜的破碎等原因导致饱和蒸汽会夹带一定量的液滴。这些液滴会随着饱和蒸汽向上流动经过各级汽水分离到达出口并进入二回路。
流体工质自传热管束区产生之后首先进入初级汽水分离器也称旋叶式汽水分离器区。单个旋叶式汽水分离器包括四个固定旋流叶片从而使液滴和蒸汽产生旋转流动,分离机制为利用液滴的离心效应的惯性分离机制。次级汽水分离器是带有一些钩子的波型叶片分离器,在其内部可以去除惯性较大的液滴。随后,蒸汽会离开蒸汽发生器进入二回路被送往汽轮机做功。目前欧美核电通用的标准要求蒸汽湿度小于0.25%,否则高速蒸汽流夹杂着这些液滴进入二回路会造成管道、关闭件的腐蚀,以及汽轮机叶片的气蚀,影响汽轮机寿命。因此,汽水分离器的分离效率就关系到核电站运行经济性和安全性。
国内外研发机构在开发波形板干燥器时,其技术路线通常为:首先采用数值模拟计算或冷态试验选型的方法确定特定型式的波形板片,然后通过冷态或热态试验确定波形板片的工作性能曲线。然而,受试验条件的限制,干燥器试验件的规模往往只能为实际蒸汽发生器干燥器的缩比模型,无法通过试验直接验证所开发干燥器在实际运行时的工作性能。此时,需要利用数值模拟方法对干燥器的三维流场进行分析。
随着计算机技术的发展,CFD程序针对蒸汽发生器干燥器三维流场的模拟能力逐渐提升。但是,在进行数值模拟时依然面临干燥器处几何结构复杂、网格划分困难和网格量大等问题。因此采用合理的方法对蒸汽发生器干燥器前的复杂结构进行简化处理将大大节省计算资源。本方法针对大型压水堆核电站内的蒸汽发生器干燥器前密封板处的三维流场模拟提出可行方法。方法适用于拥有密封板的蒸汽发生器干燥器单元的三维流场模拟。
发明内容
本发明的目的是设计一种建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法,该方法能够在保证斜密封板处网格质量的前提下,运用结构化六面体网格实现对该处的三维网格划分;能够在网格量较少的情况下,实现对立式自然循环蒸汽发生器内干燥器组件的网格模型的刻画。
本发明的解决方案是:
首先,根据干燥器组件设计图纸在solidworks中建立不带密封板的干燥器组件的几何模型,在该模型基础上,在solidworks软件中利用布尔运算功能得到不包含密封板的干燥器组件的流体域模型;然后使用结构化六面体网格刻画不含斜密封板的干燥器组件模型。
之后使用用户自定义函数检查上述六面体网格的每一个单元,若单元中心坐标符合斜密封板在投影面xoy面上的直线方程,则将该单元做标记;对所有斜密封板位置处的单元标记完成后,将这些单元设置为固体,即完成了对带斜密封板的干燥器组件的结构化六面体网格模型的建立。
本发明具体包括以下步骤:
步骤1:建立典型压水堆蒸汽发生器内干燥器组件的三维几何模型并得到其内部流体域:几何建模阶段不需要建立6个斜密封板的几何,具体地,根据真实设计参数,运用solidworks软件将组成干燥器组件实体的每个面建立出来,但不包含组成密封板的各面,并将这些面按照设计图纸中的相对位置关系组合成体,从而得到不带密封板的干燥器组件的三维几何模型;另外,在solidworks软件中重新建立一个与上述干燥器组件几何模型的直径、高度均相同的实心圆柱体;运用solidworks软件的布尔运算功能对实心圆柱体及干燥器组件模型进行布尔运算,得到不带密封板的干燥器组件的流体域模型;
步骤2:利用ANSYS-ICEM软件,使用结构化六面体网格刻画步骤1中建立的不带密封板的干燥器组件的流体域模型,具体分为以下步骤:
步骤2-1:在ANSYS-ICEM软件中创建一个三维长方体拓扑块包络上述流体域模型;将三维长方体拓扑块变为O形块以便于随后的操作中贴合干燥器组件流体域模型的圆柱形外形;
步骤2-2:依据流体域模型的几何特征,将O形块切分为多个拓扑块以贴合干燥器组件流体域模型几何结构,将以上从O形块中切分出的拓扑块的几何特征与干燥器组件流体域模型的几何特征进行关联,其中具体包括上述被切分出的拓扑块的点、线和面分别与对应的流体域模型中的点、线和面的关联;
步骤2-3:对在步骤2-2中所切分出的拓扑块各边进行节点参数设置,包括节点数目以及分布率,生成结构化六面体网格;
步骤3:在ANSYS-FLUENT软件中,使用用户自定义函数工具UDF,标记斜密封板的网格单元,将斜密封板网格标记出来并将其设置为固体,形成带斜密封板的干燥器组件的网格模型,具体分为以下步骤:
步骤3-1:根据斜密封板的倾斜角度,得到斜密封板的投影直线方程:
y=ax+b
斜密封板在与高度方向垂直的面假设为xoy面上的投影为一条倾斜直线,通过该投影可以得到斜密封板的直线方程,其中a为该直线的斜率,b为该直线的截距;
步骤3-2:检查步骤2中获得的结构化六面体网格的每个单元,若单元中心坐标符合步骤3-1中建立的直线方程,则对该单元进行标记,否则跳过该单元,从而标记出所有斜密封板的网格单元;
步骤3-3:在ANSYS-FLUENT软件中将已被标记出的斜密封板网格单元设置为固体,即得到带斜密封板的蒸汽发生器干燥器组件的网格模型。
本发明带来以下有益效果:
1)能够在数量较少的六面体结构化网格的基础上完成对次级干燥器进口前密封板的建模工作,且不影响网格质量;
2)容易实现斜密封板几何参数的变化例如倾斜角、遮挡长度等参数而不用重新进行三维建模;
3)对斜密封板区域采用标准六面体网格,可以减小网格标记的不同区域交界面处网格不连续带来的误差;
4)模型独立,方法通用性强,可以适应于不同类型的流体力学计算分析程序。
附图说明
图1为中国先进压水堆中的干燥器组件真实三维几何结构模型。
图2为不包含斜密封板的干燥器组件流体域几何结构示意图。
图3为不包含斜密封板的干燥器组件的网格模型。
图4为带斜密封板的单个干燥器单元在xoy平面上的投影。
图5为被标记出的用于刻画斜密封板的网格单元的示意图。
图6为带斜密封板的干燥器组件的网格模型示意图。
图7为本发明流程图。
具体实施方式
以下结合图7所示流程图,以先进压水堆为例,对本发明作进一步的详细描述。
本发明一种建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法,包括以下步骤:
步骤1:建立典型压水堆蒸汽发生器内干燥器组件的三维几何模型并得到内部流体域:几何建模阶段不需要建立6个斜密封板的几何,具体地,根据真实设计参数,运用solidworks软件将组成干燥器组件实体的每个面建立出来但不包含组成密封板的各面,并将这些面按照设计图纸中的相对位置关系组合成体,从而得到不带密封板的干燥器组件的三维几何模型。另外,重新建立一个与上述干燥器组件几何模型的直径、高度均相同的实心圆柱体。运用solidworks软件的布尔运算功能对实心圆柱体及干燥器组件的三维几何模型进行布尔运算,得到不带密封板的干燥器组件的流体域模型。包含密封板组件的真实三维几何结构以及不带斜密封板的组件流体域模型分别如图1和2所示。
步骤2:利用ANSYS-ICEM软件,使用结构化六面体网格刻画步骤1中建立的不带密封板的干燥器组件的流体域模型,具体分为以下步骤:
步骤2-1:在ANSYS-ICEM软件中创建一个三维长方体拓扑块包络上述流体域模型;将三维长方体拓扑变为O形块以便于随后的操作中贴合干燥器组件流体域模型的圆柱形外形;
步骤2-2:依据流体域模型的几何特征,将O形块切分为多个拓扑块以贴合干燥器组件流体域模型几何结构,将以上从O形块中切分出的拓扑块的几何特征与干燥器组件流体域模型的几何特征进行关联,其中具体包括上述被切分出的拓扑块的点、线和面分别与对应的流体域模型中的点、线和面的关联;
步骤2-3:对在步骤2-2中所切分出的拓扑块各边进行节点参数设置,包括节点数目以及分布率,生成结构化六面体网格,不包含斜密封板的干燥器组件的网格模型如图3所示。
步骤3:在ANSYS-FLUENT软件中,使用用户自定义函数工具UDF,标记斜密封板的网格单元,将斜密封板网格标记出来并将其设置为固体,形成带斜密封板的干燥器组件的网格模型,具体分为以下步骤:
步骤3-1:根据斜密封板在与高度方向垂直的投影面,xoy面上的投影,得到斜密封板的投影直线方程:
y=ax+b
如图4所示为干燥器组件中单个干燥器单元在xoy平面上的投影。其中,斜密封板在xoy平面的投影为一条倾斜直线,通过该投影可以得到斜密封板的直线方程,其中a为该直线的斜率,b为该直线的截距;
步骤3-2:检查步骤2中获得的结构化六面体网格的每个单元,若单元中心坐标符合步骤3-1中建立的直线方程,则对该单元进行标记,否则跳过该单元;从而标记出所有斜密封板的网格单元,被标记出的用于刻画斜密封板的网格单元如图5所示;
步骤3-3:在ANSYS-FLUENT软件中将已被标记出的斜密封板网格单元设置为固体,即得到带斜密封板的蒸汽发生器干燥器组件的网格模型,如图6所示。
本发明未详细说明部分属本领域技术人员公知常识。

Claims (1)

1.一种建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法,其特征在于,包括如下步骤:
步骤1:建立典型压水堆蒸汽发生器内干燥器组件的三维几何模型并得到其内部流体域:几何建模阶段不需要建立6个斜密封板的几何,具体地,根据真实设计参数,运用solidworks软件将组成干燥器组件实体的每个面建立出来,不包含组成密封板的各面,并将这些面按照设计图纸中的相对位置关系组合成体,从而得到不带密封板的干燥器组件的三维几何模型;另外,在solidworks软件中重新建立一个与上述干燥器组件几何模型的直径、高度均相同的实心圆柱体;运用solidworks软件的布尔运算功能对实心圆柱体及干燥器组件的三维几何模型进行布尔运算,得到不带密封板的干燥器组件的流体域模型;
步骤2:利用ANSYS-ICEM软件,使用结构化六面体网格刻画步骤1中建立的不带密封板的干燥器组件的流体域模型,具体分为以下步骤:
步骤2-1:在ANSYS-ICEM软件中创建一个三维长方体拓扑块包络上述流体域模型;将三维长方体拓扑块变为O形块以便于随后的操作中贴合干燥器组件流体域模型的圆柱形外形;
步骤2-2:依据流体域模型的几何特征,将O形块切分为多个拓扑块以贴合干燥器组件流体域模型几何结构,将以上从O形块切分出的拓扑块的几何特征与干燥器组件流体域模型的几何特征进行关联,其中具体包括上述被切分出的拓扑块的点、线和面分别与对应的流体域模型中的点、线和面的关联;
步骤2-3:对在步骤2-2中所切分出的拓扑块各边进行节点参数设置,包括节点数目以及分布率,生成结构化六面体网格;
步骤3:在ANSYS-FLUENT软件中,使用用户自定义函数工具UDF,标记斜密封板的网格单元,将斜密封板网格标记出来并将其设置为固体,形成带斜密封板的干燥器组件的网格模型,具体分为以下步骤:
步骤3-1:根据斜密封板的倾斜角度,得到斜密封板的投影直线方程:
y=ax+b
斜密封板在与高度方向垂直的面,假设为xoy面上的投影为一条倾斜直线,通过该投影得到斜密封板的直线方程,其中a为该直线的斜率,b为该直线的截距;
步骤3-2:检查步骤2中获得的结构化六面体网格的每个单元,若单元中心坐标符合步骤3-1中建立的直线方程,则对该单元进行标记,否则跳过该单元,从而标记出所有斜密封板的网格单元;
步骤3-3:在ANSYS-FLUENT软件中将已被标记出的斜密封板网格单元设置为固体,即得到带斜密封板的蒸汽发生器干燥器组件的网格模型。
CN201811291092.8A 2018-10-31 2018-10-31 建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法 Active CN109522624B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811291092.8A CN109522624B (zh) 2018-10-31 2018-10-31 建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811291092.8A CN109522624B (zh) 2018-10-31 2018-10-31 建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法

Publications (2)

Publication Number Publication Date
CN109522624A CN109522624A (zh) 2019-03-26
CN109522624B true CN109522624B (zh) 2019-07-02

Family

ID=65773279

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811291092.8A Active CN109522624B (zh) 2018-10-31 2018-10-31 建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法

Country Status (1)

Country Link
CN (1) CN109522624B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110009746B (zh) * 2019-04-08 2020-07-10 西安交通大学 反应堆燃料组件的带有边界层的六面体网格自动生成方法
CN111222211B (zh) * 2019-12-12 2023-03-24 合肥通用机械研究院有限公司 一种微通道单元cfd流固耦合结构化六面体网格划分方法
CN111122111A (zh) * 2019-12-23 2020-05-08 湖北保乾科技有限公司 一种波纹板式汽水分离器改进结构及其性能测试方法
CN112632735B (zh) * 2020-12-30 2022-10-28 西安交通大学 一种核反应堆蒸汽发生器网格划分方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091036A (zh) * 2014-07-31 2014-10-08 中国船舶重工集团公司第七一九研究所 一种自然循环蒸汽发生器的热传导建模与计算方法
CN105137957A (zh) * 2015-06-23 2015-12-09 黄红林 一种核电蒸汽发生器控制方法
CN106345186A (zh) * 2016-08-26 2017-01-25 上海核工程研究设计院 一种压水堆蒸汽发生器汽水分离器

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10290230B2 (en) * 2014-01-27 2019-05-14 Arizona Board Of Regents On Behalf Of Arizona State University Device specific finite element models for simulating endovascular treatment

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091036A (zh) * 2014-07-31 2014-10-08 中国船舶重工集团公司第七一九研究所 一种自然循环蒸汽发生器的热传导建模与计算方法
CN105137957A (zh) * 2015-06-23 2015-12-09 黄红林 一种核电蒸汽发生器控制方法
CN106345186A (zh) * 2016-08-26 2017-01-25 上海核工程研究设计院 一种压水堆蒸汽发生器汽水分离器

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
AP1000 蒸汽发生器管板应力和变形分析;刘晓薇 等;《科技视界》;20120630;57-59页
Thermodynamic performance analysis of OPR1000;J.A.Shamim;《The 2nd IEEE conference power engineering and renewable energy》;20141231;208-213页
蒸汽发生器干燥器声疲劳试验模型仿真研究;李奇 等;《核技术》;20130430;1-4页
蒸汽发生器干燥器气动噪声分析研究;张丰收 等;《核动力工程》;20170630;27-30页

Also Published As

Publication number Publication date
CN109522624A (zh) 2019-03-26

Similar Documents

Publication Publication Date Title
CN109522624B (zh) 建立带斜密封板的蒸汽发生器干燥器组件网格模型的方法
Savill One-point closures applied to transition
CN103488850B (zh) 一种热振复合残余应力定位均化的方法
CN110955991B (zh) 一种界面双向数据交换的流固耦合计算方法
CN104091036A (zh) 一种自然循环蒸汽发生器的热传导建模与计算方法
CN109472080B (zh) 一种建立特定脱汽份额条件下汽水分离器网格模型的方法
CN112632735B (zh) 一种核反应堆蒸汽发生器网格划分方法
CN110543704A (zh) 一种在反应堆堆芯流场计算中考虑局部结构影响的修正方法
CN113673186A (zh) 一种基于stl文件的笛卡尔网格快速生成方法
CN106126822A (zh) 一种直流锅炉的汽水分离器建模方法及系统
Cooke et al. Gas-liquid flow on smooth and textured inclined planes
Jin et al. A direct ALE multi-moment finite volume scheme for the compressible Euler equations
Castro Numerical difficulties in the calculation of complex turbulent flows
CN110728030A (zh) 基于表面式凝汽器数值模拟汽水两侧耦合的方法
CN106447781A (zh) 一种基于Minkowski和面向自动装配的碰撞检测方法
CN113609732B (zh) 一种蒸汽发生器管束支承板区域数值模拟计算方法
CN115169265A (zh) 基于数值分析的搅混系数分析方法、系统、设备和介质
Carlucci et al. The effects of symmetric/asymmetric boundary conditions on the flow of an internally heated fluid
CN110009746A (zh) 反应堆燃料组件的带有边界层的六面体网格自动生成方法
CN105247622A (zh) 用于模拟核反应堆中容器内流体流动的方法以及用于计算核反应堆堆芯组件的机械形变的方法,以及相关的计算机程序产品
Takasuo et al. Applications and development of the PORFLO 3D code in nuclear power plant thermal hydraulics
CN104573151A (zh) 一种壁面网格点数据获取的循环通信方法
CN117540507B (zh) 全环网格壁面距离的计算方法及装置
CN113158369B (zh) 静压推力轴承油垫封油边的油膜流量模拟监测方法
Batayneh et al. Modelling of Triggering and Steam Explosion Pressure Propagation with Validation Against KROTOS Experiments

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20220728

Address after: 215556 rooms e1701, e1702, e1703, e1801, e1803, e1804, No. 388 Ruoshui Road, China (Jiangsu) pilot Free Trade Zone, Suzhou, Jiangsu Province

Patentee after: SUZHOU TONGYUAN SOFTWARE & CONTROL TECHNOLOGY CO.,LTD.

Address before: Beilin District Xianning West Road 710049, Shaanxi city of Xi'an province No. 28

Patentee before: XI'AN JIAOTONG University

TR01 Transfer of patent right