International Journal of Chemical Engineering

Volume 2012 (2012), Article ID 856120, 8 pages

http://dx.doi.org/10.1155/2012/856120

## Effect of Impeller Agitation on Preparation of Tetra-*n*-Butyl Ammonium Bromide Semiclathrate Hydrate Slurries

Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan

Received 18 October 2011; Revised 28 December 2011; Accepted 15 January 2012

Academic Editor: Naoto Ohmura

Copyright © 2012 Shunsuke Hashimoto et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The slurries-containing tetra-*n*-butyl ammonium bromide (TBAB) solution and its semiclathrate hydrate have attracted a lot of interest as latent heat transport media. These hydrate slurries contain some microparticles of crystal, and the size and shape of these hydrate particles could affect the mobility of slurries. Hence, it is essential to investigate the efficient hydrate-slurry preparation methods and the effect of hydrate particles on the fluid property of slurries for the application to latent heat transport media. In the present study, the effect of agitation on particle size distribution and aggregation of particles was studied to prepare easily flowing TBAB hydrate slurries that were suitable for fluid transport. First of all, the effects of impeller rotational speed and impeller type on the particle size and frequency of aggregation were investigated. The results suggested that the particle size distribution and the frequency of particle aggregation are strongly affected by the intensity of shear rate and its uniformity, which was controllable with impeller type and its rotation speed.

#### 1. Introduction

Recently, a large-scale district (regional) cooling system has attracted attention of many researchers because it works on saving energy and global sustainability due to its high performance of utilization of exhausted heat from large plants or factories [1]. In such system, the pumping energy loss due to the transportation of heat media is not small. To solve this problem, a high-density heat transportation system involving slurries with phase change materials (hereafter, PCM) having latent heat has been developed, for example, ice/water slurries [2]. The utilization of such PCM slurries can reduce flow rate of heat media, since latent heat of some materials is much larger than sensible heat. Fukushima et al. [3] reported the ability of slurries-containing tetra-*n*-butyl ammonium bromide (hereafter, TBAB) semiclathrate hydrate and its solution as latent-heat transport media, which was more favorable than ice/water system because they can be operated at relatively high temperature (*ca*. 285 K) and atmospheric pressure. Unlike ordinary gas hydrates, in these semiclathrate hydrates, the quaternary ammonium cation and anion are incorporated with the hydrogen bonds of water molecules to construct the hydrate cage [4]. In addition, four butyl groups are one-by-one encaged with a large cage separately in this semiclathrate system. There are various reports about the crystal structure and the role of bromine for the TBAB hydrates [4, 5]. The hydration numbers of these structures are 26 and 38, which are described as Type A (tetragonal) and Type B (orthorhombic), respectively [5, 6]. The concentration of aqueous TBAB solutions results in differences of the crystal structure and thermodynamic stability of TBAB hydrate. TBAB semiclathrate hydrates have become the subject of much attention as an attractive medium for not only refrigerant but also H_{2} storage [7–9].

There are a few reports about thermodynamic properties such as phase equilibria and dissociation enthalpy of TBAB hydrate [10] and about rheological property of these hydrate slurries [11] for the application to latent heat transport media. These hydrate slurries contain some microparticles of crystal, and the size and shape of these hydrate particles may affect the mobility of slurries. In addition, with the view of heat conductivity, it is necessary to make hydrate particles smaller and subglobular and to control their agglomerating property so that hydrate particles dissociate rapidly and uniformly. Therefore, it is essential to investigate the efficient preparation method of hydrate slurries and the effect of hydrate particles on the fluid property of slurries for the application to latent heat transport media. There are some reference papers that investigate the effect of impeller agitation and bottom shape of stirred tank on the particle-size distribution and/or dispersion efficiency in the suspension systems containing polymer particles [12–14]. There, the dispersion efficiency such as the agitating speed for complete floating of solid particles was mainly investigated. In addition, some unique and complicated techniques such as multistage impeller system and unsteady agitation were usually utilized for the dispersion of solid particles. However, there are few findings that discuss the basic relation between the particle size and/or shape and the intensity and configuration of shear field. Furthermore, there is no report about the control of particle property for clathrate hydrates by means of impeller agitation.

