Thermodynamic models based on Pitzer-NRTL and Pitzer-Margules

0 downloads 0 Views 370KB Size Report
Apr 20, 2015 - The thermodynamic model of the extraction of W with primary amine under near neutral conditions is reported in this paper. The activity ...
SCIENCE CHINA Technological Sciences • Article •

May 2015 Vol.58 No.5: 935–942 doi: 10.1007/s11431-015-5819-y

Thermodynamic models based on Pitzer-NRTL and Pitzer-Margules equations for the extraction of tungstic acid with primary amine N1923 LIN Xiao1,3, NING PengGe1,3*, XU WeiFeng1,3, CAO HongBin1,3 & ZHANG Yi1,2,3 1

Research Centre for Process Pollution Control of Beijing, Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100190, China; 2 National Engineering Laboratory for Hydrometallurgical Cleaner Production Technology, Chinese Academy of Sciences, Beijing 100190, China; 3 Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100190, China Received December 27, 2014; accepted March 20, 2015; published online April 20, 2015

The thermodynamic model of the extraction of W with primary amine under near neutral conditions is reported in this paper. The activity coefficients of the nonelectrolytes in organic phase are calculated by the Margules and NRTL equations which are based on previously tested liquid-liquid extraction equilibrium data in combination with mass balances and charge balance formula. The activity coefficients of the electrolytes in aqueous phase are calculated by the Pitzer equation. The thermodynamic model of the extraction of W by a primary amine is constructed from the calculated activity coefficients of electrolytes and nonelectrolytes. The extraction of W using primary amine is also predicated, and the data is compared to the calculated results of the thermodynamic model. It is concluded that Margules equation proves to be suitable and reliable for calculating the activity coefficients of nonelectrolytes in complex organic phase systems. primary amine extraction, separation of Mo from W, thermodynamic, model parameters Citation:

Lin X, Ning P G, Xu W F, et al. Thermodynamic models based on Pitzer-NRTL and Pitzer-Margules equations for the extraction of tungstic acid with primary amine N1923. Sci China Tech Sci, 2015, 58: 935942, doi: 10.1007/s11431-015-5819-y

1 Introduction Tungsten is a rare, noble metal that has certain properties applicable in many fields, especially in defense and high technology areas. Tungsten is commonly associated with some metals [1,2], specifically molybdenum. Both elements exhibit great similarity in atomic structure, ionic radius, et al., which results in great difficulty of separation due to the strong lanthanide contraction [3]. The quantity of Tungsten is limited on the Earth; thus, the separation of tungsten and molybdenum is highly important in extractive hydrometal*Corresponding author (email: [email protected].) © Science China Press and Springer-Verlag Berlin Heidelberg 2015

lurgy. These separations are usually based on aqueous solution chemistry, which relies on differences in the composition, aqueous species, and coordination ability with other ions. Separation methods, such as precipitation [4,5], ion exchange [6], adsorption methods [7], and extraction [8,9], have been explored and applied in industrial productions. Ion floatation and liquid membrane separation have also been reported as new and improved separation methods. Extraction is one of the most promising separation methods due to its great success in generating high purity products as well as being easy to operate [10,11]. Recently, the use of a primary amine to extract and separate tungsten and molybdenum under near neutral conditions has attracted much attention. Yu [12] reported on the mechanism and impact tech.scichina.com link.springer.com

936

Lin X, et al.

Sci China Tech Sci

