Thermal Cracking of Endothermic Hydrocarbon Fuel in Regenerative

7 days ago - The chemical heat sink of Endothermic Hydrocarbon Fuels (EHFs) is generally dependent on its thermal cracking in the cooling channel, whi...
1 downloads 4 Views 4MB Size
Subscriber access provided by Kaohsiung Medical University

Fossil Fuels

Thermal Cracking of Endothermic Hydrocarbon Fuel in Regenerative Cooling Channels with Different Geometric Structures Fuqiang Li, Zaizheng Li, Kai Jing, Li Wang, Xiangwen Zhang, and Guozhu Liu Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.8b00531 • Publication Date (Web): 10 May 2018 Downloaded from http://pubs.acs.org on May 16, 2018

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 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 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.

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 44 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

Energy & Fuels

Thermal Cracking of Endothermic Hydrocarbon Fuel in Regenerative Cooling Channels with Different Geometric Structures Fuqiang Lia, Zaizheng Lia, Kai Jingb, Li Wanga, b, Xiangwen Zhanga, b, Guozhu Liua, b* a

: Key Laboratory for Green Chemical Technology of Ministry of Education, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, People’s Republic of China b

: Collaborative Innovation Center of Chemical Science and Engineering, Tianjin University, Tianjin 300072, People’s Republic of China

* Corresponding Author: Guozhu Liu E-mail address: [email protected] Postal address: 92 Weijin Road, Nankai District, Tianjin, P.R.China Tel./fax: +862285356099

ACS Paragon Plus Environment

Energy & Fuels 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

Abstract The chemical heat sink of Endothermic Hydrocarbon Fuels (EHFs) is generally dependent on its thermal cracking in the cooling channel, which is accompanied and limited by the formation of carbon deposit. In this work, HF-1 (A kerosene-based EHF) was electrically heated in the rectangular, square and circular channels with the same cross-section area under 3.5 MPa to study the effect of cooling channel geometric structures on the thermal cracking and carbon deposition behaviors. It was found that on the similar conditions (inlet flow rate of fuel, pressure, outlet temperature), conversions of HF-1 in both rectangular and square channels were slightly higher than that in the circular one with high selectivity to methane but lower selectivities to the primary cracking products (such as 1-hexene and 1-heptene, etc.). In addition, more carbon deposits were formed in the rectangular and square channels, especially around the corners of channels. Based on the CFD simulation, the possible reasons should be ascribed to the difference in the gradient uniformity near the wall of different channels. The higher temperature and lower velocity in the boundary layer of the quadratic channels might cause the thermal cracking slightly severer and the rapid secondary reactions forming carbon deposit.

ACS Paragon Plus Environment

Page 2 of 44

Page 3 of 44 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

Energy & Fuels

1. Introduction Hydrocarbon fuel is one of the few practical coolants on board for the thermal management system of the advanced aircrafts1. Flowing in the cooling channel, the fuel firstly absorbs the heat through temperature rise or phase transition, and then, after the fuel temperature exceeding the initial pyrolysis temperature (Tip), it further absorbs heat via endothermic reactions, such as thermal cracking. Endothermic Hydrocarbon Fuels (EHFs) are used to define the hydrocarbon fuels with available endothermic reactions. However, the use of chemical heat sink is generally accompanied and limited by the formation of carbon deposit. Therefore, it is crucial for the optimal design and better use of EHFs to understand the complex thermal cracking and carbon deposition behavior in the cooling channels. During the past thirty years, much work has been done towards better understanding the thermal cracking and coking (carbon deposition) of EHFs using various model compounds or kerosene fuels. Huang et al.2 measured heat sinks of JP-7, JP-8+100 and n-octane in tubular reactor with resistive heating method and reported an overall heat sink of JP-7 up to 3838 J/g, which was higher than that of JP-8+100. Goel et al.3 numerically studied jet fuel decomposition in tubular reactor using thermal cracking of n-dodecane as a model compound, and found that the outlet bulk temperature and unreacted fuel fraction could be well predicted by an one-step global model. Ward et al.4-5 performed thermal cracking experiments of n-decane and n-dodecane in tubular reactor and proposed proportional product distribution method to describe the products of mildly (conversion 650 °C) grows. The relative gap of conversion is calculated according to Eq. 9. Taking square and circular channels for example, the conversion gap (700 °C) grows to 17.07% at 4 g/s comparing to 10.59% at 2 g/s, and the gas yield gap is 31.17% at 4 g/s from 22.42% at 2 g/s.

∆X=