Generally, impeller agitation has been adopted as the method for efficient crystallization of hydrate particles. However, there are few reports about the control of particle property for clathrate hydrates by means of impeller agitation. In the present study, we aim at searching for the favorable condition of agitation that minimizes the total power needed for both agitation and pumping of the prepared hydrate slurries and the effect of agitation on particle-size distribution and aggregation of particles to prepare TBAB hydrate slurries that easily flow and are suitable for fluid transport. In this paper, the effects of impeller rotational speed and impeller type on the particle size and frequency of aggregation are investigated. In addition, the key factor to control particle characteristics in TBAB hydrate slurries is briefly discussed.

#### 2. Experimental Section

##### 2.1. Materials

Research grade TBAB (mole fraction purity 0.980) was obtained from Wako Pure Chemical Industries, Ltd. TBAB was used without further purifications. In addition, deionized water was produced using water-manufacturing equipment made by Nihon Millipore K. K.

##### 2.2. Apparatus and Procedures

The experimental apparatus in the present study consisted of the following parts: acrylic cylindrical vessel (diameter, : 120 mm, height, : 200 mm, handcrafted), three types of agitating impellers (handcrafted), agitating motors (HEIDON, BLh-300), torque meter (type: YT, Shintou Kagaku, Co., Ltd.), temperature controller (TAITEC, CL-80R), and CCD camera (KEYENCE, VH5910). Figure 1 shows the details of three types of agitating impellers used in the present study: two-bladed paddle (2-BP), Maxblend (Sumitomo Heavy Industries, Ltd., MB), and rotational cylinder (RC). The 2-BP impeller was located at a height of 60 mm. The diameter of shaft was 8 mm. Four baffles (width: 12 mm) were attached in the vessel with 2-BP and MB, while no baffle was attached in the case of RC because uniform flow field was created. In addition, the characteristics of three impellers are summarized in Table 1 (details for the velocity field of three types of impeller are shown in the former part of Section 3). The average shear rate, can be calculated by use of Metzner and Otto’s equation [15] as follows: where (depending on the type of impeller) and stand for Metzner and Otto’s constant and the rotation speed of impeller, respectively. In this case, was equivalent to 10 for 2-BP and 20 for MB, respectively, which was estimated and correlated systematically by Kamei et al. [16]. In the case of RC (simple shear field), can be generally estimated by the following equation: where stands for the diameter ratio of cylinder and vessel (in the present study, ). The values are also listed in Table 1 with value.

Figure 2 shows the thermodynamic stability of TBAB hydrate (temperature, -concentration, diagram) at atmospheric pressure. The TBAB hydrate is stable at the temperatures below the line that is equivalent to the stability boundary of TBAB hydrate. As shown in Figure 2, the stoichiometric concentration of Type A TBAB hydrate is 40.5 wt% where the TBAB hydrate can exist stably at *ca*. 285 K. In the present study, the hydrate fraction in slurries, and was determined by the following equation:
where , , and stand for the initial concentration of TBAB in aqueous solution, steady concentration of TBAB in the liquid phase of slurries, and stoichiometric concentration of Type A TBAB hydrate, respectively. In general, decreases monotonically as temperature rises. In the present study, 20 wt% of TBAB in aqueous solution was adopted as a working fluid.

Firstly, the working fluid (20 wt% solution of TBAB) was poured into the vessel up to 120 mm (liquid height, ). Then, the fluid was cooled down up to 280.6 K (subcooling degree is *ca*. 0.7 K and then . 0.1) and the rotation of impeller was started. After the system reached the cyclostationary state, some seed crystals of Type A TBAB hydrate, which were prepared in advance and annealed at 280.6 K, were injected into the vessel, and then crystallization occurred. Raman spectroscopy preliminarily revealed that only Type A TBAB hydrate formed under the present experimental conditions [8]. In this case, it was confirmed that the particles of TBAB hydrate reached the steady state after an hour. Hence, appropriate amount of hydrate particles was removed up from the mother fluid and observed by use of CCD camera on the petri dish that was chilled with Peltier device. Typical picture of hydrate particles was shown in Figure 3. In the present study, 500 particles in total were used for particle size analyses of each experimental run. The particle diameter was defined as the major axis of rod-like hydrate particle. In addition, the aggregation particle was defined as the particle that consisted of more than two particles. The density of Type A TBAB hydrate particle was measured by means of falling ball method. The density of hydrate particle was comparable with the theoretical value, which could be calculated from the lattice constant of unit cell for tetragonal Type A TBAB hydrate. Additionally, it was slightly larger than that of aqueous solution. In this case, upward flows above a certain level by external force such as impeller rotation are required to float solid particles completely, which is well known as a criterion of impeller rotational speed for the complete floating of solid particles [17]. In the present study, impeller rotational speed in all experiments was more than this threshold value. Hence, it is reasonable to decide that the hydrate particles float completely and well dispersed in the vessel.

