刘江， 李元基， 王帅
哈尔滨工业大学 a.能源科学与工程学院，黑龙江 哈尔滨 150001； b.郑州研究院，河南 郑州 450001
引用格式：
刘江， 李元基， 王帅. 铝熔体在多孔预制体渗流行为的数值模拟［J］. 中国粉体技术， 2026， 32（5）： 1-12.
Citation:Liu Jiang， Li Yuanji， Wang Shuai. Numerical simulation of infiltration behavior of aluminum melt in porous preforms［J］. China Powder Science and Technology， 2026， 32（5）： 1-12.
DOI：10.13732/j.issn.1008-5548.2026.05.001
收稿日期： 2026-05-16， 修回日期： 2026-08-06， 上线日期： 2026-08-26。
基金项目： 国家自然科学基金项目，编号：52076060。
第一作者简介： 刘江（2001—），男，硕士生，研究方向为孔隙尺度模拟。E-mail：2396684454@qq.com。
通信作者简介： 王帅（1985—），男，教授，博士，博士生导师，黑龙江省自然科学优秀青年基金获得者，研究方向为多相流动与传热。E-mail：shuaiwang@hit.edu.cn。
摘要： 【目的】 揭示燃料颗粒多孔预制体内铝熔体的渗透演化规律，实现压力浸渗制备铝基弥散燃料元件的工艺调控。【方法】 构建燃料颗粒多孔介质模型，采用相场法，对铝熔体在多孔预制体中的浸渗过程进行孔隙尺度模拟，系统分析浸渗压力、熔体黏度、接触角等参数对铝熔体浸渗行为的影响规律。【结果】 入口浸渗压力越低，熔体流动优势通道现象越显著，浸渗饱和度越低； 低黏度熔体能在极短时间内完成对多孔介质孔隙的填充，且浸渗饱和度较高，但也会因缺少黏性剪切的稳定流动作用，使孔隙空间内更多气体被困于孔道内部难以排出； 毛细驱动力与入口浸渗压力呈现协同效应。【结论】 熔体浸渗过程中的钉扎效应与弯液面几何形貌演化是导致孔隙内形成孤立气泡的主要因素。
关键词： 铝熔体； 浸渗； 饱和度； 孔隙尺度
Abstract
Objective Melt infiltration is a favorable manufacturing route for fabricating aluminum matrix dispersed fuel elements. Using preform skeletons and melt infiltration densification， this method effectively overcomes the limitations of conventional manufacturing routes， such as low fuel volume fraction， uneven fuel distribution， and poor densification. Nevertheless， the transient filling mechanism of the infiltration front in porous media and the quantitative correlation between infiltration pressure and saturation remain insufficiently understood. Therefore， clarifying the intrinsic relationship between aluminum melt infiltration behavior and its key influencing factors， and revealing the effects of infiltration processes and melt physical properties on infiltration characteristics are of great significance for developing advanced fabrication technologies for high-performance aluminum matrix dispersed fuels.
Methods In this study， the phase-field method was adopted to characterize two-phase flow behaviors. The Cahn-Hilliard equation was applied to capture the evolution of the interface position， where the chemical potential was defined as a function of mixing energy density and interface thickness control parameters. A porous medium model consisting of fuel particles was established， and pore-scale simulations of aluminum melt infiltration were carried out. The feasibility of the proposed model was verified against experimental data.
Results and Discussion Under external pressure， the melt preferentially advanced toward both sides. After the displacement fronts on both sides reached a certain depth， the melt gradually infiltrated and filled the intermediate region. When the melt flowed into dead-end pores， the liquid meniscus ruptured， and newly formed liquid-gas interfaces encapsulated residual gas within pores to form discrete isolated bubbles. When two opposing menisci converged and coalesced， the convex liquid surface could not fully spread， and residual gas was trapped by pinned liquid interfaces to form permanent pore void defects. For porous architectures featuring large cavities confined by narrow pore throats， interfacial pinning readily occurred， which suppressed effective gas evacuation from internal cavities and ultimately generated isolated bubbles. In the late infiltration stage， pressure fluctuations exhibited a lower frequency. A plateau appeared on the flow velocity curve in the late stage， indicating that the melt was trapped in narrow pores and flowed steadily along dominant flow channels at a nearly constant flow velocity. Sufficient driving pressure in the early infiltration stage yielded high infiltration saturation. Saturation fluctuated in the middle and late stages， mainly caused by Haines jumps induced by driving pressure attenuation. In the final stage， the discrepancy in infiltration saturation between samples under different inlet pressures became pronounced， and higher inlet driving pressure consistently resulted in greater infiltration saturation. Low-viscosity melts could fully occupy pore voids within an extremely short period， yielding high infiltration saturation. However， the absence of flow stabilization from viscous shearing readily triggered disordered channeling during melt infiltration. This trapped a greater amount of gas inside pore channels that could not easily escape， thereby deteriorating infiltration densification. When aluminum melt exhibited poor wettability with pore wall substrates， the liquid front became pinned at narrow pore throats， and a large external pressure was required to overcome pinning. In contrast， systems with favorable melt wettability against pore walls required much lower external pressure for depinning， and the depinning process proceeded spontaneously， facilitating uniform full pore infiltration. Poor wettability resulted in larger differences in displacement front depth， which continuously expanded as the overall infiltration depth increased. For strongly wetting interfaces， the differences in displacement front depth increased in the early infiltration stage. In the later stage， the average throat radius of unfilled pores decreased， the capillary driving force increased， and the differences in displacement front depth declined.
Conclusion In this study， pore-scale numerical simulations are conducted to characterize aluminum melt infiltration within porous preforms using a constructed packed-particle porous medium model. The results show that isolated bubbles during melt infiltration originate from meniscus evolution induced by pore geometric configurations and interface pinning effects. Under low inlet infiltration pressure， the dominant seepage channel is preferentially formed， which reduces overall infiltration saturation and increases the occurrence probability of Haines jumps. Low-viscosity melts can rapidly fill porous voids and reach high saturation. Nevertheless， the lack of flow stabilization from viscous shear induces chaotic channeling. A large amount of residual gas remains trapped inside pore channels and cannot escape， impairing the uniformity of pore filling. When capillary force acts as flow resistance in melt-solid systems， the liquid front is prone to pinning at narrow pore throats， and high external pressure is required to achieve depinning， resulting in non-uniform infiltration. At strongly wetting interfaces， capillary force cooperates with external pressure to promote infiltration， and less external force is needed for depinning. Smaller contact angles correspond to improved infiltration densification. However， isolated gas bubbles easily form when melt flows through large cavities surrounded by tiny pore throats.
Keywords： aluminum melt； infiltration； saturation； pore scale
参考文献（References）
［1］Ginting A B， Prasetyo R S E B， Boybul， et al. Post-irradiation examination of silicide fuel （U3Si2-Al） with 4.8 gU/cm3 for RSG-GAS［J］. Progress in Nuclear Energy， 2025， 183： 105673.
［2］Sihotang J C， Ajiriyanto M K， Putri A S D， et al. Study of fission gas bubbles and interaction layer on irradiated U3Si2-Al density of 4.8 gU/cm3［J］. Urania： Jurnal Ilmiah Daur Bahan Bakar Nuklir， 2022， 28（2）： 69-78.
［3］Holmström S， Baumeister B， Wight J. Qualification of advanced LEU fuels for high-power research reactor conversion designs［J］. EPJ Nuclear Sciences Technologies， 2025， 11： 18.
［4］Song Y， Lv M， Xie H. Effects of the curvatures and creep on the thermo-mechanical behaviors in curved U-Si/Al dispersion fuel plates under irradiation condition［J］. Annals of Nuclear Energy， 2025， 216： 111261.
［5］Li X， Xue W， Yin C， et al. Numerical simulation of the infiltration behavior of molten metal in a ceramic porous preform［J］.Physics of Fluids， 2025， 37： 033385.
［6］Zhao L， Xiang Y， Ma W. Numerical simulation of melt penetration in debris beds using MPS method［J］. Progress in Nuclear Energy， 2024， 167： 104982.
［7］Chen L， Fang D， Villanueva W， et al. Microscopic insights into the effect of surface wettability on melt infiltration in porous media： Experimental and simulation studies［J］. Physics of Fluids， 2026， 38： 012017.
［8］Xiao L， Zhu G， Zhang L， et al. Effects of pore-size disorder and wettability on forced imbibition in porous media［J］. Journal of Petroleum Science and Engineering， 2021， 201： 108485.
［9］Yun K， Zhou J， Wu H， et al. Modeling and experimental analysis of infiltration in graphene preforms using a level-set appr-oach［J］. Composites Part B： Engineering， 2026， 312： 113347.
［10］Song R， Tang Y， Wang Y， et al. Pore-scale numerical simulation of CO2-oil two-phase flow： a multiple-parameter analysis based on phase-field method［J］. Energies， 2023， 16（1）： 82.
［11］Wang Z， Guo C， Liu N， et al. Effect of wettability and permeability on pore-scale of CH4-water two-phase displacement behavior in the phase field model［J］. Applied Sciences， 2024， 14（15）： 6815.
［12］杨柳， 赵子恒， 张记刚， 等. 基于相场理论的含裂隙多孔介质自发渗吸规律［J］. 石油科学通报， 2025， 10（2）： 269-282.
Yang Liu， Zhao Ziheng， Zhang Jigang， et al. Spontaneous imbibition laws of fracture-containing porous media based on phase field theory［J］. Petroleum Science Bulletin， 2025， 10（2）： 269-282.
［13］Zeng C， Zhang Y， Lu H， et al. Simulation of CO2-water two-phase fluid displacement characteristics based on the phase field method［J］. Deep Underground Science and Engineering， 2025， 4（4）： 725-738.
［14］Smith C A， Keiser D D， Miller B D， et al. Microstructural dependence on fuel matrix composition in irradiated U-Mo dispersion fuels［J］. Journal of Nuclear Materials， 2021， 550： 152943.
［15］Badalassi V E， Ceniceros H D， Banerjee S. Computation of multiphase systems with phase field models［J］. Journal of Computational Physics， 2003， 190： 371-397.
［16］赵海军， 魏爱华， 彭昊， 等. 非均质多孔介质超临界CO2非混相驱替特征及机理研究［J］. 工程地质学报， 2024， 32（5）： 1825-1837.
Zhao Haijun， Wei Aihua， Peng Hao， et al. Study on characteristics and mechanism of immiscible dis-placement of supercritical CO2 in heterogeneous porous media［J］. Journal of Engineering Geology， 2024， 32（5）： 1825-1837.
［17］Rokhforouz M R， Akhlaghi Amiri H A. Effects of grain size and shape distribution on pore-scale numerical simulation of two-phase flow in a heterogeneous porous medium［J］. Advances in Water Resources， 2019， 124： 84-95.
［18］Akhlaghi Amiri H A， Hamouda A A. Pore-scale modeling of non-isothermal two phase flow in 2D porous media： Influences of viscosity， capillarity， wettability and heterogeneity［J］. International Journal of Multiphase Flow， 2014， 61： 14-27.
［19］Yi T， Yang G， Wang B， et al. Dynamics of a gas bubble penetrating through porous media. Physics of Fluids， 2022， 34： 012103.
［20］Yue P， Zhou C， Feng J J， et al. Phase-field simulations of interfacial dynamics in viscoelastic fluids using finite elements with adaptive meshing［J］. Journal of Computational Physics， 2006， 219： 47-67.
［21］Shi H， Zhu Q， Chen Z， et al. Pore-scale modeling of water–gas flow in heterogeneous porous media［J］. Physics of Fluids， 2023， 35： 072114.
［22］Yin X， Zarikos I， Karadimitriou N K， et al. Direct simulations of two-phase flow experiments of different geometry complexities using Volume-of-Fluid （VOF） method［J］. Chemical Engineering Science， 2019， 195： 820-827.
［23］Wang K， Matsushita S， Sin S， et al. Characteristics of pore-scale events and their impact on transport in porous media［J］. Physics of Fluids， 2025， 37： 032118.
［24］Yan L， Müller J C， van Noorden T L， et al. Wettability-driven pore-filling instabilities： Microfluidic and numerical insights［J］. Journal of Colloid and Interface Science， 2025， 696： 137884.
［25］刘先珊， 孙梦， 曹伊婷， 等. 基于相场法孔隙介质两相流体流动特性研究［J］. 岩石力学与工程学报， 2026， 45（5）： 1461-1476.
Liu Xianshan， Sun Meng， Cao Yiting， et al. Two-phase flow mechanism at pore scale using the phase field method［J］. Chinese Journal of Rock Mechanics and Engineering， 2026， 45（5）： 1461-1476.
［26］ 李培生， 王昭太， 赵万东， 等. 液态法制备复合材料的多组分伪势LBM模拟［J］. 华南理工大学学报（自然科学版）， 2018， 46（7）： 7-15.
Li Peisheng， Wang Zhaotai， Zhao Wandong， et al. Simulation of liquid infiltration process for reparation of composite materials by multi-component pseudopotential LBM［J］. Journal of South China University of Technology （Natural Science Edition）， 2018， 46（7）： 7-15.