Xsqaure-Xcircular Xcircular

(9)

The variations of products distributions as a function of flow rates are also obtained from the thermal cracking of HF-1. The selectivities to methane at 1 g/s and 2 g/s of inlet flow rates are presented in Figure 12(a) and (b). Other molecules like ethane, ethylene, propane and propene present no obvious variations among those three geometries no matter what the flow rates are. When the flow rate is at 1 g/s, methane seems to share equal selectivity among those three geometries. While with the flow rate increasing, the value in circular geometry falls behind that of the other two channels. Also, as shown in Figure12 (c-f), the selectivities to 1-hexene and 1-heptene in circular channels grow to higher values than that

ACS Paragon Plus Environment

Page 19 of 44 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

Energy & Fuels

of quadratic channels as flow rate increases. As for 1-octene, 1- nonene and aromatics, notable changes are not observed. The gaps of conversion and selectivity can be also influenced by temperature and velocity distribution. The distributions of temperature and velocity in the outlet cross-section at both 1 g/s and 4g/s, are shown in Figure 13 and Figure 6. From 1 g/s to 4 g/s, temperature of core fluid is decreased and that near the wall is increased, with the average velocity also enhanced. This is because an increase of inlet flow rate leads to higher fluid velocity and heat load, and in result the wall temperature rises. The variations of fluid temperature and velocity along different lines at 1 g/s and 4 g/s are presented in Figure 14 and Figure 7. Comparing rectangular and circular channels, the fluid temperature in the long centerline direction (in the range of 0%-20% and 80%-100%) is eventually beyond that along the diameter of the circular channel from 1 g/s to 4 g/s, meaning that fluid temperature in boundary layer near the short side wall is higher at 4 g/s than 1 g/s. As for the square and circular channels, the fluid temperature along the centerline is initially lower than that along diameter at 1 g/s but reverses this condition as the flow rate rises up to 4 g/s, revealing that boundary layer has much higher temperature. Considering the velocity, for square and circular channels at 1 g/s, they have the same fluid velocity at 24.4% both in the diagonal and diameter direction. While at 4 g/s, the percentage is decreased to 16.0%. And comparing rectangular and circular channels at 1 g/s, the velocities at 20% of the diagonal and diameter are 3.54 m/s and 3.75 m/s indicating the gap is 6%. When flow rate is enhanced to 4g/s, the values are increased to 14.1 m/s and 14.8

ACS Paragon Plus Environment

Energy & Fuels 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

m/s at the exit, with the gap 5%. That is, the boundary layer with low velocity in quadratic channels is narrowed compared with that of circular geometry, which is not beneficial for widening the gap of conversion. As mentioned in the introduction, the conversion is co-decided by temperature and velocity and the former one takes the leading role. In our results, though the boundary layer of low velocity in quadratic channels is narrowed compared with circular geometry, the fluid temperature especially at corners is much higher than that in circular geometry. And the integrated effect is that the gap of conversion is widened as well as selectivities to certain products.

4. Conclusions The geometric structure effect of cooling channels on thermal cracking of HF-1 was investigated in rectangular, square and circular channels. It demonstrated that at the same outlet temperature, HF-1 conversions in both rectangular and square channels were higher than that in the circular one. For example, the relatively gap of conversions between square and circular channels at 700 °C, 4 g/s of inlet flow rate was 17.2%. Meanwhile, higher selectivity to methane but lower selectivities to the primary pyrolysis products (such as 1-hexene and 1-heptene, etc.) were also found in quadratic channels compared with circular geometry. In addition, more carbon deposit was observed in quadratic channels especially at their corners than that in circular channels, with carbon thickness at corners of rectangular channel nearly 67% more than that at circular channel. With the help of simplified CFD simulation of fluid temperature and velocity distribution

ACS Paragon Plus Environment

Page 20 of 44

Page 21 of 44 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

Energy & Fuels

in the outlet cross-section, it is clear that gradient uniformity of fluid temperature and velocity shows difference in three channels. In rectangular and square channels, the fluid temperature in the boundary layer (mainly below 20% along the centerline) is higher and the fluid velocity is lower than that of the circular channel. This difference causes further extent of pyrolysis conversion and secondary reactions, and correspondingly the carbon deposit is formed more severely in quadratic channels. Our work indicates that structure design of cooling channel should pay attention to the shape of its circumference to make a balance between promoting thermal cracking and avoiding carbon deposition.

Supporting Information Grids of three channels (Figure S1-S3); Simulation results with chemical model (Figure S4-S5); Additional SEM images of the outlet cross-sections of three channels (Figure S6).

