CN107563030A - 一种针对两种流体传热交混破碎相变过程的无网格模拟方法 - Google Patents
一种针对两种流体传热交混破碎相变过程的无网格模拟方法 Download PDFInfo
- Publication number
- CN107563030A CN107563030A CN201710725986.2A CN201710725986A CN107563030A CN 107563030 A CN107563030 A CN 107563030A CN 201710725986 A CN201710725986 A CN 201710725986A CN 107563030 A CN107563030 A CN 107563030A
- Authority
- CN
- China
- Prior art keywords
- particle
- mrow
- msub
- fragment
- moment
- 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.)
- Pending
Links
Abstract
本发明提供的是一种针对两种流体传热交混破碎相变过程的无网格模拟方法。包括凝固相变模型、物性参数变化模型、表面张力模型、碎片结合判据以及流固混合模型。采用基于焓方法的凝固相变模型,对相变问题进行了简化;物性参数变化模型对熔融物在凝固破碎过程中的传热及流动物性参数进行修正,使计算更加准确;采用粒子间相互作用的表面张力模型,能够更简单地与无网格法进行结合;采用合理的碎片结合判据确定破碎产生碎片的粒子组成;采用流固混合模型对凝固后碎片的运动进行模拟。本方法基于移动粒子半隐式法,提出一种改进的无网格方法,综合各个模型可实现两种流体传热交混破碎相变过程,尤其针对两相流动、凝固传热以及流固混合等复杂问题的模拟计算求解。
Description
技术领域
本发明涉及的是多相流传热相变模拟方法,具体地说是一种针对两种流体传热交混破碎相变过程的无网格模拟方法。
背景技术
金属燃料堆芯破损事故及其导致的熔融金属材料与冷却剂接触的相互作用FCI(Fuel-Coolant Interaction),是国际上金属燃料钠冷快堆堆芯安全研究面临的一个重点和难点问题。在FCI过程中,两种或两种以上的热熔融物质与冷却剂液态钠之间的巨大温差,强化了冷热流体间的多相传热传质,使得流体接触界面发生剧烈相变,并伴随能量的瞬间释放,以及熔融物质的变形、破碎和凝固相变等,使得FCI过程极为复杂。
目前国内关于金属燃料与冷却剂钠相互作用的模拟程序比较缺乏,而国外的模拟程序主要关注FCI过程中,温度变化以及产生的压力脉冲,对于熔融金属材料的凝固破碎尺寸并没有太多的模拟研究。传统网格面对此类存在自由界面的流体传热交混破碎相变问题时会出现网格大变形以及数值扩散等问题,而无网格方法可以很好的解决这一问题。
发明内容
本发明的目的在于提供一种可实现两种流体传热交混破碎相变过程,尤其针对两相流动、凝固传热以及流固混合等复杂问题的模拟的针对两种流体传热交混破碎相变过程的无网格模拟方法。
本发明的目的是这样实现的:
步骤一,设定时间步长Δt、粒子间初始距离d0、粒子的控制半径Re以及总的模拟时间ttotal,设定tn=n·Δt,则在t0时刻输入初始化参数,包括熔融金属流体的初始温度Th、速度V0、形态,冷却流体的初始温度Tc;
步骤二,在tn时刻采用基于焓方法的凝固相变模型进行传热计算,n>0,根据tn-1时刻的粒子温度Tn-1、焓值hn-1、位置rn-1、密度ρn-1以及热导率Kn-1,计算得到tn时刻粒子的温度Tn、焓值hn以及液相分数αn;
步骤三,采用物性参数变化模型对处于相变区域粒子的物性参数进行修正,根据粒子的液相分数αn,更新得到tn时刻粒子的密度ρn、热导率Kn以及粘性系数vn;
步骤四,采用基于表面自由能的表面张力模型计算tn时刻粒子所受的表面张力
步骤五,根据步骤三和步骤四所得的粒子物性参数ρn、vn和表面张力采用原始MPS方法计算出tn时刻所有粒子的位置rn和速度un,分为以下几个小步骤:
1,不考虑压力梯度项,显式求解动量方程中的粘性项、表面张力项以及重力项,计算得出tn时刻粒子的中间速度u*和中间位置r*;
2,计算压力泊松方程,得出tn时刻粒子的压力值Pn;
3,利用粒子的压力值Pn,对粒子速度进行修正,得到tn时刻粒子的速度un和位置rn;
步骤六,根据步骤五中计算所得粒子位置rn,采用碎片结合判据确定破碎相变产生碎片的粒子组成;包括液相分数判据和粒子间距离判据,只有两个判据同时满足,才能判定两个粒子结合成碎片;液相分数判据为,两个熔融金属粒子的液相分数均需满足α1<α<α2,其中α1、α2是两个可变参数且0≤α1<α2≤1,根据具体模拟情况进行调整;粒子间距离判据为,两个熔融金属粒子之间的距离需小于1.2l0;
步骤七,采用流固混合模型(PMS方法)修正tn时刻碎片中组成粒子的位置rn和速度un,分为以下几个步骤:
1,计算碎片的速度和位置
2,计算碎片的转动惯量
3,基于角动量守恒定律,计算碎片的角速度
4,计算tn时刻碎片中粒子的位置rn和速度un;
步骤八,判断是否达到设定的总的模拟计算时间ttotal,若tn<ttotal,则跳至步骤二,进行tn+1时刻的计算;若tn=ttotal,则完成模拟。
本发明还可以包括:
1、步骤二中,焓方法的传热方程为:
2、步骤二中,粒子的液相分数定义如下:
其中下标s和l分别代表固态以及液态。
3、步骤三中,处于相变区域的粒子的密度、热导率根据液相分数加权计算,表示如下:
ρ=ρs*(1-α)+ρl*α
K=Ks*(1-α)+Kl*α。
4、步骤三中,处于相变区域的粒子的粘性系数根据液相分数以指数形式变化,表示如下:
v=vl·e5(1-α)。
5、步骤四中,表面张力模型是基于表面自由能,且以粒子间相互作用力的形式呈现,当两个粒子间的距离在(0,Re)范围内,则粒子间相互作用力表示如下:
f=C·(r-l0)·(r-Re)/m。
6、步骤六中,碎片结合判据是用来确定破碎凝固产生碎片的粒子组成,包括液相分数判据,以及粒子间距离判据,只有两个判据同时满足才能判定两个粒子结合成碎片;液相分数判据为,两个熔融金属粒子的液相分数均需满足α1<α<α2,其中α1、α2是两个可变参数且0≤α1<α2≤1,根据具体模拟情况进行调整;粒子间距离判据为,相邻的两个金属粒子之间的距离需小于ε·l0,其中l0是粒子间初始距离,ε是一个系数,根据初始布置中一个粒子i和其周围8个粒子之间距离的平均值,推算出ε≈1.2。
7、步骤七中,流固混合模型是基于角动量守恒定律,对尺寸不断变化的碎片的位置和速度进行计算。
本发明针对堆芯熔融金属破碎凝固的模拟方法的不足,提供了一种可实现两种流体传热交混破碎相变过程,尤其针对两相流动、凝固传热以及流固混合等复杂问题的模拟方法。本发明基于无网格移动粒子半隐式法提出了一种针对两种流体传热交混破碎相变过程的无网格模拟方法(Multiphase Moving Particle Semi-implicit method forFragmentation and Phase Transition simulation,MMPS-FPT)。
采用基于焓方法的凝固相变模型,对相变问题进行了简化;物性参数变化模型对熔融物在凝固破碎过程中的传热及流动物性参数进行修正,使计算更加准确;采用粒子间相互作用的表面张力模型,能够更简单地与无网格法进行结合;采用合理的碎片结合判据确定破碎产生碎片的粒子组成;采用流固混合模型对凝固后碎片的运动进行模拟。
本方法基于移动粒子半隐式法(MPS),提出一种改进的无网格方法,综合各个模型可实现两种流体传热交混破碎相变过程,尤其针对两相流动、凝固传热以及流固混合等复杂问题的模拟计算求解。
附图说明
图1本发明两种流体传热交混破碎相变过程的无网格模拟方法流程图;
图2本发明金属粒子温度、焓值及液相分数之间对应关系图。
图3本发明中碎片结合判据中的粒子间距离判据原理。
具体实施方式
下面结合附图对本发明及其实施方式作进一步详细描述。
如图1所示,本发明实现如下:
步骤一,对于一个特定工况的模拟计算,设定计算中使用的时间步长Δt、粒子间初始距离l0、粒子的控制半径Re以及总的模拟计算时间ttotal。规定tn=n·Δt,则在t0时刻输入初始化参数,包括熔融金属流体的初始温度Th、速度V0、形态,冷却流体的初始温度Tc。
步骤二:在tn时刻(n>0)采用基于焓方法的凝固相变模型进行传热计算。根据tn-1时刻的粒子温度Tn-1、焓值hn-1、位置rn-1、密度ρn-1以及热导率Kn-1等,计算得到tn时刻粒子的温度Tn、焓值hn以及液相分数αn。具体计算表示如下:
其中Cps,Cpl,L,Tmp分别代表固态定压比热、液态定压比热、潜热以及熔点温度,下标s和l分别代表固态以及液态。
图2所示为本发明金属粒子温度、焓值及液相分数之间对应关系图。焓值和温度的对应关系表示如下:
步骤三,采用物性参数变化模型对处于相变区域粒子的物性参数进行修正,根据粒子的液相分数αn,更新得到tn时刻粒子的密度ρn、热导率Kn以及粘性系数vn。表示如下:
ρn=ρs*(1-αn)+ρl*αn
Kn=Ks*(1-αn)+Kl*αn
步骤四,采用基于表面自由能的表面张力模型计算tn时刻粒子所受的表面张力以粒子间相互作用的形式呈现。当两个粒子间的距离在(0,Re)范围内,则粒子间相互作用了表示如下:
其中C为系数,r,m,σ分别为粒子间距离、粒子的质量和表面张力系数。
步骤五,根据步骤三和步骤四所得的粒子物性参数ρn、vn和表面张力采用原始MPS方法计算出tn时刻所有粒子的位置rn和速度un。可分为以下几个小步骤:1,不考虑压力梯度项,显式求解动量方程中的粘性项表面张力项以及重力项g,计算得出tn时刻粒子的中间速度u*和中间位置r*;2,计算压力泊松方程,得出tn时刻粒子的压力值Pn;3,利用粒子的压力值Pn,对粒子速度进行修正,得到tn时刻粒子的速度un和位置rn。具体计算流程如下:
r*=rn-1+Δt·u*
rn=rn-1+Δt·un
步骤六,根据步骤五中计算所得粒子位置rn,采用碎片结合判据确定破碎相变产生碎片的粒子组成。具体包括液相分数判据和粒子间距离判据,只有两个判据同时满足才能判定两个粒子结合成碎片。液相分数判据为,两个熔融金属粒子的液相分数均需满足α1<α<α2,其中α1、α2是两个可变参数且0≤α1<α2≤1,可根据具体模拟情况进行调整。图3所示为本发明中碎片结合判据中的粒子间距离判据原理,相邻的两个金属粒子之间的距离需小于ε·l0。其中l0是粒子间初始距离,ε是一个系数。根据初始布置中一个粒子i和其周围8个粒子之间距离的平均值,可以推算出
步骤七,采用流固混合模型(PMS方法)修正tn时刻碎片中组成粒子的位置rn和速度un。可分为以下几个步骤:1,计算碎片的速度和位置2,计算碎片的转动惯量3,基于角动量守恒定律,计算碎片的角速度4,计算tn时刻碎片中粒子的位置rn和速度un。对于一个包含N个粒子的碎片,具体计算表示如下:
其中,下标f代表碎片,下标i代表粒子编号,M为旋转矩阵。
步骤八,判断是否达到设定的总的模拟计算时间ttotal。若tn<ttotal,则跳至步骤二,进行tn+1时刻的计算;若tn=ttotal,则完成计算。
需要说明的是,以上介绍的是本发明的实施方案而并非限制。任何对本发明技术法案的修改或者等同替代都不脱离本发明技术法案的精神和范围,均应该涵盖在本发明的权利要求范围之内。
Claims (8)
1.一种针对两种流体传热交混破碎相变过程的无网格模拟方法,其特征是:
步骤一,设定时间步长Δt、粒子间初始距离d0、粒子的控制半径Re以及总的模拟时间ttotal,设定tn=n·Δt,则在t0时刻输入初始化参数,包括熔融金属流体的初始温度Th、速度V0、形态,冷却流体的初始温度Tc;
步骤二,在tn时刻采用基于焓方法的凝固相变模型进行传热计算,n>0,根据tn-1时刻的粒子温度Tn-1、焓值hn-1、位置rn-1、密度ρn-1以及热导率Kn-1,计算得到tn时刻粒子的温度Tn、焓值hn以及液相分数αn;
步骤三,采用物性参数变化模型对处于相变区域粒子的物性参数进行修正,根据粒子的液相分数αn,更新得到tn时刻粒子的密度ρn、热导率Kn以及粘性系数vn;
步骤四,采用基于表面自由能的表面张力模型计算tn时刻粒子所受的表面张力
步骤五,根据步骤三和步骤四所得的粒子物性参数ρn、vn和表面张力采用原始MPS方法计算出tn时刻所有粒子的位置rn和速度un,分为以下几个小步骤:
1,不考虑压力梯度项,显式求解动量方程中的粘性项、表面张力项以及重力项,计算得出tn时刻粒子的中间速度u*和中间位置r*;
2,计算压力泊松方程,得出tn时刻粒子的压力值Pn;
3,利用粒子的压力值Pn,对粒子速度进行修正,得到tn时刻粒子的速度un和位置rn;
步骤六,根据步骤五中计算所得粒子位置rn,采用碎片结合判据确定破碎相变产生碎片的粒子组成;包括液相分数判据和粒子间距离判据,只有两个判据同时满足,才能判定两个粒子结合成碎片;液相分数判据为,两个熔融金属粒子的液相分数均需满足α1<α<α2,其中α1、α2是两个可变参数且0≤α1<α2≤1,根据具体模拟情况进行调整;粒子间距离判据为,两个熔融金属粒子之间的距离需小于1.2l0;
步骤七,采用流固混合模型修正tn时刻碎片中组成粒子的位置rn和速度un,分为以下几个步骤:
1,计算碎片的速度和位置
2,计算碎片的转动惯量
3,基于角动量守恒定律,计算碎片的角速度
4,计算tn时刻碎片中粒子的位置rn和速度un;
步骤八,判断是否达到设定的总的模拟计算时间ttotal,若tn<ttotal,则跳至步骤二,进行tn+1时刻的计算;若tn=ttotal,则完成模拟。
2.根据权利要求1所述的针对两种流体传热交混破碎相变过程的无网格模拟方法,其特征是,步骤二中,焓方法的传热方程为:
<mrow>
<mfrac>
<mrow>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>D</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>&dtri;</mo>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>&dtri;</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
3.根据权利要求2所述的针对两种流体传热交混破碎相变过程的无网格模拟方法,其特征是,步骤二中,粒子的液相分数定义如下:
<mrow>
<mi>&alpha;</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo><</mo>
<msub>
<mi>h</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>h</mi>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>s</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mi>l</mi>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>h</mi>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>s</mi>
</msub>
</mrow>
<mi>L</mi>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mi>s</mi>
</msub>
<mo><</mo>
<mi>h</mi>
<mo><</mo>
<msub>
<mi>h</mi>
<mi>l</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mi>l</mi>
</msub>
<mo><</mo>
<mi>h</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中下标s和l分别代表固态以及液态。
4.根据权利要求3所述的针对两种流体传热交混破碎相变过程的无网格模拟方法,其特征是,步骤三中,处于相变区域的粒子的密度、热导率根据液相分数加权计算,表示如下:
ρ=ρs*(1-α)+ρl*α
K=Ks*(1-α)+Kl*α。
5.根据权利要求4所述的针对两种流体传热交混破碎相变过程的无网格模拟方法,其特征是,步骤三中,处于相变区域的粒子的粘性系数根据液相分数以指数形式变化,表示如下:
v=vl·e5(1-α)。
6.根据权利要求5所述的针对两种流体传热交混破碎相变过程的无网格模拟方法,其特征是,步骤四中,表面张力模型是基于表面自由能,且以粒子间相互作用力的形式呈现,当两个粒子间的距离在(0,Re)范围内,则粒子间相互作用力表示如下:
f=C·(r-l0)·(r-Re)/m。
7.根据权利要求6所述的针对两种流体传热交混破碎相变过程的无网格模拟方法,其特征是,步骤六中,碎片结合判据是用来确定破碎凝固产生碎片的粒子组成,包括液相分数判据,以及粒子间距离判据,只有两个判据同时满足才能判定两个粒子结合成碎片;液相分数判据为,两个熔融金属粒子的液相分数均需满足α1<α<α2,其中α1、α2是两个可变参数且0≤α1<α2≤1,根据具体模拟情况进行调整;粒子间距离判据为,相邻的两个金属粒子之间的距离需小于ε·l0,其中l0是粒子间初始距离,ε是一个系数,根据初始布置中一个粒子i和其周围8个粒子之间距离的平均值,推算出ε≈1.2。
8.根据权利要求7所述的针对两种流体传热交混破碎相变过程的无网格模拟方法,其特征是,步骤七中,流固混合模型是基于角动量守恒定律,对尺寸不断变化的碎片的位置和速度进行计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710725986.2A CN107563030A (zh) | 2017-08-22 | 2017-08-22 | 一种针对两种流体传热交混破碎相变过程的无网格模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710725986.2A CN107563030A (zh) | 2017-08-22 | 2017-08-22 | 一种针对两种流体传热交混破碎相变过程的无网格模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107563030A true CN107563030A (zh) | 2018-01-09 |
Family
ID=60976307
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710725986.2A Pending CN107563030A (zh) | 2017-08-22 | 2017-08-22 | 一种针对两种流体传热交混破碎相变过程的无网格模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107563030A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109215100A (zh) * | 2018-08-08 | 2019-01-15 | 天津大学 | 一种混合流体相变动画生成方法及装置 |
CN110044959A (zh) * | 2019-05-13 | 2019-07-23 | 西安交通大学 | 利用移动粒子有限容积法研究高瑞利数熔融池换热特性的方法 |
CN110070919A (zh) * | 2019-04-12 | 2019-07-30 | 上海交通大学 | 一种涉及晶相反应的熔化模型及其数值模拟方法 |
CN110321641A (zh) * | 2019-07-08 | 2019-10-11 | 西安交通大学 | 基于粒子法的熔融物与混凝土相互作用分析方法 |
CN111832214A (zh) * | 2020-06-29 | 2020-10-27 | 西安交通大学 | 基于粒子法的核反应堆严重事故碎片床熔化过程模拟方法 |
CN113191066A (zh) * | 2021-04-30 | 2021-07-30 | 西安交通大学 | 基于无网格法的核反应堆燃料元件失效分析方法 |
CN115048848A (zh) * | 2022-05-31 | 2022-09-13 | 西安交通大学 | 一种钠冷快堆熔融物碎片化评价方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101482475A (zh) * | 2008-01-11 | 2009-07-15 | 屈雪利 | 一种采集微通道中影响流动和传质因素相关参数的方法 |
CN102495968A (zh) * | 2011-12-14 | 2012-06-13 | 中国科学院工程热物理研究所 | 重力条件下高温蓄热容器内固/液相变的数值模拟方法 |
CN103279630A (zh) * | 2013-06-19 | 2013-09-04 | 天津大学 | 激光深熔点焊匙孔动态成形数值模拟方法 |
CN104200015A (zh) * | 2014-08-20 | 2014-12-10 | 清华大学 | 一种流体模拟方法及装置 |
CN104461677A (zh) * | 2014-10-30 | 2015-03-25 | 中国运载火箭技术研究院 | 一种基于cfd和fem技术的虚拟热试验方法 |
-
2017
- 2017-08-22 CN CN201710725986.2A patent/CN107563030A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101482475A (zh) * | 2008-01-11 | 2009-07-15 | 屈雪利 | 一种采集微通道中影响流动和传质因素相关参数的方法 |
CN102495968A (zh) * | 2011-12-14 | 2012-06-13 | 中国科学院工程热物理研究所 | 重力条件下高温蓄热容器内固/液相变的数值模拟方法 |
CN103279630A (zh) * | 2013-06-19 | 2013-09-04 | 天津大学 | 激光深熔点焊匙孔动态成形数值模拟方法 |
CN104200015A (zh) * | 2014-08-20 | 2014-12-10 | 清华大学 | 一种流体模拟方法及装置 |
CN104461677A (zh) * | 2014-10-30 | 2015-03-25 | 中国运载火箭技术研究院 | 一种基于cfd和fem技术的虚拟热试验方法 |
Non-Patent Citations (2)
Title |
---|
YANG ZHI 等: "Numerical simulation of a single molten 316SS droplet penetrating into sodium pool with MPS method", 《NUCLEAR SAFETY AND SIMULATION》 * |
ZHI-GANG ZHANG: "Experimental study on combustion characteristics of sodium fire in a columnar flow", 《JOURNAL OF NUCLEAR SCIENCE AND TECHNOLOGY》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109215100A (zh) * | 2018-08-08 | 2019-01-15 | 天津大学 | 一种混合流体相变动画生成方法及装置 |
CN110070919A (zh) * | 2019-04-12 | 2019-07-30 | 上海交通大学 | 一种涉及晶相反应的熔化模型及其数值模拟方法 |
CN110070919B (zh) * | 2019-04-12 | 2023-02-17 | 上海交通大学 | 一种涉及晶相反应的熔化模型及其数值模拟方法 |
CN110044959A (zh) * | 2019-05-13 | 2019-07-23 | 西安交通大学 | 利用移动粒子有限容积法研究高瑞利数熔融池换热特性的方法 |
CN110321641A (zh) * | 2019-07-08 | 2019-10-11 | 西安交通大学 | 基于粒子法的熔融物与混凝土相互作用分析方法 |
CN110321641B (zh) * | 2019-07-08 | 2020-08-04 | 西安交通大学 | 基于粒子法的熔融物与混凝土相互作用分析方法 |
CN111832214A (zh) * | 2020-06-29 | 2020-10-27 | 西安交通大学 | 基于粒子法的核反应堆严重事故碎片床熔化过程模拟方法 |
CN113191066A (zh) * | 2021-04-30 | 2021-07-30 | 西安交通大学 | 基于无网格法的核反应堆燃料元件失效分析方法 |
CN113191066B (zh) * | 2021-04-30 | 2022-12-09 | 西安交通大学 | 基于无网格法的核反应堆燃料元件失效分析方法 |
CN115048848A (zh) * | 2022-05-31 | 2022-09-13 | 西安交通大学 | 一种钠冷快堆熔融物碎片化评价方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107563030A (zh) | 一种针对两种流体传热交混破碎相变过程的无网格模拟方法 | |
Xiao et al. | Effect of metal foam on improving solid–liquid phase change in a multi-channel thermal storage tank | |
Kelkar et al. | Computational modeling of electroslag remelting (ESR) process used for the production of high-performance alloys | |
CN105302985B (zh) | 一种基于fluent软件的合金微铸造成形过程的仿真方法 | |
Zhang et al. | Atomic-scale structural evolution in selective laser melting of Cu50Zr50 metallic glass | |
CN110321641A (zh) | 基于粒子法的熔融物与混凝土相互作用分析方法 | |
CN109817284A (zh) | 一种钢液中枝晶移动的预测方法 | |
CN104881588B (zh) | 铸锭宏观偏析数值模拟方法 | |
Zhu et al. | Modelling of dendritic growth during alloy solidification under natural convection | |
Ehlen et al. | Simulation of time-dependent pool shape during laser spot welding: Transient effects | |
Chang et al. | Modeling on dry centrifugal granulation process of molten blast furnace slag | |
Liu et al. | A cellular automaton-lattice Boltzmann method for modeling growth and settlement of the dendrites for Al-4.7% Cu solidification | |
CN107092754A (zh) | 一种合金晶粒组织数值预测方法 | |
CN105665684B (zh) | 一种铸件晶粒组织数值预测的方法 | |
Kaminski et al. | Early accretion of planetesimals unraveled by the thermal evolution of the parent bodies of magmatic iron meteorites | |
Arid et al. | Numerical simulation of ice melting near the density inversion point under periodic thermal boundary conditions | |
CN107301291A (zh) | 一种基于电渣重熔过程的电磁场、温度场和流动场耦合有限元计算方法 | |
Xu et al. | Multiscale Modeling and Simulation of Directional Solidification Process of Turbine Blade Casting with MCA Method | |
Yi et al. | Gas entrapment and pore formation in metal droplet-based 3D printing | |
CN106944607B (zh) | 一种孕育合金晶粒组织数值预测方法 | |
CN105868461A (zh) | 一种激光熔覆多组分合金涂层的成分计算方法 | |
Phillips et al. | Structure and function of window glass and Pyrex | |
CN102236726B (zh) | 一种金属液凝固过程中预测缩松的方法及缩松连续预测方法 | |
Zhai et al. | Numerical simulation studies of jet rapid solidification technology for magnetic materials: a review | |
Chen et al. | Finite element numerical simulation of temperature field in metal pattern casting system and “reverse method” of defining the thermal physical coefficient |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180109 |
|
RJ01 | Rejection of invention patent application after publication |