Unconditionally Energy Stable Linear Schemes for the Diffuse Interface Model with Peng–Robinson Equation of State

Hongwei Li, Lili Ju, Chenfei Zhang, Qiujin Peng

Research output: Journal PublicationArticlepeer-review

40 Citations (Scopus)


In this paper, we investigate numerical solution of the diffuse interface model with Peng–Robinson equation of state, that describes real states of hydrocarbon fluids in the petroleum industry. Due to the strong nonlinearity of the source terms in this model, how to design appropriate time discretizations to preserve the energy dissipation law of the system at the discrete level is a major challenge. Based on the “Invariant Energy Quadratization” approach and the penalty formulation, we develop efficient first and second order time stepping schemes for solving the single-component two-phase fluid problem. In both schemes the resulted temporal semi-discretizations lead to linear systems with symmetric positive definite spatial operators at each time step. We rigorously prove their unconditional energy stabilities in the time discrete sense. Various numerical simulations in 2D and 3D spaces are also presented to validate accuracy and stability of the proposed linear schemes and to investigate physical reliability of the target model by comparisons with laboratory data.

Original languageEnglish
Pages (from-to)993-1015
Number of pages23
JournalJournal of Scientific Computing
Issue number2
Publication statusPublished - 1 May 2018
Externally publishedYes


  • Diffuse interface
  • Energy stability
  • Invariant energy quadratization
  • Linear scheme
  • Penalty formulation
  • Peng–Robinson equation of state

ASJC Scopus subject areas

  • Software
  • Theoretical Computer Science
  • Numerical Analysis
  • Engineering (all)
  • Computational Theory and Mathematics
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'Unconditionally Energy Stable Linear Schemes for the Diffuse Interface Model with Peng–Robinson Equation of State'. Together they form a unique fingerprint.

Cite this