A novel numerical simulation approach for cryogenic CO2 frosting in binary mixture gas by integrating desublimation and gas-solid phase equilibrium models

Cryogenic CO₂ separation from natural gas is a promising technology due to its environmental benefits and smaller spatial requirement. Nonetheless, the dynamics of CO₂ desublimation and subsequent frost layer formation under cryogenic conditions for natural gas pre-treatment remain largely unexplored. Addressing this gap, this paper presents a novel numerical simulation approach by integrating desublimation phase change model and vapor-solid phase equilibrium model to investigate CO₂ frost layer formation on the cold wall, revealing that both the cold wall temperature and the initial CO₂ concentration significantly influence the thickness of the CO₂ frost layer. Moreover, the study finds that the density of frost layer intensifies over time, with the highest density observed near the cold surface. Furthermore, increasing the cold wall temperature and the flow rate not only increases the average density of the CO₂ frost layer but also affects the outlet CO₂ concentration. To adhere to specific outlet CO₂ concentration targets and desublimation duration, it is critical to adjust the cold wall temperature and the inlet flow velocity dynamically. Therefore, cryogenic CO₂ removal methods, by reducing the spatial requirements of CO₂ pre-treatment units, offer a practical and efficient solution for small-scale natural gas liquefaction plants. This approach not only facilitates operational efficiency but also contributes to the broader effort of mitigating greenhouse gas emissions, aligning with environmental sustainability objectives.