{"id":2608,"date":"2018-11-01T13:48:55","date_gmt":"2018-11-01T18:48:55","guid":{"rendered":"http:\/\/www.jmcampbell.com\/tip-of-the-month\/?p=2608"},"modified":"2019-03-07T09:02:58","modified_gmt":"2019-03-07T15:02:58","slug":"how-to-generate-isothermal-and-isobaric-retrograde-regions","status":"publish","type":"post","link":"http:\/\/www.jmcampbell.com\/tip-of-the-month\/2018\/11\/how-to-generate-isothermal-and-isobaric-retrograde-regions\/","title":{"rendered":"How to Generate Isothermal and Isobaric Retrograde Regions"},"content":{"rendered":"

In facilities\u00a0operations, the understanding of where the process is on a phase diagram can often help the engineer and operator to avoid extremely embarrassing design and operating mistakes. A phase diagram is a function of only mixture compositions and components. Though it is expensive and\u00a0time consuming, a phase envelope can be determined experimentally by measuring a series of\u00a0bubblepoints\u00a0and dewpoints. Typically, a phase envelope is generated with a process simulator for specified mixture compositions. For a sensitive process, which requires an accurate phase envelope, some lab or field measured VLE (Vapor-Liquid-Equilibria) data may be required for tuning the equation of state (EOS).<\/p>\n

Process simulation programs typically don\u2019t have a utility or analysis tool to add isothermal and isobaric retrograde curves to a phase diagram. This tip presents two methods for how a process simulation program can be used to quantitatively determine the isothermal and isobaric retrograde curves. A phase diagram, including the retrograde curves for rich gas, is generated and presented.<\/p>\n

Figure 1 from Volume 1 of Gas Conditioning and Processing [1] shows the qualitative phase envelope for a natural gas mixture. Reference [1] also presents the terms used to define the location of various points on the phase envelope. A selection of these terms is listed here.<\/p>\n

 <\/p>\n

\"\"<\/p>\n

 <\/p>\n

 <\/p>\n

<\/p>\n

Figure 1.<\/strong>\u00a0<\/strong>Typical phase diagram for a multi-component natural gas mixture [1]<\/em><\/p>\n

 <\/p>\n

The blue area\u00a0represents\u00a0isothermal retrograde, where liquid condenses when pressure decreases. This is what happens in the path A to B to D.\u00a0 The condensation begins at point B and the maximum liquid occurs at point\u00a0D well\u00a0below the pressure that liquids first start to form. This phenomenon occurs in most naturally occurring natural gas mixtures. We make use of this phenomenon to choose the best operating pressures for our production equipment.<\/p>\n

The green shaded area represents isobaric retrograde vaporization. In this\u00a0area,\u00a0liquid vaporization occurs as temperature decreases or the\u00a0amount\u00a0of liquid decreases as we decrease temperature. If we come into the phase envelope from the other side, the\u00a0amount\u00a0of liquid increases as we increase temperature.<\/p>\n

In facilities operations, the understanding of where the process is on a phase diagram can often help the engineer and operator to avoid extremely embarrassing design and operating mistakes [2]. A phase diagram is a function of only mixture compositions and components. Though it is expensive and time consuming a phase envelope can be determined experimentally by measuring a series of\u00a0bubblepointsand dewpoints. Typically a phase envelope is generated with a process simulator for specified mixture compositions. For a sensitive process, which requires an accurate phase envelope, some lab or field measured VLE (Vapor-Liquid-Equilibria) data may be required for<\/p>\n

Process simulation programs typically don\u2019t have a utility or analysis tool to add isothermal and isobaric retrograde curves to a phase diagram. This tip presents two methods on how to process a simulation program can be used to determine\u00a0quantitatively<\/u>\u00a0the isothermal and isobaric retrograde curves. To demonstrate the procedure the phase diagram including the retrograde curves for a rich gas is generated and presented.<\/p>\n

 <\/p>\n

 <\/p>\n

Determination of Retrograde Curves<\/strong><\/p>\n

Table 1 presents a rich natural gas mixture composition used in this tip. The PR EOS (Peng-Robinson) equation of state [4] with its ProMax [5] default binary interaction parameters\u00a0were\u00a0used.<\/p>\n

 <\/p>\n

Table 1<\/strong>. Natural gas composition<\/p>\n

<\/p>\n

 <\/p>\n

Table 2 and Figure 2 present the calculated properties and generated phase envelope for the mixture in Table 1 by ProMax [5].<\/p>\n

 <\/p>\n

Table 2.<\/strong>\u00a0ProMax (PR EOS) calculated properties for the natural gas mixture of Table 1<\/p>\n

\"\"<\/p>\n

 <\/p>\n

 <\/p>\n

Note Figure 2\u00a0does not illustrate the two retrograde regions. This is true for other process simulation programs. The step-by-step procedures to calculate the isothermal retrograde condensation and isobaric retrograde vaporization curve are presented in the following section.<\/p>\n

<\/p>\n

Figure 2.<\/strong>\u00a0Phase diagram by ProMax (PR EOS) for the natural gas mixture of Table 1<\/em>
\nMETHOD 1: Isothermal Retrograde Condensation (isoT RG) Curve<\/strong><\/p>\n

The curve connecting points M, D, and C in the blue region of Figure 1 is the isothermal retrograde curve. This curve presents the locus of maximum liquid formation for each isotherm between the critical (Point C) and cricondentherm (Point M) temperatures. For a specified temperature and an estimate of pressure as an initial guess, the maximum liquid formation is calculated by adjusting the pressure using the Solver tool in ProMax or similar tool in other software.<\/p>\n

1.<\/strong>\u00a0Calculate a step change for temperature, \u0394T = (Cricondentherm T- Critical T)\/N, N= ~ 20<\/p>\n

2.<\/strong>\u00a0Select the first temperature, e.g. T1<\/sub>= (Cricondentherm T) \u2013 \u0394T<\/p>\n

3.<\/strong>\u00a0Use\u00a0Solver\u00a0tool to find the corresponding pressure, which maximizes the liquid formation<\/p>\n

4.<\/strong>\u00a0Record the specified temperature and calculated pressure corresponding to the maximum liquid formation<\/p>\n

5.<\/strong>\u00a0Select subsequent temperature I, e.g. TI\u00a0<\/sub>= TI-1<\/sub>\u00a0\u2013 \u0394T, and repeat steps 3 and 4 till the critical temperature is reached.<\/p>\n

6.<\/strong>\u00a0Plot the recorded pressures and temperatures on the phase diagram<\/p>\n

 <\/p>\n

 <\/p>\n

The above procedure was executed by defining a simple process flow diagram presented in Figure 3A (isoT RG). The feed composition, flow rate (100 moles), specified temperature, and the initial estimate of pressure are defined in stream 1. The pressure drop in XCHG-100 was set equal to zero. Stream 2 temperature was set equal to stream 1 temperature. Using the Solver tool, stream 2 pressure was adjusted to maximize the liquid formation (this was done by minimizing negative of liquid formation fraction).\u00a0 Table 3 presents the results for a sample calculation.<\/p>\n

 <\/p>\n

<\/p>\n

Figure 3.<\/strong>\u00a0Simple process flow diagrams for isothermal retrograde (A. isoT RG) and isobaric retrograde (B. isoB RG) curve calculations<\/em><\/p>\n

 <\/p>\n

 <\/p>\n

Table 3.\u00a0<\/strong>Single point sample calculation for isoT RG curve.<\/p>\n

\"\"<\/p>\n

Note: * Initial Estimates, ** Calculated Values<\/p>\n

 <\/p>\n

The calculated isothermal retrograde condensation (isoT RG) curve is illustrated\u00a0on\u00a0the phase diagram of Figure 4. Note as the critical point was approached\u00a0better\u00a0estimate of pressures were needed to achieve convergence. Therefore, a few point calculations were done manually.<\/p>\n

 <\/p>\n

 <\/p>\n

<\/p>\n

Figure 4.<\/strong>\u00a0Phase diagram with isothermal retrograde condensation (isoT RG) curve for the natural gas mixture of Table 1<\/em><\/p>\n

 <\/p>\n

METHOD 1: Isobaric Retrograde (isoB RG) Curve<\/strong><\/p>\n

The curve connecting points M and C in the green region of Figure 1 is the isobaric retrograde curve. This curve presents the locus of maximum liquid formation for each isobar (P) between the critical (Point C) and cricondenbar (Point N) pressures. For a specified P and an estimate of\u00a0T as\u00a0initial\u00a0guess, the maximum liquid formation is calculated by adjusting the T.<\/p>\n

1.<\/strong>\u00a0Calculate a step change for P, \u0394P = (Cricondenbar P- Critical P)\/N, N= ~ 20<\/p>\n

2.<\/strong>\u00a0Select the first pressure, e.g. P1<\/sub>= (Cricondentherm P) \u2013 \u0394P<\/p>\n

3.<\/strong>\u00a0Use the Solver tool to find the temperature which maximizes the liquid formation fraction<\/p>\n

4.<\/strong>\u00a0Record the specified pressure and the calculated temperature<\/p>\n

5.<\/strong>\u00a0Select subsequent pressure for\u00a0point\u00a0I, e.g. PI\u00a0<\/sub>= PI-1<\/sub>\u00a0\u2013 \u0394P, and repeat steps 3 and 4 till the critical pressure is reached<\/p>\n

6.<\/strong>\u00a0Plot the recorded pressure and temperatures on the phase diagram<\/p>\n

The above procedure was executed by defining a simple process flow diagram presented in Figure 3B (isoB RG). The feed composition, flow rate (100 moles), specified pressure, and the initial estimate of temperature\u00a0are\u00a0defined in stream 3. The pressure drop in XCHG-101 was set equal to zero, which resulted in the same pressure for streams 3 and 4. Using the Solver tool, stream 4 temperature was adjusted to maximize liquid formation (this was done by minimizing negative of liquid formation fraction). Table 4 presents a sample calculation results. The calculated isobaric retrograde vaporization (isoB RG) curve is illustrated\u00a0on\u00a0the phase diagram of Figure 5.<\/p>\n

 <\/p>\n

Table 4.<\/strong>\u00a0Single point sample calculation for isoB RG curve.<\/p>\n

\"\"<\/p>\n

Note: * Initial Estimates, ** Calculated Values<\/p>\n

 <\/p>\n

 <\/p>\n

<\/p>\n

Figure 5.<\/strong>\u00a0Phase diagram with isothermal (isoT RG) and isobaric (isoB RG) curve for the natural gas mixture of Table 1<\/em><\/p>\n

 <\/p>\n

 <\/p>\n

METHOD 2: Isothermal Condensation and Isobaric Retrograde Vaporization Curves<\/strong><\/p>\n

In this\u00a0method,\u00a0the process simulation program can be used to generate a series of constant liquid (or vapor) quality curves. Figure 5 presents the phase envelope generated by the SRK (Soave-Redlich-Kwong) [6] EOS of EzThermo program [7] and the constant liquid quality curves for 30%, 25%, \u2026.5%, and 2.5% L (liquid) curves. Each constant quality curve has its own cricondentherm and cricondenbar points. Note the isothermal retrograde (diamond symbol) curve passes through each constant quality curve cricondentherm and the isobaric retrograde (triangle symbol) curve passes through each constant quality curve cricondenbar. The properties of the critical point, cricondentherm, and cricondenbar by EzThermo are presented in Table 5. Note that these properties are slightly different from those calculated by the PR EOS of ProMax presented in Table 2 (due to different programs, EOS, and the binary interaction parameters).<\/p>\n

 <\/p>\n

<\/p>\n

Figure 6.\u00a0<\/strong>Phase diagram with isothermal (isoT RG) and isobaric (isoB RG) curve for the natural gas mixture of Table 1 [7]<\/em><\/p>\n

 <\/p>\n

 <\/p>\n

Table 5.<\/strong>\u00a0EzThermo (SRK EOS) calculated properties for the natural gas mixture of Table 1<\/p>\n

<\/p>\n

The generated phase diagram including the isothermal retrograde curve with a single run by EzThermo (SRK EOS) is presented in Figure 7. Note the retrograde curve does not approach the critical point correctly.<\/p>\n

 <\/p>\n

<\/p>\n

Figure 7.<\/strong>\u00a0Phase diagram by EzThermo (SRK EOS) with iso T RG for the natural gas mixture of Table 1<\/em><\/p>\n

 <\/p>\n

 <\/p>\n

Conclusions<\/strong><\/p>\n

Process simulation programs typically\u00a0generate\u00a0phase envelopes without the option to add the isothermal and isobaric retrograde curves.\u00a0 This tip presented two methods of utilizing a process simulation program for quantitative calculations of the isothermal and isobaric retrograde curves. To demonstrate the procedure two diagrams (Figures 5 and 6) including the retrograde curves for a rich gas were presented. The proposed procedure requires several steps but can be incorporated in an excel file linked to the process simulator.<\/p>\n

Ideally, it will be helpful if process simulation programs add these features to their phase envelope utilities\/analysis tools.<\/p>\n

To learn more about similar cases and how to minimize operational problems, we suggest attending ourG4 (<\/strong>Gas Conditioning and Processing)<\/strong><\/a>,<\/strong>\u00a0and\u00a0G5<\/strong>\u00a0(Practical Computer Simulation Applications in Gas Processing)<\/a>\u00a0courses.<\/p>\n

By: Dr. Mahmood Moshfeghian<\/em><\/h4>\n

<\/h1>\n
\n

Sign up to receive Tip of the Month email updates below!<\/i><\/span><\/p>\n