The required power for mixing, *,* was obtained from rotational torque, , of impeller shaft by the following equation:

#### 3. Results and Discussion

Prior to experiments, the flow velocity vector patterns of three types of impeller were obtained from the numerical simulation by the use of commercial CFD software, Rflow (RFLOW, Co., Ltd.). Figure 4 shows the flow velocity vector patterns of (a) 2-BP (with 4 baffles), (b) MB (with 4 baffles), and (c) RC (without baffles). The density and viscosity of fluid were set equal to those of TBAB aqueous solution (1080 kg/m^{3} and 0.003 Pa·s, resp.), and general model was adopted as a turbulent flow model. In calculations, the was set at 100 rpm. As it can be seen in Figure 4, 2-BP generates nonuniform velocity field in both vertical and horizontal directions, where flow velocity is locally high near the impeller and low near the wall. On the other hand, in the case of RC, uniform circulating flow (a large vortex) is vertically generated, and, in particular, uniform rotational flow is circumferentially created. In the case of MB, although there are locally high and/or low flow regions, overall relatively uniform circulating flow is vertically generated. At rpm, the flow is somewhat intensified, while the flow pattern remains the same in all impellers.

*The Experimental Data*

Particle size and distribution, aggregation frequency of particles, and required power for the agitation of slurries obtained in the present study are summarized in Table 2. Hereafter, the detailed discussion about the difference in results for the three impellers is given in each section.

##### 3.1. Particle Size and Its Distribution

Figure 5 shows the particle diameter of TBAB hydrates and its distribution for three types of impeller. In the case with 2-BP and RC, the particle diameter and its distribution do not change with the increase in shear rate. On the other hand, they change drastically as shear rate increases in the case with MB. This indicates that the growth of hydrate particle is suppressed considerably under the shear conditions above a certain level. In addition, even though the average shear rate of RC is lowest among three impellers at the equal speed of impeller rotation, the particle diameter and its distribution are relatively small in RC. Hence, the uniformity of shear rate is one of dominant factors for the control of particle diameter and its distribution. In fact, the 2-BP that generates nonuniform flow field exhibits broad size distribution characteristics of hydrate particles. Although the shear rate of MB is inhomogeneous compared with that of RC, the average shear rate of MB at rpm is relatively very high (Table 1). Consequently, it is suggested that the local region of low shear rate can be kept away in the MB due to its large area of blade and particle diameter would become small.

##### 3.2. Aggregation Frequency of Hydrate Particle

The open symbols in Figure 6 stand for the aggregation frequency of hydrate particles. The aggregation frequency becomes smaller in the order of RC, 2-BP, and MB. In addition, it becomes small as the shear rate increases. This indicates that the aggregation of hydrate particles is suppressed by the collision between hydrate particles and impeller blades and/or baffles, and the presence of field of locally high shear rate.

##### 3.3. Required Power for the Agitation of Hydrate Slurries

The closed symbols in Figure 6 represent the required power to agitate hydrate slurries. As mentioned previously, MB can create the flow field that has a distinct advantage for the control of hydrate particles and limit their aggregation. As expected, however, the required power for agitation in MB is extremely large compared to the other impellers, which is due to the relatively large area of its impeller blade. Incidentally, it is experimentally confirmed that the required power for the agitation of hydrate slurries is almost same as (slightly larger than) that of aqueous solution under the present conditions of hydrate fraction, where the difference between them is at most *ca*. 0.002 W.

##### 3.4. Conceptual Impeller for Efficient Agitation of TBAB Hydrate Slurries

