A low Mach enthalpy method to model non-isothermal gas-liquid-solid flows with melting and solidification

Ramakrishnan Thirumalaisamy,Amneet Pal Singh Bhalla

DOI: https://doi.org/10.1016/j.ijmultiphaseflow.2023.104605

2023-07-28

Abstract:Modeling phase change problems numerically is vital for understanding many natural (e.g., ice formation, steam generation) and engineering processes (e.g., casting, welding, additive manufacturing). Almost all phase change materials (PCMs) exhibit density/volume changes during melting, solidification, boiling, or condensation, causing additional fluid flow during this transition. Most numerical works consider only two phase flows (either solid-liquid or liquid-gas) for modeling phase change phenomena and some also neglect volume/density change of PCMs in the models. This paper presents a novel low Mach enthalpy method for simulating solidification and melting problems with variable thermophysical properties, including density. Additionally, this formulation allows coupling a solid-liquid PCM with a gas phase in order to simulate the free surface dynamics of PCMs undergoing melting and solidification. We revisit the two-phase Stefan problem involving a density jump between two material phases. We propose a possible means to include the kinetic energy jump in the Stefan condition while still allowing for an analytical solution. The new low Mach enthalpy method is validated against analytical solutions for a PCM undergoing a large density change during its phase transition. Additionally, a few simple sanity checks are proposed to benchmark computational fluid dynamics (CFD) algorithms that aim to capture the volume change effects of PCMs.

Fluid Dynamics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the deficiencies of existing numerical simulation methods in dealing with volume and density changes during the process of phase - change materials (PCMs). Specifically:

1. **Multiphase - flow problems**: Most of the existing numerical works only consider two - phase flow (such as solid - liquid or liquid - gas), and some also ignore the volume/density changes of phase - change materials during the phase - change process. This leads to an incomplete description of the actual phase - change process, especially when it comes to melting, solidification, boiling or condensation.

2. **Flow caused by density changes**: Many studies ignore the additional fluid flow caused by density changes. For example, in the process of metal additive manufacturing (such as selective laser melting, SLM), the gas plume in the molten pool is entrapped, resulting in porosity and keyhole defects in the finished product. These phenomena require more accurate simulations to predict and control.

3. **High - density - ratio flow**: In some applications, such as the selective laser melting of metal powders, the density ratio between gas and metal is as high as \(10^4\), which poses a challenge to the stability of numerical methods.

4. **Multiphysics coupling**: Existing methods usually cannot handle multiple phase - change phenomena such as evaporation, condensation, melting and solidification simultaneously. This limits their application in complex manufacturing processes, such as metal additive manufacturing.

To solve these problems, this paper proposes a new low Mach enthalpy method, which is able to:

- Simulate solid - liquid phase - change and melting problems with variable thermophysical properties (including density).
- Couple solid - liquid phase - change materials with the gas phase to simulate the dynamics of the free surface of phase - change materials during melting and solidification.
- Handle high - density - ratio flow and ensure the stability and accuracy of numerical methods.
- Verify the effectiveness of the new method through analytical solutions and verify the ability of the computational fluid dynamics (CFD) algorithm to capture the volume - change effects of phase - change materials through simple benchmark tests.

### Formula display

1. **Mass conservation equation**:
\[
\frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{u}) = 0
\]

2. **Momentum conservation equation**:
\[
\frac{\partial(\rho\mathbf{u})}{\partial t}+\nabla\cdot(\rho\mathbf{u}\otimes\mathbf{u}+p\mathbf{I})-\sigma C\delta(\mathbf{x}-\mathbf{s})\mathbf{n} = 0
\]
where \(\sigma\) is the surface tension coefficient between the two phases, \(C\) is the average local curvature of the interface, \(\mathbf{s}\) is the position of the interface, \(\delta\) is the Dirac delta distribution, and \(\mathbf{n}\) is the outer unit normal vector of the interface.

3. **Energy conservation equation**:
\[
\frac{\partial}{\partial t}\left(\rho\left(e+\frac{1}{2}|\mathbf{u}|^{2}\right)\right)+\nabla\cdot\left(\rho\left(e+\frac{1}{2}|\mathbf{u}|^{2}\right)\mathbf{u}+\mathbf{q}+p\mathbf{u}\right)-\sigma C\delta(\mathbf{x}-\mathbf{s})\mathbf{u}\cdot\mathbf{n} = 0
\]
where \(e\) is the internal energy and \(\mathbf{q}=-\kappa\nabla T\) is the conductive heat flux.

4. **Stefan condition**:
\[
\rho_S\left((C