Solving the UVN-flash problem in TVN-space

We investigate the phase equilibrium problem for multicomponent mixtures under specified internal energy (U), volume (V), and mole numbers (N1,N2, . . . ,Nn), commonly known as the UVN-flash problem. While conventional phase equilibrium calculations typically use pressure-temperature-mole number (PTN) specifications, the UVN formulation is essential for dynamic simulations of closed systems and energy balance computations. Existing approaches, including those based on iterative pressure-temperature updates and direct entropy maximization, suffer from computational inefficiencies due to nested iterations and reliance on inner Newton solvers. In this work, we present a novel reformulation of the UVN-flash problem as a direct entropy maximization problem that eliminates the need for inner Newton iterations, addressing key computational bottlenecks. We derive two new novel formulations: 1) a formulation based on entropy and internal energy and (2) an alternative formulation based on Helmholtz free energy. We begin with a stability analysis framework, followed by a reformulation of the UVN flash problem in natural variables. We then introduce our novel approach and discuss the numerical methods used, including gradient and Hessian computations. The proposed method is validated against benchmark cases, demonstrating improved efficiency and robustness.