Optimization Study on Fluids for the Gravity-Driven Organic Power Cycle

: The organic Rankine cycle (ORC) is efﬁcient in using low-grade heat energy, while low pump efﬁciency and high pump leakage are usually serious problems. A gravity-driven organic power cycle (GDOPC) uses gravity instead of a pump to pressurize working ﬂuid and has the potential to avoid problems associated with the pump. A theoretical method is used to study the performance and suitability of several ﬂuids for GDOPC. The results show that the ﬂow efﬁciency in high vertical pipes and the pump efﬁciency determine whether GDOPC gives better performance than ORC or not. When R245fa is selected as working ﬂuid and evaporating temperature is 62 ◦ C, speciﬁc energy of GDOPC (ﬂow efﬁciency is 80%) is 2.5% higher than that of ORC (pump efﬁciency is 60%). The improvement degree of speciﬁc energy and the liquid column height increase with increasing evaporating temperature. R1234yf and R227ea give good performance with speciﬁc energy of 4.84 kJ/kg and 4.82 kJ/kg, respectively, while they need a liquid column as much as 76.55 m and 45.65 m, respectively. Although R365mfc and cyclopentane do not give the most excellent cycle performance, they need liquid column height as low as 9.04 m and 10.88 m, respectively. Fluid with low saturated pressure and high density may need low liquid column height and has the advantage to be used in practical applications.


Introduction
Nowadays, the energy crisis is becoming more and more serious. However, the reserve of low-grade heat energy is huge in nature and the society. It is helpful for alleviating an energy shortage and protecting the environment to use low-grade heat energy. A power cycle using unconventional working fluid has good application prospects in using low-grade heat energy [1,2].
A lot of investigation was carried out on optimization of working fluid and operating conditions for different heat sources. Shen et al. [3] developed an equation of state named crossover volume translation Soave-Redlich-Kwong to predict the thermodynamic properties of several hydrocarbons. Generally, the dry and isentropic fluids are most suitable for the organic Rankine cycle (ORC) [4]. Additionally, extensive research was executed to optimize working fluid for different heat sources [5][6][7]. Pan et al. [8] improved the theoretical analysis method of ORC based on radial flow turbine and gave the performance of several fluids for ORC using 90 • C hot water as a heat source. R227ea gave the largest net output power and the lowest thermal efficiency while R141b gave the least net output

Introduction of Gravity-Driven Organic Power Cycle (GDOPC)
In a subcritical power cycle system, fluid is at either the liquid phase state or vapor phase state. Fluid density at the liquid phase state is much higher than that at the vapor phase state. Because of gravity, pressure at the liquid column bottom is much higher than that at the vapor column bottom. This pressure difference between liquid column bottom and vapor column bottom is used to drive fluid in a GDOPC. As shown in Figure 2, the process from state point 1 to state point 2 represents that vapor phase fluid with high temperature and high pressure expands in the turbine and outputs shaft power. The process from point 2 to point 3 is that fluid with low temperature and low pressure flows upwards along vertical vapor pipe (pressure at the top of the vapor column is slightly lower than that at the bottom). The process from state point 3 to state point 5 is that fluid is cooled and condensed by cooling water. For some dry fluids (with positive slope for saturated gas line), there is a large superheated degree at state point 3, so a larger condenser is needed because of the low single-phase heat transfer coefficient. The dew point (point 4) appears somewhere in the condenser. Th process from state point 5 to state point 6 is that liquid fluid flows downwards and forms a high liquid column. Under the action of gravity, pressure at the bottom of the liquid column increases and then a higher pressure-under which evaporation happens later-is reached here. The evaporating process is represented by process from state point 6 to state point 1. The bubble point and dew point appear somewhere in the evaporator.
In this article, the turbine is located at the bottom of the system as shown in Figure 2. However, it also can be located at the top of the system [28]. If the turbine is located at the top, the positions of the process from point 2 to point 3 and the process from point 1 to point 2 as shown in Figure 2b

