CN115470723A - 一种获取水射流刺破熔池液面临界条件的方法 - Google Patents

一种获取水射流刺破熔池液面临界条件的方法 Download PDF

Info

Publication number
CN115470723A
CN115470723A CN202211062257.0A CN202211062257A CN115470723A CN 115470723 A CN115470723 A CN 115470723A CN 202211062257 A CN202211062257 A CN 202211062257A CN 115470723 A CN115470723 A CN 115470723A
Authority
CN
China
Prior art keywords
liquid level
liquid
jet
impact
model
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.)
Granted
Application number
CN202211062257.0A
Other languages
English (en)
Other versions
CN115470723B (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.)
Nuclear Power Institute of China
Original Assignee
Nuclear Power Institute of China
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 Nuclear Power Institute of China filed Critical Nuclear Power Institute of China
Priority to CN202211062257.0A priority Critical patent/CN115470723B/zh
Publication of CN115470723A publication Critical patent/CN115470723A/zh
Application granted granted Critical
Publication of CN115470723B publication Critical patent/CN115470723B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Fluid Mechanics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Carbon Steel Or Casting Steel Manufacturing (AREA)
  • Continuous Casting (AREA)

Abstract

本发明公开了一种获取水射流刺破熔池液面临界条件的方法,分析水射流的射流和液滴冲击自由钢液面的稳态受力情况和液面变形情况,获取冲击力总和模型和支撑力总和模型;以变形液面轮廓曲线在三相接触线处的倾斜角
Figure DDA0003826736900000011
为临界条件,构建射流冲击刺破模型和液滴冲击刺破模型以计算射流和液滴刺破熔池液面的初始速度和下落高度;针对水射流和液滴从顶部冲击熔池自由液面的物理过程,建立数理模型,确立冷却水刺破不锈钢熔池液面的临界条件,为使用和评价堆内注水来缓解金属层热聚焦效应的方案提供理论依据,为计算出各种不同工况下的临界速度奠定基础,为1400MW级超大型堆芯熔融物的滞留提供了理论基础和分析工具。

Description

