Empirical Mass and Kinetic Models for the Flash ... - ACS Publications

Apr 17, 2018 - Then, the kinetic model of flash evaporation was obtained according to the empirical mass model. Finally, the analysis of these models ...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Kinetics, Catalysis, and Reaction Engineering

Empirical Mass and Kinetic Models for the Flash Evaporation of NaCl Water Solution Yong Liu, Qiong Luo, Guodong Wang, Siyuan Zhao, Xianlong Li, and Ping Na Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.8b00402 • Publication Date (Web): 17 Apr 2018 Downloaded from http://pubs.acs.org on April 17, 2018

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

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

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

Industrial & Engineering Chemistry Research

Empirical Mass and Kinetic Models for the Flash Evaporation of NaCl Water Solution Yong Liu,*,† Qiong Luo,† Guodong Wang,† Siyuan Zhao, ‡Xianlong Li,† Ping Na‡ †

School of Environmental Science and Engineering, Tianjin University, Tianjin 300350, China

* Corresponding author E-mail: [email protected]

School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China

ABSTRACT: This paper attempts to provide empirical mass and kinetic models for the flash evaporation of sodium chloride (NaCl) aqueous solution based on experimental phenomena. The models have nine parameters and six affecting factors including initial temperature, operation pressure, NaCl mass fraction, solution depth, evaporator diameter and time. Based on a large number of flash evaporation experimental data from various literatures, the mass model parameters were optimized and validated. After optimizing the model parameters with 283 sets of literature experimental data, the average relative error between the model values and the experimental data is about 5.7%. And a statistical method proved the mass model is well posed. The verification with other 215 sets of literature experimental data showed the mass model is in good agreement with flash evaporation phenomena, and the average relative error between the model values and the experimental data is about

ACS Paragon Plus Environment

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

8.3%. Then, the kinetic model of flash evaporation was obtained according to the empirical mass model. Finally, the analysis of these models indicated that the increase of initial temperature or evaporator diameter and the decrease of operating pressure are in favor of evaporation. Although the increase of solution depth can improve the evaporated mass, the corresponding evaporation efficiency will slightly reduce. And the increase of salt content is having a detrimental effect on the evaporation of NaCl water solution. In addition, the influence of salt content on evaporation at higher operating pressure is more obvious than that at lower operating pressure. The above results show that these models proposed in our work have high accuracy, wide practicability and good rationality.

1. INTRODUCTION

Flash evaporation belongs to an evaporation process with the intense heat and mass transfer phenomena. If the temperature of solution is higher than its boiling point, flash evaporation will happen. Flash evaporation will results in partial vaporization of solution and rapid decrease of solution temperature. According to the feeding state, flash evaporation can be divided into spray flash evaporation and static flash evaporation. Flash evaporation has attracted extensive attention because of its importance in chemical separation industry, material industry and energy industry. Nowadays, flash evaporation has been widely used in the field of crystallization1-4, desalination5-7 and concentration8. At present, flash evaporation process has been widely studied. Miyatake et al.9,10

ACS Paragon Plus Environment

Page 2 of 33

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

Industrial & Engineering Chemistry Research

used water solution to study the flash evaporation process. They found that static flash evaporation could be divided into two stages, such as fast evaporation stage and gradual evaporation stage according to the temperature drop rate. Peterson et al.11 used an interferometer and high-speed photography to observe the temperature field in the static flash evaporation of Fron-11. Popov et al.12 had used a fine wire thermocouple with 25um diameter to study the temperature profile near the interface of water at low evaporation rate conditions. The authors found that there exists a temperature jump in the region near the water-vapor interface, moreover, the temperature jump at the interface increases with the increase of evaporation rate. Liu et al.13 used an infrared thermal imager and a high speed camera to investigate the flash evaporation of salt water droplets under vacuum conditions. The results showed that the evaporation rate decreases with the increase of salt content. Saury et al.14 had studied the flash evaporation process of a water film with the initial water height of 15 mm, superheat ranging from 1 to 35 K and initial temperature from 40 to 74 oC. In 2005, Saury et al.15 had investigated the influence of initial solution height and depressurization rate on the flash evaporation of water. Zhao et al.16 had carried out many flash evaporation experiments with the initial NaCl mass fraction varied from 0 to 0.264 and superheat from 1.8 to 43.4 K. They found that both evaporation and steam-carrying effects influence the change of the equilibrium concentration of NaCl. Yan et al.17 had conducted static flash evaporation experiments for NaCl water solution with initial NaCl mass fraction from 0 to 0.26, initial height from 0.1 to 0.4 m, and superheat ranging from 1.8 to 49.5

o

C. The results indicated that NEF

ACS Paragon Plus Environment

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