Introduction of Gravity-Driven Organic Power Cycle (GDOPC)
In a subcritical power cycle system, fluid is at either the liquid phase state or vapor phase state. Fluid density at the liquid phase state is much higher than that at the vapor phase state. Because of gravity, pressure at the liquid column bottom is much higher than that at the vapor column bottom. This pressure difference between liquid column bottom and vapor column bottom is used to drive fluid in a GDOPC. As shown in Figure 2, the process from state point 1 to state point 2 represents that vapor phase fluid with high temperature and high pressure expands in the turbine and outputs shaft power. The process from point 2 to point 3 is that fluid with low temperature and low pressure flows upwards along vertical vapor pipe (pressure at the top of the vapor column is slightly lower than that at the bottom). The process from state point 3 to state point 5 is that fluid is cooled and condensed by cooling water. For some dry fluids (with positive slope for saturated gas line), there is a large superheated degree at state point 3, so a larger condenser is needed because of the low single-phase heat transfer coefficient. The dew point (point 4) appears somewhere in the condenser. Th process from state point 5 to state point 6 is that liquid fluid flows downwards and forms a high liquid column. Under the action of gravity, pressure at the bottom of the liquid column increases and then a higher pressure-under which evaporation happens later-is reached here. The evaporating process is represented by process from state point 6 to state point 1. The bubble point and dew point appear somewhere in the evaporator. should be exchanged. Usually, a GDOPC system needs a very high liquid column which is determined by heat source parameters and working fluid. Because of the large height from the geothermal heat source to the ground, the GDOPC may have the potential to use geothermal energy.

Considered Fluids
Fluid containing chlorine is harmful for ozone, so only HFCs, HFOs, CFs and HCs are selected as working fluid. The ODP of all considered fluids is 0. The basic thermal properties and environmental properties of considered fluids are shown in Table 1 [29]. Also, these fluids are In this article, the turbine is located at the bottom of the system as shown in Figure 2. However, it also can be located at the top of the system [28]. If the turbine is located at the top, the positions of the process from point 2 to point 3 and the process from point 1 to point 2 as shown in Figure 2b should be exchanged. Usually, a GDOPC system needs a very high liquid column which is determined by heat source parameters and working fluid. Because of the large height from the geothermal heat source to the ground, the GDOPC may have the potential to use geothermal energy.

Considered Fluids
Fluid containing chlorine is harmful for ozone, so only HFCs, HFOs, CFs and HCs are selected as working fluid. The ODP of all considered fluids is 0. The basic thermal properties and environmental properties of considered fluids are shown in Table 1 [29]. Also, these fluids are classified by a new ACZMN sequence method [30] with 20 • C as lowermost temperature. In the method, each letter of the term ACZMN stands for one point on the saturated curve of the fluid.  Figure 3a is a whole block diagram of the theoretical analysis method for GDOPC. For general fluids, a very high liquid column (high system) is usually needed to prove the cycle pressure difference between evaporating pressure and condensing pressure. Although vapor has much less density than liquid, the pressure drop of the vapor column from the bottom (state point 2) to the top (state point 3) cannot be ignored. Additionally, the kinetic energy loss exists along the vapor column and the liquid column. The pressure drop of the vapor column and the kinetic energy loss are all taken into account in the theoretical method. The computation methods about the liquid column and the vapor column are shown in Figure 3b,c, respectively. Therefore, processes in a vertical vapor pipe and vertical liquid pipe are non-isentropic, which can be seen from Figure 2b. Input heat source data, cooling source data, efficiency data, pinch point data, state data and acceleration due to gravity.    In this article, REFPROP 9.0 [31] is used to obtain fluid thermal properties. For example, when a few state properties are provided, such as temperature and pressure for single phase, REFPROP 9.0 can output other thermal properties, such as enthalpy, entropy, density and specific heat. In the theoretical analysis method, the turbine efficiency is designated. Fluid state parameters at the turbine exit can be computed using turbine efficiency and fluid state parameters at the turbine entrance according to Equation (1). Fluid parameters like temperature and density change along the pipe, so the process from state point 5 to state point 6 is divided into several units. The liquid fluid flow efficiency is defined as Equation (2) in order to take account of the kinetic energy loss along the liquid vertical pipe. The fluid liquid column height of each unit as well as the whole process in the vertical pipe can be expressed as Equation (3). The pressure difference between evaporating pressure and condensing pressure is completely provided by the fluid liquid column. The detailed computation method is shown in Figure 3b. In this article, REFPROP 9.0 [31] is used to obtain fluid thermal properties. For example, when a few state properties are provided, such as temperature and pressure for single phase, REFPROP 9.0 can output other thermal properties, such as enthalpy, entropy, density and specific heat. In the theoretical analysis method, the turbine efficiency is designated. Fluid state parameters at the turbine exit can be computed using turbine efficiency and fluid state parameters at the turbine entrance according to Equation (1). Fluid parameters like temperature and density change along the pipe, so the process from state point 5 to state point 6 is divided into several units. The liquid fluid flow efficiency is defined as Equation (2) in order to take account of the kinetic energy loss along the liquid vertical pipe. The fluid liquid column height of each unit as well as the whole process in the vertical pipe can be expressed as Equation (3). The pressure difference between evaporating pressure and condensing pressure is completely provided by the fluid liquid column. The detailed computation method is shown in Figure 3b.