factors of the separation of tungsten and molybdenum via a primary amine. A primary amine can efficiently separate W from the high Mo-low W system under neutral and weak alkaline conditions, where the removal yield of W can reach 90% through the 4-stage cross-current extraction. The extracted complex of primary amine (RNH2) extraction of W forms a structure of (RNH2)4·(H2WO4)3 as demonstrated by the slope, saturated capacity, and molar series methods. Our report [8] discusses the detailed research of using a primary amine to extract and separate Mo and W in order to obtain high purities of sodium molybdate and sodium tungstate. The experimental results demonstrate that the primary amine is suitable for a high Mo-low W solution in which the W/Mo mass ratio in the extraction raffinate can reach 2.876×105 through the 2-stage, cross-current extraction. These results are in agreement with Chinese standards that state in which the ratio of W in sodium molybdate is lower than 0.0001. The complexity of the aqueous solution chemistry of W and Mo [13,14] usually lead to the formation of various homopolymolybdate ions or miscellaneous polyacid ions under acidic conditions [15,16]. Such conditions increase the difficulty of separating Mo and W. There are few reports concerning the thermodynamic data of W/Mo [17] and the extraction system of W/Mo using a primary amine. A previous study [18] reported on the activity coefficients of the ions and ions in aqueous and organic phase that were calculated using Pitzer equations. Various interaction energy coefficients of molecules-molecules and ions-molecules were also obtained in the study. However, a certain error occurred between the predicated and obtained results for the extraction of W via the primary amine. Although the experiment was successful, the calculation process requires much improvement. In this work, the activity coefficients of nonelectrolytes in organic phase are obtained from the Margules and NRTL equations. These equations are used in conjunction with various equilibrium data of liquid-liquid extraction of W with primary amine and also combined with mass balances and charge balance formulas. Additionally, the Pitzer equation is used to obtain the activity coefficients of electrolytes in aqueous phase. Considering a thermodynamic model of tungsten extraction with a primary amine has been established, this work is aimed to achieve improved predication of the effect a primary amine has on the extraction of W.

May (2015) Vol.58 No.5

Chinese Academy of Sciences. The molar mass of N1923 is 310.30 g/mol. Dilute sulfuric acid was prepared with concentrated sulfuric acid (95%–98%, Sinopharm Chemical Reagent Co., Ltd., China) and ultra-pure water, homemade. All reagents were used without further purification. 2.2

Experimental apparatus

The apparatus included an electronic balance ML104 (Mettler Toledo, Switzerland) with an accuracy of ±0.0001 g, a pH meter Delta 320 (Mettler, Switzerland) with an accuracy of ±0.01, Thermostat DC-0510 (Ningbo ScientzBiotechnology Poltroon Technologies Inc., China) with an accuracy of ±0.1°C, magnetic agitator 84-1A (Shanghai SileInstrument Co., Ltd. China), a 50 mL self-made jacketed equilibrium still, pipette tube, PTFE separating funnel, beaker, and pipette gun. The analytical instruments used were an OPTIMA 5300DV inductively couple plasma-optical emission spectrometer (ICP-OES, Perkin-Elmer, USA) and ion chromatograph DX-500 (Dionex, USA). The experimental setup is shown in Figure 1, and includes a thermostat, magnetic stirrer, equilibrium still, and rubber pipe. 2.3

Experimental procedure

The extraction equilibrium method is as follows. Various concentrations of Na2WO4 aqueous solutions and primary amine solution were first prepared. The initial pH of aqueous solution of tungsten was adjusted to a desired value using dilute sulfuric acid, which was recorded as the initial value. Samples containing equal volumes of 10 mL of an aqueous phase and organic phase were first placed in the equilibrium vessel at 293.15 K for 10 min by pipette, rigorously stirred for 1 hr at rates of 1030 r/min, and allowed to settle for more than 1 hr. Thereafter, samples were transferred into the PTFE separating funnel. After phase disengagement, the aqueous phase was separated, and the equilibrium pH

2 Experiments 2.1

Materials

Sodium tungstate dehydrate (analytical grade) and toluene (analytical grade) were obtained from Beijing Chemical Works, China. The primary amine N1923, (C19–C23) secondary alkyl primary amine, with a purity of 98%, was purchased from Shanghai Institute of Organic Chemistry,

Figure 1

Apparatus for the liquid-liquid extraction equilibrium.

Lin X, et al.

Sci China Tech Sci

