Enter Or Copy the Title of the Article Here

Total Page:16

File Type:pdf, Size:1020Kb

Enter Or Copy the Title of the Article Here This is the author’s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher’s web site or your institution’s library. Kirkwood–Buff integrals for ideal solutions Elizabeth A. Ploetz, Nikolaos Bentenitis, Paul E. Smith. How to cite this manuscript If you make reference to this version of the manuscript, use the following information: Ploetz, E.A., Bentenitis, N., & Smith, P.E. (2010). Kirkwood–Buff integrals for ideal solutions. Retrieved from http://krex.ksu.edu Published Version Information Citation: Ploetz, E.A., Bentenitis, N., & Smith, P.E. (2010). Kirkwood–Buff integrals for ideal solutions. The journal of chemical physics, 132(16), 9. Copyright: © 2010 American Institute of Physics. Digital Object Identifier (DOI): doi:10.1063/1.3398466 Publisher’s Link: http://jcp.aip.org/resource/1/jcpsa6/v132/i16/p164501_s1 This item was retrieved from the K-State Research Exchange (K-REx), the institutional repository of Kansas State University. K-REx is available at http://krex.ksu.edu THE JOURNAL OF CHEMICAL PHYSICS 132, 164501 ͑2010͒ Kirkwood–Buff integrals for ideal solutions ͒ Elizabeth A. Ploetz,1 Nikolaos Bentenitis,2 and Paul E. Smith1,a 1Department of Chemistry, Kansas State University, Manhattan, Kansas 66506, USA 2Department of Chemistry and Biochemistry, Southwestern University, Georgetown, Texas 78626, USA ͑Received 20 May 2009; accepted 29 March 2010; published online 22 April 2010͒ The Kirkwood–Buff ͑KB͒ theory of solutions is a rigorous theory of solution mixtures which relates the molecular distributions between the solution components to the thermodynamic properties of the mixture. Ideal solutions represent a useful reference for understanding the properties of real solutions. Here, we derive expressions for the KB integrals, the central components of KB theory, in ideal solutions of any number of components corresponding to the three main concentration scales. The results are illustrated by use of molecular dynamics simulations for two binary solutions mixtures, benzene with toluene, and methanethiol with dimethylsulfide, which closely approach ideal behavior, and a binary mixture of benzene and methanol which is nonideal. Simulations of a quaternary mixture containing benzene, toluene, methanethiol, and dimethylsulfide suggest this system displays ideal behavior and that ideal behavior is not limited to mixtures containing a small number of components. © 2010 American Institute of Physics. ͓doi:10.1063/1.3398466͔ I. INTRODUCTION solutions. In SI solutions the KBIs are neither zero nor inde- pendent of composition. The KBIs for these solutions can be The Kirkwood–Buff ͑KB͒ theory of solutions has pro- ͑␬ ͒ expressed in terms of the isothermal compressibility T and vided a wealth of data and insight into the properties of so- ͑ ͒ 1–13 molar volumes Vi of the pure components at the same T lution mixtures. KB theory can be applied to any stable and P. Expressions for binary, ternary, and quaternary solu- mixture containing any number of components. In particular, tions are available.2,20,21 Ben-Naim has also shown that SI the theory provides exact relationships from the molecular ⌬ solutions including up to four components satisfy Gij=Gii distributions between the various species in solution to the 22 +Gjj−2Gij=0 for all i, j pairs. Our recent analysis of the corresponding thermodynamics of the mixture. The central results obtained from the KB theory of solutions for n=1 to quantities of interest in KB theory are the KB integrals ͑ ͒ 1 4 components proposed the following general expression for KBIs , the KBIs:21 ϱ ␲͵ ͓ ␮VT͑ ͒ ͔ 2 ͑ ͒ n Gij = Gji =4 gij r −1 r dr, 1 SI ␬ ␳ 2 ͑ ͒ 0 Gij = RT T − Vi − Vj + Sn, Sn = ͚ kVk , 2 k=1 where g is the corresponding radial distribution function ij R ␳ N /V ͑rdf͒ between species i and j, and r is their intermolecular where is the gas constant and i = i are the number separation. In applying KB theory to understand solutions, it densities of each species. The above expression was postu- i j n is often useful to compare the properties of a real solution lated to be valid for any and combination in any com- with those of a corresponding ideal solution.2,6,14,15 In par- ponent SI solution. This can be shown to be true for n 21 ticular, deviations from ideal solution behavior provide indi- =1–4 components. Here we provide rigorous proof that cations of specific associations or affinities between the so- this expression is valid for any number of components at 16,17 T P lution components which affect the thermodynamics. constant and . Furthermore, we also determine the corre- Furthermore, assuming ideal behavior for some components sponding expressions for the KBIs which result in ideal be- may be the only possible approach to describe solutions con- havior on the molality and molarity concentration scales. The taining a large number of components, where experimental results are then illustrated using data obtained from computer data are usually very rare.18,19 Consequently, it is important simulations. to understand the nature of ideal solutions as formulated by KB theory. Ideality in solution mixtures at constant pressure ͑P͒ and II. THEORY ͑ ͒ temperature T can be expressed using a variety of concen- The chemical potential ͑␮͒ of a solute in solution can be tration scales. The most common is the mole fraction con- expressed in a variety of ways using a series of different centration scale where the corresponding activity coefficients concentration scales. There are three major expressions re- are unity for all species at all solution compositions. These lated to whether the species concentrations are represented in ͑ ͒ mixtures are also known as symmetric ideal SI or perfect terms of mole fractions ͑x͒, molalities ͑m͒, or molarities ͑c͒. It is then typical to write a͒ Author to whom correspondence should be addressed. Tel.: 785-532-5109. ␮ ␮o,x ␮o,x ␮ex FAX: 785-532-6666. Electronic mail: [email protected]. i = i + RT ln fixi = i + i + RT ln xi, 0021-9606/2010/132͑16͒/164501/9/$30.00132, 164501-1 © 2010 American Institute of Physics Downloaded 04 Aug 2011 to 129.130.37.167. Redistribution subject to AIP license or copyright; see http://jcp.aip.org/about/rights_and_permissions 164501-2 Ploetz, Bentenitis, and Smith J. Chem. Phys. 132, 164501 ͑2010͒ ␮ץ , ␮ = ␮o,m + RT ln ␥ m i i i i ␤ͩ i ͪ = ␮ , ij ץ ln mj T,P,m ␮ ␮o,c kj i = i + RT ln yici, ␮ ␮ץ i ij 3␳ ͑ ͒ ␤ͩ ͪ = , ͑7͒⌳ ء␮ ␮ ln ␳ ␦ − ␾ ץ i = i + RT ln i i, 3 j T,P,Nkj ij j ͑ ␥ ͒ where a set of activity coefficients f , ,y are used to quan- which are general for any number of components. Alterna- tify deviations from ideal behavior on the mole fraction, mo- tively, if one starts from Eq. ͑4͒ and then takes derivatives lality, and molarity concentration scales, respectively. The with respect to pressure with all Nj and T constant one finds final equation is the statistical mechanical expression for the that chemical potential in terms of the pseudochemical potential n thermal de Broglie wavelength ͑⌳͒, and number ,͒ء␮͑ ␬ ͑␦ ͒ ͑ ͒ density.2 The standard chemical potentials ͑␮o͒ can refer to RT T = ͚ ij + Nij Vj 8 j=1 an infinitely dilute solute or the pure solute at the same T and P. They are constants and independent of the concentration for any i component. of i, whereas the pseudochemical potential is composition The above equations relate the KBIs to properties of the dependent. The following discussion is focused on changes solution. It is easier, albeit less general, if we select a par- in the chemical potentials with composition, and therefore ticular species to continue. We choose i=1 and k1, and ͑ ͒ ͑ ͒ ͑ ␳ ͒ the choice of standard state is largely irrelevant. We also note then write Eqs. 5 and 8 after multiplication by 1 in an RT d ln y when T is constant. It is generally nϫn matrix form so that= ءthat d␮ known that the activity coefficients decrease when the solute ␳ V ␳ V ␳ V ␳ V 1+N ␳ RT␬ displays significant solute self-association, whereas the activ- 1 1 1 2 1 3 ¯ 1 n 11 1 T ␮ ␮ ␮ ␮ ␾ ity coefficients increase when the solute is significantly sol- 12 22 32 n2 N12 − 2 ␮ ␮ ␮ ␮ ␾ vated by the other species in solution. 13 23 33 n3 N13 = − 3 . The general approach used here to determine the re- ΄ ΅΄ ΅ ΄ ΅ ]]] ] quired KBIs involves generating a series of expressions in- ␮ ␮ ␮ ␮ N − ␾ volving derivatives of the chemical potentials which are 1n 2n 3n ¯ nn 1n n valid for real solutions, followed by application of the appro- ͑9͒ priate conditions corresponding to ideality for the various Hence, we have a set of simultaneous equations which can concentration scales. Let us consider the species number be solved quite easily for the required KBIs to give densities in the grand canonical ensemble to be functions of ␳ ␬ temperature and all the chemical potentials. One can then 1+N11 1RT T 23 ␾ write N12 − 2 n −1 ␾ ͑ ͒ N13 = Mn − 3 , 10 ␳ ␤ ͑␦ ͒ ␮ ͑ ͒ d ln i = ͚ ij + Nij d j, 4 ΄ ΅ ΄ ΅ ] ] j=1 ␾ N1n − n which is valid for changes in the number density of any where M is the matrix from Eq. ͑9͒. This is the approach we component in any multicomponent system and any ͑thermo- n used recently for outlining a general KB inversion dynamically reasonable͒ ensemble with T constant. Here, procedure.21 After the appropriate chemical potential deriva- N =␳ G , ␤=1/RT, and ␦ is the Kroenecker delta function.
Recommended publications
  • Development Team
    Paper No: 16 Environmental Chemistry Module: 01 Environmental Concentration Units Development Team Prof. R.K. Kohli Principal Investigator & Prof. V.K. Garg & Prof. Ashok Dhawan Co- Principal Investigator Central University of Punjab, Bathinda Prof. K.S. Gupta Paper Coordinator University of Rajasthan, Jaipur Prof. K.S. Gupta Content Writer University of Rajasthan, Jaipur Content Reviewer Dr. V.K. Garg Central University of Punjab, Bathinda Anchor Institute Central University of Punjab 1 Environmental Chemistry Environmental Environmental Concentration Units Sciences Description of Module Subject Name Environmental Sciences Paper Name Environmental Chemistry Module Name/Title Environmental Concentration Units Module Id EVS/EC-XVI/01 Pre-requisites A basic knowledge of concentration units 1. To define exponents, prefixes and symbols based on SI units 2. To define molarity and molality 3. To define number density and mixing ratio 4. To define parts –per notation by volume Objectives 5. To define parts-per notation by mass by mass 6. To define mass by volume unit for trace gases in air 7. To define mass by volume unit for aqueous media 8. To convert one unit into another Keywords Environmental concentrations, parts- per notations, ppm, ppb, ppt, partial pressure 2 Environmental Chemistry Environmental Environmental Concentration Units Sciences Module 1: Environmental Concentration Units Contents 1. Introduction 2. Exponents 3. Environmental Concentration Units 4. Molarity, mol/L 5. Molality, mol/kg 6. Number Density (n) 7. Mixing Ratio 8. Parts-Per Notation by Volume 9. ppmv, ppbv and pptv 10. Parts-Per Notation by Mass by Mass. 11. Mass by Volume Unit for Trace Gases in Air: Microgram per Cubic Meter, µg/m3 12.
    [Show full text]
  • Supplement Of
    Supplement of Effects of Liquid–Liquid Phase Separation and Relative Humidity on the Heterogeneous OH Oxidation of Inorganic-Organic Aerosols: Insights from Methylglutaric Acid/Ammonium Sulfate Particles Hoi Ki Lam1, Rongshuang Xu1, Jack Choczynski2, James F. Davies2, Dongwan Ham3, Mijung Song3, Andreas Zuend4, Wentao Li5, Ying-Lung Steve Tse5, Man Nin Chan 1,6 1Earth System Science Programme, Faculty of Science, The Chinese University of Hong Kong, Hong Kong, China 2Department of Chemistry, University of California Riverside, Riverside, CA, USA 3Department of Earth and Environmental Sciences, Jeonbuk National University, Jeollabuk-do, Republic of Korea 4Department of Atmospheric and Oceanic Sciences, McGill University, Montreal, Québec, Canada 5Departemnt of Chemistry, The Chinese University of Hong Kong, Hong Kong, China’ 6The Institute of Environment, Energy and Sustainability, The Chinese University of Hong Kong, Hong Kong, China Corresponding author: [email protected] Table S1. Composition, viscosity, diffusion coefficient and mixing time scale of aqueous droplets containing 3-MGA and ammonium sulfate (AS) in an organic- to-inorganic dry mass ratio (OIR) = 1 at different RH predicted by the AIOMFAC-LLE. RH (%) 55 60 65 70 75 80 85 88 Salt-rich phase (Phase α) Mass fraction of 3-MGA 0.00132 0.00372 0.00904 0.0202 / / Mass fraction of AS 0.658 0.616 0.569 0.514 / / Mass fraction of H2O 0.341 0.380 0.422 0.466 / / Organic-rich phase (Phase β) Mass fraction of 3-MGA 0.612 0.596 0.574 0.546 / / Mass fraction of AS 0.218 0.209 0.199 0.190
    [Show full text]
  • Colloidal Suspensions
    Chapter 9 Colloidal suspensions 9.1 Introduction So far we have discussed the motion of one single Brownian particle in a surrounding fluid and eventually in an extaernal potential. There are many practical applications of colloidal suspensions where several interacting Brownian particles are dissolved in a fluid. Colloid science has a long history startying with the observations by Robert Brown in 1828. The colloidal state was identified by Thomas Graham in 1861. In the first decade of last century studies of colloids played a central role in the development of statistical physics. The experiments of Perrin 1910, combined with Einstein's theory of Brownian motion from 1905, not only provided a determination of Avogadro's number but also laid to rest remaining doubts about the molecular composition of matter. An important event in the development of a quantitative description of colloidal systems was the derivation of effective pair potentials of charged colloidal particles. Much subsequent work, largely in the domain of chemistry, dealt with the stability of charged colloids and their aggregation under the influence of van der Waals attractions when the Coulombic repulsion is screened strongly by the addition of electrolyte. Synthetic colloidal spheres were first made in the 1940's. In the last twenty years the availability of several such reasonably well characterised "model" colloidal systems has attracted physicists to the field once more. The study, both theoretical and experimental, of the structure and dynamics of colloidal suspensions is now a vigorous and growing subject which spans chemistry, chemical engineering and physics. A colloidal dispersion is a heterogeneous system in which particles of solid or droplets of liquid are dispersed in a liquid medium.
    [Show full text]
  • Calculation Exercises with Answers and Solutions
    Atmospheric Chemistry and Physics Calculation Exercises Contents Exercise A, chapter 1 - 3 in Jacob …………………………………………………… 2 Exercise B, chapter 4, 6 in Jacob …………………………………………………… 6 Exercise C, chapter 7, 8 in Jacob, OH on aerosols and booklet by Heintzenberg … 11 Exercise D, chapter 9 in Jacob………………………………………………………. 16 Exercise E, chapter 10 in Jacob……………………………………………………… 20 Exercise F, chapter 11 - 13 in Jacob………………………………………………… 24 Answers and solutions …………………………………………………………………. 29 Note that approximately 40% of the written exam deals with calculations. The remainder is about understanding of the theory. Exercises marked with an asterisk (*) are for the most interested students. These exercises are more comprehensive and/or difficult than questions appearing in the written exam. 1 Atmospheric Chemistry and Physics – Exercise A, chap. 1 – 3 Recommended activity before exercise: Try to solve 1:1 – 1:5, 2:1 – 2:2 and 3:1 – 3:2. Summary: Concentration Example Advantage Number density No. molecules/m3, Useful for calculations of reaction kmol/m3 rates in the gas phase Partial pressure Useful measure on the amount of a substance that easily can be converted to mixing ratio Mixing ratio ppmv can mean e.g. Concentration relative to the mole/mole or partial concentration of air molecules. Very pressure/total pressure useful because air is compressible. Ideal gas law: PV = nRT Molar mass: M = m/n Density: ρ = m/V = PM/RT; (from the two equations above) Mixing ratio (vol): Cx = nx/na = Px/Pa ≠ mx/ma Number density: Cvol = nNav/V 26 -1 Avogadro’s
    [Show full text]
  • Alkali Metal Vapor Pressures & Number Densities for Hybrid Spin Exchange Optical Pumping
    Alkali Metal Vapor Pressures & Number Densities for Hybrid Spin Exchange Optical Pumping Jaideep Singh, Peter A. M. Dolph, & William A. Tobias University of Virginia Version 1.95 April 23, 2008 Abstract Vapor pressure curves and number density formulas for the alkali metals are listed and compared from the 1995 CRC, Nesmeyanov, and Killian. Formulas to obtain the temperature, the dimer to monomer density ratio, and the pure vapor ratio given an alkali density are derived. Considerations and formulas for making a prescribed hybrid vapor ratio of alkali to Rb at a prescribed alkali density are presented. Contents 1 Vapor Pressure Curves 2 1.1TheClausius-ClapeyronEquation................................. 2 1.2NumberDensityFormulas...................................... 2 1.3Comparisonwithotherstandardformulas............................. 3 1.4AlkaliDimers............................................. 3 2 Creating Hybrid Mixes 11 2.1Predictingthehybridvaporratio.................................. 11 2.2Findingthedesiredmolefraction.................................. 11 2.3GloveboxMethod........................................... 12 2.4ReactionMethod........................................... 14 1 1 Vapor Pressure Curves 1.1 The Clausius-Clapeyron Equation The saturated vapor pressure above a liquid (solid) is described by the Clausius-Clapeyron equation. It is a consequence of the equality between the chemical potentials of the vapor and liquid (solid). The derivation can be found in any undergraduate text on thermodynamics (e.g. Kittel & Kroemer [1]): Δv · ∂P = L · ∂T/T (1) where P is the pressure, T is the temperature, L is the latent heat of vaporization (sublimation) per particle, and Δv is given by: Vv Vl(s) Δv = vv − vl(s) = − (2) Nv Nl(s) where V is the volume occupied by the particles, N is the number of particles, and the subscripts v & l(s) refer to the vapor & liquid (solid) respectively.
    [Show full text]
  • Receive! Osti
    DOE/MC/30175-5033 (DE96000569) RECEIVE! NOV 2 11995 OSTI Portable Sensor for Hazardous Waste Topical Report October 1993 - September 1994 Dr. Lawrence G. Piper October 1994 Work Performed Under Contract No.: DE-AC21-93MC30175 For U.S. Department of Energy U.S. Department of Energy Office of Environmental Management Office of Fossil Energy Office of Technology Development Morgantown Energy Technology Center Washington, DC Morgantown, West Virginia By Physical Sciences Inc. Andover, Massachusetts IAS DISTRIBUTION DP TMi.Q nnpiiycwr » I !MI IS UTT-n, S\c^ DISCLAIMER This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manu• facturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. This report has been reproduced directly from the best available copy. Available to DOE and DOE contractors from the Office of Scientific and Technical Information, 175 Oak Ridge Turnpike, Oak Ridge, TN 37831; prices available at (615) 576-8401. Available to the public from the National Technical Information Service, U.S.
    [Show full text]
  • Density Lab B Test
    Northern Regional: January 19th, 2019 Density Lab B Test Name(s): _______________________________________________________________________ Team Name: ________________________________________ School Name: ________________________________________ Rank: ___________ Team Number: ___________ Score: ___________ Written Test 1) All of the below are types of densities except (5 points) a) Mass Density b) Number Density c) Area Density d) Volume Density 2) Which of the following exhibits Boyle’s Law? (5 points) a) P1V1= P2V2 b) P1/T1=P2/T2 c) PV= nRT d) V1/T1=V2/T2 3) How many atoms are in 3 moles of a substance? (5 points) 23 a) 24.023*10 atoms​ ​ 24 b) 1.81*10 ​ atoms ​ c) 23 atoms 23 d) 6.02310 atoms​ ​ 4) If the temperature of a container doubled, what would have happened to the pressure of the gas inside the container? (5 points) a) Halved b) Doubled c) Quadrupled d) Nothing 5) What factors affect density? (5 points) a) Pressure b) Temperature c) Both A and B d) None of the above 6) Water is more dense as a solid. (5 points) a) True b) False 7) Pascal is the SI Unit for pressure. (5 points) a) True b) False 8) The acceleration due to gravity on Earth is 9.8 m/s. (5 points) a) True b) False 9) Salt water is more dense than deionized water. (5 points) a) True b) False 10) By Archimedes' principle, the buoyant force is equal to the weight of fluid displaced. (5 points) a) True b) False 11) What is the Archimedes Principle? (5 points) 12) What is the molarity of a solution containing 18.9 grams of RuCl3 in enough water to ​ ​ make 1.00 L of solution?
    [Show full text]
  • Stellar Structure III
    7 – Stellar Structure III introduc)on to Astrophysics, C. Bertulani, Texas A&M-Commerce 1 Fundamental physical constants a radiation density constant 7.55 × 10-16 J m-3 K-4 c velocity of light 3.00 × 108 m s-1 G gravitational constant 6.67 × 10-11 N m2 kg-2 h Planck’s constant 6.62 × 10-34 J s K Boltzmann’s constant 1.38 × 10-23 J K-1 -31 me mass of electron 9.11 × 10 kg -27 mH mass of hydrogen atom 1.67 × 10 kg 23 -1 NA Avogadro’s number 6.02 × 10 mol σ Stefan Boltzmann constant 5.67 × 10-8 W m-2 K-4 (σ = ac/4) R gas constant (k/mH) 8.26 × 103 J K-1 kg-1 e charge of electron 1.60 × 10-19 C 26 L¤ luminosity of Sun 3.86 × 10 W 30 M¤ mass of Sun 1.99 × 10 kg Teff¤ effective temperature of sun 5780 K 8 R¤ radius of Sun 6.96 × 10 m Parsec (unit of distance) 3.09 × 1016 m introduc)on to Astrophysics, C. Bertulani, Texas A&M-Commerce 2 22 7.1 - The equation of radiative transport We assume for the moment that the condition for convection is not satisfied, and we will derive an expression relating the change in temperature with radius in a star assuming all energy is transported by radiation. Hence we ignore the effects of convection and conduction. The equation of radiative transport with gas conditions is a function of only one coordinate, in this case r.
    [Show full text]
  • The Evolution of Galaxy Number Density
    Draft version October 9, 2016 A Preprint typeset using L TEX style AASTeX6 v. 1.0 THE EVOLUTION OF GALAXY NUMBER DENSITY AT Z < 8 AND ITS IMPLICATIONS Christopher J. Conselice, Aaron Wilkinson, Kenneth Duncan1, Alice Mortlock2 University of Nottingham, School of Physics & Astronomy, Nottingham, NG7 2RD UK 1Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, the Netherlands 2SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh, EH9 3HJ ABSTRACT The evolution of the number density of galaxies in the universe, and thus also the total number of galaxies, is a fundamental question with implications for a host of astrophysical problems including galaxy evolution and cosmology. However there has never been a detailed study of this important measurement, nor a clear path to answer it. To address this we use observed galaxy stellar mass functions up to z ∼ 8 to determine how the number densities of galaxies changes as a function of time and mass limit. We show that the increase in the total number density of galaxies (φT), more massive 6 −1 than M∗ = 10 M⊙ , decreases as φT ∼ t , where t is the age of the universe. We further show 7 that this evolution turns-over and rather increases with time at higher mass lower limits of M∗ > 10 6 M⊙ . By using the M∗ = 10 M⊙ lower limit we further show that the total number of galaxies in the +0.7 12 universe up to z = 8 is 2.0−0.6 × 10 (two trillion), almost a factor of ten higher than would be seen in an all sky survey at Hubble Ultra-Deep Field depth.
    [Show full text]
  • Measurement of Particle Number Density and Volume Fraction in A
    Measurement of Particle Number Density and Volume Fraction in a Fluidized Bed using Shadow Sizing Method Seckin Gokaltun, Ph.D David Roelant, Ph.D Applied Research Center FLORIDA INTERNATIONAL UNIVERSITY Objective y 2006 Roadmap Task D1 & D2: y Detailed CFB data at ~.15m ID y Non-intrusive probes y Generate detailed experimental data to develop a new mathematical analysis procedure for defining cluster formation by utilizing granular temperature. y Measure void fraction y Obtain granular temperature y Demonstrate nonintrusive probe to collect void fraction data Granular Temperature Granular theory Kinetic theory y From the ideal gas law: y Assumptions: y Very large number of PV = NkBT particles (for valid statistical treatment) where kB is the Boltzmann y Distance among particles much larger than molecular constant and the absolute size temperature, thus: y Random particle motion with 2 constant speeds PV = NkBT = Nmvrms / 3 2 y Elastic particle-particle and T = mvrms / 3kB particle-walls collisions (no n loss of energy) v 2 = v 2 y Molecules obey Newton Laws rms ∑ i i =1 is the root-mean-square velocity Granular Temperature y From velocity 1 θ = ()σ 2 + σ 2 + σ 2 3 θ r z 2 2 where: σ z = ()u z − u y uz is particle velocity u is the average particle velocity y From voidage for dilute flows 2 3 θ ∝ ε s y From voidage for dense flows 2 θ ∝1 ε s 6 Glass beads imaged at FIU 2008 Volume Fraction y Solids volume fraction: nV ε = p s Ah y n : the number of particles y Vp : is the volume of single particle y A : the view area y h : the depth of view (Lecuona et al., Meas.
    [Show full text]
  • 1 Atmospheric Sciences 501. Course Notes Autumn 1998. 1. Ideal Gas
    1 Atmospheric Sciences 501. Course Notes Autumn 1998. 1. Ideal Gas Law and Hydrostatics. 1.1 Ideal Gas Law. Basic thermodynamic variables are: Pressure, p (N m-2; 1 N m-2 = 1 pascal (Pa) = 0.01 millibar (mb) = 0.01hectoPascals (hPa). Density, r (kgm-3) Temperature, T (K) . If the energy due to intermolecular forces is small relative to molecular kinetic energy, the ideal gas law is an excellent approximation to the equation of state relating these variables: pRT=r (1.1) R* where R = = gas constant; R* = universal gas constant = 8314J kmol-1 K-1; M = molar M weight (kg/kg-mole) and R = specific gas constant for dry air = 287 J kg-1 K-1. 1.2 Molecules and the ideal gas law, Maxwell's kinetic interpretation. Note that r=nmÃà (1.2) where nà = number density (number of molecules per unit volume m-3), and mà = molecular mass (kg). Also MmA= à (1.3) where A = Avogadro's Number = 6.022x1026 molecules per mole; note that, since SI units are used, moles here are kilogram-moles, not gram-moles. Substituting (1.2) and (1.3) into (1.1): 2 æ R* ö p = ç ÷ nTÃÃ= knT (1.4) è A ø * æ R ö -23 -1 where k = ç ÷ = 1.381x10 J K is Boltzmann's Constant. è A ø Now consider the pressure exerted by molecules of gas on a surface perpendicular to the x- axis in a Cartesian reference frame. This is the force per unit area and is equivalent to the average flux in the x-direction of molecules with x-momentum mvà , mvÃà nv x ()xx× where () represents averaging over molecular velocities.
    [Show full text]
  • A Theory for Calculating the Number Density Distribution of Small Particles on a Flat Wall from Pressure Between the Two Walls
    A theory for calculating the number density distribution of small particles on a flat wall from pressure between the two walls Kota Hashimotoa and Ken-ichi Amanoa aDepartment of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan. E-mail: [email protected] ABSTRACT Surface force apparatus (SFA) and atomic force microscopy (AFM) can measure a force curve between a substrate and a probe in liquid. However, the force curve had not been transformed to the number density distribution of solvent molecules (colloidal particles) on the substance due to the absence of such a transform theory. Recently, we proposed and developed the transform theories for SFA and AFM. In these theories, the force curve is transformed to the pressure between two flat walls. Next, the pressure is transformed to number density distribution of solvent molecules (colloidal particles). However, pair potential between the solvent molecule (colloidal particle) and the wall is needed as the input of the calculation and Kirkwood superposition approximation is used in the previous theories. In this letter, we propose a new theory that does not use both the pair potential and the approximation. Instead, it makes use of a structure factor between solvent molecules (colloidal particles) which can be obtained by X-ray or neutron scattering. MAIN TEXT Structure of solid-liquid interface has relationship with physical properties and reactions at interfaces. Measurement of the density distribution of solvent molecules at the interface is important to understand them. Furthermore, the density distribution of colloid particles on solids or membranes is important for study of soft matter.
    [Show full text]