The design and optimization of Organic Rankine Cycle (ORC) require knowledge concerning the thermophysical properties of the working fluids: pure components or mixtures. These properties are generally calculated by thermodynamic and transport property models (thermodynamic or equation of state or correlations). The parameters of these models are adjusted on accurate experimental data. The main experimental data of interest concern phase equilibrium properties (noncritical and critical data), volumetric properties (density and speed of sound), energetic properties (enthalpy, heat capacity), and transport properties (dynamic viscosity and thermal conductivity). In this chapter, some experimental techniques frequently used to obtain the experimental data are presented. Also, we will present some models frequently used to correlate the data and some results (comparison between experimental data and model predictions).
- working fluid design
- experimental techniques
- transport properties
- thermodynamic properties
- equations of state
The utilization of energy available at low, average, and high temperature can be one solution to reduce the energy consumption particularly the fossil energy and to reduce emission of CO2 in the atmosphere. Closed power cycle (Brayton cycle and ORC) proves to be the solution to convert heat sources into energy. Heat source can come from geothermal, solar, and biomass energy or from processes and energy systems. In general, the conditions of the heat source are fixed, and a temperature glide may be observed. In consequence, it is necessary to design the most suitable cycle architecture and to select the best working fluid in order to obtain the best performance and design of each component of the system. The selection of fluid requires thermodynamic models, and these models currently need experimental data in order to optimize their parameters. The knowledge of experimental techniques is very important to measure the thermophysical properties and to estimate their experimental uncertainties. The choice of the most suitable technique depends on the type of data to be determined, the range of pressure and temperature, the precision required, and the composition of the mixture if necessary. In the literature several studies published concerning the investigation of working fluid for ORC applications exist. We will not present and described all the working fluid investigation. In 1985, Badr et al.  have examined around 68 pure working fluids (including natural fluid and hydrofluorocarbons) and given the main characteristic of the working fluids and their impacts on heat exchangers and turbine design. We can notice that in 1985, chlorofluorocarbon were not definitively banned. More recently, Saleh et al.  made a screening of 31 pure working fluids. They concluded that the thermal efficiency is better if the critical temperature of the working fluid is higher as possible. To increase the performance of heat exchangers and so reduce their size, Maizza and Maizza  suggest to use fluid with high density and high heat of vaporization. In 2012, Liu et al. , in the context of power generation, presented a two-stage Rankine cycle for electricity power plants. Ten pure working fluids from different chemical families (aromatics, hydrofluoroolefin, hydrofluorocarbons, hydrocarbons, ammonia) were tested. They concluded that the performance was not affected by the fluid selected. They concluded that selection would be realized on installation volume, size of the different components in the cycle, environmental protection, and operator safety. In 2017, Rahbar et al.  published a complete review of ORC for small-scale applications. In their paper, they present some main characteristic of the working fluid.
2. Description of Organic Rankine Cycle
An Organic Rankine Cycle is composed of a boiler/evaporator, a condenser, a pump, and a turbine (expander). The main application of ORC concerns the transformation/utilization of heat source at low or intermediate temperature (around 80°C). Figure 1 reminds the schematic diagram of ORC. The pump compresses the liquid state working fluid until its desired pressure (and so temperature). The liquid is heated and vaporized in the boiler/evaporator which is also called the generator of vapor (the heat source). The vapor state working fluid is expanded in the turbine. During this operation electricity can be produced thanks to a generator. At low pressure, the working fluid is cooled in the condenser.
The cycle performance depends on the architecture of the system but also on the chosen working fluid and the operating conditions. The Figure 2 presents a typical T-s diagram of an ORC with a heating and cooling medium.
3. Working fluid selection: characteristics
The choice of a working fluid is very crucial in the ORC systems. The working fluid should guarantee the cycle performance and system efficiency and must be adapted to the operating conditions of the system (temperature, pressure). Also hygiene, safety, and environment (HSE) aspects of the fluid have to be taken into account: the working fluid should satisfy the environment protection standards and ensure the safety of the operators. Moreover, the thermophysical properties of the working fluids impact strongly on component’s size.
As the efficiency of the cycle is better for fluids with high critical temperature value , the working fluids are often heavy compounds with important molecular weight and low boiling point. Three types of organic fluids exist: dry, isentropic, and wet fluids. Figure 3 presents the T-S diagrams of these fluids. For dry fluid, dT/dS > 0, and for isentropic fluid, dT/dS = 0. It signifies that during expansion, there is no formation of liquid droplet. For wet fluid, dT/dS < 0 (like water). During expansion, liquid droplets are formed and can damage the equipment. So, it can be necessary to add a superheating equipment.
In general, fluorinated components are used as working fluids. It is important to note that in 2014, European F-gas directive plans the prohibition of fluorinated working fluids with GWP of 2500 or more from 2020. The 2009 F-gas regulation fixes the limits of GWP for each year (Bolaji ). In 2018, the objective was to use fluids with GWP < 1300 close to the GWP of R134A (GWP = 1300). In 2020, the objective is to use fluids with GWP < 1000. In order to reach the objectives in terms of GWP, two solutions exist: the first one consists of the development of new fluids with low GWP values, such as hydrofluoroolefins (HFOs). The second one consists in developing new blends of refrigerants, less than four components . With the existing equipment retrofit aspects are also very important. For all the cases, it is important to consider the thermophysical properties for the selection of the organic working fluid.
3.1 Chemical properties
3.2 Thermodynamic properties
The other thermodynamic properties are also useful but not crucial for the selection for the working fluid (heat capacity, surface tension).
3.3 Transport properties
4. Experimental techniques for the estimation of thermophysical properties
In this section, we will present several experimental techniques used for the experimental determination of thermodynamic and transport properties. Experimental methods for the investigation of thermophysical properties belong either to closed or open circuit methods.
4.1 Equilibrium properties and critical point
Synthetic or analytic methods can be considered for the determination of equilibrium properties. Vapor-liquid equilibrium properties can be obtained using the “static-analytic” method. Herein, the mixture is enclosed in an equilibrium cell equipped with a mixing mechanism to get fast equilibrium conditions. When the equilibrium is reached, small quantities of the phases are sampled and analyzed through chromatographic analyzers. A complete description of the setup is available in Wang et al.’s  paper. The apparatus is similar to the one present on Figure 5. Capillary samplers like ROLSI™ (Armines’s patent) can be used to take samples.
The variable volume cell technique (Figure 6) can be cited as a static-synthetic method [13, 14]. The components of the mixture are introduced separately, and the composition is known by weighing procedure or after analysis. The volume of the cell is modified with a piston to study bubble points. At fixed temperature, saturating properties (pressure and saturated molar volume) of the mixture are determined through the pressure vs. volume curve recorded that displays a break point.
Isochoric method (Figure 7) can be used also to measure the dew point of multicomponent system. The components of the mixture are introduced separately, and the composition is known by weighing procedure or after analysis. The mixture is introduced at its vapor state, and then the temperature slowly decreases. The pressure and temperature are recorded. When the first drop of liquid appears, a break in the P-T curve is observed. The break point corresponds to the dew point. This technique is identical to the technique used to determine gas hydrate dissociation points .
Concerning critical point measurement, it is necessary to use a special device. The technique is based on dynamic-synthetic method where the mixture is circulating through the equilibrium cell (Figure 8) under specific conditions of temperature and pressure. A critical point can be determined by visually observing the critical opalescence and the simultaneous disappearance and reappearance of the meniscus, i.e., of the liquid-vapor interface from the middle of the view cell.
4.2 Volumetric properties or densities
It is well known that density is required for the development of equations of state and the development of models for mass and heat transfers. Several techniques which can be used to measure the density exist. We can cite the hydrostatic balance densitometer coupling with magnetic suspension (single-sinker method from Wagner et al. ), density measurement with vibrating bodies, bellows , and isochoric method . More details concerning these techniques of measurements (and others) are detailed in the IUPAC book dedicated to experimental thermodynamics . Herein, we will only describe the technique based on vibrating tube densitometer.
A mixture with known composition can circulate through a vibrating U-tube (static or dynamic mode). Density is deduced from careful calibration using reference fluids. This apparatus can be used to obtain (PρT) data of compressed phases. A complete description of this technique is available in the papers of Coquelet et al. . This technique is not very recommended close to the critical point. In effect, vibration of the tube may provoke a phase transition. Other techniques based on isochoric method can be used also to determine the volumetric properties at the vicinity of the critical point  (Figure 9).
4.3 Speed of sound
The speed of sound data is also very important to determine the equation of state, and it is linked to other thermodynamic properties. In effect, the isothermal compressibility κT is related to the isentropic compressibility κS via Maxwell’s relations (Eq. (1)):
In Eq. (2), v [m3/mol] is the molar volume of the compound, Cp is the heat capacity [J/mol], and κS is the isentropic compressibility [Pa−1] and is determined thanks to the measurements of speed of sound c [m/s] and density ρ [kg/m3] reached in this work. Indeed, these three properties are linked in the liquid phase via . Like with density measurement, several techniques which were developed to obtain the value of speed of sound exist. The spherical resonator developed by Mehl and Moldover  (with high accuracy in gases), Trusler and Zarari , and Benedetto et al.  can be cited as example. For liquid and dense fluid, pulse-echo techniques are preferred for measuring the speed of sound particularly at high pressure. More details concerning these techniques of measurements (and others) are detailed in the IUPAC book dedicated to experimental thermodynamics . Herein, we will describe one technique used at Heriot-Watt University . Figure 10 describes the cell of measurement. A cylindrical acoustic cell with well-known dimension is considered to measure the sound speed in the fluid using through-transmission method of ultrasonic testing. In this method, a transducer is located on one side of the cell, and a detector is placed on the opposite side of the acoustic cell (electric signal is converted into ultrasound waves and vice versa). An oscilloscope is used to observe the waves. Speed of sound is obtained by dividing the period of the waves by the distance between the speed of sound transducer and detector.
4.4 Heat capacity
The determination of isobaric heat capacity is done using a differential scanning calorimeter (DSC). The equipment is composed of two cells: one is the measurement cell, and one is the reference cell. A sample is introduced into the measurement cell, and a temperature ramp is applied. Knowing the heat flux transferred (absorbed or released, ) and the ramp, it is possible to estimate the heat capacity (Eq. (2)):
In Eq. (4),
4.5 Transport properties
4.5.1 Interest of transport properties in engineering
Chemical and mechanical engineers use correlations with nondimensional numbers (Eq. (3)) to estimate heat transfer coefficient (
In Eq. (1), Nu is the Nusselt number , Re is the Reynolds number , and Pr is the Prandtl number with h being the heat transfer coefficient [J.m−2], L a characteristic length [m], λ the thermal conductivity [J.m−1], ρ the density [kg.m3], μ the viscosity [Pa.s], Cp the heat capacity [J.mol−1.K−1], and v the speed of the fluid [m.s−1]. The correlation can be modified with the utilization of Weber number (, γ is the surface tension) taking into account the effect of interfacial tension during the formation of the bubble of gas in the evaporator or drop of liquids in the condenser.
We remind that for a heat exchanger, global heat transfer coefficient (U) depends on the local heat transfers (h) of the two fluids and the thermal conductivity of the material of the heat exchanger.
4.5.2 Dynamic viscosity
Like density measurements, several techniques to determine viscosity of fluids exist. For example, we can cite the falling ball technique , capillary technique, vibration quartz , and vibrating bodies [29, 30]. We invite the reader to have a look in the paper from Le Neindre  for a complete overview of the techniques of measurement. Herein we will only describe the viscometer used to measure dynamic viscosities using the capillary tube viscosity method. A schematic view of the setup used at Heriot-Watt University is shown in Figure 12. The apparatus is comprised of two small cylinders connected to each other through a capillary tube and a temperature-dependent calibrated internal diameter. A complete description of this technique is available in the paper of Kashefi et al. . The temperature of the system was set to the desired condition, and the desired pressure was set using the hand pump. Poiseuille equation (Eq. (4)) (in laminar flow conditions) can relate the pressure drop across the capillary tube to the viscosity, tube characteristics, and also flow rate for laminar flow:
In Eq. (3), ΔP is the differential pressure across the capillary tube viscometer in psi, Q represents the flow rate in cm3/sec, L is the length of the capillary tube in cm, D refers to the internal diameter of the capillary tube in cm, μ is the dynamic viscosity of the flown fluid in cP, and C is the unit conversion factor equal to 6,894,757 if the above units are used.
4.5.3 Thermal conductivity
Several techniques exist to measure the thermal conductivity of fluids . Among those methods, the transient hot-wire method is considered to be the most used technique and to be a very accurate and reliable technique to measure this thermophysical property. The basic theory of the transient hot-wire method is presented in Healy’s paper , and procedure of measurement is fully described in the paper of Marsh et al. . Generally, a transient hot-wire apparatus consists of one highly pure platinum wire, a current source, a voltage meter, a data acquisition system, and a computer (Figure 13). The current source provides a constant current for the platinum wire, which is embedded in the tested fluid. Then the temperature of the wire will rise because of the Joule effect. Subsequently, the temperature of the fluid will also change as a result of the heat conduction between the wire and the fluid (heat radiation and convection are in general neglected). The temperature rise of the platinum wire as a function of the duration t is given by Eq. (5):
In Eq. (5), λ is the thermal conductivity of the sample, q is the constant power provided to the wire, W/m. κ is the thermal diffusivity of the fluid, r is the radius of the hot wire, and C is Euler’s constant, whose value is 1.781. As the relation between the ΔT and ln(t) can be determined through the experiments, the thermal conductivity of the tested fluid can be calculated using this equation. Compared with other methods, the convenience, accuracy, and the short duration of transient hot-wire method make it a widely used method nowadays.
5. Data treatment
Process or system simulator required models or correlation in order to define the best operating conditions (T, P, flow, compositions) with the maximum of efficiency and the best coefficient of performance (COP) but also to design each component (heat exchangers, expanders, valves, expander). In this section we will just present few models which can be used with the experimental data in order to predict the phase diagrams and thermodynamic and transport properties.
5.1 Thermodynamic models
In this section, we propose a presentation of three types of thermodynamic models. These models are briefly described. The main difference concerns the number of parameters we have to adjust and so the quantity and type of experimental data we have to acquire in order to adjust the parameters.
The thermodynamic properties are obtained by using equations of state (EoSs). In the process simulators, the most famous EoSs are of cubic type, such as Eq. (6):
In Eq. (6), R is the ideal gas constant, a and b are the parameters of the EoS calculated using the critical temperature (Tc) and pressure (Pc) of each component, and α(T) is a function of temperature, acentric factor, Tc, and Pc. u and w are other parameters.
Another type of equation of state of the molecular type can be used. The Helmholtz energy is calculated by considering all the molecular interactions like dispersion, polarity, H bonding (association), etc. Equation (7) describes the method of calculation of the Helmholtz energy A:
In Eq. (7), kb is the Boltzmann constant, T is the temperature, and N is the mole number. The most known molecular EoSs of this type are SAFT type. Based on the Wertheim’s statistical theory of associative fluids (1984), Chapman et al.  developed the first EoS SAFT (statistical associating fluid theory) called SAFT-0. Many versions exist today, such as SAFT-VR , PC-SAFT , and SAFT Mie . The various versions differ mainly in the choice of the reference fluid, the radial distribution function, and explicit expressions of the terms of perturbation.
The last type of equation of state which can be used to estimate the thermodynamic properties is based on multi-fluid approximation. It is well known that from the knowledge of Helmholtz energy, all thermodynamic properties can be calculated. This equation of state is explained in terms of reduced molar Helmholtz free energy (Eq. (8)). Temperature and density are expressed in the dimensionless variables and :
In Eq. (8), the exponent “id” stands for the ideal gas contribution, and exponent “res” is the residual contribution. and are the reduced density and temperature, respectively. Concerning the development of equation of state for the mixtures, the first possibility is to consider mixing rules for each parameter like in the cubic equation of state. The best approach is to consider the multi-fluid approximation. This approach was introduced by Tillner-Roth in 1993 . It applies mixing rules to the Helmholtz free energy of the mixture of components (Eq. (9)):
In Eq. (9), superscript “E” is for excess properties, and subscripts “p,” “q,” and “j” are the component index. An excess property is defined to calculate the deviation from ideal mixture. is called departure function from the ideal solution. It is an empirical function, like Eq. (10), fitted to experimental binary mixture data, mainly densities, speed of sound, or heat of mixing. In this departure function, the Fpq parameters take into account the behavior of one binary pair with another. If only vapor-liquid equilibrium properties are available, can be considered to be equal to zero :
Equation (10) contains several adjustable parameters αi, αk, tk, dk, and lk. Experimental data is required to adjust these parameters. Moreover, if lk = 0, then γ = 0 and if lk ≠ 0, then γ =1.
With the multi-fluid approximation, it is important to calculate the new critical properties corresponding to the mixture studied (superscript “mix”) because reduced parameters are used in Eq. (8).
For example, we can use and where and are adjustable binary interaction parameters. VLE data should be used to fit parameters, and if experimental data concerning densities of mixture are available, it is possible to fit .
5.2 Transport property models
Concerning the transport properties, different approaches exist. One consists in using the corresponding state method. The most famous approach is the TRAPP method developed by the NIST. Huber et al.  and Klein et al.  have developed a series of equations adapted to the prediction of viscosities and thermal conductivities of pure components and mixtures. The approach consists in modifying the transport properties in the ideal dilute gas state taking into account the molecular interactions (and so the density of the fluid with temperature and pressure).
The viscosity of a dilute gas can be determined as a function of temperature by Eq. (11). A dilute gas is composed by noninteractive rigid spheres of diameter σ:
In Eq. (11), T is temperature in K, M is the molar mass in g/mol, σ is the diameter of the rigid sphere in nm, and C is a constant depending of the molecule considered. Chapman-Enskog cited in Poling and Prausnitz  have introduced an additional term called “collision integral” which takes into account the collision between the molecules (Eq. (12)):
In Eq. (12), , and 𝜖/k is an empirical factor link to the potential of interaction.
Concerning thermal conductivity, λ0(T) represents the dilute gas thermal conductivity and can be calculated by Eq. (13):
When the pressure and so the density of the fluid increases, molecular interaction between molecules cannot be neglected, and so it is necessary to apply a correction (like a residual term for equation of state). Equations (14) and (15) lead to calculate dynamic viscosity and thermal conductivity:
In Eq. (14), is the residual viscosity:
In Eq. (15), Δλr (ρ, T) is the residual thermal conductivity, and ΔλC(ρ, T) is the empirical critical enhancement.
The R134a is the reference fluid for the refrigerants , but for the hydrocarbons, it is better to consider propane or methane . The viscosity and thermal conductivity of the other fluids are calculated from the properties of reference fluids. Critical properties of the fluid are required. Equation (16) presents the equation for the dynamic viscosity:
6. Presentation of some results
In this section, we will present some results obtained for pure component (transport properties) and multicomponent systems (equilibrium properties) on working fluids already published.
6.1 Phase diagram
We will present the results obtained for three binary systems. The first one concerns a mixture of CO2 and R1234ze(e) published by Wang et al. ). Figure 14 presents the phase diagram. The system presents no azeotropic behavior. We can notice that with the experimental data, we can also predict a critical point using asymptotic laws of behavior . Figure 15 presents a comparison with the critical point measured by Juntaratchat et al.  using a critical point setup similar to the equipment already presented in Section 4.
The binary system R245fa + isopentane is presented on Figure 16 at 392.87 K. Measurement was done using static-analytic method . We have used REFPROP 10.0 to correlate the data (REFPROP 10.0  uses Fundamental Helmholtz equation). We can observe a good agreement between REFPROP prediction and the experimental data (we have one comment concerning this point: the reference of the data used in the data treatment by REFPROP is not clearly mentioned).
6.2 Bubble point
Variable volume cell was used to determine bubble pressure of the ternary system composed with R32 + R290 + R227ea . Figure 17 presents the results obtained. The data were correlated by a cubic equation of state (Redlich-Kwong-Soave EoS , MHV1  mixing rules, and NRTL  activity coefficient model).
Vibrating tube densitometer technique was used to measure the densities of the R1216 . The results are presented on Figure 18. The same technique is used to obtain density data of R1234yf (Figure 19) at saturation . Density at the vicinity of the critical point was obtained by an isochoric method used by Tanaka and Higashi .
6.4 Dynamic viscosity
Capillary viscometer was used by Laesecke et al.  to determine the dynamic viscosity of liquid R245fa at saturation. Figure 20 presents the results and comparison with prediction using REFPROP 10.0. A good agreement is observed.
6.5 Thermal conductivity
Marrucho et al.  have used the transient hot-wire technique to measure the thermal conductivity of R365mfc. Figure 21 presents the results obtained at 336.85 and 377.4 K for several pressures from 1 to 4.5 bar. Comparison was done with prediction from REFPROP 10.0. A good agreement is observed.
The optimization of ORC depends strongly on the capability of models to predict the thermophysical properties of working fluids for their selection and retrofit in existing ORC equipment. The main thermophysical properties include phase equilibria (and so critical point), densities, speed of sound, dynamic viscosity, and thermal conductivity. Through this book chapter, the reader can easily understand that several experimental techniques developed to measure the thermophysical properties exist. Some of them were presented and described. In general, these techniques are very accurate, and we recommend to use the experimental technique with which the operator/engineer is the most familiar with. It is obvious to see that we have the capability to determine all the thermophysical properties of interest for the design and optimization of ORC. In our opinion, there are two main challenges for the future. The first one concerns the development of equipment which permits the determination of several thermophysical properties together and particularly thermophysical properties at equilibrium. The second one concerns the development of equipment which require a small quantity of sample and so the utilization of in situ analysis techniques. In effect, it is difficult to have important quantity of new synthesis molecules (with high purity) to realize the measurement of accurate thermophysical properties of pure components and mixtures. The experimental data will be used to adjust parameters on models or correlations. Data acquisition is an essential step for a good working fluid design. In the literature an optimal selection (Figure 22) of working fluids as a function of the temperature of the heating medium exists . This classification does not concern mixtures, and in the case of the utilization of mixtures, experimental characterization has to be done.