was measured using a pH meter. Analytical methods: the aqueous solution and raffinate were first diluted followed by analysis of the W concentration by ICP-OES with power of 1.3 kW, assistant gas flow of 0.2 L/min, cooling air flow of 15 L/min, and carrier gas flow of 0.8 L/min. The active wavelength used for tungsten analysis was 207.912 nm. The concentration of SO42 in aqueous solution was determined by ion chromatography with eluent conditions of 1.1 g/L NaOH, eluent flow rate of 1 mL/min, AC11 detection column, ED50 detector, and 100 mA suppressor. 2.4

Models

2.4.1 Chemical equilibrium Ning et al. [8] reported that when the initial mass concentration of W is less than 1.0 g/L and the equilibrium pH of the aqueous phase is approximately 7.00, the composition proportion of a complex molecule can be determined by the slope method. The chemical equilibrium reaction is described as follows:   H 2 WO 4  RNH 2org. 2H + +WO 24   RNH org.  2  K

(1)

where RNH2 is primary amine N1923, and org. is the organic phase. The corresponding thermodynamic equilibrium constant, K, is expressed as:

K

aHorg. 2 WO4  RNH 2 org. aH2  aWO 2 aRNH 2 4

=

mHorg.  Horg. 2 WO 4  RNH 2 2 WO 4  RNH 2 org. org.  RNH aH2   mWO 2  WO 2  mRNH 2 2 4

,

(2)

4

where ai is activity of species i, mi is molality of species i, and i is activity coefficient of species i. 2.4.2 Mass and charge balance Before the extraction process, dilute sulfuric acid was added to the feed solution to adjust the initial pH of the aqueous phase. The entire process is described by the following equilibria reactions.

937

May (2015) Vol.58 No.5

When T=293.15 K, Ks is 0.01234. The equilibrium constant Ks can be described according to eq. (2):

Ks 

aH aSO2

=

4

aHSO

aH mSO2  SO2 4

4

mHSO  HSO

4

4

mSO2  SO2 4

4

mHSO  HSO 4

=

4

Ks . aH

(10)

Because the equilibrium pH of the raffinate is close to 7.00 and the aqueous phase is relatively dilute, the value of  SO2  HSO can be considered 1 after phase equilibrium 4

4

under neutral conditions. The molality of SO42 is nearly 1234500 times as much as that of HSO4 after calculating eq. (10). It is worth noting that the anion HSO4 can be ignored after phase equilibrium since the molality of SO42 is far greater than the molality of HSO4. Similarly, the primary (K1) and secondary (K2) ionization constants of tungstic acid obtained from the literature are 104.6 and 103.5, respectively. The concentrations of HWO4 are far less than that of WO42 under phase equilibrium; therefore, WO4 is considered the dominant ion. Additionally, the ions of H+ and OH− are minimal and can also be ignored when the pH value of the raffinate is close to 7.00. As discussed above, only the ions of Na+, SO42, and WO42 should be considered in the equilibrium aqueous phase. As for the organic phase, the complex H2WO4·RNH2, extractant (N1923), and diluent (toluene) need to be acknowledged. The mass balances and charge balance equations are as follows:



mNa+  2 mWO2



4

mH2 WO4 RNH2  C7 H8  mWO2





initial

,

H 2 SO 4  H + +HSO 4

(4)

mH2 WO4 RNH2  mRNH2   mRNH2 

Ks   H + +SO 42  HSO 4  

(5)

2 mSO2  2 mWO2  mNa+ .

K1   H + +HWO 4 H 2 WO 4  

(6)

K2   H + +WO 24  HWO 4 

(7)

4

2825.2 . T

(8)

4

4

initial

(11)



(12)

,

(13)

 mWO2  H2 O ,

(3)

ln K s =  14.0321+

(9)

Eq. (9) transforms into:

Na 2 WO 4  2Na + +WO24 

