Simulations of Ga melting based on multiple-relaxation time lattice Boltzmann method performed with CUDA in Python

A new solver, via the enthalpy multiple-relaxation lattice Boltzmann method, is developed to simulate the Ga melting (considering Ga as a phase change material) for different settings. At first, the phase change simulation of a simple bar is performed, this case is implemented to validate the heat transfer in our model via the analytical solution. Second, the solid–liquid phase change simulation with convection driven by gravity of Ga immersed in a 2D non-Darcy heterogeneous porous media, obtained from an image, is provided. Here, we present a novel process where Kozeny law is used globally for the homogeneous porosity and locally for the heterogeneous porosity. The second case is validated by experimental data already published in the literature. The fact that our solver is enforced in a single Nvidia GPU device with CUDA technology in Phyton 3.8 is a new feature introduced in this paper. Our results are given as a different approach to phase transitions embedded in porous media within an acceptable error margin from analytical and experimental results.