Based on the present experimental results, to prepare the hydrate particles that have small diameter and sharp distribution of diameter, three factors are essential. They are (1) uniformity of shear rate, (2) absence of the local region of low shear rate, and (3) presence of field that has locally high shear rate (or the collision between hydrate particles and impeller blade). Herein, a conceptual impeller that can satisfy above key factors was developed. Figure 7 shows the schematic illustration of this conceptual impeller, velocity vector of its pattern of flow, and its characteristics of generated hydrate particle. As shown in Figure 7(a), this conceptual impeller consists of 2-bladed paddle blade (bottom) and cylinder (top). Figure 7(b) indicates that this conceptual impeller can generate both the region of locally high shear rate (bottom side) and that of uniform shear rate (top side) at the same time. As a result, the uniform and nearly steady circulating flow can be generated. Needless to say, this impeller can prevent the region of low shear rate. In addition, it is expected that the part of impeller blade can prevent the aggregation of hydrate particles. In fact, as shown in Figure 7(c), this conceptual impeller can prepare TBAB hydrate particles that have uniform diameter, and their average diameter is relatively small (0.804 mm). Additionally, the aggregation frequency of TBAB hydrate for this impeller is relatively low (5.6%) and the required power is also low (0.180 W). After all, this conceptual impeller manufactured based on some experimentally obtained dominant factors has good performance for the preparation of TBAB hydrate particles. The findings obtained in the present study are useful for the development of new impeller, which is suitable for the preparation of TBAB hydrate slurries that easily flow and are suitable for fluid transport.

#### 4. Conclusion

To prepare TBAB hydrate slurries that easily flow and are suitable for fluid transport, the effects of impeller rotational speed and impeller type on the particle size and frequency of aggregation have been investigated. The experimental findings are simply summarized in Table 3. For the preparation of hydrate particles that have small diameter and sharp distribution of diameter, three factors and their balances are essential as follows: (1) uniformity of shear rate, (2) absence of the region of low shear rate, and (3) presence of field that has locally high shear rate. Based on these factors, a new conceptual impeller has been proposed, which consists of 2-bladed paddle blade (bottom) and cylinder (top). This conceptual impeller has succeeded in the preparation of smaller and uniform TBAB hydrate particles with relatively low agitation power.

#### Notations

: | Concentration, wt% |

: | Diameter of vessel, m |

: | Height of vessel, m |

: | Height of working fluid, m |

: | Hydrate fraction in slurries |

: | Rotation speed of impeller, rpm |

: | Required power for agitation, W |

: | Temperature, K |

: | Torque, N·m |

: | Diameter ratio of cylinder and vessel |

: | Coordinate of -direction (vertical direction) |

: | Shear rate, s^{−1}. |

#### Acknowledgments

This work was financially supported by Grant-in-Aid for Scientific Research to S. Hashimoto (Young Research, B, 22760581). One of the authors (S. Hashimoto) expresses his special thanks to Dr. H. Oka (RFLOW, Co., Ltd.) for the numerical simulation by use of commercial CFD software, Rflow.

#### References