一种获取水射流刺破熔池液面临界条件的方法
技术领域
本发明涉及反应堆热工水力技术领域,具体涉及一种获取水射流刺破熔池液面临界条件的方法。
背景技术
轻水反应堆在严重事故条件下在下封头内形成两层或三层的熔池结构,中间为主要发热源的氧化物层,顶部为轻质的不锈钢和少量锆的金属层,整个熔池结构使得下封头外壁面在两个关键角度位置的热流密度较高:靠近氧化物熔融池顶部的位置和面向金属层的位置。前者是因为氧化物熔融池自然对流导致的,或者是因为金属层热聚焦效应导致的,两个位置相邻,这使得反应堆下封头容易在中上部角度位置失效。
业内为了防止熔融物堆内滞留的失效,采用了两个方面的手段在应对这种致命缺陷。①压力容器下封头外壁面建立保温层内流道,建立下封头外表面水冷系统,通过自然循环换热或强迫循环增强外壁面的排热能力。②增加顶部轻金属层的厚度,在设计环节增加堆芯围板、支撑孔板等的不锈钢质量,使得熔化形成的轻金属层厚度增大,增加金属层与压力容器接触处的传热面积,减缓甚至消除热聚焦效应。
然而在面对1400MW以上的大功率堆芯时,提升外部冷却能力和增加熔融态不锈钢量的方法已到极限,业内提出向下封头熔池顶部注水的方式在缓解热聚焦效应,但这将会使得熔融钢相与冷却水直接接触,存在触发蒸汽爆炸的风险;传统方法无法准确获取射流和液滴触发蒸汽爆炸的临界条件;也没有针对大功率堆芯射流和液滴触发蒸汽爆炸过程的具体分析模型以便获取蒸汽爆炸的临界阈值。
发明内容
本发明所要解决的技术问题是:对于1400MW以上的大功率堆芯,传统的方法无法准确获取射流和液滴触发蒸汽爆炸的临界条件;也没有针对大功率堆芯射流和液滴触发蒸汽爆炸过程的具体分析模型以便获取蒸汽爆炸的临界阈值;鉴于此本发明提供一种获取水射流刺破熔池液面临界条件的方法,通过水射流冲击自由钢液面的稳态受力情况和液面变形情况进行微观分析构建射流冲击刺破模型和液滴冲击刺破模型,针对水射流和液滴从顶部冲击熔池自由液面的物理过程,建立数理模型,确立冷却水刺破不锈钢熔池液面的临界条件,解决了上述技术问题。
本发明通过下述技术方案实现:
本方案提供一种获取水射流刺破熔池液面临界条件的方法,包括步骤:
步骤一:分析水射流的射流冲击自由钢液面的稳态受力情况和液面变形情况,获取冲击力总和模型和支撑力总和模型;
步骤二:分析水射流的液滴冲击自由钢液面的稳态受力情况和液面变形情况,获取支撑力模型;
步骤三:以变形液面轮廓曲线在三相接触线处的倾斜角
Figure BDA0003826736880000021
为临界条件,构建射流冲击刺破模型和液滴冲击刺破模型;
步骤四:根据射流冲击刺破模型计算射流刺破熔池液面的初始速度和下落高度;根据液滴冲击刺破模型计算液滴刺破熔池液面的初始速度和下落高度。
本方案工作原理:对于1400MW以上的大功率堆芯,传统的方法无法准确获取射流和液滴触发蒸汽爆炸的临界条件;也没有针对大功率堆芯射流和液滴触发蒸汽爆炸过程的分析模型以便获取蒸汽爆炸的临界阈值;鉴于此本发明提供一种获取水射流刺破熔池液面临界条件的方法,通过水射流冲击自由钢液面的稳态受力情况和液面变形情况进行微观分析构建射流冲击刺破模型和液滴冲击刺破模型,针对水射流和液滴从顶部冲击熔池自由液面的物理过程,建立数理模型,确立冷却水刺破不锈钢熔池液面的临界条件,解决了上述技术问题。
目前APR1400标准版反应堆就采用顶部注水方案来应对金属层热聚焦效应;国外CEA对此过程也进行了实验研究,得出结果显示:射流和液滴均无法触发蒸汽爆炸,顶部注水无蒸汽爆炸风险。而发明人开展了实验研究,发现10m/s的Φ10mm射流和Φ15mm射流在0.5m的喷嘴高度均触发了剧烈的蒸汽爆炸,爆炸冲击波达6MPa。目前尚无针对顶部注水冲击金属熔池自由液面的过程的理论模型,来给出触发蒸汽爆炸的临界阈值。本方案针对射流、液团、液滴冲击金属熔池自由液面的物理过程,建立数理模型,确立冷却水刺破不锈钢熔池液面的临界条件。反应堆严重事故下的下封头熔融池冷却,单独的外部冷却不足以满足1400MW以上的大功率堆芯的熔池冷却需求,特别是薄金属层诱发的热流集中效应使得下封头在中上部失效的概率升高,降低了IVR策略的安全裕度,向熔池顶部的轻金属层注水成为了IVR技术继续推行的重要解决方案,但同时存在触发蒸汽爆炸并加速下封头破裂的风险。建立水射流冲击熔池自由液面的数理模型,分析触发蒸汽爆炸的临界条件,结合实验研究数据验证模型的精确性,指导堆内注水策略。
使用和评价堆内注水来缓解金属层热聚焦效应的方案终于有了一款理论工具,可以计算出各种不同工况下的临界速度,堆内注水的缓解措施中流速的设计必须不超过此阈值。此外,此模型亦可以用于严重事故分析软件中,评价各种事故序列下顶部注水后的堆内滞留效果,助力新型严重事故缓解措施的研发。同时,对于1400MW级超大型堆芯熔融物的滞留提供了理论基础和分析工具。
进一步优化方案为,步骤一包括:
S11、分析水射流的射流冲击自由钢液面的冲击压力,对冲击区内的径向冲击力积分构建出冲击力总和模型;
S12、分析自由钢液面对射流的支撑力和液面变形情况构建出支撑力总和模型。
进一步优化方案为,所述冲击力总和模型为:
Figure BDA0003826736880000031
其中,vj为射流中心的最大流速;H为射流下落高度;ρl为射流的密度。
进一步优化方案为,所述支撑力总和模型为:
Figure BDA0003826736880000032
其中,Fσ为自由钢液面熔池表面张力,Fg为水射流没入钢液面以下排开的钢液重力;
Figure BDA0003826736880000033
为变形液面轮廓曲线三相接触线处的倾斜角,所述三相为钢液、空气和冷却水三相;ρ为钢液的密度;rm为冲击区最大直径;g为重力加速度;σ为钢液面的表面张力;hm为射流冲击区的高度。
进一步优化方案为,所述射流冲击刺破模型获取方法为:
在射流达稳态条件下,射流冲击自由钢液面的冲击力总和等于自由钢液面对射流的支撑力总和;
Figure BDA0003826736880000034
的临界条件下,自由钢液面对射流的支撑力总和达到最大值得到射流冲击刺破模型。
进一步优化方案为,所述射流冲击刺破模型为:
Figure BDA0003826736880000035
进一步优化方案为,步骤二包括以下子步骤:
T21、将水射流的液滴等效为等体积、等质量的球形固态物体,分析球形固态物体冲击自由钢液面的受力情况和液面变形情况;
T22、通过自由钢液面的表面张力沿竖直方向分力与水射流的液滴没入自由钢液面以下排开的钢液重力之和,推导出自由钢液面对液滴的支撑力模型。
进一步优化方案为,所述自由钢液面对液滴的支撑力模型为:
Figure BDA0003826736880000036
β为浸没角,rt为球形固态物体的球半径;ρ为钢液的密度;g为重力加速度;hlr为三相接触线与熔池自由液面之间的距离;
Figure BDA0003826736880000041
为变形液面轮廓曲线在三相接触线处的倾斜角。
进一步优化方案为,所述液滴冲击刺破模型为:
Figure BDA0003826736880000042
其中,hcr为球形固态物体刺破自由钢液面的下陷深度;θc为自由钢液面在球形固态物体处的静态接触角;z为液滴的下落高度;ρ表示液滴的密度;v0为液滴从下落处最高点的初始速度;
Figure BDA0003826736880000043
式中A=2πσrt
Figure BDA0003826736880000044
v0为球形固态物体从下落高处最高点的初始速度。
进一步优化方案为,根据液滴冲击刺破模型,在给定的rt,z条件下,得到液滴刺破熔池液面的临界初速vcr
Figure BDA0003826736880000045
的临界条件下,根据式:
Figure BDA0003826736880000046
计算球形固态物体刺破熔池自由液面的临界下陷深度hcr
本发明与现有技术相比,具有如下的优点和有益效果:
本发明提供一种获取水射流刺破熔池液面临界条件的方法,通过水射流冲击自由钢液面的稳态受力情况和液面变形情况进行微观分析构建射流冲击刺破模型和液滴冲击刺破模型,针对水射流和液滴从顶部冲击熔池自由液面的物理过程,建立数理模型,确立冷却水刺破不锈钢熔池液面的临界条件,为使用和评价堆内注水来缓解金属层热聚焦效应的方案提供理论依据,为计算出各种不同工况下的临界速度奠定基础,为1400MW级超大型堆芯熔融物的滞留提供了理论基础和分析工具。
附图说明
为了更清楚地说明本发明示例性实施方式的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。在附图中:
图1为获取水射流刺破熔池液面临界条件的方法流程示意图;
图2为冲击射流结构示意图;
图3为射流冲击自由钢液面受力示意图;
图4为液滴冲击自由钢液面受力示意图;
图5为不同喷嘴高度下的临界射流最大速度;
图6为不同液滴直径下的临界下落初始速度。
在附图中:
A-喷嘴,B-核心射流区,C-自由射流区,D-滞止区,E-壁面射流区,F-壁面,G-射流边界,H-壁面驻点,I-冲击区。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
本实施例提供一种获取水射流刺破熔池液面临界条件的方法,如图1所示,包括步骤:
步骤一:分析水射流的射流冲击自由钢液面的稳态受力情况和液面变形情况,获取冲击力总和模型和支撑力总和模型;
具体包括子步骤:
S11、分析水射流的射流冲击自由钢液面的冲击压力,对冲击区内的径向冲击力积分构建出冲击力总和模型;
水射流冲击熔池自由钢液面,会在熔池液面形成冲击区,自由钢液面凹陷,给射流冲击区给予支撑,当射流冲击从初始接触到达到稳态之后,射流冲击区对钢液面的冲击力达到最大值。分析水射流冲击钢液面的稳态受力情况,以此建立水射流冲击金属熔池自由液面的数理模型。
圆形紊动垂直射流的结构如图2所示,射流通过喷嘴A喷出,在射流边界G以内构成核心射流区B和自由射流区C,射流与壁面F之间形成滞止区D、壁面射流区E;核心射流区B的中心撞击点为壁面驻点H。
射流对于壁面F的冲击力只分布于冲击区I内,且呈正态分布,则在冲击区I内沿径向的冲击力表达式为:
Figure BDA0003826736880000061
式中pm为冲击区滞止点处压力,也即最大时均压力,bp为冲击压力等于
Figure BDA0003826736880000062
处的r值。冲击区范围约为:2.2≤r/bp≤2.2,冲击射流理论也给出,冲击区I的范围主要跟射流下落高度相关联,冲击区具有高度x/H≈0.86,径向r/H≈0.22的尺寸,其中H为射流下落高度。
由冲击区范围可知,
Figure BDA0003826736880000063
冲击区最大直径rm=0.22H,则可知射流的冲击力总和为:
Figure BDA0003826736880000064
式中
Figure BDA0003826736880000065
vj为射流中心的最大流速,出现于射流进入冲击区前;积分之后得到冲击力总和模型:
Figure BDA0003826736880000066
其中,vj为射流中心的最大流速;H为射流下落高度;ρl为射流的密度。由冲击力总和模型可知总的冲击力只跟射流的密度、最大流速、下落高度呈正相关。
S12、分析自由钢液面对射流的支撑力和液面变形情况构建出支撑力总和模型。
射流形成的冲击区压迫熔池液面凹陷导致整个自由钢液面弯曲变形,熔池液面张力产生一个整体向上的支撑力托住冲击区的液相。
由广义阿基米德浮力定律可知,这个支撑力由表面张力沿竖直向上的分力,再加上水射流没入钢液面以下排开的钢液重力组成,三相接触线以下冲击区的熔池液面为二次曲线,如示意图3所示;支撑力总和模型表示为:
Figure BDA0003826736880000067
其中,Fσ为自由钢液面熔池表面张力,Fg为水射流没入钢液面以下排开的钢液重力;
Figure BDA0003826736880000068
为变形液面轮廓曲线三相接触线处的倾斜角,所述三相为钢液、空气和冷却水三相;ρ为钢液的密度;rm为冲击区最大直径;g为重力加速度;σ为钢液面的表面张力;hm为射流冲击区的高度。
步骤三之以变形液面轮廓曲线在三相接触线处的倾斜角
Figure BDA0003826736880000069
为临界条件,构建射流冲击刺破模型具体包括以下子步骤:
在射流达稳态条件下,射流冲击自由钢液面的冲击力总和等于自由钢液面对射流的支撑力总和;
Figure BDA0003826736880000071
的临界条件下,自由钢液面对射流的支撑力总和达到最大值得到射流冲击刺破模型:
F冲击力=F支撑力
Figure BDA0003826736880000072
式中的vj为下落高度H下的临界射流最大速度。
步骤二:分析水射流的液滴冲击自由钢液面的稳态受力情况和液面变形情况,获取支撑力模型;
在下落冲击自由钢液面过程中,液滴或液团的形状时刻变化,向“扁平化”发展过程中,无法准确获知其粒径,更无法得到跟熔池表面发生作用的相接触线。同时液滴直径大于6mm的液滴在下落过程中容易分裂。
故做如下假设:将液滴或液团等效为等体积、等质量的球形固态物体来接触熔池液面,且球形固态物体下落过程中忽略空气阻力的影响。在这种条件下,液滴或液团受到的熔池液面支撑力≤实际液滴变形状态下的支撑力,球形固态物体的假设具有保守性。
步骤二包括以下子步骤:
T21、将水射流的液滴等效为等体积、等质量的球形固态物体,分析球形固态物体的液滴冲击自由钢液面的受力情况和液面变形情况;
液滴以“球形固态物体”接触熔池表面后,自由钢液面变形及受力情况如图4所示,通过分析表面张力沿竖直方向分力与静水压力之和,来推导液滴总的支撑力计算公式。
由图4几何关系可得:
Figure BDA0003826736880000073
式中β—浸没角;rt—球形固态物体的球半径;rlr—三相接触线与液滴中心沿水平方向距离;hce—球形固态物体的球心距离自由水平面位移;θc—熔池液面在球形固态物体处的静态接触角。
T22、通过自由钢液面的表面张力沿竖直方向分力与水射流的液滴没入自由钢液面以下排开的钢液重力之和,推导出自由钢液面对液滴的的支撑力模型。
则自由钢液面对液滴的的支撑力模型的支撑力模型为:
Figure BDA0003826736880000081
β为浸没角,rt为球形固态物体的球半径;ρ为钢液的密度;g为重力加速度;hlr为三相接触线与熔池自由液面之间的距离;
Figure BDA0003826736880000082
为变形液面轮廓曲线在三相接触线处的倾斜角。
步骤三之以变形液面轮廓曲线在三相接触线处的倾斜角
Figure BDA0003826736880000083
为临界条件,构建液滴冲击刺破模型;
如图4所示,液滴(球形固态物体)冲击熔池自由液面的过程中,随着球形固态物体压迫熔池液面不断下陷,浸没角不断增大,变形液面轮廓曲线在三相接触线处的倾斜角
Figure BDA0003826736880000084
也不断增大,所以表面张力竖直方向的分力也不断增大,直到倾斜角
Figure BDA0003826736880000085
为90°时达到最大,此时熔池表面对液球形固态物体的支撑力也达到最大值。经过此位置之后,球形固态物体的下陷将更加容易,且此时的浸没角也已大于90°(钢液静态接触角大于90°),球形固态物体处于即将被完全浸没的状态。在准静态条件下,变形液面轮廓曲线在三相接触线处的倾斜角
Figure BDA0003826736880000086
达到90°即为刺破自由液面的临界条件。液滴刺破熔池自由液面临界条件下的下陷深度为:
Figure BDA0003826736880000087
以高速液滴、液团从一定高处z以初速v0落入到熔池自由液面处,并冲击液面致使变形下陷,直到液滴冲击速度降为零,此时恰好熔池液面对液滴的支撑力达到最大值,则此时的液滴初速v0即为高度z下的临界速度。液滴冲击熔池自由液面过程中,熔池液面对液滴的支撑力所做的负功即为液滴初始所具有的机械能,则得到液滴冲击刺破模型为:
Figure BDA0003826736880000088
其中,,hcr为球形固态物体刺破自由钢液面的下陷深度;θc为自由钢液面在球形固态物体处的静态接触角;z为液滴的下落高度;ρ表示液滴的密度;v0为液滴从下落处最高点的初始速度;
Figure BDA0003826736880000089
式中A=2πσrt
Figure BDA0003826736880000091
v0为液滴从下落高处最高点的初始速度。
根据液滴冲击刺破模型,在给定的rt,z条件下,得到液滴刺破熔池液面的临界初速vcr
Figure BDA0003826736880000092
的临界条件下,根据式:
Figure BDA0003826736880000093
计算球形固态物体刺破熔池自由液面的临界下陷深度hcr
步骤四:根据射流冲击刺破模型计算射流刺破熔池液面的初始速度和下落高度;根据液滴冲击刺破模型计算液滴刺破熔池液面的初始速度和下落高度。
针对水射流的射流和液滴从顶部冲击熔池自由液面的物理过程,建立数理模型,确立冷却水刺破不锈钢熔池液面的临界条件,且经过实验研究数据验证了模型的精确性。使用和评价堆内注水来缓解金属层热聚焦效应的方案终于有了一款理论工具,可以计算出各种不同工况下的临界速度,堆内注水的缓解措施中流速的设计必须不超过此阈值。此外,此模型亦可以用于严重事故分析软件中,评价各种事故序列下顶部注水后的堆内滞留效果,助力新型严重事故缓解措施的研发。
实施例2
金属熔池的自由钢液面也具有一些特殊性,金属熔池较大的液相密度和表面张力使得其液面不容易被刺破,相比于水蝇在水面自由滑动而不沉,金属熔池液面有过之而无不及。钢液的密度为ρ=8523-0.8358t(kg/m3),是水的密度的7~8倍;金属键物质的表面张力最大,钢液的表面张力在1823K时为1.8~1.88N/m,是水的表面张力(20℃时为0.0728N/m)的25倍,故而刺破钢液面不会像刺破水面那么容易。
本实施例使用本发明提供的获取水射流刺破熔池液面临界条件的方法进行计算常温射流冲击1600℃自由波动钢液面的过程中,不同喷嘴高度下的临界流速,其结果如图5所示。
本实施例还在展开的不锈钢熔池顶部注水实验中,在0.5m的喷嘴高度下,以Φ10mm,Φ15mm的喷嘴形成约10m/s初速的射流去冲击自由波动的钢液面,均刺破钢液面,发生了剧烈的蒸汽爆炸。
以本发明方案中的计算值来看,水射流冲击钢液面的情景,0.5m注水喷嘴高度下,射流临界速度为11.9m/s,对实验值给出了较好的预测。
实施例3
本实施例使用本发明提供的获取水射流刺破熔池液面临界条件的方法进行计算水射流的液滴冲击钢液面过程。
已知炼钢炉中,氧化铝夹杂物和钢液的接触角为144°,在下落高度为0.5m的条件下,使用5.2节中的公式1-10到1-12进行计算常温液滴冲击1600℃自由波动钢液面的过程中,不同液滴直径下的临界下落初速,计算结果如图6所示。
从图6可以看出,10mm以下的不同粒径的液滴,其临界下落初速变化很大,10mm以上的液团,其临界下落初速变化趋缓,逐渐达到恒值。在中国核动力院钢熔池顶部注水冷却实验中,液滴喷淋、边缘射流注水方式产生的液滴粒径均为6mm,在边缘水膜中产生的液团粒径也在5~10mm范围内,液滴和液团下落初速均不超过10m/s。在实验过程中均未发现冷却水液滴或液团刺破钢液表面,也未发生爆炸,凝固的钢液表面平整。由图6看出,此粒径范围的液滴,其下落初速需达到25m/s以上才会刺破钢液表面,对于实验中10m/s以下初速的液滴无法刺破钢液表面给出了很好的预测。
实际液滴下落受到空气阻力的很大影响,暴雨最大液滴粒径6mm,落地速度8~9m/s,粒径0.5mm的毛毛细雨的落地速度为2m/s。可见,液滴0~6mm的范围内,不管高度多高,初速为零的自然下落液滴喷淋均达不到刺破钢液面的临界速度。而6~20mm的液滴,需要高于25m/s的初速才能刺破钢液面。所以,使用喷淋液滴进行钢熔池表面的排热冷却是合理安全的,此模型从理论上给出了推断。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种获取水射流刺破熔池液面临界条件的方法,其特征在于,包括步骤:
步骤一:分析水射流的射流冲击自由钢液面的稳态受力情况和液面变形情况,获取冲击力总和模型和支撑力总和模型;
步骤二:分析水射流的液滴冲击自由钢液面的稳态受力情况和液面变形情况,获取支撑力模型;
步骤三:以变形液面轮廓曲线在三相接触线处的倾斜角
Figure FDA0003826736870000011
为临界条件,构建射流冲击刺破模型和液滴冲击刺破模型;
步骤四:根据射流冲击刺破模型计算射流刺破熔池液面的初始速度和下落高度;根据液滴冲击刺破模型计算液滴刺破熔池液面的初始速度和下落高度。
2.根据权利要求1所述的一种获取水射流刺破熔池液面临界条件的方法,其特征在于,步骤一包括:
S11、分析水射流的射流冲击自由钢液面的冲击压力,对冲击区内的径向冲击力积分构建出冲击力总和模型;
S12、分析自由钢液面对射流的支撑力和液面变形情况构建出支撑力总和模型。
3.根据权利要求2所述的一种获取水射流刺破熔池液面临界条件的方法,其特征在于,所述冲击力总和模型为:
Figure FDA0003826736870000012
其中,vj为射流中心的最大流速;H为射流下落高度;ρl为射流的密度。
4.根据权利要求3所述的一种获取水射流刺破熔池液面临界条件的方法,其特征在于,所述支撑力总和模型为:
Figure FDA0003826736870000013
其中,Fσ为自由钢液面熔池表面张力,Fg为水射流没入钢液面以下排开的钢液重力;
Figure FDA0003826736870000014
为变形液面轮廓曲线三相接触线处的倾斜角,所述三相为钢液、空气和冷却水三相;ρ为钢液的密度;rm为冲击区最大直径;g为重力加速度;σ为钢液面的表面张力;hm为射流冲击区的高度。
5.根据权利要求4所述的一种获取水射流刺破熔池液面临界条件的方法,其特征在于,所述射流冲击刺破模型获取方法为:
在射流达稳态条件下,射流冲击自由钢液面的冲击力总和等于自由钢液面对射流的支撑力总和;
Figure FDA0003826736870000021
的临界条件下,自由钢液面对射流的支撑力总和达到最大值得到射流冲击刺破模型。
6.根据权利要求5所述的一种获取水射流刺破熔池液面临界条件的方法,其特征在于,所述射流冲击刺破模型为:
Figure FDA0003826736870000022
7.根据权利要求1所述的一种获取水射流刺破熔池液面临界条件的方法,其特征在于,步骤二包括以下子步骤:
T21、将水射流的液滴等效为等体积、等质量的球形固态物体,分析球形固态物体的液滴冲击自由钢液面的受力情况和液面变形情况;
T22、通过自由钢液面的表面张力沿竖直方向分力与水射流的液滴没入自由钢液面以下排开的钢液重力之和,推导出自由钢液面对液滴的支撑力模型。
8.根据权利要求7所述的一种获取水射流刺破熔池液面临界条件的方法,其特征在于,所述支撑力模型为:
Figure FDA0003826736870000023
β为浸没角,rt为球形固态物体的球半径;ρ为钢液的密度;g为重力加速度;hlr为三相接触线与熔池自由液面之间的距离;
Figure FDA0003826736870000024
为变形液面轮廓曲线在三相接触线处的倾斜角。
9.根据权利要求8所述的一种获取水射流刺破熔池液面临界条件的方法,其特征在于,所述液滴冲击刺破模型为:
Figure FDA0003826736870000025
其中,hcr为球形固态物体刺破自由钢液面的下陷深度;θc为自由钢液面在球形固态物体处的静态接触角;z为液滴的下落高度;ρ表示液滴的密度;v0为液滴从下落处最高点的初始速度;
Figure FDA0003826736870000026
式中A=2πσrt
Figure FDA0003826736870000031
10.根据权利要求9所述的一种获取水射流刺破熔池液面临界条件的方法,其特征在于,根据液滴冲击刺破模型,在给定的rt,z条件下,得到液滴刺破熔池液面的临界初速vcr
Figure FDA0003826736870000032
的临界条件下,根据式:
Figure FDA0003826736870000033
计算球形固态物体刺破熔池自由液面的临界下陷深度hcr
CN202211062257.0A 2022-08-31 2022-08-31 一种获取水射流刺破熔池液面临界条件的方法 Active CN115470723B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211062257.0A CN115470723B (zh) 2022-08-31 2022-08-31 一种获取水射流刺破熔池液面临界条件的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211062257.0A CN115470723B (zh) 2022-08-31 2022-08-31 一种获取水射流刺破熔池液面临界条件的方法

