Kinetics of Trifluoromethane Clathrate Hydrate Formation from CHF3

Sep 11, 2017 - We report the formation kinetics of trifluoromethane clathrate hydrate (CH) from less than 75 μm diameter ice particles and CHF3 gas. ...
1 downloads 6 Views 2MB Size
Subscriber access provided by Columbia University Libraries

Article

Kinetics of Trifluoromethane Clathrate Hydrate Formation from CHF Gas and Ice Particles 3

Jaruwan Amtawong, Suvrajit Sengupta, Michael T Nguyen, Nicole C Carrejo, Jin Guo, Everly B. Fleischer, Rachel W. Martin, and Kenneth C. Janda J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.7b08730 • Publication Date (Web): 11 Sep 2017 Downloaded from http://pubs.acs.org on September 17, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1: Schematic of the aaparatus 65x54mm (600 x 600 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. (a) Pressure versus time curves for a reference run and experimental runs at an initial flow rate of 30 mmol/h and four temperatures: 268 K, 263 K, 258 K, and 253 K. (b) Uptake rate versus time. (c) Percent yield versus time. 82x136mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 2 of 41

Page 3 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3. Percent yield versus time curves of experiments at 268 and 253 K, run over 55 hours. 82x49mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. (a) Pressure versus time curves for a reference run and experimental runs at an initial flow rate of 15 mmol/h and four temperatures: 268 K, 263 K, 258 K, and 253 K. (b) Uptake rate versus time, and (c) Percent yield versus time. 82x125mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 4 of 41

Page 5 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5. (a) Pressure versus time curves for a reference run and experimental runs at an initial flow rate of 48 mmol/h and four temperatures: 268 K, 263 K, 258 K, and 253 K. (b) Uptake rate versus time, and (c) Percent yield versus time. 82x125mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. Percent yield curves for experiments at 258 K with three different flow rates: 15, 30, and 48 mmol/h. 82x49mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 6 of 41

Page 7 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 7. Arrhenius plots for (a) 15 mmol/h, (b) 30 mmol/h, and (c) 48 mmol/h flow rate experiments. 82x136mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Graphic 81x43mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 8 of 41

Page 9 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Kinetics of Trifluoromethane Clathrate Hydrate Formation from CHF3 Gas and Ice Particles Jaruwan Amtawonga, Suvrajit Sengupta, Michael T. Nguyenb, Nicole C. Carrejoc, Jin Guo, Everly B. Fleischer, Rachel W. Martin, Kenneth C. Janda* Department of Chemistry, University of California, Irvine, California 92697, United States Email: [email protected]. Tel.: (949) 824-6022.

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 41

Abstract. We report the formation kinetics of trifluoromethane clathrate hydrate (CH) from 75 μm) were transferred under liquid nitrogen to a precooled reaction cell (310 mL) which is then submerged in a pre-cooled bath and evacuated to ~0.05 MPa. The system is equilibrated to a preset temperature. The trifluoromethane is added via a needle valve adjusted to give the desired gas flow rate into the reaction cell. Gas is

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 41

continuously added until the cell pressure reaches the final pressure fixed by the gas regulator. As the pressure in the reaction cell increases, the gas flow rate into the cell slows down because the pressure drop across the needle valve decreases. The actual gas flow into the cell as a function of the pressure drop across the needle valve is calibrated by gas flow into a cell that contains sand instead of ice. To characterize the gas flow rate into the reaction cell in a particular run, we use the amount of gas delivered during the first hour of the experiment. Refrigerant grade trifluoromethane is used as purchased from Airgas. The actual mass of ice particles used in any experiment was determined by measuring the mass of water left in the reaction vessel after decomposing the hydrate sample. All numerical analysis is based on this mass. The Peng-Robinson equation of state, PV  ZnRT

(1)

where, Z is the compressibility factor as a function of temperature and pressure, is used to calculate the number of moles of gas (n) in the cell at any given instant, from the measured pressure (P) and temperature (T), and the known volume (V) of the cell (adjusted for the volume of ice). The uptake rate of trifluoromethane by ice

dnice can be calculated from: dt

