Jisheng Kou

Jisheng Kou, Shuyu Sun, Xiuhua Wang

In this paper, for the first time we propose two linear, decoupled, energy-stable numerical schemes for multi-component two-phase compressible flow with a realistic equation of state (e.g. Peng-Robinson equation of state). The methods are constructed based on the scalar auxiliary variable (SAV) approaches for Helmholtz free energy and the intermediate velocities that are designed to decouple the tight relationship between velocity and molar densities. The intermediate velocities are also involved in the discrete momentum equation to ensure the consistency with the mass balance equations. Moreover, we propose a component-wise SAV approach for a multi-component fluid, which requires solving a sequence of linear, separate mass balance equations. We prove that the methods preserve the unconditional energy-dissipation feature. Numerical results are presented to verify the effectiveness of the proposed methods.

Jisheng Kou, Shuyu Sun, Xiuhua Wang

The Peng-Robinson equation of state (PR-EoS) has become one of the most extensively applied equations of state in chemical engineering and petroleum industry due to its excellent accuracy in predicting the thermodynamic properties of a wide variety of materials, especially hydrocarbons. Although great efforts have been made to construct efficient numerical methods for the diffuse interface models with PR-EoS, there is still not a linear numerical scheme that can be proved to preserve the original energy dissipation law. In order to pursue such a numerical scheme, we propose a novel energy factorization (EF) approach, which first factorizes an energy function into a product of several factors and then treats the factors using their properties to obtain the semi-implicit linear schemes. We apply the EF approach to deal with the Helmholtz free energy density determined by PR-EoS, and then propose a linear semi-implicit numerical scheme that inherits the original energy dissipation law. Moreover, the proposed scheme is proved to satisfy the maximum principle in both the time semi-discrete form and the cell-centered finite difference fully discrete form under certain conditions. Numerical results are presented to demonstrate the stability and efficiency of the proposed scheme.

Jisheng Kou, Shuyu Sun

A general diffuse interface model with a realistic equation of state (e.g. Peng-Robinson equation of state) is proposed to describe the multi-component two-phase fluid flow based on the principles of the NVT-based framework which is a latest alternative over the NPT-based framework to model the realistic fluids. The proposed model uses the Helmholtz free energy rather than Gibbs free energy in the NPT-based framework. Different from the classical routines, we combine the first law of thermodynamics and related thermodynamical relations to derive the entropy balance equation, and then we derive a transport equation of the Helmholtz free energy density. Furthermore, by using the second law of thermodynamics, we derive a set of unified equations for both interfaces and bulk phases that can describe the partial miscibility of two fluids. A relation between the pressure gradient and chemical potential gradients is established, and this relation leads to a new formulation of the momentum balance equation, which demonstrates that chemical potential gradients become the primary driving force of fluid motion. Moreover, we prove that the proposed model satisfies the total (free) energy dissipation with time. For numerical simulation of the proposed model, the key difficulties result from the strong nonlinearity of Helmholtz free energy density and tight coupling relations between molar densities and velocity. To resolve these problems, we propose a novel convex-concave splitting of Helmholtz free energy density and deal well with the coupling relations between molar densities and velocity through very careful physical observations with a mathematical rigor. We prove that the proposed numerical scheme can preserve the discrete (free) energy dissipation. Numerical tests are carried out to verify the effectiveness of the proposed method.

Jisheng Kou, Shuyu Sun

In this paper, we consider a diffuse-interface gas-liquid two-phase flow model with inhomogeneous temperatures, in which we employ the Peng-Robinson equation of state and the temperature-dependent influence parameter instead of the van der Waals equation of state and the constant influence parameter used in the existing models. As a result, our model can characterize accurately the physical behaviors of numerous realistic gas-liquid fluids, especially hydrocarbons. Furthermore, we prove a relation associating the pressure gradient with the gradients of temperature and chemical potential, and thereby derive a new formulation of the momentum balance equation, which shows that gradients of the chemical potential and temperature become the primary driving force of the fluid motion. It is rigorously proved that the new formulations of the model obey the first and second laws of thermodynamics. To design efficient numerical methods, we prove that Helmholtz free energy density is a concave function with respect to the temperature under certain physical conditions. Based on the proposed modeling formulations and the convex-concave splitting of Helmholtz free energy density, we propose a novel thermodynamically stable numerical scheme. We rigorously prove that the proposed method satisfies the first and second laws of thermodynamics. Finally, numerical tests are carried out to verify the effectiveness of the proposed simulation method.

Jisheng Kou, Shuyu Sun

In this paper, we consider mathematical modeling and numerical simulation of non-isothermal compressible multi-component diffuse-interface two-phase flows with realistic equations of state. A general model with general reference velocity is derived rigorously through thermodynamical laws and Onsager's reciprocal principle, and it is capable of characterizing compressibility and partial miscibility between multiple fluids. We prove a novel relation among the pressure, temperature and chemical potentials, which results in a new formulation of the momentum conservation equation indicating that the gradients of chemical potentials and temperature become the primary driving force of the fluid motion except for the external forces. A key challenge in numerical simulation is to develop entropy stable numerical schemes preserving the laws of thermodynamics. Based on the convex-concave splitting of Helmholtz free energy density with respect to molar densities and temperature, we propose an entropy stable numerical method, which solves the total energy balance equation directly, and thus, naturally satisfies the first law of thermodynamics. Unconditional entropy stability (the second law of thermodynamics) of the proposed method is proved by estimating the variations of Helmholtz free energy and kinetic energy with time steps. Numerical results validate the proposed method.

Xiuhua Wang, Jisheng Kou, Jianchao Cai

The Allen--Cahn equation is one of fundamental equations of phase-field models, while the logarithmic Flory--Huggins potential is one of the most useful energy potentials in various phase-field models. In this paper, we consider numerical schemes for solving the Allen--Cahn equation with logarithmic Flory--Huggins potential. The main challenge is how to design efficient numerical schemes that preserve the maximum principle and energy dissipation law due to the strong nonlinearity of the energy potential function. We propose a novel energy factorization approach with the stability technique, which is called stabilized energy factorization approach, to deal with the Flory--Huggins potential. One advantage of the proposed approach is that all nonlinear terms can be treated semi-implicitly and the resultant numerical scheme is purely linear and easy to implement. Moreover, the discrete maximum principle and unconditional energy stability of the proposed scheme are rigorously proved using the discrete variational principle. Numerical results are presented to demonstrate the stability and effectiveness of the proposed scheme.