Prediction of Vapour Pressures of Dioxin Congeners. R & D Center, Baoshan Iron & Steel Co., Ltd., Fujin Rd., Baoshan District, Shanghai , China

Similar documents
dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential

δq T = nr ln(v B/V A )

Phase transition. Asaf Pe er Background

The Second Law: The Machinery

FUGACITY. It is simply a measure of molar Gibbs energy of a real gas.

Institute of Multidisciplinary Research for Advanced Materials (IMRAM), Tohoku University, Sendai , Japan

3. High Temperature Gases FPK1 2009/MZ 1/23

/ p) TA,. Returning to the

Chemistry 531 Spring 2009 Problem Set 6 Solutions

Phase Equilibrium Calculations by Equation of State v2

Liquid water static energy page 1/8

4. A Brief Review of Thermodynamics, Part 2

I affirm that I have never given nor received aid on this examination. I understand that cheating in the exam will result in a grade F for the class.

602 ZHANG Zhi and CHEN Li-Rong Vol Gibbs Free Energy From the thermodynamics, the Gibbs free energy of the system can be written as G = U ; TS +

General Physical Chemistry I

Phase transitions. Lectures in Physical Chemistry 4. Tamás Turányi Institute of Chemistry, ELTE. Phases

AT 25 C! CH10090 Thermodynamics (part 2) Enthalpy changes during reactions. Let s remember what we did in CH10089

F = U TS G = H TS. Availability, Irreversibility S K Mondal s Chapter 5. max. actual. univ. Ns =

Speed of sound measurements in liquid Methane at cryogenic temperature and for pressure up to 10 MPa

The Role of Water Vapor. atmosphere (we will ignore the solid phase here) Refer to the phase diagram in the web notes.

arxiv:cond-mat/ v2 25 Sep 2002

Chapter 2 Solutions. 2.1 (D) Pressure and temperature are dependent during phase change and independent when in a single phase.

Carrying out thermodynamic calculations and definition of the main reactions of decomposition of vapours of ethyl alcohol

The extreme case of the anisothermal calorimeter when there is no heat exchange is the adiabatic calorimeter.

Homework #11. (Due December 5 at the beginning of the class.) Numerical Method Series #6: Engineering applications of Newton-Raphson Method to

THERMODYNAMICS. Contents

Comparison of R744 and R134a heat transfer coefficients during flow boiling in a horizontal circular smooth tube

Mathematics as the Language of Physics.

Thermodynamics in combustion

1 Entropy 1. 3 Extensivity 4. 5 Convexity 5

ANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE

whether a process will be spontaneous, it is necessary to know the entropy change in both the

not to be republished NCERT STATES OF MATTER UNIT 5 INTRODUCTION

Internal Energy in terms of Properties

Chapter III: Statistical mechanics

Chemistry 420/523 Chemical Thermodynamics (Spring ) Examination 1

Standard partial molal properties of aqueous alkylphenols and alkylanilines over a wide range of temperatures and pressures

I have not proofread these notes; so please watch out for typos, anything misleading or just plain wrong.

Formation and transformation of volatile nanoparticles from a diesel engine during exhaust dilution

STATES OF MATTER UNIT 5

Chapter 9 Practical cycles

The temperature dependence of the isothermal bulk modulus at 1 bar pressure

Mixture Homogeneous Mixtures (air, sea water ) same composition, no chemical bond components are NOT distinguishable

Eric M. Yezdimer, and Robert H. Wood

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Evaporation

Casimir Force Between the Two Moving Conductive Plates.

Chapter 20: Exercises: 3, 7, 11, 22, 28, 34 EOC: 40, 43, 46, 58

Compressible Flow Introduction. Afshin J. Ghajar

The International Association for the Properties of Water and Steam

Calculators are permitted. Computers, PDAs, and other electronic devices with a keyboard are not. Cell phones may not be used as calculators.

Effect of hole-trapping on mass/charge transport properties in acceptor-doped BaTiO 3

PHYS1001 PHYSICS 1 REGULAR Module 2 Thermal Physics Chapter 17 First Law of Thermodynamics

ONE. The Earth-atmosphere system CHAPTER

First law of thermodynamics (Jan 12, 2016) page 1/7. Here are some comments on the material in Thompkins Chapter 1

Quasi-particle Contribution in Thermal Expansion and Thermal Conductivity in Metals

16. CHARACTERISTICS OF SHOCK-WAVE UNDER LORENTZ FORCE AND ENERGY EXCHANGE

