Numerical Solutions of Population-Balance Models in Particulate Systems
|
|
- George Booker
- 5 years ago
- Views:
Transcription
1 Outline Numerical Solutions of Population-Balance Models in Particulate Systems Shamsul Qamar Gerald Warnecke Institute for Analysis and Numerics Otto-von-Guericke University Magdeburg, Germany In collaboration with Max-Planck Institute for Dynamical Systems, Magdeburg Harrachov, August 19 5, 007
2 Outline Outline 1 Aim Application Areas 3 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 4
3 Outline Aim Applications 1 Aim Application Areas 3 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 4
4 Aim Applications Aim To model and simulate nucleation, growth, aggregation and Breakage phenomena in processes engineering by solving population balance equations (PBEs). Numerical Methods To solve population balance models we use the high resolution finite volume schemes as well as their combination with the method of characteristics
5 Outline Aim Applications 1 Aim Application Areas 3 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 4
6 Industrial Applications Aim Applications Applications Pharmaceutical Chemical industries Biomedical science Aerosol formation Atmospheric physics Food industries
7 Outline 1 Aim Application Areas 3 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 4
8 f(t, x) t [G(t, x)f(t, x)] + = Q ± agg x (t, x) + Q± break (t, x) + Q+ nuc (t, x) f(0, x) = f 0, x R + :=]0, + [, t 0 1 f(t, x) is the number density function, t denotes the time and x is an internal coordinate 3 G(t, x) is the growth/dissolution rate along x, 4 Q ± α (t, x) are the aggregation, breakage and nucleation terms for α = {agg, break, nuc}. 5 The entities in the population density can be crystals, droplets, molecules, cells, and so on.
9 Q out phase boundary Q + nuc Q ± agg [Gf] x Q dis Q ± break Q in Figure: A schematic representation of different particulate processes
10 Q ± agg(t, x) = 1 x 0 0 β(t, x, x x ) f(t, x ) f(t, x x)dx β(t, x, x ) f(t, x) f(t, x )dx. Where: β = β(t, x, x ) is the rate at which the aggregation of two particles with respective volumes x and x produces a particle of volume x + x and is a nonnegative symmetric function, 0 β(t, x, x ) = β(t, x, x), x ]0, x[, (x, x ) R +.
11 Q ± break (t, x) = b(t, x, x ) S(x ) f(t, x )dx S(x) f(t, x). x b := b(t, x, x ) is the probability density function for the formation of particles of size x from particle of size x. The selection function S(x ) describes the rate at which particles are selected to break. Moments: µ i (t) = 0 x i f(t, x)dx, i = 0, 1,,, µ 0 (T) =total number of particles, µ 1 (t) =total volume of particles
12 Outline 1 Aim Application Areas 3 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 4
13 Multiply the original PBE with x and re-arrange the terms, we get f(t, x) t + [(G f)(t, x)] x (G f)(t, x) = F agg(t, x) + F break(t, x) + Q nuc, x x x f(0, x) = f0, x R +, t 0, where f(t, x) := xf(t, x), Q nuc = xq + nuc and x F agg (t, x) = 0 F break (t, x) = x u x 0 u β(t, u, v) f(t, u) f(t, v) dvdu (Filbet & Laurencot, 004) x u b(t, u, v) S(v) f(t, v) dvdu.
14 Outline 1 Aim Application Areas 3 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 4
15 Amino acid enantiomers
16 Ternary Phase Diagram Metastable zone Equilibrium point Solvent Seeding with E 1 seeds T cryst E Solubility curves M A T cryst + T E 1 E Real trajectories after seeding with E 1
17 Preferential Batch Crystallizer With Fines Dissolution preferred crystals (seeded) counter crystals (unseeded) F (k) (x, t) m (k) liq,1, mw,1 τ 1, V 1 V Fines Dissolution Loop m (k) liq,, mw, Annular Settling Zone τ, V Tank A
18 Model for Preferential Crystallization Balance for solid phase f (k) (t, x) t = G (k) (t) f (k) (t, x) x 1 τ 1 h(x)f (k) (t, x), k [p, c]. Mass balance for liquid phase in cyrstallizer dm (k) (t) dt = ṁ (k) in (t) ṁ(k) out (t) 3ρk vg (k) (t) 0 x f (k) (t, x)dx. f (k) (t, 0) = B(k) (t) G (k) (t), w(k) (t) = S (k) (t) = w(k) (t) w (k) eq m (k) (t) m (p) (t) + m (c) (t) + m W (t) 1, G (k) (t) = k g [S (k) (t)] α, k g 0, α 1.
19 B (p) 0 B (c) 0 ( ) b (t) = k(p) b S (p) (p) (t) µ (p) 3 (t) (t) = k(c) b b (c) e ln(s (c) (t)) ṁ (k) out (t) = w(k) (t)ρ liq (T) ṁ (k) in (t) = ṁ(k) out (t τ ) + k vρ τ 1 0 x 3 h(x)f (k) (t τ, x) dx.
20 Outline Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 1 Aim Application Areas 3 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 4
21 Domain Discretization Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes Regular/Irregular grid: Let N be a large integer and denote by ) i {1,,N+1} a mesh of [x min, x max ]. We set (x i 1 x 1/ = x min, x N+1/ = x max, x i+1/ = x min + i x i, i = 1,, N 1. x i 1 x i+ 1 x i 1 x i x i+1 x Here x i = (x i 1/ + x i+1/ )/, x i = x i+1/ x i 1/.
22 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes Geometric grid: x 1/ = x min, x i+1/ = x min + (i N)/q (x max x min ), i = 1,,, N where the parameter q is any positive integer. Let Ω i = [ ] x i 1/, x i+1/ for i 0. We approximate the initial data f 0 (x) in each grid cell by f i = 1 f 0 (x)dx. x i Ω i
23 Outline Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 1 Aim Application Areas 3 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 4
24 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes Method 1: Combination of MOC and FVS Let us substitute the growth rate G(t, x) by Then we have to solve: d f i dt = 1 [ (F agg ) x i (t) i+ 1 + G i+ fi ( 1 x i (t) G i+ 1 dx dt := ẋ(t) = G(t, x). (F agg ) i 1 G i 1 ] + 1 [ (F break ) x i (t) i+ 1 ) fi x i (t) + Q i (F break ) i 1 ] dx i+ 1 dt where = G i+ 1, i = 1,,, N with i.c. f(0, xi ) = f 0 (x i )
25 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes i N β(x α i,k 1/, x k ) (F agg ) i+1/ = x k (t) f k x k=0 dx fj β(x, x k ) + x dx fαi,k 1, j=α i,k Λ h x i+1/ x k j i (F break ) i+1/ = x N S(x ) fj x dx dx + O( x 3 ). k=0 Ω k j=i+1 Ω j b(x, x ) Here, the integer α i,k corresponds to the index of the cell such that x i+1/ (t) x k (t) Ω αi,k 1(t). A standard ODE-solver can be used to solve the above ODEs.
26 Outline Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 1 Aim Application Areas 3 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes 4
27 Method : Semidiscrete HR-schemes [ Integration of PBE over the control volume Ω i = Ω i f(t, x) t dx+ Ω i [G(t, x) f(t, x)] x = Ω i F agg (t, x) x Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes dx Ω i dx + Ω i x i 1 G(t, x) f(t, x) x F break (t, x) x, x i+ 1 dx ] implies dx + Let f i = f i (t) and Q i = Q i (t) be the averaged values, then we have f i t = 1 x + [ F i+ 1 [ (F break ) i+ 1 where F i+ 1 = (G f) i+ 1 F i 1 ] 1 x (F break ) i 1 [ (F agg ) i+ 1 (F agg ) i 1 ] + G i+ fi 1 + Q i, x i Ω i Q(t, x) dx. and (F agg ) & (F break ) are as given in Method 1. ]
28 Domain Discretization Numerical Method 1: Combination of MOC and FVS Numerical method : Semi-Discrete HR-Schemes The flux F i+ 1 at the right cell interface is given as (Koren, 1993): (F i + 1 ( Φ F i+ 1 and Φ is defined as: = Φ(r i+ 1 ) = max ( ( 0, min r i+ 1 r i+ 1, min ) ) (F i F i 1 ) The argument r i+ 1 of the function Φ is given as r i+ 1 = F i+1 F i + ε F i F i 1 + ε. ( ))) 3 r i+,. 1 Analogously, one can formulate the flux F i 1. Here, ε = There are several other limiting functions, namely, minmod, superbee and MC limiters, etc. Each of them leeds to a different HR-scheme (LeVeque 00, Koren 1993).
29 Example 1: All Processes Further Reading The initial data: f(0, x) = { 100 for 0.4 x 0.6, 0.01 elsewhere. B.C.: f(t, 0) = exp( 10 4 (t 0.15) ). G = 1.0, β = , b(t, x, x ) = x and S(x) = x. The exact solution in growth and nucleation case is: exp( 10 4 ((G t x) 0.15) ) for 0 x Gt, f(t, x) = 10 for 0.4 x Gt 0.6, 0.01 elsewhere. t max = 0.5 and N = 00.
30 Further Reading Results of Method 1: MOC+FVM 10 6 Pure Growth and Nucleation MOC 10 6 Nucleation + Growth + Aggregation MOC+FVM Analytical Pure Growth and Nucleation Analytical Pure Growth and Nucleation number density 10 number density volume volume Nucleation+Growth + Breakage MOC+FVM Analytical Pure Growth and Nucleation Nucleation + Growth + Aggregation + Breakage 10 6 MOC+FVM Anaytical Pure Growth and Nucleation 10 4 number density number density volume volume
31 Results of Method : FVM Further Reading 10 6 Pure Growth and Nucleation FVM 10 6 Nucleation + Growth + Aggregation FVM Analytical Pure Growth and Nucleation Analytical Pure Growth and Nucleation number density 10 number density volume volume Nucleation + Growth + Breakage FVM Analytical Pure Growth and Nucleation Nucleation +Growth + Aggregation + Breakage 10 6 FVM Analytical Pure growth and Nucleation 10 4 number density number density volume volume
32 Example : Pure Growth Further Reading The initial data are: f(0, x) = with G = 1.0, N = 100, t = 60. particle density x 10 9 Exact First Order LeVeque Koren Uniform Mesh { if 10 < x < 0, 0 elsewhere. particle density x 10 9 Exact First Order LeVeque Koren Adaptive Mesh crystal size crystal size For mesh adaptation we have used a moving mesh technique of T. Tang et al. (003)
33 Preferential Crystallization Further Reading Isothermal Case Temperature = 33C o. Non-isothermal Case T(t)[C o ] = e 7 t e 5 t t+33.
34 Further Reading Example 3: Preferential Crystallization The initial data: f (p) (0, x) = [ 1 1 πσia x exp 1 ( ) ] ln(x) µ, σ f (c) (0, x) = 0, with I a = k V ρ s µ (p) M 3 (0). s Here M s = kg is the mass of initial seeds. The maximum crystal size is x max = m with N = 500 for t = 600 min.
35 Further Reading Example 3: Preferential Crystallization particle density [#/m] x Without Fines Dissolution Initial Data First Order HR κ= 1 HR κ=1/3 particle density [#/m] 10 x 108 With Fines Dissolution Initial Data First Order 8 HR κ= 1 HR κ=1/3 6 4 particle density [#/m] crystal size [m] x 10 3 x Without Fines Dissolution Initial Data First Order HR κ= 1 HR κ=1/3 particle density [#/m] crystal size [m] x x 109 With Fines Dissolution Initial Data First Order 5 HR κ= 1 HR κ=1/ crystal size [m] x crystal size [m] x 10 3
36 Further Reading Example 3: Preferential Crystallization growth rate [m/min] 6 x 10 6 Isothermal Case p without fines c without fines 5 p with fines c with fines growth rate [m/min] 6 x 10 6 Non isothermal Case p without fines c without fines 5 p with fines c with fines nucleation rate [m/min] time [min] Isothermal Case p without fines c without fines p with fines c with fines time [min] nucleation rate [m/min] time [min] 4 x 104 Non isothermal Case p without fines c without fines p with fines c with fines time [min]
37 Further Reading Example 3: Mass Preservation in the Schemes Table: Percentage errors in mass preservation without fines dissolution. Method Isothermal Non-isothermal CPU time (s) (isothermal) N=500 N=1000 N=500 N=1000 N=500 N=1000 First order HR-κ = HR-κ = 1/ MOC Table: Percentage errors in mass preservation with fines dissolution. Method Isothermal Non-isothermal CPU time (s) (isothermal) N=500 N=1000 N=500 N=1000 N=500 N=1000 First order HR-κ = HR-κ = 1/ MOC
38 Current Project Further Reading Results of isothermal (30 C) seeded growth experiments with mandelic acid in water. Left:without counter enantiomer; Right: with counter-enantiomer (Lorenz et al., 006).
39 For Further Reading Further Reading S. Qamar, M. Elsner, I. Angelov, G. Warnecke and A. Seidel-Morgenstern A comparative study of high resolution schemes for solving population balances in crystallization. Compt. & Chem. Eng., Vol. 30, , 006. S. Qamar, and G. Warnecke Solving population balance equations for two-component aggregation by a finite volume scheme. Chem. Sci. Eng., 6, , 006. S. Qamar, and G. Warnecke Numerical solution of population balance equations for nucleation growth and aggregation processes. Compt. & Chem. Eng. (in press), 007.
40 Further Reading Thanks for your Attention
Stochastic Modeling of Chemical Reactions
Stochastic Modeling of Chemical Reactions Eric Haseltine Department of Chemical Engineering University of Wisconsin-Madison TWMCC 27 September Texas-Wisconsin Modeling and Control Consortium 22TWMCC Outline
More informationChapter 17. Fundamentals of Atmospheric Modeling
Overhead Slides for Chapter 17 of Fundamentals of Atmospheric Modeling by Mark Z. Jacobson Department of Civil & Environmental Engineering Stanford University Stanford, CA 94305-4020 August 21, 1998 Mass
More informationStabilization of continuous fluidized bed spray granulation - a Lyapunov approach
8th IFAC Symposium on Nonlinear Control Systems University of Bologna, Italy, September 1-3, 21 Stabilization of continuous fluidized bed spray granulation - a Lyapunov approach Stefan Palis 1, Achim Kienle
More informationDiscrepancy based control of systems of population balances
1st IFAC Workshop on Control of Systems Governed by Partial Differential Equations Discrepancy based control of systems of population balances Stefan Palis 1,, Andreas Bück 1, Achim Kienle 1, 1 Otto-von-Guericke-Universität,
More informationPopulation Balance Modeling
Population Balance Modeling JMBC-Particle Technology 2015 Pieter Vonk Advanced Chemical Engineering Solutions March 6 th, 2015 Overview Introduction & Industrial challenges Granulation processes Crystallization
More informationODE Math 3331 (Summer 2014) June 16, 2014
Page 1 of 12 Please go to the next page... Sample Midterm 1 ODE Math 3331 (Summer 2014) June 16, 2014 50 points 1. Find the solution of the following initial-value problem 1. Solution (S.O.V) dt = ty2,
More informationRecap: Introduction 12/1/2015. EVE 402 Air Pollution Generation and Control. Adsorption
EVE 402 Air Pollution Generation and Control Chapter #6 Lectures Adsorption Recap: Solubility: the extent of absorption into the bulk liquid after the gas has diffused through the interface An internal
More informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationHybrid CFD-Multizonal Modelling of Polymorphs and Agglomeration Phenomena in Crystallisation Processes
Hybrid CFD-Multizonal Modelling of Polymorphs and Agglomeration Phenomena in Crystallisation Processes Sean K. Bermingham 1, Stefania Trivellato 1,2, Maddalena Vernier 1,2, Fabrizio Bezzo 2,* 1 Process
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationModels for Nucleation and Condensation on Nanoparticles
PFAU X Models for Nucleation and Condensation on Nanoparticles Stefan Radl, 1 Mario Schriefl 2,3 1 TU Graz, Institute of Process and Particle Engineering 2 TU Graz, Institute of Electronics Tristan Reinisch,
More informationReview Session. The concentration of substrate that results in 50% of the maximal growth rate is 100 mg/ml.
ReviewSession Problem1 Acultureofcellsisinitiatedwith0.2mg/mlofcellmass.Theinitialglucose concentrationis1mg/ml.whatisthevalueofthemonodconstantassuming externalgrowthcontrolafter7daysofgrowth?howmuchloweristhemonod
More informationA Comparative Study of Coupled Preferential Crystallizers for the Efficient Resolution of Conglomerate-Forming Enantiomers
pharmaceutics Article A Comparative Study of Coupled Preferential Crystallizers for the Efficient Resolution of Conglomerate-Forming Enantiomers Aniruddha Majumder 1,2, * and Zoltan K. Nagy 2,3 1 School
More informationSimulation of Condensing Compressible Flows
Simulation of Condensing Compressible Flows Maximilian Wendenburg Outline Physical Aspects Transonic Flows and Experiments Condensation Fundamentals Practical Effects Modeling and Simulation Equations,
More informationErrata for SOLIDIFICATION (First Edition, 2009)
Errata for SOLIDIFICATION (First Edition, 29) J. A. Dantzig and M. Rappaz March 2, 217 Chapter 1: Overview Page 2, first line: the definition of r should have a square root. i.e., r = (ξ 2 + y 2 + z 2
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationChemical Reaction Engineering. Lecture 2
hemical Reaction Engineering Lecture 2 General algorithm of hemical Reaction Engineering Mole balance Rate laws Stoichiometry Energy balance ombine and Solve lassification of reactions Phases involved:
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationControl of Crystal Size Distribution in Batch Crystallization Processes with Size-dependent Growth using Reachability Analysis
Control of Crystal Size Distribution in Batch Crystallization Processes with Size-dependent Growth using Reachability Analysis K. Zhang, M. Nadri, C.Z. Xu, B. Uccheddu and G. Févotte Laboratoire d Automatique
More informationDevelopment of Dynamic Models. Chapter 2. Illustrative Example: A Blending Process
Development of Dynamic Models Illustrative Example: A Blending Process An unsteady-state mass balance for the blending system: rate of accumulation rate of rate of = of mass in the tank mass in mass out
More informationMath 3313: Differential Equations Second-order ordinary differential equations
Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties
More informationDynamic simulation and Control of a CO 2 Compression and Purification Unit for Oxy-Coal-Fired Power Plants
Dynamic simulation and Control of a CO 2 Compression and Purification Unit for Oxy-Coal-Fired Power Plants Authors A. Chansomwong, K.E. Zanganeh, A. Shafeen, P.L. Douglas,E. Croiset, L.A. Ricardez-Sandoval,
More informationModeling and Optimization of Semi- Interpenetrating Polymer Network (SIPN) Particle Process
Modeling and Optimization of Semi- Interpenetrating Polymer Network (SIPN) Particle Process Weijie Lin, Lorenz T. Biegler, Annette Jacobson Chemical Engineering Department Carnegie Mellon University, Pittsburgh,
More informationWell-Balanced Schemes for the Euler Equations with Gravity
Well-Balanced Schemes for the Euler Equations with Gravity Alina Chertock North Carolina State University chertock@math.ncsu.edu joint work with S. Cui, A. Kurganov, S.N. Özcan and E. Tadmor supported
More informationA particle-in-cell method with adaptive phase-space remapping for kinetic plasmas
A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas Bei Wang 1 Greg Miller 2 Phil Colella 3 1 Princeton Institute of Computational Science and Engineering Princeton University
More informationNon-equilibrium mixtures of gases: modeling and computation
Non-equilibrium mixtures of gases: modeling and computation Lecture 1: and kinetic theory of gases Srboljub Simić University of Novi Sad, Serbia Aim and outline of the course Aim of the course To present
More informationResolution of the ethanolamine salt of (±)mandelic acid by using the AS3PC method: Principle, application and results
J. Phys. IV France 122 (2004) 157-162 EDP Sciences, Les Ulis DOI: 10.1051/jp4:2004122024 Resolution of the ethanolamine salt of (±)mandelic acid by using the AS3PC method: Principle, application and results
More informationModeling Industrial Crystallizers
Modeling Industrial Crystallizers Kumar Dhanasekharan, Ph.D. Fluent Inc., 10 Cavendish Court, Lebanon, NH 03766 Phone: (603) 643-2600 x323; Email: kd@fluent.com One of the main challenges in industrial
More informationMATHCAD SOLUTIONS TO THE CHEMICAL ENGINEERING PROBLEM SET
MATHCA SOLUTIONS TO THE CHEMICAL ENGINEERING PROBLEM SET Mathematical Software Session John J. Hwalek, epartment of Chemical Engineering University of Maine, Orono, Me 4469-5737 (hwalek@maine.maine.edu)
More informationModel-based analysis of continuous. crystallization/reaction processes separating. conglomerate forming enantiomers
Model-based analysis of continuous crystallization/reaction processes separating conglomerate forming enantiomers Till Köllges and Thomas Vetter School of Chemical Engineering and Analytical Science, The
More informationDIFFUSION COEFFICIENT MEASUREMENTS USING STEFAN TUBE AT SUPERCRITICAL CONDITIONS
DIFFUSION COEFFICIENT MESUREMENTS USING STEFN TUBE T SUPERCRITICL CONDITIONS Irmak E. Ozguler, Sermin G. Sunol and ydin K. Sunol* Department of Chemical Engineering University of South Florida 4 East Fowler
More informationTHE VERIFICATION OF THE TAYLOR-EXPANSION MOMENT METHOD IN SOLVING AEROSOL BREAKAGE
144 THERMAL SCIENCE, Year 1, Vol. 16, No. 5, pp. 144-148 THE VERIFICATION OF THE TAYLOR-EXPANSION MOMENT METHOD IN SOLVING AEROSOL BREAKAGE by Ming-Zhou YU a,b * and Kai ZHANG a a College of Science, China
More informationSub -T g Relaxation in Thin Glass
Sub -T g Relaxation in Thin Glass Prabhat Gupta The Ohio State University ( Go Bucks! ) Kyoto (January 7, 2008) 2008/01/07 PK Gupta(Kyoto) 1 Outline 1. Phenomenology (Review). A. Liquid to Glass Transition
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationRadiative transfer in planetary atmospheres
Chapter 9 Radiative transfer in planetary atmospheres Warning: This chapter is just a very rough introduction to planetary atmospheres. Since the discovery of planets around other stars than the sun, the
More informationComputing phase diagrams of model liquids & self-assembly of (bio)membranes
Computing phase diagrams of model liquids & self-assembly of (bio)membranes Ulf Rørbæk Pedersen (ulf@urp.dk) Crystal structure of a simplistic model of a molecule. What is the melting temperature? A typical
More information1. Strahlungsgesetze, Ableitung der Planck-schen Strahlungsformel, Einstein-Koeffizienten, Extinktinskoeffizient, Oszillatorenstärke
1. Strahlungsgesetze, Ableitung der Planck-schen Strahlungsformel, Einstein-Koeffizienten, Extinktinskoeffizient, Oszillatorenstärke Einheiten in diesem Kapitel: diesmal cgs. Energy volume density of blackbody
More informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS 1D-time averaged models
A SHORT INTRODUCTION TO TWO-PHASE FLOWS 1D-time averaged models Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +3304 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
More informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
More information2 D. Terminal velocity can be solved for by equating Fd and Fg Fg = 1/6πd 3 g ρ LIQ = 1/8 Cd π d 2 ρ air u
ecture 8 Collection Growth of iquid Hydrometeors. I. Terminal velocity The terminal velocity of an object is the maximum speed an object will accelerate to by gravity. At terminal velocity, a balance of
More informationThermodynamic stability and transformation of pharmaceutical polymorphs*
Pure Appl. Chem., Vol. 77, No. 3, pp. 581 591, 2005. DOI: 10.1351/pac200577030581 2005 IUPAC Thermodynamic stability and transformation of pharmaceutical polymorphs* Mitsutaka Kitamura Department of Mechanical
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationChemical Engineering Applications in Scilab
Chemical Engineering Applications in Scilab Prashant Dave Indian Institute of Technology Bombay (IIT Bombay) Introduction In Chemical Engineering the type of problems that occur are Modeling and Simulation
More informationScale-up problems are often perceived as difficult. Here the reaction calorimetry has proven to be
APPLICATION OF REACTION CALORIMETRY FOR THE SOLUTION OF SCALE-UP PROBLEMS A paper from the RC User Forum Europe, Interlaken, 1995 Francis Stoessel, Ciba AG, Basel, Switzerland. Scale-up problems are often
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationESTIMATION OF CRYSTAL SIZE DISTRIBUTION OF A BATCH CRYSTALLIZATION PROCESS USING A GROWTH SIZE DEPENDENT MODEL
8th International IFAC Symposium on Dynamics and Control of Process Systems Preprints Vol, June 6-8, 7, Cancún, Mexico ESTIMATION OF CRYSTAL SIZE DISTRIBUTION OF A BATCH CRYSTALLIZATION PROCESS USING A
More informationTurbulence modulation by fully resolved particles using Immersed Boundary Methods
Turbulence modulation by fully resolved particles using Immersed Boundary Methods Abouelmagd Abdelsamie and Dominique Thévenin Lab. of Fluid Dynamics & Technical Flows University of Magdeburg Otto von
More informationStochastic Differential Equations
CHAPTER 1 Stochastic Differential Equations Consider a stochastic process X t satisfying dx t = bt, X t,w t dt + σt, X t,w t dw t. 1.1 Question. 1 Can we obtain the existence and uniqueness theorem for
More information6 Physical transformations of pure substances
6 Physical transformations of pure substances E6.b E6.2b E6.3b E6.4b Solutions to exercises Discussion questions Refer to Fig. 6.8. The white lines represent the regions of superheating and supercooling.
More informationPhysics and Thermodynamics of Water and Ice. Ottmar Möhler
Forschungszentrum Karlsruhe in der Helmholtz-Gemeinschaft Physics and Thermodynamics of Water and Ice Ottmar Möhler Institute for Meteorology and Climate Research (IMK-AAF) ESF workshop on Microbiological
More informationGROWTH RATE OF HYDROXYAPATITE CRYSTALS OBTAINED BY PRECIPITATION
STUDIA UNIVERSITATIS BABEŞ-BOLYAI, CHEMIA, XLV,, TOM I, 1 Dedicated to the memory of Prof. dr. Ioan Silaghi-Dumitrescu marking 6 years from his birth GROWTH RATE OF HYDROXYAPATITE CRYSTALS OBTAINED BY
More informationME Thermodynamics I
Homework - Week 01 HW-01 (25 points) Given: 5 Schematic of the solar cell/solar panel Find: 5 Identify the system and the heat/work interactions associated with it. Show the direction of the interactions.
More informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
More informationMicro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling
Micro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2, Giacomo Dimarco 3 et Mohammed Lemou 4 Saint-Malo, 14 décembre 2017 1 Université de
More informationPart I.
Part I bblee@unimp . Introduction to Mass Transfer and Diffusion 2. Molecular Diffusion in Gasses 3. Molecular Diffusion in Liquids Part I 4. Molecular Diffusion in Biological Solutions and Gels 5. Molecular
More informationSpeeding up path integral simulations
Speeding up path integral simulations Thomas Markland and David Manolopoulos Department of Chemistry University of Oxford Funded by the US Office of Naval Research and the UK EPSRC Outline 1. Ring polymer
More informationPopulation balances in case of crossing characteristic curves: Application to T-cells immune response
Population balances in case of crossing characteristic curves: Application to T-cells immune response Qasim Ali 1,3, Ali Elkamel 2, Frédéric Gruy 3, Claude Lambert 4, and Eric Touboul 5 1 Department of
More information1. Summary of the work and key results. Elena Simone, School of Food Science and Nutrition, University of Leeds, Leeds, UK
1 Application of process analytical technology (PAT) tools for the better understanding and control of the crystallization of polymorphic and impure systems Elena Simone, School of Food Science and Nutrition,
More informationChapter 5. On-line resource
Chapter 5 The water-air heterogeneous system On-line resource on-line analytical system that portrays the thermodynamic properties of water vapor and many other gases http://webbook.nist.gov/chemistry/fluid/
More informationOn the Modeling and Simulation of a Laser-Induced Caviation Bubble
1 On the Modeling and Simulation of a Laser-Induced Caviation Bubble Ali Zein Supervisor: Prof. Dr. Gerald Warnecke Co-supervisor: Dr. Maren Hantke Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationMath 266 Midterm Exam 2
Math 266 Midterm Exam 2 March 2st 26 Name: Ground Rules. Calculator is NOT allowed. 2. Show your work for every problem unless otherwise stated (partial credits are available). 3. You may use one 4-by-6
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Pseudo-Time Integration 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 10 Outline 1
More informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate Theory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationModelling and Simulation of a Batch Poly(Vinyl Chloride) Reactor
769 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 32, 23 Chief Editors: Sauro Pierucci, Jiří J. Klemeš Copyright 23, AIDIC Servizi S.r.l., ISBN 978-88-968-23-; ISSN 974-979 The Italian Association
More informationAppendix 4. Appendix 4A Heat Capacity of Ideal Gases
Appendix 4 W-143 Appendix 4A Heat Capacity of Ideal Gases We can determine the heat capacity from the energy content of materials as a function of temperature. The simplest material to model is an ideal
More informationA population balance approach for continuous fluidized bed dryers
A population balance approach for continuous fluidized bed dryers M. Peglow, U. Cunäus, C. Kettner, T. Metzger, E. Tsotsas, Thermal Process Engineering, University Magdeburg, 396 Magdeburg, ermany Abstract
More informationEE 143 Microfabrication Technology Fall 2014
EE 143 Microfabrication Technology Fall 2014 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 EE 143: Microfabrication
More informationSimulation of Naturally Fractured Reservoirs with Dual Porosity Models
Simulation of Naturally Fractured Reservoirs ith Dual Porosity Models Pallav Sarma Prof. Khalid Aziz Stanford University SUPRI-HW Motivations! NFRs represent at least 0%of orld reserves, but difficult
More informationA general well-balanced finite volume scheme for Euler equations with gravity
A general well-balanced finite volume scheme for Euler equations with gravity Jonas P. Berberich, Praveen Chandrashekar, Christian Klingenberg Abstract We present a second order well-balanced Godunov-type
More informationInvestigation of Packing Effect on Mass Transfer Coefficient in a Single Drop Liquid Extraction Column
Iranian Journal of Chemical Engineering Vol. 7, No. 4 (Autumn), 2010, IAChE Investigation of Packing Effect on Mass Transfer Coefficient Z. Azizi, A. Rahbar, H. Bahmanyar Engineering College, Chemical
More informationShear Flow of a Nematic Liquid Crystal near a Charged Surface
Physics of the Solid State, Vol. 45, No. 6, 00, pp. 9 96. Translated from Fizika Tverdogo Tela, Vol. 45, No. 6, 00, pp. 5 40. Original Russian Text Copyright 00 by Zakharov, Vakulenko. POLYMERS AND LIQUID
More informationAdaptive discrepancy based control of continuous fluidized bed spray granulation with internal classification
9th IFAC Symposium on Nonlinear Control Systems Toulouse, France, September 4-6, 13 ThA3.3 Adaptive discrepancy based control of continuous fluidized bed spray granulation with internal classification
More informationPhysics 208, Spring 2016 Exam #3
Physics 208, Spring 206 Exam #3 A Name (Last, First): ID #: Section #: You have 75 minutes to complete the exam. Formulae are provided on an attached sheet. You may NOT use any other formula sheet. You
More informationLiquid and solid bridges during agglomeration in spray fluidized beds
Liquid and solid bridges during agglomeration in spray fluidized beds Evangelos Tsotsas Thermal Process Engineering Otto von Guericke University Magdeburg PIKO-Workshop, 15-16 July 2014, Paderborn 1 SFB
More informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More informationThe Material Balance for Chemical Reactors
The Material Balance for Chemical Reactors Copyright c 2015 by Nob Hill Publishing, LLC 1 General Mole Balance V R j Q 0 c j0 Q 1 c j1 Conservation of mass rate of accumulation of component j = + { rate
More informationThe Material Balance for Chemical Reactors. Copyright c 2015 by Nob Hill Publishing, LLC
The Material Balance for Chemical Reactors Copyright c 2015 by Nob Hill Publishing, LLC 1 General Mole Balance V R j Q 0 c j0 Q 1 c j1 Conservation of mass rate of accumulation of component j = + { rate
More informationGradient Flows: Qualitative Properties & Numerical Schemes
Gradient Flows: Qualitative Properties & Numerical Schemes J. A. Carrillo Imperial College London RICAM, December 2014 Outline 1 Gradient Flows Models Gradient flows Evolving diffeomorphisms 2 Numerical
More informationPeculiarities of first-order phase transitions in the presence of an electric field
Peculiarities of first-order phase transitions in the presence of an electric field Yu. Dolinsky* and T. Elperin The Pearlstone Center for Aeronautical Engineering Studies Department of Mechanical Engineering
More informationFormulas for probability theory and linear models SF2941
Formulas for probability theory and linear models SF2941 These pages + Appendix 2 of Gut) are permitted as assistance at the exam. 11 maj 2008 Selected formulae of probability Bivariate probability Transforms
More informationBeer-Lambert (cont.)
The Beer-Lambert Law: Optical Depth Consider the following process: F(x) Absorbed flux df abs F(x + dx) Scattered flux df scat x x + dx The absorption or scattering of radiation by an optically active
More informationElmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid
Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.
More informationCFD Models for Polydisperse Solids Based on the Direct Quadrature Method of Moments
CFD Models for Polydisperse Solids Based on the Direct Quadrature Method of Moments Rodney O. Fox H. L. Stiles Professor Department of Chemical Engineering Iowa State University Acknowledgements: US National
More informationExam 2: Cloud Physics April 16, 2008 Physical Meteorology Questions 1-10 are worth 5 points each. Questions are worth 10 points each.
Exam : Cloud Physics April, 8 Physical Meteorology 344 Name Questions - are worth 5 points each. Questions -5 are worth points each.. Rank the concentrations of the following from lowest () to highest
More informationLecture 9. Systems of Two First Order Linear ODEs
Math 245 - Mathematics of Physics and Engineering I Lecture 9. Systems of Two First Order Linear ODEs January 30, 2012 Konstantin Zuev (USC) Math 245, Lecture 9 January 30, 2012 1 / 15 Agenda General Form
More informationNewtonian Mechanics. Chapter Classical space-time
Chapter 1 Newtonian Mechanics In these notes classical mechanics will be viewed as a mathematical model for the description of physical systems consisting of a certain (generally finite) number of particles
More informationIntroduction to Molecular Dynamics
Introduction to Molecular Dynamics Dr. Kasra Momeni www.knanosys.com Overview of the MD Classical Dynamics Outline Basics and Terminology Pairwise interacting objects Interatomic potentials (short-range
More informationAn efficient approach to stochastic optimal control. Bert Kappen SNN Radboud University Nijmegen the Netherlands
An efficient approach to stochastic optimal control Bert Kappen SNN Radboud University Nijmegen the Netherlands Bert Kappen Examples of control tasks Motor control Bert Kappen Pascal workshop, 27-29 May
More informationMelting of ice particles:
Melting of ice particles: When ice particles fall below 0 C they begin to melt, but the process takes some time since heat transfer needs to occur (heat from ambient environment has to supply the latent
More informationA Multi-Continuum Multi-Component Model for Simultaneous Enhanced Gas Recovery and CO 2 Storage in Stimulated Fractured Shale Gas Reservoirs Jiamin
A Multi-Continuum Multi-Component Model for Simultaneous Enhanced Gas Recovery and CO 2 Storage in Stimulated Fractured Shale Gas Reservoirs Jiamin Jiang M.S. Candidate Joined Fall 2013 1 Main Points Advanced
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationPower Series. Part 2 Differentiation & Integration; Multiplication of Power Series. J. Gonzalez-Zugasti, University of Massachusetts - Lowell
Power Series Part 2 Differentiation & Integration; Multiplication of Power Series 1 Theorem 1 If a n x n converges absolutely for x < R, then a n f x n converges absolutely for any continuous function
More informationTopics covered so far:
Topics covered so far: Chap 1: The kinetic theory of gases P, T, and the Ideal Gas Law Chap 2: The principles of statistical mechanics 2.1, The Boltzmann law (spatial distribution) 2.2, The distribution
More informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations
A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
More informationTime-optimal scheduling for high throughput screening processes using cyclic discrete event models
Mathematics and Computers in Simulation 66 2004 181 191 ime-optimal scheduling for high throughput screening processes using cyclic discrete event models E. Mayer a,, J. Raisch a,b a Fachgruppe System
More informationObservational constraints on mixed-phase clouds imply higher climate sensitivity
Observational constraints on mixed-phase clouds imply higher climate sensitivity TRUDE STORELVMO AND I. TAN (YALE UNIVERSITY) COLLABORATORS: M. KOMURCU (U. OF NEW HAMPSHIRE), M. ZELINKA (PNNL) TAKE-HOME
More informationexp ( κh/ cos θ) whereas as that of the diffuse source is never zero (expect as h ).
Homework 3: Due Feb 4 1. 2.11 Solution done in class 2. 2.8 The transmissivity along any dection is exp ( κh/ cos θ) where h is the slab thickness and θ is the angle between that dection and the normal
More informationEnsemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors
Ensemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors Eugenia Kalnay Lecture 2 Alghero, May 2008 Elements of Ensemble Forecasting
More information