(non-equilibrium function) with crystallization will reach a low equilibrium value. Wang et al.18 used aqueous NaCl solution as working fluid in their continuous flash evaporation experiment with superheat from 0.3 oC to 11.5 oC to study the bubble morphology and evaporation rate. Yang et al.19 had used salt water with the initial mass fraction of 0.05 to 0.26 to explore the flash evaporation phenomena under different conditions. They believed that, when superheat and initial liquid level are same, the effect of salt mass fraction on evaporation mass is very weak, and the increase of initial liquid level can promote total evaporation mass, but will slightly reduce evaporation efficiency, and evaporation mass will linearly increases with superheat.

For the desalination of seawater or wastewater and the industrial crystallization of NaCl, the accurate control of evaporation process seriously affects the process efficiency and the product quality of salt. Therefore, evaporation kinetic model plays an important role in the process design of desalination and NaCl crystallization. However, nowadays, the evaporation kinetics with wide applicability and good accuracy is very rare. The well-known kinetic theory of gases (KTG) was often used to analyze the evaporation process for decades at molecular level.20 Hertz-Knudsen equation, based on the kinetic theory of gases, can be used to estimate evaporation flux with the interfacial temperatures and evaporation coefficient. Based on Hertz-Knudsen equation, in ANSYS Fluent software, an evaporation-condensation model with the coefficients of 0.1 was given. However, in the actual simulation process, the evaporation coefficient needs to be corrected in order to match the

ACS Paragon Plus Environment

Page 4 of 33

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

Industrial & Engineering Chemistry Research

simulation results with the experimental phenomena. In fact, the evaporation coefficient of water from a large number of literature reports is different and ranges from 0.01 to 1.21 The reliability of the evaporation-condensation model for process simulation and design is often not high unless there is an accurate evaporation coefficient obtained from experiments. In order to avoid the uncertainty of evaporation coefficient, based on the transition probability concept of quantum mechanics, called as a statistical rate theory (SRT), Ward et al.22 had develop an expression to predict evaporation flux. The authors believed the new model values are better than those obtained by the KTG compared with experimental observations. But this model requires accurate vapor temperature and liquid temperature near the interface. Up to now, it is difficult to determine the temperature of interfacial gas-liquid phase.11,12 Therefore, the model is quite inconvenient for engineering application. Furthermore, at present the model is mainly used for single component case such as fresh water or carbon monoxide. For a multi component system such as NaCl water solution, seawater and wastewater, it should be difficult to truly predict evaporation flux using this new model. Some authors had used flash evaporation experimental data to present a correlation. Based on the experiments of fresh water and saline water (3.5% NaCl solution), under the conditions of initial temperature from 25oC to 80oC, superheat from 0.5oC to l0oC, and solution depth from 165 mm to 457 mm, Gopalakrihna et al.23 had develop an empirical correlation of the total mass of vapor released during flash evaporation. Owing to the narrow range of experimental conditions, the applicability of Gopalakrishna’s model should be not

ACS Paragon Plus Environment

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

high. Through many flash evaporation experiments, Saury et al.14 and Yang et al.19 had respectively proposed the formula of final evaporation mass. But these models can't describe flash evaporation process.

In order to provide a more accurate and widely applicable evaporation model for the flash evaporation of NaCl water solution, based on the typical change curve of steam mass versus time, this paper tries to present empirical mass and kinetic models. These models include six influencing factor, such as initial temperature, operation pressure, solution depth, evaporator diameter, salt content and time. Compared with Gopalakrishna’s model or the evaporation-condensation model, the empirical mass and kinetic models presented in this work are more consistent with the experimental phenomena. Finally, in this paper the influence of various factors on evaporation mass and rate will be analyzed under a certain condition. These new evaporation models should be an important tool for the simulation and design of the flash crystallization process of NaCl water solution and the desalination process of seawater or wastewater.

2. EMPIRICAL MASS MODEL

For comparison of the evaporated mass in various literatures and model optimization, the evaporated mass per unit volume of solution will be used in this study. The evaporated mass per unit volume of solution can be calculated using the following formula.

ACS Paragon Plus Environment

Page 6 of 33

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

Industrial & Engineering Chemistry Research

mv =

m AH

(1)

Based on the flash evaporation experiments in the literatures14,24, the typical change curves of vapor mass versus time are shown in Figure 1. According to these phenomena, the relationship between vapor mass and time is in good agreement with exponential function (eq 2). m v = m vf (1-exp(-p 2 t ))

(2)

Based on fresh water and saline water flash evaporation experiments, Gopalakrihna et al.23 had proposed a dimensionless equation of final evaporation mass like eq 3, which is mainly related to solution depth, superheat, salt content and solution properties including density, specific heat, viscosity, thermal conductivity, vaporization heat. m vf =0.8867 ρ l JaJa T 0. 05 Pr − 0.05 ( ∆ P / H) − 0.05 (1 + c) 0.06

(3)

Through various experimental studies on pool brine flash evaporation, Yang et al.19 had also obtained a more simple correlation formula of final evaporation mass, which can be written as eq 4. mvf =0.914ρl Ja0.95 Pr −0.24 (H/D)- 0. 06

(4)