The secondary ionization thermodynamic equilibrium constant of sulfuric acid, Ks, can be calculated by the semi-empirical formula (eq. (8)).

.

4

4

initial

(14)

The density of water and toluene at 273.15 K in eqs. (2)–(12) is 0.9982 and 0.8670 kg/L respectively. 2.4.3 Pitzer model in the electrolyte aqueous solution In this work, the activity coefficients in the electrolyte aqueous solution were obtained from the Pitzer model [19–21]. The activity coefficient, or activity formula, of WO42, SO42, Na+, H+, and activity calculation formula of solvent H2O are expressed as:

938

Lin X, et al.



Sci China Tech Sci

May (2015) Vol.58 No.5

 I 2  ln 1  b I  1  b I b 



ln  WO  4 A  2 4



(0)  2mNa   WO +



 mWO mNa 2 4

2 + 4 ,Na





(1) 2  WO 2 ,Na + 4

 2I

(1) 4  WO 2 ,Na + 4

+

 mSO2 mNa + 4

 2I 2 (1) 4 SO ,Na 2 4

 2I 2

+

 1  1   I exp  I     



 





 



1  1   I   2 I 2 exp  I   



 



1  1   I   2 I 2 exp  I  ,  

(15) where A is the Debye-Hückel constant of water, I represents the ionic strength and b is Pitzer parameter which is a fixed constant (b=1.2).

ln aH2 O 



 mNa+  mSO2  mWO2 4

1000

4

M

,

H2 O

(16)

(17)

2.4.4 Nonelectrolyte organic phase model The activity coefficients of H2WO4·RNH2, N1923 (RNH2), and toluene (C7H8) are described by f1, f2, and f3 respectively. (I) Margules models [22] 2-Suffix formula ln f1  A12 x22  A13 x32 +(A12  A13  A23 ) x2 x3 ,

(18)

ln f 2  A x  A x +(A12  A23  A13 ) x1 x3 ,

(19)

ln f3  A x  A x +(A13  A23  A12 ) x1 x2 ,

(20)

2 13 1

2 23 3

2 23 2

3-Suffix formula

ln f1  x22  A12  2 x1  A21  A12    x32  A13  2 x1  A31  A13    x2 x3  A21  A13  A32  2 x1  A31  A13  2 x3  A32  A23   C (1  2 x1 )  ,

(21)

ln f 2  x32  A23  2 x2  A32  A23    x12  A21  2 x2  A12  A21    x1 x3  A32  A21  A12  2 x2  A12  A21  2 x1  A13  A31   C (1  2 x2 )  ,

(22)

ln f3  x12  A31  2 x3  A13  A31    x22  A32  2 x3  A23  A32    x1 x2  A13  A32  A21  2 x3  A23  A32  2 x2  A21  A12   C (1  2 x3 )  ,

m1 , m1  m2  m3

(25)

x3 

m3 . m1  m2  m3

(26)

(II) NRTL model [22] The activity coefficients of the organic complexes including H2WO4·RNH2, N1923(RNH2) and toluene (C7H8) are described as follows, and the i, j, k represent any compound among H2WO4·RNH2, RNH2 and C7H8. ln  i 



ji

G ji x j

j

G

x

 j

k

 l  li Gli xl Gij x j    k Gki xk  ij k Gkj xk

G ji  exp   ji ji  ,

(23)

(24)

  ,  

(27)

(28)

where Gji is the excess Gibbs free energy and x stands for mole fraction.

 ji   ij  0.2 .

(29)

The mole fractions of each substance in organic phases are expressed as:

xi 

mi .  mi

(30)

i

2.4.5 Methods of model solution The unknown parameters are solved by the nonlinear optimization method. GAMS (General Algebraic Modeling System) is used to program and accurately solve the models. This system is a computer numerical analysis commercial software, or a high-level modeling system designed for mathematical programming and optimization, solving linear, nonlinear, and mixed-integer optimization problems that are tailored for complex, large-scale modeling applications. Target function F is defined as eq. (31) during the modeling process. AARDWO2 is defined as the model error of eq. 4