Analysis Method of GDOPC
The process from state point 2 to state point 3 is also divided into several units and the vapor fluid flow efficiency is defined as Equation (4). The vapor column height in each unit as well as the whole vertical pipe can be expressed as Equation (5). The pressure values of state point 1 and state point 3 can be easily obtained by computing the saturated pressure at evaporating temperature and condensing temperature. However, when the pressure difference of the vapor column and the kinetic energy losses are all taken into account, state point 2 parameters must be obtained by iterative method which is detailedly shown in Figure 3c.
The heat capacity in the evaporator can be expressed as Equation (6). Power is not inputted into the GDOPC system, so the net output power is the same as the turbine output power. Therefore, the thermal efficiency can be expressed as Equation (7a). The ORC system is usually integrated with a fluid pump which consumes power, so its thermal efficiency is expressed as Equation (7b).
Power output per unit mass flow rate of the heat source fluid can represent the benefits from the heat source and is usually selected as the objective function. In this article, the specific energy (E specific ) is defined as Equation (8a) for GDOPC while it is defined as Equation (8b) for conventional ORC. The unit of specific energy is kJ/kg. The higher the specific energy, the more the benefit can be obtained from the heat source.

Comparison between GDOPC and Organic Rankine Cycle (ORC)
The significant advantage of GDOPC is that the fluid pump is omitted in the system. Although a pump consumes much lower power than turbine outputs, the seal problem of the pump is solved and pressurizing efficiency can be increased slightly. These improvements can increase the thermal efficiency and decrease the system's cost. The flow efficiency along high vertical pipes determines the thermal efficiency and the vertical pipes' impacts on the system cost. In addition, the new system is a little more complicated than ORC system, because of the high vertical pipes. The pressurizing process can be self-adapting. In other words, the liquid column height varies with the difference between evaporating pressure to condensing pressure if the vertical pipes are high enough.
In this section, a comparative study between GDOPC and ORC is carried out using R245fa as working fluid in both cycles. Specified conditions are shown in Table 2. As shown in Figure 4, for both cycles, there is a maximum specific energy appearing in the condition with evaporating temperature of about 62 • C. For ORC, specific energy increases with pump efficiency, especially near the optimal condition. In considered conditions, the specific energy of GDOPC is a little higher than that of ORC. The average absorbing heat temperature in the evaporator increases with the increase of evaporating temperature. Consequently, the cycle thermal efficiency increases with increasing evaporating temperature, which is helpful for enhancing specific energy. The hot water outlet temperature increases with the increase of evaporating temperature. Accordingly, heat capacity in the evaporator decreases with increasing evaporating temperature, which is harmful for increasing specific energy. Therefore, the maximum specific energy appears under the influence of the two factors. In the conventional ORC system, a pump is used to pressurize and transport the fluid. The net power of the system is lower than the turbine output power because of electric power consumed by the fluid pump. The net power of GDOPC is equal to the turbine output power because input power is not consumed by other instruments. But there is kinetic energy loss along the high vertical pipes, which decreases the turbine output power. The essence of GDOPC is pressurizing fluid directly by heat energy. As shown in Figure 5, it is obviously that the specific energy is improved by using GDOPC in considered conditions. The lower the pump efficiency in the ORC, the higher the improvement degree. The advantage becomes more and more obvious with increasing evaporating temperature. At the optimal condition with evaporating temperature of 62°C, the improvement degree is 2.5% with the pump efficiency of 60% in the ORC. In a practical application, the pump efficiency is usually about 60%. An appropriate efficiency of 80% is used for the flow efficiency in the high vertical pipes. In considered conditions, the higher the evaporating temperature, the higher the improvement degree. However, a higher fluid liquid column is needed in the conditions with higher evaporating temperature, as shown in Figure 6. For example, in the condition with heat source temperature of 150 °C, evaporating temperature of 130 °C and working fluid of R245fa, the improvement degree of specific energy and cycle thermal efficiency can reach 7.1% and 6.8% while the fluid liquid column As shown in Figure 5, it is obviously that the specific energy is improved by using GDOPC in considered conditions. The lower the pump efficiency in the ORC, the higher the improvement degree. The advantage becomes more and more obvious with increasing evaporating temperature. At the optimal condition with evaporating temperature of 62 • C, the improvement degree is 2.5% with the pump efficiency of 60% in the ORC. In a practical application, the pump efficiency is usually about 60%. An appropriate efficiency of 80% is used for the flow efficiency in the high vertical pipes. As shown in Figure 5, it is obviously that the specific energy is improved by using GDOPC in considered conditions. The lower the pump efficiency in the ORC, the higher the improvement degree. The advantage becomes more and more obvious with increasing evaporating temperature. At the optimal condition with evaporating temperature of 62°C, the improvement degree is 2.5% with the pump efficiency of 60% in the ORC. In a practical application, the pump efficiency is usually about 60%. An appropriate efficiency of 80% is used for the flow efficiency in the high vertical pipes. In considered conditions, the higher the evaporating temperature, the higher the improvement degree. However, a higher fluid liquid column is needed in the conditions with higher evaporating temperature, as shown in Figure 6. For example, in the condition with heat source temperature of 150 °C, evaporating temperature of 130 °C and working fluid of R245fa, the improvement degree of specific energy and cycle thermal efficiency can reach 7.1% and 6.8% while the fluid liquid column In considered conditions, the higher the evaporating temperature, the higher the improvement degree. However, a higher fluid liquid column is needed in the conditions with higher evaporating temperature, as shown in Figure 6. For example, in the condition with heat source temperature of 150 • C, evaporating temperature of 130 • C and working fluid of R245fa, the improvement degree of specific energy and cycle thermal efficiency can reach 7.1% and 6.8% while the fluid liquid column reaches 164.6 m. That means the system must be higher under higher evaporating temperature. Then, the construction difficulty and the system cost increase with increasing the system height. reaches 164.6 m. That means the system must be higher under higher evaporating temperature. Then, the construction difficulty and the system cost increase with increasing the system height.