Acknowledgments The authors sincerely acknowledge financial support from the National Natural Science Foundation of China (21522605 and 21776210)

References (1). Gascoin, N.; Abraham, G.; Gillard, P., Synthetic and jet fuels pyrolysis for cooling and combustion applications. Journal of Analytical and Applied Pyrolysis, 2010, 89, (2), 294-306. (2). Huang, H.; Spadaccini, L. J.; Sobel, D. R., Fuel-cooled thermal management for advanced aeroengines. Journal of Engineering for Gas Turbines and Power-Transactions of

ACS Paragon Plus Environment

Energy & Fuels 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 Asme, 2004, 126, (2), 284-293. (3). Goel, P.; Boehman, A. L., Numerical Simulation of Jet Fuel Degradation in Flow Reactors. Energy & Fuels, 2000, 14, (5), 953-962. (4). Ward, T.; Ervin, J. S.; Striebich, R. C.; Zabarnick, S., Simulations of Flowing Mildly-Cracked Normal Alkanes Incorporating Proportional Product Distributions. Journal

of Propulsion and Power, 2004, 20, (3), 394-402. (5). Ward, T. A.; Ervin, J. S.; Zabarnick, S.; Shafer, L., Pressure Effects on Flowing Mildly-Cracked n-Decane. Journal of Propulsion & Power, 2005, 21, (2), 344-355. (6). Jiang, R. P.; Liu, G. Z.; Zhang, X. W., Thermal Cracking of Hydrocarbon Aviation Fuels in Regenerative Cooling Microchannels. Energy & Fuels, 2013, 27, (5), 2563-2577. (7). Wang, Z.; Guo, Y. S.; Lin, R. S., Pyrolysis of hydrocarbon fuel ZH-100 under different pressures. Journal of Analytical and Applied Pyrolysis, 2009, 85, (1-2), 534-538. (8). Jin, B.; Jing, K.; Liu, J.; Zhang, X.; Liu, G., Pyrolysis and coking of endothermic hydrocarbon fuel in regenerative cooling channel under different pressures. Journal of

Analytical and Applied Pyrolysis, 2017, 125, (Supplement C), 117-126. (9). Zhou, W.; Jia, Z.; Qin, J.; Bao, W.; Yu, B., Experimental study on effect of pressure on heat sink of n-decane. Chemical Engineering Journal, 2014, 243, 127-136. (10). Eser, S.; Altin, O.; Pradhan, B. K., Formation of carbon nanotubes from jet fuel on superalloys at moderate temperature and high pressure. Carbon, 2000, 38, (10), 1512-1515. (11). Altin, O.; Eser, S., Analysis of Carboneceous Deposits from Thermal Stressing of a JP-8 Fuel on Superalloy Foils in a Flow Reactor. Industrial & Engineering Chemistry

ACS Paragon Plus Environment

Page 22 of 44

Page 23 of 44 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

Energy & Fuels

Research, 2001, 40, (2), 589-595. (12). Altin, O.; Eser, S., Analysis of Solid Deposits from Thermal Stressing of a JP-8 Fuel on Different Tube Surfaces in a Flow Reactor. Industrial & Engineering Chemistry

Research, 2001, 40, (2), 596-603. (13). Gascoin, N.; Gillard, P.; Bernard, S.; Bouchez, M., Characterisation of coking activity during supercritical hydrocarbon pyrolysis. Fuel Processing Technology, 2008, 89, (12), 1416-1428. (14). Xie, W.; Fang, W.; Li, D.; Yan, X.; Guo, Y.; Lin, R., Coking of Model Hydrocarbon Fuels under Supercritical Condition. Energy & Fuels, 2009, 23, (6), 2997-3001. (15). DeWitt, M. J.; Edwards, T.; Shafer, L.; Brooks, D.; Striebich, R.; Bagley, S. P.; Wornat, M. J., Effect of Aviation Fuel Type on Pyrolytic Reactivity and Deposition Propensity under Supercritical Conditions. Industrial & Engineering Chemistry Research, 2011, 50, (18), 10434-10451. (16). Stiegemeier, B.; Meyer, M.; Taghavi, R., A Thermal Stability and Heat Transfer Investigation of Five Hydrocarbon Fuels. 38th AIAA/ASME/SAE/ASEE Joint Propulsion