(32), which is the average absolute relative deviation of WO4 in raffinate . cal. exp.  mWO 2 – m WO24 4 Target function: min F =  exp. mWO i 1  2  4 n

where fi represents the activity coefficient of species I and Aij is the terminal parameter in Margules model. The mole fractions of organic phases are expressed as:

x1 

m2 , m1  m2  m3

ki k

aH+  10 – pH .

2 12 1

x2 

2

  , 

cal. exp. 1 n mWO24 – mWO24 , AARD WO2 =  exp. 4 n i 1 mWO 2

(31)

(32)

4

where exp and cal correspond to the experimental data and calculated data respectively.

Lin X, et al.

3

Sci China Tech Sci

illustrate, the extraction percentages of W decrease by reducing pH value, which is consistent with previous literature research. The experiments are completed and represented by dots in the Figure 2, which are used to verify the prediction results from models. The lines in the figures are the results of the predictions. The experimental results agree well with the predicated results at low equilibrium pH range, but large deviation occurs at high equilibrium pH range. Since the model was established at low equilibrium pH range.

Results and discussion

3.1

939

May (2015) Vol.58 No.5

Liquid-Liquid phase equilibrium data

The data agrees with previous literature [14] and listed in Table 1. 3.2 Pitzer-Margules model of primary amine extraction of W

3.2.1 Results of model solution The parameters concluded by the 2-suffix Margules model are listed in Table 2. The parameters calculated by the 3-suffix Margules model are listed in Table 3. As for the table results, the AARD from the 2-suffix Margules model is more larger than the results from the 3-suffix Margules model. The binary interaction parameters, Aij from the 3-suffix Margules model, are double the values from the 2-suffix Margules model and up to 106 larger, which is not reasonable. The 2-suffix Margules model is taken into consideration in the model analyzing process while ignoring the 3-suffix Margules model.

3.3 Pitzer-NRTL model of primary amine extraction of W

3.3.1 Solution of model The parameters of the model calculated by nonlinear solver in GAMS are listed in Table 4. The calculated value from the Pitzer-NRTL model and experiments of AARDWO2 is 4

13.38%, which is the same as the results from the PitzerMargules (2-suffix) model. The calculated thermodynamic equilibrium constant is also consistent with the PitzerMargules (2-suffix) model.

3.2.2 Model predicated results New Pitzer-Margules models are used to predict the effect of equilibrium pH on the extraction of W with the fixed initial concentration of extract N1923 and W. As Figure 2

3.3.2 Model predicated results A line represents the predicated results in Figure 2, whereas the experiments results are represented as dots. The predicated

Table 1 Experimental data for the extraction of W by primary amine N1923 diluted in toluene at 293.15 K (phase ratio: O/A=1, equilibrium time: 1 h, standing time: 1 h)

m  WO24 

Note:

initial

(103 mol kg1)

m  org. RNH 2

(103 mol kg1)

initial

Initial pH

m  exp. WO24 

res.

(103 mol kg1)

Equilibrium pH

4.3513

5.8172

5.47

1.9113

6.73

4.3513

5.8172

5.47

1.8950

6.76

5.5324

5.8172

5.53

2.3282

6.70

5.5324

5.8172

5.53

2.3596

6.70

5.5324

5.8172

4.01

1.8065

6.19

5.5324

5.8172

4.01

1.7754

6.19

5.5324

12.2220

5.51

2.3173

6.93

5.5324

12.2220

5.51

2.3214

6.97

4.8254

12.2220

5.54

2.0953

6.88

4.8254

12.2220

5.54

2.1157

6.82

4.3513

12.2220

5.51

1.7969

6.93

4.3513

12.2220

5.51

1.8065

6.91

3.2737

12.2220

5.52

1.5204

6.84

3.2737

12.2220

5.52

1.4972

6.80