According to pool water flash evaporation experiments under various initial temperature and operation pressure, Saury et al.14,15 had found final evaporation mass can be well expressed as the following equation.

ACS Paragon Plus Environment

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

mvf =ρ lC p ∆T/h

Page 8 of 33

(5)

When solution depth was 15mm, Saury et al.14 found the final evaporation mass under various conditions is proportional to superheat, and the value of ρ l C p / h is a constant of about 1.558kg.m-3.K-1, indicating that solution temperature and operating pressure have weak effect on the value of ρ l C p / h even though the initial temperature ranges from about 40oC to about 74oC, and the operation pressure ranges from 5000Pa to 20000Pa. Later, Saury et al.15 examined the effect of solution depth, ranging from 25mm to 275mm, on the final evaporation mass through many flash evaporation experiments of water at 5000Pa. They found that d m f / d H > 0 and d 2 m f / dH 2 < 0 . For the flash evaporation of NaCl water solution, Yang et al.19 had

also found the similar phenomena mentioned by Saury et al.14,15 According to the law of mass transfer, it is clear that the cross section area of evaporator has an important influence on evaporation. According to the study by Yang et al.19, the effect of salt concentration on evaporation is mainly reflected in superheat. Therefore, based on the analysis above, an empirical formula of final evaporation mass per unit volume of solution is mainly related to superheat, solution depth and evaporator dimeter, and can be simply expressed as a powder function like eq 6. m vf = a 1 ∆ T a 2 H a 3 D a 4

(6)

Relaxation time, defined as the time from non-equilibrium state to equilibrium state, also known as flashing time or bubble time, plays an important role in tuning the boiling regimes.25 Through a large number of flash experiments with water in an

ACS Paragon Plus Environment

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

Industrial & Engineering Chemistry Research

evaporator, Guo et al.26 believed flashing time is mainly affected by superheat and solution depth. Miyatake et al.9,10 thought that flashing time is related to superheat and equilibrium temperature. Gopalakrihna et al.23 believed flashing time is related to superheat, solution depth and solution properties. As sated earlier, the cross section area of evaporator affects evaporation rate, and must affect flashing time. Based on the above viewpoint about flashing time, p2 in eq 2 can be empirically represented as a power function like eq 7. (7)

p 2 = a 5 Tea 6 ∆ T a 7 H a 8 D a 9

