Lecture 16. Kinetics and Mass Transfer in Crystallization

Similar documents
PID CONTROL. Presentation kindly provided by Dr Andy Clegg. Advanced Control Technology Consortium (ACTC)

Chapter 6 Control Systems Design by Root-Locus Method. Lag-Lead Compensation. Lag lead Compensation Techniques Based on the Root-Locus Approach.

Period #8: Axial Load/Deformation in Indeterminate Members

Notes on Implementation of Colloid Transport and Colloid-Facilitated Solute Transport into HYDRUS-1D

ES 247 Fracture Mechanics Zhigang Suo. Applications of Fracture Mechanics

DISCHARGE MEASUREMENT IN TRAPEZOIDAL LINED CANALS UTILIZING HORIZONTAL AND VERTICAL TRANSITIONS

Heat transfer and absorption of SO 2 of wet flue gas in a tube cooled

Chapter 4. Simulations. 4.1 Introduction

Thermochemistry and Calorimetry

Inverse Kinematics 1 1/21/2018

TENSILE STRENGTH MODELING OF GLASS FIBER-POLYMER COMPOSITES AND SANDWICH MATERIALS IN FIRE

KINETICS OF IRON OXIDE DIRECT REDUCTION BY COAL E.R. ABRIL 1

@(; t) p(;,b t) +; t), (; t)) (( whih lat line follow from denition partial derivative. in relation quoted in leture. Th derive wave equation for ound

Bio-Systems Modeling and Control

The Model Predictive Control System for the Fluid Catalytic Cracking Unit

Deepak Rajput

UVa Course on Physics of Particle Accelerators

Compensation Techniques

IMPACT OF PRESSURE EQUALIZATION SLOT IN FLOW CHANNEL INSERT ON TRITIUM TRANSPORT IN A DCLL-TYPE POLOIDAL DUCT. H. Zhang, A. Ying, M.

Where Standard Physics Runs into Infinite Challenges, Atomism Predicts Exact Limits

Lecture 4. Kinetics and Transport in Sorption (1)

THE SOLAR SYSTEM. We begin with an inertial system and locate the planet and the sun with respect to it. Then. F m. Then

Lecture 24: Spinodal Decomposition: Part 3: kinetics of the

crown/cap crest α b core/base bedding and/or filter

Supplementary Materials for

Department of Pesticide Regulation

MODULE 3: MASS TRANSFER COEFFICIENTS

Establishment of Model of Damping Mechanism for the Hard-coating Cantilever Plate

Analysis of Feedback Control Systems

Perfusion-enhanced growth of. tissue-engineered cartilage in a. bioreactor: A finite element study

Line Radiative Transfer

Velocity distributions and correlations in homogeneously heated granular media

MIKE 21 & MIKE 3 Flow Model FM. Mud Transport Module. Scientific Documentation

Homework Set 4. gas B open end

Transition from phase slips to the Josephson effect in a superfluid 4 He weak link

IV Transport Phenomena. Lecture 21: Solids and Concentrated Solutions

F O R SOCI AL WORK RESE ARCH

The elementary spike produced by a pure e + e pair-electromagnetic pulse from a Black Hole: The PEM Pulse

COMM 602: Digital Signal Processing. Lecture 8. Digital Filter Design

5.2.6 COMPARISON OF QUALITY CONTROL AND VERIFICATION TESTS

INDIVIDUAL OVERTOPPING EVENTS AT DIKES

Nonlinear Dynamics of Single Bunch Instability in Accelerators *

Torsional resistance of high-strength concrete beams

The elementary spike produced by a pure e + e pair-electromagnetic pulse from a Black Hole: The PEM Pulse

Steady-state response of systems with fractional dampers

Determination of the reaction order

Cogging torque reduction of Interior Permanent Magnet Synchronous Motor (IPMSM)

NUMERICAL SIMULATION ON THE FIREPROOF BEHAVIOR OF RC BEAM STRENGTHENED WITH STRANDED MESH AND POLYMER MORTAR

Design Manual to EC2. LinkStudPSR. Version 3.1 BS EN : Specialists in Punching Shear Reinforcement.

Critical Percolation Probabilities for the Next-Nearest-Neighboring Site Problems on Sierpinski Carpets

Wave Phenomena Physics 15c

To determine the biasing conditions needed to obtain a specific gain each stage must be considered.

Duct Acoustics. Chap.4 Duct Acoustics. Plane wave

A consistent beam element formulation considering shear lag effect

ANALYSIS OF A REDUNDANT SYSTEM WITH COMMON CAUSE FAILURES

Modeling and Simulation of Buck-Boost Converter with Voltage Feedback Control

Feedback Control Systems (FCS)

Convective Heat Transfer

DETERMINATION OF THE POWER SPECTRAL DENSITY IN CAPACITIVE DIGITAL ACCELEROMETERS USING THEORY OF LIMIT CYCLES

Developing Transfer Functions from Heat & Material Balances

Mass Transfer (Stoffaustausch) Fall 2012

NONISOTHERMAL OPERATION OF IDEAL REACTORS Plug Flow Reactor

PREDICTION OF PRESSURE DROP IN HELICAL COIL WITH SINGLE PHASE FLOW OF NON-NEWTONIAN FLUID

Mass Transfer 2. Diffusion in Dilute Solutions

MOLECULAR ORBITAL THEORY- PART I

Study of Upper Critical Magnetic Field of Superconducting HoMo6Se8

Chemistry (Physical chemistry) Lecture 10.

EE 333 Electricity and Magnetism, Fall 2009 Homework #11 solution

Subject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (

1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #11: Solute Transport in Groundwater

Today s Lecture. Block Diagrams. Block Diagrams: Examples. Block Diagrams: Examples. Closed Loop System 06/03/2017

CEE 670 TRANSPORT PROCESSES IN ENVIRONMENTAL AND WATER RESOURCES ENGINEERING. Kinetics Lecture #1

A Queueing Model for Call Blending in Call Centers

Sound Propagation through Circular Ducts with Spiral Element Inside

Energy-Work Connection Integration Scheme for Nonholonomic Hamiltonian Systems

A Simplified Steady-State Analysis of the PWM Zeta Converter

Einstein's Energy Formula Must Be Revised

12th International Congress on the Deterioration and Conservation of Stone Columbia University, New York, 2012

OLIGONUCLEOTIDE microarrays are widely used

On the Stationary Convection of Thermohaline Problems of Veronis and Stern Types

Flexural Reinforced Concrete Elements Residual Stresses of Fibres with Limited Reinforcement Bond

New Directions in Interconnect Performance Optimization

SOLUTIONS

Control Engineering An introduction with the use of Matlab

Activation energies of the In Si -Si i defect transitions obtained by carrier lifetime measurements

Econ 455 Answers - Problem Set 4. where. ch ch ch ch ch ch ( ) ( ) us us ch ch us ch. (world price). Combining the above two equations implies: 40P

General Equilibrium. What happens to cause a reaction to come to equilibrium?

Creep and Shrinkage Analysis of Curved Composite Beams Including the Effects of Partial Interaction

Lecture 24. Ideal Gas Law and Kinetic Theory

ANALYSIS OF FLEXIBLE CLAMPING IN TENSILE TESTS OF MULTIDIRECTIONAL LAMINATES

Sidelobe-Suppression Technique Applied To Binary Phase Barker Codes

Friction-induced Vibration in Lead Screw Systems

CART: Classification and Regression Trees

IMPACT PERFORATION OF COMPOSITE SANDWICH PANELS

The longevity of water ice on Ganymedes and Europas around migrated giant planets Owen R. Lehmer 1, David C. Catling 1, Kevin J.

A generic method for modeling the behavior of. anisotropic metallic materials : application to. recrystallized zirconium alloys

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

Chapter 14. The Concept of Equilibrium and the Equilibrium Constant. We have for the most part depicted reactions as going one way.

Advanced Computational Fluid Dynamics AA215A Lecture 4

REVIEW QUESTIONS Chapter 15

Transcription:

Leture 16. Kineti and Ma Tranfer in Crytallization Crytallization Kineti Superaturation Nuleation - Primary nuleation - Seondary nuleation Crytal Growth - Diffuion-reation theory - Srew-diloation theory

Crytallization Kineti Crytallization i a omplex phenomenon involving three tep: (1) nuleation, (2) ma tranfer of olute to the rytal urfae, and (3) inorporation of olute into the rytal lattie Þ Crytallization kineti Superaturation i the driving fore for rytallization kineti A rytal ize dereae, olubility notieably inreae, making it poible to uperaturate a olution if it i ooled lowly without agitation The olubility of very mall rytal an fall in the metatable region Point a: the olution i underaturated; rytal of all ize diolve Point b: equilibrium between a aturated olution and rytal that an be een by naked eye

Superaturation Point : metatable region; rytal an grow but annot nuleate Point d: pontaneou nuleation of very mall rytal, that are inviible to the naked eye, our Limiting uperaturation, D limit m Superaturation, D Limit of uperaturation temperature differene Superaturation temperature differene Kelvin equation: a relationhip between olubility and rytal ize æ ö ln ç è ø 4v, L vrtd Relative uperaturation - p v : molar volume of the rytal,l : interfaial tenion v : number of ion/moleule of olute / : uperaturation ratio (S) In pratie, i uually le than 2% - 1 S - 1

Nuleation (1) To determine volume or reidene time for magma in a rytallizer, the rate of nuleation (birth) of rytal and their rate of growth mut be etablihed Relative rate of nuleation and growth are important beaue they determine rytal ize and ize ditribution Primary nuleation - Superaturated olution i free of rytalline urfae - Requre high uperaturation and i the prinipal mehanim in preipitation

Nuleation (2) - Homogeneou nuleation our with uperaturated olution in the abene of foreign matter, uh a dut moleule in the olution firt aoiate to form a luter, whih may dioiate or grow if a luter get large enough to take on the appearane of a lattie truture, it beome an embryo table rytalline nuleu rate of homogeneou nuleation B o é ê -16p v Aexp ê 2 3 ê 3v ( RT) éln ë ë N 2 3, L a ( ) B o : rate of homogeneou primary nuleation, number of nulei/m 3 A : frequeny fator N a : Avogadro number ù û 2 ù ú ú úû

Nuleation (3) Seondary nuleation - Superaturated olution ontain rytal - Key in ommerial rytallizer, where rytalline urfae are preent and large rytal are deired - Initiated by (1) fluid hear pat rytal urfae that weep away nulei, (2) olliion of rytal with eah other, and (3) olliion of rytal with metal urfae (rytallizer veel wall or agitator blade) - (2) and (3) mehanim are ontat nuleation: mot ommon ine they happen at the low value of relative uperatuation,, typial of indutrial appliation - Empirial power-law funtion B k M N N b j r T B : rate of eondary nuleation M T : ma of rytal per volume of magma N : agitation rate (e.g., rpm of an impeller) The ontant k N, b, j, and r are determined from experiment

Crytal Growth (1) Ma tranfer theory of rytal growth baed on equilibrium at the rytal olution interfae [Noye and Whitney, 1897] dm dt k A( - ) dm/dt : rate of ma depoited on the rytal urfae k : ma-tranfer oeffiient A : urfae area of the rytal : ma olute onentration in the bulk uperaturated olution : olute ma onentration in the olution at aturation Thin, tagnant film of olution adjaent to the rytal fae through whih olute moleular diffuion take plae [Nernt, 1904] k D d D : diffuivity d : film thikne Two-tep theory of rytal growth, diffuion-reation theory [Valeton, 1924] - Firt tep: ma tranfer of olute from the bulk of the olution to the rytal-olution interfae our dm dt k A( - ) i : uperaturated onentration at the interfae i

Crytal Growth (2) - Seond tep: a firt-order reation i aumed to our at the rytal-olution interfae, in whih olute moleule are integrated into the rytal-lattie truture (kineti tep) dm dt k A( - ) dm dt i i A( - ) 1 k + 1 k i - At low veloitie, growth rate may be ontrolled by the firt tep (ma tranfer) - The eond tep (kineti) i important when olution veloity pat the rytal urfa i high: k > k i - The ma-tranfer oeffiient, k, for the firt tep i independent of the rytallization proe - The kineti oeffiient, k i, i unique to the rytallization proe

Crytal Growth (3) - Srew-diloation theory a diloation i an imperfetion in the rytal truture predit a growth rate proprotional to ( i - ) 2 at low uperaturation and to ( i - ) at high uperaturation

Crytal Growth (4) - Although rytal do not grow a phere, an equation an be derived for the diameter of a pherial rytal dm dt K A( - ) K : overall ma-tranfer oeffiient ddp 2K ( - ) dt r p D A p D 2 p and m 6 If growth rate i ontrolled by k i, whih i aumed to be independent ( ) DDp 2ki D of D p Dt r Crytal-ize inreae i linear in time for a ontant uperaturation If growth rate i ontrolled by k at a low veloity K k 2D D D : olute diffuivity p ddp 4D( D) dt D r p 2K ( D) r 3 p r

Crytal Growth (5) Integrating from D po to D p 2 2 Dp - Dp 4D( D) o t 2 r If D po D p, D p æ 8D( D) t ö ç è r ø 1 2 In thi ae, the inreae in rytal diameter low with time At higher olution veloitie where k till ontrol, K k C D 1 2 1 p For thi ae, the inreae in rytal diameter alo low with time, but not a rapidly a predited in the previou ae