基于GPU并行的烧结-干道复合吸液芯渗透率格子Boltzmann模拟研究

Lattice Boltzmann Simulation on Permeability of Sintered-arterial Composite Wick Using GPU Parallel

  • 摘要: 针对先进热管堆对高性能吸液芯的迫切需求,基于GPU并行对850 °C液态钠在吸液芯中的渗透率开展了格子Boltzmann模拟。利用四参数随机生长方法生成烧结多孔介质,分析表明其渗透率随孔隙率(ε)的增大而增大,随初始种子分布(S_\mathrmd)的增大而减小。研究发现,烧结-干道复合吸液芯的渗透率随干道直径增大显著提高,流动机制会逐渐从多孔介质渗流主导转变为干道管流主导。基于流动独立性假设建立了渗透率预测模型,预测值与模拟结果吻合良好,最大相对误差为11.4%。同时发现由于干道与烧结介质之间的流体交换,模拟值会系统性高于预测值。这体现为等效干道直径大于实际直径,该差值随ε的增大而增加,随S_\mathrmd的增大而减小。本文研究为复合吸液芯的优化设计提供了理论见解与模型支持。

     

    Abstract: Permeability is the core transport parameter of heat pipe wicks, and its accurate prediction is essential for the design of high-temperature heat pipes used in advanced heat pipe cooled reactors. Sintered-arterial composite wicks balance high capillary force and high permeability, but their cross-scale microstructure makes traditional computational fluid dynamics methods difficult to perform fine simulations. The lattice Boltzmann method (LBM) has unique advantages in this field. It can accurately handle complex porous media geometries, has natural parallelism suitable for large-scale computing, and provides clear physical images. This study aims to systematically investigate the permeability characteristics of a sintered-arterial composite wick and establish a predictive model, thereby providing crucial insights for the design of high-performance wicks. In this paper, a GPU-accelerated LBM solver based on the D3Q19 discrete velocity model and half-step bounce-back boundary condition was developed. Large-scale simulations of liquid sodium flow at 850 °C were performed, achieving a maximum grid scale of over 100 million on a single NVIDIA RTX A6000 graphics card. The sintered porous matrix was reconstructed using the quartet structure generation set (QSGS) method. Firstly, a comprehensive parametric study on the pure sintered wick was conducted to analyze the effects of porosity (ε), initial seed distribution (S_\mathrmd), and reference growth probability (G_\mathrmref). Subsequently, circular channel of varying diameters (D) were introduced into the sintered matrix to form the composite wick. The flow dynamics were analyzed through velocity contours and streamlines. Finally, a simplified permeability prediction model for the composite wick was developed based on the superposition of independent Poiseuille flow in the arterial and Darcy flow in the sintered matrix. For the pure sintered wick, permeability increases significantly with ε but decreases with a higher S_\mathrmd, while remaining largely insensitive to G_\mathrmref. In the composite wick, the overall permeability increases with the channel diameter, transitioning through three distinct regimes: sintered-matrix-dominated, mixed-flow, and arterial-dominated flow. The flow in the arterial exhibits a parabolic profile, similar to laminar pipe flow, while fluid exchange occurs at the arterial-matrix interface. The proposed prediction model agrees well with the LBM simulation results, with a maximum relative deviation of 11.4%. However, the model systematically underestimates the permeability because it neglects the interfacial fluid exchange, which is quantified by the positive difference between the equivalent arterial diameter and the actual diameter. The magnitude of this deviation increases with ε but decreases with S_\mathrmd, showing negligible dependence on G_\mathrmref. This work demonstrates the superior performance of LBM in cross-scale porous media flow simulation, clarifies the key factors governing composite wick permeability, and provides a practical predictive tool for the optimized design of high-performance wicks in heat pipe cooled reactors

     

/

返回文章
返回