OpenFOAM is a leading open-source CFD platform widely used for simulating complex multiphase and phase-change flows. It employs a mass-conserving Volume of Fluid (VOF) method to capture fluid interfaces through scalar advection, enabling robust handling of large deformations and accurate interface reconstruction. Building on this foundation, the present work focuses on developing a multicomponent phase-change solver in OpenFOAM capable of handling the coupled transport of mass, momentum, energy, and species across interfaces. The solver incorporates interfacial heat and mass transfer models, thermodynamic equilibrium relations, and variable thermophysical properties for multiple components. Advanced numerical techniques are implemented to maintain stability and resolve sharp interfacial gradients. The framework is verified against benchmark cases and designed to remain fully compatible with OpenFOAM’s modular solver architecture. This integrated approach enhances OpenFOAM’s capability for high-fidelity simulations of realistic phase-change processes relevant to advanced thermal management systems.

​

​

​

​

​

​

​

​

​

​

​

​

​

 

​Efficient heat transfer underpins emerging technologies such as electronics cooling, microfluidics, and advanced phase change systems, where interfacial dynamics critically affect thermal performance. Modeling these processes demands accurate representation of interfacial motion, deformation, and phase transformation under coupled thermal and flow conditions. To address these challenges, this research employs the OpenFOAM computational framework for modeling multiphase and phase-change flows. The work utilizes and customizes the advanced geometric Volume of Fluid (VOF) formulations available in the TwoPhaseFlow framework, which employ interface reconstruction techniques to maintain sharp phase boundaries and accurately capture surface tension effects. These capabilities enable high fidelity prediction of droplet impact, film evolution, and evaporation -condensation behaviour under diverse thermal and flow conditions. Building upon this framework, electrohydrodynamic (EHD) effects are investigated to understand how electric fields can influence interfacial motion and augment heat transfer. The interplay between EHD forces, Marangoni stresses, and natural convection reveals mechanisms for active control of fluid flow and thermal transport in both single phase and phase change regimes. Such coupled analyses deepen the understanding of field driven interfacial phenomena and contribute to the design of next generation cooling and energy conversion systems.