A consistent, volume preserving, and adaptive mesh refinement-based framework for modeling non-isothermal gas-liquid-solid flows with phase change

Ramakrishnan Thirumalaisamy,Amneet Pal Singh Bhalla

2024-07-08

Abstract:This work expands on our recently introduced low Mach enthalpy method [1] for simulating the melting and solidification of a phase change material (PCM) alongside (or without) an ambient gas phase. The method captures PCM's volume change (shrinkage or expansion) by accounting for density change-induced flows. We present several improvements to the original work. First, we introduce consistent time integration schemes for the mass, momentum, and enthalpy equations, which enhance the method stability. Demonstrating the effectiveness of this scheme, we show that a system free of external forces and heat sources can conserve its initial mass, momentum, enthalpy, and phase composition. This allows the system to transition from a non-isothermal, non-equilibrium, phase-changing state to an isothermal, equilibrium state without exhibiting unrealistic behavior. Furthermore, we show that the low Mach enthalpy method accurately simulates thermocapillary flows without introducing spurious phase changes. We propose an analytical model to validate advanced CFD codes for simulating metal manufacturing processes like welding and 3D printing. These processes involve a heat source melting metal or alloy in an inert gas environment. Traditionally, validation relied on manipulating material properties to match complex experiments. Our model uses the Stefan problem with a density jump to provide a straightforward method for validating multiphysics simulations involving heat sources and phase changes in three-phase flows.Lastly, we demonstrate the practical utility of the method in modeling porosity defects (gas bubble trapping) during metal solidification. A field extension technique is used to accurately apply surface tension forces in a three-phase flow situation. This is where part of the bubble surface is trapped within the (moving) solidification front.

Fluid Dynamics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to accurately simulate non - isothermal gas - liquid - solid three - phase flows involving phase transitions (such as melting and solidification) in the metal manufacturing process. Specifically, the paper proposes an improved framework based on the low Mach enthalpy method, which has the following characteristics:

1. **Accurate capture of volume change**: In traditional methods, when simulating the melting and solidification processes of phase - change materials (PCM), the volume change (contraction or expansion) caused by density change is often ignored. The method proposed in this paper can accurately capture the volume change of PCM by considering the flow caused by density change.

2. **Consistent time - integration scheme**: In order to improve the stability of the method, this paper introduces a consistent time - integration scheme to solve the mass, momentum and enthalpy equations. This enables the system to maintain the initial mass, momentum, enthalpy and phase composition in the absence of external forces and heat sources, so that it can smoothly transition from a non - isothermal, non - equilibrium, phase - transition state to an isothermal, equilibrium state without unrealistic behavior.

3. **Adaptive mesh refinement (AMR)**: In order to reduce the computational cost, especially in three - dimensional simulations, this paper adopts the adaptive mesh refinement technique. By comparing two cell marking / unmarking criteria, it is found that the gradient - based method is more effective. This method ensures that the moving thin mushy region is always captured by the fine mesh, even if it temporarily falls within the sub - grid level.

4. **Accurate simulation of thermocapillary flow**: The method in this paper can accurately simulate the thermocapillary flow caused by temperature - dependent surface tension without introducing false phase - transition phenomena. This is particularly important for simulations in the additive manufacturing (such as 3D printing) process.

5. **Multi - physics field simulation verification**: This paper proposes an analytical model for verifying multi - physics field simulations involving heat sources (such as lasers) and phase - transition phenomena. This model utilizes the analytical solution of the Stefan problem with density jumps, providing a simple method to verify three - phase flow simulations involving heat sources and phase - transition phenomena.

6. **Modeling of pore defects**: This paper also shows the practical application of this method in simulating pore defects (such as gas bubble capture) in the metal solidification process. A partial differential equation (PDE) extension technique is adopted, which can accurately apply surface tension forces in the case of three - phase flow. Part of the gas bubble surface is captured by the (moving) solidification front, while the rest resists deformation due to the action of surface tension.

In summary, the main objective of this paper is to provide a more stable, efficient framework that can handle complex physical phenomena by improving the low Mach enthalpy method, in order to solve the multi - physics field simulation problems involving phase transitions in the metal manufacturing process.