1. Introduction
Liquefying gases is a very energy demanding process. The further away the temperature and pressure of a stream or system are from ambient temperature and pressure (), the higher its exergy. In a liquefaction process, a substance in gaseous state is taken from ambient conditions to a thermodynamically distant state in which it becomes liquid. Therefore, the exergy of the substance at the end of the process will be substantially high.
In this sense, the process of liquefaction is an injection of exergy. Starting from a gaseous state, the substance is given exergy by the liquefaction facility, a compressor-driven thermodynamic cycle designed to remove its sensible and latent heat. Very low temperatures must be reached in order to produce liquid, which in the case of nitrogen can be in the order of 77 K, approximately C.
Removing heat at temperatures so far from the ambient is necessarily exergy consuming. It can be observed in Equation (1) that the input of exergy,
, required per unit of heat removed,
Q, increases steeply the further temperature
T is from the ambient,
:
A thermodynamic cycle for gas liquefaction can be designed around a basic principle based on two flows of gas. The flow of gas which is going to be partially liquefied, “forward flow”, must be at a high pressure. It is cooled down and then expanded at a Joule–Thomson valve coming out as a two-phase stream. The liquid fraction is the product and is extracted from the process. The remaining vapour, which is cold and at a low pressure after expansion, is used to cool down the forward flow. This flow constitutes the “return flow”.
The most direct realisation of this principle is the Linde–Hampson process, which is used for low-scale applications. Liquid air energy storage (LAES) systems consist of an air liquefaction unit for charging a liquid air reservoir and a power unit for discharging it. An analysis of a LAES system based on a modified Linde–Hampson cycle is presented in [
1]. It includes a pre-cooling unit capable of extracting heat out of the process. The system is analysed operating with different degrees of pre-cooling: from the original configuration without pre-cooling to up to 200 K pre-cooling. An increasing trend in the liquid yield is clearly observed. This topology operates at high pressures reaching an optimum liquid yield of approximately 11% at 260 bar. It is concluded that the efficiency of the process can be improved with a Claude topology.
A Claude cycle [
2] includes an expander connecting the forward and return flows of the Linde–Hampson cycle at an intermediate stage of the cooling. Part of the mass flow is kept on the forward flow and part is diverted to the expander. In the expander the pressure and temperature of the expanded mass flow will drop producing work, thus allowing to recover a fraction of the compression duty. The expander outlet is mixed with the returning mass flow of the vapour fraction after expansion, increasing the overall cooling capacity from this point on. With the addition of the expander the Claude cycle operates at significantly lower pressure and produces a higher liquid yield than the Linde–Hampson. A Claude LAES system is analysed in [
3] operating at 40 bar reaching a liquid yield of 22%.
Minor modifications [
4] soon become better options when specific power consumption requirement (per unit mass of liquefied gas) and specific investment costs are considered even at medium scale. The Kapitza topology is a modification of the Claude by which the expander is connected in parallel to the expansion valve [
5], simplifying the design. In [
6], several liquefaction topologies for medium scale cryogenic energy storage (CES) are compared technically and economically: Linde–Hampson, pre-cooled Linde–Hampson, Claude and pre-cooled Kapitza cycles. The pre-cooled Kapitza topology shows very high performance in terms of liquid yield and results the better option from an economic perspective, presenting the lower investment cost and being the most cost effective process. A sensitivity analysis shows a strong dependence of the power consumption per unit of liquid product (specific power consumption (SPC)) on the ratio of the total mass flow diverted to the expander in the Claude and Kapitza cycles. This last shows an optimum at approximately 38% of expanded mass flow.
Large-scale plants for gas liquefaction are based on Kapitza or Collins topologies as indicated in [
7]. Instead of combining a pre-cooling unit with an expander, the Collins topology is a modification of the Claude cycle using two expanders [
8].
A diagram of the standard Collins cycle is shown in
Figure 1. If no mass flow was diverted to the expanders the cycle would operate as a Linde–Hampson. If only one of the expanders was used it would operate as a Claude. Calculating the ratios of mass flow diverted to each of the expanders that minimise SPC is a core design problem.
The gas is first compressed isothermally in C (in a series of adiabatic compression–after-cooling stages), then pre-cooled in E1. After E1 the flow is split. A fraction is diverted to expander Ex1. The remaining flow is further cooled in E2 and E3. At this point it is split again and a fraction is sent to the downstream expander Ex2. The remaining forward flow continues cooling down in E4 and E5 and then flashed in the Joule–Thomson expansion valve. The flow exits the valve in two-phase into the vessel V where the liquid fraction is extracted. The vapour fraction, at low pressure and cold, is recycled into E5 after which it is mixed with the flow expanded in Ex2, heated in E4 and E3, then mixed again with Ex1 outlet and finally heated in E2 and E1. Finally the feed gas is introduced previous to recompression in C.
The temperature of the gas at the inlet of the J–T valve has been specifically indicated in
Figure 1 by
. The liquid yield will be greater the cooler
, so the objective of an industrial cycle will be to achieve the lowest
possible with the lowest SPC.
The compromise between
and SPC is not straightforward. Two aspects must be taken into consideration: first, that the net power consumption of the cycle is the balance between the power consumption of the isothermal compression C, and the power produced at expanders Ex1 and Ex2. Second, that the efficiency of the heat exchange will determine
. While the power consumption at C is defined by the pressure ratio between forward and return flows, the efficiency of the heat exchange is greatly dependent on the share of the total mass flow that is expanded at Ex1 and Ex2. This is because the changes in the mass flow ratio between the forward and return streams caused by diverting mass to the expanders determine their heat capacity ratios, and thus the temperature difference between the streams. In [
9], the heat transfer process is studied for a helium liquefaction cycle. The study analyses the efficiency of the heat transfer process in depth, evaluating the effect of the expanded fractions in Ex1 and Ex2 (see
Figure 1) on the efficiency of the heat exchange. It concludes that efficient heat exchange in E1 and E5 through optimum expanded mass flows must be ensured. The expanded mass flow ratio at Ex1 mainly affects E1, while the expanded mass flow ratio at Ex2 mainly affects E5. In [
7], a discussion of the types of heat transfer in cryogenic cycles is provided. As commented in the text, the heat exchange at E5 may take place in the vicinity of the critical point when liquefying substances such as nitrogen. The specific heat will show large variations during the process, leading to an irregular temperature profile of the forward flow, so that the minimum temperature approach (MITA) in E5 will take place in between both ends.
Several studies analyse how the share of expanded flow affect the performance of a standard Collins cycle. In [
10], the expanded mass flow ratios for minimum SPC are calculated for helium liquefaction. The optimum expanded mass flow ratios are 45% for Ex1 and 35% for Ex2, adding a total expanded mass flow of 80% of the total flow at compressor C. The total expanded mass flow ratio can be reduced if the heat exchange is increased, either by more effective heat exchangers or external pre-cooling, down to 75–79%. Similar results are obtained in a more recent study for an analogous Collins helium liquefaction cycle in [
11]. Expanded flow ratios for maximum liquid yield are calculated, with similar assumptions to those of the previous study. The cycle is simulated with Aspen HYSYS across a range of expanded mass flow ratios. Although the main conclusion of the study is similar, that the optimum expanded mass flow is 80%, the resulting share between Ex1 and Ex2 is 50%, differing from the larger Ex1 ratio of the previous study. The article points out that results are similar for uneven distributions between the expanders up to 40–60% or 60–40%, and that distributions with higher Ex1 flow provide more greater liquid yield. The analysis of the heat process in this study, by the same author as [
9], is coherent with the observations about the importance of E1 and E5 for the efficiency of the cycle concluded in [
9].
As it has been mentioned, nitrogen liquefaction presents some particularities due to the large variations of the specific heat in the vicinity of the critical point. However, few specific assessments exist in the literature. In [
12], a Kapitza cycle is simulated. The influence of diverting mass to the expander on the heat exchange process is identified, along the lines commented above and in [
7]. A small-scale pre-cooled Kapitza cycle is analysed in [
13]. Due to the scope of the study, the compressor discharge pressure is limited to 8 bar. Liquid yield results slightly over 6%. As the study mentions, this exceptionally low value would increase with higher compression ratios. The study provides the temperature profiles of the forward and return flows. It is interesting to observe that at 8 bar, well below the critical for nitrogen (33.9 bar), no irregularities due to large variations of specific heat can be observed. The large exergy losses in the aftercooler are identified as a major cause of exergy destruction, 26.03% of the total. The study [
14] analyses an industrial pre-cooled Kapitza cycle for nitrogen liquefaction, operating at a compressor discharge pressure of 43.27 bar, sufficiently close to the critical point to observe the large variations in specific heat. The relation between heat transfer at E5 and the expanded mass flow ratio at Ex2 is assessed, analysing the variation of the temperature profiles for off-design operation.
In this article, an industrial Collins-based nitrogen liquefaction process is analysed in detail in order to find the optimum expanded flows that minimise SPC. The cycle does not correspond to the standard cycle of
Figure 1. As any industrial cycle, it must be adapted to the requirements of the product and integrate modifications for large-scale liquefaction [
15]. The high pressure level of the cycle is 36 bar, close to the critical 33.9 bar, so in contrast to the helium case, the large variations in specific heat mentioned earlier will affect the last stage of the heat exchange and appear clearly on the temperature profiles. Although the underlying Collins topology explains a similar behaviour of the studied cycle to the standard, the resulting SPC and optimum expanded mass flow show differences.
A detailed explanation of the cycle and its modifications over the standard Collins cycle will be developed in
Section 2.1. The last two compression stages have been coupled to the two expansions in the cycle for better efficiency, as would happen in a real industrial cycle, using
companders (explained in
Section 2.1) instead of independent compressor stages and expanders [
16]. The topology also integrates a subcooling unit at the end of the cycle, necessary to prepare the liquid product for transportation and distribution (to reduce boil-off). It consists in an additional heat exchanger and a second expansion valve, where the previously liquefied nitrogen is flashed again and subcooled. The low pressure gaseous fraction is recycled, which requires an additional compressor. The cycle has four pressure levels distributed in two forward and two return flows, instead of the two of the standard cycle. The pressure of the nitrogen gas feed to the process is also adapted to the typical industrial value of 4 bar, instead of the usual value employed in other studies of approximately 1 bar. Global assumptions and design parameters are described in
Section 2.2.
Section 3 describes the series of simulations that have been carried out in Unisim Design R451. The expanded mass flows that minimise SPC have been calculated and a global optimum has been identified
Section 3.1. This case has been analysed in depth and a brief sensitivity analysis to component efficiencies developed in
Section 3.2.
An exergy characterisation of the process has been carried out in
Section 3.3. The rational exergy efficiency of the process and the exergy destruction in each component have been calculated. A Sankey diagram illustrating the exergy flows in the process have been generated.
2. Methodology
Figure 2 shows the component diagram of the industrial N
liquefaction cycle, which will be developed in detail in
Section 2.1. The cycle will be simulated across a range of different mass flow ratios,
and
, in order to find the lowest SPC configuration for a fixed set of process assumptions presented in
Section 2.2. It is the same design problem as in [
10,
11] for helium but with a more complex topology in which the highest pressure level, HHP, is actually determined by
and
, due to the use of
companders (see
Section 2.1).
The optimisation problem for this cycle and the role of the
companders will be summarised in
Section 2.3, in which the similarities with the underlying Collins cycle can be observed. The process is described in terms of exergy in
Section 2.4, and a
diagram is shown in
Figure 3.
2.1. Thermodynamic Cycle
The diagram of the cycle is shown in
Figure 2. It can be observed that it is a modified Collins cycle with four pressure levels: LP, MP, HP and HHP. The main forward flow is ‘c’-‘p’, cooled by the MP and LP return flows (“q”–“a” and “t”–“z”, respectively). The flows diverted to expander T-101,
, and to T-102,
, increase the cooling power of the MP return flow. The main expansion valve is V-101, flashing two-state N
into the F-101 vessel. A subcooler after F-101 is required to adapt the product for transportation. This is done with E-102 and a second expansion valve V-102.
A main N compressor K-101 with an aftercooler C-101 compresses the MP Nitrogen to HP level. It is then split into two fractions: The first one is cooled at the first cold box section, E-101.1 and expanded in T-101 to MP level into the MP return flow.
The second one is further compressed to HHP and cooled in K-102, C-102, K-103 and C-103. These compressors are powered by the expanders downstream, in companders, as will be explained later. The compressed gas leaving C-103 at state “j” is then cooled at E-101.1-3 and then split again, diverting ratio to expander T-102, where it expands to MP level into the F-101 vessel. The remaining HHP flow is further cooled at E-101.4 and finally flashed at the Joule-Thomson valve V-101 into the F-101 vessel.
The vessel F-101 is at MP level. There the liquid and vapour phases are separated. Vapour is returned to the start through the E-101 series of exchangers. The liquid fraction is subcooled at E-102, then split into two, one part is the actual dispatched product. The other stream is flashed to LP level at V-102, producing the temperature drop required for the subcooling, so after flashing it is returned to E-102 and from there, through the series of E-101 exchangers to the start. Logically, before it is integrated with the feed and MP return flows, it has to be compressed back to MP at K-104 and cooled at C-104.
In a
compander set, compressor and expander blades are fixed to the same shaft forcing them to rotate together (there can exist two shafts coupled with a gearbox to adapt rotation speeds). The expanding flow will generate power through the expander section which in turn, via the shaft, will drive the associated compression stage. Both streams are independent so they can be connected to different processes altogether, or to different points within the same process. In this cycle, in order to use the power produced by the expanded mass flow compressor K-102 and expander T-101 are integrated into a
compander. The same happens with K-103 and T-102, integrated into a second
compander. In this way, the compression duty is provided by the power produced at a different point of the cycle directly, instead of passing through an intermediate conversion to electricity, making the process more efficient and avoiding electric motors and generators, thus reducing cost. An example of a nitrogen cryogenic cycle using companders can be found in [
17].
From the point of view of analysing the cycle, however, companders introduce an additional restriction that does not exist in the original Collins cycle: the power of the downstream expansions equal the power of the corresponding compression stages upstream. Therefore, increasing the maximum pressure of the cycle cannot be done independently, as it will necessarily be associated with greater amount of gas expansion downstream. This will alter mass flows and operating temperatures affecting the heat transfer between forward and return flows, the key factor for cycle performance.
2.2. Assumptions and Cycle Parameters
The following assumptions have been made.
The system is in steady state.
Component efficiencies do not vary with pressure, temperature and mass flow rate.
The cold box is adiabatic.
Pressure drops in pipelines are negligible.
These assumptions are used in order to isolate the behaviour of the cycle from the particularities of any specific equipment, and they are common in directly related studies [
10,
11]. Constant component efficiencies allow designing the optimum cycle configuration, which will be the target of this work (analysing a given facility across its operating range would require variable component efficiency).
Table 1 details the main component specifications. Polytropic efficiency was used because it is more representative of the actual technology level and its value does not depend on the pressure ratio. The efficiency of compander compressor stages is lower than in the main compressor, because the coupling with the expander will limit performance. Performance of heat exchangers is defined by manufacturer values of Minimum Temperature Approach (MITA), within the 0 to 5.36
C range specified in [
18].
The feed is assumed to be 100% N at medium pressure, MP = 4 bar. Pressure levels correspond to LP = 1 bar, MP = 4 bar, HP = 18 bar and HHP resultant.
The simulations have been carried out with Honeywell Unisim R451. Nitrogen properties have been calculated using a modified Peng–Robinson equation of state. The original Peng–Robinson equation [
19] showed a remarkable accuracy for nitrogen. A number of improvements have been developed since [
20]. It is of widespread use for pure nitrogen [
13,
21] and nitrogen mixture [
22,
23,
24,
25] cryogenic processes. Unisim Design R451 uses a modified Peng-Robinson EOS specifically recommended for all components of air in the range of pressure and temperature conditions of this article [
26].
2.3. Optimisation Problem
The specific power consumption of the process is given by the duty of the main compressor, K-101, and the LP return compressor, K-104. As it will be explained in
Section 2.1, the other two compressors K-102 and K-103 do not influence SPC directly:
where
and
indicate main and low pressure cycle compressor power consumption, K-101 and K-104, respectively, as will be developed in
Section 2.1, and
indicates total liquefied mass flow.
As it has been introduced in
Section 1, there are two opposing effects in the Collins cycle (it would equally apply to Kapitza and other cycles). They also exist in the industrial cycle analysed here, although some details must be taken into account.
The temperature at the inlet of the expansion valve V-101 (state “o” in
Figure 2) will be indicated by
for clarity. It has to be as low as possible in order to maximise the liquid fraction,
, after flashing (state “p”).
This temperature is achieved by cooling the gas from state “c” (after the main aftercooler, C-101) in the series of heat exchangers, in which the cold, return flow extracts heat from the main, forward flow from “c” to “o”. The return flow is formed by two major streams. At medium pressure (MP), the non-liquefied fraction (which should be as small as possible) into which the outlet from the two expanders are injected. At low pressure (LP), the outlet from the subcooler (E-102 + V-102).
The greater the return mass flow, the greater the cooling prior to the valve (state ‘o’ in
Figure 2), which will lower
. The return flow can be increased by diverting as much nitrogen as possible to the expanders. Due to the
companders linking T-101 and T-102 with K-102 and K-103 (see
Section 2.1), changes in expansion power will also lead to a higher HHP pressure. This effect will tend to move state “o” further left toward the liquid region in
Figure 3, thus increasing the liquid yield,
.
The optimisation problem consists in calculating the values of and which minimise SPC. In appearance, this would mean that as much mass flow as possible should be diverted to the expanders, which equals large values of and . However, this has an opposite effect that will become significant over a certain limit. The greater the expanded fraction is, the less nitrogen mass flow at point ‘o’ left to flash in V-101. When it is flashed, the liquid fraction, , will be high but the actual mass flow small. So there must exist a certain ratio of the mass flows (, ) that, when sent to the expanders, will generate enough cooling while maintaining a sufficient mass flow through the expansion valve and minimises SPC.
2.4. Exergy Fundamentals
The general exergy balance for a given open volume
j at rest, with streams
i entering or leaving through the boundary, during a time differential
can be formulated [
27]:
where
is the change in total exergy (exergy+mechanical exergy) of the control volume,
the exergy content of the heat received during the process, being
T the temperature of the control volume and
that of the environment,
the exergy content of the work performed by the control volume and
the exergy destruction in the control volume.
is the balance of exergy entering and exiting the control volume through all inlets and outlet flows, where
is the exergy of flow
i:
where
are specific enthalpy and entropy at flow
i, respectively.
The liquefaction cycle analysed here can be taken as a control volume with one nitrogen input, MP GAN, and one output, MP LIN, referring to
Figure 2. The minimum exergy required for the process would be the difference in flow exergy between these two states: 173.4 kWh/ton. This is achievable by the blue path [
7] in the
chart for nitrogen shown in
Figure 3, which consists of an isothermal compression and an isentropic expansion, both reversible thus having zero exergy destruction, but technologically impossible. The real process has also been represented in the figure. The exergy analysis developed in
Section 3.3 will compare both processes.
In order to obtain the exergy breakdown of the process, Equation (
2) has to be applied to all components, with the assumptions of
Section 2.2. The resulting expressions for exergy destruction at each component are summarised in
Table 2.