- S. Kilkis and B. Kilkis, “An analytical, multi-objective an optimization algorithm for energy efficiency in district cooling systems and green buildings,”
*ASHRAE Transactions*, vol. 117, pp. a13/1–a13/8, 2011. View at Google Scholar - S. Fukusako, Y. Kozawa, M. Yamada, and M. Tanino, “Fundamental researches and developments on fine crystallineice slurries,”
*Transactions of the Japan Society of Refrigerating and Air Conditioning Engineers*, vol. 17, no. 4, pp. 413–437, 2000. View at Google Scholar - S. Fukushima, S. Takao, H. Ogoshi et al., “Development of High-density Cold Latent Heat with Clathrate Hydrate,”
*NKK Technical Report*, vol. 166, pp. 65–70, 1999 (Japanese). View at Google Scholar - L. S. Aladko, Y. A. Dyadin, T. V. Rodionova, and I. S. Terekhova, “Clathrate hydrates of tetrabutylammonium and tetraisoamylammonium halides,”
*Journal of Structural Chemistry*, vol. 43, no. 6, pp. 990–994, 2002. View at Publisher · View at Google Scholar · View at Scopus - W. Shimada, T. Ebinuma, H. Oyama et al., “Separation of gas molecule using tetra-
*n*-butyl ammonium bromide semi-clathrate hydrate crystals,”*Japanese Journal of Applied Physics*, vol. 42, no. 2A, pp. L129–L131, 2003. View at Publisher · View at Google Scholar · View at Scopus - W. Shimada, M. Shiro, H. Kondo et al., “Tetra-
*n*-butylammonium bromide-water (1/38),”*Acta Crystallographica Section C*, vol. 61, no. 2, pp. o65–o66, 2005. View at Publisher · View at Google Scholar · View at Scopus - S. Hashimoto, S. Murayama, T. Sugahara, H. Sato, and K. Ohgaki, “Thermodynamic and Raman spectroscopic studies on H
_{2}+*tetrahydrofuran*+*water*and H_{2}+tetra-*n*-butyl ammonium*bromide*+*water*mixtures containing gas hydrates,”*Chemical Engineering Science*, vol. 61, no. 24, pp. 7884–7888, 2006. View at Publisher · View at Google Scholar · View at Scopus - S. Hashimoto, T. Sugahara, M. Moritoki, H. Sato, and K. Ohgaki, “Thermodynamic stability of hydrogen + tetra-
*n*-butyl ammonium bromide mixed gas hydrate in nonstoichiometric aqueous solutions,”*Chemical Engineering Science*, vol. 63, no. 4, pp. 1092–1097, 2008. View at Publisher · View at Google Scholar · View at Scopus - S. Hashimoto, T. Tsuda, K. Ogata, T. Sugahara, Y. Inoue, and K. Ohgaki, “Thermodynamic properties of hydrogen + tetra-
*n*-butyl ammonium bromide semi-clathrate hydrate,”*Journal of Thermodynamics*, vol. 2010, Article ID 170819, 5 pages, 2010. View at Publisher · View at Google Scholar - H. Oyama, W. Shimada, T. Ebinuma et al., “Phase diagram, latent heat, and specific heat of TBAB semiclathrate hydrate crystals,”
*Fluid Phase Equilibria*, vol. 234, no. 1-2, pp. 131–135, 2005. View at Publisher · View at Google Scholar · View at Scopus - M. Darbouret, M. Cournil, and J. M. Herri, “Rheological study of TBAB hydrate slurries as secondary two-phase refrigerants,”
*International Journal of Refrigeration*, vol. 28, no. 5, pp. 663–671, 2005. View at Publisher · View at Google Scholar · View at Scopus - S. I. Kondo, M. Motoda, K. Takahashi, and H. Horiguchi, “The influence of the bottom shape of an agitated vessel stirred by dual impellers on the distribution of solid concentration,”
*Journal of Chemical Engineering of Japan*, vol. 40, no. 8, pp. 617–621, 2007. View at Publisher · View at Google Scholar · View at Scopus - K. Shimizua, K. Takahashi, E. Suzukia, and T. Nomura, “Effect of baffle geometries on crystal size distribution of aluminum potassium sulfate in a seeded batch crystallizer,”
*Journal of Crystal Growth*, vol. 197, no. 4, pp. 921–926, 1999. View at Google Scholar · View at Scopus - A. Barresi and G. Baldi, “Solid dispersion in an agitated vessel,”
*Chemical Engineering Science*, vol. 42, no. 12, pp. 2949–2956, 1987. View at Google Scholar · View at Scopus - A. B. Metzner and R. E. Otto, “Agitation of non-Newtonian fluids,”
*AIChE Journal*, vol. 3, no. 1, pp. 3–10, 1957. View at Publisher · View at Google Scholar - N. Kamei, S. Hiraoka, Y. Kato et al., “Development of power correlations for various impellers in an agitated vessel in laminar region based on a numerical analysis of 2-D flow,”
*Kagaku-kogaku Ronbunshu*, vol. 20, no. 5, pp. 595–603, 1994 (Japanese). View at Google Scholar - T. N. Zwietering, “Suspending of solid particles in liquid by agitators,”
*Chemical Engineering Science*, vol. 8, no. 3-4, pp. 244–253, 1958. View at Google Scholar · View at Scopus