Analysis on Pressure Drop and Heat Transfer of GDOPC
There are two important factors that impact the performance of GDOPC, namely, heat transfer in heat exchangers and resistance pressure drop in the system. Total pressure difference in the system is produced by the fluid liquid column and is proportional to its height. In the cycle, the total pressure difference is consumed by turbine and flow resistance in the system. The local resistance in the system can be expressed as ( ) . The on-way resistance coefficient λ is determined by the pipe roughness and the Reynolds number according to the flow pattern. Therefore, the local resistance is proportional to 2 u while the on-way resistance is not in some flow patterns. The rough turbulence zone is selected in this analysis because the on-way resistance coefficient λ is only determined by the pipe roughness in this flow pattern. Then in the analysis, not only the local resistance, but also the on-way resistance, are proportional to 2 u . Hence the total pressure difference can be expressed as, There are three heat transfer types in the heat exchangers in the system, such as single phase heat transfer, boiling heat transfer and condensation heat transfer. Many factors, namely, mass flow rate, flow area, flow pattern, heat flux density and temperature difference impact on the boiling heat transfer and the condensation heat transfer. Single phase heat transfer, by contrast, is simpler, as shown in Equation (10), if the shell and tube heat exchanger is selected and the working fluid flows in the tube pass. The coefficient n is specified as 0.4 when the fluid is heated and is specified as 0.3 when the fluid is cooled. From Equation (10), it can be concluded that the heat transfer coefficient is proportional to 8 . 0 u . The higher the velocity, the higher the heat transfer coefficient is. In this analysis, single phase heat transfer is selected to study the variation of heat transfer coefficient.

