In Situ Determination of Solid Fraction from the Measured Hydrate Slurry Flow Rate and Pressure Drop across Oriﬁce

: Two-phase ﬂow is encountered in various engineering areas, including the pharmaceutical, chemical, and food industries, desalination facilities, and thermal energy storage systems. Cost-e ﬀ ective and non-invasive monitoring of the solid volume fraction, which governs the thermos-physical properties of two-phase medium, is important for ﬂow assurance. The ﬂow loop having an inner diameter of 21.5 mm and length of about 12.2 m was equipped with square-edged oriﬁce and slash plate pump. Tetraﬂuroethane (R134a) hydrate slurry of the speciﬁed solid volume fraction could be formed within the ﬂow loop by removing an appropriate amount of water, and simultaneously injecting the pertinent amount of R134a while chilled at 275 K. The uncertainty in the thus-obtained solid volume fraction was smaller than 9%, with the largest contribution originating from the uncertain hydration number. The near power-law relationship between the oriﬁce pressure loss coe ﬃ cient and Metzner–Reed Reynolds number was recognized. However, the nonlinear nature of the Reynolds number with respect to the solid volume fraction inevitably makes the solution procedure iterative. The short span pressure di ﬀ erences across the oriﬁce were regressed to yield empirical correlation, with which the solid volume fraction of R134a slurry could be determined from the measured pressure drop across the oriﬁce and the ﬂow rate. The uncertainty was less than 12% of the thus determined solid volume fraction.


Introduction
Over the past decade, the use of phase change material (PCM) as thermal energy storage (TES) medium has received significant research attention. Due to the larger equivalent heat capacity compared to single phase medium and better heat transfer characteristics compared to gas phase, slurries are best suited as a heat transfer fluid and TES medium [1][2][3][4][5][6]. Among numerous phase change materials, ice slurries are one of the commonly available phase change material. Ice slurries are recognized as suitable materials for TES and they have already been used in many practical processes [7,8].
Interest in gas hydrates technology has also grown drastically, because of its diverse disciplines, including the potential utilization as effective phase change material for TES medium [9][10][11]. Clathrate hydrates are ice-like crystalline inclusion compounds that are formed even at above zero temperature through the hydrogen bonding of water molecules, with the assistance of gas molecules if the pressure is high enough [12,13]. Hydrate slurry is a relatively new type of phase change medium, in which the solid phase is clathrate hydrate particles, and the liquid phase is liquid water, which serves as the