dnm dncell dnice   dt dt dt

(2)

where,

dnm dncell is the rate of change of the number of moles in the gas phase, is the flow rate dt dt

the gas phase calculated from the reference run. The percent yield of the experiment is calculated by integrating the uptake rate as a function of time.

ACS Paragon Plus Environment

6

Page 15 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1. Schematic of the apparatus. III. Results and Data Analysis III.A. Temperature Dependence of the Reaction Rate. A typical set of data, including a reference run and four runs at 253, 258, 263, and 268 K, is shown in Figure 2a. These data are for runs in which the initial gas flow rate into the reaction cell is about 30 mmol/h. Because there is no clathrate hydrate formation in the reference run (top blue curve), the pressure rises monotonically and reaches the maximum pressure of 0.46 MPa. This calibration allows us to determine the gas flow rate into the cell as a function of cell pressure. For each experiment with ice, the pressure initially rises monotonically because there is no hydrate formation at low pressure. Once the initiation pressure is reached, more gas is absorbed into the ice (due to hydrate formation) than leaked into the reaction cell, causing a decrease in pressure. We refer to this reaction regime as stage I, in analogy to the shrinking core model.21,22 For the highest temperature (268 K), the cell pressure reaches 0.37 MPa before the reaction starts. For the lowest temperature (253 K), the onset of reaction is at 0.25 MPa, which occurs nearly an hour sooner than for 268 K. Intermediate temperatures yield results between these two limits. For each temperature, upon initiation, the cell pressure drops to a minimum, at which the reaction rate equals

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 41

the flow rate into the cell, and then gradually increases toward the limiting value set by the gauge of the gas source. Empirically, we refer to the region for which the reaction rate is about equal to the gas flow rate as stage II, and the long-term, slower reaction rate as stage III in analogy to the shrinking core model. We note that our definition of the three reaction stages is purely empirical, as the transition between stages in the shrinking core model is somewhat ill defined. For the data set reported in Figure 2, the duration of stage II increases slightly with temperature. The experiments performed here are too time consuming to repeat them often enough to determine accurate error bars for the results. However, each run was repeated two or more times to ensure that the data is reproducible. The reproducibility of the methods has been more completely examined in Refs. 15-17. In the data tables, we report the results to two significant figures, and they are reliable to 5% of the stated value. As described previously,20 the raw data is converted to uptake rate versus time and percent yield versus time, shown in Figures 2b and 2c, respectively. Figure 2b shows that the peak uptake rate is also inversely proportional to temperature, even though it occurs at a lower cell pressure for lower temperatures. For instance, the maximum uptake rate at 253 K is over two times higher than that of the 268 K run, even though the pressures in the cell are 0.25 and 0.37 MPa, respectively. The thermodynamic driving force for the enclathration process might be understood from the degree of undercooling from the ice-hydrate-vapor equilibrium line. At 0.25 MPa and 0.37 MPa the ice-hydrate vapor equilibrium is approximately at 267 K and 272 K respectively.23 Thus, in the 253 K experiment, the degree of undercooling at the initiation of enclathration is about 14 K, while for the 268 K experiment it is about 4 K. The higher uptake rates at lower temperatures also translate to higher percent yields for stages I and II of the reaction: e.g. 40% for 268 K and 50% for 253 K. Our experimental technique is not sensitive enough to measure the percent yield after

ACS Paragon Plus Environment

8

Page 17 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

many days since the errors in integrating the slow reaction rate over long times will be cumulative. A summary of all data from this study is presented in Table 1.

Figure 2. (a) Pressure versus time curves for a reference run and experimental runs at an initial flow rate of 30 mmol/h and four temperatures: 268 K, 263 K, 258 K, and 253 K. (b) Uptake rate versus time. (c) Percent yield versus time.

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 41

Table 1. Kinetic Data Reported for Trifluoromethane Clathrate Hydrate Formation. temperature (K)

particle size (μm)

flow rate (mmol/h)

initiation pressure (MPa)

maximum uptake rate (μmol/s)

percent yielda

stage III rateb (%/104 s)

duration of stages I and II (s/104)

268 K