Exponential Time Differencing Schemes for the Peng-Robinson Equation of State with Preservation of Maximum Bound Principle

The Peng-Robison equation of state,one of the most extensively applied equations of state in the petroleum industry and chemical engineering,has an excel-lent appearance in predicting the thermodynamic properties of a wide variety of ma-terials.It has been a great challenge on how to design numerical schemes with preser-vation of mass conservation and energy dissipation law.Based on the exponential time difference combined with the stabilizing technique and added Lagrange multi-plier enforcing the mass conservation,we develop the efficientfirst-and second-order numerical schemes with preservation of maximum bound principle(MBP)to solve the single-component two-phase diffuse interface model with Peng-Robison equation of state.Convergence analyses as well as energy stability are also proven.Several two-dimensional and three-dimensional experiments are performed to verify these theo-retical results.