3.2737

12.2220

4.03

0.7793

6.68

3.2737

12.2220

4.03

0.7902

6.63

2.1307

12.2220

3.92

0.4738

6.54

2.1307

12.2220

3.92

0.4761

6.58

m  WO24 

initial

2 4

is the molality of WO

tration unit in mol/kg toluene.

in the residue solution. For aqueous phase, the concentration unit in mol/kg. For organic phase, the concen-

940

Lin X, et al.

Sci China Tech Sci

May (2015) Vol.58 No.5

Table 2 Pitzer-Margules (2-suffix) parameters for Na2WO4-H2SO4-primary amine N1923-toluene system at 293.15 K

i: H2WO4·RNH2 Aij 0 

 Na

+

,WO24 

1

 Na

+

,WO24 

j: RNH3

409.5658

i: RNH2; j: toluene 1.1108

i: H2WO4·RNH2

j: toluene

6.4714

15.0000 48.8769

ln(K)

28.3921

AARD

13.38 %

Table 3 Pitzer-Margules (3-suffix) parameters for Na2WO4-H2SO4-primary amine N1923-toluene system at 293.15 K

i: H2WO4·RNH2

j: RNH3

i: RNH2; j: toluene

i: toluene j: H2WO4·RNH2

Aij

20328.2060

1.1243

7905.3705

Aji

1948729.994

1557.2018

0.9682

C

1950916.467

0

 Na

+

,WO 24 

1

 Na

+

,WO24 

15.0000 48.9129

ln(K)

36.0340

AARD

12.86 %

Figure 2 (a) Effect of equilibrium pH on the extraction of W, (0.6398 g/L initial W, 1.5650 g/L initial N1923); (b) Effect of initial molality of N1923 on the extraction of W, (the 0.6738 g/L initial W and constant equilibrium pH 6.54). Temperature: 293.15 K, O/A=1, equilibrium time: 1 h, standing time: 1 h. Line M1, M2——calculated using the Pitzer-Margules (2-suffix) model, line N1, N2——calculated using the Pitzer-NRTL model, dot: experimental data.

Lin X, et al.

Sci China Tech Sci

941

May (2015) Vol.58 No.5

Table 4 Pitzer-NRTL parameters for Na2WO4-H2SO4-primary amine N1923 -toluene system at 293.15 K

i: RNH2; j: toluene

i: toluene j: H2WO4·RNH2

ij

i: H2WO4·RNH2 6.0671

133.0648

6.0767

ji

22.3217

1.1111

57.5248

0

 Na

+

15.000

,WO 24 

1

 Na

+

j: RNH3

48.8773

,WO 24 

ln(K)

28.7830

AARD

13.38%

results are fairly consistent with Pitzer-Margules results. 3.4

Summary and comparison of different model results

The two established models in this paper are compared to the models in previous literature to test accuracy. The equilibrium constant of each experimental group is recalculated by the model parameters and defined as lnKcal.. The equilibrium constant of the model regression is defined as lnKreg., and the error between them is expressed as average absolute relative deviation as shown in Table 5. A summary of the organic phase model parameters, equilibrium constants of different models, model error AARDWO2 , 4

1 lnK cal.  lnK reg. lnK reg. is listed in Table 5. The n differences between 1 AARDWO2 and lnK cal.  lnK reg. lnK reg. 4 n of the three models are negligible according to the results. There are still large differences between these models based on the numbers of parameters and physical meanings. For example, organic Pitzer model is a semi-empirical model with non-obvious physical meaning of Pitzer parameters  ; thus, it is difficult to carry deep theory analysis based on these data. The Pitzer- Margules (2-suffix) model is derived from the Scatchard-Hildebrand solution theory with obvious physical meaning of terminal parameter Aij, which can qualitatively explain the strength of interaction among various components in the organic phase and the relationship between molecule structures. The repulsive force between the different components is greater when Aij>0, which the attractive force is dominant when Aij