EXAMINATION OF PHYSICAL PROPERTIES OF FUELS AND MIXTURES WITH ALTERNATIVE FUELS By Anne Lauren Lown A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree of Chemical Engineering – Doctor of Philosophy 2015 UNCLASSIFIED: Distribution Statement A. Approved for public release. ABSTRACT EXAMINATION OF PHYSICAL PROPERTIES OF FUELS AND MIXTURES WITH ALTERNATIVE FUELS By Anne Lauren Lown The diversity of alternative fuels is increasing due to new second generation biofuels. By modeling alternative fuels and fuel mixtures, types of fuels can be selected based on their properties, without producing and testing large batches. A number of potential alternative fuels have been tested and modeled to determine their impact when blended with traditional diesel and jet fuels. The properties evaluated include cloud point and pour point temperature, cetane number, distillation curve, and speed of sound. This work represents a novel approach to evaluating the properties of alternative fuels and their mixtures with petroleum fuels. Low temperature properties were evaluated for twelve potential biofuel compounds in mixtures with three diesel fuels and one jet fuel. Functional groups tested included diesters, esters, ketones, and ethers, and alkanes were used for comparison. Alkanes, ethers, esters, and ketones with a low melting point temperature were found to decrease the fuel cloud point temperature. Diesters added to fuels display an upper critical solution temperature, and multiple methods were used to confirm the presence of liquid-liquid immiscibility. These behaviors are independent of chain length and branching, as long as the melting point temperature of the additive is not significantly higher than the cloud point temperature of the fuel. Physical properties were estimated for several potential fuel additive molecules using group contribution methods. Quantum chemical calculations were used for ideal gas heat capacities. UNCLASSIFIED UNCLASSIFIED Fuel surrogates for three petroleum based fuels and six alternative fuels were developed. The cloud point temperature, distillation curve, cetane number, and average molecular weight for different fuel surrogates were simultaneously represented. The proposed surrogates use the experimental mass fractions of paraffins, and the experimental concentrations of mono- and diaromatics, isoparaffins, and naphthenics. The surrogates represent both low and high temperature properties better than most surrogates in the literature. Three different methods were developed to predict the cetane number of alternative fuels and their mixtures with JP-8, a military jet fuel. The same six alternative fuels were distilled, as well as blended with JP-8, and the cetane numbers measured. The Ghosh and Jaffe model represented the neat fuels with pseudocomponents to predict the cetane numbers of blends. This model worked well for the neat fuels, but the mixture behavior was predicted with incorrect curvature. The second and third methods used near infrared (NIR) and Fourier transform infrared (FTIR) spectroscopy to correlate the cetane number. The correlation provides prediction of the cetane numbers of the blends based on spectral measurements. Both the FTIR and NIR correlations are able to predict mixture cetane numbers within experimental error, but the NIR model was found to be the most reliable of all three methods. Finally, the SAFT-BACK and ESD equations of state were used to model the density and speed of sound for hydrocarbons at elevated pressures. The SAFT-BACK equation was found to be more accurate, and the model was extended to predicting the speed of sound. Mixtures of hydrocarbons were also predicted, but the SAFT-BACK is limited in capability for representing compressed alkanes heavier than octane. UNCLASSIFIED Copyright by ANNE LAUREN LOWN 2015 ACKNOWLEDGEMENTS I would first like to thank my advisor, Dr. Carl Lira, for all of his guidance during my time as his PhD student. He has been a tremendous mentor to me. I would also like to thank my committee, Dr. Dennis Miller, Dr. Andre Lee, and Dr. James Jackson for all of their support and assistance. I want to thank Dr. Lars Peereboom for all his mentoring and training. I also want to thank Ravinder Singh and Jacob Anibal for their assistance with experiments and coding. In addition, I would like to thank my collaborators: James Anderson and Sherry Mueller at Ford Motor Company; Nichole Hubble and Eric Sattler at Tank Automotive Research, Development, and Engineering Center; and Robert Morris, James Edwards, and Linda Shafer at Air Force Research Laboratory. This thesis was also supported through funding from the Defense Logistics Agency (contract numbers SP4701-09-C-0037 and SP4701-11-C-0011), and the United States Army (contract number W56HZV-13-C-0340). Finally, I would like to thank my friends and family for their extensive support and love throughout my time in graduate school. In particular, I would like to thank my husband, without whose support this never would have been possible. v TABLE OF CONTENTS LIST OF TABLES ......................................................................................................................... ix LIST OF FIGURES ..................................................................................................................... xiv CHAPTER 1 ................................................................................................................................... 1 Introduction and Background ......................................................................................................... 1 1.1 Introduction ................................................................................................................... 2 1.2 Background ................................................................................................................... 5 CHAPTER 2 ................................................................................................................................... 7 Cold Flow Properties for Blends of Biofuels with Diesel and Jet Fuels ........................................ 7 2.1 Introduction ................................................................................................................... 8 2.2 Materials and Methods .................................................................................................. 9 2.2.1 Materials ................................................................................................................. 9 2.2.2 Experimental methods .......................................................................................... 13 2.2.2.1 Gas chromatography method .......................................................................... 13 2.2.2.2 Cloud point testing (ASTM D7683) ............................................................... 13 2.2.2.3 Cloud point testing (ASTM D2500) ............................................................... 16 2.2.2.4 Cold filter plugging point (ASTM D6371) ..................................................... 17 2.3 Results and Discussion ............................................................................................... 17 2.3.1 Pure compounds .................................................................................................... 17 2.3.2 Functional group effect on cloud point temperature ............................................. 18 2.3.3 Effect of branching and chain length on cloud point temperature ........................ 22 2.3.4 Dibutyl succinate blends and liquid-liquid phase formation ................................ 24 2.3.4.1 Visual observation .......................................................................................... 24 2.3.4.2 Testing a system with a known liquid-liquid phase separation ...................... 26 2.3.4.3 Comparison of cloud point and cold filter plugging point temperatures ........ 27 2.3.5 Discussion summary ............................................................................................. 28 2.4 Conclusions ................................................................................................................. 29 CHAPTER 3 ................................................................................................................................. 30 Physical Property Prediction Methods for Fuels and Fuel Components ...................................... 30 3.1 Introduction ................................................................................................................. 31 3.2 Property development and extension .......................................................................... 32 3.2.1 Overview ............................................................................................................... 32 3.2.2 Gas phase properties ............................................................................................. 33 3.2.3 Liquid phase properties ......................................................................................... 35 3.3 Summary ..................................................................................................................... 46 CHAPTER 4 ................................................................................................................................. 47 Development of an Adaptable and Widely Applicable Fuel Surrogate ........................................ 47 4.1 Introduction ................................................................................................................. 48 vi 4.2 Materials and Methods ................................................................................................ 50 4.2.1 Materials ............................................................................................................... 50 4.2.1.1 Group A fuels .................................................................................................. 50 4.2.1.2 Group B fuels .................................................................................................. 51 4.2.2 Analytical testing .................................................................................................. 52 4.2.2.1 GC-MS for n-paraffin content ........................................................................ 52 4.2.2.2 HPLC for aromatic content ............................................................................. 54 4.2.2.3 Average molecular weight .............................................................................. 55 4.2.2.4 GCxGC for isoparaffin and naphthenes content ............................................. 58 4.2.2.5 Cloud and freezing point testing ..................................................................... 63 4.2.2.6 Distillation testing ........................................................................................... 63 4.2.3 Models................................................................................................................... 65 4.2.3.1 Cloud point model........................................................................................... 65 4.2.3.2 Cloud point modeling of literature surrogates ................................................ 67 4.2.3.3 Distillation curve model .................................................................................. 69 4.2.3.4 Cetane number model ..................................................................................... 70 4.2.4 Surrogate optimization .......................................................................................... 71 4.3 Results and Discussion ............................................................................................... 74 4.4 Conclusions ................................................................................................................. 87 CHAPTER 5 ................................................................................................................................. 88 Prediction of Cetane Number for Fuels and Fuel Mixtures .......................................................... 88 5.1 Introduction ...................................................................................................... 89 5.2 Materials and Methods ................................................................................................ 91 5.2.1 Materials ............................................................................................................... 91 5.2.1.1 Fuels .................................................................................................... 91 5.2.1.2 Pure compounds .............................................................................................. 92 5.2.2 Experimental Methods .......................................................................................... 92 5.2.2.1 Distillation....................................................................................................... 92 5.2.2.2 Derived Cetane Number ................................................................................. 93 5.2.2.3 Two Dimensional Gas Chromatography ........................................................ 94 5.2.2.4 Near-Infrared Spectroscopy (NIR) ................................................................. 94 5.2.2.5 Fourier Transform Infrared Spectroscopy (FTIR) ........................................ 100 5.2.3 Calculation Methods ........................................................................................... 106 5.2.3.1 Cetane Number by Molecular Class via the method of Ghosh and Jaffe ..... 106 5.2.3.2 Cetane Number Using Partial Least Squares Regression ............................. 110 5.3 Results and Discussion ............................................................................................. 111 5.3.1 Ghosh and Jaffe Method ..................................................................................... 111 5.3.2 NIR and FTIR Predictions .................................................................................. 118 5.3.2.1 Fuel Prediction .............................................................................................. 118 5.3.2.2 Model Validation Using Pure Components .................................................. 126 5.4 Conclusions .................................................................................................... 127 CHAPTER 6 ............................................................................................................................... 130 Isothermal Compressibility of Fuels and Fuel Mixtures ............................................................ 130 6.1 Introduction ............................................................................................................... 131 vii 6.2 Calculation Methods ................................................................................................. 134 6.2.1 ESD equation of state.......................................................................................... 134 6.2.2 SAFT-BACK equation of state ........................................................................... 135 6.2.3 Calculation method ............................................................................................. 137 6.2.4 Extension to mixtures ......................................................................................... 139 6.3 Results and Discussion ............................................................................................. 139 6.3.1 Model selection ................................................................................................... 139 6.3.2 Density prediction for pure fluids ....................................................................... 142 6.3.3 Speed of sound prediction for pure fluids ........................................................... 144 6.3.4 Speed of sound prediction for mixtures .............................................................. 153 6.4 Conclusions ............................................................................................................... 155 APPENDICES ............................................................................................................................ 157 APPENDIX A ............................................................................................................................. 158 Cloud point and pour point data for all mixtures ........................................................................ 158 APPENDIX B ............................................................................................................................. 165 Tabulated data for the prediction of physical properties ............................................................ 165 APPENDIX C ............................................................................................................................. 173 Distillation Curve Data ............................................................................................................... 173 APPENDIX D ............................................................................................................................. 181 Database of Potential Surrogate Components............................................................................. 181 APPENDIX E ............................................................................................................................. 185 Near-IR and Fourier Transform IR Spectra for all Fuel Distillations and Mixtures .................. 185 APPENDIX F.............................................................................................................................. 199 Proof of Pseudocomponent Mixtures for Cetane Number Prediction and Tabulated Cetane Numbers and Beta Values ........................................................................................................... 199 APPENDIX G ............................................................................................................................. 205 Regression Code and Coefficients .............................................................................................. 205 APPENDIX H ............................................................................................................................. 229 Derivatives for the ESD and SAFT-BACK Equations of State .................................................. 229 REFERENCES ........................................................................................................................... 240 viii LIST OF TABLES Table 2.1. Base petroleum fuel characteristics. USD – US standard #2 diesel; ESD – European standard diesel; HAD – high aromatic diesel. .................................................................. 10 Table 2.2. Evaluated compounds, their structures, and melting points. a) FAME is a mixture of multiple esters, represented here with methyl oleate, the main component in canola FAME. For FAME only, the cloud point is given, rather than the melting temperature. . 11 Table 3.1. Availability of literature data for the requested properties. Abbreviations: x – no data available in the literature. M – measured in our lab at ambient temperature and pressure. D – data available through DIPPR database [54]. P – data available from a predicted method. E – data available from experimental methods. R – data available though Reaxys [66]. #T – number of temperature points available. RT – range of temperature points available. CT – critical temperature. ................................................................................. 33 Table 4.1: Base petroleum fuel characteristics. USD – US standard #2 diesel; HAD – high aromatic diesel. ................................................................................................................. 51 Table 4.2. Table of fuel name, supplier, source, and production method for each fuel of interest. ........................................................................................................................................... 52 Table 4.3. Alternative fuels characteristics. All tests were performed at TARDEC in Warren, MI. Note a) No wax formation was found for IPK down to -78 ˚C. ....................................... 52 Table 4.4. Average molecular weights for all fuels as determined by melting point depression and GCxGC when available. GCxGC testing was not available for the USD or HAD fuels. .................................................................................................................................. 58 Table 4.5. Distillation vapor temperatures for all fuels. JP-8 through HAD have been adjusted to the normal boiling point using the Sydney Young equation. Note a) HAD temperatures from [15]; b) from D86. .................................................................................................... 65 Table 4.6. Weight fractions for the fixed components of the proposed surrogates for all fuels, split by type of testing. ...................................................................................................... 72 Table 4.7. Total weight fractions for the remaining functional classes in the surrogates. For the Group A fuels, the naphthenic and isoparaffin components were treated as one class. ... 72 Table 4.8. Summary of surrogate predictions and errors. ............................................................. 74 Table 4.9. Adjusted components for the proposed JP-8 surrogate developed using the GC-MS and HPLC data. ................................................................................................................. 82 Table 4.10. Adjusted components for the proposed HAD surrogate. ........................................... 82 ix Table 4.11. Adjusted components for the proposed USD surrogate............................................. 83 Table 4.12. Adjusted components for the proposed JP-8 surrogate developed using the GCxGC data. ................................................................................................................................... 83 Table 4.13. Adjusted components for the proposed IPK surrogate. ............................................. 84 Table 4.14. Adjusted components for the proposed HRJ surrogate. ............................................ 84 Table 4.15. Adjusted components for the proposed SPK2 surrogate. .......................................... 85 Table 4.16. Adjusted components for the proposed HRJ-8 surrogate. ......................................... 85 Table 4.17. Adjusted components for the proposed SPK surrogate. ............................................ 86 Table 4.18. Adjusted components for the proposed HRD surrogate. ........................................... 86 Table 5.1. Branched compounds used for comparison, with their cetane numbers and the method used to find the cetane number. Cetane numbers and measurement methods from the Murphy Compendium [35]. .............................................................................................. 92 Table 5.2. Derived cetane numbers for all neat fuels, distillation fractions, and fuel mixtures. .. 93 Table 5.3. Surrogate concentrations and components as determined by averaging the GCxGC results (JP-8, IPK, HRJ, and SPK2). ............................................................................... 116 Table 5.4. Surrogate concentrations and components as determined by averaging the GCxGC results (HRJ-8, SPK, and HRD) ..................................................................................... 116 Table 5.5 Average absolute error for the prediction of cetane number using NIR and FTIR data, with varying regression conditions. The testing set and training set are evaluated separately. ....................................................................................................................... 121 Table 5.6. Average absolute error for the prediction of cetane number using NIR and FTIR data, with varying regression conditions and the baseline removed. The testing set and training set are evaluated separately. ............................................................................................ 124 Table 5.7. Errors for the regression of various branched compounds using the neat fuels and distillation fractions as the training set. .......................................................................... 126 Table A.1. Cloud point temperatures for mixtures of USD and various additives (Part 1). ....... 159 Table A.2. Cloud point temperatures for mixtures of USD and various additives (Part 2). ....... 160 Table A.3. Cloud point temperatures for mixtures of HAD and various additives (Part 1). ...... 161 Table A.4. Cloud point temperatures for mixtures of HAD and various additives (Part 2). ...... 162 x Table A.5. Cloud point (CP), adjusted cloud point (CP*), and cold filter plugging point (CFPP) temperatures for mixtures of canola FAME in the HAD+ fuel. The cloud point data was collected using ASTM D2500, and CFFP data collected by ASTM D6371. The cloud point data is adjusted for comparison with ASTM D7683 using Cloud point (D2500) = Cloud point (D7683) + 1.68. ........................................................................................... 163 Table A.6. Cloud point (CP) and adjusted cloud point (CP*) temperatures for mixtures of dibutyl succinate in the HAD* fuel. The cloud point data was collected using ASTM D2500. The cloud point data is adjusted for comparison with ASTM D7683 using Cloud point (D2500) = Cloud point (D7683) + 1.68. ......................................................................... 163 Table A.7. Cloud point (CP), adjusted cloud point (CP*), and cold filter plugging point (CFPP) temperatures for mixtures of dibutyl succinate in the HAD+ fuel. The cloud point data was collected using ASTM D2500. The cloud point data is adjusted for comparison with ASTM D7683 using Cloud point (D2500) = Cloud point (D7683) + 1.68. ................... 164 Table A.8. Cold filter plugging point (CFPP) temperatures for mixtures of butyl nonanoate in the HAD fuel. ........................................................................................................................ 164 Table B.1. Constants for the regression of Cp/R for the lower temperature range. .................... 166 Table B.2. Constants for the regression of Cp/R for the high temperature range. ...................... 166 Table B.3. Constants for the regression of H/RT for the low temperature range. ...................... 166 Table B.4. Constants for the regression of H/RT for the high temperature range. ..................... 166 Table B.5. Constants for the regression of S/T for the low temperature range. ......................... 167 Table B.6. Constants for the regression of S/T for the high temperature range. ........................ 167 Table B.7. Predicted data for the viscosity of various compounds extrapolated to the critical point (Part 1). ........................................................................................................................... 168 Table B.8. Predicted data for the viscosity of various compounds extrapolated to the critical point (Part 2). ........................................................................................................................... 169 Table B.9. Predicted data for the heat capacity of various compounds extrapolated to the critical point (Part 1). .................................................................................................................. 170 Table B.10. Predicted data for the heat capacity of various compounds extrapolated to the critical point (Part 2). .................................................................................................................. 171 Table B.11. Predicted gas phase heat capacities for various compounds. Units are J/kgK. .......... 172 Table C.1. Uncorrected and corrected distillation data for JP-8. Data taken with an atmospheric pressure of 737 mmHg. ................................................................................................... 174 xi Table C.2. Uncorrected and corrected distillation data for IPK. Data taken with an atmospheric pressure of 746 mmHg. ................................................................................................... 175 Table C.3. Uncorrected and corrected distillation data for HRJ. Data taken with an atmospheric pressure of 746 mmHg. ................................................................................................... 176 Table C.4. Uncorrected and corrected distillation data for SPK2. Data taken with an atmospheric pressure of 758 mmHg. ................................................................................................... 177 Table C.5. Uncorrected and corrected distillation data for HRJ-8. Data taken with an atmospheric pressure of 753 mmHg. ................................................................................................... 178 Table C.6. Uncorrected and corrected distillation data for HRJ-8. Data taken with an atmospheric pressure of 749.6 mmHg. ................................................................................................ 179 Table C.7. Uncorrected and corrected distillation data for HRD. Data taken with an atmospheric pressure of 743.8 mmHg. ................................................................................................ 180 Table D.1. Selected properties of potential surrogate components – n-paraffins, mono-aromatics, di-aromatics, and naphthenes. ......................................................................................... 183 Table D.2. Selected properties of potential surrogate components – isoparaffins...................... 184 Table F.1. Cetane numbers and β values for n-paraffins and mono-isoparaffins used for cetane number surrogates. .......................................................................................................... 202 Table F.2. Cetane numbers and β values for multi-isoparaffins and mono-cylcoparaffins used for cetane number surrogates. ............................................................................................... 203 Table F.3. Cetane numbers and β values for mono- and di-aromatics used for cetane number surrogates. ....................................................................................................................... 204 Table G.1. Regression coefficients (RC) for the full NIR dataset (1 of 2). WL – wavelength (cm). ......................................................................................................................................... 210 Table G.2. Regression coefficients for the full NIR dataset (2 of 2). ......................................... 211 Table G.3. Regression coefficients for the background removed NIR dataset. .......................... 212 Table G.4. Regression coefficients for the full FTIR dataset (1 of 12). WN – wavenumbers (cm1 ). ..................................................................................................................................... 213 Table G.5. Regression coefficients for the full FTIR dataset (2 of 12). ..................................... 214 Table G.6. Regression coefficients for the full FTIR dataset (3 of 12). ..................................... 215 Table G.7. Regression coefficients for the full FTIR dataset (4 of 12). ..................................... 216 Table G.8. Regression coefficients for the full FTIR dataset (5 of 12). ..................................... 217 xii Table G.9. Regression coefficients for the full FTIR dataset (6 of 12). ..................................... 218 Table G.10. Regression coefficients for the full FTIR dataset (7 of 12). ................................... 219 Table G.11. Regression coefficients for the full FTIR dataset (8 of 12). ................................... 220 Table G.12. Regression coefficients for the full FTIR dataset (9 of 12). ................................... 221 Table G.13. Regression coefficients for the full FTIR dataset (10 of 12). ................................. 222 Table G.14. Regression coefficients for the full FTIR dataset (11 of 12). ................................. 223 Table G.15. Regression coefficients for the full FTIR dataset (12 of 12). ................................. 224 Table G.16. Regression coefficients for the background removed FTIR dataset (1 of 3). ......... 225 Table G.17. Regression coefficients for the background removed FTIR dataset (2 of 3). ......... 226 Table G.18. Regression coefficients for the background removed FTIR dataset (3 of 3). ......... 227 Table G.19. Regression intercepts (C) corresponding to the coefficients presented in this appendix. ......................................................................................................................... 228 xiii LIST OF FIGURES Figure 1.1. GC-MS spectra for JP-8 jet fuel. .................................................................................. 4 Figure 2.1. Examples of experimental cloud point runs. Run 1 shows interference, while run 2 show typical cloud point behavior. ................................................................................... 15 Figure 2.2. Cloud point temperature for USD in mixtures with various low molecular weight compounds. The uncertainty in these measurements is ±1.2 ˚C. ...................................... 18 Figure 2.3. Cloud point temperatures for HAD in mixtures with various low molecular weight compounds. The uncertainty in these measurements is ±1.2 ˚C. ...................................... 19 Figure 2.4. Cloud point temperatures for JP-8 in mixtures with various low molecular weight compounds. The uncertainty in these measurements is ±1.2 ˚C. ...................................... 19 Figure 2.5. Cloud point temperatures for mixtures of various diesel fuels and dibutyl succinate. Differences between HAD, HAD* and HAD+ are described in section 2.1. HAD* and HAD+ measured by ASTM D2500 and shifted as described in the text; all other cloud points were measured by ASTM D7683. The uncertainty in these measurements from ASTM D7683 is ±1.2 ˚C and for ASTM D2500 is ±2.0 ˚C. ............................................ 22 Figure 2.6. Cloud point temperatures for USD in mixtures with high molecular weight compounds. The uncertainty in these measurements is ±1.2 ˚C. ...................................... 23 Figure 2.7. Cloud point temperatures for HAD in mixtures with high molecular weight compounds. The uncertainty in these measurements is ±1.2 ˚C. ...................................... 23 Figure 2.8. Liquid-liquid phase separation in mixtures of diesel fuel and DBS, initially at -25˚C. (a) 75% DBS and 25% HAD. (b) 75% DBS and 25% HAD, 10 seconds after picture (a) was taken. (c) 75% DBS and 25% HAD. ......................................................................... 25 Figure 2.9. Comparison of cloud point experimental data and literature data [65] for dimethyl succinate in n-heptane. ...................................................................................................... 26 Figure 2.10. Cloud point (CP) and cold filter plug point (CFPP) temperatures for HAD in mixtures with: (a) fatty acid methyl esters (FAME) of Canola oil; (b) dibutyl succinate (DBS); and (c) butyl nonanoate (BN). CP, BN measured by ASTM D7683 (uncertainty = ±1.2 ˚C); CP, FAMES and CP, DBS measured by ASTM D2500 and shifted as described in the text........................................................................................................................... 28 Figure 3.1. Gas phase heat capacity for all molecules. DIPPR data is shown as points, Gaussian predictions are shown as lines........................................................................................... 34 Figure 3.2. Surface tension for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. ................................................................................................. 37 xiv Figure 3.3. Viscosity for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. ............................................................................................................ 38 Figure 3.4. Heat of vaporization for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. ................................................................................. 39 Figure 3.5. Vapor pressure for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. ................................................................................................. 40 Figure 3.6. Density for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. ............................................................................................................ 41 Figure 3.7. Liquid phase heat capacity for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. ................................................................................. 42 Figure 3.8. Thermal conductivity prediction using two methods, and comparison to the DIPPR [54] values for methyl stearate. ......................................................................................... 43 Figure 3.9. Thermal conductivity prediction using two methods, and comparison to the DIPPR [54] values for heptamethylnonane. .................................................................................. 44 Figure 3.10. Thermal conductivity prediction using two methods, and comparison to the DIPPR [54] values for TGME. ...................................................................................................... 44 Figure 3.11. Thermal conductivity for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. ................................................................................. 45 Figure 4.1. GC-MS of JP-8, with the n-paraffins labeled. ............................................................ 53 Figure 4.2. Distribution of n-paraffin content in various fuels. Coutinho sample shown for comparison [17]. Open points indicated interpolated values due to being unable to get baseline resolution on individual peaks. ........................................................................... 54 Figure 4.3. Aromatic distribution for the group A fuels. .............................................................. 55 Figure 4.4. Distribution of calculated average molecular weight of JP-8 in p-xylene. The red and green symbols represent repeat runs taken a year apart. The red and green lines represent the average of the melting point depression runs, and the blue line represents the average molecular weight found by GCxGC testing...................................................................... 57 Figure 4.5. Composition of JP-8 broken down by carbon number and functional class. ............. 60 Figure 4.6. Composition of IPK broken down by carbon number and functional class. .............. 60 Figure 4.7. Composition of HRJ broken down by carbon number and functional class. ............. 61 Figure 4.8. Composition of SPK2 broken down by carbon number and functional class. ........... 61 Figure 4.9. Composition of HRJ-8 broken down by carbon number and functional class. .......... 62 xv Figure 4.10. Composition of SPK broken down by carbon number and functional class. ........... 62 Figure 4.11. Composition of HRD broken down by carbon number and functional class. .......... 63 Figure 4.12. Deviations from experimental cloud point temperature for various jet fuel surrogates. Jet A surrogates 1 and 2 from [14]. ................................................................ 67 Figure 4.13. Deviations from experimental cloud point temperature for various high aromatic diesel (HAD) fuel surrogates. HAD surrogates 1 and 2 from [15]. .................................. 68 Figure 4.14. Deviations from experimental cloud point temperature for various US diesel (USD) fuel surrogates. Diesel surrogate 1 from [24] and surrogate 2 from [23]. ........................ 68 Figure 4.15. Experimental and predicted distillation curve for the JP-8 fuel and proposed surrogate developed using the GC-MS and HPLC data. .................................................. 77 Figure 4.16. Experimental and predicted distillation curve for the HAD fuel and proposed surrogate. ........................................................................................................................... 77 Figure 4.17. Experimental and predicted distillation curve for the USD fuel and proposed surrogate. ........................................................................................................................... 78 Figure 4.18. Experimental and predicted distillation curve for the JP-8 fuel and proposed surrogate developed using the GCxGC data. .................................................................... 78 Figure 4.19. Experimental and predicted distillation curve for the IPK fuel and proposed surrogate. ........................................................................................................................... 79 Figure 4.20. Experimental and predicted distillation curve for the HRJ fuel and proposed surrogate. ........................................................................................................................... 79 Figure 4.21. Experimental and predicted distillation curve for the SPK2 fuel and proposed surrogate. ........................................................................................................................... 80 Figure 4.22. Experimental and predicted distillation curve for the HRJ-8 fuel and proposed surrogate. ........................................................................................................................... 80 Figure 4.23. Experimental and predicted distillation curve for the SPK fuel and proposed surrogate. ........................................................................................................................... 81 Figure 4.24. Experimental and predicted distillation curve for the HRD fuel and proposed surrogate. ........................................................................................................................... 81 Figure 5.1. Cetane numbers for all alternative fuels in mixtures with JP-8.................................. 93 Figure 5.2. NIR spectra for all neat fuels. ..................................................................................... 95 Figure 5.3. NIR spectra of the neat fuels, with focus on the third overtone region. ..................... 95 xvi Figure 5.4. NIR spectra of the neat fuels, with focus on the second combination overtone region. ........................................................................................................................................... 96 Figure 5.5. NIR spectra of the neat fuels, with focus on the first combination overtone region. . 97 Figure 5.6. NIR spectra for neat IPK and its distillation fractions. Differences between spectra a difficult to see. .................................................................................................................. 98 Figure 5.7. NIR spectra for neat HRJ-8 and its distillation fractions. .......................................... 98 Figure 5.8. NIR spectra for mixtures of JP-8 and IPK.................................................................. 99 Figure 5.9. NIR spectra for mixtures of JP-8 and HRD.............................................................. 100 Figure 5.10. FTIR spectra for all neat fuels. ............................................................................... 101 Figure 5.11. FTIR spectra for all neat fuels, with a focus on the aromatic out of plane bending region. ............................................................................................................................. 102 Figure 5.12. FTIR spectra for all neat fuels, with a focus on the carbon-carbon stretch region. 103 Figure 5.13. FTIR spectra for all neat fuels, with a focus on the carbon-hydrogen stretch region. ......................................................................................................................................... 103 Figure 5.14. FTIR spectra for neat JP-8 and its distillation fractions. ........................................ 104 Figure 5.15. FTIR spectra for neat SPK2 and its distillation fractions. ...................................... 105 Figure 5.16. FTIR spectra for mixtures of JP-8 and IPK. ........................................................... 105 Figure 5.17. FTIR spectra for mixtures of JP-8 and SPK. .......................................................... 106 Figure 5.18. Plot of mono- and multi-isoparaffins, with representative lines for the values used in the cetane number prediction. ......................................................................................... 109 Figure 5.19. Predictions of the cetane number of mixtures of JP-8 and IPK compared to experimental results. Both the JP-8 and IPK fuels are varied to find the correct isoparaffin content and carbon number. ............................................................................................ 112 Figure 5.20. Predictions of the cetane number of mixtures of JP-8 and HRJ compared to experimental results. ....................................................................................................... 113 Figure 5.21. Predictions of the cetane number of mixtures of JP-8 and SPK2 compared to experimental results. ....................................................................................................... 113 Figure 5.22. Predictions of the cetane number of mixtures of JP-8 and HRJ-8 compared to experimental results. ....................................................................................................... 114 xvii Figure 5.23. Predictions of the cetane number of mixtures of JP-8 and SPK compared to experimental results. ....................................................................................................... 114 Figure 5.24. Predictions of the cetane number of mixtures of JP-8 and HRD compared to experimental results. ....................................................................................................... 115 Figure 5.25. Predicted vs experimental cetane number for all fuels, distillation fractions, and mixtures. Prediction done using the NIR data and 25 regression components. .............. 118 Figure 5.26. Predicted vs experimental cetane number for all fuels, distillation fractions, and mixtures. Prediction done using the NIR data and seven regression components.......... 119 Figure 5.27. Predicted vs experimental cetane number for all fuels, distillation fractions, and mixtures. Prediction done using the FTIR data and 26 regression components. ............ 119 Figure 5.28. Predicted vs experimental cetane number for all fuels, distillation fractions, and mixtures. Prediction done using the FTIR data and seven regression components. ....... 120 Figure 5.29. Regression coefficients plotted against the NIR spectra for the JP-8, IPK, and HRD fuels for the regression using 7 components. .................................................................. 122 Figure 5.30. Regression coefficients plotted against the FTIR spectra for the JP-8, IPK, and HRD fuels for the regression using 26 components. ................................................................ 122 Figure 5.31. Regression weights plotted against the NIR spectra for the JP-8, IPK, and HRD fuels for the regression using 7 components. .................................................................. 123 Figure 5.32. Regression weights plotted against the FTIR spectra for the JP-8, IPK, and HRD fuels for the regression using 26 components. ................................................................ 123 Figure 5.33. Regression weights plotted against the NIR spectra of JP-8, IPK, and HRD fuels with the baseline removed. ............................................................................................. 125 Figure 5.34. Regression weights plotted against the FTIR spectra for the JP-8, IPK, and HRD fuels for the regression using 7 components, with the baseline removed during regression. ......................................................................................................................................... 125 Figure 5.35. Predicted vs experimental cetane numbers of various branched compounds regressed using the neat fuels and distillation fractions as the training set. ................... 127 Figure 6.1. Prediction of the density of n-butane at various temperatures and pressures using the ESD EOS. Correlated values from the NIST Webbook [123] shown as points, ESD prediction shown as lines. ............................................................................................... 140 Figure 6.2. Prediction of the density of n-decane at various temperatures and pressures using the ESD EOS. Correlated values from the NIST Webbook [123] shown as points, ESD prediction shown as lines. ............................................................................................... 141 xviii Figure 6.3. Prediction of the density of n-butane at various temperatures and pressures using the SAFT-BACK EOS. Correlated values from the NIST Webbook [123] shown as points, SAFT-BACK prediction shown as lines. ........................................................................ 141 Figure 6.4. Prediction of the density of n-decane at various temperatures and pressures using the SAFT-BACK EOS. Correlated values from the NIST Webbook [123] shown as points, SAFT-BACK prediction shown as lines. ........................................................................ 142 Figure 6.5. Density and pressure for various small alkanes predicted using the SAFT-BACK equation at multiple temperatures. GERG correlation values shown as points [123] and predictions shown as lines. Individual molecules are (a) methane, (b) ethane, (c) propane, and (d) butane. ................................................................................................................ 143 Figure 6.6. Density and pressure for various larger alkanes predicted using the SAFT-BACK equation at multiple temperatures. GERG correlation values shown as points [123] and predictions shown as lines. Individual molecules are (a) hexane, (b) octane, (c) decane, and (d) dodecane. ............................................................................................................ 144 Figure 6.7. Speed of sound for various small alkanes predicted using the SAFT-BACK equation at multiple temperatures. GERG correlation values shown as points [116] and predictions shown as lines. Individual molecules are (a) methane, (b) ethane, (c) propane, and (d) butane. ............................................................................................................................. 146 Figure 6.8. Speed of sound for various larger alkanes predicted using the SAFT-BACK equation at multiple temperatures. GERG correlation values shown as points [116] and predictions shown as lines. Individual molecules are (a) hexane, (b) octane, (c) decane, and (d) dodecane. ........................................................................................................................ 147 Figure 6.9. Predicted heat capacities for methane at multiple temperatures. GERG correlation values from [116]. ........................................................................................................... 148 Figure 6.10. Predicted heat capacities for butane at multiple temperatures. GERG correlation values from [116]. ........................................................................................................... 149 Figure 6.11. Predicted heat capacities for octane at multiple temperatures. GERG correlation values from [116]. ........................................................................................................... 149 Figure 6.12. Predicted heat capacities for dodecane at multiple temperatures. GERG correlation values from [116]. ........................................................................................................... 150 Figure 6.13. Pressure predicted by SAFT-BACK at multiple temperatures for dodecane at various densities. The segments at high molar volume are connected through a minimum that is off scale. Red dotted line represents the close-packed density. ........................... 151 Figure 6.14. Pressure predicted by SAFT-BACK at 298.15 K for multiple n-alkanes at various densities. Curves for all molecules except propane have been truncated after the 4th root, to improve visibility. The two segments at high molar volume for propane are connected through a minimum that is off scale. .............................................................................. 152 xix Figure 6.15. Zoomed in version of Figure 6.14 to show detail in the pressure range of interest for the three largest alkanes. All roots are shown in this version of the figure, and segments are connected through a minimum which is off scale..................................................... 153 Figure 6.16. Speed of sound predicted for mixtures of ethane and butane at various temperatures, compared to GERG correlation values from [116]. (a) 200 K, (b) 300 K, and (c) 400 K. ......................................................................................................................................... 154 Figure 6.17. Speed of sound predicted for mixtures of propane and hexane at various temperatures, compared to GERG correlation values from [116]. (a) 300 K, (b) 400 K, and (c) 500 K. ................................................................................................................. 155 Figure E.1. NIR spectra for the JP-8 fuel and its distillation fractions. ...................................... 186 Figure E.2. NIR spectra for the IPK fuel and its distillation fractions........................................ 186 Figure E.3. NIR spectra for the HRJ fuel and its distillation fractions ....................................... 187 Figure E.4. NIR spectra for the SPK2 fuel and its distillation fractions ..................................... 187 Figure E.5. NIR spectra for the HRJ-8 fuel and its distillation fractions.................................... 188 Figure E.6. NIR spectra for the SPK fuel and its distillation fractions ....................................... 188 Figure E.7. NIR spectra for the HRD fuel and its distillation fractions...................................... 189 Figure E.8. NIR spectra for mixtures of JP-8 and IPK. .............................................................. 189 Figure E.9. NIR spectra for mixtures of JP-8 and HRJ. ............................................................. 190 Figure E.10. NIR spectra for mixtures of JP-8 and SPK2. ......................................................... 190 Figure E.11. NIR spectra for mixtures of JP-8 and HRJ-8. ........................................................ 191 Figure E.12. NIR spectra for mixtures of JP-8 and SPK. ........................................................... 191 Figure E.13. NIR spectra for mixtures of JP-8 and HRD. .......................................................... 192 Figure E.14. FTIR spectra for the JP-8 fuel and its distillation fractions. .................................. 192 Figure E.15. FTIR spectra for the IPK fuel and its distillation fractions. ................................... 193 Figure E.16. FTIR spectra for the HRJ fuel and its distillation fractions. .................................. 193 Figure E.17. FTIR spectra for the SPK2 fuel and its distillation fractions. ................................ 194 Figure E.18. FTIR spectra for the HRJ-8 fuel and its distillation fractions. ............................... 194 Figure E.19. FTIR spectra for the SPK fuel and its distillation fractions. .................................. 195 xx Figure E.20. FTIR spectra for the HRD fuel and its distillation fractions. ................................. 195 Figure E.21. FTIR spectra for mixtures of JP-8 and IPK. .......................................................... 196 Figure E.22. FTIR spectra for mixtures of JP-8 and HRJ........................................................... 196 Figure E.23. FTIR spectra for mixtures of JP-8 and SPK2. ....................................................... 197 Figure E.24. FTIR spectra for mixtures of JP-8 and HRJ-8. ...................................................... 197 Figure E.25. FTIR spectra for mixtures of JP-8 and SPK. ......................................................... 198 Figure E.26. FTIR spectra for mixtures of JP-8 and HRD. ........................................................ 198 xxi CHAPTER 1 Introduction and Background 1 1.1 Introduction The production of biofuels and alternative fuels is a rapidly expanding area of research. Wide varieties of molecules and production schemes are being developed for many different applications. These fuels range from traditional first generation biofuels like ethanol to second generation fuels like Fischer-Tropsch isoparaffins. Many logistical issues can arise during the consideration of fuels and fuel blends ranging from fuel standards to acceptability in distribution and combustion systems. Fuel properties which need be addressed for biofuels to be considered compatible include energy content, vapor pressure, octane/cetane number, materials compatibility, and cold flow properties. First generation biofuels are made primarily from starches, sugars, or vegetable oils. These include ethanol from corn and fatty acid methyl esters (FAMES) produced from soybean oil or animal fats. While FAMES from soybean oil are the most commonly used form of biodiesel in the United States [1], the use of this fuel is limited in cold climates due to its high cloud point temperature, even when in mixtures with petroleum diesel [1-3]. Traditional FAMES biodiesel is unacceptable as a bio-based jet fuel, due to the low temperatures encountered at high flight altitudes [4]. Without additional processing, such as distillation [5] or the use of urea to precipitate the saturated esters [6], FAME fuel is not able to achieve the required low cloud point temperatures [7-9]. Second generation biofuels are biofuels made from non-food feedstocks. This classification includes ethanol from corn stover and other lignocellulosic sources, FischerTropsch liquids, and butanol or mixed alcohols from renewable sources [10]. In addition, these fuels can be reacted together or upgraded to produce fuels with more desirable qualities [11, 12]. 2 These processes create molecules with a wide variety of functional groups and oxygen content, and can usually be tailored to produce a specific molecule type. In order to streamline the production of these fuels, and to identify which types of molecules would be desirable as fuels, evaluations of the molecule’s effects on cloud points and other fuel properties are needed. Often, these molecules are fairly new to large scale production, and as such have limited physical property data. Thermodynamic property models can be used to extrapolate property data from a limited experimental data [13]. Because of the wide variety of possible fuel components, physical testing of every possible molecule would be labor and cost intensive. Models for the physical properties of these biofuels in mixtures with traditional fuels facilitates the screening of biofuel components without the need for large production batches or extensive testing. Traditional petroleum fuels contain hundreds of different components. Figure 1.1 shows a GC-MS spectrum for JP-8, a military jet fuel. The vast majority of components are not baseline resolved and are too similar to be identified easily. In particular, isoparaffins and naphthenic species can not be resolved from one another in the analysis below, due to very similar fragmentation patterns. GC-MS is also unable to distiguish between a carbon fragment that has broken off from an aromatic ring and a fragment from another isoparaffin, thus these species cannot be resolved. GCxGC testing can determine the differences between isoparaffin and naphthenic species, but can only give a carbon number, not the specific structure of each component. 3 Figure 1.1. GC-MS spectra for JP-8 jet fuel. It is uncommon to identify and quantify all the components in a traditional fuel. Instead, surrogate mixtures are developed for petroleum fuels. A fuel surrogate uses a limited number of compounds with known physical properties that when mixed mimic the properties of a specific fuel. This surrogate can then be used to calculate properties such as cetane number, distillation curve, cloud point temperature, and compressibility. This work represents a novel approach to predicting the properties of bio- or alternative fuels based on their properties in mixtures with petroleum fuels. By modeling fuel mixtures, potential biofuels can be identified based on their properties without producing and testing large batches. The objective of this research is to evaluate different methods of developing a fuel surrogate which can predict both low and high temperature fuel properties. The fuel surrogates could then be used to predict the properties of alternative fuel mixtures. A number of potential biofuels have been tested to determine their impact on the cloud point temperature of different fuels. A model was developed to simultaneously optimize the surrogate composition to match the cloud point temperature, distillation curve, cetane number, and average molecular weight. A surrogate model is used to relate the composition of the fuel to cetane number, and those surrogates are used to predict the cetane numbers of mixtures of fuels. 4 Near-infrared and Fourier transform infrared spectroscopy is also used to predict the cetane number of fuels. Finally, the speed of sound of pure component and mixtures are predicted using an equation of state. 1.2 Background The methods for designing a surrogate can vary depending on the type of property that will be predicted and the model used to perform the predictions. Most surrogates are designed to match one or two related fuel properties, such as boiling curve [14, 15], cloud point [16-22], spray formation [23], combustion kinetics [24], or cetane number [25-32]. Cetane number does not need to meet a minimum value for jet fuels, but due to the one fuel for the battlefield program, JP-8 is being used in diesel engines as well [33-35]. When these surrogates are used to model a property for which they weren’t designed, they give very inaccurate results, with errors usually in the range of 150-250%. For example, a diesel fuel surrogate with a focus on cetane number and combustion was developed by simply adjusting the ratio of cetane to methyl naphthalene until the cetane number was correct [24]. This surrogate, however, exhibits a cloud point temperature of -38˚C, when the actual temperature is -12˚C. These deviations are likely due to the difference in approach and focus of each model. In some cases, components are picked from a database guided by a small amount of compositional information and the concentrations are adjusted until the prediction is correct [14, 15]. Other surrogates use partial compositional data for the paraffin content and only use one or two components to represent the rest of the fuel [17]. For cloud point modeling, some surrogates use a solid solution model [17, 18, 20] or a pure solid model [19], while others use a sequential precipitation technique [16, 21, 22]. For cetane number surrogates, either compositional data is used [25, 26] or components of the fuel spectrum are fit to an empirical equation [27-31]. 5 Compressibility is another fuel property which can affect the performance of biofuels in mixtures. If a fuel is too compressible, the mass of fuel injected for a given injection volume can be different based on the injection pressure. These differences can cause incomplete or insufficient combustion, which can lead to an increase in wear on the engine and higher emissions. The compressibility of a fuel is usually calculated by experimentally determining the speed of sound and density of the fuel at a wide variety of temperatures and pressures [36-38]. For some mixtures, the speed of sound can be predicted by assuming a direct relationship between the pure component speed of sound and the mole fraction [37, 39-41]. An equation of state model can also be used to predict the speed of sound of a pure fluid. This has been done for pure alkanes and binary mixtures using various modifications of the SAFT equation, as well as the Peng-Robinson equation [42-46]. To use an equation of state model, a surrogate would also need to be used. Thermodynamic models are preferable to empirical models because they represent the chemical interactions in the surrogates rather than other correlations which require a large number of tests to fit or are not accurate if the composition changes. By using thermodynamic models, all predictions are sensitive to changes in composition and the optimized surrogate is more likely to be widely applicable. Different production batches of fuels can meet the required specifications but can have different compositions [47, 48]. The surrogates implemented in this work take elements of some of the above methods and combine them. This approach allows for better representation of the fuel across both low and high temperature properties. 6 CHAPTER 2 Cold Flow Properties for Blends of Biofuels with Diesel and Jet Fuels Portions of this chapter have been previously published as Lown A.L., Peereboom L., Mueller S.A., Anderson J.E., Miller D.J., Lira C.T., Cold flow properties for blends of biofuels with diesel and jet fuels, Fuel, 117, Part A (2014) 544-51. 7 2.1 Introduction By the year 2022, it is projected that 60 billion gallons of biofuel will be required worldwide to meet government mandates [49]. There are many logistical issues for consideration of fuel blends ranging from fuel standards to acceptability in distribution and combustion systems. Various fuel properties must be evaluated for biofuels to be considered compatible, including energy content, octane/cetane number, materials compatibility, volatility and cold flow properties. In this chapter, an expanded view of possible biofuels is taken, and compounds that contain oxygen in different functional groups are considered. First-generation biofuels are made primarily from starches, sugars, or vegetable oils. This includes ethanol from corn and fatty acid methyl esters (FAME) produced from vegetable oil or animal fats. While FAME from soybean oil is the most commonly used form of biodiesel in the United States [1], the use of this fuel is limited in cold climates due to its high cloud point temperature, even when mixed with petroleum diesel [1-3]. FAME biodiesel is more problematic as a bio-jet fuel, due to the low temperatures encountered at high flight altitudes [4]. Without additional processing, such as distillation [5] or the use of urea to precipitate the saturated esters [6], FAME is not able to achieve the required low cloud point temperatures. There are many papers in the literature which describe the effects of biodiesel on the low temperature and various fuel properties of diesel fuels [7-9]. Second-generation biofuels are often defined as biofuels made from non-food feedstocks, including ethanol from corn stover and other lignocellulosic sources, Fischer-Tropsch oils from wood, and butanol or mixed alcohols from renewable sources [10]. In addition, these fuels can be reacted together or upgraded to produce fuels with more desirable qualities [11, 12, 50]. Chemical reactions create compounds with a wide variety of functional groups and oxygen 8 content, and can usually be tailored to produce a specific molecule type. Dibutyl succinate has previously been considered as a diesel-range oxygenate [51], as well as other ethers [52]. Evaluations of the functional group’s effects on cloud points and other cold flow properties are needed to streamline the development of these fuels, and to identify which types of molecules would be desirable as fuels. Twelve compounds representing different chemical classes (alkanes, ethers, esters, ketones, and diesters) were examined in fuel blends, to develop an understanding of the impact of molecular functionality on cold flow properties. ASTM D7683 and D2500 were to determine cloud point temperatures, and ASTM D6371 was used for cold filter plugging point temperatures. Multiple types of both diesel and jet fuels were used to determine the effects of the molecule types across fuel specifications. The compounds evaluated have the potential to be produced through bio-derived pathways. Most are usually not produced directly by fermentation, but are envisioned as products of the upgrading of fermentation products. 2.2 Materials and Methods 2.2.1 Materials Three petroleum diesel fuels and one jet fuel with properties shown in Table 2.1 were used, covering a range of intended application temperatures. These fuels include a representative US standard #2 diesel (USD), a representative European standard diesel (ESD), and a diesel with comparatively high aromatic content, denoted high aromatic diesel (HAD). Three samples of HAD fuel were used, obtained from three fuel batches prepared to the same specifications, but with slightly different measured properties. HAD was the initial sample, HAD+ was obtained 2 months later, and HAD* was obtained 4 months after the first sample. The same batch of HAD used in these cloud point studies was previously characterized by Windom et. al. [15]. The JP-8 9 jet fuel was donated by the US Air Force. JP-8 is the petroleum-based fuel used by the US military and it is very similar in specifications to the commercially available fuel Jet A1. However, JP-8 contains an additional additive package including a corrosion inhibitor, a lubricity enhancer, an icing inhibitor, and a static dissipater. Table 2.1. Base petroleum fuel characteristics. USD – US standard #2 diesel; ESD – European standard diesel; HAD – high aromatic diesel. Cloud point (˚C) Cold filter plugging point (˚C) Specific gravity at 15.6 ˚C Distillation T10 (˚C) Distillation T50 (˚C) Distillation T90 (˚C) Aromatics (%vol) Saturates (%vol) Olefins (%vol) Cetane Number USD ESD HAD JP-8 Test Method -15 n/a 0.847 207 256 313 -23 -33 0.834 209 263 332 -27 -36 0.859 214 253 312 -52 n/a 0.798 177 203 241 ASTM D2500 ASTM D6371 ASTM D4052 29 68 3 44 22 74 4 52 34 63 3 41 15.7 n/a 0.8 n/a ASTM D86 ASTM D1319 ASTM D5291 ASTM D613 Chemical classes for the selected biofuel compounds include diesters, esters, ketones, ethers, and alkanes. Table 2.2 summarizes all compounds tested, their structures, and their melting temperatures. For FAME, the structure of methyl oleate is shown as a representative component. All melting points are either from NIST Chemistry Webbook [53], or the DIPPR database [54] and for FAME the cloud point is tabulated. The compounds can be divided into two classes based on molecular weight. The lighter components contain 7 or 8 carbons, and the heavier components contain from 12 to 17 carbons. 10 Table 2.2. Evaluated compounds, their structures, and melting points. a) FAME is a mixture of multiple esters, represented here with methyl oleate, the main component in canola FAME. For FAME only, the cloud point is given, rather than the melting temperature. Compound Name Structure O Dihexyl ketone FAME (see note a) Melting Point (Tm) (˚C) Chemical Class 30 Ketone 3.0 Ester -29 Diester -32.5 Ketone -38 Ester -43 Ether -48 Ester -55 Diester -92 Ester -95 Alkane -97.9 Ether < -65 Ester < -65 Ester O Dibutyl succinate O O O O 4-Heptanone O Butyl nonanoate O Dihexyl ether O O Isobutyl nonanoate O O Diisobutyl succinate O O O O Butyl butyrate O Hexanes (mixture of isomers) Dibutyl ether O O Butyl ethyl hexanoate O O Ethyl hexyl nonanoate O Hexanes (99.8% purity), butyl butyrate (98% purity), dibutyl ether (99.3% purity), dihexyl ketone (97% purity), and dihexyl ether (97% purity) were obtained from Sigma Aldrich. All other chemical compounds were synthesized in the MSU laboratories and purities were 11 confirmed using gas chromatography. Water contents were analyzed using Karl-Fisher titration. All compounds had ≥97% purity and <0.1% water. Dibutyl succinate (DBS) is a possible biofuel component because it can be made from succinic acid and butanol, which are both potential products of fermentation. Dibutyl succinate and diisobutyl succinate were produced adapting the method of Kolah et al. [55], using n-butanol or isobutanol, respectively, instead of ethanol. Dibutyl ether, 4-heptanone, and butyl butyrate were selected because they can be made from butanol and butyric acid. Butyric acid is also a fermentation product [56]. Dibutyl ether can be made by acid-catalyzed dehydration of nbutanol, 4-heptanone by ketonization of butanol [57, 58], and butyl butyrate can be produced by esterification of n-butanol with butyric acid [59]. The cetane numbers (CN) of certain compounds such as DBS (CN = 15-20 [60]) or butyl butyrate (CN = 17.5 [61]) are too low to be viable as a diesel fuel in their own, but could be envisioned as blending components to provide other beneficial properties (e.g., oxygen content, reduction in exhaust particulates). The DBS CN = 15 is based on extrapolation of ASTM D613 testing DBS/FAME blends by Paragon Laboratories, Livonia, MI. The experimental cetane numbers are listed here in the format of (v/v% DBS, CN): (0%, 59.3), (20%, 46.5); (40%, 40.1), (50%, 35.6). Larger compounds (CNs measured by Paragon Laboratories, ASTM D613) such as butyl nonanoate (CN = 51), isobutyl nonanoate (CN = 45), and ethylhexyl nonanoate (CN = 59) have higher cetane numbers that are more typical of diesel fuel and thus are more desirable for diesel blending. Dihexyl ketone and dihexyl ether are included to evaluate whether the trends exhibited by esters, ethers, and ketones continue into the higher molecular weight range. Butyl ethyl hexanoate and ethyl hexyl nonanoate can be made from other biofuel intermediates and are 12 included to evaluate the effect of side chain branching for higher molecular weight compounds. 2-Ethylhexanol can be produced by the Guerbet reaction of n-butanol [62], which can then be esterified with nonanoic acid to produce ethyl hexyl nonanoate. Nonanoic acid can be produced by ozonolysis of oleic acid [63]. 2-Ethylhexanoic acid can be made by oxidizing 2-ethylhexanol, which can then be esterified with n-butanol to produce butyl ethyl hexanoate. Fatty acid methyl esters (FAME) are currently blended with diesel and are included here to compare the trends in cloud point temperatures with cold filter plugging point (CFPP) temperatures. The FAME was produced from canola oil using a standard acid-catalyzed esterification. The resulting fuel had >97% purity and was not distilled to decrease the cloud point temperature. 2.2.2 Experimental methods 2.2.2.1 Gas chromatography method A Varian 450 gas chromatograph was used with helium carrier gas at 10 mL/min. The column was a J&W Scientific DB-WAX 30 m x 0.54 mm with a film thickness of 1 μm, with on-column injection of 0.2 μL at 250 ˚C. The oven temperature program consisted of a hold at 40 ˚C for 2 min, a ramp of 10 ˚C/min to 150 ˚C, a ramp of 30 ˚C/min to 230 ˚C, and a hold for 2 min. Detection was by thermal conductivity with a 240 ˚C operating temperature. 2.2.2.2 Cloud point testing (ASTM D7683) Cloud point (CP) testing was conducted at Michigan State University (MSU), and Ford Motor Company. Data in all figures except some series in Figure 2.5 and Figure 2.9 were collected using ASTM D7683 in the MSU laboratory, as indicated in the captions. Data in Figure 2.5 and Figure 2.9 were collected at Ford Motor Company using ASTM D2500 as outlined in section 2.2.2.3. The two methods were used for convenience due to the equipment availability at 13 the physical locations of the collaborating groups. Comparison of the methods using round-robin testing shows a small systematic difference in cloud point results of 1.68 °C as reported in ASTM D7683. Because the majority of the data in this study was collected using ASTM D7683, the data from ASTM D2500 have been shifted to ASTM D7683 values using CP(D7683) = CP(D2500) - 1.68 to facilitate comparison in the figures. The original measurements are tabulated in Appendix A. Cloud point tests at MSU were performed using ASTM method D7683 with a Tanaka Mini-Pour/Cloud Point Tester (Model MPC-102A/102L). The cloud point temperature is defined as the temperature at which the first solid appears in the liquid. The commercial unit includes a sample chamber into which a cylindrical vial is inserted. The vial holds a 3.5 mL sample with 3 cm of head space above the sample. After loading the sample, the temperature decreases at a programmed rate. The unit monitors light reflected from a source in the top of the sample via a reflector on the bottom of the sample chamber. The cloud point is identified as the temperature at which the formation of a second phase in the solution scatters the light and decreases the reflected light by a predefined amount. The commercial apparatus was modified using a Varian Star 800 module interface to log data every 0.27 seconds. Data included the temperature of the bath, sample temperatures at both the lower and middle positions in the sample chamber, and the light sensor voltage response. The lower temperature sensor is used to determine the cloud point temperature, and the middle position is used to determine the pour point temperature. The absence of stirring in the apparatus causes a small temperature gradient. The solution at the bottom of the chamber becomes colder slightly faster than in the center of the chamber, thus the cloud point is visible at the bottom of 14 the chamber first. The temperature difference between the two locations was on average 4 ˚C when the cooling rate was 1.0 ˚C/minute. Figure 2.1 shows the light detection for both a successful run and a run which displays interference. The unmodified commercial apparatus self-determines and displays the cloud and pour point temperatures of the sample. However, the self-reporting was susceptible to erroneous results. Logging of the light scattering showed that several downward steps occasionally occurred in the plateau of the light scattering trace, and a step was occasionally misinterpreted by the instrument as the cloud point. The sources of the errors were suspected to be electrical interference from other instruments on the same electrical circuit, electrical noise, issues with reflectors, or physical interference, but they were not reproducible. By logging the signals and comparing results from multiple runs, it was determined that runs with this behavior did result in accurate cloud points when the trace was manually interpreted. Figure 2.1. Examples of experimental cloud point runs. Run 1 shows interference, while run 2 show typical cloud point behavior. 15 All cloud points reported were determined using the logged data from the light sensor and the lower sample temperature sensor. Data were collected at a uniform cooling rate of 1.0 ˚C/min. To decrease the run time for some samples, a variable cooling rate was used, usually 5 ˚C/min until the sample temperature was 30 ˚C warmer than the expected cloud point, and then 1 ˚C/min afterwards. A minimum of 3 runs were performed to determine the cloud point at each concentration. Before the cloud point was reached, the light detection reached a plateau. The cloud point temperature was taken to be the temperature at which the light detection was at 90% of the plateau value, and only when the decrease was a smooth curve. Figure 2.2-Figure 2.7 and Figure 2.9 display data measured using this method. When there was no evidence of supercooling, the repeatability of the measurements was ±0.22 ˚C, and the error of the thermocouple was less than ±1 ˚C. Supercooling was observed in some runs, predominately with nearly pure, or pure compounds. Supercooling is caused by a lack of nucleation sites for the formation of a solid phase, and the sample cools as a liquid below the freezing point. When nucleation occurs, the exothermic heat of fusion warms the sample back to the cloud point/freezing temperature. If supercooling was observed in a test reported here, and the solution warmed by less than 1 ˚C after nucleation occurred, the cloud point was reported as normal. When behavior other than supercooling was observed, additional runs were performed. A consistent cloud point could not be found for compositions of >85% butyl butyrate in JP-8 and therefore, these runs were excluded from the analysis. 2.2.2.3 Cloud point testing (ASTM D2500) At Ford Motor Company, cloud point was determined using an automated cloud point analyzer (Lawler model DR-14L automated CP and CFPP analyzer, Lawler Manufacturing 16 Corp., Edison, NJ) per ASTM D2500 specifications. The fuel sample was cooled at the specified rate and examined in 1 ˚C increments for the formation of insoluble materials as measured by light scattering. The cloud point was identified as the temperature at which the insoluble materials were first detected at the bottom of the sample vial. The repeatability for these measurements is ±2.0 ˚C. 2.2.2.4 Cold filter plugging point (ASTM D6371) Cold filter plugging point (CFPP) testing was performed at Ford Motor Company. The CFPP temperature was determined by an automated analyzer (Lawler model DR-14L automated CP and CFPP analyzer, Lawler Manufacturing Corp., Edison, NJ) following ASTM D6371. The CFPP temperature is defined as the temperature at which the fuel has solidified or gelled enough to plug the fuel filter. As the fuel sample was cooled at intervals of 1 ˚C, an attempt was made to draw 20 mL of the fuel sample into a pipette under a controlled vacuum through a standard wire mesh filter. If successful, the procedure was repeated at the next lower temperature. As the sample cools, insoluble materials (e.g., wax crystals) form and begin to inhibit flow through the filter. The cold filter plugging point is defined as the temperature for which the time taken to fill the pipette exceeds 60 seconds. The repeatability for these measurements is ±1.76 ˚C. 2.3 Results and Discussion 2.3.1 Pure compounds Pure components do not display a cloud point temperature. There is often significant supercooling when a pure component freezes. Supercooling is detected when the temperature increases after the first solid forms. For Figure 2.2-Figure 2.7, the cloud point temperature shown for the pure added compound is represented by that compound’s melting point temperature as listed in Table 2.2. 17 2.3.2 Functional group effect on cloud point temperature For most mixtures tested, the cloud point decreases with addition of esters, ethers, and alkanes. For diesters, a maximum is seen in the cloud point curve, with cloud points that are higher than the cloud point of the base fuel or the melting point of the pure component, which is indicative of a liquid-liquid phase separation. Figure 2.2 - Figure 2.4 compare the low molecular weight compounds with different compounds added to USD, HAD, and JP-8, respectively. When hexanes are added to USD, HAD, and JP-8, the cloud point decreases steadily with the increasing mass fraction of hexanes from the base fuel’s cloud point, though remaining considerably higher than the melting temperature for the hexanes, even at high mass fractions. The hexanes are expected to exhibit approximately ideal solution behavior in the petroleum base fuel, with the cloud point decreasing due to dilution of the high molecular weight n-paraffin content of the base fuel. Figure 2.2. Cloud point temperature for USD in mixtures with various low molecular weight compounds. The uncertainty in these measurements is ±1.2 ˚C. 18 Figure 2.3. Cloud point temperatures for HAD in mixtures with various low molecular weight compounds. The uncertainty in these measurements is ±1.2 ˚C. Figure 2.4. Cloud point temperatures for JP-8 in mixtures with various low molecular weight compounds. The uncertainty in these measurements is ±1.2 ˚C. 19 In Figure 2.2 and Figure 2.3, dibutyl ether also decreases the cloud point temperature, but to a lesser extent than hexane in blends. Comparison of the behavior indicates that the activity coefficient of the fuel paraffins in dibutyl ether is greater than in hexanes. Butyl butyrate shows a greater departure from ideal behavior than dibutyl ether, possibly due to the higher oxygen-tocarbon ratio and/or the greater polarity. For the two ketones tested (4-heptanone and dihexyl ketone), the behavior depends on the relative properties of the ketone and the base fuel. The ketones tested display a eutectic point which shifts based on the difference between the compound’s melting point temperature and the base fuel’s cloud point temperature. When the ketone’s melting point temperature was lower than or close to the cloud point temperature of the base fuel, the ketone in small additions behaves the same as the esters, as shown by 4-heptanone in Figure 2.2 - Figure 2.4. In Figure 2.2, the eutectic is expected to be very near pure 4-heptanone composition and melting temperature, though the composition is not determined explicitly by the experimental data. In Figure 2.3, the eutectic is at concentrations greater than or equal to 85 wt% 4-heptanone. For 4-heptanone in mixtures with JP-8 (Figure 2.4), the eutectic point has shifted to approximately 25 mass% 4heptanone. To the left of the eutectic point, in the absence of solid solutions, the base fuel’s nparaffins are expected to be the precipitating component, whereas to the right of the eutectic the 4-heptanone is precipitating. In general, the eutectic point of the mixture shifts towards the lower concentrations of the added compound as the cloud point of the base fuel decreases. Dibutyl succinate (DBS), a di-ester, in mixtures with the three base fuels (Figure 2.2Figure 2.4) shows a maximum cloud point that is greater than the cloud point of the base fuel, despite the fact that the melting point temperature of DBS is lower than the three base fuels. Typically, when a compound is added to a fuel, the cloud point temperature of the resulting 20 mixture is between that of the pure fuel and the melting temperature of the compound for a solid solution, or below the cloud point temperature of the two when the system displays eutectic behavior. The observed maximum in cloud point is attributed to a liquid-liquid phase separation as discussed in detail in section 2.3.4. This behavior was seen for DBS mixtures with five different base diesel fuels, including three different HAD fuel batches (HAD, HAD*, and HAD+), as shown in Figure 2.5. Fuel distillation and blending can produce products with multiple compositions which all meet specifications, but with different cloud point results [48, 64]. Each of these fuels displayed the same tendency of increased cloud point temperatures with addition of DBS, despite the low melting temperature of DBS. As shown in Figure 2.5, the differences in the cloud points of these three base fuel samples (-27 ˚C for HAD, -26 ˚C for HAD+, and -24 ˚C for HAD*) were reflected as differences in the cloud points of their mixtures with DBS. 21 Figure 2.5. Cloud point temperatures for mixtures of various diesel fuels and dibutyl succinate. Differences between HAD, HAD* and HAD+ are described in section 2.1. HAD* and HAD+ measured by ASTM D2500 and shifted as described in the text; all other cloud points were measured by ASTM D7683. The uncertainty in these measurements from ASTM D7683 is ±1.2 ˚C and for ASTM D2500 is ±2.0 ˚C. 2.3.3 Effect of branching and chain length on cloud point temperature Figure 2.6 and Figure 2.7 show the cloud points of higher molecular weight esters, ketones, and ethers in mixtures with USD and HAD respectively. With the exception of dihexyl ketone which has a considerably higher melting point temperature than the others, the behavior of the high molecular weight compounds is similar, showing a gradual reduction in cloud point with increasing concentration of the added compound. In Figure 2.6 and Figure 2.7, the eutectic point for the dihexyl ketone system (Tm = +30 ˚C) has shifted almost to the neat base fuel. Because the melting temperature of the ketone is so high, it precipitates at the cloud point instead of fuel paraffins. 22 Figure 2.6. Cloud point temperatures for USD in mixtures with high molecular weight compounds. The uncertainty in these measurements is ±1.2 ˚C. Figure 2.7. Cloud point temperatures for HAD in mixtures with high molecular weight compounds. The uncertainty in these measurements is ±1.2 ˚C. 23 Diisobutyl succinate (DIBS) exhibits behavior almost identical to DBS (Figure 2.2Figure 2.4), despite a difference of ~30 ˚C in melting temperatures between the succinates. The monoesters (butyl nonanoate, isobutyl nonanoate, butyl ethyl hexanoate, and ethyl hexyl nonanoate) all show very similar cloud point temperatures in mixtures with USD and HAD (Figure 2.6 and Figure 2.7). The cloud point temperature of butyl butyrate in mixtures with USD and HAD (Figure 2.2 and Figure 2.4) are also indistinguishable (within experimental error) from the mixtures with the higher molecular weight esters. The data suggests that the degree of branching and the length of the chain in the alkane side groups have less of an impact on the cloud point than the oxygen-containing functional group. 2.3.4 Dibutyl succinate blends and liquid-liquid phase formation When DBS or DIBS was added to the base fuel, the measured cloud points of the mixtures at intermediate compositions are higher than the melting temperature of each component and the cloud point of the base fuels. Figure 2.5 shows DBS in mixtures with a variety of base fuels with different properties and hydrocarbon compositions. This elevated cloud point for mixtures is indicative of a liquid-liquid phase separation rather than the normal solidliquid phase separation expected at a cloud point. The measurement is reported as a cloud point because the liquid-liquid phase behavior creates light scattering. The behavior is different from other cloud point data reported in this study where solids form at the cloud point temperature Multiple methods were pursued to confirm the existence of two liquid phases. 2.3.4.1 Visual observation Mixtures were placed in a test tube and allowed to equilibrate in a low temperature bath at varying temperatures to determine the phase behavior for these systems. The HAD, ESD, and USD fuels were mixed with DBS:fuel volume ratios of 50:50 and 25:75. The samples were 24 cycled from -10 ˚C to -35 ˚C in 2.5 ˚C intervals with 2 to 48 hours at each temperature, then warmed to -10 ˚C and then cooled again to -35 ˚C as before. When equilibrating for 12-24 hours at temperatures lower than the instrument-reported cloud point temperature, the mixtures separated into two liquid phases, with a cloudy phase on top and clear phase below (Figure 2.8a). The phases had similar densities and required long settling times for complete phase separation. This two-phase behavior was highly temperature-dependent; when removed from the bath, the phases homogenized within 10-15 seconds starting from the outside edges of the test tube (Figure 2.8b). As the equilibrium temperature was decreased, the upper cloudy phase became more opaque and increased in volume (Figure 2.8c) with each step until the entire test tube volume was a single cloudy phase. Once the entire solution was cloudy, the solution gelled and then solidified within an additional 5 ˚C cooling as the temperature decrease continued. Figure 2.8. Liquid-liquid phase separation in mixtures of diesel fuel and DBS, initially at -25˚C. (a) 75% DBS and 25% HAD. (b) 75% DBS and 25% HAD, 10 seconds after picture (a) was taken. (c) 75% DBS and 25% HAD. 25 2.3.4.2 Testing a system with a known liquid-liquid phase separation Dimethyl succinate in n-heptane, a mixture that exhibits a known liquid-liquid phase separation [65], was tested to experimentally demonstrate that the cloud point instrument could detect a liquid-liquid separation as a cloud point. The phase envelope [65] was determined by placing samples in a bath and chilling or heating the well-stirred sample at a rate of 0.1 ˚C/min. The experimental data shown in Figure 2.9 are in good agreement with the data of Manzini and Crescenzl [65], but exhibit some scatter, which is attributed to the faster cooling rate (1.0 ˚C/min) and lack of stirring in the cloud point instrument used in the present study. Overall, the similarity in behavior demonstrates that the cloud point instrument used in this study can identify a liquid-liquid phase separation and report it as a cloud point. Figure 2.9. Comparison of cloud point experimental data and literature data [65] for dimethyl succinate in nheptane. 26 2.3.4.3 Comparison of cloud point and cold filter plugging point temperatures When canola oil FAME was mixed with the base diesel fuel, both the cloud point and CFPP increased (Figure 2.10a), consistent with the cloud point of canola oil FAME being higher than the cloud point of the diesel fuel. When butyl nonanoate (having a lower melting point) was added, the opposite trend was observed (Figure 2.10b). In the case of canola oil FAME-diesel fuel blend, a component of the FAME precipitated at the cloud point rather than the paraffins present in the diesel fuel because the melting temperatures for saturated FAME components (in the range of 30-38 ˚C) are very high compared to the base diesel fuel. For butyl nonanoate, the cloud point and CFPP decreased because the paraffins in the diesel precipitate on cooling, and were diluted by the added compound. When DBS was added to the same fuel, however, the cloud point increased (Figure 2.10c), but the CFPP remained approximately constant over the entire concentration range. This data supports the hypothesis that the instrument-detected cloud point for these DBS-diesel mixtures was actually a liquid-liquid phase separation. 27 Figure 2.10. Cloud point (CP) and cold filter plug point (CFPP) temperatures for HAD in mixtures with: (a) fatty acid methyl esters (FAME) of Canola oil; (b) dibutyl succinate (DBS); and (c) butyl nonanoate (BN). CP, BN measured by ASTM D7683 (uncertainty = ±1.2 ˚C); CP, FAMES and CP, DBS measured by ASTM D2500 and shifted as described in the text. 2.3.5 Discussion summary Compounds with melting point temperatures lower than the petroleum fuel to which they were added generally decreased the cloud point temperature of the mixtures by dilution of the relatively high-melting temperature paraffins in the fuel. Addition of lower molecular weight compounds resulted in cloud point depressions that were dependent on functionality. Cloud point depressions at a given weight percent increased in the following order: ketone < ester < ether < small alkane. When the molecular weight of the added compound was greater, all behaved similarly. When the melting temperature of the added compound was similar to or higher than the cloud point of the base fuel, eutectic behavior sometimes interfered with the precipitation of the paraffins and/or the added compound precipitated instead. Diesters exhibited liquid-liquid 28 phase immiscibility and the cloud points of their mixtures with these diesel fuels were higher than both the cloud point of the neat diesel fuel and the melting point of the diester. 2.4 Conclusions Potential biofuel compounds have been studied in blends with petroleum diesel and jet fuels to increase understanding of blend behavior. For mixtures of diesel fuels and diesters, the diester functionality caused liquid-liquid immiscibility over the small range of carbon chain length tested here. When diesel fuels were mixed with compounds of low molecular weight, the oxygen moieties present had more effect on miscibility than the degree of branching of the carbon chains. Functional groups in order from most effective to least effective at decreasing cloud point were alkanes, ethers, esters, ketones, and diesters. Higher molecular weight compounds exhibited more consistent effects, despite up to 60 ˚C differences in melting temperatures, over the composition range until the eutectic temperature was encountered. 29 CHAPTER 3 Physical Property Prediction Methods for Fuels and Fuel Components 30 3.1 Introduction The diversity of the molecules used as second generation biofuels is continually expanding. Second generation biofuels are often defined as biofuels made from non-food feedstocks. This classification includes ethanol from corn stover and other lignocellulosic sources, Fischer-Tropsch hydrocarbons from wood, and butanol or mixed alcohols from renewable sources [10]. In addition, fuel quality can be improved by blending or reacting these fuels to produce fuels with more desirable qualities [11, 12]. These processes create molecules with a wide variety of functional groups and oxygen content, and can usually be tailored to produce a specific molecule type. In order to streamline the production of these fuels, and to identify which types of molecules would be desirable as fuels, evaluations of the molecule’s physical properties such as freezing point, vapor pressure, and density are needed. Often, these molecules are fairly new to large scale production, and as such have limited physical property data. Testing of these properties at a number of temperatures and pressures can be time and resource consuming, especially when there are multiple potential biomolecules. Thermodynamic property evaluation can be used to extrapolate property data from a limited experimental data [13]. Then, these properties can be used to predict the behavior of the new biomolecules in mixtures with traditional fuels. This chapter will apply predictive methods for fuel components. Section 3.2 discusses the method for the prediction of the properties of dibutyl succinate, as well as the extension of the properties for other fuel surrogate components to the critical point. Second 3.2.2 discusses gas phase properties, while section 3.2.3 discusses liquid phase and critical properties. 31 3.2 Property development and extension 3.2.1 Overview Collaborators at Ford Motor Company requested properties to be used in predicting combustion and spray formation of alternative fuels. Property data was requested for dibutyl succinate (DBS), methyl oleate, methyl stearate, cetane, 2,2,4,4,6,8,8-heptamethylnonane (HMN), and tripropylene glycol monomethyl ether (TGME). Heptamethylnonane is used in mixtures with cetane to create the scale used for determination of cetane numbers in diesel fuels. Ideal gas phase heat capacity, enthalpy, and entropy were requested in the temperature range of 300-5,000 K, and fitted based on a 4th degree heat capacity polynomial equation for all compounds. Liquid phase viscosity, surface tension, heat of vaporization, vapor pressure, thermal conductivity, density, and heat capacity data were requested from 0 K to the critical point for each compound. For methyl oleate, methyl stearate, cetane, HMN, and TGME, data was available for all liquid phase properties in the DIPPR database [54]. DBS has very limited literature data and is not included in the DIPPR database. Reaxys was used to find additional literature data for DBS [66]. The properties available in the literature were either experimental or predicted values, or some combination thereof. Table 3.1 details which properties are experimental and predicted, and the source for each property. For the liquid phase, the properties are given for the temperature range that is given in DIPPR. If the data did not go all the way to the critical temperature, it was because there either was no experimental data in that range, or the predictive equations used are not valid above the temperature given. In these cases, properties were extended to the critical point using the same prediction methods as those for DBS. In addition, 32 gas phase data for heat capacity for these compounds was available up to ~1200 K. Quantum calculations from Gaussian [67] extended the gas phase heat capacity to the requested 5000 K. Table 3.1. Availability of literature data for the requested properties. Abbreviations: x – no data available in the literature. M – measured in our lab at ambient temperature and pressure. D – data available through DIPPR database [54]. P – data available from a predicted method. E – data available from experimental methods. R – data available though Reaxys [66]. #T – number of temperature points available. RT – range of temperature points available. CT – critical temperature. Dibutyl Methyl Methyl 2,2,4,4,6,8,8Tripropylene glycol Cetane succinate Oleate Stearate Heptamethylnonane monomethyl ether x D-P D-P D-P D-P D-P Gas phase Cp x x x x x x Gas phase H x x x x x x Gas phase S x D-E D-E D-E D-E D-P Viscosity R - RT D - E&P D-E D-E D - E&P D - E&P Surface Tension Heat of R - 2T D-P D-P D - E&P D-P D-P Vaporization R - RT D - E&P D-E D-E D - E&P D - E&P Vapor pressure Thermal x D-P D-P D-E D-P D - E&P Conductivity M D-E D-E D - E&P D-E D-E Density Liquid heat x D-P D-E D - E&P D-P D-P capacity 3.2.2 Gas phase properties For the high temperature ideal gas heat capacity data, Gaussian was used [67]. The structure of each molecule was minimized using the BY3LP method and the 6-31G(d) basis set. Then the minimum energy structure was used to find the frequencies and thermochemical properties. The equation to calculate the statistical-mechanical heat capacity for a linear molecule [68] is 2 5 𝜃𝑣,𝐾 /𝑇 𝜃𝑣,𝐾 /𝑇 𝐶𝑉 = 𝐶𝑡 + 𝐶𝑟 + 𝐶𝑣 = 𝑅 + 𝑅 ∑ 𝑒 ( −𝜃 /𝑇 ) 2 𝑒 𝑣,𝐾 − 1 ( 3.1 ) 𝐾 where CV is the total constant volume heat capacity (J/kg), Ct is the translational contribution to the heat capacity (J/kg), Cr is the rotational contribution to the heat capacity (J/kg), Cv is the vibrational contribution to the heat capacity (J/kg), R is the ideal gas constant (J/kg K), K is the number of vibrational modes, θv,k is the characteristic vibrational temperature (K), and T is the 33 temperature of the system (K). For equation ( 3.1 ), the value of 5/2 is replaced with 3 for a nonlinear molecule. The relationship between CV (constant volume heat capacity) and CP (constant pressure heat capacity) for an ideal gas is CP = CV + R. The built-in Gaussian utility “freqchk” was used to calculate heat capacity values at multiple temperatures. The deviations in the overlapping region between the low temperature DIPPR data and the high temperature Gaussian was less than 5% for all molecules. The Gaussian data could be scaled to allow for a better alignment. Figure 3.1. Gas phase heat capacity for all molecules. DIPPR data is shown as points, Gaussian predictions are shown as lines. Entropy and enthalpy constants were found using their standard relationship to heat capacity [69]. The heat capacity values were fitted to the following equation: 𝐶𝑃 = 𝐴1 + 𝐴2 𝑇 + 𝐴3 𝑇 2 + 𝐴4 𝑇 3 + 𝐴5 𝑇 4 𝑅 where CP is the constant pressure heat capacity (J/kg K) and A1-5 are the fitted constants. 34 ( 3.2 ) Ideal gas enthalpy was calculated relative to 298.15 K. The equation for enthalpy: 𝐻 1 𝑇 𝐴2 𝐴3 𝐴4 𝐴5 𝐴6 = ∫ 𝐶𝑃 𝑑𝑇 = 𝐴1 + 𝑇 + 𝑇 2 + 𝑇 3 + 𝑇 4 + 𝑅𝑇 𝑇 𝑇𝑟𝑒𝑓 2 3 4 5 𝑇 𝐴2 2 𝐴3 3 𝐴4 4 𝐴5 5 𝐴6 = −𝐴1 𝑇𝑟𝑒𝑓 − 𝑇𝑟𝑒𝑓 − 𝑇𝑟𝑒𝑓 − 𝑇𝑟𝑒𝑓 − 𝑇𝑟𝑒𝑓 2 3 4 5 ( 3.3 ) ( 3.4 ) where H is the enthalpy (J/kg), A1-5 are the same constants as equation ( 3.2 ). Tref is the reference temperature (298.15 K), A6 is the integration constant, which has the form shown in equation ( 3.4 ). Ideal gas entropy is also calculated using a reference state of 298.15 K. The equation for entropy: 𝑇 𝑆 𝐶𝑃 𝐴3 𝐴4 𝐴5 =∫ 𝑑𝑇 = 𝐴1 ln 𝑇 + 𝐴2 𝑇 + 𝑇 2 + 𝑇 3 + 𝑇 4 + 𝐴7 𝑅 2 3 4 𝑇𝑟𝑒𝑓 𝑇 𝐴7 = −𝐴1 ln 𝑇𝑟𝑒𝑓 − 𝐴2 𝑇𝑟𝑒𝑓 − 𝐴3 2 𝐴4 3 𝐴5 4 𝑇𝑟𝑒𝑓 − 𝑇𝑟𝑒𝑓 − 𝑇𝑟𝑒𝑓 2 3 4 ( 3.5 ) ( 3.6 ) where S is the entropy (J/kgK), and A7 is the integration constant with the form shown in equation ( 3.6 ). 3.2.3 Liquid phase properties Very few property data are available for DBS in the literature. Data for a small range of temperatures was found for liquid phase heat of vaporization, surface tension, and vapor pressure. Density was also measured at one temperature, with an experimental error of less than 1%. Group contribution methods were used to find data for all liquid phase properties. When literature values were available, they were compared to the group contribution method to estimate the percent error in the estimation method. Data was also not available for all compounds for viscosity, liquid heat capacity, and thermal conductivity all the way to the critical point. The same methods that are described below are used to extend these properties to the critical point. 35 Critical properties for DBS were estimated using the method of Constantinou and Gani [13]. For all the other compounds, the values presented in the DIPPR database are used [54]. This is a group contribution method with errors generally <10%. While the acentric factor is not a critical property, it is included in this section because the same method was used. The critical temperature for DBS was found to be 715.4 K, the critical pressure 17.7 bar, the critical volume 753 cm3/mol, and the acentric factor 0.8179. The equations for each property are: ( 3.7 ) 𝑇𝐶 = 181.128 ln [∑ 𝑁𝑘 (𝑡𝑐1𝑘) + 𝑊 ∑ 𝑀𝑗 (𝑡𝑐2𝑗)] 𝑘 𝑗 where TC is the critical temperature (K), k is the type of first order component, Nk is the number of k-type first order components, tc1k is the first order constant, W is set to 0 for first order only calculations, and 1 to include second order calculations, j is the type of second order component, and tc2j is the second order constant. −2 𝑃𝐶 = [∑ 𝑁𝑘 (𝑝𝑐1𝑘) + 𝑊 ∑ 𝑀𝑗 (𝑝𝑐2𝑗) + 0.10022] 𝑘 + 1.3705 ( 3.8 ) 𝑗 where Pc is the critical pressure (bar), pc1k is the first order constant, and pc2j is the second order constant. ( 3.9 ) 𝑉𝐶 = −0.00435 + [∑ 𝑁𝑘 (𝑣𝑐1𝑘) + 𝑊 ∑ 𝑀𝑗 (𝑣𝑐2𝑗)] 𝑘 𝑗 where VC is the critical volume (cm3/mol), vc1k is the first order constant, and vc2j is the second order constant. 1⁄0.5050 𝜔 = 0.4085 {ln [∑ 𝑁𝑘 (𝑤1𝑘) + 𝑊 ∑ 𝑀𝑗 (𝑤2𝑗) + 1.1507]} 𝑘 𝑗 36 ( 3.10 ) where ω is the acentric factor, w1k is the first order constant, and w2j is the second order constant. The Brock and Bird method [13] was used to calculate surface tension. This method can give higher errors for polar molecules than for non-polar. Comparison of the predicted values to the experimental data from the literature shows an approximate error of ~3% for DBS. The trend and shape of the predicted data also matches the trends for the experimental data for the other molecules (Figure 3.2). Figure 3.2. Surface tension for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. The correlation is: ⁄ ⁄ 𝜎 = 𝑃𝐶2 3 𝑇𝐶1 3 (0.132𝛼𝑐 − 0.279)(1 − 𝑇𝑟 )11⁄9 𝑇𝑏𝑟 ln(𝑃𝐶 ⁄1.01325) 𝛼𝑐 = 0.9076 [1 + ] 1 − 𝑇𝑏𝑟 37 ( 3.11 ) ( 3.12 ) where σ is the surface tension (dyn/cm), Tr is the reduced temperature (T/Tc), and Tbr is the reduced boiling temperature (Tb/Tc). Viscosity was calculated using the Orrick and Erbar [13] method. For monoesters, the deviation between the predicted values and the literature show an approximate error of 15%. The prediction values for the highly branched compounds deviate more from experimental data than the components with linear chains (Figure 3.3). Figure 3.3. Viscosity for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. The equation for viscosity is: 𝐵 𝜂𝐿 = 𝜌𝐿 𝑀𝑒 𝐴+𝑇 ( 3.13 ) where ηL is the liquid phase viscosity (cP), ρL is the liquid phase density at 20 ˚C (g/cm3), M is the molecular weight (g/mol), and A and B are constants based off of group contribution constants. 38 The Watson relation was used to calculate the heat of vaporization, with an approximate error of 10%. The predictions for DBS again match the shape and trend of the experimental data (Figure 3.4). Figure 3.4. Heat of vaporization for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. There were two data points in the literature for DBS. The first is used as the reference point for the relation below, and the second used to evaluate the error. The relation is: 𝑛 ∆𝐻𝑣,2 1 − 𝑇𝑟,2 = ∆𝐻𝑣,1 ( ) 1 − 𝑇𝑟,1 ( 3.14 ) where ΔHv,2 is the heat of vaporization at temperature 2, ΔHv,1 is the heat of vaporization at temperature 1 (the reference temperature), Tr,2 is the reduced temperature 2, Tr,1 is the reduced temperature 1, and n is a constant (0.375). A shortcut estimation was used to find the vapor pressure [69], because no data was available for the heat of vaporization, so the Clausius-Clapeyron equation could not be used. The 39 shortcut method assumes there is a linear relationship between the natural log of the vapor pressure and the inverse of the temperature, at low temperatures. The available data has that linear trend, but the model is only a good approximation when the reduced temperature is greater than 0.5 and the vapor pressure is greater than 2 bar or 200,000 Pa. The accuracy of the correlation also depends upon prediction of the critical conditions and the acentric factor. The literature data is outside these limits, and the predictions exhibit an error of 30-300% for DBS in the low pressure range (Figure 3.5). Figure 3.5. Vapor pressure for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. The relation is: log10 𝑃𝑟𝑠𝑎𝑡 = 7 1 (𝜔 + 1) (1 − ) 3 𝑇𝑟 where Prsat is the reduced vapor pressure (Psat/PC). 40 ( 3.15 ) The modified Rackett equation was used to estimate the liquid phase density [13]. Because one density was measured, this density could be used as a reference point. The equation gives the density as the molar volume, which has units of cm3/mol. To get density in the same units as the literature data, the inverse of the molar volume is taken, and converted from moles to mass. The literature indicates the error in this method to be generally less than 1%, but no compounds from their set were as large as DBS, so the error could be larger. The shape and trend of the DBS curve matches the trend of the literature data, however (Figure 3.6). Figure 3.6. Density for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. The Rackett equation is: 𝑟𝑒𝑓 (0.29056 − 0.08775𝜔)𝜙 𝑇𝑟𝑒𝑓 2⁄7 𝑇 2⁄7 𝜙 = (1 − ) − (1 − ) 𝑇𝐶 𝑇𝐶 𝑉𝑠 = 𝑉𝑠 ( 3.16 ) ( 3.17 ) where VS is the molar volume at the specified temperature (cm3/mol), and VSref is the molar volume at Tref. 41 A modified Bondi method [13] was used for the liquid phase heat capacity. This method is good for liquids which do not associate, and is better for larger molecules. An estimated error would be <5% based on alcohols and monoesters found in the literature. However, the HMN shows errors of ~12% at higher temperatures, while the other molecules have errors less than 5%. The shapes of the curves are also consistent with the literature data (Figure 3.7). Figure 3.7. Liquid phase heat capacity for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. The Bondi method is: 𝐶𝑃 − 𝐶𝑃° 0.49 6.3(1 − 𝑇𝑟 )1⁄3 0.4355 = 1.586 + + 𝜔 [4.2775 + + ] 𝑅 1 − 𝑇𝑟 𝑇𝑟 1 − 𝑇𝑟 ( 3.18 ) where CP˚ is the ideal gas heat capacity at the calculation temperature. Thermal conductivity was estimated using both the Latini method and the Sastri method [13]. According to the literature, both methods can have errors as high as 15% for organic molecules. For methyl oleate and cetane, the Latini method and the DIPPR data are in good 42 agreement. For methyl stearate, HMN, and TGME, there is considerable uncertainty in the predictions from DIPPR, and neither estimation method provides good agreement. For methyl stearate, the Latini and Sastri methods bracket the DIPPR prediction (Figure 3.8). The Sastri method is closer to the DIPPR prediction for HMN, but is still outside DIPPR’s uncertainty (Figure 3.9). For TGME, the Latini method is within DIPPR’s stated uncertainty (Figure 3.10). Figure 3.8. Thermal conductivity prediction using two methods, and comparison to the DIPPR [54] values for methyl stearate. 43 Figure 3.9. Thermal conductivity prediction using two methods, and comparison to the DIPPR [54] values for heptamethylnonane. Figure 3.10. Thermal conductivity prediction using two methods, and comparison to the DIPPR [54] values for TGME. 44 There is significant uncertainty inherent in the prediction of the thermal conductivity, so the uncertainty could be as high as 30%. The values shown in Figure 3.11 were found using the prediction method that was closest to the experimental values. Again, the errors for the highly branched compounds were worse than the linear compounds (Figure 3.11). Figure 3.11. Thermal conductivity for all molecules. DIPPR [54] data is shown as points, and predicted values are shown as lines. The Latini method is: 𝜆𝐿 = 𝐴(1 − 𝑇𝑟 )0.38 𝐴= ⁄ 𝑇𝑟1 6 𝐴∗ 𝑇𝑏𝛼 ( 3.19 ) ( 3.20 ) 𝛾 𝑀𝛽 𝑇𝐶 where λL is the liquid phase thermal conductivity (W/m K), Tb is the normal boiling point (K), and A*, α, β, γ are all group contribution constants. The Sastri method is: 45 𝜆𝐿 = 𝜆𝑏 𝑎 𝑚 1 − 𝑇𝑟 𝑛 𝑚 =1−( ) 1 − 𝑇𝑏𝑟 where λb is the thermal conductivity at the normal boiling point (determined by group ( 3.21 ) ( 3.22 ) contribution methods), a is a structure dependent constant (0.16), and n is a structure dependent constant (0.2). 3.3 Summary Multiple methods were used to predict the physical properties of various molecules for use in fuel modeling. Most of the compounds evaluated had substantial literature data and only needed extending of the data to the critical point. DBS has very limited literature data available, and more extensive methods were required to predict its properties, including critical properties. Properties were estimated for the linear and monoester compounds with errors of less than 15%. Branched and diester compounds had larger errors, usually of less than 30%. 46 UNCLASSIFIED: Distribution Statement A. Approved for public release. CHAPTER 4 Development of an Adaptable and Widely Applicable Fuel Surrogate UNCLASSIFIED 47 UNCLASSIFIED 4.1 Introduction A wide variety of different alternative fuels are currently being produced. These can range in composition from esters from oils and fats to paraffins and isoparaffins from FischerTropsch synthesis. These processes create molecules with a wide variety of functional groups and oxygen content, and usually are tailored to produce specific molecules. Often, these molecules are fairly new to large scale production and characterization, and thus physical property data are limited. Because of the wide variety of possible fuel components, physical testing of every possible blend of new molecules would be labor and cost intensive. Creating models for the physical properties of these biofuels in mixtures with traditional fuels would allow for potential biofuel components to be identified without large production batches or extensive testing. It is uncommon to both identify and quantify all the components in a traditional fuel, making rigorous computer modeling of these fuels difficult [70]. Instead, petroleum fuels are represented with either surrogate mixtures or pseudocomponents. A fuel surrogate uses a limited number of compounds with known physical properties that when mixed approximate the properties of the target fuel. The pseudocomponent method differs from the surrogate method in that the components used to make up the surrogate are not real compounds. They are instead theoretical compounds which are given properties to meet the desired results. The surrogate method is used here because it more closely relates to the real fuel by using real components. The methods for designing a surrogate can vary depending on the type of property that will be predicted and the model used to perform the predictions [14, 15, 17-21, 23-31]. Most surrogates are designed to match one or two related fuel properties, such as boiling curve and viscosity [14, 15], cloud point [17-21], spray formation [23], combustion kinetics [24], or cetane UNCLASSIFIED 48 UNCLASSIFIED number [25-32]. Cetane number is required to be reported for jet fuels, but does not need to meet a certain value. Due to the one fuel for the battlefield program, JP-8 is being used in diesel engines as well [33-35], thus the cetane value is important in the field. When surrogates are used to model a property for which they weren’t designed, they give very inaccurate results, with errors usually in the range of 150-250%. For example, a diesel fuel surrogate with a focus on cetane number and combustion was developed by adjusting the ratio of cetane to methyl naphthalene until the cetane number was correct [24]. This surrogate, however, gives a cloud point temperature of -38 ˚C, when the actual temperature is -12 ˚C. Surrogates are developed using a variety of methods for the selection of surrogate components and compositions. In some cases, components are picked from a database guided by a small amount of compositional information until the property is matched [14, 15]. Other surrogates use partial compositional data for the paraffin content and use physical prediction methods but only use one or two components to represent the remainder of the fuel [17]. The method of property prediction can also impact the surrogate design. For cloud point modeling, some surrogates use a solid solution model [17, 18, 20] or a pure solid model [19], while others use a sequential precipitation technique [21]. For cetane number surrogates, either compositional data is used [25, 26] or components of the fuel spectrum are fit to an empirical equation [27-31]. This work combines established methods to develop a surrogate which correctly predicts multiple fuel properties simultaneously. Physical property calculations are combined with empirical adjustment of certain component concentrations to optimize a surrogate composition. Both low and high temperature properties are used for fitting. A single JP-8 fuel as well as two fuel groups were assessed at different times during research. The first group, Group A, are diesel fuels: a high aromatic diesel (HAD), and a US standard diesel (USD). The second group, Group UNCLASSIFIED 49 UNCLASSIFIED B, is made up of alternative jet fuels. Group A fuels were tested with GC-MS to determine paraffin content and HPLC for aromatic content. Group B fuel surrogate design was guided by GCxGC testing. The same JP-8 fuel has been tested with both Group A and Group B methods and JP-8 surrogates were developed based on compositions for each test. The resulting two JP-8 surrogates developed in this work are extremely similar. 4.2 Materials and Methods 4.2.1 Materials 4.2.1.1 Group A fuels The properties of two petroleum diesel fuels and one JP-8 fuel with properties shown in Table 4.1 were used for surrogate development. These fuels include a representative US standard #2 diesel (USD), and a diesel with comparatively high aromatic content, denoted high aromatic diesel (HAD). The USD and HAD fuels were both donated by Ford Motor Company, Dearborn, MI, and are the same fuels that are discussed in Chapter 2, section 2.2.1. The same batch of HAD used in these cloud point studies was previously characterized by Windom et. al. [15]. The JP-8 jet fuel was donated by the US Air Force Research Lab at Wright Patterson Air Force Base. The same batch of JP-8 (POSF 6169, also designated POSF 11732 during a repeat test) was used for all testing and characterization. JP-8 is a traditional petroleum-based fuel used by the US military and it is very similar in specifications to the commercially available fuel Jet A1. However, JP-8 contains an additional additive package including a corrosion inhibitor, a lubricity enhancer, an icing inhibitor, and a static dissipater. UNCLASSIFIED 50 UNCLASSIFIED Table 4.1: Base petroleum fuel characteristics. USD – US standard #2 diesel; HAD – high aromatic diesel. Cloud point (˚C) Distillation T10 (˚C) Distillation T50 (˚C) Distillation T90 (˚C) Aromatics (%v) Saturates (%v) Olefins (%v) Cetane Number JP-8 USD HAD Test Method -52 177 203 241 15.5 83.2 1.3 45.78 -15 207 256 313 29 68 3 44 -27 214 253 312 34 63 3 41 ASTM D2500 ASTM D86 ASTM D1319 ASTM D5291 ASTM D613 4.2.1.2 Group B fuels Six different group B fuels were evaluated. They were all given designated names which relate to their composition or manufacturing process. The IPK (POSF 7629 and POSF 11736) and HRJ (POSF 7635 and POSF 11734) fuels were donated by Wright Patterson Air Force Base, Dayton, Ohio, and the SPK2 (POSF 11737), HRJ-8 (POSF 11735), SPK (POSF 11738), and HRD (POSF 11733) fuels were provided by the US Army Tank Automotive Research, Development, and Engineering Center (TARDEC) (Table 4.2 and Table 4.3). The HRJ fuel is a mixture of multiple batches of fuel, one of which is the HRJ-8 fuel. Select fuel characteristics are listed in Table 4.3. Dynamic Fuels and Syntroleum are part of the Renewable Energy Group at the time of this writing. UNCLASSIFIED 51 UNCLASSIFIED Table 4.2. Table of fuel name, supplier, source, and production method for each fuel of interest. Fuel Isoparaffin Kerosene Fuel Designation IPK Supplier Production Method and Source Sasol Fischer-Tropsch, natural gas Hydrotreated Renewable Jet HRJ Dynamic Fuels Hydrotreating and isomerization, mixture of sources Synthetic Paraffin Kerosene SPK2 Shell Fischer-Tropsch, natural gas Hydrotreated Renewable Jet HRJ-8 UOP Hydrotreating, camelina oil Synthetic Paraffin Kerosene SPK Syntroleum Fischer-Tropsch, natural gas Hydrotreated Renewable Diesel HRD Dynamic Fuels Hydrotreating and isomerization, waste fat and oil Table 4.3. Alternative fuels characteristics. All tests were performed at TARDEC in Warren, MI. Note a) No wax formation was found for IPK down to -78 ˚C. Freeze point (˚C) Distillation T10 (˚C) Distillation T50 (˚C) Distillation T90 (˚C) Aromatics (%v) Saturates (%v) Olefins (%v) Derived Cetane Number IPK HRJ SPK2 HRJ-8 SPK HRD Test Method <-78a 164.6 177.4 201.2 0.5 99.1 0.4 -49.2 180.2 222.2 258.5 1.4 95.8 2.8 -53.5 161.6 169.0 184.5 0.6 98.6 0.8 -53.6 164.5 216.4 271.8 0.4 98.8 0.8 -58.9 169.6 205.1 256.6 0.6 98.8 0.6 -10.2 ASTM D7153 228.0 279.0 ASTM D86 297.7 5.1 ASTM D1319 93.6 ASTM D5291 1.3 31.59 57.99 58.11 56.88 60.05 64.65 ASTM D6890 4.2.2 Analytical testing 4.2.2.1 GC-MS for n-paraffin content Gas chromatography coupled with mass spectroscopy (GC-MS) was used to quantify the paraffin content. Linear paraffins create sharp, well resolved peaks on GC-MS. Standards can be used to determine response factors for various paraffins, which allows the paraffin content to be quantified, both for mass fraction and carbon number. A Varian 3800 GC was used, coupled to a Varian 2000 MS using an ion trap detector. The column used in the GC was a Supelco SLB-5ms UNCLASSIFIED 52 UNCLASSIFIED 30 m x 0.25 mm with a film thickness of 0.25 μm. The GC temperature program was: 40 ˚C for 5 min, ramp to 190 ˚C at 2 ˚C/min, ramp to 325 ˚C at 30 ˚C/min, and hold at 325 ˚C for 5.5 min. Three paraffin standards were obtained from Sigma Aldrich. One standard contained n-pentane through n-octane, one contained n-heptane through n-decane, and the final contained n-decane, n-dodecane, n-tetradecane, and n-hexadecane. These standards had a linear relationship for the response factors, and the values were extrapolated to the higher carbon numbers. The GC-MS for JP-8 is shown in Figure 4.1. The middle range of the spectrum has so many different peaks that they are unable to be baseline resolved. For the paraffins in that section, valley-to-valley Tetradecane integration is used, so potentially the paraffins are underestimated. Figure 4.1. GC-MS of JP-8, with the n-paraffins labeled. UNCLASSIFIED 53 UNCLASSIFIED Figure 4.2. Distribution of n-paraffin content in various fuels. Coutinho sample shown for comparison [17]. Open points indicated interpolated values due to being unable to get baseline resolution on individual peaks. Distilled fuels display a lognormal distribution of n-paraffins [71]. The JP-8 fuel does display the expected lognormal distribution (Figure 4.2), as does the Coutinho fuel [17]. The two diesel fuels, however, have a biased distribution towards the higher molecular weight side of the bell curve. This could be caused by the mixing of multiple distilled batches. 4.2.2.2 HPLC for aromatic content High performance liquid chromatography (HPLC) was used to quantify the aromatic content for the group A fuels using a similar method to that of Segudovic et al. [72]. This method can only distinguish between mono-, di-, and tri-aromatics. A Shimadzu LC-2010A HT pump system was used with a Shimadzu RIT-10A refractive index detector. The column was a Restek Pinnacle II Cyano 5 μm 250 x 3.2 mm. N-heptane was used for the mobile phase. Unlike the method of Segudovic et al., there was no backflow. Cyclohexane, o-xylene, 1UNCLASSIFIED 54 UNCLASSIFIED methylnaphthalene, and anthracene were used to determine the relative response for each group of aromatics, which allowed for the mass fraction of each group to be quantified. Figure 4.3. Aromatic distribution for the group A fuels. Four peaks are found through HPLC testing, though the fourth peak is very small and only present for the USD (Figure 4.3). The first peak at about 3.3 min represents all the components without aromatic rings. This peak includes paraffins as well as cyclohexanes. The second peak, at about 3.6 min, represents the mono-aromatics compounds. The second peak at 3.8 min represents the di-aromatic compounds, and the very small peak at 4.2 min is the triaromatics. 4.2.2.3 Average molecular weight The average molecular weight was found by using melting point depression in p-xylene for both group A and B fuels. Small concentrations of diesel fuels were added to p-xylene, and UNCLASSIFIED 55 UNCLASSIFIED the sample was placed in a constant temperature bath 3-5 ˚C below the expected cloud point. The temperature of the solution was recorded using a thermocouple data recorder. Before a sample starts to freeze, the temperature drops with a constant rate. Once freezing begins, the slope of the sample temperature curve changes. By plotting lines through each section and finding the intersection point, the temperature of the onset of freezing can be found. The difference between the fuel/p-xylene mixture’s freezing temperature and that of the pure p-xylene relates to the average molecular weight of the added fuel. Equation ( 4.1 ) relates the melting point calculated for the mixture to that of the pure pxylene in order to calculate the mole fraction of p-xylene in the mixture. Then, the equation for mole fraction is rearranged in equation ( 4.2 ) to use the calculated mole fraction of p-xylene to find the corresponding molecular weight of the fuel. 𝑓𝑢𝑠 ln[𝑥𝑝𝑋 ] = − ∆𝐻𝑝𝑋 𝑀𝑊𝑓𝑢𝑒𝑙 1 1 ( − ) 𝑅 𝑇𝑚,𝑚𝑖𝑥 𝑇𝑚,𝑝𝑋 𝑀𝑓𝑢𝑒𝑙 𝑀𝑊𝑝𝑋 = 1 𝑀𝑝𝑋 (𝑥 − 1) 𝑝𝑋 ( 4.1 ) ( 4.2 ) where xpX is the mole fraction of p-xylene, ΔHfuspX is the heat of fusion of p-xylene, R is the ideal gas constant, Tm,mix is the melting point of the fuel/p-xylene mixture, Tm,pX is the melting point of the pure p-xylene, MWfuel is the average molecular weight of the fuel, Mfuel is the mass of the fuel in the mixture, MWpX is the molecular weight of p-xylene, and MpX is the mass of p-xylene in the mixture. Best practices for using this method included allowing the fuel mixture to supercool by no more than 1-2 ˚C during operation. The sides of the test tube were too smooth to provide many nucleation sites, and sometimes nucleation would need to be induced by sharply tapping UNCLASSIFIED 56 UNCLASSIFIED the test tube or thermocouple. This allowed for a longer linear region during freezing, which in turn made for a more repeatable calculation of the melting temperature. Figure 4.4. Distribution of calculated average molecular weight of JP-8 in p-xylene. The red and green symbols represent repeat runs taken a year apart. The red and green lines represent the average of the melting point depression runs, and the blue line represents the average molecular weight found by GCxGC testing. Figure 4.4 shows a typical distribution of average molecular weights found by melting point depression for JP-8 for two separate runs. The bars represent the maximum and minimum calculated values at each composition, rather than the error. The symbols represent the average at each composition. The horizontal lines represent the average of values given by the symbols of corresponding color. While the distribution is large, the overall average gives a good prediction. When butyl butyrate was tested using this method, the error was 3.13%, and the absolute average error between the GCxGC results (explained in section 4.2.2.4) and the melting point depression UNCLASSIFIED 57 UNCLASSIFIED method is 1.81%. The melting point results are presented with a confidence of ±20 g/mol (Table 4.4). Table 4.4. Average molecular weights for all fuels as determined by melting point depression and GCxGC when available. GCxGC testing was not available for the USD or HAD fuels. Fuel JP-8 USD HAD IPK HRJ SPK2 HRJ-8 SPK HRD Average molecular weight (g/mol) Melting Point Depression GCxGC 162.6 159 193.3 191.0 151.0 154 170.6 176 144.7 143.4 171.9 172 163.0 166 207.8 213.3 4.2.2.4 GCxGC for isoparaffin and naphthenes content Two dimensional gas chromatography (GCxGC) was performed on the neat fuels at the Air Force Research Laboratory, Fuels and Energy Branch, located at Wright-Patterson Air Force Base in Dayton, OH. The method is discussed in detail in Striebich et al. [73, 74]. The sample is separated in the first dimension by volatility using a non-polar column. Six-second fractions of the initial separation are then directed into a polar column for separation by polarity in the second dimension. A combination of flame ionization and mass spectral detection (FID and MS) is used. This method allows for the separation and quantitation of hydrocarbon functional groups by carbon number, as well as identification of individual components that are present above a certain level. It was found that for the HRJ, HRJ-8, SPK2, and HRD fuels, there were significant amounts of isoparaffins co-eluting with the n-paraffins. For these fuels, GC-MS was employed to determine the n-paraffin content using molecular ion quantitation. The amount each n-paraffin UNCLASSIFIED 58 UNCLASSIFIED was over-predicted by GCxGC-FID was added to the iso-paraffin content of one carbon number larger. The results for the fuels are presented as stacked bar graphs in Figure 4.5-Figure 4.11. JP8 displays the expected normal distribution for a distilled fuel. The JP-8 also has a large diversity of components in it, since it is petroleum based. The IPK and SPK2 fuels (Figure 4.6 and Figure 4.8) both have a very narrow distribution. The IPK fuel is almost entirely isoparaffins, while the SPK2 fuel is split between n-paraffins and isoparaffins. The HRJ fuel has a broader distribution than either IPK or SPK2, but it does not have the same bell-curve shape (Figure 4.7). The HRJ-8 fuel has a bimodal distribution, which is very unusual for a fuel (Figure 4.9). The HRJ fuel does not have as pronounced of a bimodal distribution as the HRJ-8, likely because it is a mixture of multiple fuel batches (one of which is the HRJ-8 fuel). The SPK fuel has close to the normal distribution that the JP-8 and other distilled fuels display, but it does not have the same diversity of composition (Figure 4.10). The composition of the HRD fuel is concentrated around C16-C18, which is expected for a fuel made from a plant based oil (Figure 4.11). UNCLASSIFIED 59 UNCLASSIFIED Figure 4.5. Composition of JP-8 broken down by carbon number and functional class. Figure 4.6. Composition of IPK broken down by carbon number and functional class. UNCLASSIFIED 60 UNCLASSIFIED Figure 4.7. Composition of HRJ broken down by carbon number and functional class. Figure 4.8. Composition of SPK2 broken down by carbon number and functional class. UNCLASSIFIED 61 UNCLASSIFIED Figure 4.9. Composition of HRJ-8 broken down by carbon number and functional class. Figure 4.10. Composition of SPK broken down by carbon number and functional class. UNCLASSIFIED 62 UNCLASSIFIED Figure 4.11. Composition of HRD broken down by carbon number and functional class. 4.2.2.5 Cloud and freezing point testing Cloud point temperatures were found for the group A using the same method as outlined in Chapter 2, section 2.2.2.2. The freezing point of the group B fuels were tested at TARDEC in Warren, MI, following ASTM D7153. The difference between the two methods is that the cloud point is detected when crystals from upon cooling, and the freezing point is detected when the crystals disappear upon warming. The cloud point method can suffer from supercooling effects, while the freezing point method would not. 4.2.2.6 Distillation testing Distillation curves were collected for the group B fuels and JP-8. Two liter distillations were performed with a method scaled to mimic that of ASTM D86. Two liters of each fuel was placed into a 3 L three neck round bottom flask. Nitrogen was fed through one of the necks of the flask to purge the headspace of oxygen prior to heating. The flow rate of nitrogen was 40 UNCLASSIFIED 63 UNCLASSIFIED mL/min, and it was kept at this flow rate throughout distillation to ensure an oxygen-free environment inside the vessel. Thermocouples were placed in the pot liquid and at the intake to the condenser. The flask and distillation head were insulated using fiberglass insulation up to the condenser intake. Both a Liebig and Friedrichs condenser were used in series due to the high vapor flow rate. Heating was adjusted until the distillation rate was approximately 100 mL per 4-5 min, which corresponds to approximately 4 L/min of vapor flow. This distillation rate corresponds volumetrically to a scaled ASTM D86 distillation rate. Both the vapor and pot temperatures were recorded at each 100 mL distilled. The pot was usually 9-10 ˚C warmer than the vapor temperature. The vapor temperature corresponds more closely to what the equilibrium boiling temperature of the fuel would be, since the pot is usually superheated, and for this reason the vapor temperatures were used for all distillation curves. The distillation was ended after 95% of the total volume was distilled. The sample that was collected between first boiling and the first 100 mL (5% of the total sample volume) was designated as the 5% sample, the next section from 5%-10% was designated as the 10%, etc. The atmospheric pressure was recorded prior to all runs. This allowed the vapor temperatures to be adjusted using the Sydney Young equation to the correct boiling temperature at standard atmospheric pressure [75, 76]. The uncorrected data for the group B fuels can be found in Appendix C. Two of the fuels were not distilled in this work. The values used for the HAD fuel distillation curve are from Windom et. al [15], and the USD values are from the D86 fuel specification sheet from Ford Motor Company. UNCLASSIFIED 64 UNCLASSIFIED Table 4.5. Distillation vapor temperatures for all fuels. JP-8 through HAD have been adjusted to the normal boiling point using the Sydney Young equation. Note a) HAD temperatures from [15]; b) from D86. Volume % Distilled 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 Boiling Range JP-8 167.7 183.8 187.2 189.8 192.5 196.1 197.8 200.4 202.3 206.3 209.3 212.3 215.3 219.4 223.4 228.4 233.4 239.4 245.4 256.5 88.7 IPK 162.5 170.3 171.5 172.6 173.5 175.0 176.0 177.4 178.8 180.2 181.8 183.6 185.5 187.6 190.0 192.7 195.9 199.7 205.8 214.8 52.3 HRJ 160.9 184.4 192.2 197.8 201.8 207.8 211.8 215.8 220.8 224.8 229.8 232.8 237.9 241.9 245.9 250.9 255.9 259.9 264.9 270.9 110.0 Vapor Temperature (˚C) SPK2 HRJ-8 SPK 159.8 147.8 146.7 165.2 173.8 176.8 165.9 177.6 180.5 166.2 181.4 183.9 166.9 185.3 187.1 167.5 189.8 191.0 168.2 195.6 194.7 168.9 200.4 199.0 169.6 208.4 203.6 170.5 215.4 206.6 171.6 222.4 211.6 172.5 230.4 216.6 173.8 240.4 221.6 174.9 249.4 226.6 176.8 256.4 231.6 178.2 263.5 237.6 180.8 270.5 243.6 182.6 275.5 251.7 186.2 279.5 260.7 191.5 283.5 270.7 31.7 135.7 124.0 HRD HADa 164.8 209.1 231.0 220.3 243.0 225.6 255.0 229.3 263.0 233.4 269.1 237.6 274.1 242.1 279.1 246.8 283.1 251.8 286.1 256.8 289.1 262.1 291.1 267.5 294.1 274.3 296.1 280.8 298.1 287.0 300.1 293.6 302.1 301.3 305.1 309.6 308.1 320.1 314.1 342.0 149.3 132.9 USDb 177.8 197.8 208.9 224.4 236.7 247.2 256.7 266.7 277.2 290.6 311.1 332.2 349.4 171.7 4.2.3 Models 4.2.3.1 Cloud point model A variety of methods have been described in the literature for the determination of cloud point temperature or wax appearance temperature for petroleum based fuels. Two different models were evaluated in this work for the prediction of cloud point temperature. These predictions require the use of surrogates to represent the composition of the fuel. The first assumes that only the component with the warmest melting point will solidify [77, 78]. The second allows a solid solution to form [18, 22, 79]. UNCLASSIFIED 65 UNCLASSIFIED If only one surrogate component is allowed to solidify, the equation takes the form: 𝑓𝑢𝑠 ln[𝑥2 𝛾2 ] = − ∆𝐻2 𝑅 1 1 ( − ) 𝑇 𝑇𝑚,2 ( 4.3 ) where x2 is the mole fraction of component 2 in liquid phase, γ2 is the activity coefficient of component 2 in liquid phase, H2fus (J/mol) is the heat of fusion of component 2, Tm,2 (K) is the melting temperature of component 2, and T is the cloud point temperature (K). An activity coefficient model can be used, or the liquid phase can be assumed to be ideal. A more accurate representation of the solid phase is to allow the solid phase to be a mixture of components: 𝑓𝑢𝑠 𝑥𝑖𝐿 𝛾𝑖 ∆𝐻 𝐾𝑖 = 𝑆 = 𝑒𝑥𝑝 [ 𝑖 𝑅 𝑥𝑖 1 1 ( − )] 𝑇 𝑇𝑚,𝑖 ( 4.4 ) where Ki is the equilibrium ratio of component i between the liquid and solid phases, xiS is the mole fraction of i in the solid phase, and xiL is the mole fraction of i in the liquid phase. The cloud point temperature is found by iteration: T is changed until the sum of xiS is equal to 1. This also allows us to manipulate which classes of components are allowed to precipitate by forcing the mole fraction of the other components to be zero. In order to evaluate which model and conditions to use for surrogate modeling, five sets of assumptions were evaluated using 6 different fuel surrogates. The first two calculations were done using equation ( 4.3 ) and assuming either ideal behavior (designated “Pure 1 Component Solid (Ideal)”), or using the UNIFAC model for the activity coefficient (designated “Pure 1 Component Solid (UNIFAC)”). The remaining three conditions used equation ( 4.4 ) with the UNIFAC model, and allow either all components (designated “Solid Solution”), only the alkane components (designated “Solid Solution of Alkanes Only”), or only the aromatic components to UNCLASSIFIED 66 UNCLASSIFIED precipitate (designated “Solid Solution of Aromatics Only”). Each group A fuel is represented by two surrogates from the literature. 4.2.3.2 Cloud point modeling of literature surrogates For Figure 4.12-Figure 4.14, all predicted cloud point temperatures are represented using their deviation from the experimental cloud point for the fuel. The experimental cloud point is shown as the horizontal line in all figures. The predicted cloud point is shown above or below the bar, and the bar is used to show the deviation relative to the experimental cloud point. Figure 4.12. Deviations from experimental cloud point temperature for various jet fuel surrogates. Jet A surrogates 1 and 2 from [14]. UNCLASSIFIED 67 UNCLASSIFIED Figure 4.13. Deviations from experimental cloud point temperature for various high aromatic diesel (HAD) fuel surrogates. HAD surrogates 1 and 2 from [15]. Figure 4.14. Deviations from experimental cloud point temperature for various US diesel (USD) fuel surrogates. Diesel surrogate 1 from [24] and surrogate 2 from [23]. UNCLASSIFIED 68 UNCLASSIFIED The calculation can be constrained so that only certain surrogate components are modeled in the precipitate. Since the paraffins have been found to be the main component in the precipitate [17], a solid solution of only n-paraffins was used for the optimization predictions. The study was configured to allow all components to be in the solid phase, or limit it to only certain hydrocarbon classes, such as paraffins or aromatics. Coutinho [17] and Dirand [71] state that the components which come out of solution during an experimental cloud point of a fuel are the paraffins, and that the distribution of paraffins needs to be accurately represented for an accurate cloud point temperature. To simplify optimization, only the paraffins are allowed to precipitate during iterations and a solid solution of alkanes only is assumed. The solid phase is also assumed to be an ideal solution. 4.2.3.3 Distillation curve model A single-stage batch distillation model was used to model the distillation curve. The bubble point temperature of the mixture is calculated using Raoult’s law [69]. A fixed, small volume of liquid is evaporated using the equilibrium vapor phase concentrations at the bubble point temperature (equation ( 4.5 )). To simplify calculations, the K-ratio is used (equation ( 4.6 )), which is a rearrangement of Raoult’s law [69]. The new mole fractions in the liquid phase are calculated by subtracting the change in the liquid phase mole fraction (dxi in equation ( 4.5 )), and are used to calculate the new bubble point temperature. The iterations continue until a component mole fraction in the liquid phase reaches 0.9999, meaning that there is then only one component remaining in the liquid phase. 𝑑𝑥𝑖 = (𝐾𝑖 − 1)𝑥𝑖 𝑑(𝑙𝑛[𝑛𝐿 ]) 𝑦𝑖 𝑃𝑖𝑠𝑎𝑡 𝐾𝑖 = = 𝑥𝑖 𝑃 UNCLASSIFIED 69 ( 4.5 ) ( 4.6 ) UNCLASSIFIED here dxi is the change in liquid phase mole fraction of component i, Ki is the ratio of gas to liquid mole fraction of component i, xi is the equilibrium liquid phase mole fraction of component i., d(ln[nL]) is the amount of liquid allowed to evaporate during each iteration step, yi is the equilibrium gas phase mole fraction of component i, Pisat is the saturation pressure of component i at the bubble point temperature, and P is the system pressure, which is set to be atmospheric. ASTM D86 is the accepted standard experimental method for determining the distillation curve of a fuel. However, this method can be inconsistent and the products of distillation are not necessarily at equilibrium during operation [70, 80]. The initial boiling temperature can vary by up to 30 ˚C depending on how well the heating rate is controlled, and the temperatures are almost always higher than the equilibrium boiling temperature. Because of these inconsistencies, it is difficult to model a D86 distillation using an equilibrium model. By using a larger distillation volume, a stirred distillation flask, and insulating the still head up to the thermocouple, the distillation method used for the group B fuels will give results that are closer to equilibrium temperatures. 4.2.3.4 Cetane number model The cetane number of each surrogate mixture was calculated using the method of Ghosh and Jaffe [26]. The calculation is done using the cetane number, volume fraction, and β value of each component (equation ( 4.7 )). The β values are fitted parameters which represent the impact that a specific class of fuel component has on the cetane number. For example, n-paraffins have a low β value and isoparaffins have a high β value. 𝐶𝑁 = ∑ 𝑖 𝑣𝑖 𝛽𝑖 𝐶𝑁𝑖 𝑣𝑖 𝛽𝑖 UNCLASSIFIED 70 ( 4.7 ) UNCLASSIFIED where CN is the cetane number of the mixture, vi is the volume fraction of component i, βi is the blending coefficient of component i, and CNi is the cetane number of component i. This method requires a cetane number for each potential surrogate component. Ghosh and Jaffe [26] give cetane numbers for components with different carbon numbers in each functional class. The beta values and cetane numbers for each potential surrogate component are taken directly from this paper, unless cetane values for the specific isomer were specified in the Murphy Compendium [61]. Each surrogate is made up of n-paraffins (β = 0.5212), isoparaffins (β = 7.3717), naphthenes (β = 0.0727), mono-aromatics (β = 3.1967), and di-aromatics (β = 3.1967). The cetane method will be discussed in more detail in Chapter 5. 4.2.4 Surrogate optimization To develop the surrogate for each fuel, the compositional data was used as a basis. The distribution of the paraffins was fixed based on the amount found in the GC-MS for the group A fuels, and the GCxGC for the group B fuels and JP-8. Methyl naphthalene was used for the diaromatic component for all fuels. The total mass fractions of mono-, di-, and tri-aromatics were fixed to match the compositional analysis, using the HPLC data for the group A fuels and the GCxGC data for the group B fuels. For the group A fuels, the difference between the naphthenes and isoparaffins could not be identified using GC-MS, so the total amount for both groups was found by subtracting the composition of the other classes from 1. The total mass fraction of naphthenes was fixed using the GCxGC data for the group B fuels, and the isoparaffin fraction was set as the balance of the other components. The exact paraffin distributions and di-aromatic mass fractions are presented in Table 4.6, and the remaining mass fractions are presented in Table 4.7 for all fuel surrogates. Then, components within each class were selected, and the properties of the surrogate were calculated as outlined in section 4.2.3. The boiling curve, cloud UNCLASSIFIED 71 UNCLASSIFIED point temperature, average molecular weight, and cetane number of the neat fuel were used to optimize the composition of each component within the overall constraint on each class. Table 4.6. Weight fractions for the fixed components of the proposed surrogates for all fuels, split by type of testing. GC-MS / HPLC GCxGC n-Paraffin Carbon JP-8 HAD USD JP-8 IPK HRJ SPK2 HRJ-8 SPK HRD Number C6 0.0008 C7 0.0014 0.0004 0.0017 0.0000 0.0040 0.0003 0.0004 0.0005 0.0010 C8 0.0058 0.0004 0.0023 0.0056 0.0001 0.0064 0.0151 0.0118 0.0162 0.0020 C9 0.0172 0.0012 0.0038 0.0188 0.0006 0.0095 0.1875 0.0171 0.0266 0.0031 C10 0.0361 0.0031 0.0055 0.0426 0.0006 0.0137 0.2171 0.0110 0.0313 0.0065 C11 0.0392 0.0038 0.0072 0.0463 0.0006 0.0138 0.0375 0.0079 0.0284 0.0074 C12 0.0321 0.0044 0.0089 0.0407 0.0005 0.0124 0.0031 0.0060 0.0214 0.0065 C13 0.0250 0.0051 0.0106 0.0312 0.0005 0.0092 0.0001 0.0054 0.0141 0.0063 C14 0.0130 0.0057 0.0123 0.0212 0.0003 0.0140 0.0026 0.0082 0.0083 C15 0.0057 0.0064 0.0110 0.0090 0.0002 0.0038 0.0038 0.0036 0.0095 C16 0.0020 0.0066 0.0090 0.0031 0.0001 0.0070 0.0022 0.0014 0.0278 C17 0.0052 0.0060 0.0011 0.0003 0.0016 0.0005 0.0171 C18 0.0028 0.0043 0.0003 0.0002 0.0001 0.0490 C19 0.0017 0.0032 0.0034 C20 0.0003 0.0015 0.0024 C21 0.0002 0.0005 0.0014 C22 0.0009 C23 0.0005 C24 0.0003 C25 0.0002 C26 0.0001 Subtotal weight 0.1784 0.0866 0.0469 0.2218 0.0035 0.0942 0.4609 0.0700 0.1524 0.1538 fraction paraffins Weight fraction di0.0120 0.0590 0.0350 0.0158 0.0002 0.0002 0.0001 0.0001 0.0001 0.0107 aromatics Table 4.7. Total weight fractions for the remaining functional classes in the surrogates. For the Group A fuels, the naphthenic and isoparaffin components were treated as one class. GC-MS / HPLC GCxGC JP-8 HAD USD JP-8 IPK HRJ SPK2 HRJ-8 SPK HRD Mono0.148 0.223 0.237 0.1229 0.0044 0.0004 0.0006 0.0002 0.0002 0.0225 Aromatic Naphthenic 0.2152 0.0013 0.0171 0.0064 0.0137 0.0207 0.0439 0.6616 0.6712 0.6434 Isoparaffin 0.4243 0.9904 0.8881 0.5320 0.9160 0.8265 0.7691 A database of potential surrogate components was developed. For each component, the tabulated properties were: melting point temperature, freezing temperature, liquid density at 25 UNCLASSIFIED 72 UNCLASSIFIED ˚C, heat of fusion, UNIFAC groups, vapor pressure equation coefficients, cetane numbers, and β values. The vapor pressure equation is the same one used in the DIPPR database [54]. For all potential components, the DIPPR database [54] and the NIST Webbook [53] were used to find the required properties, with the exception of cetane number and β values, which were found as discussed in section 4.2.3.4. A table of all the potential surrogate components and select properties can be found in Appendix D. Surrogates were developed using the Matlab function “fgoalattain” to adjust composition while constraining the weight fractions in each functional class to minimize a composite objective function including the cetane number, distillation curve, cloud point, and average molecular weight. Each optimization used an initial composition based on component freezing temperatures and boiling temperatures. Multiple initial guesses were given during the course of optimization to rule out local minima, as well as to help select the optimum components. Each property that was being used as a metric for the optimization was also given a weight. The cloud point temperature and distillation curve objective functions were weighted with a value of 100, while the cetane number and average molecular weight were assigned a weight of 10. The errors for each property were calculated using the square of the difference between the calculated property and the desired property, with the exception of the distillation curve. For the distillation curve, an error was calculated at each experimental point using a point in the calculated curve with the same volume distilled. The square of the differences at each point in the experimental curve was calculated and summed, then divided by the total number of points. These errors were set as the values of a vector and the minimization of the values was the objective function for fgoalattain. UNCLASSIFIED 73 UNCLASSIFIED 4.3 Results and Discussion The optimized surrogates had differing degrees of ability to correctly predict both low and high temperature properties. Two JP-8 surrogates have been proposed. The first uses the GCMS and HPLC data to define the paraffin and aromatic content, and allowed the program adjust the naphthenic and isoparaffinic components directly. This was the method used for the group A fuels as well. The second JP-8 surrogate uses the GCxGC information to fix the total content of all groups, like the group B fuels. A summary of the surrogates, their predictions, and percent errors are shown in Table 4.8. The errors in cloud point prediction temperature are represented as error in the absolute temperature measurement (K), while the values are shown in degrees Celsius. Table 4.8. Summary of surrogate predictions and errors. Fuel JP-8 GC-MS / HPLC HAD USD JP-8 GCxGC IPK HRJ SPK2 HRJ-8 SPK HRD Number of Adjusted Components Predicted Cloud Point (˚C) % Error Predicted Cetane Number % Error Surrogate Ave. Mol. Weight % Error 4 -35.76 7.42% 44.40 -3.01% 156.92 -3.47% 7 6 -23.90 -17.30 1.32% -0.83% 44.16 44.16 0.37% 7.71% 187.84 190.14 -1.65% -1.63% 6 -30.14 9.96% 47.73 4.25% 161.47 -0.67% 6 9 5 9 9 6 -46.48 -32.67 -47.95 -40.23 -38.15 -1.30 16.24% 7.45% 2.59% 6.16% 9.76% 3.44% 37.73 50.45 48.07 50.01 50.93 57.37 19.44% -13.00% -17.28% -12.08% -15.19% -11.26% 152.40 170.11 140.45 168.58 170.67 220.12 0.93% -0.29% -2.94% -1.96% 4.71% 5.94% The predicted cetane numbers have a low percent error for the petroleum based group A fuels, with errors less than 7%, while the alternative fuels have errors of up to 19.5%. The larger errors with the alternative fuels are attributed to the high isoparaffin content. The scatter of UNCLASSIFIED 74 UNCLASSIFIED isoparaffin cetane numbers and the shortcomings of the Ghosh and Jaffe method create large uncertainties, and the issues will be addressed further in Chapter 5. All surrogates do a good job of matching the average molecular weight, with errors of less than 5%. For a typical fuel, an error of about 8% represents a change in composition by one carbon number. The exception is the HRD surrogate. The HRD surrogate has larger errors due to the DIPPR database not containing branched isoparaffins in the boiling range required to fit the distillation curve. The cloud point temperatures are predicted within 8% error for the group A fuels and 16.5% error for the group B fuels. The two JP-8 surrogates developed here differ in the distribution of the paraffin content (Table 4.6) and the limits set by the greater resolution obtained by the GCxGC testing. The cloud point temperature prediction is significantly dependent on the paraffin content and distribution. Because the GCxGC testing detects almost 5 wt% more paraffins than the GC-MS, the cloud point predicted for the surrogates differs by 5.5 ˚C. Part of the differences between the two methods could be due to the valley-to-valley integration employed by the GC-MS method. The GCxGC testing projects the peaks to a theoretical baseline, which has been shown to over predict the mass of paraffins in the fuel [22]. This was corrected in part by using GC-MS, and for the corrected fuels the error is lower than 7.5%. While the JP-8 GC-MS surrogate still over-predicts the cloud point by approximately 17 ˚C, this is still an improvement over the literature surrogates by about 6 ˚C, and the distillation curves for the proposed surrogates are as good as the literature surrogates [14]. The distillation curves for all the predicted fuel surrogates are shown in Figure 4.15Figure 4.24.The surrogates which give the best fits are the HAD, both JP-8s, IPK, and SPK fuels. UNCLASSIFIED 75 UNCLASSIFIED In some cases, like the USD, HRJ, and SPK2, the general trend of prediction is correct, but the values are shifted by 5-20 ˚C. This shift is caused by the optimization model adding more of a compound to help match the cetane number or average molecular weight. The two JP-8 surrogates have slightly different predicted distillation curves. The GC-MS surrogate generally under predicts the boiling temperature by about 3 ˚C, while the GCxGC surrogate over predicts by about 5 ˚C during the latter half of distillation. The differences between the two surrogates are primarily due to the fact that the GCxGC testing can determine the difference between the isoparaffin and naphthenic components, while the GC-MS testing cannot. This results in different restrictions of the components that can be used in the surrogates, and how much of each component that can be used. The HRD, HRJ, and HRJ-8 fuels have particularly poor predictions. The bimodal distribution of the HRJ-8 fuel makes it harder to predict a distillation curve with two plateaus. In addition, the database does not include isoparaffins that boil between 250-350 ˚C, which compromises the ability to create high boiling surrogates comprised of high isoparaffin content. The only alkane component in the current DIPPR database that boils between 250-350 ˚C is a naphthenic compound. Because the GCxGC testing shows the HRD, HRJ, and HRJ-8 fuels are at least 70% isoparaffins, a surrogate built from the database cannot represent the high temperature portion of the distillation curve correctly. The initial boiling point is usually over predicted by 15-20 ˚C. UNCLASSIFIED 76 UNCLASSIFIED Figure 4.15. Experimental and predicted distillation curve for the JP-8 fuel and proposed surrogate developed using the GC-MS and HPLC data. Figure 4.16. Experimental and predicted distillation curve for the HAD fuel and proposed surrogate. UNCLASSIFIED 77 UNCLASSIFIED Figure 4.17. Experimental and predicted distillation curve for the USD fuel and proposed surrogate. Figure 4.18. Experimental and predicted distillation curve for the JP-8 fuel and proposed surrogate developed using the GCxGC data. UNCLASSIFIED 78 UNCLASSIFIED Figure 4.19. Experimental and predicted distillation curve for the IPK fuel and proposed surrogate. Figure 4.20. Experimental and predicted distillation curve for the HRJ fuel and proposed surrogate. UNCLASSIFIED 79 UNCLASSIFIED Figure 4.21. Experimental and predicted distillation curve for the SPK2 fuel and proposed surrogate. Figure 4.22. Experimental and predicted distillation curve for the HRJ-8 fuel and proposed surrogate. UNCLASSIFIED 80 UNCLASSIFIED Figure 4.23. Experimental and predicted distillation curve for the SPK fuel and proposed surrogate. Figure 4.24. Experimental and predicted distillation curve for the HRD fuel and proposed surrogate. UNCLASSIFIED 81 UNCLASSIFIED Detailed surrogate composition information is presented in Table 4.9-Table 4.18. Each table outlines the individual components in each surrogate. An interesting trend that can be noted by looking at the Group B surrogates is that increasing the number of potential components does not increase the model accuracy. Table 4.9. Adjusted components for the proposed JP-8 surrogate developed using the GC-MS and HPLC data. Melting Temperature (˚C) Boiling Temperature (˚C) Mass Fraction in Surrogate n-Octylbenzene -36.00 264.55 0.1480 nPropylcyclohexane -74.73 181.13 0.1945 2,2,5,5Tetramethylhexane -12.60 137.61 0.0776 3-Methylundecane -58.00 210.95 0.3896 Name Structure Table 4.10. Adjusted components for the proposed HAD surrogate. Name Melting Boiling Temperature Temperature (˚C) (˚C) Structure Mass Fraction in Surrogate n-Butylbenzene -87.85 183.46 0.1143 n-Tridecylbenzene 10.00 341.43 0.1087 Ethylcyclohexane -111.31 131.95 0.0225 n-Butylcyclohexane -74.73 181.13 0.0572 3-Methylundecane -58.00 210.95 0.1166 n-Decylcyclohexane -1.73 297.75 0.2366 2,2,4,4,6,8,8Heptamethylnonane -110.15 246.50 0.2382 UNCLASSIFIED 82 UNCLASSIFIED Table 4.11. Adjusted components for the proposed USD surrogate. Name Melting Boiling Temperature Temperature (˚C) (˚C) Structure n-Butylbenzene Mass Fraction in Surrogate -87.85 183.46 0.11672 10 341.43 0.12028 2,5Dimethylhexane -91.15 109.26 0.05842 n-Butylcyclohexane -74.73 181.131 0.046535 -1.73 297.75 0.28896 -110.15 246.5 0.24946 n-Tridecylbenzene nDecylcyclohexane 2,2,4,4,6,8,8Heptamethylnonane Table 4.12. Adjusted components for the proposed JP-8 surrogate developed using the GCxGC data. Name Melting Boiling Temperature Temperature (˚C) (˚C) Structure Mass Fraction in Surrogate n-Hexylbenzene -61.15 226.26 0.1229 nPropylcyclohexane -111.31 131.95 0.1355 n-Butylcyclohexane -74.73 181.131 0.1366 2,2,5,5Tetramethylhexane -12.60 137.61 0.0601 3-Methylundecane -58.00 210.95 0.1537 2,2,4,4,6,8,8Heptamethylnonane -110.15 246.5 0.1537 UNCLASSIFIED 83 UNCLASSIFIED Table 4.13. Adjusted components for the proposed IPK surrogate. Melting Temperature (˚C) Boiling Temperature (˚C) Mass Fraction in Surrogate n-Butylbenzene -87.85 183.46 0.0044 n-Butylcyclohexane -74.73 181.13 0.0013 3,3,5Trimethylheptane -108.15 155.83 0.2618 2,2,3,3Tetramethylhexane -54.00 160.46 0.3298 3-Methylundecane -58.00 210.95 0.3957 2,2,4,4,6,8,8Heptamethylnonane -110.15 246.50 0.0031 Melting Temperature (˚C) Boiling Temperature (˚C) Mass Fraction in Surrogate n-Butylbenzene -87.85 183.46 0.0004 n-Butylcyclohexane -74.73 181.13 0.0086 n-Decylcyclohexane -1.73 297.75 0.0086 2,5-Dimethylhexane -91.15 109.26 0.0786 4-Methyloctane -113.20 142.59 0.0786 4-Methylnonane -98.70 165.85 0.1572 3-Methylundecane -58.00 210.95 0.1965 -110.15 246.50 0.3537 -38.00 447.00 0.0236 Name Structure Table 4.14. Adjusted components for the proposed HRJ surrogate. Name Structure 2,2,4,4,6,8,8Heptamethylnonane Squalane UNCLASSIFIED 84 UNCLASSIFIED Table 4.15. Adjusted components for the proposed SPK2 surrogate. Melting Temperature (˚C) Boiling Temperature (˚C) Mass Fraction in Surrogate n-Butylbenzene -87.85 183.46 0.0006 n-Butylcyclohexane -74.73 181.13 0.0064 2-Methyloctane -80.10 143.45 0.0859 4-Methylnonane -98.70 165.85 0.3435 3-Methylundecane -58.00 210.95 0.1026 Name Structure Table 4.16. Adjusted components for the proposed HRJ-8 surrogate. Name Melting Boiling Temperature Temperature (˚C) (˚C) Structure Mass Fraction in Surrogate n-Butylbenzene -87.85 183.46 0. 0179 n-Butylcyclohexane -74.73 181.13 0. 0069 n-Decylcyclohexane -1.73 297.75 0. 0069 2,5-Dimethylhexane -91.15 109.26 0. 084 4-Methyloctane -113.20 142.59 0. 084 4-Methylnonane -98.70 165.85 0. 1681 3-Methylundecane -58.00 210.95 0. 2101 -110.15 246.50 0. 3361 -38.00 447.00 0. 0336 2,2,4,4,6,8,8Heptamethylnonane Squalane UNCLASSIFIED 85 UNCLASSIFIED Table 4.17. Adjusted components for the proposed SPK surrogate. Name Melting Boiling Temperature Temperature (˚C) (˚C) Structure Mass Fraction in Surrogate n-Butylbenzene -87.85 183.46 0.0002 n-Butylcyclohexane -74.73 181.13 0.0103 n-Decylcyclohexane -1.73 297.75 0.0104 2,5-Dimethylhexane -91.15 109.26 0.0844 4-Methylnonane -98.70 165.85 0. 1684 3-Methylundecane -58.00 210.95 0. 2109 -110.15 246.50 0. 3375 -38.00 447.00 0. 0253 Melting Temperature (˚C) Boiling Temperature (˚C) Mass Fraction in Surrogate n-Octylbenzene -61.15 226.26 0.0059 n-Decylbenzene -14.38 298.04 0.0166 nDecylcyclohexane -1.73 297.75 0.0439 2,7-Dimethyloctane -54.00 160.02 0.0576 -110.15 246.50 0.6603 -38.00 447.00 0.0512 2,2,4,4,6,8,8Heptamethylnonane Squalane Table 4.18. Adjusted components for the proposed HRD surrogate. Name Structure 2,2,4,4,6,8,8Heptamethylnonane Squalane UNCLASSIFIED 86 UNCLASSIFIED 4.4 Conclusions Ten different fuel surrogates were developed to fit both traditional petroleum fuels and a variety of alternative fuels. All surrogates used 10 or less adjusted components. Paraffins were constrained to match experimental compositions. Two different surrogates were developed for the JP-8 fuel based on two different sets of composition information. The difference in the paraffin distributions between the GC-MS and GCxGC predictions causes a difference of 5 ˚C in the cloud point prediction temperatures, with the GC-MS prediction being much closer to the experimental value. The JP-8 surrogates represent an improvement of 6˚C over the literature Jet A surrogates at predicting the cloud point temperature, and predict the distillation curve equally well [14]. The USD and HAD fuel surrogates predict all properties within 10% for the cetane number, average molecular weight, and cloud point temperature, and show an improvement of 15 and 20 ˚C respectively over the literature surrogates cloud points [15, 23, 24], while modeling the distillation curves simultaneously. The alternative fuel surrogates for IPK and SPK predict the distillation curves within 5 ˚C and cloud point temperatures within 16%. The SPK2 surrogate predicts the cloud point temperature with 12% error and gets the general shape of the distillation curve correct, but shifted by 3-5 ˚C. The HRD, HRJ and HRJ-8 fuels were difficult to represent with surrogates. These higher boiling fuels were unable to be fitted due to a lack of high molecular weight isoparaffin components in the database. It is likely that a lack of diversity in the set of potential surrogate components hindered optimization for the higher boiling fuels. UNCLASSIFIED 87 UNCLASSIFIED: Distribution Statement A. Approved for public release. CHAPTER 5 Prediction of Cetane Number for Fuels and Fuel Mixtures UNCLASSIFIED 88 UNCLASSIFIED 5.1 Introduction The cetane number of a fuel is an important metric for the evaluation of diesel fuels. The cetane number is a measure of the ignition delay in a compression-ignition engine. It is a property that is usually evaluated for jet fuels, but since jet engines do not rely on autoignition there is no minimum or maximum value for the fuel specification. However, the one fuel for the battlefield program [33-35, 81] calls for the use of JP-8 fuel in all military vehicles, including tanks and trucks. Because cetane number is not a specification for jet fuel, different batches of JP-8 as well as biofuels from different sources can have largely different cetane numbers, while still being within spec for all other properties. Knowledge of blend behavior is thus important in the field. Cetane number prediction is challenging. Ignition delay depends on fuel structural, thermodynamic, and transport properties, because the process involves droplet breakup dependent on surface tension, vaporization, mixing, and the free radical process of combustion. Common cetane number predictive correlations use other fuel properties [82]. These can include aniline point [83-87], API gravity [83-86], density [83, 85], distillation curve [32, 83-86], and viscosity [83-87]. The majority of these equations are empirical relations, with fitted parameters. They work well for other petroleum based fuels from outside the training set. However, when alternative fuels with different compositions or biofuels with oxygenated components are used, the empirical relationships can fail [25, 88]. Quantitative structure-property relationship (QSPR) methods have also been used to relate the structure of fuel compounds to the compound cetane number [89-92]. The QSPR method relates the different structural elements of individual molecules to a selected property. This can be done using either basic structural information about the molecule [89], or using UNCLASSIFIED 89 UNCLASSIFIED NMR spectra to relate peaks to the desired property [90]. This method works very well for nalkanes, but has the same amount of scatter as other methods for branched compounds. It has also only been applied to potential fuel compounds, and not to fuel surrogates or actual fuel mixtures. An alternative to comparative property prediction is to use spectroscopy to predict cetane numbers. Nuclear magnetic resonance (NMR) spectroscopy has been used to relate the ratio of aromatic and paraffin carbons to the cetane number [88]. Other groups have used carbon NMR to relate the density of the CH2 and CH3 groups to the cetane number [31]. The DeFries method works well for n-paraffin cetane number prediction, but has a large amount of scatter for isoparaffin prediction. Being able to relate the cetane number to common fuel components allows for the model to be applied to multiple fuel types. Chemometric methods can be used to relate the IR spectra to cetane number [27-30, 89, 93-95]. Chemometrics is defined by using data analysis methods to extract results from large volumes of data. In the case of cetane number, the peaks in the spectra are correlated to the cetane number using various types of regression. This work focuses on using the partial least squares (PLS) method of regression. In PLS regression, the absorbance at each wavelength or wavenumber is linearly related to the cetane number or other fuel property [89, 94, 95]. This differs from principle component analysis (PCA) which finds links and correlations between sets of data [96, 97]. In PLS, a vector of linear coefficients are developed, which can then be used to predict cetane numbers for fuels that weren’t used to create the coefficients. In this chapter, two different methods for predicting cetane number are evaluated. The first method takes a strictly composition-based approach to predicting cetane number. Two dimensional gas chromatography (GCxGC) testing is used to analyze the fuel in terms of each UNCLASSIFIED 90 UNCLASSIFIED molecular class (paraffin, isoparaffin, naphthenes, aromatic, ect.) by carbon number. The composition of each class is then used to find the volume average carbon number for each class, and a pseudocomponent of that carbon number is used to represent the individual molecular class. Then, the method of Ghosh and Jaffe [26] is used to predict the cetane number for the fuel based on the pseudocomponents. Because of the form of the cetane formula, the pseudocomponent method used here is mathematically equivalent to including each surrogate component because each class shares the same blending number. The second method for cetane prediction uses partial least squares (PLS) regression to relate the cetane number to the infrared spectra of the fuel [27-30, 93]. Both near-infrared (NIR) and Fourier Transform infrared (FTIR) spectra were used to create separate regressions. NIR and FTIR were selected for this work due to their greater portability than NMR. Neat fuels and the distillation fraction samples were used as a training set, and mixtures of the fuels were used as the testing set. A chemometric algorithm was developed that can be rapidly applied to predict cetane number based on spectroscopy. 5.2 Materials and Methods 5.2.1 Materials 5.2.1.1 Fuels The seven fuels evaluated here are previously described in Chapter 4, section 4.2.1. This chapter focuses on the jet fuels: the six group B alternative fuels plus JP-8. Two different techniques were used to develop a wider variety of cetane numbers from the test set. The fuels were distilled, and the T10, T50, and T90 fractions were tested, which resulted in different compositions and cetane numbers for each fraction. The second technique analyzed blends of UNCLASSIFIED 91 UNCLASSIFIED 15%, 35%, 50%, 65% and 85% by volume of each alternative fuel in JP-8. Derived cetane numbers for each distillation fraction and blend were measured by TARDEC in Warren, MI. 5.2.1.2 Pure compounds Six pure isoparaffin compounds were used for validation of the PLS model. These compounds were: 2-methylpentane, 3-methylpentane, 2,3-dimethylpentane, 2,4dimethylpentane, 2,2,4-trimethylpentane, and 2,2,4,4,6,8,8-heptamethylnonane. The compounds were selected because the specific isomer had a cetane number reported in the Murphy Compendium of Cetane Numbers [35], and were available from Sigma Aldrich. The cetane numbers and testing methods for the selected compounds are recorded in Table 5.1. Table 5.1. Branched compounds used for comparison, with their cetane numbers and the method used to find the cetane number. Cetane numbers and measurement methods from the Murphy Compendium [35]. Compound 2-Methylpentane 3-Methylpentane 2,3-Dimethylpentane 2,4-Dimethylpentane 2,2,4-Trimethylpentane 2,2,4,4,6,8,8Heptamethylnonane Measurement Method Blend or Correlated from Octane number Blend IQT CFR Engine Blend Direct Cetane Range (If multiple values given) Cetane number used for benchmarking 23-34 33 12-17.5 30 21 29 14 15 5.2.2 Experimental Methods 5.2.2.1 Distillation The distillation method is discussed in detail in Chapter 4, section 4.2.2.6. A two liter sample was distilled using a scaled up method similar to ASTM D86. One-hundred mL samples were collected, and the temperature was recorded after each sample collection. The sample referred to as T10 corresponds to the mixture of the 10% and 15% volume distilled samples, T50 corresponds to the mixture of the 45% and 50%, and T90 to the mixture of 85% and 90%. UNCLASSIFIED 92 UNCLASSIFIED 5.2.2.2 Derived Cetane Number Tests for derived cetane number were performed at TARDEC in Warren, Michigan. An ignition quality tester (IQT) from Advanced Engine Technologies (AET) was used. ASTM D6890 was followed for all tests (Table 5.2). The cetane numbers for all fuels are shown in Figure 5.1, and are linear with composition in JP-8. Table 5.2. Derived cetane numbers for all neat fuels, distillation fractions, and fuel mixtures. Neat Fuel Distillation Fraction T10 Distillation Fraction T50 Distillation Fraction T90 15:85 JP-8:Fuel Mixture 35:65 JP-8:Fuel Mixture 50:50 JP-8:Fuel Mixture 65:35 JP-8:Fuel Mixture 85:15 JP-8:Fuel Mixture JP-8 45.78 41.95 44.45 49.91 - IPK 31.59 31.2 30.82 30.77 33.35 36.66 39.22 41.52 43.77 HRJ 57.99 53.11 57.87 64.53 55.88 54.00 52.13 50.02 48.01 SPK2 58.11 57.30 57.80 58.38 56.21 54.07 51.73 50.01 47.68 HRJ-8 56.88 49.88 55.84 64.89 54.96 52.70 51.05 49.53 47.76 SPK 60.05 52.96 57.38 63.83 56.36 53.98 52.52 49.62 47.99 Figure 5.1. Cetane numbers for all alternative fuels in mixtures with JP-8. UNCLASSIFIED 93 HRD 64.65 55.08 67.59 73.90 62.16 58.04 55.62 52.97 48.97 UNCLASSIFIED 5.2.2.3 Two Dimensional Gas Chromatography The GCxGC method is discussed in detail in Chapter 4, section 4.2.2.4. The GCxGC results were used to select the simple psuedocomponents for cetane number predction. 5.2.2.4 Near-Infrared Spectroscopy (NIR) NIR spectroscopy was used on all fuels, distillation fractions, and fuel mixtures. A Perkin Elmer Lambda 900 UV/VIS/NIR spectrometer was used. Light absorbance was measured from 880–1580 nm in 2 nm intervals. The slit width was set to 5.00 nm with a gain of 1 and an integration time of 0.68 s. The baseline was assessed at both normal transmittance and zero transmittance. A quartz cuvette was used with a path length of 1 cm, and a blank cuvette was used as the reference. Simpson [93] was used for all NIR peak identification discussed in this section. The NIR spectra for the neat fuels are presented in Figure 5.2-Figure 5.5. Figure 5.2 shows the entire region, while Figure 5.3-Figure 5.5 focus in on the third overtone, second combination overtone, and first combination overtone respectively. There is a double peak in the range of 880–980 nm, as well as a broad peak from 980–1100 nm. The double peak at 880–980 nm represents the ratio of CH3 to CH2 groups. The peak at 980–1100 nm has not been identified. UNCLASSIFIED 94 UNCLASSIFIED Figure 5.2. NIR spectra for all neat fuels. Figure 5.3. NIR spectra of the neat fuels, with focus on the third overtone region. UNCLASSIFIED 95 UNCLASSIFIED Figure 5.4 shows a triplet peak from 1120–1280 nm. The shoulder at 1150–1160 nm is likely the aromatic CH group. The peak at 1160–1200 nm is likely caused by the paraffin CH3 group, and the 1200–1240 nm peak the paraffin CH2 group. Figure 5.5 shows the final overtone section, with a triplet peak with a shoulder. It is unclear what contributes to the shoulder at 1340–1380 nm. The peak at 1380–1400 nm is most likely the CH3 group, the peak from 1400–1420 nm the CH2 group, and the peak at 1420–1460 nm the paraffin CH group. Figure 5.4. NIR spectra of the neat fuels, with focus on the second combination overtone region. UNCLASSIFIED 96 UNCLASSIFIED Figure 5.5. NIR spectra of the neat fuels, with focus on the first combination overtone region. A few examples for the distillation fraction spectra and mixture spectra are shown here. The spectra for all fuels and mixtures can be found in Appendix E. The distillation fractions of IPK show very little difference in absorbance relative to the neat fuel (Figure 5.6). The IPK fuel has a narrower mass distribution than the other fuels. Because of this, the distillation fractions have a smaller differences in composition than other fuels, and also very little difference in derived cetane numbers as shown in Table 5.2. The HRJ-8 fuel has a bimodal distribution of molecular weight components, spread out over a larger range than the IPK fuel. Consequently, the distillation fractions exhibit a larger difference in spectra (Figure 5.7). For all fuels, the more diverse the fuel composition, the larger the difference between the T10 and the T90 spectra. UNCLASSIFIED 97 UNCLASSIFIED Figure 5.6. NIR spectra for neat IPK and its distillation fractions. Differences between spectra a difficult to see. Figure 5.7. NIR spectra for neat HRJ-8 and its distillation fractions. UNCLASSIFIED 98 UNCLASSIFIED The IPK and the JP-8 fuels have largely different compositions, so the spectra of the blends in Figure 5.8 show a strong composition dependence. The largest changes are in the 1200–1240 nm peak, and the shoulder at 1380 nm. The HRD and JP-8 are closer together in composition and the spectral differences between the mixtures can mainly be seen in the 1200– 1240 nm peak and the triplet peak at 1340–1460 nm (Figure 5.9). Figure 5.8. NIR spectra for mixtures of JP-8 and IPK. UNCLASSIFIED 99 UNCLASSIFIED Figure 5.9. NIR spectra for mixtures of JP-8 and HRD. 5.2.2.5 Fourier Transform Infrared Spectroscopy (FTIR) FTIR was performed on all fuels, distillation fractions, and fuel mixtures. A Perkin Elmer Spectrum One FTIR spectrometer was used with a universal ATR sampling unit. 2-3 drops of each sample were placed on the ATR plate, and a Teflon lined plastic vial cap was placed over the sample to reduce evaporation during scans. Four scans were taken from 4000–650 cm-1 and an ATR correction was applied. Absorbance values were reported at each wavenumber. The spectra for the neat fuels are presented in Figure 5.10-Figure 5.13. The FTIR peak identification in this section comes from Socrates [98], as well as comparison with known compounds through NIST Webbook [53]. The full range is shown in Figure 5.10, and Figure 5.11-Figure 5.13 show the three main peak absorbance regions. Figure 5.11 shows a close up view of the aromatic UNCLASSIFIED 100 UNCLASSIFIED region, Figure 5.12 shows a view of the carbon-carbon stretch region, and Figure 5.13 shows the carbon-hydrogen stretch region. Figure 5.10. FTIR spectra for all neat fuels. The first region is from 650–950 cm-1 and has a quintuple or double peak, depending on the fuel (Figure 5.11). This region shows the out of plane bending for the aromatic ring. The number of peaks reflects the substitution pattern of the aromatic rings. Most of the fuels are mono-substituted, with a peak at about 725 cm-1 and a small shoulder at 745 cm-1. The exception is JP-8, which has a quintuple peak in this region. It is possible this is a sextuplet or more split peak, but the signal-to-noise ratio in this region makes it hard to distinguish the peaks from the background noise. UNCLASSIFIED 101 UNCLASSIFIED Figure 5.11. FTIR spectra for all neat fuels, with a focus on the aromatic out of plane bending region. Figure 5.12 shows the carbon-carbon stretch region. There are two peaks in this region, a double peak from 1350–1400 cm-1, and what might be a non-resolved double peak at 1425–1500 cm-1. For most fuels except the IPK, the double peak at 1350–1400 cm-1 appears more as a single peak with a shoulder. Since the IPK likely has a much higher ratio of multi-branched to monobranched isoparaffins than the other fuels, the peak at 1350–1400 cm-1 most likely represents the side chain carbon stretch. The peak from 1425–1500 cm-1 likely represents the stretching of the backbone carbon chain. Some of the fuels show the peak “leaning” to one side or the other of the main peak. This “lean” might reflect the position of the side chains on the paraffin backbone. UNCLASSIFIED 102 UNCLASSIFIED Figure 5.12. FTIR spectra for all neat fuels, with a focus on the carbon-carbon stretch region. Figure 5.13. FTIR spectra for all neat fuels, with a focus on the carbon-hydrogen stretch region. UNCLASSIFIED 103 UNCLASSIFIED Figure 5.13 shows a close up view of the carbon-hydrogen stretching region. There is a quadruplet peak from 2825–3000 cm-1. This region contains the stretches for all the carbonhydrogen bonds, including the CH3, CH2, and CH groups in paraffins, isoparaffins, and aromatic side chains. The FTIR spectra for distillation fractions and mixtures of select fuels are presented in Figure 5.14-Figure 5.17. Figures for all fuel distillation fractions and mixtures can be found in Appendix E. As with the NIR spectra, the more diverse the composition of the neat fuel, the more differences there are in the distillation fraction and mixture spectra. However, the peaks in the FTIR spectra are sharper than the NIR spectra, and the FTIR differences are primarily in peak height and more difficult to discern with the eye than the NIR differences. Figure 5.14. FTIR spectra for neat JP-8 and its distillation fractions. UNCLASSIFIED 104 UNCLASSIFIED Figure 5.15. FTIR spectra for neat SPK2 and its distillation fractions. Figure 5.16. FTIR spectra for mixtures of JP-8 and IPK. UNCLASSIFIED 105 UNCLASSIFIED Figure 5.17. FTIR spectra for mixtures of JP-8 and SPK. 5.2.3 Calculation Methods 5.2.3.1 Cetane Number by Molecular Class via the method of Ghosh and Jaffe Cetane numbers of fuel mixtures were estimated by the method of Ghosh and Jaffe [26]. This model relates the volume fraction, vi, pure component cetane number CNi, and an adjustable parameter, i, to the mixture cetane number, CN. The general form of this equation is: 𝐶𝑁 = ∑𝑖 𝑣𝑖 𝛽𝑖 𝐶𝑁𝑖 ∑𝑖 𝑣𝑖 𝛽𝑖 ( 5.1 ) Equation ( 5.1 ) permits the lumping of fuel classes into a single pseudocomponent. This can both simplify the individual fuel surrogates, and also allow a mixture of two fuels to be represented as a binary mixture of two pseudocomponents. The strategy is justified by the form of Equation 5.1 and explained here. The pseudocomponent model is valid if the UNCLASSIFIED 106 UNCLASSIFIED pseudocomponent β value, βx, is defined as the volume average of the surrogate component βs such that: 𝑣𝑥 = ∑ 𝑣𝑖 𝑣𝑖 𝛽𝑥 = ∑ 𝛽𝑖 𝑣𝑥 𝑣𝑖 𝑣𝑥 𝛽𝑖 𝐶𝑁𝑖 𝐶𝑁𝑥 = ∑ 𝛽𝑥 ( 5.2 ) ( 5.3 ) ( 5.4 ) 𝑖 where vx is the volume fraction of the summed components assuming excess volume is zero. The Ghosh and Jaffe method is consistent with equations ( 5.2 ), ( 5.3 ), and ( 5.4 ) as proven in Appendix F, where pure component βs and cetane numbers are also tabulated. According to equations ( 5.1 ) through ( 5.4 ), all components sharing a common  value can be lumped as a single psuedocomponent using algebraic association. For example, because all paraffins have the same β, the paraffin portion of the fuel surrogate can be represented by one representative paraffin pseudocomponent calculated by the volume fraction weighted cetane number. The cetane number for compounds other than n-paraffins can depend on the isomerization of the compounds, and there can be multiple components with the same carbon number and class that have significantly different cetane numbers. For this reason, the pseudocomponents used for these surrogates have cetane numbers that are purely based on the carbon number. This work makes a distinction between mono-isoparaffins (single-branched) and multi-isoparaffins (multi-branched) for the determination of the surrogate cetane number because the cetane number tends to decrease with branching for a given carbon number. A monoisoparaffin has one tertiary carbon, and a multi-isoparaffin has multiple tertiary or quaternary carbons, but even this distinction does not result in clear trends. Cetane numbers for various UNCLASSIFIED 107 UNCLASSIFIED compounds from the Murphy Compendium of Experimental Cetane Number Data [61] are plotted against their carbon numbers in Figure 5.18. The isoparaffin cetane numbers scatter significantly, with an overall trend of increasing with carbon number. Selecting an exact cetane number for a given carbon number is impossible because cetane numbers for a given carbon number can vary by nearly 40 cetane numbers. The Ghosh and Jaffe paper has a discrepancy between the cetane numbers listed for their mono- and multi-isoparaffins in the table and those shown on their figures. Because a fuel is a composite of many branched compounds, some averaging is justified. An arbitrary line was placed through both the mono- and multi-isoparaffins, as well as an average line between the two. For the fuel surrogates presented here, either the plotted line for the mono-isoparaffins are used, or a 50-50 mixture of the mono-iso- and multi-isoparaffins are used, depending on what gives the best fit for the neat fuel. Thus, the neat fuel cetane numbers were correlated instead of predicted. Blends were then predicted based on the fitted neat fuel cetane numbers. UNCLASSIFIED 108 UNCLASSIFIED Figure 5.18. Plot of mono- and multi-isoparaffins, with representative lines for the values used in the cetane number prediction. The cetane number that is used to represent each lump could either be calculated using a volume average of the cetane numbers of the individual paraffin components, or the average carbon number can be calculated using the volume fractions, and that carbon number component can represent the lump. For the isoparaffin lump, because a linear approximation is used, the two methods give the same representative cetane number. For the paraffin lump, changing the method resulted in a 0.02-0.2% change in the prediction of the neat fuel cetane number, depending on the total paraffin content of the fuel. Because these changes were negligible, the simpler method of using the volume average carbon number was used for the pseudocomponent selection. An overall fuel pseudocomponent is then built from pseudocomponents for each of the 7 molecular classes: n-paraffins, isoparaffins, monocycloparaffins, dicycloparaffins, mononuclear UNCLASSIFIED 109 UNCLASSIFIED aromatics, polynuclear aromatics, and monocycloaromatics. The method is algebraically consistent with writing the sum for each component using equation ( 5.1 ), but provides a simple way to represent the fuel as a single component when fuels are blended assuming excess volume is zero. 5.2.3.2 Cetane Number Using Partial Least Squares Regression Experimental cetane values were correlated with NIR and FTIR spectra using Partial Least Squares (PLS) regression. To provide wide dissemination of the resulting correlation, the program R was used to do the regression [99]. R was selected because it is an open source program that can be downloaded and run on any computer, and has a PLS regression package built into the program. Code was written for R to process and regress the data in a consistent manner, and is presented in Appendix G. The cetane number and absorbance for each sample was fed to the program, and the data was broken up into a training set and a testing set. The training set was regressed to develop the regression coefficients, and those coefficients were used to calculate the cetane numbers for the data not used in the development of the coefficients. This allows the model to be tested with fuels that were not directly used in the regression. The neat fuels and distillation fractions were used as the training set, and the blends were used as the predicted (testing) set. A section of the code permits the removal of baseline sections if desired. An important factor in the regressions is the number of regression components. Within the PLS regression, the term “components” refers to the number of times that the error is reregressed. The R statistical package “pls” combines these regression components into one set of regression coefficients. The coefficients are then multiplied against the absorbance values at each wavelength or wavenumber and summed to give the predicted cetane number. The resulting equation looks like: UNCLASSIFIED 110 UNCLASSIFIED 𝐶𝑁𝑝𝑟𝑒𝑑𝑖𝑐𝑡𝑒𝑑 = ∑ 𝑏𝑛 𝐴𝑛 + 𝐶 ( 5.5 ) 𝑛 where CNpredicted is the predicted cetane number, n is the number of wavelengths or wavenumbers at which there is a measured absorbance, bn is the regressed coefficient at the nth wavelength or wavenumber, and An is the absorbance at the nth wavelength or wavenumber, C is the intercept. The data was mean centered during processing. Mean centering of the data allows for the intercept C to be calculated. Mean centering is performed on the training and testing sets separately. To mean center the data, the absorbance at each wavelength is averaged, then the average is subtracted from the absorbance of each sample at that wavelength. The experimental cetane numbers are mean centered as well. The intercept is calculated using the values from the training set, where the experimental cetane numbers are known. This allows the model to be extended to datasets that were not used in the training set. The intercept is calculated using the average values for the cetane number and the absorbance at each wavelength such that: 𝐶 = 𝐶𝑁𝑡𝑟𝑎𝑖𝑛,𝑚𝑒𝑎𝑛 − ∑ 𝑏𝑛 𝐴𝑡𝑟𝑎𝑖𝑛,𝑚𝑒𝑎𝑛,𝑛 ( 5.6 ) 𝑛 where CNtrain,mean is the average of the cetane numbers in the training set, and Atrain,mean,n is the average of the absorbance values of the training set at the nth wavelength or wavenumber. For each type of spectroscopy, the testing set used the same coefficients and intercepts developed from the training set. 5.3 Results and Discussion 5.3.1 Ghosh and Jaffe Method The predicted cetane numbers are compared to the experimental data in Figure 5.19- Figure 5.24. In these figures, the experimental data is presented as points with error bars for the UNCLASSIFIED 111 UNCLASSIFIED experimental error, and the predicted cetane numbers are presented as curves. The curves are named based on the assumption that was made about the branching of the isoparaffin pseudocomponent that was used. In Figure 5.19, the selected isoparaffin component for both the JP-8 and IPK fuels are adjusted simultaneously. For Figure 5.20-Figure 5.24, the optimum JP-8 surrogate was used, and only the isoparaffin component of the alternative fuel was adjusted. The cetane numbers used for both groups are presented in Appendix F. Figure 5.19. Predictions of the cetane number of mixtures of JP-8 and IPK compared to experimental results. Both the JP-8 and IPK fuels are varied to find the correct isoparaffin content and carbon number. UNCLASSIFIED 112 UNCLASSIFIED Figure 5.20. Predictions of the cetane number of mixtures of JP-8 and HRJ compared to experimental results. Figure 5.21. Predictions of the cetane number of mixtures of JP-8 and SPK2 compared to experimental results. UNCLASSIFIED 113 UNCLASSIFIED Figure 5.22. Predictions of the cetane number of mixtures of JP-8 and HRJ-8 compared to experimental results. Figure 5.23. Predictions of the cetane number of mixtures of JP-8 and SPK compared to experimental results. UNCLASSIFIED 114 UNCLASSIFIED Figure 5.24. Predictions of the cetane number of mixtures of JP-8 and HRD compared to experimental results. In some cases, the composition of the neat alternative fuel needed to be adjusted to get a better fit. This is noted by, for example, “+2C”, which indicates that a compound is selected for the surrogate that two carbons larger than the average determined by the GCxGC testing. The initial surrogates and the concentrations of each functional class are presented in Table 5.3 and Table 5.4. UNCLASSIFIED 115 UNCLASSIFIED Table 5.3. Surrogate concentrations and components as determined by averaging the GCxGC results (JP-8, IPK, HRJ, and SPK2). Class n-Paraffin Iso-Paraffin MonoCyclohexane DiCyclohexane Alkylbenzene Diaromatic Cycloaromatic JP-8 Carbon Vol Number % 12 23.44 12 34.65 IPK Carbon Vol Number % 11 0.36 11 98.97 HRJ Carbon Vol Number % 12 12.46 13 85.87 SPK2 Carbon Vol Number % 10 44.1 11 55.25 12 20.99 10 0.11 10 1.53 10 0.52 12 5.07 14 0.04 12 0 10 0.05 11 12 12 11.09 1.24 3.29 10 12 10 0.38 0.02 0.12 10 12 10 0.04 0.01 0.02 10 10 10 0.05 0 0 Table 5.4. Surrogate concentrations and components as determined by averaging the GCxGC results (HRJ-8, SPK, and HRD) HRJ-8 SPK HRD Carbon Number Vol % Carbon Number Vol % Carbon Number Vol % n-Paraffin Iso-Paraffin 12 13 10.67 88.01 11 13 10.67 88.01 16 16 19.78 74.05 MonoCyclohexane 9 1.21 9 1.21 14 2.66 DiCyclohexane 11 0.08 12 0.08 13 1.49 Alkylbenzene Diaromatic Cycloaromatic 10 10 10 0 0 0 10 10 10 0 0 0 14 14 13 2.03 2.03 0.85 Class For all fuels except the IPK fuel, the mono-isoparaffin numbers give better pseudocomponent results than the multi-isoparaffins. For all fuels, the method is able to model the cetane number of the neat fuels within experimental error, once an isoparaffin component with a higher carbon number is substituted for the neat fuels. The adjustment however, compromises the intention of providing a predictive method. For all fuels except the IPK (Figure 5.19), the adjusted carbon number pseudocomponents can predict the cetane number across all mixture concentrations to within experimental error as well. Because there is high uncertainty in UNCLASSIFIED 116 UNCLASSIFIED the correct values to use for the pure component cetane numbers, the effect of a shift of one or two carbon numbers in the surrogate compositions was explored. The only fuel that needs a shift of more than two carbon numbers is the SPK2 fuel (Figure 5.21). The model always predicts curvature not present in the experiments. Representation of curvature was intended by the model authors, but the fuels considered in this work do not exhibit curvature. The curvature occurs when beta values for the components are different. Curvature is dependent on the difference in  values for the pseudocomponents. The larger the difference between the two pseudocomponent β values, the larger the curvature. A possible cause of the curvature is incorrect beta values for the isoparaffins. The fuels in this study had a much higher isoparaffin content than most petroleum fuels, and thus may be quite different than refinery fuels used to develop the  values. The GCxGC results shown in Chapter 4, section 4.2.2.4, show a distinct difference in the diversity of components represented in alternative fuels versus petroleum fuels. Because of these large composition differences, it is possible that re-regressing the beta values with more alternative fuels would result in a better prediction for the alternative fuels. The surrogates presented here have significantly less error than the surrogates developed in Chapter 4. The error in the cetane number for those surrogates was approximately 20%, while the neat fuels here have <4% error for the final surrogates. However, the surrogates developed here have highly inaccurate cloud point temperatures and distillation curves. These surrogates work well for the prediction of neat fuel cetane number, and the cetane number of blends, but do not extend to any additional fuel properties. UNCLASSIFIED 117 UNCLASSIFIED 5.3.2 NIR and FTIR Predictions 5.3.2.1 Fuel Prediction The results of the PLS regressions are shown in Figure 5.25-Figure 5.28. In the case of the NIR results, reducing the number of regression components from 25 to seven results in an increase for the training set error from 0.00027% to 1.0% (Table 5.5). However, the reduction improves the prediction in the testing set from 7.3% error to 0.96% (Table 5.5). This is due to the regression coefficients becoming overly specific to the training set when too many regression components are used. The improvement in the training set prediction outweighs the small increase in error for the training set. Figure 5.25. Predicted vs experimental cetane number for all fuels, distillation fractions, and mixtures. Prediction done using the NIR data and 25 regression components. UNCLASSIFIED 118 UNCLASSIFIED Figure 5.26. Predicted vs experimental cetane number for all fuels, distillation fractions, and mixtures. Prediction done using the NIR data and seven regression components. Figure 5.27. Predicted vs experimental cetane number for all fuels, distillation fractions, and mixtures. Prediction done using the FTIR data and 26 regression components. UNCLASSIFIED 119 UNCLASSIFIED Figure 5.28. Predicted vs experimental cetane number for all fuels, distillation fractions, and mixtures. Prediction done using the FTIR data and seven regression components. For the FTIR data, both the training and testing set errors are larger for seven regression components compared to 25 but the correlation is still within the experimental uncertainty of the measurements, which is typically ±2-3% of the cetane number. The training set error goes from essentially 0% to 0.55%, and the testing set error increases by an almost negligible amount (0.078%, Table 5.5). This is a small increase in error, (Figure 5.27 and Figure 5.28). The larger amount of noise in the baseline of the FTIR data could explain this increase. Another explanation for the increase of error in the testing set could be that there are not enough data sets in the training set, to allow for a decrease in regression components. UNCLASSIFIED 120 UNCLASSIFIED Table 5.5 Average absolute error for the prediction of cetane number using NIR and FTIR data, with varying regression conditions. The testing set and training set are evaluated separately. Spectra Type Number of Regression Components Data Set Training set Testing set Training set Testing set Training set Testing set Training set Testing set 25 NIR 7 26 FTIR 7 Average Absolute % Error 0.00027% 7.3% 1.0% 0.96% 8.3E-13% 0.85% 0.55% 0.92% The regression coefficients for the NIR regression with seven regression components and both FTIR regressions are plotted against the spectra in Figure 5.29-Figure 5.30. The NIR coefficients appear to relate directly to the spectra peaks, while the FTIR coefficients don’t appear to relate to the spectra. If the weights (coefficients multiplied by absorbance) are inspected instead (Figure 5.31 and Figure 5.32), both the NIR and FTIR coefficient peaks line up well with the spectral peaks. Fodor [27, 29] and DeFries [31] both discuss removal of the baseline, but determine it has no effect on the result beyond decreasing computing time. They do not, however, report the relationship between the regression coefficients and peaks. While the weighted coefficients have a good agreement with the peaks, there are also areas with higher weights in the background region. For this reason, portions of the NIR and FTIR baselines were manually removed, and the data was re-regressed evaluate the effect on cetane prediction. UNCLASSIFIED 121 UNCLASSIFIED Figure 5.29. Regression coefficients plotted against the NIR spectra for the JP-8, IPK, and HRD fuels for the regression using 7 components. Figure 5.30. Regression coefficients plotted against the FTIR spectra for the JP-8, IPK, and HRD fuels for the regression using 26 components. UNCLASSIFIED 122 UNCLASSIFIED Figure 5.31. Regression weights plotted against the NIR spectra for the JP-8, IPK, and HRD fuels for the regression using 7 components. Figure 5.32. Regression weights plotted against the FTIR spectra for the JP-8, IPK, and HRD fuels for the regression using 26 components. UNCLASSIFIED 123 UNCLASSIFIED Visual inspection was used to remove sections of the baseline for the data sets. The NIR sets use the data in the range of 1100-1680nm, and the FTIR data uses the ranges of 650900cm-1, 1340-1520cm-1, 2800-3000cm-1. Removing the baseline results in increased errors for all training sets except the seven component NIR set, both FTIR testing sets and the NIR testing set using seven regression components (Table 5.6). In the two cases where the baseline removal improves the error, the improvement is by less than 0.1%. Table 5.6. Average absolute error for the prediction of cetane number using NIR and FTIR data, with varying regression conditions and the baseline removed. The testing set and training set are evaluated separately. Spectra Type Number of Regression Components Data Set Training set Testing set Training set Testing set Training set Testing set Training set Testing set 25 NIR 7 26 FTIR 7 Average Absolute % Error 0.00025% 6.4% 1.0% 1.4% 1.2E-9% 1.4% 0.73% 1.1% The NIR weights did not changed dramatically; while the absolute values were changed, the positions and the peaks which were emphasized did not (Figure 5.33). Even with the baseline removed, the FTIR weights were still noisy (Figure 5.34). The removal of the baseline did increase the absolute values of the weights, and helped to separate the aromatic region from the baseline. But, because the errors were increased overall, the removal of the baseline only helped to decrease processing time, not model accuracy. The coefficients and intercepts for all PLS regressions with 7 components can be found in Appendix G. UNCLASSIFIED 124 UNCLASSIFIED Figure 5.33. Regression weights plotted against the NIR spectra of JP-8, IPK, and HRD fuels with the baseline removed. Figure 5.34. Regression weights plotted against the FTIR spectra for the JP-8, IPK, and HRD fuels for the regression using 7 components, with the baseline removed during regression. UNCLASSIFIED 125 UNCLASSIFIED 5.3.2.2 Model Validation Using Pure Components To check the accuracy of the NIR and FTIR models in prediction the cetane number of fuels or compounds outside of the fuel set, six pure isoparaffin components were tested. These six compounds, which were detailed in Table 5.1, were tested with both the NIR and FTIR methods, and the spectra were predicted using the same coefficients as developed during previous regressions. The average absolute errors are shown in Table 5.7. The average error of the compounds is very high for all models, with no regression method giving an average error lower than 33% Table 5.7. Errors for the regression of various branched compounds using the neat fuels and distillation fractions as the training set. Spectra Type Number of Regression Components 25 NIR 7 26 FTIR 7 Background Sections Removed No Yes No Yes No Yes No Yes Average Absolute % Error 112% 141% 34% 40% 41% 33% 45% 61% The individual predictions are shown along with the training set in Figure 5.35. Three of the compounds are predicted well, while the other three compounds are predicted poorly. 2,2,4,4,6,8,8-Heptamethylnonane (CN = 15) is predicted almost exactly. This compound is used as a calibration standard in the cetane number testing method, so it is a set value. The experimental value for 2,4-dimethylpentane (CN = 29) was found using direct measurement, and it is predicted within experimental error. 2-Methylpentane (CN = 33) was also predicted within experimental error. UNCLASSIFIED 126 UNCLASSIFIED Part of the error for the poor predictions could come from the literature value. For both 2,2,4-trimethylpentane (CN = 14) and 3-methylpentane (CN = 30), the reported cetane number is from the experimental testing of a blend, and the pure cetane number is extrapolated from the blending data (Table 5.1). The extrapolation of a neat compound’s cetane number from blending data can be unpredictable, depending on the curvature or lack of curvature displayed, so it is possible that the reported cetane value is incorrect. For 2,3-dimethylpentane (CN = 21), the Murphy Compendium states the value is from an IQT test, but gives no information about the purity of the sample, or a literature reference. Figure 5.35. Predicted vs experimental cetane numbers of various branched compounds regressed using the neat fuels and distillation fractions as the training set. 5.4 Conclusions Six different alternative fuels were examined in mixtures with JP-8 to evaluate methods for predicting their cetane numbers. The mixture cetane numbers were found experimentally to UNCLASSIFIED 127 UNCLASSIFIED be linear with composition. The neat fuels were also distilled to create a larger training set for the IR prediction methods. Two different methods of predicting the cetane number of mixtures were examined. The first used the method of Ghosh and Jaffe. GCxGC data were averaged to develop a surrogate composition, which was in turn used to create a pseudocomponent for the neat fuel. The pseudocomponent was then used to predict the cetane number of mixtures. If the carbon number of the isoparaffins was increased, the model could represent the paraffin number to within experimental error for all neat fuels, and all fuel mixtures except for the JP-8:IPK mixture. The Ghosh and Jaffe method did not work well for the fuels with high isoparaffin content because the model predicts incorrect curvatures. It is likely that the prediction would improve if the beta values were to be re-regressed with a more diverse range of fuels. This kind of re-regression requires a very large data set, and the data set for this study was too limited for good optimization. The second method used partial least squares regression to relate the NIR or FTIR spectra to the cetane number of the fuel mixture. The FTIR spectra give a smaller overall error for the training sets than the NIR spectra. However, the FTIR regressions do not improve when the number of regression components is decreased. When the baseline was removed, the error increased for all data sets except for the NIR regression with 25 components, while still remaining within the experimental uncertainty. The NIR spectra gave coefficients and weights that were more correlated to the spectral peaks than the FTIR regression weights. While the FTIR regression gives an average error of 0.42%, and the NIR one of 1.0%, the NIR coefficients are more clearly related to the spectra. The NIR regression is the best cetane model out of those evaluated for these fuels. The models were evaluated using six pure isoparaffins with data from UNCLASSIFIED 128 UNCLASSIFIED the literature. The three isoparaffins were predicted well, including heptamethylnonane, which is used as the calibration standard for cetane number testing. UNCLASSIFIED 129 CHAPTER 6 Isothermal Compressibility of Fuels and Fuel Mixtures 130 6.1 Introduction The compressibility of a fuel affects its behavior during the injection. As the fuel is subjected to high pressures during injection, the compressibility determines how much the volume will change. This can affect the mass injected, and is related to the optimum injection timing [100, 101]. In addition, biofuels can have a substantially different (usually lower) compressibility than petroleum fuels [101, 102]. Prediction of compressibility for different fuels based on a representative surrogate composition will allow the engine operation to be optimized. Prediction of fuel compressibility using surrogate compositions is unexplored in the literature. The adiabatic compressibility is more relevant than the isothermal compressibility because the compression process is so fast that there is little time for heat transfer. In general adiabatic and isothermal compressibility are defined as: 1 𝜕𝑉 𝜅𝑆 = − ( ) 𝑉 𝜕𝑃 𝑆 1 𝜕𝑉 𝜅𝑇 = − ( ) 𝑉 𝜕𝑃 𝑇 ( 6.1 ) ( 6.2 ) where κS is adiabatic compressibility, P is pressure, V is the molar volume, S is entropy, κT is the isothermal compressibility, and T is temperature. The adiabatic compressibility is related to the isothermal compressibility and heat capacity as shown by Tyrer [103] and Liebenberg, Mills, and Bronson [104]: 𝑇𝑉𝛼𝑃2 𝜅𝑇 − 𝜅𝑆 = 𝐶𝑃 ( 6.3 ) where CP is the constant pressure heat capacity, CV is the constant volume heat capacity, T is the temperature, and αP is the thermal expansion coefficient. For calculation of the adiabatic compressibility, this equation requires accurate isothermal compressibility, thermal expansion coefficient, and heat capacity. Experimental data for either fuel components or the fuels 131 themselves are usually available for all of these properties. Kleppa [105] gives a similar equation to ( 6.3 ) with a slightly different form. However, it requires iteration on κT. It is impractical to measure adiabatic compressibility directly using P-V-T experiments, and difficult to measure isothermal compressibility accurately. The adiabatic compressibility is most accessible via measurement of speed of sound. The relationship between adiabatic compressibility and speed of sound, v, is: 𝜈2 = 1 𝜌𝜅𝑆 ( 6.4 ) Isothermal compressibility can be related to the speed of sound via equations ( 6.3 ) and ( 6.4 ) and can be rearranged as the isothermal derivative of density with pressure: 1 𝜕𝑉 1 𝜕𝜌 1 𝑇𝑉𝛼𝑃2 𝜅𝑇 = − ( ) = ( ) = + 𝑉 𝑑𝑃 𝑇 𝑉𝜌2 𝜕𝑃 𝑇 𝜌𝑣 2 𝐶𝑃 𝜕𝜌 1 𝑇𝛼𝑃2 ( ) = 2+ 𝜕𝑃 𝑇 𝑣 𝐶𝑃 ( 6.5 ) Some groups ([106] and [38]) integrate equation ( 6.5 ) to represent the density change as a function of the speed of sound. 𝑃 𝑃 2 𝑑𝑃 𝛼𝑃 + 𝑇 ∫ 𝑑𝑃 2 𝑃0 𝜈 𝑃0 𝐶𝑃 𝜌(𝑃, 𝑇) = 𝜌(𝑃0 , 𝑇) + ∫ ( 6.6 ) where P0 is the reference pressure. Dzida et al. [106] and Payri [38] et al. ignore the second term, claiming it only accounts for a small error in density. Both have been used to predict the density of biodiesel or diesel systems. When predicting adiabatic compressibility of mixtures, there are two routes that can be taken: empirical constants can be fitted to the experimental data, or an equation of state model can be used. Both Fort and Moore [40] and Jain et al. [39] use very simple empirical parameters to fit the excess compressibility of binary mixtures to concentration: 132 𝜅𝑆𝐸 = 𝐴 + 𝐵(𝑥1 − 𝑥2 ) + 𝐶(𝑥1 − 𝑥2 )2 𝑥1 𝑥2 ( 6.7 ) where κSE is the excess adiabatic compressibility, x1 is the mole fraction of component 1, x2 is the mole fraction of component 2, and A, B, and C are fitted empirical constants. Wang and Nur [107] and Berryman [108] use a volumetric average and assumption of ideal mixing to fit the speed of sound. Multiple research groups fit the density data of mixtures using the Tait equation [36, 37, 106]. They can obtain fits with <3% error this way, however, each mixture composition is fit with different parameters, so the method is not predictive to other types or batches of fuels. Instead of using empirical relations, an equation of state (EOS) model can be used to predict the density and speed of sound. An EOS model directly relates the density of the fluid to the physical properties of the solution, which allows for mixtures to be predicted from pure component properties using mixing rules that have already been developed. Application of an EOS to fuels requires the development of a surrogate to represent the fuel. Surrogates are mixtures of components with well-known properties that, when mixed with the correct ratios, give good predictions of the physical properties of the fuel. Many equations of state have been used to model speed of sound. The Tao-Mason EOS is similar to the van der Waals EOS with a slight modification to the attractive term [109-111]. This equation can produce predictions with <3% in density for heavy alkanes, but it still has four parameters which would need to be fitted. Another option is the Goharshadi-Morsali-Abbaspour (GMA) EOS [112-114]. It has higher accuracy than Tao-Mason, with errors of <1% for alkanes, but has six fitted parameters, with no mixing rules. The Peng-Robinson equation has been used to predict the speed of sound in the gas phase, but it gives large errors for liquid densities [46]. The GERG equation can fit the data to within experimental uncertainties for multiple physical 133 properties [115]. It is a highly complicated equation that requires a large number of data sets to accurately predict mixture properties. While it is too complicated for easy use for surrogate mixtures, it is used in the REFPROP program [116] available from NIST and can be used to predict properties of binary mixtures of limited components for model evaluation. The statistical associating fluid theory (SAFT) model has many different modifications, some of which have been used to predict speed of sound and compressibility. The original SAFT overpredicts the speed of sound by approximately 30%, and PC-SAFT underpredicts speed of sound by 10-15% [42, 43]. PC-SAFT can be improved by re-regressing the constants for alkanes using the speed of sound as well as density and saturated pressure data, and the errors are reduced to an average of 3.26% [43]. Soft-SAFT gives variable predictions depending on the type of fluid. Perfluoroalkanes were modeled very well [117], but alkanes and alcohols had higher errors [44]. Finally, the SAFT-BACK equation is reported to give very good prediction of speed of sound data for pure fluids [45, 46, 118, 119]. The BACK modification uses a hard convex body form for the radial distribution function. The SAFT-BACK equation is also referred to as the SAFT-CP equation [120], but the two equations are identical except for the parameter values used. Two equations, the SAFT-BACK equation, and the Elliott, Suresh, Donohue (ESD) equation, have been chosen for evaluation. 6.2 Calculation Methods 6.2.1 ESD equation of state The ESD equation of state (EOS) is a cubic equation which more accurately predicts the repulsive forces than the van der Waals form used in the popular Peng-Robinson equation. The 134 ESD equation has been used to predict phase equilibria of pure fluids and mixtures, including polymers, but has not been used for speed of sound prediction. The ESD is defined as: 𝑍= 𝑃𝑉 2.1𝑞𝜂 1.9𝜂 9.5𝑞𝑌𝜂 = 1+ + − 𝑅𝑇 1 − 1.9𝜂 1 − 1.9𝜂 1 + 1.7745𝑌𝜂 ( 6.8 ) Where Z is compressibility factor, V is volume, P is pressure, R is the gas constant, T is temperature, and q, η, and Y are all variables defined as: 𝑞 = 1 + 1.90476(𝑐 − 1) 𝑏 𝜂= 𝑉 𝑌 = 𝑒 𝜖/(𝑘𝑇) + 1.0617 ( 6.9 ) ( 6.10 ) ( 6.11 ) where b is the packing fraction, c is a shape parameter which gives an indication of how spherical the molecule is, ε is the attractive term, and k is the Boltzmann constant. The ESD is a cubic equation; it can be re-arranged to have the form of a third degree polynomial in Z. This is beneficial because it is simpler to calculate the roots of a polynomial, rather than solve a more complex equation. When the ESD is rearranged is has the form: 𝑍 3 + 𝑎2 𝑍 2 + 𝑎1 𝑍 + 𝑎0 = 0 𝑏𝑃 𝑏𝑃 𝑎2 = −1.9 + 1.7745𝑌 −1 𝑅𝑇 𝑅𝑇 𝑏𝑃 2 𝑏𝑃 𝑏𝑃 𝑏𝑃 𝑎1 = −3.37155𝑌 ( ) − 1.7745𝑌 − 2.1𝑞 + 9.5𝑞𝑌 𝑅𝑇 𝑅𝑇 𝑅𝑇 𝑅𝑇 𝑏𝑃 2 𝑎0 = −21.77645𝑞𝑌 ( ) 𝑅𝑇 ( 6.12 ) ( 6.13 ) ( 6.14 ) ( 6.15 ) 6.2.2 SAFT-BACK equation of state The SAFT-BACK EOS [118] is a modification of the statistical associating fluid theory (SAFT) model. The SAFT equation [121] calculates the Helmholtz energy of a fluid by considering individual segments of molecules. Molecules are broken up into spherical sites, whose energy is independently calculated. These sites are then assembled into chains, and the energy of the assembly is calculated. The sites are then allowed to interact with both other sites 135 on the same molecule, and sites on other molecules, and this contribution is added to the energy. SAFT-BACK differs from SAFT in the way that the individual spheres are modeled, and this modification allows good prediction all the way to the critical point for pure n-alkanes [45, 118]. It is defined in terms of the Helmholtz energy: 𝐴𝑟𝑒𝑠 = 𝐴 − 𝐴𝑖𝑑𝑒𝑎𝑙 = 𝐴ℎ𝑐𝑏 + 𝐴𝑐ℎ𝑎𝑖𝑛,ℎ𝑐𝑏 + 𝐴𝑑𝑖𝑠 + 𝐴𝑐ℎ𝑎𝑖𝑛,𝑑𝑖𝑠 𝐴ℎ𝑐𝑏 𝛼2 𝛼 2 − 3𝛼 = 𝑚[ − − (1 − 𝛼 2 ) ln(1 − 𝜂) − 3𝛼] (1 − 𝜂)2 𝑅𝑇 1−𝜂 𝐴𝑐ℎ𝑎𝑖𝑛,ℎ𝑐𝑏 = (1 − 𝑚) ln[𝑔ℎ𝑐𝑏 (𝑑)] 𝑅𝑇 𝐴𝑑𝑖𝑠 𝑢 𝑖 𝜂 𝑗 = 𝑚 ∑ ∑ 𝐷𝑖𝑗 ( ) ( ) 𝑅𝑇 𝑘𝑇 𝜏 𝑖 𝐴 𝑐ℎ𝑎𝑖𝑛,𝑑𝑖𝑠 𝑅𝑇 =𝑚 𝜆𝐴 𝑐ℎ𝑎𝑖𝑛,ℎ𝑐𝑏 𝐴ℎ𝑐𝑏 ( 6.16 ) ( 6.17 ) ( 6.18 ) ( 6.19 ) 𝑗 𝐴𝑐ℎ𝑎𝑖𝑛,ℎ𝑐𝑏 𝐴𝑑𝑖𝑠 𝑢 𝜂 𝑅𝑇 ) ( 𝑅𝑇 ) ∑ ∑ 𝐷𝑖𝑗 ( ) ( ) = 𝜆 𝐴ℎ𝑐𝑏 𝑘𝑇 𝜏 ( 𝑅𝑇 ) 𝑖 𝑗 𝑖 𝑗 ( ( 6.20 ) where A is Helmholtz energy, k is the Boltzmann constant, m is the number of segments per molecule, α is a shape parameter, η is the packing fraction, ghcb(d) is the hard convex body radial distribution function, Dij is a power law expansion coefficient [122], u is the attractive energy, and τ is a constant. η, ghcb(d), u, and τ are defined as: 1 𝑏 𝜂 = 𝜋𝑁𝐴𝑉 𝑚𝜌𝑑 3 = 𝑏𝜌 = 6 𝑉 1 3𝛼(1 + 𝛼)𝜂 2𝛼 2 𝜂2 ℎ𝑐𝑏 (𝑑) 𝑔 = + + 1 − 𝜂 (1 − 𝜂)2 (1 + 3𝛼) (1 − 𝜂)3 (1 + 3𝛼) 𝜀 𝑢 = 𝑢0 (1 + ) 𝑘𝑇 √2𝜋 𝜏= 6 ( 6.21 ) ( 6.22 ) ( 6.23 ) ( 6.24 ) where ρ is the molar density, b is a constant which lumps together the other constants in equation ( 6.21 ), NAV is Avogadro’s number, d is the temperature dependent segment diameter, and ε is a constant which determines the temperature dependence of u. ε/k is set to 1 for small molecules 136 (methane in this work) and 10 for larger molecules. σ is the temperature independent segment diameter. d and σ are defined as: 𝑑 = 𝜎[1 − 0.12𝑒^(−3𝑢0 /(𝑘𝑇)] 1 𝑣 00 = 𝜋𝑁𝐴𝑉 𝜎 3 6 ( 6.25 ) ( 6.26 ) Helmholtz energy is related to pressure by: 𝐴 𝜕 (𝑅𝑇) 𝑍 = 𝜌( ) 𝜕𝜌 𝐴 𝜕 (𝑅𝑇) = 𝜂( ) 𝜕𝜂 𝑁,𝑇 ( 6.27 ) 𝑁,𝑇 This makes Z a function of each of the five terms used to make up the Helmholtz energy such that: 𝑍 = 𝑍 𝑖𝑑𝑒𝑎𝑙 + 𝑍 ℎ𝑐𝑏 + 𝑍 𝑐ℎ𝑎𝑖𝑛,ℎ𝑐𝑏 + 𝑍 𝑑𝑖𝑠 + 𝑍 𝑐ℎ𝑎𝑖𝑛,𝑑𝑖𝑠 𝑍 𝑖𝑑𝑒𝑎𝑙 = 1 (1 + 3𝛼)𝜂 + (3𝛼 2 − 3𝛼 − 2)𝜂2 + (1 − 𝛼 2 )𝜂3 𝑍 ℎ𝑐𝑏 = 𝑚 [ ] (1 − 𝜂)3 (1 − 𝑚)𝜂 3𝛼(1 + 𝛼) 𝑍 𝑐ℎ𝑎𝑖𝑛,ℎ𝑐𝑏 = [1 + (1 − 𝜂) (1 + 3𝛼)(1 − 𝜂)2 𝑔ℎ𝑐𝑏 (𝑑) 4𝛼 2 𝜂 + ] (1 + 3𝛼)(1 − 𝜂)3 𝑔ℎ𝑐𝑏 (𝑑) 𝑢 𝑖 𝜂 𝑗 𝑍 𝑑𝑖𝑠 = 𝑚 ∑ ∑ 𝑗𝐷𝑖𝑗 [ ] [ ] 𝑘𝑇 𝜏 𝑖 𝑍 𝑐ℎ𝑎𝑖𝑛,𝑑𝑖𝑠 ( 6.28 ) ( 6.29 ) ( 6.30 ) ( 6.31 ) ( 6.32 ) 𝑗 𝑑𝑖𝑠 𝐴𝑐ℎ𝑎𝑖𝑛,ℎ𝑐𝑏 𝑑𝑖𝑠 𝐴 𝐴𝑐ℎ𝑎𝑖𝑛,ℎ𝑐𝑏 𝐴𝑑𝑖𝑠 ℎ𝑐𝑏 𝑐ℎ𝑎𝑖𝑛,ℎ𝑐𝑏 = 𝜆[ 𝑍 + ℎ𝑐𝑏 𝑍 − 𝑍 ] (𝐴ℎ𝑐𝑏 )2 𝐴ℎ𝑐𝑏 𝐴 ( 6.33 ) 6.2.3 Calculation method The relation between adiabatic compressibility and speed of sound is given in equation ( 6.4 ). Equations ( 6.3 ) is used to relate the adiabatic to the isothermal compressibility using heat capacity: κ 𝑇 𝐶𝑃 = κ𝑆 𝐶𝑉 137 ( 6.3 ) where κT is the isothermal compressibility, CP is the constant pressure heat capacity, and CV is the constant volume heat capacity. The relationship between CP and CV is: 𝐶𝑃 − 𝐶𝑉 = 𝑇 ( 𝜕𝑃 𝜕𝑉 ) ( ) 𝜕𝑇 𝑉 𝜕𝑇 𝑃 ( 6.34 ) Thus, equations ( 6.3 ) and ( 6.4 ) can be rearranged: 𝐶𝑃 1 𝑣2 = √ ( ) 𝐶𝑉 𝜌𝜅𝑇 ( 6.35 ) The speed of sound propagation through a fluid depends highly on the density of the fluid. The density of the fluid is also dependent on the temperature and pressure of the system. The steep slope of the isotherm on the P-V diagram in the liquid region relates to the isothermal compressibility, which has a significant impact on the speed of sound (equation ( 6.35 )). Values for CP for pure liquids are readily available at atmospheric pressure in the literature. Values for both CP and CV as an ideal gas can also be found. To calculate values of CP and CV at elevated pressures, however, an adjustment needs to be applied to the heat capacity values. There are three options for this step. If there are experimental values for CP available at the correct pressure, these values can be used and CV can be found using equation ( 6.34 ). For ranges where there is no literature data available, the route can use CV instead. In the case of the two models that have been evaluated, one path is easier than the other for each model. 𝜌 𝜕𝐶𝑉 𝜕 2 𝑃 𝜕𝜌 𝐶𝑉 − 𝐶𝑉𝐼𝐺 = ( ) = ∫ 𝑇 ( 2) 2 𝜕𝑉 𝑇 𝜕𝑇 𝑉 𝜌 0 2 𝑟𝑒𝑠 𝜕 𝐴 𝐶𝑉 − 𝐶𝑉𝐼𝐺 = −𝑇 ( ) 𝜕𝑇 2 ( 6.36 ) ( 6.37 ) where CVIG is the ideal gas constant volume heat capacity, and Ares is the residual Helmholtz energy as shown in equation ( 6.16 ). All the required derivatives are presented in Appendix H. 138 6.2.4 Extension to mixtures The mixing rules from Maghari and Hamzehloo [120] were used for the prediction of the speed of sound of mixtures of hydrocarbons. The segment terms and the chain terms are treated differently. The different treatment allows for mixtures of molecules of different size and chain length. For the hard convex body and dispersion terms, the mixing rules are: 𝑚 = ∑ 𝑥𝑖 𝑚𝑖 ( 6.38 ) 𝑖 𝛼 = ∑ 𝑥𝑖 𝛼𝑖 𝑖 𝜂= 𝜌 ∑ 𝑥𝑖 𝑚𝑖 𝑣𝑖𝑂𝑂 𝑖 ( 6.39 ) 3 −3𝑢𝑖0 [1 − 0.12 ( )] 𝑘𝑇 1/2 𝑢𝑖𝑗 = (𝑢𝑖 𝑢𝑗 ) ∑𝑖 ∑𝑗 𝑥𝑖 𝑥𝑗 𝑚𝑖 𝑚𝑗 𝑢𝑖𝑗 𝑢= ∑𝑖 ∑𝑗 𝑥𝑖 𝑥𝑗 𝑚𝑖 𝑚𝑗 ( 6.40 ) ( 6.41 ) ( 6.42 ) For the hard convex body chain term, the mixing rule and the radial distribution function are given by: 𝐴𝑐ℎ𝑎𝑖𝑛,ℎ𝑐𝑏 = ∑ 𝑥𝑖 (1 − 𝑚𝑖 ) ln[𝑔𝑖𝑖ℎ𝑐𝑏 (𝑑𝑖𝑖 )] 𝑅𝑇 ( 6.43 ) 1 3𝛼𝑖 𝜂(1 + 𝛼𝑖 ) 4𝛼𝑖2 𝜂2 𝑔𝑖𝑖ℎ𝑐𝑏 (𝑑𝑖𝑖 ) = + + 1 − 𝜂 (1 + 3𝛼𝑖 )(1 − 𝜂)2 (1 + 3𝛼𝑖 )(1 − 𝜂)3 ( 6.44 ) 𝑖 6.3 Results and Discussion 6.3.1 Model selection Both the ESD equation and the SAFT-BACK equation were used to predict the density for n-butane and n-decane at various temperatures and pressures. At room temperature, the ESD equation does fairly well at predicting the density (Figure 6.1 and Figure 6.2). However, as pressure and temperature increase, the ESD fails to predict the density accurately. The slope of 139 the isotherm at room temperature is to low, which leads to incorrect values for the derivative properties, and the errors will propagate and lead to an underprediction of the speed of sound. The SAFT-BACK equation is better at predicting density, and matches both the values and the slope particularly well. For this reason, only the SAFT-BACK EOS will be used here to calculate the speed of sound. Figure 6.1. Prediction of the density of n-butane at various temperatures and pressures using the ESD EOS. Correlated values from the NIST Webbook [123] shown as points, ESD prediction shown as lines. 140 Figure 6.2. Prediction of the density of n-decane at various temperatures and pressures using the ESD EOS. Correlated values from the NIST Webbook [123] shown as points, ESD prediction shown as lines. Figure 6.3. Prediction of the density of n-butane at various temperatures and pressures using the SAFT-BACK EOS. Correlated values from the NIST Webbook [123] shown as points, SAFT-BACK prediction shown as lines. 141 Figure 6.4. Prediction of the density of n-decane at various temperatures and pressures using the SAFT-BACK EOS. Correlated values from the NIST Webbook [123] shown as points, SAFT-BACK prediction shown as lines. 6.3.2 Density prediction for pure fluids The SAFT-BACK equation was used to predict the density of 8 pure n-alkanes, ranging in size from methane to n-dodecane. The smaller alkanes are shown in Figure 6.5, and the larger alkanes in Figure 6.6. For the smaller molecules, methane through butane, the prediction of the density is highly accurate (Figure 6.5), with the exception of ethane (Figure 6.5 (b)). The errors were usually smaller than the size of the markers, but ethane is slightly more inaccurate. In particular, the slope of the density at the lower temperatures for these molecules are very accurate (circle markers), with the exception of ethane. However, for the larger molecules (decane and dodecane), the prediction exhibits large errors (Figure 6.6 (c) and (d)). Both the slope and values predicted are incorrect at the lower temperatures for the two largest molecules, and for the case of dodecane, the predictions are incorrect across all temperatures. The fluid is predicted as being too incompressible for the liquid phase. 142 Figure 6.5. Density and pressure for various small alkanes predicted using the SAFT-BACK equation at multiple temperatures. GERG correlation values shown as points [123] and predictions shown as lines. Individual molecules are (a) methane, (b) ethane, (c) propane, and (d) butane. 143 Figure 6.6. Density and pressure for various larger alkanes predicted using the SAFT-BACK equation at multiple temperatures. GERG correlation values shown as points [123] and predictions shown as lines. Individual molecules are (a) hexane, (b) octane, (c) decane, and (d) dodecane. 6.3.3 Speed of sound prediction for pure fluids The speed of sound for the same molecules were also predicted. For the small molecules, the predictions are accurate for all but ethane (Figure 6.7). Ethane has larger errors at the lower temperatures, due to a combination of the inaccurate density predictions. But again, as the molecules get larger, the error significantly increases (Figure 6.8). For the density, both hexane 144 and octane were predicted well, but for speed of sound, only the hexane is predicted well for only the higher temperatures. The results reported by Maghari [45] appear to have a much higher accuracy than those reported here, but percent absolute deviations were not reported in that paper. However, we were unable to reproduce Maghari’s results given the coefficients and equations stated in his paper. Because of this discrepancy, we carefully evaluated our derivatives term-by-term against numerical results and found excellent agreement to five significant figures. Our results were close to that of Chen [118] for density (1.98% AAD compared to Chen’s 2.05% AAD for decane), but Chen does not evaluate speed of sound. Literature data for dodecane was not available at the time of the Chen publishing, and so the dodecane coefficients were not optimized with data. In addition, Chen only evaluates temperatures down to 447 K, and the more significant speed of sound errors are found at the lower temperatures. 145 Figure 6.7. Speed of sound for various small alkanes predicted using the SAFT-BACK equation at multiple temperatures. GERG correlation values shown as points [116] and predictions shown as lines. Individual molecules are (a) methane, (b) ethane, (c) propane, and (d) butane. 146 Figure 6.8. Speed of sound for various larger alkanes predicted using the SAFT-BACK equation at multiple temperatures. GERG correlation values shown as points [116] and predictions shown as lines. Individual molecules are (a) hexane, (b) octane, (c) decane, and (d) dodecane. One of the potential sources of the inaccuracies in the speed of sound is the prediction of CV and CP. The calculation process uses the ideal gas CV, which is corrected using the second derivative of Helmholtz energy with respect to temperature (equation ( 6.37 )). The predicted CV is then converted to CP using the derivatives of volume and pressure with respect to temperature (equation ( 6.34 )). In the case of methane and butane, the prediction of CV becomes much more 147 inaccurate as temperature decreases and pressure increases (Figure 6.9 (a) and Figure 6.10 (a)). The derivatives used to convert CV to CP fortuitously cancel the error in CV in the case of methane and butane. However, as the molecules become bigger, the error in the second derivative becomes much larger, especially at lower temperatures (Figure 6.11 and Figure 6.12 (a)). In the cases of octane and dodecane, the low temperature CV has the completely wrong slope as well. Because the error is so large, the cancellation of errors is no longer able to compensate, and the prediction of speed of sound deteriorates. Figure 6.9. Predicted heat capacities for methane at multiple temperatures. GERG correlation values from [116]. 148 Figure 6.10. Predicted heat capacities for butane at multiple temperatures. GERG correlation values from [116]. Figure 6.11. Predicted heat capacities for octane at multiple temperatures. GERG correlation values from [116]. 149 Figure 6.12. Predicted heat capacities for dodecane at multiple temperatures. GERG correlation values from [116]. In an attempt to find the source of these errors, the roots of the SAFT-BACK equation were examined. Normally, the equation is used to find the correct density for a given temperature and pressure, but it can also be used to calculate the pressure at various theoretical densities (Figure 6.13 and Figure 6.14). Plotting an isotherm provides visualization of the theoretical density roots at a given pressure. In the case of a cubic equation, such as the ESD equation, there can be one or three unique real solutions or “roots” at a given pressure. For the SAFT-BACK equation, there can be up to five real roots at a given pressure. The behavior of the SAFT EOS roots has been discussed in detail in the literature [124, 125]. 150 Figure 6.13. Pressure predicted by SAFT-BACK at multiple temperatures for dodecane at various densities. The segments at high molar volume are connected through a minimum that is off scale. Red dotted line represents the close-packed density. The liquid densities of interest for dodecane shown in Figure 6.13 are to the left of the red dotted line, which represents the close-packed density. The hump that is seen around 0.0044 m3/mol represents the saturated liquid (3rd and 4th) roots, and the height of it decreases with decreasing temperature (Figure 6.13). For this work, we seek to model liquids up to 50 MPa. When the height of this hump decreases below the pressure of interest, these roots disappear as solutions, and the only root remaining is the nonphysical 5th root. This helps to explain the why the predictions become spurious as the temperature decreases. A similar trend can be seen for multiple components at room temperature (Figure 6.14). The liquid densities of interest for propane shown in Figure 6.14 are to the left of the maximum pressure hump at 700 MPa. This pressure high enough for propane that experimental conditions for propane can be represented by the SAFT-BACK as shown in Figure 6.14. In order to represent compressed liquids, this pressure maximum must be above the experimental pressures. 151 However, the 5th root is above the close-packed density for propane (0.024 m3/mol), and as such should not be considered as a realistic root. In the case of all molecules shown in Figure 6.14, this 5th root is above the close-packed density for each molecule. Figure 6.14. Pressure predicted by SAFT-BACK at 298.15 K for multiple n-alkanes at various densities. Curves for all molecules except propane have been truncated after the 4th root, to improve visibility. The two segments at high molar volume for propane are connected through a minimum that is off scale. A distinct trend in the height of the stable liquid region “hump” at 298.15 K can been seen in Figure 6.14. An enlarged image is shown in Figure 6.15. As the molecule size increases, the height of the hump decreases strongly. As the size of the molecule increases, the maximum pressure that be represented at 298.15 K decreases. For hexadecane, a molecule representative of a jet/diesel fuel component, the liquid molar volume below 10 MPa is in the range of 0.0016 m3/mol, but above about 10 MPa, the only molar volume root jumps to the range of 0.0058 m3/mol (the nonphysical 5th root) and is a nonphysical solution. 152 Figure 6.15. Zoomed in version of Figure 6.14 to show detail in the pressure range of interest for the three largest alkanes. All roots are shown in this version of the figure, and segments are connected through a minimum which is off scale. 6.3.4 Speed of sound prediction for mixtures The speed of sound for binary mixtures was predicted using the compounds for which the SAFT-BACK worked moderately well. Mixtures included ethane/butane and propane/hexane, and were evaluated at multiple compositions and temperatures. It was found that the mixing rules work moderately well. As long as the individual compounds were predicted well at the selected temperatures, the mixtures were predicted accurately (Figure 6.16 and Figure 6.17). These results are similar to those shown by Maghari [120]. Maghari’s results are again more accurate than we were able to reproduce at the higher molecular weights, but we have very similar results at the low molecular weights. 153 Figure 6.16. Speed of sound predicted for mixtures of ethane and butane at various temperatures, compared to GERG correlation values from [116]. (a) 200 K, (b) 300 K, and (c) 400 K. 154 Figure 6.17. Speed of sound predicted for mixtures of propane and hexane at various temperatures, compared to GERG correlation values from [116]. (a) 300 K, (b) 400 K, and (c) 500 K. 6.4 Conclusions Density and speed of sound were predicted for multiple n-alkanes and some binary mixtures. Both the ESD and SAFT-BACK equations were evaluated for density, and the SAFTBACK equation was found to work better for all molecules selected. The SAFT-BACK equation was then used to predict the density, speed of sound, and heat capacities for eight n-alkanes, ranging in size from methane to n-dodecane. Density was predicted well for molecules smaller than octane, and speed of sound for molecules smaller than hexane. The slope of CV is predicted incorrectly, which causes the fluid to be predicted as too compressible. This leads to the speed of 155 sound being overpredicted for the larger molecules. The mixtures are well predicted in the temperature and pressure ranges where the pure fluids are well predicted. The SAFT-BACK equation works well for small molecules and high temperatures, but fails for larger molecules near room temperature because the only real volume root at high pressure is the nonphysical root that is smaller than the close-packed volume. 156 APPENDICES 157 APPENDIX A Cloud point and pour point data for all mixtures 158 This appendix gives all the raw data and adjusted data for the cloud point and pour point temperature measurements. Table A.1. Cloud point temperatures for mixtures of USD and various additives (Part 1). Additive USD only Dihexyl ketone Dibutyl succinate Cloud Point Mass % Temperature Additive (˚C) 0.0 -27.1 0.0 -26.8 9.7 -4.7 9.7 -5.2 9.7 -8.1 14.6 0.3 14.6 0 14.6 0.7 24.4 9.6 24.4 8.9 24.4 9.9 16.5 -14.4 16.5 -14.2 16.5 -14.2 27.1 -13.3 27.1 -13.4 27.1 -13.4 27.1 -13.9 27.1 -13.8 27.1 -13.9 52.8 -11 52.8 -11.3 52.8 -11.8 52.8 -11.8 52.8 -12.2 77.0 -11 77.0 -10 77.0 -11.5 86.4 -11.3 86.4 -11.6 86.4 -12.1 91.0 -17.7 91.0 -15.9 91.0 -16.2 91.0 -15.9 159 Additive 4-Heptanone Butyl nonanoate Dihexyl ether Cloud Point Mass % Temperature Additive (˚C) 24.2 -16.7 24.2 -16.6 24.2 -16.7 48.9 -17 48.9 -17.1 48.9 -16.9 74.1 -19.3 74.1 -19 74.1 -19.1 84.4 -22 84.4 -21.5 84.4 -21.9 5.0 -15.5 5.0 -15.4 5.0 -15 15.0 -16 15.0 -15.8 15.0 -15.8 24.9 -16.4 24.9 -16 24.9 -16.1 49.9 -17.6 49.9 -18.3 49.9 -18.3 74.9 -23.2 74.9 -22.8 74.9 -23.1 85.0 -27.8 85.0 -25.5 85.0 -26 23.6 -17.6 23.6 -17.7 23.6 -17.6 73.6 -27.4 73.6 -27.2 73.6 -27.3 Table A.2. Cloud point temperatures for mixtures of USD and various additives (Part 2). Additive Isobutyl nonanoate Diisobutyl succinate Butyl butyrate Mass % Additive 25.4 25.4 25.4 50.5 50.5 50.5 75.4 75.4 75.4 27.1 27.1 27.1 52.8 52.8 52.8 77.0 77.0 77.0 15.3 15.3 15.3 25.4 25.4 25.4 50.5 50.5 50.5 75.4 75.4 75.4 85.3 85.3 85.3 Cloud Point Temperature (˚C) -16.1 -16.4 -16.2 -18 -18 -18.2 -21.5 -22.1 -22.4 -13.5 -13.6 -13.4 -11.8 -11.2 -11.9 -10.6 -11.5 -10.5 -15.9 -16 -15.9 -16.7 -16.4 -16.6 -18 -17.9 -18 -21.1 -21 -20.9 -23.9 -24.2 -24.2 160 Additive Hexanes Dibutyl ether Ethyl hexyl nonanoate Mass % Additive 20.6 20.6 20.6 43.7 43.7 43.7 70.0 70.0 70.0 81.5 81.5 23.1 23.1 23.1 73.0 73.0 73.0 24.9 24.9 24.9 74.9 74.9 74.9 Cloud Point Temperature (˚C) -21.6 -21.6 -21.6 -26.7 -26.4 -26.8 -35.1 -35 -34.6 -41.1 -39.5 -18.3 -18.3 -18.2 -29.1 -29.2 -29.1 -15.8 -15.7 -15.9 -22.7 -23.1 -22.7 Table A.3. Cloud point temperatures for mixtures of HAD and various additives (Part 1). Additive HAD only Dihexyl ketone Dibutyl succinate Mass % Additive 0 0 9.7 9.7 9.7 24.4 24.4 24.4 16.5 16.5 16.5 27.1 27.1 27.1 52.8 52.8 52.8 77.0 77.0 77.0 77.0 77.0 86.4 86.4 86.4 91.0 91.0 91.0 95.5 95.5 95.5 Cloud Point Temperature (˚C) -14.3 -15.2 -7 -7.6 -9 6.8 7.2 6.1 -25.9 -26.1 -26.1 -26 -25.9 -25.6 -24.1 -24.9 -24.2 -23.6 -23.6 -23.7 -23.4 -23.3 -25.6 -25.5 -25.6 -28 -27.9 -28 -34.6 -33.5 -34.2 161 Additive 4-Heptanone Butyl nonanoate Dihexyl ether Mass % Additive 48.9 48.9 48.9 74.1 74.1 79.3 79.3 79.3 5.0 5.0 5.0 15.0 15.0 15.0 24.9 24.9 24.9 49.9 49.9 49.9 74.9 74.9 74.9 85.0 85.0 85.0 23.6 23.6 23.6 73.6 73.6 73.6 Cloud Point Temperature (˚C) -29.8 -29.8 -29.7 -33 -33 -34 -34.2 -34.1 -27.6 -27.6 -27.7 -28.5 -28 -28.2 -28.5 -28.4 -28.5 -30.5 -31.1 -30.4 -35.7 -35.6 -35.8 -38.5 -38.5 -38 -30 -30.1 -29.9 -38.9 -38.5 -38.7 Table A.4. Cloud point temperatures for mixtures of HAD and various additives (Part 2). Cloud Point Mass % Temperature Additive (˚C) 25.4 -28.8 25.4 -28.6 25.4 -28.6 50.5 -30.5 Isobutyl nonanoate 50.5 -30.4 50.5 -30.7 75.4 -35.4 75.4 -35.9 75.4 -35.7 27.1 -25.6 27.1 -25.6 27.1 -25.7 52.8 -23.7 Diisobutyl succinate 52.8 -24.2 52.8 -23.8 77.0 -23 77.0 -22.9 77.0 -23.3 15.3 -28.6 15.3 -28.6 15.3 -28.5 25.4 -28.8 25.4 -28.8 25.4 -28.9 50.5 -30.9 Butyl butyrate 50.5 -31 50.5 -31 75.4 -35.5 75.4 -35.4 75.4 -34.9 85.3 -38.2 85.3 -37.8 85.3 -37 Additive 162 Additive Hexanes Dibutyl ether Ethyl hexyl nonanoate Butyl ethyl hexanoate Mass % Additive 20.6 20.6 20.6 43.7 43.7 43.7 70.0 70.0 70.0 81.5 81.5 81.5 23.1 23.1 23.1 73.0 73.0 73.0 24.9 24.9 24.9 49.9 49.9 49.9 74.9 74.9 74.9 24.9 24.9 24.9 49.9 49.9 49.9 74.9 74.9 74.9 Cloud Point Temperature (˚C) -34.1 -34 -33.7 -39.8 -40.1 -39.5 -48.9 -48.7 -48.4 -54.2 -53.8 -54.3 -31.8 -31.5 -31.3 -42.9 -42.3 -42.8 -28.2 -28.4 -28.8 -30.5 -31.1 -30.9 -36.2 -36.5 -35.9 -28.9 -29.2 -29.3 -31 -31.3 -31.5 -35.9 -36 -35.8 Table A.5. Cloud point (CP), adjusted cloud point (CP*), and cold filter plugging point (CFPP) temperatures for mixtures of canola FAME in the HAD+ fuel. The cloud point data was collected using ASTM D2500, and CFFP data collected by ASTM D6371. The cloud point data is adjusted for comparison with ASTM D7683 using Cloud point (D2500) = Cloud point (D7683) + 1.68. Mass% Canola FAME CP (˚C) CP* (˚C) CFPP (˚C) 0.0 0.0 10.3 10.3 15.4 20.5 20.5 30.7 30.7 40.8 40.8 60.8 80.5 100.0 -24 -24 -20 -20 -20 -19 -19 -16 -16 -15 -15 -9 -7 -6 -25.7 -25.7 -21.7 -21.7 -21.7 -20.7 -20.7 -17.7 -17.7 -16.7 -16.7 -10.7 -8.7 -7.7 -38 -36 -35 -34 -33 -31 -31 -28 -28 -21 -21 -19 -14 -10 Table A.6. Cloud point (CP) and adjusted cloud point (CP*) temperatures for mixtures of dibutyl succinate in the HAD* fuel. The cloud point data was collected using ASTM D2500. The cloud point data is adjusted for comparison with ASTM D7683 using Cloud point (D2500) = Cloud point (D7683) + 1.68. Mass % Dibutyl Succinate CP (˚C) CP* (˚C) 0 0.0 11.3 22.3 22.3 33.0 43.3 53.4 63.2 72.8 82.1 91.2 91.2 100.0 -26 -25 -25 -24 -24 -24 -23 -22 -21 -21 -21 -24 -22 -29 -27.7 -26.7 -26.7 -25.7 -25.7 -25.7 -24.7 -23.7 -22.7 -22.7 -22.7 -25.7 -23.7 -29 163 Table A.7. Cloud point (CP), adjusted cloud point (CP*), and cold filter plugging point (CFPP) temperatures for mixtures of dibutyl succinate in the HAD+ fuel. The cloud point data was collected using ASTM D2500. The cloud point data is adjusted for comparison with ASTM D7683 using Cloud point (D2500) = Cloud point (D7683) + 1.68. Mass % Dibutyl Succinate CP (˚C) CP* (˚C) CFPP (˚C) 0.0 0.0 11.3 22.3 33.0 43.3 53.4 63.2 72.8 82.1 91.2 100.0 -24 -24 -23 -23 -22 -21 -20 -20 -19 -20 -23 -29 -25.7 -25.7 -24.7 -24.7 -23.7 -22.7 -21.7 -21.7 -20.7 -21.7 -24.7 -29 -38 -36 -41 -41 -41 -42 -43 -42 -42 -41 -41 -41 Table A.8. Cold filter plugging point (CFPP) temperatures for mixtures of butyl nonanoate in the HAD fuel. Mass % Butyl Nonanoate CFPP (˚C) 0.0 0.0 10.0 10.0 29.9 29.9 49.9 49.9 69.9 69.9 89.9 89.9 -38 -36 -32 -32 -33 -33 -35 -35 -47 -42 -50 -49 164 APPENDIX B Tabulated data for the prediction of physical properties 165 This appendix tabulates the values of all the constants given in Chapter 3, as well as the predicted data for all compounds. Table B.1. Constants for the regression of Cp/R for the lower temperature range. A1 A2 A3 A4 A5 Dibutyl Succinate 32.619 6.004 -0.187 -6.91E-04 1.40E-04 Methyl Oleate -5.455 0.245 -1.83E-04 7.46E-08 -1.31E-11 Methyl Stearate -2.430 0.239 -1.75E-04 7.36E-08 -1.43E-11 Cetane -10.191 0.234 -1.98E-04 9.58E-08 -1.97E-11 Heptamethylnonane -12.176 0.241 -1.92E-04 8.04E-08 -1.40E-11 TGME -9.185 0.198 -2.07E-04 1.26E-07 -3.25E-11 Table B.2. Constants for the regression of Cp/R for the high temperature range. A1 A2 A3 A4 A5 Dibutyl Succinate 36.700 7.31E-02 -2.91E-05 5.28E-09 -3.59E-13 Methyl Oleate 47.343 0.115 -4.58E-05 8.30E-09 -5.64E-13 Methyl Stearate 49.352 0.121 -4.81E-05 8.72E-09 -5.93E-13 Cetane 36.630 0.109 -4.39E-05 8.01E-09 -5.48E-13 Heptamethylnonane 40.981 0.103 -4.07E-05 7.35E-09 -4.98E-13 TGME 33.656 6.99E-02 -2.77E-05 5.01E-09 -3.40E-13 Table B.3. Constants for the regression of H/RT for the low temperature range. A1 A2 A3 A4 A5 A6 Dibutyl Succinate 32.619 6.004 -0.187 -6.91E-04 1.40E-04 -6.31E+07 Methyl Oleate -5.455 0.245 -1.83E-04 7.46E-08 -1.31E-11 -7.78E+03 Methyl Stearate -2.430 2.39E-01 -1.75E-04 7.36E-08 -1.43E-11 -8.51E+03 Heptamethylnonane -10.191 -12.176 2.34E-01 0.241 -1.98E-04 -1.92E-04 9.58E-08 8.04E-08 -1.97E-11 -1.40E-11 -5.79E+03 -5.55E+03 Cetane TGME -9.185 0.198 -2.07E-04 1.26E-07 -3.25E-11 -4.45E+03 Table B.4. Constants for the regression of H/RT for the high temperature range. A1 A2 A3 A4 A5 A6 Dibutyl Methyl Succinate Oleate 36.700 47.343 0.073 0.115 -2.91E-05 -4.58E-05 5.28E-09 8.30E-09 -3.59E-13 -5.64E-13 -1.39E+04 -1.88E+04 Methyl Stearate 49.352 0.121 -4.81E-05 8.72E-09 -5.93E-13 -1.97E+04 166 Heptamethylnonane 36.630 40.981 0.109 0.103 -4.39E-05 -4.07E-05 8.01E-09 7.35E-09 -5.48E-13 -4.98E-13 -1.54E+04 -1.64E+04 Cetane TGME 33.656 0.070 -2.77E-05 5.01E-09 -3.40E-13 -1.29E+04 Table B.5. Constants for the regression of S/T for the low temperature range. A1 A2 A3 A4 A5 A7 Dibutyl Succinate 32.619 6.004 -0.187 -6.91E-04 1.40E-04 -2.64E+05 Methyl Oleate -5.455 0.245 -1.83E-04 7.46E-08 -1.31E-11 -34.423 Methyl Stearate -2.430 0.239 -1.75E-04 7.36E-08 -1.43E-11 -50.396 Cetane -10.191 2.34E-01 -1.98E-04 9.58E-08 -1.97E-11 -3.703 HeptamethylTGME nonane -12.176 -9.185 0.241 0.198 -1.92E-04 -2.07E-04 8.04E-08 1.26E-07 -1.40E-11 -3.25E-11 5.274 1.538 Table B.6. Constants for the regression of S/T for the high temperature range. A1 A2 A3 A4 A5 A7 Dibutyl Methyl Methyl Succinate Oleate Stearate 36.700 47.343 49.352 0.073 0.115 0.121 -2.91E-05 -4.58E-05 -4.81E-05 5.28E-09 8.30E-09 8.72E-09 -3.59E-13 -5.64E-13 -5.93E-13 -2.30E+02 -3.02E+02 -3.15E+02 167 Cetane 36.630 0.109 -4.39E-05 8.01E-09 -5.48E-13 -2.39E+02 HeptamethylTGME nonane 40.981 33.656 0.103 0.070 -4.07E-05 -2.77E-05 7.35E-09 5.01E-09 -4.98E-13 -3.40E-13 -2.62E+02 -2.11E+02 Table B.7. Predicted data for the viscosity of various compounds extrapolated to the critical point (Part 1). Dibutyl Succinate Methyl Oleate Methyl Stearate Temperature Viscosity Temperature Viscosity Temperature Viscosity (K) (Pa*s) (K) (Pa*s) (K) (Pa*s) 250 2.34E-02 293 1.04E-02 310 6.73E-03 260 1.64E-02 300 8.53E-03 320 5.19E-03 270 1.18E-02 310 6.53E-03 330 4.07E-03 280 8.73E-03 320 5.08E-03 340 3.23E-03 290 6.57E-03 330 4.02E-03 350 2.60E-03 300 5.04E-03 340 3.22E-03 360 2.12E-03 310 3.94E-03 350 2.61E-03 370 1.75E-03 320 3.12E-03 360 2.14E-03 380 1.46E-03 330 2.51E-03 370 1.78E-03 390 1.22E-03 340 2.04E-03 380 1.49E-03 400 1.04E-03 350 1.68E-03 390 1.26E-03 410 8.86E-04 360 1.40E-03 400 1.07E-03 420 7.63E-04 370 1.18E-03 410 9.24E-04 430 6.62E-04 380 1.00E-03 420 7.99E-04 440 5.77E-04 390 8.58E-04 430 6.97E-04 450 5.07E-04 400 7.40E-04 440 6.11E-04 460 4.48E-04 410 6.43E-04 450 5.39E-04 470 3.97E-04 420 5.63E-04 460 4.78E-04 480 3.55E-04 430 4.95E-04 470 4.26E-04 490 3.18E-04 440 4.39E-04 480 3.81E-04 500 2.86E-04 450 3.90E-04 490 3.43E-04 510 2.59E-04 460 3.49E-04 500 3.10E-04 520 2.35E-04 470 3.14E-04 510 2.81E-04 530 2.14E-04 480 2.84E-04 520 2.56E-04 540 1.95E-04 490 2.57E-04 530 2.34E-04 550 1.79E-04 500 2.34E-04 540 2.15E-04 560 1.65E-04 510 2.14E-04 550 1.97E-04 570 1.52E-04 520 1.96E-04 560 1.82E-04 580 1.41E-04 530 1.80E-04 570 1.68E-04 590 1.30E-04 540 1.66E-04 580 1.56E-04 600 1.21E-04 550 1.54E-04 590 1.45E-04 610 1.13E-04 560 1.43E-04 600 1.35E-04 620 1.06E-04 570 1.33E-04 610 1.27E-04 630 9.88E-05 580 1.24E-04 620 1.18E-04 640 9.27E-05 590 1.16E-04 630 1.11E-04 650 8.71E-05 600 1.09E-04 640 1.05E-04 660 8.20E-05 610 1.02E-04 650 9.84E-05 670 7.74E-05 620 9.60E-05 660 9.29E-05 680 7.31E-05 630 9.05E-05 670 8.78E-05 690 6.92E-05 640 8.54E-05 680 8.32E-05 700 6.56E-05 650 8.08E-05 690 7.89E-05 710 6.23E-05 660 7.66E-05 700 7.49E-05 720 5.93E-05 670 7.27E-05 710 7.13E-05 730 5.64E-05 680 6.91E-05 720 6.79E-05 740 5.38E-05 690 6.58E-05 730 6.47E-05 750 5.14E-05 700 6.28E-05 740 6.18E-05 760 4.91E-05 710 5.99E-05 750 5.91E-05 770 4.70E-05 760 5.66E-05 780 4.50E-05 764 5.56E-05 781 4.48E-05 168 Table B.8. Predicted data for the viscosity of various compounds extrapolated to the critical point (Part 2). Cetane Heptamethylnonane TGME Temperature Viscosity Temperature Viscosity Temperature Viscosity (K) (Pa*s) (K) (Pa*s) (K) (Pa*s) 291 3.47E-03 163 2.44E+01 231 5.04E-02 300 2.86E-03 170 1.12E+01 240 3.24E-02 310 2.34E-03 180 4.06E+00 250 2.06E-02 320 1.94E-03 190 1.64E+00 260 1.35E-02 330 1.63E-03 200 7.28E-01 270 9.16E-03 340 1.38E-03 210 3.48E-01 280 6.38E-03 350 1.18E-03 220 1.78E-01 290 4.56E-03 360 1.02E-03 230 9.68E-02 300 3.33E-03 370 8.87E-04 240 5.52E-02 310 2.49E-03 380 7.77E-04 260 2.05E-02 320 1.89E-03 390 6.85E-04 270 1.32E-02 330 1.46E-03 400 6.08E-04 280 8.76E-03 340 1.14E-03 410 5.43E-04 290 5.98E-03 350 9.10E-04 420 4.88E-04 300 4.19E-03 360 7.33E-04 430 4.40E-04 310 3.01E-03 370 5.97E-04 440 3.99E-04 320 2.20E-03 380 4.92E-04 450 3.63E-04 330 1.64E-03 390 4.09E-04 460 3.32E-04 340 1.25E-03 400 3.43E-04 470 3.04E-04 360 7.52E-04 410 2.91E-04 480 2.80E-04 370 5.96E-04 420 2.48E-04 490 2.59E-04 380 4.78E-04 430 2.13E-04 500 2.40E-04 390 3.88E-04 440 1.85E-04 510 2.23E-04 410 2.64E-04 450 1.61E-04 520 2.08E-04 420 2.20E-04 460 1.41E-04 530 1.95E-04 430 1.86E-04 470 1.24E-04 540 1.82E-04 440 1.58E-04 480 1.10E-04 550 1.71E-04 460 1.16E-04 490 9.81E-05 560 1.61E-04 470 1.01E-04 500 8.78E-05 570 1.52E-04 480 8.77E-05 510 7.89E-05 580 1.44E-04 490 7.69E-05 520 7.12E-05 590 1.36E-04 500 6.78E-05 530 6.45E-05 600 1.29E-04 510 6.00E-05 540 5.86E-05 610 1.23E-04 520 5.34E-05 550 5.34E-05 620 1.17E-04 530 4.77E-05 560 4.89E-05 630 1.11E-04 540 4.28E-05 570 4.49E-05 640 1.06E-04 560 3.49E-05 580 4.14E-05 650 1.02E-04 570 3.17E-05 590 3.82E-05 660 9.75E-05 580 2.89E-05 600 3.54E-05 670 9.35E-05 590 2.64E-05 610 3.28E-05 680 8.97E-05 600 2.42E-05 620 3.05E-05 690 8.63E-05 610 2.22E-05 630 2.85E-05 700 8.30E-05 620 2.05E-05 640 2.66E-05 710 8.00E-05 630 1.89E-05 650 2.49E-05 720 7.71E-05 640 1.75E-05 657 2.38E-05 723 7.63E-05 660 1.51E-05 670 1.41E-05 680 1.32E-05 690 1.23E-05 692 1.22E-05 169 Table B.9. Predicted data for the heat capacity of various compounds extrapolated to the critical point (Part 1). Dibutyl Succinate Methyl Oleate Methyl Stearate Temperature Heat Capacity Temperature Heat Capacity Temperature Heat Capacity (K) (J/kgK) (K) (J/kgK) (K) (J/kgK) 250 1.87E+03 293 2.09E+03 310 2.10E+03 260 1.89E+03 300 2.11E+03 320 2.13E+03 270 1.92E+03 310 2.14E+03 330 2.16E+03 280 1.94E+03 320 2.17E+03 340 2.19E+03 290 1.97E+03 330 2.20E+03 350 2.22E+03 300 1.99E+03 340 2.23E+03 360 2.25E+03 310 2.02E+03 350 2.25E+03 370 2.28E+03 320 2.05E+03 360 2.28E+03 380 2.31E+03 330 2.08E+03 370 2.31E+03 390 2.34E+03 340 2.10E+03 380 2.34E+03 400 2.37E+03 350 2.13E+03 390 2.37E+03 410 2.40E+03 360 2.16E+03 400 2.40E+03 420 2.43E+03 370 2.19E+03 410 2.43E+03 430 2.46E+03 380 2.22E+03 420 2.46E+03 440 2.49E+03 390 2.24E+03 430 2.48E+03 450 2.52E+03 400 2.27E+03 440 2.51E+03 460 2.55E+03 410 2.30E+03 450 2.54E+03 470 2.58E+03 420 2.33E+03 460 2.57E+03 480 2.61E+03 430 2.36E+03 470 2.60E+03 490 2.63E+03 440 2.39E+03 480 2.62E+03 500 2.66E+03 450 2.41E+03 490 2.65E+03 510 2.69E+03 460 2.44E+03 500 2.68E+03 520 2.72E+03 470 2.47E+03 510 2.71E+03 530 2.75E+03 480 2.50E+03 520 2.74E+03 540 2.78E+03 490 2.53E+03 530 2.76E+03 550 2.80E+03 500 2.56E+03 540 2.79E+03 560 2.83E+03 510 2.58E+03 550 2.82E+03 570 2.86E+03 520 2.61E+03 560 2.84E+03 580 2.89E+03 530 2.64E+03 570 2.87E+03 590 2.92E+03 540 2.67E+03 580 2.90E+03 600 2.94E+03 550 2.70E+03 590 2.93E+03 610 2.97E+03 560 2.73E+03 600 2.96E+03 620 3.00E+03 570 2.76E+03 610 2.99E+03 630 3.03E+03 580 2.79E+03 620 3.02E+03 640 3.06E+03 590 2.83E+03 630 3.05E+03 650 3.09E+03 600 2.86E+03 640 3.08E+03 660 3.13E+03 610 2.90E+03 650 3.11E+03 670 3.16E+03 620 2.94E+03 660 3.15E+03 680 3.20E+03 630 2.98E+03 670 3.19E+03 690 3.24E+03 640 3.04E+03 680 3.23E+03 700 3.28E+03 650 3.10E+03 690 3.28E+03 710 3.33E+03 660 3.17E+03 700 3.34E+03 720 3.39E+03 670 3.28E+03 710 3.42E+03 730 3.47E+03 680 3.43E+03 720 3.52E+03 740 3.58E+03 690 3.68E+03 730 3.67E+03 750 3.74E+03 700 4.25E+03 740 3.93E+03 760 4.04E+03 710 6.85E+03 750 4.54E+03 770 4.86E+03 760 8.17E+03 775 6.25E+03 170 Table B.10. Predicted data for the heat capacity of various compounds extrapolated to the critical point (Part 2). Cetane Heptamethylnonane TGME Temperature Heat Capacity Temperature Heat Capacity Temperature Heat Capacity (K) (J/kgK) (K) (J/kgK) (K) (J/kgK) 291 2.16E+03 163 1.49E+03 231 2.11E+03 300 2.19E+03 170 1.52E+03 240 2.12E+03 310 2.23E+03 180 1.56E+03 250 2.14E+03 320 2.27E+03 190 1.60E+03 260 2.15E+03 330 2.31E+03 200 1.64E+03 270 2.17E+03 340 2.35E+03 210 1.68E+03 280 2.19E+03 350 2.38E+03 220 1.73E+03 290 2.21E+03 360 2.42E+03 230 1.77E+03 300 2.24E+03 370 2.46E+03 240 1.82E+03 310 2.26E+03 380 2.50E+03 260 1.91E+03 320 2.28E+03 390 2.53E+03 270 1.95E+03 330 2.30E+03 400 2.57E+03 280 2.00E+03 340 2.33E+03 410 2.61E+03 290 2.04E+03 350 2.35E+03 420 2.64E+03 300 2.09E+03 360 2.38E+03 430 2.68E+03 310 2.13E+03 370 2.40E+03 440 2.71E+03 320 2.18E+03 380 2.42E+03 450 2.75E+03 330 2.22E+03 390 2.45E+03 460 2.79E+03 340 2.26E+03 400 2.47E+03 470 2.82E+03 360 2.35E+03 410 2.50E+03 480 2.86E+03 370 2.39E+03 420 2.52E+03 490 2.89E+03 380 2.44E+03 430 2.55E+03 500 2.92E+03 390 2.48E+03 440 2.57E+03 510 2.96E+03 400 2.52E+03 450 2.60E+03 520 2.99E+03 410 2.56E+03 460 2.63E+03 530 3.03E+03 420 2.60E+03 470 2.65E+03 540 3.06E+03 430 2.64E+03 480 2.68E+03 550 3.10E+03 440 2.68E+03 490 2.71E+03 560 3.13E+03 460 2.76E+03 500 2.73E+03 570 3.17E+03 470 2.80E+03 510 2.76E+03 580 3.20E+03 480 2.84E+03 520 2.79E+03 590 3.24E+03 490 2.88E+03 530 2.82E+03 600 3.28E+03 500 2.92E+03 540 2.86E+03 610 3.32E+03 510 2.95E+03 550 2.90E+03 620 3.36E+03 520 2.99E+03 560 2.94E+03 630 3.41E+03 530 3.03E+03 570 2.98E+03 640 3.46E+03 540 3.07E+03 580 3.03E+03 650 3.52E+03 560 3.15E+03 590 3.10E+03 660 3.59E+03 570 3.19E+03 600 3.18E+03 670 3.67E+03 580 3.23E+03 610 3.29E+03 680 3.79E+03 590 3.28E+03 620 3.45E+03 690 3.96E+03 600 3.33E+03 630 3.72E+03 700 4.26E+03 610 3.38E+03 640 4.29E+03 710 4.98E+03 620 3.44E+03 650 6.47E+03 720 1.05E+04 630 3.50E+03 640 3.59E+03 660 3.86E+03 670 4.14E+03 680 4.86E+03 690 1.26E+04 171 Table B.11. Predicted gas phase heat capacities for various compounds. Units are J/kgK. Temperature Dibutyl Methyl Methyl HeptamethylCetane TGME (K) Succinate Oleate Stearate nonane 1000 3079.94 3474.89 3615.14 3991.84 4026.39 3231.88 1250 3335.16 3787.03 3940.80 4376.07 4392.90 3504.67 1500 3507.66 3998.34 4161.31 4631.37 4643.60 3690.17 1750 3627.26 4144.93 4314.26 4805.96 4818.62 3819.31 2000 3712.44 4249.38 4423.24 4929.11 4943.87 3911.57 2250 3774.78 4325.78 4502.95 5018.52 5035.74 3979.21 2500 3821.49 4383.05 4562.70 5085.17 5104.75 4029.98 2750 3857.26 4426.91 4608.45 5136.00 5157.69 4068.92 3000 3885.20 4461.16 4644.18 5175.56 5199.06 4099.36 3250 3907.40 4488.37 4672.57 5206.90 5231.97 4123.56 3500 3925.29 4510.30 4695.44 5232.12 5258.52 4143.07 3750 3939.92 4528.23 4714.15 5252.70 5280.23 4159.05 4000 3952.02 4543.06 4729.63 5269.70 5298.21 4172.26 4250 3962.14 4555.47 4742.56 5283.91 5313.23 4183.31 4500 3970.70 4565.94 4753.48 5295.88 5325.93 4192.64 4750 3977.96 4574.86 4762.79 5306.08 5336.74 4200.59 5000 3984.21 4582.51 4770.77 5314.82 5346.03 4207.43 172 UNCLASSIFIED: Distribution Statement A. Approved for public release. APPENDIX C Distillation Curve Data UNCLASSIFIED 173 UNCLASSIFIED This appendix provides supporting data for Chapter 4. The measured temperature is corrected to atmospheric pressure using the Sydney Young equation as cited in Chapter 4. Table C.1. Uncorrected and corrected distillation data for JP-8. Data taken with an atmospheric pressure of 737 mmHg. Volume Distilled (mL) Measured T (˚C) Pressure Correction (˚C) Corrected T (˚C) 0 166.5 1.213 167.7 100 182.5 1.257 183.8 200 185.9 1.267 187.2 300 188.5 1.274 189.8 400 191.2 1.281 192.5 500 194.8 1.291 196.1 600 196.5 1.296 197.8 700 199.1 1.303 200.4 800 201 1.308 202.3 900 205 1.319 206.3 1000 208 1.328 209.3 1100 211 1.336 212.3 1200 214 1.344 215.3 1300 218 1.355 219.4 1400 222 1.366 223.4 1500 227 1.380 228.4 1600 232 1.394 233.4 1700 238 1.410 239.4 1800 244 1.427 245.4 1900 255 1.457 256.5 UNCLASSIFIED 174 UNCLASSIFIED Table C.2. Uncorrected and corrected distillation data for IPK. Data taken with an atmospheric pressure of 746 mmHg. Volume Distilled (mL) Measured T (˚C) Pressure Correction (˚C) Corrected T (˚C) 0 161.8 0.7305 162.5 100 169.6 0.7436 170.3 200 170.8 0.7456 171.5 300 171.9 0.7474 172.6 400 172.8 0.7489 173.5 500 174.2 0.7513 175.0 600 175.2 0.7530 176.0 700 176.6 0.7553 177.4 800 178 0.7577 178.8 900 179.4 0.7600 180.2 1000 181 0.7627 181.8 1100 182.8 0.7657 183.6 1200 184.7 0.7689 185.5 1300 186.8 0.7725 187.6 1400 189.2 0.7765 190.0 1500 191.9 0.7810 192.7 1600 195.1 0.7864 195.9 1700 198.9 0.7928 199.7 1800 205 0.8030 205.8 1900 214 0.8182 214.8 UNCLASSIFIED 175 UNCLASSIFIED Table C.3. Uncorrected and corrected distillation data for HRJ. Data taken with an atmospheric pressure of 746 mmHg. Volume Distilled (mL) Measured T (˚C) Pressure Correction (˚C) Corrected T (˚C) 0 160.2 0.7278 160.9 100 183.6 0.7671 184.4 200 191.4 0.7802 192.2 300 197 0.7896 197.8 400 201 0.7963 201.8 500 207 0.8064 207.8 600 211 0.8131 211.8 700 215 0.8198 215.8 800 220 0.8282 220.8 900 224 0.8350 224.8 1000 229 0.8434 229.8 1100 232 0.8484 232.8 1200 237 0.8568 237.9 1300 241 0.8635 241.9 1400 245 0.8702 245.9 1500 250 0.8786 250.9 1600 255 0.8870 255.9 1700 259 0.8938 259.9 1800 264 0.9022 264.9 1900 270 0.9122 270.9 UNCLASSIFIED 176 UNCLASSIFIED Table C.4. Uncorrected and corrected distillation data for SPK2. Data taken with an atmospheric pressure of 758 mmHg. Volume Distilled (mL) Measured T (˚C) Pressure Correction (˚C) Corrected T (˚C) 0 159.7 0.0519 159.8 100 165.1 0.0526 165.2 200 165.8 0.0527 165.9 300 166.1 0.0527 166.2 400 166.8 0.0528 166.9 500 167.4 0.0528 167.5 600 168.1 0.0529 168.2 700 168.8 0.0530 168.9 800 169.5 0.0531 169.6 900 170.4 0.0532 170.5 1000 171.5 0.0533 171.6 1100 172.4 0.0534 172.5 1200 173.7 0.0536 173.8 1300 174.8 0.0537 174.9 1400 176.7 0.0540 176.8 1500 178.1 0.0541 178.2 1600 180.7 0.0544 180.8 1700 182.5 0.0547 182.6 1800 186.1 0.0551 186.2 1900 191.4 0.0557 191.5 UNCLASSIFIED 177 UNCLASSIFIED Table C.5. Uncorrected and corrected distillation data for HRJ-8. Data taken with an atmospheric pressure of 753 mmHg. Volume Distilled (mL) Measured T (˚C) Pressure Correction (˚C) Corrected T (˚C) 0 147.4 0.3531 147.8 100 173.4 0.3750 173.8 200 177.2 0.3782 177.6 300 181 0.3814 181.4 400 184.9 0.3846 185.3 500 189.4 0.3884 189.8 600 195.2 0.3933 195.6 700 200 0.3973 200.4 800 208 0.4040 208.4 900 215 0.4099 215.4 1000 222 0.4158 222.4 1100 230 0.4225 230.4 1200 240 0.4309 240.4 1300 249 0.4385 249.4 1400 256 0.4444 256.4 1500 263 0.4502 263.5 1600 270 0.4561 270.5 1700 275 0.4603 275.5 1800 279 0.4637 279.5 1900 283 0.4670 283.5 UNCLASSIFIED 178 UNCLASSIFIED Table C.6. Uncorrected and corrected distillation data for HRJ-8. Data taken with an atmospheric pressure of 749.6 mmHg. Volume Distilled (mL) Measured T (˚C) Pressure Correction (˚C) Corrected T (˚C) 0 146.2 0.5232 146.7 100 176.2 0.5606 176.8 200 179.9 0.5652 180.5 300 183.3 0.5695 183.9 400 186.5 0.5735 187.1 500 190.4 0.5783 191.0 600 194.1 0.5829 194.7 700 198.4 0.5883 199.0 800 203 0.5940 203.6 900 206 0.5978 206.6 1000 211 0.6040 211.6 1100 216 0.6103 216.6 1200 221 0.6165 221.6 1300 226 0.6228 226.6 1400 231 0.6290 231.6 1500 237 0.6365 237.6 1600 243 0.6440 243.6 1700 251 0.6540 251.7 1800 260 0.6652 260.7 1900 270 0.6777 270.7 UNCLASSIFIED 179 UNCLASSIFIED Table C.7. Uncorrected and corrected distillation data for HRD. Data taken with an atmospheric pressure of 743.8 mmHg. Volume Distilled (mL) Measured T (˚C) Pressure Correction (˚C) Corrected T (˚C) 0 164 0.850 164.8 100 230 0.978 231.0 200 242 1.001 243.0 300 254 1.024 255.0 400 262 1.040 263.0 500 268 1.052 269.1 600 273 1.061 274.1 700 278 1.071 279.1 800 282 1.079 283.1 900 285 1.085 286.1 1000 288 1.091 289.1 1100 290 1.094 291.1 1200 293 1.100 294.1 1300 295 1.104 296.1 1400 297 1.108 298.1 1500 299 1.112 300.1 1600 301 1.116 302.1 1700 304 1.122 305.1 1800 307 1.128 308.1 1900 313 1.139 314.1 UNCLASSIFIED 180 APPENDIX D Database of Potential Surrogate Components 181 This appendix provides the database of compounds used for surrogate development in Chapter 4. Property information was collected from DIPPR and the NIST Webbook as cited in Chapter 4. Most cetane values are from the Murphy Compendium of Cetane Numbers as cited in Chapter 4, except the compendium did not report values for many isoparaffins. For the isoparaffins only, the values reported in Ghosh and Jaffe were used. 182 Table D.1. Selected properties of potential surrogate components – n-paraffins, mono-aromatics, di-aromatics, and naphthenes. Name Molecular class n-Hexane n-Heptane n-Octane n-Nonane n-Decane n-Undecane n-Dodecane n-Tridecane n-Tetradecane n-Pentadecane n-Hexadecane n-Heptadecane n-Octadecane n-Nonadecane n-Eicosane n-Heneicosane Benzene Toluene Ethylbenzene p-Xylene o-Xylene n-Propylbenzene Mesitylene n-Butylbenzene m-Diethylbenzene n-Pentylbenzene Hexylbenzene Cyclohexylbenzene n-Heptylbenzene n-Octylbenzene n-Decylbenzene Methyl Napthalene Ethylcyclohexane n-Propylcyclohexane trans-Decalin n-Butylcyclohexane 1,1-Diethylcyclohexane iso-Butylcyclohexane sec-Butylcyclohexane Bicyclohexyl n-Decylcyclohexane n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic mono-Aromatic di-Aromatic Naphthene Naphthene Naphthene Naphthene Naphthene Naphthene Naphthene Naphthene Naphthene Melting Boiling Molecular Cetane Temperature Temperature Weight Number (g/mol) (˚C) (˚C) 86.18 -95.32 68.88 43.9 100.2 -90.58 98.58 52.5 114.23 -56.77 125.83 64.4 128.26 -53.49 150.97 73 142.28 -29.64 174.31 76.7 156.31 -25.58 196.08 81.1 170.33 -9.58 216.47 82.5 184.36 -5.39 235.62 89.5 198.39 5.86 253.73 95 212.41 9.92 270.84 96 226.44 18.16 287.01 100 240.47 21.98 302.3 105 254.49 28.16 316.86 106.3 268.52 31.89 330.05 110 282.55 36.43 343.93 110 296.57 40.2 356.65 100 78.11 5.53 80.24 6.8 92.14 -94.97 110.78 7.4 106.17 -94.95 136.35 8 106.17 13.26 138.51 1 106.17 -25.17 144.58 9 120.19 -99.6 159.39 7.6 120.19 -44.73 159.39 7.6 134.22 -87.85 183.46 10 134.22 -83.89 181.29 0 148.24 -75 205.61 9 162.27 -61.15 226.26 26 160.26 6.99 240.27 9 176.3 -48 246.25 34.7 190.32 -36 264.55 35.3 218.38 -14.38 298.04 56.8 142.2 -30.48 244.83 0 112.21 -111.31 131.95 34.1 126.24 -94.9 156.9 40 138.25 -30.36 187.46 40 140.27 -74.73 181.13 61.8 140.27 -100 176.82 61.8 140.27 -94.8 171.45 61.8 140.27 -110 179.49 61.8 166.3 3.63 239.19 50 224.43 -1.73 297.75 90 183 Table D.2. Selected properties of potential surrogate components – isoparaffins. Name Isobutane Isopentane 3-Methylpentane 2,3-Dimethylbutane 3-Methylhexane 3-Ethylpentane 2,4-Dimethylpentane 2,2,3-Trimethylbutane 2,5-Dimethylhexane 2,3,4-Trimethylpentane 2-Methyl-3-Ethylpentane 2-Methyloctane 3-Methyloctane 4-Methyloctane 2,4-Dimethyl-3-ethylpentane 3-Ethylheptane 2,7-Dimethyloctane 3,3,5-Trimethylheptane 2,2,5,5-Tetramethylhexane 2,2,3,3-Tetramethylhexane 2-Methylnonane 3-Methylnonane 4-Methylnonane 5-Methylnonane 3-Methylundecane 2,2,4,4,6,8,8Heptamethylnonane Squalane Molecular class Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Isoparaffin Melting Boiling Molecular Cetane Temperature Temperature Weight Number (g/mol) (˚C) (˚C) 58.12 -159.61 -11.57 20 72.15 -159.9 27.99 24.28 86.18 -162.9 63.42 29.09 86.18 -127.96 58.13 11.8 100.2 -119.4 92 32.83 100.2 -118.6 93.62 32.83 100.2 -119.24 80.64 13.94 100.2 -24.58 81.03 13.94 114.23 -91.15 109.26 15.37 114.23 -109.2 113.62 15.37 114.23 -114.95 115.8 15.37 128.26 -80.1 143.45 41.38 128.26 -107.6 144.38 41.38 128.26 -113.2 142.59 41.38 128.26 -122.36 136.87 25.17 128.26 -114.9 143.35 41.38 142.28 -54 160.02 29.8 142.28 -108.15 155.83 29.8 142.28 -12.6 137.61 29.8 142.28 -54 160.46 29.8 142.28 -74.65 167.15 46.73 142.28 -84.8 167.95 46.73 142.28 -98.7 165.85 46.73 142.28 -87.7 165.3 46.73 170.33 -58 210.95 49.58 Isoparaffin 226.44 -110.15 246.5 60 Isoparaffin 422.81 -38 447 60 184 UNCLASSIFIED: Distribution Statement A. Approved for public release. APPENDIX E Near-IR and Fourier Transform IR Spectra for all Fuel Distillations and Mixtures UNCLASSIFIED 185 UNCLASSIFIED Figure E.1. NIR spectra for the JP-8 fuel and its distillation fractions. Figure E.2. NIR spectra for the IPK fuel and its distillation fractions UNCLASSIFIED 186 UNCLASSIFIED Figure E.3. NIR spectra for the HRJ fuel and its distillation fractions Figure E.4. NIR spectra for the SPK2 fuel and its distillation fractions UNCLASSIFIED 187 UNCLASSIFIED Figure E.5. NIR spectra for the HRJ-8 fuel and its distillation fractions Figure E.6. NIR spectra for the SPK fuel and its distillation fractions UNCLASSIFIED 188 UNCLASSIFIED Figure E.7. NIR spectra for the HRD fuel and its distillation fractions Figure E.8. NIR spectra for mixtures of JP-8 and IPK. UNCLASSIFIED 189 UNCLASSIFIED Figure E.9. NIR spectra for mixtures of JP-8 and HRJ. Figure E.10. NIR spectra for mixtures of JP-8 and SPK2. UNCLASSIFIED 190 UNCLASSIFIED Figure E.11. NIR spectra for mixtures of JP-8 and HRJ-8. Figure E.12. NIR spectra for mixtures of JP-8 and SPK. UNCLASSIFIED 191 UNCLASSIFIED Figure E.13. NIR spectra for mixtures of JP-8 and HRD. Figure E.14. FTIR spectra for the JP-8 fuel and its distillation fractions. UNCLASSIFIED 192 UNCLASSIFIED Figure E.15. FTIR spectra for the IPK fuel and its distillation fractions. Figure E.16. FTIR spectra for the HRJ fuel and its distillation fractions. UNCLASSIFIED 193 UNCLASSIFIED Figure E.17. FTIR spectra for the SPK2 fuel and its distillation fractions. Figure E.18. FTIR spectra for the HRJ-8 fuel and its distillation fractions. UNCLASSIFIED 194 UNCLASSIFIED Figure E.19. FTIR spectra for the SPK fuel and its distillation fractions. Figure E.20. FTIR spectra for the HRD fuel and its distillation fractions. UNCLASSIFIED 195 UNCLASSIFIED Figure E.21. FTIR spectra for mixtures of JP-8 and IPK. Figure E.22. FTIR spectra for mixtures of JP-8 and HRJ. UNCLASSIFIED 196 UNCLASSIFIED Figure E.23. FTIR spectra for mixtures of JP-8 and SPK2. Figure E.24. FTIR spectra for mixtures of JP-8 and HRJ-8. UNCLASSIFIED 197 UNCLASSIFIED Figure E.25. FTIR spectra for mixtures of JP-8 and SPK. Figure E.26. FTIR spectra for mixtures of JP-8 and HRD. UNCLASSIFIED 198 APPENDIX F Proof of Pseudocomponent Mixtures for Cetane Number Prediction and Tabulated Cetane Numbers and Beta Values 199 F.1 Proof of Pseudocomponent Model To represent a fuel as one component rather than multiple, it first needs to be shown that a mixture can be represented with a pseudocomponent and still give the same results. To reiterate, the Ghosh and Jaffe equation states that the cetane number of a mixture can be predicted using: 𝐶𝑁 = ∑𝑖 𝑣𝑖 𝛽𝑖 𝐶𝑁𝑖 ∑𝑖 𝑣𝑖 𝛽𝑖 ( F.1 ) where CN is the mixture cetane number, vi is the volume fraction of component i, CNi is the cetane number of pure i, and i is an adjustable parameter. A 4 component system is considered such that: 𝐶𝑁 = 𝑣1 𝛽1 𝐶𝑁1 + 𝑣2 𝛽2 𝐶𝑁2 + 𝑣3 𝛽3 𝐶𝑁3 + 𝑣4 𝛽4 𝐶𝑁4 𝑣1 𝛽1 + 𝑣2 𝛽2 + 𝑣3 𝛽3 + 𝑣4 𝛽4 ( F.2 ) The first 3 components could be condensed into a pseudocomponent mixture, designated by x. The resulting equation is: 𝐶𝑁 = 𝑣𝑥 𝛽𝑥 𝐶𝑁𝑥 + 𝑣4 𝛽4 𝐶𝑁4 𝑣𝑥 𝛽𝑥 + 𝑣4 𝛽4 ( F.3 ) where vx is defined as the sum of the volume fractions in the lump: 𝑣𝑥 = 1 − 𝑣4 = 𝑣1 + 𝑣2 + 𝑣3 ( F.4 ) βx is defined as the volume average of the individual βs such that: 𝛽𝑥 = 𝑣1 𝑣2 𝑣3 1 𝛽1 + 𝛽2 + 𝛽3 = (𝑣1 𝛽1 + 𝑣2 𝛽2 + 𝑣3 𝛽3 ) 𝑣𝑥 𝑣𝑥 𝑣𝑥 𝑣𝑥 The cetane number of the lump is defined using the original equation such that: 200 ( F.5 ) 𝑣1 𝑣2 𝑣3 𝛽 𝐶𝑁 + 𝛽 𝐶𝑁 + 1 1 2 2 𝑣𝑥 𝑣𝑥 𝑣𝑥 𝛽3 𝐶𝑁3 𝐶𝑁𝑥 = 𝑣1 𝑣2 𝑣3 𝑣𝑥 𝛽1 + 𝑣𝑥 𝛽2 + 𝑣𝑥 𝛽3 𝑣1 𝛽1 𝐶𝑁1 + 𝑣2 𝛽2 𝐶𝑁2 + 𝑣3 𝛽3 𝐶𝑁3 = 𝑣1 𝛽1 + 𝑣2 𝛽2 + 𝑣3 𝛽3 ( F.6 ) where in equation ( F.6 ), 1/vx is factored out of both the numerator and denominator, so it cancels out and the equation is simplified. The lump cetane number, CNx, in equation ( F.6 ) can be used in the original equation of Ghosh and Jaffe, represented by equation ( F.1 ). When equations ( F.4 ), ( F.5 ), and ( F.6 ) are plugged into equation ( F.3 ), 𝐶𝑁 1 𝑣 𝛽 𝐶𝑁 + 𝑣2 𝛽2 𝐶𝑁2 + 𝑣3 𝛽3 𝐶𝑁3 𝑣𝑥 [𝑣 (𝑣1 𝛽1 + 𝑣2 𝛽2 + 𝑣3 𝛽3 )] [ 1 1 1 ] + 𝑣4 𝛽4 𝐶𝑁4 𝑣1 𝛽1 + 𝑣2 𝛽2 + 𝑣3 𝛽3 𝑥 = 1 𝑣𝑥 [𝑣 (𝑣1 𝛽1 + 𝑣2 𝛽2 + 𝑣3 𝛽3 )] + 𝑣4 𝛽4 𝑥 𝑣1 𝛽1 𝐶𝑁1 + 𝑣2 𝛽2 𝐶𝑁2 + 𝑣3 𝛽3 𝐶𝑁3 + 𝑣4 𝛽4 𝐶𝑁4 = 𝑣1 𝛽1 + 𝑣2 𝛽2 + 𝑣3 𝛽3 + 𝑣4 𝛽4 ( F.7 ) Once everything has been canceled out, the original equation returns. Therefore, any number of components can be grouped into one pseudocomponent lump using the volume averaging method for  shown above, and the Ghosh and Jaffe equation will still be valid. This also means that a mixture of two fuels can be treated as a binary mixture of two representative pseudocomponents. These pseudocomponents take the form of: 𝑣𝑥 = ∑ 𝑣𝑖 𝑣𝑖 𝛽𝑥 = ∑ 𝛽𝑖 𝑣𝑥 𝑣𝑖 𝑣 𝛽𝑖 𝐶𝑁𝑖 𝐶𝑁𝑥 = ∑ 𝑥 𝛽𝑥 𝑖 201 ( F.8 ) ( F.9 ) ( F.10 ) F.2 Tabulated Cetane Numbers and Blending Values  values are used by each functionality class as indicated by Ghosh and Jaffe referenced in Chapter 4. Table F.1. Cetane numbers and β values for n-paraffins and mono-isoparaffins used for cetane number surrogates. Molecular Class n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin n-Paraffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Mono-Isoparaffin Carbon Number C5 C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 C5 C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 C17 C18 C19 C20 C21 C22 C23 C24 C25 202 Cetane Number 31.72 39.97 54.61 61.56 85.4 80.91 84.74 83.15 88 100.4 96.3 100 27.1135 30.17 33.2265 36.283 39.3395 42.396 45.4525 48.509 51.5655 54.622 57.6785 60.735 63.7915 66.848 69.9045 72.961 76.0175 79.074 82.1305 85.187 88.2435 β Value 0.5212 0.5212 0.5212 0.5212 0.5212 0.5212 0.5212 0.5212 0.5212 0.5212 0.5212 0.5212 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 Table F.2. Cetane numbers and β values for multi-isoparaffins and mono-cylcoparaffins used for cetane number surrogates. Molecular Class Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Multi-Isoparaffin Mono-cycloparaffin Mono-cycloparaffin Mono-cycloparaffin Mono-cycloparaffin Mono-cycloparaffin Mono-cycloparaffin Mono-cycloparaffin Mono-cycloparaffin Mono-cycloparaffin Mono-cycloparaffin Mono-cycloparaffin Carbon Number C5 C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 C17 C18 C19 C20 C21 C22 C23 C24 C25 C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 203 Cetane Number 14.88 17.75 20.61 23.48 26.34 29.21 32.07 34.94 37.8 40.67 43.53 46.4 49.26 52.13 54.99 57.86 60.72 63.59 66.45 69.32 72.18 18.5 24.9 34.1 40 61.8 62 64.5 78.5 83.7 90 90 β Value 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 7.3717 0.0727 0.0727 0.0727 0.0727 0.0727 0.0727 0.0727 0.0727 0.0727 0.0727 0.0727 Table F.3. Cetane numbers and β values for mono- and di-aromatics used for cetane number surrogates. Molecular Class Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Mono-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Di-Aromatic Carbon Number C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 C17 C18 C19 C20 C10 C11 C12 C13 C14 C15 C16 C17 C18 C19 C20 C21 C22 C23 204 Cetane Number 10.7 2.6 7.4 7.6 10 9 9 37.5 39.2 50.1 56.8 56.3 62.8 50.5 40 22.6 0 5 10 14 16 16.5 17 17.5 18 18.5 19 19.5 20 β Value 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 3.1967 UNCLASSIFIED: Distribution Statement A. Approved for public release. APPENDIX G Regression Code and Coefficients UNCLASSIFIED 205 UNCLASSIFIED G.1 R Code This section is written in the programming language used in the R statistical software. A “#” before a line indicates a comment. To use this code, copy this section into an R script, and highlight the desired sections to run. Functionality is built in to work with both FTIR and NIR data sets. #Prior to a run, the “pls” package must be loaded. This package contains all the functions used in #this code beyond the basic math set. #The dataset must be in a “.csv” file (comma separated file) with no headers. Each row contains #the data for an individual run. The first column must be the cetane numbers, with the #absorbance values for the run continuing the row. Every run needs to have the absorbance #values recorded at the same wavelengths, and the data needs to be lined up so that each column #represents the absorbance at that wavelength for all runs. So each experiment needs to be #contained in a single row, with the form “CN(1), AB(1,1), AB(2,1), AB(3,1)” with CN(1) #representing the experimental cetane number for run 1, AB(1,1) representing the absorbance #value at wavelength 1 for run 1, AB(2,1) the absorbance at wavelength 2 for run 1, ect. The #wavelengths or wavenumbers should not be included in the file as headers, the first line should #begin with the first run. #Read in dataset ALL