Experiments
The formation and circulation of aqueous tetraflouroethane (R134a) hydrate slurry were made to occur within a flow loop having an inner diameter of 21.5 mm. Figure 1 shows a schematic of the experimental apparatus. The flow loop consisted of a rectangular service flow loop section, and extended test flow loop section.
The volumetric flow rate was measured using an electromagnetic flowmeter (KTM-800), whose measuring range is a maximum 10 m/s with an uncertainty of ±0.5% for the flow velocity larger than 0.3 m/s. The correlation coefficient between the pump rotation speed and volume flow rate was 99.99% for the experimental conditions of the present study. The flow velocity can be estimated from the pump rotation speed with the same uncertainty as the flow meter after calibration. Flow visualization was enabled using sight glass assemblies installed upstream and downstream of the pump. The test flow loop section was connected to the first downstream corner from the pump of the service flow loop and was about 8200 mm long. The orifice test segment was located 1.5 m downstream from the 180-degree bend of the test flow loop. For the non-Newtonian flow it is important to know required pipe length for the fluid to become fully developed as given by Poole and Ridely et al. [51]. In this study the calculated entry length of straight pipe was smaller than 0.97 m for the test conditions employed in the discussion. The orifice test segment was a glass tube having a concentric square-edged orifice in the middle. The thickness of the orifice plate was 6 mm, and the bore diameter of the orifice was either (10 or 15) mm. SENSUS differential pressure transducers having an uncertainty of 0.025% up to 200 kPa were employed to measure the differences in pressure at 50 mm upstream (PO1) and 100 mm downstream (PO2) from the upstream face of the orifice. The pressure at 1000 mm downstream of the orifice (PO3), as well as the representative pressures of the service flow loop (PS) and test flow loop (PT), were measured with SENSUS pressure transducers having an uncertainty of 0.4%, and maximum range of 1000 kPa.
The service flow loop was immersed in an acrylic water bath. The test flow loop was immersed in a stainless steel water bath, which had a chilling water supply tube, and a chilling water overflow drain tube. The temperatures of the water baths were maintained close to 275 K by circulating the chilled water from constant temperature baths, both with the flow rate of 28 LPM. The R134a storage tank was also immersed in the water bath, whose temperature was maintained close to 295 K, to avoid condensation. The temperatures of the flow loop and (TS1, TS2, TT), chilling water baths, and R134a storage tank were measured with T-type thermocouples having an uncertainty of ±0.3 K after The service flow loop section was about 4000 mm long, and rigged with flow driving pump, bubble removal chamber, water supply/drain port, volumetric flow meter, R134a gas supply port, sight glass assembly, and tabs to measure the service loop pressure (P S ) and temperatures (T S1 and T S2 ). The flow is driven by the swash plate pump manufactured by Parker (Model F12-090-MF-IV-D), which had a displacement volume of 90 cm 3 per revolution. The rotation speed of the pump driving motor was controlled by a Rexroth (EFC 3600) inverter. The bubble removal chamber was a commercially available 1 L size sight glass assembly. The top of the bubble removal chamber was closed by a blank plate with a vent valve, and the bottom was welded to the T-type pipe fitting. A pneumatic control valve was installed between the gas supply port and the R134a storage tank to control the flow loop pressure close to the desired value engaged with a pressure transducer for (P S ). The volumetric flow rate was measured using an electromagnetic flowmeter (KTM-800), whose measuring range is a maximum 10 m/s with an uncertainty of ±0.5% for the flow velocity larger than 0.3 m/s. The correlation coefficient between the pump rotation speed and volume flow rate was 99.99% for the experimental conditions of the present study. The flow velocity can be estimated from the pump rotation speed with the same uncertainty as the flow meter after calibration. Flow visualization was enabled using sight glass assemblies installed upstream and downstream of the pump.
The test flow loop section was connected to the first downstream corner from the pump of the service flow loop and was about 8200 mm long. The orifice test segment was located 1.5 m downstream from the 180-degree bend of the test flow loop. For the non-Newtonian flow it is important to know required pipe length for the fluid to become fully developed as given by Poole and Ridely et al. [51]. In this study the calculated entry length of straight pipe was smaller than 0.97 m for the test conditions employed in the discussion. The orifice test segment was a glass tube having a concentric square-edged orifice in the middle. The thickness of the orifice plate was 6 mm, and the bore diameter of the orifice was either (10 or 15) mm. SENSUS differential pressure transducers having an uncertainty of 0.025% up to 200 kPa were employed to measure the differences in pressure at 50 mm upstream (P O1 ) and 100 mm downstream (P O2 ) from the upstream face of the orifice. The pressure at 1000 mm downstream of the orifice (P O3 ), as well as the representative pressures of the service flow loop (P S ) and test flow loop (P T ), were measured with SENSUS pressure transducers having an uncertainty of 0.4%, and maximum range of 1000 kPa.
The service flow loop was immersed in an acrylic water bath. The test flow loop was immersed in a stainless steel water bath, which had a chilling water supply tube, and a chilling water overflow drain tube. The temperatures of the water baths were maintained close to 275 K by circulating the chilled water from constant temperature baths, both with the flow rate of 28 LPM. The R134a storage tank was also immersed in the water bath, whose temperature was maintained close to 295 K, to avoid condensation. The temperatures of the flow loop and (T S1 , T S2 , T T ), chilling water baths, and R134a storage tank were measured with T-type thermocouples having an uncertainty of ±0.3 K after calibration. Pressure and temperature data were processed and recorded by data logger (Agilent, 34970A) and software program HP-VEE.
R134a hydrate is a nonstoichiometric inclusion solid compound formed from guest gas molecules and host water molecules under sufficiently high-pressure condition at a given temperature. The hydration number of gas hydrate is defined as the average number of water molecules per unit cell that contains one gas molecule. As the measured average hydration number (r) of R134a hydrate is 17 [52], the fraction of water molecule (x H 2 O ) is 73.86 wt.%, and that of R134a molecule (x R ) is 26.14 wt.%, considering that the molecular weight of water (M H 2 O ) is 18.02 kg/kmol, and that of R134a (M R ) is 102.03 kg/kmol: The rheological characteristics of two-phase flow strongly depend on the solid volume fraction, especially in the laminar regime [39].
In order to obtain the valid empirical correlation for the pressure difference across the orifice, the solid volume fraction needs to be precisely acquired or evaluated. In the present work, hydrate slurry having the prespecified solid volume fraction of (15, 20, 25, 28, or 30)% was formed within the flow loop. Considering the difference in liquid water density and hydrate density [33,34] shown in Table 1, the net amount of water to be removed from the flow loop ( H 2 O ,remove) and the amount of R134a to be injected (M R ) to attain specific hydrate solid fractions (φ s ) can be calculated as follows: where, V T ( = 5.68 L) is the total volume of the flow loop. Table 2 presents the mass of water that needs to be removed from the flow loop, and the mass of the R134a that needs to be injected into the flow loop for each specific hydrate volume fraction. In order to prepare for an experimental run, the flow loop was completely filled with water. At the same time, air bubbles were eliminated through the bubble removal chamber, while the pump was running with a volume flow rate of 15.8 LPM, which is equivalent to a flow velocity of 0.73 m/s. When the removal of air bubble was complete, the pump was paused, and the valves connecting the air removing chamber to the flow loop were closed. Then, the gas supply port was opened to pressurize the flow loop to help push the water out. The calculated amount of water in Table 2 was withdrawn. During the intermittent manual withdrawal of liquid water, the loop pressure fluctuated as shown in Figure 2 The peak pressures were less than 400 kPa, due to the dissolution of R134a gas. After water removal, the gas supply port was closed, and the pump operation was resumed to check the stabilized loop pressure. The flow loop pressure after 40 s of stabilization was about 120 kPa. When the temperature of H 2 O and R134a binary system is 275 K, hydrate formation is driven at pressures larger than 79.6 kPa (the three-phase equilibrium pressure among water, vapor R134a, and hydrate). If the pressure is larger than 314.6 kPa (saturation pressure of R134a), the condensation of injected gas occurs simultaneously [53]. The duration of stabilization was determined taking two opposite influences into account: to allow for the significant dissolution of gas, and to prevent the hydrate formation. Since the flow loop pressure was larger than 79.6 kPa, R134a hydrate can form past the induction time. The stabilized flow loop pressure smaller than 314.6 kPa indicated that the void resulted from water withdrawal was filled with gaseous R134a.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 18 was running with a volume flow rate of 15.8 LPM, which is equivalent to a flow velocity of 0.73 m/s. When the removal of air bubble was complete, the pump was paused, and the valves connecting the air removing chamber to the flow loop were closed. Then, the gas supply port was opened to pressurize the flow loop to help push the water out. The calculated amount of water in Table 2 was withdrawn. During the intermittent manual withdrawal of liquid water, the loop pressure fluctuated as shown in Figure 2 The peak pressures were less than 400 kPa, due to the dissolution of R134a gas. After water removal, the gas supply port was closed, and the pump operation was resumed to check the stabilized loop pressure. The flow loop pressure after 40 s of stabilization was about 120 kPa. When the temperature of H2O and R134a binary system is 275 K, hydrate formation is driven at pressures larger than 79.6 kPa (the three-phase equilibrium pressure among water, vapor R134a, and hydrate). If the pressure is larger than 314.6 kPa (saturation pressure of R134a), the condensation of injected gas occurs simultaneously [53]. The duration of stabilization was determined taking two opposite influences into account: to allow for the significant dissolution of gas, and to prevent the hydrate formation. Since the flow loop pressure was larger than 79.6 kPa, R134a hydrate can form past the induction time. The stabilized flow loop pressure smaller than 314.6 kPa indicated that the void resulted from water withdrawal was filled with gaseous R134a. When the gas supply port was opened again, the flow loop was pressurized with R134a up to 400 kPa, to initiate the formation of hydrate. The gauging of injected gas was also started by marking the pressure and temperature of the gas storage tank as the initial state. The gas supply under near constant pressure condition was continued while the flow of hydrate slurry was qualitatively observed with the help of the sight glass. The gas supply was ceased when the amount of injected gas, which was calculated by the Peng-Robinson equation of state, reached the predetermined value. Figure 2 shows that the flow loop pressure (Ps) decreased afterwards approaching steady state, which indicated the completion of hydrate formation. The hydrate formation took about five min for the solid fraction of 15 vol.%, and about 25 min for the solid fraction of 30 vol.%.
Finally, the pressure drops across the orifice were measured for the steady slurry flow of different flow rates. The pressure measurements for a specific combination of flow rate and solid fraction were made at least five times while the flow rates were increased and decreased, to check for the presence of any hysteresis effect, which was not found. After completing a series of pressure measurements for a given solid fraction, the entire flow loop was emptied and flushed with water, before repeating the experiment for the different hydrate fractions.
1. The uncertainty in the thus-obtained solid volume fraction is mainly attributed to four contributions: When the gas supply port was opened again, the flow loop was pressurized with R134a up to 400 kPa, to initiate the formation of hydrate. The gauging of injected gas was also started by marking the pressure and temperature of the gas storage tank as the initial state. The gas supply under near constant pressure condition was continued while the flow of hydrate slurry was qualitatively observed with the help of the sight glass. The gas supply was ceased when the amount of injected gas, which was calculated by the Peng-Robinson equation of state, reached the predetermined value. Figure 2 shows that the flow loop pressure (Ps) decreased afterwards approaching steady state, which indicated the completion of hydrate formation. The hydrate formation took about five min for the solid fraction of 15 vol.%, and about 25 min for the solid fraction of 30 vol.%.
Finally, the pressure drops across the orifice were measured for the steady slurry flow of different flow rates. The pressure measurements for a specific combination of flow rate and solid fraction were made at least five times while the flow rates were increased and decreased, to check for the presence of any hysteresis effect, which was not found. After completing a series of pressure measurements for a given solid fraction, the entire flow loop was emptied and flushed with water, before repeating the experiment for the different hydrate fractions.
1. The uncertainty in the thus-obtained solid volume fraction is mainly attributed to four contributions: 2. The amount of R134a that filled the void created by the water removal is neglected in the calculation of the required gas injection amount. For example, the void having 6.8% of total volume is created in order to form the 30% hydrate volume fraction slurry. The void at the pressure and temperature of 275 K and 120 kPa (the stabilized pressure in Figure 2) contains the R134a amount equivalent to 0.03% of the required gas injection.
3. The amount of R134a dissolved into liquid water, which does not participate in hydrate formation, is neglected. The solubility of R134a at the pressure and temperature of 275 K and 140 kPa (the steady flow condition of Figure 2) is about 570 ppm, which is equivalent to 1.9% of the required gas injection amount to form the slurry of 10% hydrate volume fraction.
4. The largest contribution originates from the uncertain hydration number. The hydration number can increase 7% due to hydrate clathrate unoccupied with guest molecule. This contribution is equivalent to an uncertainty of 5.2% [Equation (1)].
The contribution caused by the error in measuring the mass of removed water is about 0.5%, the error in measuring the pressure is about 0.4%, etc.
Therefore, it can be concluded that the uncertainty in the solid volume fraction is smaller than 9%. During the pressure drop measurements, the hydrate slurry within the flow loop was brought into near phase equilibrium, even though the pressure and temperature fluctuated within the limited extents. If the amount of water remaining after removal or the amount of injected gas was less than the calculated value, gas bubbles would be enclosed within the flow loop after hydrate formation to balance the insufficient volume, and the flow loop would be at a pressure of three-phase equilibrium among water, vapor R134a, and hydrate (79.6 kPa). In contrast, if the amount of remaining water or the amount of injected gas was more than the calculated value, liquid R134a droplets would be enclosed. The flow loop pressure would be equal to or larger than the pressure of three-phase equilibrium among liquid R134a, vapor R134a, and hydrate (314.6 kPa), owing to the comparable density values of the condensed phases (Table 1). Figure 2 shows that the steady flow loop pressure can exceed 79.6 kPa, due to the inclusion of volatile impurity gases. As the migration of bubbles could not be observed through sight glass, the R134a and impurity gas mixture is believed to be trapped at relatively higher spaces within the flow loop, or to be of negligible amount.

