Energy landscape analysis for two-phase multi-component NVT flash systems by using ETD type high-index saddle dynamics☆
Introduction
Thanks to the rapid development in computing equipment and techniques, the past decades have borne remarkable progress in the more refined modeling and simulation of multi-phase fluid flows with multiple components. Studies on phase distributions and phase transitions have been incorporated into an increasing number of engineering applications to improve the production efficiency [1], ensure the operation safety, [2] and enhance the environmental protection [3]. As a concept defining the aggregated state of homogeneous substances with the same physical and chemical properties, phases are determined and separated by a clear physical interface, crossing which a sudden change in certain thermo properties can be captured [4]. Thus, the observations on the phase behaviors are focused on finding the change of thermodynamic properties, especially modeling the non-analytical free energy function representing phase transitions [5]. The vapor phase and liquid phase are the two major fluid phases existing in the multi-phase flow problems, and the description of the phase transitions between the two phases is usually carried out by flash calculation [6]. By calculating the compositional distributions of each phase at equilibrium, the flash calculation provides the fundamental thermodynamic knowledge for the multi-phase flow modeling and simulation, including the number of phases and the initial phase distributions with physical meanings [7]. In addition, phase equilibrium investigation can also used for independent calculation, analysis, design and optimization of engineering processes and facilities, including reducing unnecessary substances (for example, CO2 and H2O) to purify the production of hydrogen, removing acidic gases (for example, CO2 and H2S) to protect pipelines and equipment from corrosion, determining the amount of inhibitors (for example, methanol and glycol) to avoid the formation of gas hydrates, controlling asphaltene precipitation to enhance flow assurance, and many others [8]. Based on the different constraints on phase equilibrium, flash calculations are usually divided into two categories: NPT flash with constant chemical compositions, pressure, and temperature, and NVT flash with constant chemical compositions, volume, and temperature. Comparing with the NPT flash framework, NVT flash model has well-posed formulation and exhibits a unique pressure-volume relation so that the root selection procedure in PT flash framework can be avoided [9], [10]. At the same time, considering the easier and more intuitive treatment on capillarity, which is the critical mechanism in a number of popular engineering applications including shale gas recovery [11] and microfluidic fuel cell [12], NVT flash has been attracting increasing attentions in the phase equilibrium studies of modern industry.
Despite the enabled prediction on phase equilibrium states by thermodynamic modeling and numerical analysis, the conventional NVT flash calculation methods are usually limited by the necessary step of the phase stability test, which is the guarantee of thermodynamic consistency. Besides, in some of the NVT flash calculation methods, only the local minima on the energy landscape of the system is reached, while the numerical efficiency is challenged in the other methods calculating the global minima of the system [13], [14], [15], [16]. The emerging deep learning algorithm has been successfully applied in accelerating flash calculations [17] with a good performance on both accuracy and efficiency, and both the two steps of phase stability tests and phase splitting calculations have been incorporated in the integrated Thermodynamics-Informed Neural Network, also known as TINN [18]. However, the expensive and limited experimental phase equilibrium data results in the selection of iterative flash calculation results as the ground truth for training and testing in the deep neural network, which places a higher requirement on the reliability of the NVT flash scheme.
A significant breakthrough was proposed in [19], which generated a new thermodynamic analysis approach studying the minimum energy path of the total energy landscape of the two-component two-phase NVT flash system. With a preliminary understanding of the saddle point information and energy landscape of this system, the global minima can be easier to be reached. However, the minimum energy path can only inform the transition route between different local minima on the energy landscape, while the information regarding more saddle points is expected to help generate the complete energy landscape of such system. By calculating the stationary points of one given system, the stable state can be obtained by the 0-order saddle points calculation, but the algorithm stability is always a shortcoming compared with the conventional gradient method converging to the stable state. A number of methods have been proposed on resolving the unstable nature of the saddle points, including the dimer method [22], nudged elastic band method [23], [24], string method [25], [26], [27], gentlest ascent dynamic method [28], [29], minimax method [30] and optimization-based shrinking dimer method [31]. Recently, inspired by the idea of optimization-based shrinking dimer method, a high-order saddle points dynamic method (HiSD) was proposed in [33] to compute saddle points at any orders on the energy landscape of a nonlinear system. This numerical approach can help us to build the solution landscape of the object system, which is a pathway map gathered by all the stationary points on the energy landscape of the system and their connecting paths. This technique has been successfully applied to various fields, such as confined nematic liquid crystal [20], [21], period crystals and quasicrystals [35], diblock copolymers [36] and other applications [37], [38], [39]. Through these applications, the HiSD approach played a very important role to help us to find steady states (local minima) and meta-stable states ( order saddle points) of the numerical model systematically and build the related solution landscape, which can be used to further understand the relevant physical laws and phenomena.
As mentioned in [19], the NVT flash systems with Peng-Robinson equation of state have ill-conditioned Hessian matrices with large condition numbers due to the strong nonlinearity of the equation of state. Normal HiSD approaches need strict step sizes to make the algorithm converge, which significantly slows down the phase equilibrium calculations and further challenges the requirement of obtaining the phase equilibrium state quickly in engineering practice. In order to build the solution landscape of the NVT flash system efficiently in a numerically-compatible manner, in this study, we need to modify the HiSD and make it possible to handle the challenges caused by the very strong stiffness terms in conventional NVT flash models. The string method can efficiently determine the energy transition pathways for complex systems with smooth energy landscape by evolving strings rather than points in the configuration space. This approach consists of two main steps, the evolution of the string by some ordinary differential equation (ODE) solvers and the reparametrization of the string by interpolation. The accuracy of this approach significantly relies on ODE solvers. Moreover, explicit numerical approaches are required because the real scale of the problem is huge. Some common ODE solvers, like Euler and Runge-Kutta methods, require a very small time step for solving ill-conditioned problems. Thus, the attentions have been focused on the exponential time differencing (ETD) scheme, which can provide a satisfactory accuracy and stability for problems with strong stiffness. The scheme involves the exact integration of the governing equation, followed by an explicit approximation of the time integration of the nonlinear term. For the nonlinear problem with strong stiffness, Hochbruck and Ostermann proposed Rosenbrock-type ETD scheme [40], [41] to treat the stiff nonlinear term without the loss of the original stability. This idea proposed a new numerical routine to establish a good linearization approach for nonlinear problems. In this paper, the idea of Rosenbrock-type ETD scheme is applied to build an ETD-type HiSD scheme to calculate the solution landscape of the two-phase two-component VT flash system.
In this paper, the high order dynamic method will be applied to calculate the saddle points in the NVT flash system, so as to generate the total energy landscape of such system to help proceed the thermodynamic analysis on phase equilibrium states, while the manual intervention in the conventional NVT flash framework is expected to be minimized. The fundamental NVT flash scheme is introduced in Section 2, followed by the construction of the high order dynamic method searching for the saddle points. Thermodynamic analysis on the two-component two-phase NVT flash system is presented in Section 3, showing the example saddle points calculations in a fluid mixture, as well as the phase equilibrium analysis based on the saddle points and energy landscape. Concluding remarks are organized in Section 4, and a simple introduction on the Peng-Robinson equation of state, a common-used model in thermodynamic analysis regarding hydrocarbons, is included in the appendix.
Section snippets
Multi-component multi-phase NVT flash scheme
Let us consider a fluid mixture of M components occupying volume V at temperature T. If the investigated mixture is unstable and splits into multiple phases, the total Helmholtz free energy of such system is given by where and denote the mole number vector and volume of phase , respectively.
As the temperature is a constant environmental condition in one defined NVT flash system, T will be dropped in the following
High order saddle point dynamic
To find the stationary points of the energy landscape of the two-phase NVT flash system (2.11), an efficient numerical algorithm is expected to resolve the optimization problem formulated in Equation (2.7) by calculating the saddle points. As explained in [33], inspired by Morse theory [30], [32], finding a k-saddle point of a gradient system can be transformed as finding the solution to the following mini-max problem where V is a k-dimensional subspace spanned by k
Conclusions and discussions
Information regarding phase equilibrium states is of essential importance in the study of multi-component multi-phase fluid flow problems to construct the physically-meaningful initial distribution and to develop the thermodynamically-consistent multi-phase flow models. In this paper, the HiSD algorithm is used for the first time to calculate the saddle points on the energy landscape of the two-component two-phase NVT flash model developed based on the Peng-Robinson equation of state, and the
CRediT authorship contribution statement
Yuze Zhang: methodology, programming, writing.
Xuguang Yang: conceptualization, visualization.
Lei Zhang: reviewing, editing.
Yiteng Li: methodology.
Tao Zhang: methodology, reviewing, supervision.
Shuyu Sun: conceptualization, editing, supervision.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgements
This work is partially supported by the National Key Research and Development Program of China (Grants No. 2021YFF1200500), the National Natural Science Foundation of China (Grants No. 12225102 and 12050002), the Natural Science Foundation of Hunan Province (Grants No. 2020JJ4235), the science research project of Hunan Provincial Education Department (Grants No. 21B0818), the Natural Science Foundation of Jiangsu Province (Grants No. BK20220368), Natural science fund for colleges and
References (43)
- et al.
A coupled Lattice Boltzmann approach to simulate gas flow and transport in shale reservoirs with dynamic sorption
Fuel
(2019) - et al.
Gray-box surrogate models for flash, distillation and compression units of chemical processes
Comput. Chem. Eng.
(2021) - et al.
Phase equilibrium in the hydrogen energy chain
Fuel
(2022) - et al.
Fast and robust algorithm for calculation of two-phase equilibria at given volume, temperature, and moles
Fluid Phase Equilib.
(2013) - et al.
A stable algorithm for calculating phase equilibria with capillarity at specified moles, volume and temperature using a dynamic model
Fluid Phase Equilib.
(2018) - et al.
Based microfluidic fuel cells and their applications: a prospective review
Energy Convers. Manag.
(2022) - et al.
Global stability analysis and phase equilibrium calculations at high pressures using the enhanced simulated annealing algorithm
Chem. Eng. Sci.
(2000) - et al.
Multiphase equilibria calculation by direct minimization of Gibbs free energy with a global optimization method
Comput. Chem. Eng.
(2002) - et al.
A self-adaptive deep learning algorithm for accelerating multi-component flash calculation
Comput. Methods Appl. Mech. Eng.
(2020) - et al.
A comprehensive review of pore scale modeling methodologies for multiphase flow in porous media
Adv. Geo-Energy Res.
(2018)
Subsurface multiphase reactive flow in geologic CO2 storage: key impact factors and characterization approaches
Adv. Geo-Energy Res.
A quantitative study on the approximation error and speed-up of the multi-scale MCMC (Monte Carlo Markov chain) method for molecular dynamics
J. Comput. Phys.
Bivariate free energy and anomalous diffusion
Europhys. Lett.
Intelligent natural gas and hydrogen pipeline dispatching using the coupled thermodynamics-informed neural network and compressor Boolean neural network
Processes
Thermodynamics-Informed Neural Network (TINN) for phase equilibrium calculations considering capillary pressure
Energies
Complex multiphase equilibrium calculations by direct minimization of Gibbs free energy by use of simulated annealing
SPE Reserv. Eval. Eng.
A review of global optimization methods for phase equilibrium modeling and calculations
Open Thermodynam. J.
Thermodynamically-consistent flash calculation in energy industry: from iterative schemes to a unified thermodynamics-informed neural network
Int. J. Energy Res.
Construction of a minimum energy path for the VT flash model by the string method coupled with the exponential time differencing scheme
Commun. Comput. Phys.
Construction of a pathway map on a complicated energy landscape
Phys. Rev. Lett.
Solution landscape of a reduced Landau–de Gennes model on a hexagon
Nonlinearity
Cited by (1)
Numerical Investigation of Natural Convention to a Pseudoplastic Fluid in a Long Channel using a Semi-Implicit Scheme
2023, Applied Sciences (Switzerland)
- ☆
This work is partially supported by the National Key Research and Development Program of China (Grants No. 2021YFF1200500), the National Natural Science Foundation of China (Grants No. 12225102 and 12050002), the Natural Science Foundation of Hunan Province (Grants No. 2020JJ4235), the science research project of Hunan Provincial Education Department (Grants No. 21B0818), the Natural Science Foundation of Jiangsu Province (Grants No. BK20220368), Natural science fund for colleges and universities in Jiangsu Province (Grants No. 22KJD110003), and King Abdullah University of Science and Technology (KAUST) through the grant BAS/1/1351-01-01.