First Name Last Name CIRCLE YOUR LECTURE BELOW: Div. 1 10:30 am Div. 2 2:30 pm Div. 3 4:30 pm Prof. Gore Prof. Udupa Prof. Chen

Modelling heat transfer and fluid flow inside a pressure cooker

ENTROPY GENERATION AND THERMAL PERFORMANCE OF A PULSATING HEAT PIPE. A Thesis. presented to. the Faculty of the Graduate School

Chemical Kinetics and Equilibrium - An Overview - Key

MODULE 2: DIFFUSION LECTURE NO DIFFUSION COEFFICIENT: MEASUREMENT AND PREDICTION

Chapter 1 Fundamentals

Thermodynamic Properties are Measurements p,t,v, u,h,s - measure directly -measure by change

Study of Current Decay Time during Disruption in JT-60U Tokamak

Entransy analysis of open thermodynamic systems

Equilibrium Thermodynamics

Lab 4 Module α 2. Miscibility Gaps

HEAT, WORK, AND THE FIRST LAW OF THERMODYNAMICS

Optical self-energy of superconducting Pb in the terahertz region

df da df = force on one side of da due to pressure

THERMODYNAMICS. Prepared by Sibaprasad Maity Asst. Prof. in Chemistry For any queries contact at

Minimum Fluidization Velocities of Binary- Solid Mixtures: Model Comparison

Pressure-sensitivity Effects on Toughness Measurements of Compact Tension Specimens for Strain-hardening Solids

What is Physical Chemistry?

Dioxins & PCBs concerns

Phase Diagrams. NC State University

Thermodynamic equations of state of polymers and conversion processing

Foundations of thermodynamics Fundamental thermodynamical concepts

The International Association for the Properties of Water and Steam

The Second Law of Thermodynamics (Chapter 4)

CHAPTER 2 THERMODYNAMICS

Thickness and refractive index measurements using multiple beam interference fringes (FECO)

Actual exergy intake to perform the same task

We now turn to considerations of mixtures. To keep our discussion reasonably simple,

Distribution of populations in excited states of electrodeless discharge lamp of Rb atoms

Evaluation of New Dual-Layer Carbon Reversible Column in Analysis of PCDDs/PCDFs and Co-Planar PCBs

On the q-deformed Thermodynamics and q-deformed Fermi Level in Intrinsic Semiconductor

Thermodynamic Modeling and Analysis of an Optical Electric-Field Sensor

A Simulation Model for the Application of Nanofluids as Condenser Coolants in Vapor Compression Heat Pumps

Chapter-6: Entropy. 1 Clausius Inequality. 2 Entropy - A Property

Determination of Pressure Losses in Hydraulic Pipeline Systems by Considering Temperature and Pressure

A numerical tool for plasma spraying. Part II: Model of statistic distribution of alumina multi particle powder.

Modeling Volume Changes in Porous Electrodes

NOTICE: This is the author's version of a work that was accepted for publication in Chemical Engineering Research and Design. Changes resulting from

The characteristic curves of water

Journal of Thermal Science and Technology

12.808: Some Physical Properties of Sea Water or, More than you ever wanted to know about the basic state variables of the ocean

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation

Transcription:

Prediction of Vaour Pressures of Dioxin Congeners Xian-Wei Li, 1 Etsuro Shibata, 2 Eiki Kasai, 2 Takashi Nakamura 2 1 R & D Center, Baoshan Iron & Steel Co., Ltd., Fujin Rd., Baoshan District, Shanghai 201900, China 2 Institute of Multidiscilinary Research for Advanced Material, Tohoku University, 1,1 Katahira, 2-Chome, Aobaku, Sendai 980-8577, Jaan Abstract: The vaour ressures and thermal roerties of non-measured olychlorinated dibenzo--dioxins and furans (PCDDs/PCDFs) are redicted by using a correlation method based on the exerimental vaor ressures obtained by the Knudsen effusion method. The vaor ressures of all the 59 PCDDs and 131 PCDFs redicted here are more or less higher than the data redicted by Rordorf, although the calculation method is the same. For the most toxic 2,3,7,8-TeCDD, the vaour ressure at 298 K redicted in this study is 6.2 10-6 Pa, 31 times higher than the value rovided by Rordorf. Keywords: PCDDs/PCDFs, Dioxins, Vaour ressure 1 Introduction Vaour ressures are fundamental roerties of ersistent organic ollutants and are imortant in determining their distributions and results in the environment. Due to large number of dioxin congeners and of other POPs, exerimental measurements can only cover a minor fraction of them, so that the rediction of vaour ressures of dioxin congeners by semi emirical correlations is useful and necessary. Prediction methods thus lay an imortant role in environmental research, roviding data otherwise inaccessible, and those more doubtful in testing. Numerous equations and correlations for estimating vaour ressure are resented in literature. Most methods are emirical and the derived results are very rough. [1-4] Rordorf [5-7] suggested a good correlation method to redict vaour ressure, which derives from thermodynamic theory. Based on exerimental vaour 1.