Results and Discussion
The measurements of pressure drop across the orifice were made at different flow velocities, for different solid volume fractions of (15, 20, 25, 28, and 30)%, and for orifice bores of (10 and 15) mm. Figure 3 exhibits the pressure difference (∆P 13 ) between the pressures measured at 50 mm upstream (P O1 ) and 1100 mm downstream (P O3 ) from the upstream face of the orifice. For a given solid fraction, ∆P 13 increases almost linearly with the increase of flow velocity for both orifice bores of (10 and 15) mm.  The net pressure loss (ΔP) caused by the orifice is calculated by subtracting from ΔP13 the pressure loss (ΔPLd) caused by the straight pipe having the same length between pressure difference measuring locations: 13 Ld where, LLd is the pipe length, and D is the pipe diameter. The pressure loss coefficient (K) is calculated according to the common definition: The net pressure loss (∆P) caused by the orifice is calculated by subtracting from ∆P 13 the pressure loss (∆P Ld ) caused by the straight pipe having the same length between pressure difference measuring locations: The empirical equation reported by Zabd et al. [39] is utilized to evaluate ∆P Ld for the pipe length of 1150 mm: where, L Ld is the pipe length, and D is the pipe diameter. The pressure loss coefficient (K) is calculated according to the common definition: The Metzner-Reed Reynolds number (Re MR ), frequently employed in non-Newtonian fluid analysis [54], is introduced: where,ρ HS is the density of hydrate slurry. Empirical correlations for the consistency index k and behavior index n of Equation (8) were also available [39], which revealed the pseudoplastic non-Newtonian behavior of aqueous R134a hydrate slurry: The exponent m is a constant having the same value for both orifice bores. A model equation for the multiplier C in Equation (11) can be further constructed as a function of s φ , which may contain different sets of coefficient estimates for the respective orifice bore.   The exponent m is a constant having the same value for both orifice bores. A model equation for the multiplier C in Equation (11) can be further constructed as a function of φ s , which may contain different sets of coefficient estimates for the respective orifice bore.
However, the highly nonlinear nature of Equation (11), especially Re MR with respect to φ s [Equations (8)-(10)], causes a serious disadvantage when used for in situ determination of the solid volume fraction. The solution procedure of Equation (11) for φ s from the measured pair of net pressure loss and flow velocity data is inevitably iterative. In order to determine the solid volume fraction of slurry without iteration, the net pressure loss (∆P) is regressed as the simplest polynomial function of φ s : The results of regression shown in Figure 5 imply that Equation (12) well characterizes the linear relation between ∆P and φ s for the given pair of U and orifice bore (d), despite the average value of slope m in Equation (11) being −0.616, which is noticeably different from −1, as shown in Figure 4. The minimum R2 value is larger than (0.99 and 0.97), when the orifice bore is (10 and 15) mm, respectively. The linear relation between ∆P and φ s is better validated for the smaller orifice bore. The results of regression shown in Figure 5 imply that Equation (12) well characterizes the linear relation between ΔP and s φ for the given pair of U and orifice bore (d), despite the average value of slope m in Equation (11) being −0.616, which is noticeably different from −1, as shown in Figure 4. The minimum R2 value is larger than (0.99 and 0.97), when the orifice bore is (10 and 15) mm, respectively. The linear relation between ΔP and s φ is better validated for the smaller orifice bore.   (12) for the 15 mm orifice bore, using the net pressure loss.
Empirical correlations to determine the solid volume fraction of slurry from the measured values of net pressure loss and velocity are developed separately for the (10 and 15) mm orifice bores from Equation (12) (12) for the 15 mm orifice bore, using the net pressure loss.
Empirical correlations to determine the solid volume fraction of slurry from the measured values of net pressure loss and velocity are developed separately for the (10 and 15) mm orifice bores from Equation (12) Figure 9 presents the reproduced solid volume fractions using Equation (14) for the 15 mm bore orifice. The maximum error is 8.45 vol.%, and the standard error is 2.6 vol.%. The relatively larger errors with the 15 mm bore orifice originate from the smaller R 2 values ( Figure 5). Figure 9. Reproduced solid volume fraction using Equation (14) for the 15 mm bore orifice. Figure 9. Reproduced solid volume fraction using Equation (14) for the 15 mm bore orifice.
The application of Equations (13) and (14) still requires iterative solution procedure, because the evaluation of ∆P Ld using Equation (6) needs a priori knowledge of φ s . Moreover, the pressure measurement location at far downstream distance of 1100 mm may not be acceptable for the compact desalination facility design. Figure 10 presents the short span pressure difference (∆P 12 ) between the pressures measured at 50 mm upstream (P O1 ) and 100 mm downstream (P O2 ) from the upstream face of the orifice.
The ∆P 12 data were processed in the same way as the net pressure loss (∆P), to yield the following empirical correlations; The application of Equations (13) and (14) still requires iterative solution procedure, because the evaluation of ΔPLd using Equation (6) needs a priori knowledge of s φ . Moreover, the pressure measurement location at far downstream distance of 1100 mm may not be acceptable for the compact desalination facility design. Figure 10 presents the short span pressure difference (ΔP12) between the pressures measured at 50 mm upstream (PO1) and 100 mm downstream (PO2) from the upstream face of the orifice.
(a) The ΔP12 data were processed in the same way as the net pressure loss (ΔP), to yield the following empirical correlations;  Figure 11 presents the reproduced solid volume fractions ( ' s φ ) using Equation (15). The maximum error is 2.33 vol.%, and the standard error is 0.68 vol.% (less than 4.6% of the determined solid volume fraction). Figure 11. Reproduced solid volume fractions using Equation (15) for the 10 mm bore orifice. Figure 12 presents the reproduced solid volume fractions using Equation (16). The maximum error is 6.45 vol.%, and the standard error is 2.5 vol.%. The performance of Equation (15) is not inferior to that of Equation (13).  Figure 11 presents the reproduced solid volume fractions (φ s ) using Equation (15). The maximum error is 2.33 vol.%, and the standard error is 0.68 vol.% (less than 4.6% of the determined solid volume fraction).
Appl. Sci. 2020, 10, x FOR PEER REVIEW 14 of 18 Figure 11 presents the reproduced solid volume fractions ( ' s φ ) using Equation (15). The maximum error is 2.33 vol.%, and the standard error is 0.68 vol.% (less than 4.6% of the determined solid volume fraction). Figure 11. Reproduced solid volume fractions using Equation (15) for the 10 mm bore orifice. Figure 12 presents the reproduced solid volume fractions using Equation (16). The maximum error is 6.45 vol.%, and the standard error is 2.5 vol.%. The performance of Equation (15) is not inferior to that of Equation (13).  Figure 12 presents the reproduced solid volume fractions using Equation (16). The maximum error is 6.45 vol.%, and the standard error is 2.5 vol.%. The performance of Equation (15) is not inferior to that of Equation (13). Figure 11. Reproduced solid volume fractions using Equation (15) for the 10 mm bore orifice. Figure 12 presents the reproduced solid volume fractions using Equation (16). The maximum error is 6.45 vol.%, and the standard error is 2.5 vol.%. The performance of Equation (15) is not inferior to that of Equation (13). It is concluded that the solid volume fraction of aqueous R134a hydrate slurry can be monitored using Equation (15) from the measured flow rate and pressure difference between 50 mm upstream and 100 mm downstream from the upstream face of the square-edged orifice having a bore of 10 mm. The uncertainty is less than 12% of the determined solid volume fraction, in which about 7.2% is the biased uncertainty in obtaining raw solid volume fraction data [contributions (1)-(3) to the uncertainty in the solid volume fraction, as discussed in the "Experiments"], while about 4.6% is the random error resulting from the regression process, in addition to contribution (4). A short span between pressure measuring locations is advisable to ease the complexity of facility layout design and cost constraints. Comparing to the other technologies, this method is non-destructive, cost effective and does not require frequent calibration. It is concluded that the solid volume fraction of aqueous R134a hydrate slurry can be monitored using Equation (15) from the measured flow rate and pressure difference between 50 mm upstream and 100 mm downstream from the upstream face of the square-edged orifice having a bore of 10 mm. The uncertainty is less than 12% of the determined solid volume fraction, in which about 7.2% is the biased uncertainty in obtaining raw solid volume fraction data [contributions (1)-(3) to the uncertainty in the solid volume fraction, as discussed in the "Experiments"], while about 4.6% is the random error resulting from the regression process, in addition to contribution (4). A short span between pressure measuring locations is advisable to ease the complexity of facility layout design and cost constraints. Comparing to the other technologies, this method is non-destructive, cost effective and does not require frequent calibration.

Conclusions
The empirical correlations were developed for in situ determination of the solid volume fraction from the measured pressure drop across the orifice and flow rate of R134a hydrate slurry. The conclusions drawn from the experimental results are as follows: 1. R134a hydrate slurry of the specified solid volume fraction can be formed within the flow loop chilled at 275 K by removing an appropriate amount of water, and simultaneously injecting the pertinent amount of R134a. The uncertainty in the thus-obtained solid volume fraction is smaller than 9%, with the largest contribution originating from the uncertain hydration number.
2. The velocity of slurry flow can be estimated from the rotation speed of the slash plate pump with the same uncertainty as employed for the volume flow meter after calibration.
The near linear relationship between ln(K) and ln(Re MR ) was recognized with the average value of slope m to be −0.616. However, the nonlinear nature Re MR with respect to φ s inevitably makes the solution procedure for φ s iterative. The evaluation of pressure loss (∆P Ld ) caused by the straight pipe also requires a priori knowledge of φ s .
3. The short span pressure difference (∆P 12 ) between the pressures measured at 50 mm upstream and 100 mm downstream of the 10 mm bore orifice were regressed as linear functions of φ s to yield the empirical correlation, Equation (15), which showed excellent reproducibility.
4. The uncertainty is less than 12% of the thus-determined solid volume fraction, which includes about 7.2% of biased uncertainty in obtaining the raw solid volume fraction data, and about 4.6% of random error.

Conflicts of Interest:
The authors declare no conflict of interest.