Diffuse-interface modeling and simulation of the freezing of binary fluids with the Marangoni effect

TL;DR

This paper develops a comprehensive diffuse-interface model using a phase-field and lattice Boltzmann method to simulate complex multiphase freezing processes involving solute transport, volume change, and Marangoni effects, validated against experiments.

Contribution

The paper introduces a novel diffuse-interface model that integrates phase-field, enthalpy, and solute transport methods with a lattice Boltzmann solver for multiphase freezing with Marangoni effects.

Findings

Model accurately captures freezing dynamics with impurities.

Numerical results agree well with experimental data.

The approach preserves mass and volume during phase change.

Abstract

This paper proposes a diffuse-interface model for simulating gas-liquid-solid multiphase flows involving solid-liquid phase change, solute transport, and the Marangoni effect. In this model, a phase-field method is employed to capture the evolution of fluid-fluid interfaces, while an enthalpy-based approach is used to describe the temperature field and implicitly track the solid-liquid interface. Solute transport is modeled using a constrained scalar-transport model combined with a pseudo-potential concentration approach. The proposed diffuse-interface model satisfies the reduction consistency, and can degenerate to the conservative phase-field method for incompressible two-phase flow and the classical enthalpy method for binary material solidification in an appropriate way. Furthermore, the model not only can preserve the mass conservation, but also can capture the volume change induced by phase change. To solve the diffuse-interface model, a lattice Boltzmann (LB) method is then developed, and the numerical tests demonstrate that the method has a good performance in the study of the freezing process coupled with Marangoni flow, phase-change-induced volume change, and solute transport. Finally, the model is applied to investigate the freezing dynamics of a system containing an insoluble impurity, revealing the complex interaction between the advancing freezing front and the impurity. It is found that the numerical results are in good agreement with experimental data.