Based on eqs 2, 6 and 7, the evaporation mass per unit volume of solution at t moment can be written as: 4 m v = a 1 ∆ T a 2 H a 3 D a( 1 − exp( -a 5 T ea 6 ∆ T a 7 H a 8 D a 9 t)

(8)

3. EVAPORATION MASS CALCULATION METHODS

Saury et al.14,15 estimated the evaporated mass under various time using a simple formula like eq 9.

m = ρ l AH(1 −

1 + Cp T / h 1 + C p T0 / h



(9)

or eq 10   C  m = ρ l AH 1 − exp  − v (T0 -T)    h  

(10)

The above two equations both were neglected the second order terms during the

ACS Paragon Plus Environment

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

Page 10 of 33

formula deduction, which leads to some errors. In order to estimate evaporation mass between i and i+1 moments more accurately, difference equation of eq 11 was used in this work. m i C p Ti -( m i − ∆ m i +1 ) C p Ti+1 = ∆ m i +1 h

i = 0,1, 2 L L

(11)

Rearranging eq 11, we can get

∆ m i +1 =

m i C p (Ti -Ti+1 ) h − C p Ti+1

i = 0,1, 2 LL

(12)

Therefore, the total evaporation mass at i+1 moments can be determined using eq 13. i

mi +1 = ∑ ∆m j+1

(13)

j=0

For the evaporated data in the literatures14,19, eq 9 was firstly used to obtain temperature values, and then eqs 11 ~13 were used to calculated the evaporated mass at various time. According to eq 10, the evaporated mass in the literature15 was corrected using the similar methods. Based on the temperature at various time, the evaporated mass in the literature24,27,28 was calculated using eqs 11 ~13. The experimental data in the literatures14,15,19,24,27,28 were obtained using scale reading method with the relative error of less than 5%. Finally, the evaporated mass per unit volume of solution was calculated according to eq 1. These literature experimental data are shown in Table S1~S4.

ACS Paragon Plus Environment

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

Industrial & Engineering Chemistry Research

4. RESULTS AND DISCUSSIONS

4.1Model Optimization and Verification

4.1.1 Model parameter optimization

Using the methods mentioned earlier, 283 sets of literature experimental data (shown in Table S1~S3) under various conditions are shown in Figure 1 (a) ~ (d) and Figure 2, respectively. Among them, the literature14 provided 156 sets of experimental data for fresh water, 49 sets of experimental data for aqueous salt solution were got from the literature19 and other 78 sets of experimental data for fresh water were obtained from the literature24.

Based on eq 8 and 283 sets of literature experimental data, the nine parameters were obtained using nonlinear optimization techniques. According to the physical meaning of parameters, the values of a1 and a5 must be greater than 0. Based on the basic law of evaporation,9,10,14,15,19,24 the values of a2, a4, a6 and a9 should be larger than 0, the values of a3, a7 and a8 should be less than 0, moreover, a3 should be in the range of -1 to 0 according the flash evaporation phenomena from literatures.15,19 With these restrictive conditions of the model parameters, after optimization, the optimal values of a1~a9 are 2.05, 9.72e-1, -5.54e-2, 8.61e-2, 1.00e+4, 1.95, -2.16e-1, -1.09e-2 and 7.85, respectively.

ACS Paragon Plus Environment

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

Page 12 of 33

(a)

(b)

(c)

(d)

Figure 1. The change of the evaporated mass of experiments and eq 8 with time during water flash evaporation. (a) H=0.06m;24 (b) ~ (d) H=0.015m.14

The values calculated by eq 8 are shown in Figures 1 and 2, respectively. The

ACS Paragon Plus Environment

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

Industrial & Engineering Chemistry Research

overall fitting effect of the experimental data14,19,24 is shown in Figure 3. mv,E and mv,C in Figure 3 represent the experimental data and the model values, respectively. There is a good linear relationship between the experimental data and the model values, and the slope and the regression index (R2) are 0.997 and 0.989(seen in Figure 3), respectively. The average relative error between 283 sets of experimental values and the calculated data is about 5.7%. The results indicated that the values calculated by eq 8 are in very good agreement with the experimental data.

Figure 2. Experimental data from literature19 and the values estimated by eq 8 under various conditions for NaCl water solution when t=20s.

ACS Paragon Plus Environment

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

Figure 3. The overall fitting effect of the experimental data.14,19,24

In addition, a statistical method has been used to analyze the well posedness of eq 8, and the model statistics are shown in Table 1. Mp and M are the number of the model parameters and the number of experimental data sets, respectively, and p2 is equal to (1- sum of residuals squares/ sum of experimental data squares), F can be calculated with the expression of p2/Mp/((1-p2)/(M-Mp)), and FT is the F value under the corresponding freedom degree with a significant level of 5%.29 It is generally believed that the model is well posed when p2 is greater than 0.9 and F is greater than 10 times of FT.29 According to Table 1, it is clear that the model statistics meet the above requirements, so the mass model of eq 8 and its parameters are quite nice. Table 1. The statistics of eq 8 Mp M-Mp 9

274

Re

p2

F

10×FT

765 0.997 9695

19.1

*Re =Sum of residuals squares

4.1.2 Model verification

The previous analysis shows that the proposed mass model has good regression effects. In order to show its applicability, this section will use other experimental data to verify the model. The flash evaporation experiments in the literature27 and the literature15 both were carried out using fresh water, but the former was square flash chamber of 0.125m×0.13m, the initial solution height was 0.05m, the later was cylindrical chamber with a diameter of 0.299m, whose solution level was range from

ACS Paragon Plus Environment

Page 14 of 33

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

Industrial & Engineering Chemistry Research

about 0.025m to about 0.275m. The flash evaporation experiments for NaCl water solution in the literatures19, 28 were run 20s and the experimental systems were same, but the experimental conditions are different. All of the experimental data used here were not included in the previous optimization process. Among them, 59 sets of experimental data were got from the literature15, 57 sets of experimental data were got from the literature27 and 30 sets of experimental data were obtained from the literature19. 69 sets of experimental data were obtained from the literature28. The number of the sets of experimental data is 215 in total (shown in Table S4).

Figure 4 (a) shows the experimental data and the calculated data by eq 8. The average relative error between the experimental data and the model values calculated by eq 8 is about 8.3%. The slope of the linear function and R2 are 1.027 and 0.967 (Figure 4(a)), respectively. This result proves that the proposed mass model is very nice, and has good practicability.

In order to compare with the effects of other existing evaporation models, Gopalakrihna’s formula23 and the evaporation-condensation model with the coefficient of 0.1, which has been used in Fluent software, both were used to estimate the

evaporation

mass

under

the

literature

conditions.

For

the

evaporation-condensation model, numerical method was used to calculate the evaporation mass per unit volume of solution with the time step of 1.0e-4s. It was found that the errors of the calculated results can be ignored when the time step is less than or equal to 0.0001s. Specially noted, the condensation effect was ignored during the numerical calculation of the evaporation-condensation model because that the ACS Paragon Plus Environment

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

Page 16 of 33

evaporation rate under those conditions was far greater than the condensation rate.

The corresponding results calculated by the two evaporation models are shown in Figure 4 (b) and (c), respectively. Obviously, the linear relationships between the experimental data and the model values are very poor compared with Figure 4 (a). Figure 4(b) indicates the slope and R2 are 0.778 and 0.722, respectively. The average relative error between the experimental data and the calculated values by Gopalakrihna’s model is about 38.5%. Figure 4(c) shows the slope and R2 are 0.845 and 0.851, respectively. The average relative error between the experimental data and the calculated values by the evaporation-condensation model is about 27.3%. Obviously, the average relative errors of the two evaporation models are largely exceed the allowable maximum error of about 20% in chemical engineering design industry. However, the empirical mass model proposed in this work fully meets the error requirements.

(a)

(b)

(c)

Figure 4. The experimental data15,19,27,28 and the model values. (a) the model of eq 8; (b) Gopalakrihna’s formula; (c) the evaporation-condensation model.

ACS Paragon Plus Environment

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

Industrial & Engineering Chemistry Research

In comparison, the effect of the evaporative-condensation model seems to be better than that of Gopalakrihna’s model. The reason may be related to the narrower scope of Gopalakrihna’s experiment conditions23. In the process of model calculation, it was found that, for Gopalakrihna’s model, the effect of liquid level on evaporation is contrary to the literature experimental phenomena.15,19

4.2 Analysis of Mass and Kinetic Models

According to eq 8 and its parameter values, the empirical mass model for the flash evaporation of water or NaCl water solution can be expressed as follows: m v =2.05∆T 0.972 H -0.0554 D 0.0861 (1 − exp(-10 4 Te1.95 ∆ T -0.216 H -0.109 D 7.85 t)

(14)

Then, the evaporation kinetic model can be easily got after taking the time derivative of eq 14, which can be written as follows: r =2.05 × 10 4 T e 1.95 ∆ T 0.756 H -0.0663 D 7.94 exp(-10 4 T e 1.95 ∆ T -0.216 H -0.0109 D 7.85 t )

(15)

Based on eqs 14 and 15, the influence of various affecting factors on evaporation mass and rate can be analyzed. In order to investigate the effect of operating pressure and salt content on superheat, a formula for calculating the boiling point under various conditions must be needed. Based on the boiling point diagram of NaCl water solution under different salt concentration and pressure conditions,30 the Antoine equation of water,31 the boiling point of salt water solution at atmospheric pressure32 and Duhring’s rule, the boiling point of NaCl water solution can be expressed as follows:

ACS Paragon Plus Environment

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

Te =

Page 18 of 33

1657.46(1.00+0.32c-1.60c 2 +5.79c 3 ) -227.02-93.37c+571.27c 2 -1836.59c 3 (16) 7.07406-lg(P/1000)

Then the superheat can be calculated using eq.17. ∆ T=T0 -Te

(17)

Based on eqs 14~17, the effects of initial temperature, operation pressure, NaCl mass fraction, solution depth, evaporator diameter and time on the evaporation mass and rate of NaCl solution under various conditions can be simulated. The results are shown in Figures 5~9, respectively.

Based on Figure 5, it is clear that, when other factors are fixed, the higher the initial temperature is, the greater the final evaporation mass is (Figure 5(a)). In addition, at the same operating time, evaporation rate increases with initial temperature (Figure 5(b)). This result is consistent with the experimental phenomena of the literature14,19,24. The reason is that when other conditions are fixed, the increase of initial temperature directly leads to the increase of superheat.

(a)

ACS Paragon Plus Environment

(b)

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

Industrial & Engineering Chemistry Research

Figure 5. The effect of temperature on evaporation mass and rate under the conditions of P=15000Pa, c=0.1, H=0.1m and D=0.15m.

(a)

(b)

Figure 6. The effect of pressure on evaporation mass and rate under the conditions of T0=75oC, c=0.1, H=0.1m and D=0.15m According to Figure 6, when other influencing factors are fixed, the decrease of operating pressure will lead to a larger final evaporation mass (Figure 6 (a)). The reason is that the decrease of pressure leads to the decrease of the saturation temperature of solution, which causes the increase of the superheat. This phenomenon is why the evaporation of brine solution is usually operated at lower pressure, and is also consistent with a large number of literature experiments.14,15 For evaporation rate, the effect of pressure is shown in Figure 6(b). In the earlier stage of evaporation, the evaporation rate curve for 10000Pa is intersected with the rate curve of 15000Pa (Figure 6(b)). However, when operating pressure is 20000Pa, the rate curve is always below those for P=10000Pa and P=15000Pa in the time range of 0 to 1s (Figure 6 (b)).

ACS Paragon Plus Environment

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

(b)

(a)

Figure 7. The effect of solution depth on evaporation mass and rate under the conditions of T0=75oC, P=15000Pa, c=0.1 and D=0.15m

(a)

(b)

Figure 8. The effect of evaporator dimeter on evaporation mass and rate under the conditions of T0=75oC, P=15000Pa, c=0.1 and H=0.1m From Figure 7(a) and (b), the increase of solution depth is beneficial to the final evaporation mass and the evaporation rate when other factors keep at constant. This phenomenon is in agreement with the results of the literature experiments.15,19,24 The fundamental reason is that higher solution depth leads to more heat energy, which can

ACS Paragon Plus Environment

Page 20 of 33

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

Industrial & Engineering Chemistry Research

be used to vaporize more of solution at a certain saturation temperature. When solution depth is 0.1m, 0.2m and 0.3m, respectively, the corresponding final evaporation mass is 0.0654kg, 0.1259kg and 0.1846kg. The ratios of the final evaporation mass for 0.2m and 0.3m to that of 0.1m are 1.925 (less than 2) and 2.823 (less than 3), respectively. Clearly, the increase of solution depth will slightly reduce evaporation efficiency.

(a-1)

(b-1)

ACS Paragon Plus Environment

(a-2)

(b-2)

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

(c-2)

(c-1)

Figure 9. The effect of salt content on evaporation mass and rate at various pressure when T0=75oC, H=0.1m and D=0.15m. (a) P=5000Pa, (b) P=15000Pa, (c) P=30000Pa

Figure 8 (a) indicates that, when other conditions are constant, if evaporator diameter is large, the final evaporation mass will be high. The reason is that big evaporator diameter will result in large solution volume, similar to the effect of liquid depth, it must have more energy for the vaporization of solution. Compared with Figure 7(a), the increase of evaporator diameter can not only increase evaporation mass, but also shorten flashing time. That's why people often want to increase the diameter of evaporator. There is a cross phenomenon between the corresponding evaporation rate curves under different diameters (Figure 8(b)). This is because the

ACS Paragon Plus Environment

Page 22 of 33

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

Industrial & Engineering Chemistry Research

larger the diameter is, the higher the evaporation efficiency is. The higher evaporation efficiency leads to a rapid decrease in evaporation rate, eventually resulting in the cross phenomenon among the rate curves.

Figure 9 shows that salt mass fraction has influence on evaporation mass and rate. Moreover, operating pressure will obviously affect the influence of salt mass fraction on evaporation. When P=5000Pa, salt mass fraction, such as 0, 0.1 and 0.2, has a weak impact on evaporation mass and rate. After five second, the final evaporation mass is about 0.1322kg, 0.1317kg and 0.1282kg, respectively. The ratios of final evaporation mass for the salt mass fraction of 0 and 0.1 to that of the salt mass fraction of 0.2 are 1.031 and 1.027, respectively. And the evaporation rate curves almost coincides with each other. However, when operating pressure is 30000Pa, the final evaporation mass is distinctly different. After five second, the corresponding final evaporation mass for the salt mass fraction of 0, 0.1 and 0.2 is 0.0223kg, 0.0187kg and 0.0112kg, respectively. The corresponding ratios calculated using the similar method mentioned above are 1.991 and 1.670, respectively. It is very clear that salt concentration has a great influence on evaporation amount when operating pressure is higher. Of course, high operation pressure is not conducive to solution evaporation (Figure 6). The reason for this phenomenon should be related to the effect of pressure and salt concentration on the boiling point of NaCl water solution.

5. CONCLUSIONS

Based on the typical flash evaporation phenomena and experimental data of fresh

ACS Paragon Plus Environment

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

water and NaCl water solution in many literatures, this paper presents an empirical evaporation mass and kinetic models for the flash evaporation of NaCl water solution. The results proved that this mass model and its parameters are well posed, and can be used to accurately simulate the evaporation mass of aqueous NaCl solution. For 498 sets of literature flash evaporation experimental data, the average relative error between the model values and the experimental data is less than about 8.3%. The calculation accuracy of the models proposed in our work is much better than that of Gopalakrihna’s model or the evaporative-condensation model. The model analysis shows that the model is well consistent with the basic laws of the flash evaporation of water and NaCl water solution. The increase of initial solution temperature, solution depth and evaporator diameter can improve vapor mass. But the increase of liquid depth will reduce evaporation efficiency. However, the increase of the diameter of evaporator is beneficial to the rapid and highly efficient flash evaporation of solution. The increase of operating pressure and salt content will be detrimental to solution flash evaporation. In addition, the operating pressure have some influence on the effect of salt concentration on evaporation mass and rate.

ASSOCIATED CONTENT

Supporting Information Available: the literature evaporation mass under various conditions are available free of charge via the internet at http://pubs.acs.org.

■AUTHOR INFORMATION

ACS Paragon Plus Environment

Page 24 of 33

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

Industrial & Engineering Chemistry Research

Corresponding Author

*E-mail: [email protected].

ORCID

Yong Liu: 0000-0002-6551-1608

■Notes

The authors declare no competing financial interest.

■ACKNOWLEDGMENTS

This work was supported by the National Natural Science Foundation of China (No.51478308), the Scientific Research Special Fund of Marine Public Welfare Industry (No. 201405008) and the Natural Science Foundation of Tianjin (No.14JCYBJC23300).

■NOMENCLATURE A=horizontal cross-sectional area, m2

a1~a9=correlation parameters c=salt mass fraction Cp=isobaric specific heat capacity, J·kg−1·oC−1 Cv=equivalent specific heat capacity, J·kg−1·oC−1 D=evaporator diameter, m

ACS Paragon Plus Environment

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

H=initial height of liquid, m h=latent heat of vaporization, J·kg-1

Ja=Jakob number (=Cp∆T/h) JaT= ρl Cp∆T/ρvh Pr= Prandtl number of the liquid

m=mass, kg mv=evaporated mass per unit volume, kg·m-3 mf=final evaporation mass, kg mvf=final evaporation mass, kg·m-3 P=operation pressure, Pa T=solution temperature, oC T0=initial solution temperature, oC Te= saturation temperature, oC t=time, s r=evaporation rate per unit volume of solution, kg·m-3·s-1 R= evaporation rate, kg·s-1

Greek Symbols

ACS Paragon Plus Environment

Page 26 of 33

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

Industrial & Engineering Chemistry Research

ρl = solution density, kg·m-3 ρv= vapor density, kg·m-3 ∆mi+1= steam mass between i and i+1 moments, kg ∆P= pressure drop, Pa ∆T= initial superheat, oC

■REFERENCES

(1) Kiew, M. K.; Moses, O. T. Nonlinear control of a simulated industrial evaporation system using a feedback linearization technique with a state observer. Ind. Eng. Chem. Res. 1999, 38, 2995–3006.

(2) Kiew, M. K.; Moses, O. T.; Gade, P. R.; Yu, C. T. Strategies for enhancing geometric nonlinear control of an industrial evaporator system. Ind. Eng. Chem. Res. 2001, 40, 656–667.

(3) Kim, J.W.; Shin, M.S.; Kim, J.K.; Kim, H.S.; Koo, K. K. Evaporation crystallization of RDX by ultrasonic spray. Ind. Eng. Chem. Res. 2011, 50, 12186– 12193.

(4) Zijlema, T.G.; Geertman R.M.; Witkamp, G.J.; Rosmalen, G.M.; Graauw, J. Antisolvent crystallization as an alternative to evaporative crystallization for the production of sodium chloride. Ind. Eng. Chem. Res. 2000, 39, 1330–1337.

(5) Ramon, G.; JoséMaría, P.; Mahmoud, M. Optimal design of water desalination

ACS Paragon Plus Environment

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

systems involving waste heat recovery. Ind. Eng. Chem. Res. 2017, 56, 1834−1847.

(6) Jin, Z.; Ye, H.; Wang, H.; Li, H.; Qian, J. Thermodynamic analysis of siphon flash evaporation desalination system using ocean thermal energy. Energy Convers. Manage. 2017, 136, 66–77.

(7) Paden, J.J.R.; Tewari, P.K.; Venkatram, D.; Barnabas, D. Spray flash evaporator for low-temperature saline water desalination application. Int. J. Nuclear Desalination 2007, 2, 393–405.

(8) Macdonald, R.J.; Middlewood, P.G.; MacManus,R.L. Development and characterisation of a vacuum flash evaporator for concentrating a heat sensitive aqueous peptide stream. Desalination 2008, 218, 238-247.

(9) Miyatake, O.; Murakami, K.; Kawata, Y. Fundamental experiments with flash evaporation. Heat Transfer Jpn. Res. 1973, 2, 89–100.

(10) Miyatake, O.; Fujii, T.; Hashimoto, T. An experimental study of multi-stage flash evaporation phenomena. Heat Transfer Jpn. Res. 1977, 6, 25–35.

(11) Peterson, R.J.; Grewal, S.S.: El-Wakil, M.M. Investigations of liquid flashing and evaporation due to sudden depressurization. Int. J. Heat Mass Transfer 1984, 27, 301– 310.

(12)Popov, S.; Melling, A.; Durst, F.; Ward, C.A. Apparatus for investigation of evaporation at free liquid–vapour interfaces. Int. J. Heat Mass Transfer 2005, 48, 2299–2309.

ACS Paragon Plus Environment

Page 28 of 33

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

Industrial & Engineering Chemistry Research

(13) Liu, L.; Bi, Q.; Li, H. Experimental investigation on flash evaporation of saltwater droplets released into vacuum. Microgravity Sci. Technol 2009, 21, S255– S260.

(14) Saury, D.; Harmand, S.; Siroux, M. Experimental study of flash evaporation of a water film. Int. J. Heat Mass Transfer 2002, 45, 3447–3457.

(15) Saury, D.; Harmand, S.; Siroux, M. Flash evaporation from a water pool: Influence of the liquid height and of the depressurization rate. Int. J. Heat Mass Transfer 2005, 44,953–965.

(16) Zhao, B.; Zhang, D.; Yang, Q.; Yan, J. Experimental study on equilibrium water film concentration in static flash evaporation of aqueous NaCl solution. Desalination 2014, 353, 109–117.

(17) Yang, Q.; Zhao, B.; Zhang, D.; Wang, Y.; Yan, J. Experimental study on heat transfer characteristics in static flash. Int. J. Heat Mass Transfer 2016, 102, 1093– 1099.

(18) Wang, Y.; He, Q.; Yu, L.; Yan, J. Experimental investigation of flash speed of NaCl solution in circulatory flash evaporation. Journal of Engineering Thermophysics 2017, 38, 1024–1028(in Chinses).

(19) Yang, Q.; Liu, G.; Zhang, D.; Zhao, B.; Yan. J. Evaporation static flash

evaporation

of aqueous

characteristics of

NaCl solution. CIESC Journal 2013, 64,

4068–4073(in Chinese).

ACS Paragon Plus Environment

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

(20) Pao, Y. P. Application of kinetic theory to the problem of evaporation and condensation. Phys. Fluids 1971, 14,306 − 312.

(21) Eames, I. W.; Marr, N. J.; Sabir, H. The evaporation coefficient of water: a review. Int. J. Heat Mass Transfer 1997, 40, 2963 − 2973.

(22) Ward, C.A.; Fang, G. Expression for predicting liquid evaporation flux: statistical rate theory approach. Phys. Rev. E 1999, 59, 429 − 440.

(23) Gopalakrihna, S.; Purshothaman, V. M.: Lior, N. An experimental study of flash evaporation from liquid pools. Desalinution 1987, 65, 139−151.

(24) Guo, Y.; Wei, D.; Yan, J.; Wang, G. Influence of initial conditions on transient flash process. Journal of Engineering Thermophysics 2008, 29, 1335-1338(in Chinese).

(25) Liu, Y.; Olewski, T.; Vechot, L. N. Modeling of a cryogenic liquid pool boiling by CFD simulation.Journal of Loss Prevention in the Process Industries 2015, 35, 125-134.

(26) Guo, Y.L.; Li, Yan, J.J.; Deng, W. Liu, J.P.; Chong, D.T. Experimental investigation on flash duration time. Journal of Engineering Thermophysics 2009, 30, 1339-1342.

(27) Yan, J.; Zhang, D.; Wei, D.; Guo, Y. Experimental investigation of the instantaneous mass transfer coefficient at steam-liquid interface during water film flash evaporation in closed chamber. Journal of Xian Jiao Tong University 2008, 42,

ACS Paragon Plus Environment

Page 30 of 33

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

Industrial & Engineering Chemistry Research

515-518 (in Chinese).

(28) Zhang, D.; Chong, D.; Yan, J.; Zhao, B. Experimental study on static flash evaporation of aqueous NaCl solution. International Journal of Heat and Mass Transfer 2012, 55, 7199-7206.

(29) Song, W.; Zhu, B.; Wang, H.; Zhu, M.; Sun, Q.; Zhang, J. Reaction kinetics of methanol synthesis in the presence of C301 Cu-based catalyst (I): model of intrinsic kinetics. Journal of Chemical Industry and Engineering 1988, (4), 401- 408(in Chinese).

(30) Liu, G.; Ma, L.; Liu, J. Handbook of Physical Property Data of Chemistry and Chemical Engineering. Beijing: Chemical Industry Press, 2002(in Chinese).

(31) Wang, S.; Cheng, H. Calculation of saturated vapor pressure of water. Henan Chemical Industry 1999, (11), 29-34 (in Chinese).

(32) Xi, H. Physical property relationship of sodium chloride solution. Journal of Tianjin Institute of Light Industry 1997, (2), 72-74(in Chinese).

ACS Paragon Plus Environment

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

Figure 1. The change of the evaporated mass of experiments and eq 8 with time during water flash evaporation. (a) H=0.06m;24 (b) ~ (d) H=0.015m.14 Figure 2. Experimental data from literature19 and the values estimated by eq 8 under various conditions for NaCl water solution when t=20s. Figure 3. The overall fitting effect of the experimental data.14,19,24 Figure 4. The experimental data15,19,27,28 and the model’s values. (a) the model of eq 8; (b) Gopalakrihna’s formula; (c) the evaporation-condensation model with the coefficient of 0.1.

Figure 5. The effect of temperature on evaporation mass and rate under the conditions of P=15000Pa, c=0.1, H=0.1m and D=0.15m.

Figure 6. The effect of pressure on evaporation mass and rate under the conditions of T0=75oC, c=0.1, H=0.1m and D=0.15m Figure 7. The effect of solution depth on evaporation mass and rate under the conditions of T0=75oC, P=15000Pa, c=0.1 and D=0.15m Figure 8. The effect of evaporator dimeter on evaporation mass and rate under the conditions of T0=75oC, P=15000Pa, c=0.1 and H=0.1m Figure 9. The effect of salt content on evaporation mass and rate at various pressure when T0=75oC, H=0.1m and D=0.15m. (a) P=5000Pa, (b) P=15000Pa, (c) P=30000Pa Table 1. The statistics of eq 8

ACS Paragon Plus Environment

Page 32 of 33

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

Industrial & Engineering Chemistry Research

ACS Paragon Plus Environment