ressure data, this method gives not only the enthaly and entroy of sublimation but also the enthaly of melting and boiling oint. Based on our exerimental data reorted in the recently aer [8], this research emloys the correlation method to redict the vaour ressure of all PCDDs/PCDFs. Correlation methods are of articular imortance when dealing with dioxin congeners since there are thousands of ossible homologues and isomers. Correlation methods also allow testing of exerimental data set for self-consistency. Calculation of related substance roerties (i.e. boiling oints and enthalies of fusion) is ossible starting from vaour ressure data. 2 Prediction methods Equations that relate vaour ressure to temerature are commonly derived by integration of the Clausius-Claeyron equation: [1] d ln /dt= H v /( Z R T 2 ) (I) where is the vaour ressure, H v is the heat of vaorization, R is the gas constant, T is the temerature, and Z is a comressibility factor, given by Z = V/(R T) (II) where V is the volume difference between vaour and liquid. The correlation method uses the liquid hase as a reference state, as has been the case for most redictive vaour ressure methods. Suersaturated liquid hases are usually defined at the temerature of interest. A full use is made of the measured temerature deendence of the vaour ressure of the solid. The data are extraolated to the melting oint and the liquid hase is brought in at this oint. Enthalies ( H m ) and entroies of fusion ( S m ) are defined for the melting temerature (T m ) and unknown temerature deendencies of these functions are not involved. The two arameters are determined for the different comounds and are found to vary significantly between these comounds. The redicted method is not limited by assumtion of either temerature or substance indeendence of the entroies of fusion. Entroies of fusion are affected by symmetry relations between molecule and crystal lattice. Figure 1 illustrates the calculation rocedure of vaour ressure correlation method suggested by Rordorf. 2.