Publications (2)

Publication Number Publication Date
CN115470723A true CN115470723A (zh) 2022-12-13
CN115470723B CN115470723B (zh) 2023-04-28

Family

ID=84368502

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211062257.0A Active CN115470723B (zh) 2022-08-31 2022-08-31 一种获取水射流刺破熔池液面临界条件的方法

Country Status (1)

Country Link
CN (1) CN115470723B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107609288A (zh) * 2017-09-21 2018-01-19 中国科学院力学研究所 模拟激光电弧复合焊接中熔滴对熔池冲击作用的计算方法
CN111386577A (zh) * 2018-11-01 2020-07-07 原子能设计股份公司 核反应堆堆芯熔化物冷却方法及核反应堆堆芯熔化物冷却控制系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107609288A (zh) * 2017-09-21 2018-01-19 中国科学院力学研究所 模拟激光电弧复合焊接中熔滴对熔池冲击作用的计算方法
CN111386577A (zh) * 2018-11-01 2020-07-07 原子能设计股份公司 核反应堆堆芯熔化物冷却方法及核反应堆堆芯熔化物冷却控制系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
崔杰等: "《射流冲击模型在水下爆炸实验中的应用研究》" *
李明磊: "《层流液体射流碰撞非融合现象及其临界特性的研究》" *
罗朝霞等: "《液滴冲击无限大页面过程的边界元模拟》" *
胡钰文等: "《熔融池相变热传特性的大涡模拟数值研究》" *