Analysis on Pressure Drop and Heat Transfer of GDOPC
There are two important factors that impact the performance of GDOPC, namely, heat transfer in heat exchangers and resistance pressure drop in the system. Total pressure difference in the system is produced by the fluid liquid column and is proportional to its height. In the cycle, the total pressure difference is consumed by turbine and flow resistance in the system. The local resistance in the system can be expressed as ς(ρu 2 /2) and the on-way resistance can be expressed as (λl/d) · (ρu 2 /2). The on-way resistance coefficient λ is determined by the pipe roughness and the Reynolds number according to the flow pattern. Therefore, the local resistance is proportional to u 2 while the on-way resistance is not in some flow patterns. The rough turbulence zone is selected in this analysis because the on-way resistance coefficient λ is only determined by the pipe roughness in this flow pattern. Then in the analysis, not only the local resistance, but also the on-way resistance, are proportional to u 2 . Hence the total pressure difference can be expressed as, There are three heat transfer types in the heat exchangers in the system, such as single phase heat transfer, boiling heat transfer and condensation heat transfer. Many factors, namely, mass flow rate, flow area, flow pattern, heat flux density and temperature difference impact on the boiling heat transfer and the condensation heat transfer. Single phase heat transfer, by contrast, is simpler, as shown in Equation (10), if the shell and tube heat exchanger is selected and the working fluid flows in the tube pass. The coefficient n is specified as 0.4 when the fluid is heated and is specified as 0.3 when the fluid is cooled. From Equation (10), it can be concluded that the heat transfer coefficient is proportional to u 0.8 . The higher the velocity, the higher the heat transfer coefficient is. In this analysis, single phase heat transfer is selected to study the variation of heat transfer coefficient.
For R245fa, the optimal condition with evaporating temperature of 62 • C is selected as the standard condition to analyze the pressure drop and the heat transfer. In the standard condition, the turbine pressure drop is 0.2366 MPa and the resistance pressure drop is 0.0025 MPa. The resistance pressure drop is as much as 1.05% of the total pressure difference. Figure 7 shows the variations of cycle performance, namely, fluid velocity and heat transfer coefficient with the increase of resistance pressure drop. As shown in the figure, velocity in the pipes increases with the increase of resistance pressure drop, as well as the heat transfer coefficient. The resistance pressure drop is proportional to u 2 , which is stated above. Therefore, the velocity increases with increasing the resistance pressure drop ratio. From the above analysis about heat transfer, heat transfer coefficient is proportional to u 0.8 . Therefore, the heat transfer coefficient also increases with increasing resistance pressure drop ratio. For R245fa, the optimal condition with evaporating temperature of 62°C is selected as the standard condition to analyze the pressure drop and the heat transfer. In the standard condition, the turbine pressure drop is 0.2366 MPa and the resistance pressure drop is 0.0025 MPa. The resistance pressure drop is as much as 1.05% of the total pressure difference. Figure 7 shows the variations of cycle performance, namely, fluid velocity and heat transfer coefficient with the increase of resistance pressure drop. As shown in the figure, velocity in the pipes increases with the increase of resistance pressure drop, as well as the heat transfer coefficient. The resistance pressure drop is proportional to 2 u , which is stated above. Therefore, the velocity increases with increasing the resistance pressure drop ratio. From the above analysis about heat transfer, heat transfer coefficient is proportional to

The Performance of Several Organic Fluids for GDOPC
The performance of 15 organic fluids for GDOPC is investigated in this section. The specified conditions are the same as in section 3.1, as shown in Table 2. As shown in Figure 8, each fluid gives different specific energy from each other, especially in the range from 60 °C to 66 °C, where each fluid has a maximum specific energy. In considered conditions, R1234yf and R227ea give higher specific energy than the other fluids while cyclopentane gives the lowest specific energy.

