搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

枝晶生长和气泡形成的数值模拟

吴伟 孙东科 戴挺 朱鸣芳

引用本文:
Citation:

枝晶生长和气泡形成的数值模拟

吴伟, 孙东科, 戴挺, 朱鸣芳

Modeling of dendritic growth and bubble formation

Wu Wei, Sun Dong-Ke, Dai Ting, Zhu Ming-Fang
PDF
导出引用
  • 本文建立了二维的格子玻尔兹曼方法-元胞自动机(lattice Boltzmann method-cellular automaton, LBM-CA)耦合模型, 对凝固过程中枝晶生长和气泡形成进行模拟研究. 本模型采用CA方法模拟枝晶的生长, 根据界面溶质平衡法计算枝晶生长的驱动力. 采用基于Shan-Chen多相流的LBM模拟气泡在液相中的生长和运动. 在LBM-CA的耦合模型中包含了固-液-气三相之间的相互作用. 应用Laplace定理和模拟气-液-固三相之间的润湿现象对模型进行了验证. 应用所建立的LBM-CA耦合模型模拟研究了气-液相互作用系数对单气泡生长的影响. 发现单气泡的生长速度和平衡半径随气-液相互作用系数的增大而增大. 定向凝固过程中枝晶和气泡生长的模拟结果再现了枝晶的择优生长、 气泡的优先形核位置、气泡的长大、合并、在枝晶间受挤变形以及在液相通道中的运动等物理现象, 与实验结果符合良好. 此外, 初始气体含量越高, 凝固结束时气泡的体积分数也相对较高. 本模型的模拟结果可以揭示在凝固过程中气泡形核、 生长和运动演化以及与枝晶生长相互作用的物理机理.
    A two-dimensional (2D) lattice Boltzmann method (LBM)-cellular automaton (CA) coupled model is developed for the simulation of dendritic growth and bubble formation during solidification. In the model, the dendritic growth is simulated by a CA approach. The driving force of dendritic growth is determined by a local solute equilibrium approach. The LBM based on the Shan-Chen multiphase flow scheme is adopted to simulate the growth and the motion of bubbles in liquid. The interaction mechanism between dendrites and bubbles is embedded in the model. Model validation is carried out by comparing the simulations with the Laplace law, and by simulating the wettability of a bubble on a smooth solid surface. The proposed model is used to study the effect of gas-liquid interaction coefficient on single bubble growth. It is found that the growth velocity and the equilibrium radius of bubble increase with the gas-liquid interaction coefficient. The simulations of the dendritic growth and bubble formation during directional solidification reproduce the physical phenomena, including dendritic competitive growth, the preferential nucleation locations of bubbles, and bubble growth, coalescence, deformation due to the squeeze of neighboring dendrites, as well as bubble motion in the liquid channels. The simulation results are compared reasonably well with the experimental results. In addition, gas pore volume fraction increases with the initial gas content. The simulations of the present LBM-CA model provide an insight into the physical mechanism of bubble nucleation, growth, and motion, as well as the interaction between the dendritic growth and bubble formation during solidification.
    • 基金项目: 国家自然科学基金(批准号: 50971042)资助的课题.
    • Funds: Project supported by the National Natural Science Foudation of China (Grant No. 50971042).
    [1]

    Han Q Y 2006 Scripta Mater. 55 871

    [2]

    Zhao L, Liao H C, Pan Y, Wang L, Wang Q G 2011 Scripta Mater. 65 795

    [3]

    Atwood R C, Lee P D 2003 Acta Mater. 51 5447

    [4]

    Dong S Y, Xiong S M, Liu B C 2004 Mater Sci Technol. 20 23

    [5]

    Raabe D 2004 Modelling Simul. Mater. Sci. Eng. 12 R13

    [6]

    Guo Z L, Zheng C G 2009 Theory and Applications of Lattice Boltzmann Method (Beijing: Science Press) p9-10, 61-63 (in Chinese) [郭照立, 郑楚光 2009 格子Boltzmann方法的原理和应用(第一版) (北京: 科学出版出版社)第9—10, 61—63页]

    [7]

    Zeng J B, Li R J, Liao Q, Jiang F M 2011 Acta Phys.Sin. 60 066401 (in Chinese) [曾建邦, 李隆键, 廖全, 蒋方明 2011 物理学报 60 066401]

    [8]

    Wang W X, Shi J, Qiu B, Li H B 2010 Acta Phys. Sin. 59 8371 (in Chinese) [王文霞, 施娟, 邱冰, 李华兵 2010 物理学报 59 8371]

    [9]

    Wang J F, Liu Y, Xu Y S 2009 Biomed Microdevices 11 351

    [10]

    Guo Z L, Asinari, Pietro, Zheng C G 2009 Phys. Rev. E 79 026702

    [11]

    Sun D K, Zhu M F, Pan S Y, Raabe D 2009 Acta Mater. 6 1755

    [12]

    Yin H, Felicelli SD, Wang L 2011 Acta Mater. 59 3124

    [13]

    Selzer M, Jainta M, Nestler B 2009 Phys. Status Solid B 246 1197

    [14]

    Zhou F M, Sun D K, Zhu M F 2010 Acta Phys. Sin. 59 3394 (in Chinese) [周丰茂, 孙东科, 朱鸣芳 2010 物理学报 59 3394]

    [15]

    Li Q, MaY C, Liu K, Kang X H, Li D Z 2007 Acta Metall. Sin. 43 217 (in Chinese) [李强, 马颖澈, 刘奎, 康秀红, 李殿中 2007 金属学报 43 217]

    [16]

    Shan B W, Lin X, Wei L, Huang W D 2009 Acta Phys. Sin. 58 1132 (in Chinese) [单博炜, 林鑫, 魏雷, 黄卫东 2009 物理学报 58 1132]

    [17]

    Huang H B, Thorne D T, Jr, Schaap M G, Sukop M C 2007 Phys. Rev. E 76 066701

    [18]

    Shan X W 2006 Phys. Rev. E 73 047701

    [19]

    Zhu M F, Stefanescu D M 2007 Acta Mater. 55 1741

    [20]

    Sukop M C, Thorne D T, Jr 2005 Lattice Boltzmann Modeling-An Introduction for Geoscientists and Engineers (2nd ed) (New York: Springer) p114

  • [1]

    Han Q Y 2006 Scripta Mater. 55 871

    [2]

    Zhao L, Liao H C, Pan Y, Wang L, Wang Q G 2011 Scripta Mater. 65 795

    [3]

    Atwood R C, Lee P D 2003 Acta Mater. 51 5447

    [4]

    Dong S Y, Xiong S M, Liu B C 2004 Mater Sci Technol. 20 23

    [5]

    Raabe D 2004 Modelling Simul. Mater. Sci. Eng. 12 R13

    [6]

    Guo Z L, Zheng C G 2009 Theory and Applications of Lattice Boltzmann Method (Beijing: Science Press) p9-10, 61-63 (in Chinese) [郭照立, 郑楚光 2009 格子Boltzmann方法的原理和应用(第一版) (北京: 科学出版出版社)第9—10, 61—63页]

    [7]

    Zeng J B, Li R J, Liao Q, Jiang F M 2011 Acta Phys.Sin. 60 066401 (in Chinese) [曾建邦, 李隆键, 廖全, 蒋方明 2011 物理学报 60 066401]

    [8]

    Wang W X, Shi J, Qiu B, Li H B 2010 Acta Phys. Sin. 59 8371 (in Chinese) [王文霞, 施娟, 邱冰, 李华兵 2010 物理学报 59 8371]

    [9]

    Wang J F, Liu Y, Xu Y S 2009 Biomed Microdevices 11 351

    [10]

    Guo Z L, Asinari, Pietro, Zheng C G 2009 Phys. Rev. E 79 026702

    [11]

    Sun D K, Zhu M F, Pan S Y, Raabe D 2009 Acta Mater. 6 1755

    [12]

    Yin H, Felicelli SD, Wang L 2011 Acta Mater. 59 3124

    [13]

    Selzer M, Jainta M, Nestler B 2009 Phys. Status Solid B 246 1197

    [14]

    Zhou F M, Sun D K, Zhu M F 2010 Acta Phys. Sin. 59 3394 (in Chinese) [周丰茂, 孙东科, 朱鸣芳 2010 物理学报 59 3394]

    [15]

    Li Q, MaY C, Liu K, Kang X H, Li D Z 2007 Acta Metall. Sin. 43 217 (in Chinese) [李强, 马颖澈, 刘奎, 康秀红, 李殿中 2007 金属学报 43 217]

    [16]

    Shan B W, Lin X, Wei L, Huang W D 2009 Acta Phys. Sin. 58 1132 (in Chinese) [单博炜, 林鑫, 魏雷, 黄卫东 2009 物理学报 58 1132]

    [17]

    Huang H B, Thorne D T, Jr, Schaap M G, Sukop M C 2007 Phys. Rev. E 76 066701

    [18]

    Shan X W 2006 Phys. Rev. E 73 047701

    [19]

    Zhu M F, Stefanescu D M 2007 Acta Mater. 55 1741

    [20]

    Sukop M C, Thorne D T, Jr 2005 Lattice Boltzmann Modeling-An Introduction for Geoscientists and Engineers (2nd ed) (New York: Springer) p114

  • [1] 张士杰, 王颖明, 王琦, 李晨宇, 李日. 基于元胞自动机-格子玻尔兹曼模型的枝晶碰撞行为模拟. 物理学报, 2021, 70(23): 238101. doi: 10.7498/aps.70.20211292
    [2] 胡梦丹, 张庆宇, 孙东科, 朱鸣芳. 纳米结构超疏水表面冷凝现象的三维格子玻尔兹曼方法模拟. 物理学报, 2019, 68(3): 030501. doi: 10.7498/aps.68.20181665
    [3] 夏瑱超, 王伟丽, 罗盛宝, 魏炳波. 三元等原子比Fe33.3Cu33.3Sn33.3合金的快速凝固机理与室温组织磁性研究. 物理学报, 2016, 65(15): 158101. doi: 10.7498/aps.65.158101
    [4] 冯黛丽, 冯妍卉, 石珺. 介孔复合材料声子输运的格子玻尔兹曼模拟. 物理学报, 2016, 65(24): 244401. doi: 10.7498/aps.65.244401
    [5] 陈瑞, 许庆彦, 柳百成. 基于元胞自动机方法的定向凝固枝晶竞争生长数值模拟. 物理学报, 2014, 63(18): 188102. doi: 10.7498/aps.63.188102
    [6] 永贵, 黄海军, 许岩. 菱形网格的行人疏散元胞自动机模型. 物理学报, 2013, 62(1): 010506. doi: 10.7498/aps.62.010506
    [7] 梁善勇, 王江安, 宗思光, 吴荣华, 马治国, 王晓宇, 王乐东. 基于多重散射强度和偏振特征的舰船尾流气泡激光探测方法. 物理学报, 2013, 62(6): 060704. doi: 10.7498/aps.62.060704
    [8] 刘云龙, 张阿漫, 王诗平, 田昭丽. 基于边界元法的近平板圆孔气泡动力学行为研究. 物理学报, 2013, 62(14): 144703. doi: 10.7498/aps.62.144703
    [9] 孙东科, 项楠, 陈科, 倪中华. 格子玻尔兹曼方法模拟弯流道中粒子的惯性迁移行为. 物理学报, 2013, 62(2): 024703. doi: 10.7498/aps.62.024703
    [10] 张云鹏, 林鑫, 魏雷, 彭东剑, 王猛, 黄卫东. 界面能各向异性对定向凝固枝晶生长的影响. 物理学报, 2013, 62(17): 178105. doi: 10.7498/aps.62.178105
    [11] 陈海楠, 孙东科, 戴挺, 朱鸣芳. 凝固前沿和气泡相互作用的大密度比格子玻尔兹曼方法模拟. 物理学报, 2013, 62(12): 120502. doi: 10.7498/aps.62.120502
    [12] 张云鹏, 林鑫, 魏雷, 王猛, 彭东剑, 黄卫东. 用CA方法模拟界面能各向异性对胞晶生长形态的影响 . 物理学报, 2012, 61(22): 228106. doi: 10.7498/aps.61.228106
    [13] 刘云龙, 张阿漫, 王诗平, 田昭丽. 基于边界元法的气泡同波浪相互作用研究 . 物理学报, 2012, 61(22): 224702. doi: 10.7498/aps.61.224702
    [14] 魏雷, 林鑫, 王猛, 黄卫东. 基于MeshTV界面重构算法的二元合金自由枝晶生长元胞自动机模型. 物理学报, 2012, 61(9): 098104. doi: 10.7498/aps.61.098104
    [15] 周丰茂, 孙东科, 朱鸣芳. 偏晶合金液-液相分离的格子玻尔兹曼方法模拟. 物理学报, 2010, 59(5): 3394-3401. doi: 10.7498/aps.59.3394
    [16] 孙东科, 朱鸣芳, 杨朝蓉, 潘诗琰, 戴挺. 强制对流和自然对流作用下枝晶生长的数值模拟. 物理学报, 2009, 58(13): 285-S291. doi: 10.7498/aps.58.285
    [17] 黄锋, 邸洪双, 王广山. 用元胞自动机方法模拟镁合金薄带双辊铸轧过程凝固组织. 物理学报, 2009, 58(13): 313-S318. doi: 10.7498/aps.58.313
    [18] 单博炜, 林鑫, 魏雷, 黄卫东. 纯物质枝晶凝固的元胞自动机模型. 物理学报, 2009, 58(2): 1132-1138. doi: 10.7498/aps.58.1132
    [19] 周文平, 万松明, 张 霞, 张庆礼, 孙敦陆, 仇怀利, 尤静林, 殷绍唐. PbMoO4晶体生长基元和生长习性的高温拉曼光谱研究. 物理学报, 2008, 57(11): 7305-7309. doi: 10.7498/aps.57.7305
    [20] 赵代平, 荆 涛, 柳百成. 相场方法模拟铝合金三维枝晶生长. 物理学报, 2003, 52(7): 1737-1742. doi: 10.7498/aps.52.1737