Conference & Exhibit, American Institute of Aeronautics and Astronautics: 2002. (17). Linne, D.; Meyer, M.; Edwards, T.; Eitman, D.; Linne, D.; Meyer, M.; Edwards, T.; Eitman, D., Evaluation of heat transfer and thermal stability of supercritical JP-7 fuel.

Proceedings of the 33rd Joint Propulsion Conference and Exhibit, American Institute of Aeronautics and Astronautics: 1997. (18). Liu, G.; Wang, X.; Zhang, X., Pyrolytic depositions of hydrocarbon aviation fuels in

ACS Paragon Plus Environment

Energy & Fuels 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

regenerative cooling channels. Journal of Analytical and Applied Pyrolysis, 2013, 104, 384-395. (19). Liu, P.; Zhou, H.; Gao, X.; Zhu, Q.; Wang, J.; Li, X., An experimental and numerical investigation on thermal cracking of n-decane in the microchannel. Petroleum Science and

Technology, 2016, 34, (6), 555-561. (20). Bao, W.; Zhang, S.; Qin, J.; Zhou, W.; Xie, K., Numerical analysis of flowing cracked hydrocarbon fuel inside cooling channels in view of thermal management. Energy, 2014, 67, 149-161. (21). Zhang, S.; Qin, J.; Xie, K.; Feng, Y.; Bao, W., Thermal behavior inside scramjet cooling channels at different channel aspect ratios. Journal of Propulsion and Power, 2016,

32, (1), 57-70. (22). Zhang, S. L.; Feng, Y.; Jiang, Y. G.; Qin, J.; Bao, W.; Han, J. C.; Haidn, O. J., Thermal behavior in the cracking reaction zone of scramjet cooling channels at different channel aspect ratios. Acta Astronautica, 2016, 127, 41-56. (23). Xu, K.; Meng, H., Modeling and Simulation of Supercritical-Pressure Turbulent Heat Transfer of Aviation Kerosene with Detailed Pyrolytic Chemical Reactions. Energy & Fuels,

2015, 29, (7), 4137-4149. (24). Ely, J. F.; Hanley, H. J. M., Prediction of transport properties. 2. Thermal conductivity of pure fluids and mixtures. Industrial & Engineering Chemistry Fundamentals, 1983, 22, (1), 90-97. (25). Meng, H.; Yang, V., A unified treatment of general fluid thermodynamics and its

ACS Paragon Plus Environment

Page 24 of 44

Page 25 of 44 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

Energy & Fuels

application to a preconditioning scheme. J. Comput. Phys., 2003, 189, (1), 277-304. (26). Zhu, Y. H.; Liu, B.; Jiang, P. X., Experimental and Numerical Investigations on n-Decane Thermal Cracking at Supercritical Pressures in a Vertical Tube. Energy & Fuels,

2014, 28, (1), 466-474. (27). Kossiakoff, A.; Rice, F. O., Thermal Decomposition of Hydrocarbons, Resonance Stabilization and Isomerization of Free Radicals. Journal of the American Chemical Society,

1943, 65, (4), 590-595. (28). Savage, P. E., Mechanisms and kinetics models for hydrocarbon pyrolysis. Journal of

Analytical and Applied Pyrolysis, 2000, 54, (1-2), 109-126. (29). Kunzru, D.; Shah, Y. T.; Stuart, E. B., Thermal Cracking of n-Nonane. Industrial &

Engineering Chemistry Process Design and Development, 1972, 11, (4), 605-612. (30). Yu, J.; Eser, S., Thermal Decomposition of C10−C14 Normal Alkanes in Near-Critical and Supercritical Regions:  Product Distributions and Reaction Mechanisms. Industrial &

Engineering Chemistry Research, 1997, 36, (3), 574-584.

ACS Paragon Plus Environment

Energy & Fuels 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

Table captions Table 1. Properties and major components of HF-1. Table 2. Parameters of channel geometries.

ACS Paragon Plus Environment

Page 26 of 44

Page 27 of 44 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

Energy & Fuels

Table 1. Properties and major components of HF-1.

Parameters

Value

Mean molecular weight

149.9

Density (20 °C), g/cm3

0.7915

ASTM distillation, °C

Method

GB/T 1884 ASTM D2887

IBP, °C

190.9

10%

202.3

20%

203.8

50%

204.4

90%

205.1

FBP, °C

217.9

PONA, wt%

SH/T 0606-2005

Normal paraffin

46.5

One-ring cycloparaffin

34.5

Two-ring cycloparaffin

18.0

Three-ring cycloparaffin

0.4

Aromatics

0.6

Flash point, °C

69.0

ASTM D93

Total sulfur, ppm