Peng–Robinson equation

Peng–Robinson equation
Name	Peng–Robinson equation
State	solid/liquid/gas
Developer	Ding-Yu Peng, Donald B. Robinson
Year	1976
Field	Thermodynamics
Applications	Petroleum industry, Chemical engineering, Natural gas processing, Phase equilibria

Peng–Robinson equation The Peng–Robinson equation is a cubic equation of state developed to model the P–V–T behavior of nonpolar and mildly polar fluids for applications in Petroleum industry, Chemical engineering, Natural gas processing, Refining, and Cryogenics. It was introduced by Ding-Yu Peng and Donald B. Robinson in 1976 and has become a standard tool alongside other cubic equations of state used by practitioners at institutions such as ExxonMobil, Shell plc, Chevron Corporation, BP, and academic groups at Massachusetts Institute of Technology, University of Cambridge, and Imperial College London. The model balances accuracy and computational efficiency, enabling its adoption in process simulators from vendors like AspenTech, Honeywell Process Solutions, and Schneider Electric.

Introduction

The equation was proposed in response to limitations observed with earlier cubic models developed by engineers at Rudolf Clausius-era research and later formalized by developers of the van der Waals equation, Redlich–Kwong equation, and Soave modification. Peng and Robinson presented an empirical attractive term calibrated using critical properties and acentric factors from data compiled by organizations such as National Institute of Standards and Technology and industry consortia in the American Petroleum Institute. Its uptake was reinforced by validation studies published in journals like AIChE Journal and Chemical Engineering Science, and by its use in textbooks authored by academics at University of California, Berkeley and University of Texas at Austin.

Equation and Parameters

The Peng–Robinson model expresses pressure as a cubic polynomial in molar volume, with parameters tied to critical temperature and pressure data that are routinely reported by agencies such as NIST, Society of Petroleum Engineers, and national laboratories including Sandia National Laboratories. The formulation introduced coefficients derived from empirical fits similar in spirit to approaches used by John D. Van der Waals and later adapted by G. Soave in the Soave–Redlich–Kwong work. The equation’s temperature-dependent attraction parameter uses the acentric factor concept developed from datasets maintained by International Union of Pure and Applied Chemistry and experimental compilations by ASTM International. Practitioners at DOW Chemical Company, BASF, and DuPont commonly use the reported critical constants and acentric factors to compute the Peng–Robinson parameters for hydrocarbons and light components encountered in operations overseen by regulators like United States Environmental Protection Agency and Norwegian Petroleum Directorate.

Thermodynamic Properties and Applications

Peng–Robinson predictions feed calculations of phase equilibria, fugacity coefficients, enthalpy, entropy, and speed of sound used in process design at firms such as Fluor Corporation and Bechtel. It supports vapor–liquid equilibrium computations in separation units designed by engineering consultants affiliated with KBR, Inc. and is implemented in reservoir simulation workflows used by operators like ConocoPhillips and TotalEnergies. Academic investigations comparing Peng–Robinson results against experimental data from facilities at National Research Council (Canada) and Argonne National Laboratory often reference benchmarking efforts coordinated by International Energy Agency. The model underpins safety analyses conducted under standards from American Society of Mechanical Engineers and informs hydrate formation studies relevant to Offshore Petroleum Engineering.

Comparison with Other Cubic Equations of State

Peng–Robinson is frequently compared to the Redlich–Kwong equation, the Soave–Redlich–Kwong equation, and modifications of the van der Waals equation, with comparative assessments published by researchers at University of Alberta, University of Texas at Austin, University of Calgary, and Texas A&M University. Industry evaluations by Shell plc and ExxonMobil have highlighted trade-offs in liquid density prediction, critical region behavior, and computational robustness versus models like the Benedict–Webb–Rubin equation and multi-parameter empirical correlations developed at National Physical Laboratory (UK). International standards bodies, including ISO committees and API task groups, often recommend model selection guidelines that reference these comparative studies.

Limitations and Modifications

Despite broad utility, the Peng–Robinson equation shows systematic deviations for highly polar substances and complex associating fluids studied by research groups at ETH Zurich, University of Manchester, and Caltech. To address such limitations, extensions incorporate volume translation techniques proposed by researchers at University of Strasbourg and activity coefficient coupling methods from researchers affiliated with University of Barcelona and Peking University. Hybrid models combining Peng–Robinson with equations like the Soave–Redlich–Kwong equation variants, or with association models such as Statistical Associating Fluid Theory, have been developed in collaborations involving Chevron Research and university labs at University of Alberta and National University of Singapore.

Numerical Implementation and Solution Methods

Solution of the Peng–Robinson cubic involves root-finding algorithms used in simulation packages from AspenTech, Honeywell Process Solutions, and open-source projects maintained by contributors from MIT and University of California, Berkeley. Common numerical methods include analytic cubic solvers traced to techniques taught at Princeton University and iterative methods like Newton–Raphson used in process control codes at Siemens, ABB Group, and Emerson Electric. Robust phase-envelope tracing and flash calculations rely on algorithms validated by interlaboratory comparisons organized by NIST and computational standards developed at IEEE conferences. High-performance implementations for reservoir-scale simulations use parallel solvers and preconditioning strategies studied at Lawrence Berkeley National Laboratory and deployed by companies such as Schlumberger and Halliburton.

Category:Equations of state