Study on Interparticle Interaction Force Model to Correct Saturation Density of Real Cryogenic Fluid for LBM Simulation

Abstract

Cryogenic liquefaction energy storage is an important form of storage for sustainable energy liquid hydrogen and other gases. The weighting parameter A in the parameter-adjusted two-phase LBM model is important for the deviation of simulation results. The aim of this paper is to discover the appropriate parameter to eliminate the deviation, and to solve the problem of large deviation between the theoretical solution and the simulated value that is caused by using different equations of state in LBM simulation. The modified PT equation of state, which is suitable for cryogenic fluids, is combined with the parameter-adjustable two-phase model to simulate the saturation density at different temperatures. Four typical cryogenic fluids—nitrogen, hydrogen, oxygen, and helium—are exploratively simulated to find the suitable parameters to eliminate errors by analyzing the results with theoretical solutions. This is an efficient solution to the deviation between the simulated value and the theoretical solutions, which is caused by the different equation of state in LBM. The optimal A-value of the model based on the PT equation of state was obtained as −0.21, while droplets and bubbles were set into the calculation region, and an inverse relationship between the interface density gradient and temperature was analyzed. The analysis and comparison of the simulation results under the optimal value and the experimental values have laid an important foundation for the phase change simulation of the real cryogenic fluids at the mesoscopic scale.