Skip to main content
Engineering LibreTexts

10.1: Multiple Roots and Cubic Behavior

  • Page ID
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \) \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\) \(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)

    It comes as no surprise that cubic equations of state yield three different roots for volume and compressibility factor. This is simply because they are algebraic equations, and any nth order algebraic equation will always yield “n” roots. However, those “n” roots are not required to be distinct, and that is not all: they are not required be real numbers, either. A quadratic expression (n = 2) may have zero real roots (e.g., \(x^2 + 1 = 0\)); this is because those roots are complex numbers. In the case of cubic expressions (n = 3), we will either have one or three real roots; this is because complex roots always show up in pairs (i.e., once you have a complex root, its conjugate must also be a solution.) In our case, and because we are dealing with physical quantities (densities, volumes, compressibility factors), only real roots are of interest. More specifically, we look for real, positive roots such that \(\bar{V} > b\) in the case of molar volume and \( Z > Pb/RT} in the case of compressibility factor.

    In a cubic equation of state, the possibility of three real roots is restricted to the case of sub-critical conditions (\(T < T_c\)), because the S-shaped behavior, which represents the vapor-liquid transition, takes place only at temperatures below critical. This restriction is mathematically imposed by the criticality conditions. Anywhere else, beyond the S-shaped curve, we will only get one real root of the type \(\bar{V} > b\). Figure \(\PageIndex{1}\) illustrates this point.

    Multiple Roots in Cubic EOS showing single-phase, two-phase. TTc Supercritical, T=Tc Critical, T
    Figure \(\PageIndex{1}\): Multiple Roots in Cubic EOS

    Let us examine the three cases presented in Figure \(\PageIndex{1}\):

    1. Supercritical isotherms (\(T > T_c\): At temperatures beyond critical, the cubic equation will have only one real root (the other two are imaginary complex conjugates). In this case, there is no ambiguity in the assignment of the volume root since we have single-phase conditions. The occurrence of a unique real root remains valid at any pressure: any horizontal (isobaric) line cuts the supercritical isotherm just once in Figure \(\PageIndex{1}\).
    2. Critical isotherm (\(T = T_c\)): At the critical point (\(P = P_c\)), vapor and liquid properties are the same. Consequently, the cubic equation predicts three real and equal roots at this special and particular point. However, for any other pressure along the critical isotherm (\(P < P_c\) or \(P > P_c\),) the cubic equation gives a unique real root with two complex conjugates.
    3. Subcritical isotherm (\(T < T_c\)): Predictions for pressures within the pressure range for metastability (\(P_A’ < P < P_B’\)) or for the saturation condition (\(P = P^{sat}\)) will always yield three real, different roots. In fact, this is the only region in Figure \(\PageIndex{1}\) where an isobar cuts the same isotherm more than once. The smallest root is taken as the specific volume of the liquid phase; the largest is the specific volume of the vapor phase; the intermediate root is not computed as it is physically meaningless. However, do not get carried away. Subcritical conditions will not always yield three real roots of the type \(\bar{v} > b\). If the pressure is higher than the maximum of the S-shaped curve, \(P_B\), we will only have one (liquid) real root that satisfies \(\bar{v} > b\). By the same token, pressures between \(0 < P < P_A’\) yield only one (vapor) root. In the case of \(P_A’\) being a negative number, three real roots are to be found even for very low pressures when the ideal gas law applies. This can be seen in Figure 10.1 as well. The largest root is always the correct choice for the gas phase molar volume of pure components.

    Most of these considerations apply to the cubic equation of state in Z (compressibility factor). The most common graphical representation of compressibility factor is the well-known chart of Standing and Katz (Figure \(\PageIndex{2}\)), where Z is plotted against pressure. Standing and Katz presented their chart for the compressibility factor (Z) of sweet natural gases in 1942. This chart was based on experimental data. Graphical determination of properties was widespread until the advent of computers, and thus the Standing and Katz Z-chart became very popular in the natural gas industry. Typical Standing and Katz charts are given for high temperature conditions (\(T > T_c\) or \(T_r > 1\)). Figure \(\PageIndex{2}\), using a cubic equation of state for a pure gas, presents the qualitative behavior of the solution of \(Z\) versus pressure. Isotherms (T > Tc) show the typical qualitative behavior we are accustomed to seeing in the Standing and Katz chart. Cases T < Tc (Tr < 1) are not as familiar to us, as they were not considered by Standing and Katz. For such isotherms, it is clear that you come up with two values of Z (liquid, gas) at saturation conditions.

    Figure \(\PageIndex{2}\): Compressibility Factor versus Pressure

    Contributors and Attributions

    • Prof. Michael Adewumi (The Pennsylvania State University). Some or all of the content of this module was taken from Penn State's College of Earth and Mineral Sciences' OER Initiative.

    This page titled 10.1: Multiple Roots and Cubic Behavior is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Michael Adewumi (John A. Dutton: e-Education Institute) via source content that was edited to the style and standards of the LibreTexts platform; a detailed edit history is available upon request.

    • Was this article helpful?