The Performance of Several Organic Fluids for GDOPC
The performance of 15 organic fluids for GDOPC is investigated in this section. The specified conditions are the same as in Section 3.1, as shown in Table 2. As shown in Figure 8, each fluid gives different specific energy from each other, especially in the range from 60 • C to 66 • C, where each fluid has a maximum specific energy. In considered conditions, R1234yf and R227ea give higher specific energy than the other fluids while cyclopentane gives the lowest specific energy. Specific energy is determined by cycle thermal efficiency and hot water outlet temperature, which is explained in section 3.1. In considered conditions, deviations for different fluids exist for not only cycle thermal efficiency but also hot water outlet temperature. Deviation of hot water outlet temperature tends to zero with increasing evaporating temperature, as shown in Figure 9, while deviation of cycle thermal efficiency tends to zero with decreasing evaporating temperature, which is indicated in Figure 10. Therefore, the maximum deviations of specific energy for different fluids appear in considered conditions. The deviation tends to zero with increasing or decreasing evaporating temperature. Specific energy is determined by cycle thermal efficiency and hot water outlet temperature, which is explained in Section 3.1. In considered conditions, deviations for different fluids exist for not only cycle thermal efficiency but also hot water outlet temperature. Deviation of hot water outlet temperature tends to zero with increasing evaporating temperature, as shown in Figure 9, while deviation of cycle thermal efficiency tends to zero with decreasing evaporating temperature, which is indicated in Figure 10. Therefore, the maximum deviations of specific energy for different fluids appear in considered conditions. The deviation tends to zero with increasing or decreasing evaporating temperature. Specific energy is determined by cycle thermal efficiency and hot water outlet temperature, which is explained in section 3.1. In considered conditions, deviations for different fluids exist for not only cycle thermal efficiency but also hot water outlet temperature. Deviation of hot water outlet temperature tends to zero with increasing evaporating temperature, as shown in Figure 9, while deviation of cycle thermal efficiency tends to zero with decreasing evaporating temperature, which is indicated in Figure 10. Therefore, the maximum deviations of specific energy for different fluids appear in considered conditions. The deviation tends to zero with increasing or decreasing evaporating temperature. As shown in Figure 10, for each considered fluid, cycle thermal efficiency increases with increasing evaporating temperature. In considered conditions, cyclopentane gives the highest cycle thermal efficiency while R227ea and RC318 give lower cycle thermal efficiency. Cycle thermal efficiency depends on average absorbing heat temperature in the evaporator and average releasing heat temperature in the condenser. Cycle thermal efficiency increases with increasing absorbing heat temperature, as well as decreasing releasing heat temperature. In the considered conditions, average absorbing temperature increases with increasing evaporating temperature while average releasing temperature keeps constant in general. Therefore, thermal efficiency increases with the increase of evaporating temperature, which is similar to that of ORC. As shown in Figure 10, for each considered fluid, cycle thermal efficiency increases with increasing evaporating temperature. In considered conditions, cyclopentane gives the highest cycle thermal efficiency while R227ea and RC318 give lower cycle thermal efficiency. Cycle thermal efficiency depends on average absorbing heat temperature in the evaporator and average releasing heat temperature in the condenser. Cycle thermal efficiency increases with increasing absorbing heat temperature, as well as decreasing releasing heat temperature. In the considered conditions, average absorbing temperature increases with increasing evaporating temperature while average releasing temperature keeps constant in general. Therefore, thermal efficiency increases with the increase of evaporating temperature, which is similar to that of ORC. The optimal condition for each fluid is defined as the condition in which specific energy appears. Figure 11 shows specific energy, cycle thermal efficiency and liquid column height in the optimal condition for each fluid. Among considered fluids, R1234yf and R227ea give higher specific energy than the others and the values are 4.84 kJ/kg and 4.82 kJ/kg, respectively. Their cycle thermal efficiencies are 4.09% and 4.03%, respectively. However, they need high fluid liquid columns and the values are 76.55 m and 45.65 m, respectively. Although R1234yf and R227ea have good performance for GDOPC, the disadvantage of high construction difficulty and system cost limits the use of the two fluids. R365mfc and cyclopentane do not have excellent specific energy and cycle thermal efficiency, but they give very low liquid column height. In the optimal condition, R365mfc gives a specific energy of 4.46 kJ/kg and a cycle thermal efficiency of 3.93%. Its liquid column height is as low as 9.04 m. For cyclopentane, a liquid column height of 10.88 m is needed while the corresponding specific energy and cycle thermal efficiency are 4.35 kJ/kg and 4.06%, respectively. Therefore, they have potential to be used with high evaporating temperature for GDOPC. The optimal condition for each fluid is defined as the condition in which specific energy appears. Figure 11 shows specific energy, cycle thermal efficiency and liquid column height in the optimal condition for each fluid. Among considered fluids, R1234yf and R227ea give higher specific energy than the others and the values are 4.84 kJ/kg and 4.82 kJ/kg, respectively. Their cycle thermal efficiencies are 4.09% and 4.03%, respectively. However, they need high fluid liquid columns and the values are 76.55 m and 45.65 m, respectively. Although R1234yf and R227ea have good performance for GDOPC, the disadvantage of high construction difficulty and system cost limits the use of the two fluids. R365mfc and cyclopentane do not have excellent specific energy and cycle thermal efficiency, but they give very low liquid column height. In the optimal condition, R365mfc gives a specific energy of 4.46 kJ/kg and a cycle thermal efficiency of 3.93%. Its liquid column height is as low as 9.04 m. For cyclopentane, a liquid column height of 10.88 m is needed while the corresponding specific energy and cycle thermal efficiency are 4.35 kJ/kg and 4.06%, respectively. Therefore, they have potential to be used with high evaporating temperature for GDOPC.  Figure 12 shows three parameters, namely, difference between evaporating pressure (saturated pressure at 62°C) and condensing pressure (saturated pressure at 40°C), saturated liquid density at condensing temperature (40°C), and liquid column height with evaporating temperature of 62°C. The liquid column height is correlated with pressure difference and saturated liquid density at the condensing temperature. For fluids from number 1 to number 10, the saturated liquid density are similar, so the difference of liquid column height is caused generally by the pressure difference. The spider diagrams of pressure difference and liquid column height are similar for fluids from number 1 to number 10. For fluids from number 11 to number 15, the saturated liquid densities are very different from each other and the spider diagram shape of the liquid column height is determined by the pressure difference and the saturated liquid density.
It can be noted that the low liquid column height of R365mfc benefits from a low pressure difference while cyclopentane also benefits from high density. Therefore, except for good cycle performance, low saturated pressure (or difference between evaporating pressure and condensing pressure) and high density should also be determined as the principle for selecting fluid for GDOPC.  Figure 12 shows three parameters, namely, difference between evaporating pressure (saturated pressure at 62 • C) and condensing pressure (saturated pressure at 40 • C), saturated liquid density at condensing temperature (40 • C), and liquid column height with evaporating temperature of 62 • C. The liquid column height is correlated with pressure difference and saturated liquid density at the condensing temperature. For fluids from number 1 to number 10, the saturated liquid density are similar, so the difference of liquid column height is caused generally by the pressure difference. The spider diagrams of pressure difference and liquid column height are similar for fluids from number 1 to number 10. For fluids from number 11 to number 15, the saturated liquid densities are very different from each other and the spider diagram shape of the liquid column height is determined by the pressure difference and the saturated liquid density.  Figure 12 shows three parameters, namely, difference between evaporating pressure (saturated pressure at 62°C) and condensing pressure (saturated pressure at 40°C), saturated liquid density at condensing temperature (40°C), and liquid column height with evaporating temperature of 62°C. The liquid column height is correlated with pressure difference and saturated liquid density at the condensing temperature. For fluids from number 1 to number 10, the saturated liquid density are similar, so the difference of liquid column height is caused generally by the pressure difference. The spider diagrams of pressure difference and liquid column height are similar for fluids from number 1 to number 10. For fluids from number 11 to number 15, the saturated liquid densities are very different from each other and the spider diagram shape of the liquid column height is determined by the pressure difference and the saturated liquid density.
It can be noted that the low liquid column height of R365mfc benefits from a low pressure difference while cyclopentane also benefits from high density. Therefore, except for good cycle performance, low saturated pressure (or difference between evaporating pressure and condensing pressure) and high density should also be determined as the principle for selecting fluid for GDOPC. It can be noted that the low liquid column height of R365mfc benefits from a low pressure difference while cyclopentane also benefits from high density. Therefore, except for good cycle