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

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.

References in corpus (10)

• A robust incompressible Navier-Stokes solver for high density ratio multiphase flows
• Numerical Study of Molten Metal Melt Pool Behaviour during Conduction-mode Laser Spot Melting
• Direct numerical simulation of variable surface tension flows using a Volume-of-Fluid method
• A DLM immersed boundary method based wave-structure interaction solver for high density ratio multiphase flows
• Simulating water-entry/exit problems using Eulerian-Lagrangian and fully-Eulerian fictitious domain methods within the open-source IBAMR library
• The inertial sea wave energy converter (ISWEC) technology: device-physics, multiphase modeling and simulations
• Breaching the capillary time-step constraint using a coupled VOF method with implicit surface tension
• A low Mach enthalpy method to model non-isothermal gas-liquid-solid flows with melting and solidification
• Preventing mass loss in the standard level set method: New insights from variational analyses
• An effective preconditioning strategy for volume penalized incompressible/low Mach multiphase flow solvers

Cited by in corpus (4)

• Consistent continuum equations and numerical benchmarks for a perturbation-based, variable-coefficient acoustofluidic solver
• The influence of volumetric shrinkage on the metal solidification process under localized energy deposition
• Simulating acoustically-actuated flows in complex microchannels using the volume penalization technique
• Analytical and numerical solutions to the three-phase Stefan problem with simultaneous occurrences of melting, solidification, boiling, and condensation phenomena