计量
  • 文章访问数:  5730
  • PDF下载量:  993
  • 被引次数: 0
出版历程
  • 收稿日期:  2011-10-27
  • 修回日期:  2011-12-30
  • 刊出日期:  2012-08-05

枝晶生长和气泡形成的数值模拟

  • 1. 东南大学江苏省先进金属材料高技术研究重点实验室, 南京 211189
    基金项目: 国家自然科学基金(批准号: 50971042)资助的课题.

摘要: 本文建立了二维的格子玻尔兹曼方法-元胞自动机(lattice Boltzmann method-cellular automaton, LBM-CA)耦合模型, 对凝固过程中枝晶生长和气泡形成进行模拟研究. 本模型采用CA方法模拟枝晶的生长, 根据界面溶质平衡法计算枝晶生长的驱动力. 采用基于Shan-Chen多相流的LBM模拟气泡在液相中的生长和运动. 在LBM-CA的耦合模型中包含了固-液-气三相之间的相互作用. 应用Laplace定理和模拟气-液-固三相之间的润湿现象对模型进行了验证. 应用所建立的LBM-CA耦合模型模拟研究了气-液相互作用系数对单气泡生长的影响. 发现单气泡的生长速度和平衡半径随气-液相互作用系数的增大而增大. 定向凝固过程中枝晶和气泡生长的模拟结果再现了枝晶的择优生长、 气泡的优先形核位置、气泡的长大、合并、在枝晶间受挤变形以及在液相通道中的运动等物理现象, 与实验结果符合良好. 此外, 初始气体含量越高, 凝固结束时气泡的体积分数也相对较高. 本模型的模拟结果可以揭示在凝固过程中气泡形核、 生长和运动演化以及与枝晶生长相互作用的物理机理.

English Abstract

参考文献 (20)

目录

    /

    返回文章
    返回