These correlations ermit the estimation of boiling oints and enthalies of evaoration for dioxins and the solid-state vaour ressures can be redicted by the equations of Fig 1 for dioxin congeners of known melting oints. ln (/Pa) T b (T b )=101,325 Pa Liquid T b = boiling oint G = molar Gibbs energy T m = melting oint T m liqud = solid H = molar enthaly S = molar entroy Subcooled liquid C = molar heat caacity R = gaas constant Solid (T) = vaour ressure at T Subscrit s : sublimation Subscrit v : vaorization K/T Subscrit m : fusion At the boiling oint: H v (T b ) = T b K F (36.61+R ln T b ) (1) where K F =1.01 S v (T b ) = H v (T b )/T b + R ln (2) where =101,325 Pa In the liquid hase: H ( T v 2 ) = H v ( T ) + C 1 T2 T1 gas liq. dt gives after integration for T 1 =T b and T 2 =T: H v (T) = H v (T b ) [1 + K ( 1 T/T b )] (4) with K = -T b C gas-liq. / H v (T b ) (5) Setting formula (4) into the Clausius-Claeyron equation, d ln ( T ) = T b T H v ( T b )[1 + K(1 T / T b (3) 2 )]/ RT dt gives after integration from the melting to boiling oint: ln (T m ) = ln - ( H v (T b )/R T b )[(1+K)(T b /T m -1)-K ln(t b /T m ) (7) At the melting oint: From G s = G m + G v (8) H m = H s - H v (9) and G = H - T S (10) follows S m = S s - S v (11) For T=T m : G m (T m )= 0 (12) and S m (T m ) = H m (T m )/ T m (13) From G v = - RT ln liq. = H v - T S v (14) follows for T = T m : S v (T m ) = R ln liq (T m ). + H v (T m )/T m (15) From equation (4), for T=T m : H v (T m ) = H v (T b ) [1 + K ( 1 T m /T b )] (4) In the solid hase: G s (T) = -RT ln solid (T) = H s (T) T S s (T) (16) R ln solid ( T ) = [ S s ( T m ) + T Tm C g s / TdT ] [ H R ln solid (T) = S s (T m ) - H s (T m )/T (18) s ( T m ) / T + Fig 1 The calculation rocedure of correlation method suggested by Rordorf [5]. (Relationshi between the vaour ressures of the solid and the liquid hase) (6) T Tm C g s dt ] (17) 3.

3 Predicted vaour ressure results for PCDDs/PCDFs Enthalies and entroies of sublimation of the 22 PCDDs/PCDFs were obtained by linear regression of the exerimental vaour ressures over the indicated temerature intervals in the aer ublished recently [8]. They hold for the mid temerature, T mid =2(T max T min )/(T max + T min ), of the investigated temerature ranges, and are recalculated for the melting oints and room temerature by means of bracket exressions of equation (17). The melting oint (T m ) data come from literature [9]. The boiling oint (T b ) aears in several exressions and is determined in an iterative rocedure. EXCEL is used for the solution of the equations shown in Fig 1 by an iterative rocedure. Figures 2 and 3 show that the enthaly of sublimation and derived boiling oint of PCDDs and PCDFs correlate well with the substitution number of chlorine, while the melting oint and entroy of sublimation show a great deal of scatter. 160 150 PCDDs 350 PCDDs 140 330 H sub(298k), kj/mol 130 120 110 S sub, J/mol/K 310 290 100 90 y = 6.5821x + 95.907 R 2 = 0.9174 270 y = 2.3454x + 297.57 R 2 = 0.1676 80 250 (a) (b) 160 350 150 PCDFs PCDFs 140 330 H sub(298k), kj/mol 130 120 110 S sub(298k), J/mol/K 310 290 100 90 y = 7.3232x + 88.561 R 2 = 0.9444 80 4. 270 y = 4.8032x + 281.83 R 2 = 0.6587 250

(c) (d) Fig 2 Chlorine correlations of thermal roerties, (a) enthaly of sublimation of PCDDs, (b) entroy of sublimation of PCDDs, (c) enthaly of sublimation of PCDFs, (d) entroy of sublimation of PCDFs 850 850 800 PCDDs 800 PCDFs 750 750 Tb (K) 700 650 Tb (K) 700 650 600 600 550 y = 29.105x + 564.001 R 2 = 0.949 550 y = 34.452x + 554.597 R 2 = 0.998 500 (a) 500 Cl substituion number (b) Fig 3 Chlorine correlations of boiling oint of PCDDs and PCDFs The determined boiling oints are then set for correct rediction of the vaour ressures ( liquid (T m )= solid (T m )) at the melting oints. The determined boiling oints and enthalies of fusion of the dioxin congeners are also correlated with the degree of chlorine substitution. Figures 3 and 4 show the correlations of the redicted boiling oints and enthalies of fusion with the degree of chlorination (x). The relation is as the following. For PCDDs: T b = 29.10 x + 564.00 ( K) H m (T m ) = 3.808x + 22.105 ( kj mol -1 ) For PCDFs: T b = 34.45 x + 554.60 (K) H m (T m ) = 3.553x + 20.801 (kj mol -1 ) Using these formulas, boiling oints and enthalies of fusion are redicted for unknown dioxin congeners, and then the solid-state vaour ressures are calculated by the equations of Fig 1. 5.

The thermal roerties and vaour ressures for 59 PCDDs and 131 PCDFs are redicted. 100 PCDDs 100 PCDFs 80 80 H (T m), kj/mol 60 40 H (T m), kj/mol 60 40 20 Hv(Tm) Hm(Tm) 0 20 Hv(Tm) Hm(Tm) 0 (a) (b) Fig 4 Chlorine correlations of enthalies of evaoration and fusion of PCDDs and PCDFs at melting oint All the vaour ressures of PCDDs/PCDFs redicted here are more or less higher than the data redicted by Rordorf [7] for the same comounds. Figure 5 shows the vaour ressure differences between Rordorf s study and this study for the toxic congeners, which are all substituted at the ositions 2, 3, 7, and 8. For the most toxic 2,3,7,8-TeCDD, the redicted vaour ressure at 25 by this study is 6.2 10-6 Pa, which is 31 times higher than the value given by Rordorf. Therefore, the concentration of air saturated with 2,3,7,8-TeCDD at 25 can arrive at 805 ng/m 3. 6.

0-2 -5-7 ln ( /Pa) -10 ln ( /Pa) -12-15 2,3,7,8-TeCDD,This study 2,3,7,8-TeCDD, Rordorf 1,2,3,7,8-PeCDD,This study 1,2,3,7,8-PeCDD, Rordorf 1,2,3,4,7,8-HxCDD,This study 1,2,3,4,7,8-HxCDD, Rordorf -20 2.4 2.6 2.8 3 3.2 3.4 1000 (K/T ) 1,2,3,7,8,9-HxCDD,This study -17 1,2,3,7,8,9-HxCDD, Rordorf 1,2,3,6,7,8-HxCDD,This study 1,2,3,6,7,8-HxCDD, Rordorf 1,2,3,4,6,7,8-HCDD,This study 1,2,3,4,6,7,8-HCDD, Rordorf -22 2.4 2.6 2.8 3 3.2 3.4 1000 (K/T ) (a) (b) 0 0-5 -5 ln ( /Pa) -10 ln ( /Pa) -10-15 2,3,7,8-TeCDF, This study 2,3,7,8-TeCDF, Rordorf 1,2,3,7,8-PeCDF, This study 1,2,3,7,8-PeCDF, Rordorf 1,2,3,4,7,8-HxCDF, This study 1,2,3,4,7,8-HxCDF, Rordorf -20 2.4 2.6 2.8 3 3.2 3.4-15 2,3,4,7,8-PeCDF, This study 2,3,4,7,8-PeCDF, Rordorf 2,3,4,6,7,8-HxCDF, This study 2,3,4,6,7,8-HxCDF, Rordorf 1,2,3,4,7,8,9-HCDF, This study 1,2,3,4,7,8,9-HCDF, Rordorf -20 2.4 2.6 2.8 3 3.2 3.4 1000 (K/T ) 1000 (K/T ) (c) (d) Fig 5 Comarison of vaour ressure of high chlorinated PCDDs/PCDFs redicted by Rordorf s study and this study 4 Conclusions Based on our measured vaour ressure data set of a number of dioxin congeners, the vaour ressure and thermal roerties of non-measured PCDDs/PCDFs are redicted by using a correlation method. The enthaly of 7.

sublimation and calculated boiling oint correlate well with the degree of chlorination, while the melting oint and entroy of sublimation show a great deal of scatter. All vaour ressures of the 59 PCDDs and 131 PCDFs redicted here are more or less higher than the data redicted by Rordorf, although the calculated method is the same. For the most toxic 2,3,7,8-TeCDD, the vaour ressure at 25 redicted by this study is 6.2 10-6 Pa, which is 31 times higher than the value given by Rordorf. The concentration of air saturated with 2,3,7,8-TeCDD at 25 can hence attain 805 ng/m 3. References 1 Lyman WJ, Reehl WF, Rosenblatt DH. Handbook of Chemical Proerty Estimation Methods, American Chemical Society, Washington, DC, 1990. 2 Miller DG. Estimating vaor ressures a comarison of equations. Ind. Eng. Chem., 56 (1964), 46. 3 Miller DG. Derivation of two equations for the estimation of vaor ressures, J. Phys. Chem., 68 (1964), 1399. 4 Partington JR. An Advanced Treatise on Physical Chemistry( Vol. 2). Longmans, Green and Co., London, 1951. 5 Rordorf BF. Thermal roerties of dioxins, furans and related comounds. Chemoshere, 15 (1986), 1325 1332. 6 Rordorf BF. Prediction of vaor ressures, boiling oints and enthalies of fusion for twenty-nine halogenated dibenzo--dioxines. Thermochim. Acta, 112 (1987), 117 122. 7 Rordorf BF. Prediction of vaour ressures, boiling oints and enthalies of fusion for twenty-nine halogenated dibenzo--dioxins and fifty-five dibenzofurans by a vaor ressure correlation method. Chemoshere, 18 8.

(1989), 783 788. 8 Li X-W, Shibata E, Kasai E, and Nakamura T.. Vaor ressures and enthalies of sublimation of 17 olychlorinated dibenzo--dioxins and 5 olychlorinated dibenzofurans. Environmental Toxicology and Chemistry,23 (2004), 348-354. 9 Shiu W-Y, Ma K-C. Temerature deendence of hysical-chemical roerties of selected chemicals of environmental interest. II. Chlorobenzenes, olychlorinated bihenyls, olychlorinated dibenzo--dioxins, and dibenzofurans. J. Phys. Chem. Ref. Data,29 (2000) 387-462. 9.