Also Published As

Publication number Publication date
CN115470723B (zh) 2023-04-28

Similar Documents

Publication Publication Date Title
Sanchez-Perez et al. A physical model for the two-phase flow in a continuous casting mold
Cho et al. Bubble behavior and size distributions in stopper-rod nozzle and mold during continuous casting of steel slabs
Chattopadhyay et al. Physical and mathematical modelling to study the effect of ladle shroud mis-alignment on liquid metal quality in a tundish
Liang et al. A study of a single liquid drop impact on inclined wetted surfaces
Ling et al. Multiphase flow behavior in a single-strand continuous casting tundish during ladle change
CN115470723A (zh) 一种获取水射流刺破熔池液面临界条件的方法
Sibamoto et al. Visualization and measurement of subcooled water jet injection into high-temperature melt by using high-frame-rate neutron radiography
Spencer et al. Effect of boiling regime on melt stream breakup in water
Zhang et al. Effects of gas velocity and break size on steam penetration depth using gas jet into water similarity experiments
CN111341469A (zh) 一种应用于核电站压力容器冷却的浸没式喷射系统
Zhu et al. Experimental study of bubble growth and departure at the tip of capillary tubes with various wettabilities in a stagnant liquid
Ribeiro et al. Insights on bubble encapsulation after drop impact on thin liquid films
Wang et al. Dynamics of cavitation bubbles in viscous liquids in a tube during a transient process
Wendel et al. Experiments and Simulations With Large Gas Bubbles in Mercury Towards Establishing a Gas Layer to Mitigate Cavitation Damage
Senthilkumar et al. Optimization of spillway weir to mitigate gas entrainment in a fast reactor
Alam et al. Jet Formation After Droplet Impact on Microholed Hydrophilic Surfaces
Chung et al. Collapse of a hot vapor bubble in subcooled liquid
Wu et al. Numerical Investigation on the Oblique Impact of Droplets on Moving Film in Spray Cooling
CN116078760B (zh) 用于反应堆压力容器外壁面的汽泡清扫搅混装置及方法
Spencer et al. Phenomenological investigations of cavity interactions following postulated vessel meltthrough
Nagai Shock pressures exerted by breaking waves on breakwaters
Haas et al. A Review of Bubble Dynamics in Liquid Metals. Metals 2021, 11, 664
Anderson et al. Experimental Studies of High-Speed Liquid Films on Downward-Facing Surfaces for IFE Wet Wall Concepts
Van Tichelen et al. MYRRHA: Design and verification experiments for the windowless spallation target of the ADS Prototype MYRRHA
Sinha et al. Droplet Impact on a Curved Moving Liquid Film

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