CFD Models for Polydisperse Solids Based on the Direct Quadrature Method of Moments
|
|
- Eleanore Sanders
- 5 years ago
- Views:
Transcription
1 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 Science Foundation, US Dept. of Energy, Numerous Collaborators CFD in CRE IV: Barga, Italy June 19-24,
2 1. Introduction Population Balances Coupling with CFD Outline 2. Population Balances in CFD Population Balance Equation Direct Solvers Quadrature Methods 3. Implementation for Gas-Solid Flow Overview of MFIX Polydisperse Solids Model Application of DQMOM 4. Two Open Problems CFD in CRE IV: Barga, Italy June 19-24,
3 Population Balances Number density function (NDF) particle surface area time particle volume (mass) spatial location CFD provides a description of the dependence of n(v,a) on x For multiphase flows, the NDF will include the phase velocities (as in kinetic theory) CFD in CRE IV: Barga, Italy June 19-24,
4 Population Balances Moments of number density function Choice of k and l depends on what can be measured Solving for moments in CFD makes the problem tractable due to smaller number of scalars Multi-fluid model solves for moments from kinetic theory CFD in CRE IV: Barga, Italy June 19-24,
5 Population Balances Physical processes leading to size changes Nucleation J(x,t) produces new particles, coupled to local solubility, and properties of continuous phase Growth G(x,t) mass transfer to surface of existing particles, coupled to local properties of continuous phase Restructuring particle surface/volume and fractal dimension changes due to shear and/or physio-chemical processes Aggregation/Agglomeration particle-particle interactions, coupled to local shear rate, fluid/particle properties Breakage system dependent, but usually coupled to local shear rate, fluid/particle properties CFD provides a description of the local conditions CFD in CRE IV: Barga, Italy June 19-24,
6 Population Balances What can we compare to in-situ experiments? Sub-micron particles small-angle static light scattering Larger particles optical methods zero-angle intensity radius of gyration 1.8 < d f <3 length projected fractal dimension CFD model should predict measurable quantities accurately CFD in CRE IV: Barga, Italy June 19-24,
7 Coupling with CFD Do particles follow the flow? Stokes number Particle diameter Kinematic viscosity If St > 0.14, particle velocities must be found from a separate momentum equation in the CFD simulation CFD in CRE IV: Barga, Italy June 19-24,
8 Coupling with CFD Do PBE timescales overlap with flow timescales? Residence time Recirculation time Local mixing timescale Kolmogorov timescale CFD simulations w/o PBE can be used to determine timescales for a particular piece of equipment CFD in CRE IV: Barga, Italy June 19-24,
9 Outline 1. Introduction Population balances Coupling with CFD 2. Population Balances in CFD Population Balance Equation Direct Solvers Quadrature Methods 3. Implementation for Gas-Solid Flow Overview of MFIX Polydisperse Solids Model Application of DQMOM 4. Two Open Problems CFD in CRE IV: Barga, Italy June 19-24,
10 Population Balance Equation Typical NDF Transport Equation (small Stokes) Advection Diffusion Nucleation + Growth } Aggregation Breakage CFD in CRE IV: Barga, Italy June 19-24,
11 Population Balance Equation Aggregation Kernel Shear-induced Brownian Sub-micron aggregates: Brownian >> Shear-induced Breakage and restructuring determine fractal dimension d f In granular flow, particle-particle collisions must be added CFD in CRE IV: Barga, Italy June 19-24,
12 Population Balance Equation Breakage Kernels exponential power law Breakage due to fluid shear only ==> additional term due to collisions in gas-solid flows Parameters determined empirically and depend on chemical/physical properties of aggregates CFD in CRE IV: Barga, Italy June 19-24,
13 Population Balance Equation Daughter Distribution binary Equal sized: f = ½ Erosion: f << Number density (# / m 3 ) m(i) / m(1) f=0.5 f=0.1 f=0.01 f=0.001 CFD in CRE IV: Barga, Italy June 19-24,
14 Direct Solvers Sectional or Class Methods Finer grid N t=0 / N Finer grid R g / R g,p Time / min Time / min Accurate predictions for higher-order moments require finer grid (range: bins) CFD in CRE IV: Barga, Italy June 19-24,
15 Direct Solvers Difficulties encountered when coupled with CFD n(v; x, t) represented by N scalars n i (x,t) where 25<N<120 Depending on kernels, initial conditions, etc., source terms for these scalars can be stiff If particles are large (measured by Stokes number), multiphase models with N momentum equations required Extension to multi-variate distributions scales like N D accounting for morphology changes will be intractable Need methods that accurately predict experimentally observable moments, but at low computational cost CFD in CRE IV: Barga, Italy June 19-24,
16 Quadrature Methods Quadrature Method of Moments (QMOM) weights abscissas k th moment of CSD: CFD in CRE IV: Barga, Italy June 19-24,
17 Quadrature Methods Product-Difference algorithm (univariate CSD) Inverse problem solved on the fly in CFD simulation CFD in CRE IV: Barga, Italy June 19-24,
18 Quadrature Methods Transport 2N moments in CFD simulation Advection Diffusion Nucleation + Growth Aggregation Breakage CFD in CRE IV: Barga, Italy June 19-24,
19 Quadrature Methods Comparison with direct method R g / R g,p Fixed pivot QMOM N = 2 QMOM N = 3 QMOM N = 4 QMOM N = time, min Number density (# / m 3 ) m(i) / m(1) Using 2N = 8 scalars, QMOM reproduces the grid-independent moments of the direct method CFD in CRE IV: Barga, Italy June 19-24,
20 Quadrature Methods Multi-variate extension is straightforward weights abscissas (k,l) th moment of CSD: But inverse problem cannot be solved on the fly! CFD in CRE IV: Barga, Italy June 19-24,
21 Quadrature Methods Direct Quadrature Method of Moments (DQMOM) Weights Volume Source terms found from linear system on the fly Area CFD in CRE IV: Barga, Italy June 19-24,
22 Polydisperse Gas-Solid Flow DQMOM with size and momentum of solid phase Number Mass Momentum Source terms for mass and momentum can be found from kinetic theory for gas-solid flows Reduces to two-fluid model when α =1 CFD in CRE IV: Barga, Italy June 19-24,
23 Outline 1. Introduction Population Balances Coupling with CFD 2. Population Balances in CFD Population Balance Equation Direct Solvers Quadrature Methods 3. Implementation for Gas-Solid Flow Overview of MFIX Polydisperse Solids Model Application of DQMOM 4. Two Open Problems CFD in CRE IV: Barga, Italy June 19-24,
24 Overview of MFIX Gas-solid multi-fluid model Momentum equations Mass & energy equations Chemical species equations Population balance equations ISAT DQMOM Kinetic theory Mass & heat transfer models Detailed chemistry Aggregation, breakage and growth CFD in CRE IV: Barga, Italy June 19-24,
25 MFIX Governing Equations (I) Mass balances ( ε ρ ) + ( ε ρ u ) = t g g g g g gα n α = 1 n= 1 ( ε ρ ) + ( ε ρ u ) = t N N s sα sα sα sα sα gαn n= 1 N s M M Mass transfer from gas to solid phases Momentum balances N ( ερ g gug) + ( ερ g gugug) = σg + fgα + ερ g gg t N ( ε ρ u ) + ( ε ρ u u ) = σ -f + f + ε ρ g t α= 1 sα sα sα sα sα sα sα sα gα βα sα sα β= 1, β α g: Gas phase sα: Solid phases α=1, N Stress tensor Interaction with gas and other solid phases Body force CFD in CRE IV: Barga, Italy June 19-24,
26 MFIX Governing Equations (II) Thermal energy balances ε ε ρ T ρ C Tsα + u T = q t + H H N g g gcpg + ug Tg = qg Hg α Hrg + Hwall Twall Tg t α = 1 sα sα psα sα sα sα gα rsα Conductive heat flux Chemical species balances Heat transfer between phases N ( ε ρ X ) + ( ε ρ X u ) = R M t ( ε ρ X ) + ( ε ρ X u ) = R + M t g g gn g g gn g gn gαn α = 1 sα sα sαn sα sα sαn sα sαn gαn Heat of reaction ( ) Heat lost to walls g: Gas phase sα: Solid phases α=1, N Reactions Mass transfer CFD in CRE IV: Barga, Italy June 19-24,
27 Polydisperse Solids Model Population balance equation for solid phase L, u s Force acting to accelerate particles nl (, us; xt, ) t ( ) ( ) + u sn L,u s;x, t + u n L, (, ;, ) s F,u s;x t = S L us x t Joint size & velocity distribution function Aggregation, breakage and chemical reaction CFD in CRE IV: Barga, Italy June 19-24,
28 Direct Quadrature Method of Moments nl (, us;,) xt + snl ( s, t) + nl (, ) (, ;, ) N s s t = L s xt t u,u;x u F,u;x S u nl (, us; xt,) = ωα(,) xtδ( L Lα(,)) xt δ( us us α(,)) xt α = 1 Integrate out solid velocity nlxt (;,) + s L nl (, t) = ( Lxt ;,) t u ;x S Use distribution function N N ( L L ) s ( L L ) ( L; x, t) t ωδ α α + u α ωδ α α = S α = 1 α = 1 N α = 1 [ δ δ ] ( L L ) a '( L L )( b L a ) = S( L; x, t) N α= 1 α α α α α α Simplify equation Moment transform k k 1 al α α + klα ( bα La α α) = Sk ( xt, ) Ax = w 1 w 2 L 1 L 2 L 3 u s1 u s2 u s3 N w3 ( ;x, ) ω ( x, ) δ - ( x, ) n L t = α t L Lα t b α = 1 { ωα + ( uαωα) = aα t ωαl α ( u ω L ) t + = b α α α α + N k k mk ( x, t) = n( L;x, t) LdL ω α L α 0 α= 1 L CFD in CRE IV: Barga, Italy June 19-24,
29 Modifications to MFIX Relation between volume fractions and weights: ε = klω k v : volumetric shape factor 3 sα v α α Transport equations for volume fractions and lengths: ( ε sαρ sα) t ( ε L ρ ) sα α sα t + ( ε ρ u ) = 3k ρ L b 2k ρ L a 2 3 sα sα sα v sα α α v sα α α + ( ε L ρ u ) = 4k ρ L b 3k ρ L a 3 4 sα α sα sα v sα α α v sα α α CFD in CRE IV: Barga, Italy June 19-24,
30 CFD in CRE IV: Barga, Italy June 19-24, DQMOM Source Terms S S S S = Ax b L L L L L L L L b b a a Matrix A relates moments to weights and lengths Source term x is obtained by forcing moments to be exact
31 Aggregation and Breakage Death due to Birth due to a aggregation aggregation Bk (,) x t D a k (,) xt a a b b Sk(,) x t = Bk(,) x t Dk(,) x t + Bk(,) x t Dk(,) x t Birth due to breakage B b k (,) x t L Death due to breakage b Dk (,) xt a k /3 Bk (,) x t = n(;,) λ x t β(, u λ)( u λ ) n(;,) u x t dudλ a k Dk (,) xt = LnLxt (;,) β(, Lλ)(;,) nλ xtdλdl b 0 0 k Bk (,) x t = L a()( λ b L λ)(;,) n λ x t dλdl a k Dk (,) x t = La() L n( LxtdL ;, ) 0 Apply DQMOM 1 S x t L L L ab L a N N N N N N 3 3 k/3 k ( k ) k k (, ) = ωω i j ( i + j ) βij ωω i j i βij + ωi i i ωi i i 2 i= 1 j= 1 i= 1 j= 1 i= 1 i= 1 CFD in CRE IV: Barga, Italy June 19-24,
32 Aggregation and Breakage Kernels Aggregation and breakage kernels are obtained from kinetic theory Number of collisions: θ mi + m s j 2 Nij = πωiω jσ ijg ij ( s ) u σ ij π 2mm i j 3 Aggregation kernel: Breakage kernel: β ij Nij = ψ a ωω i j a i Nij = ψ b ω i i Efficiencies (ψ a and ψ b ) depend on temperature, particle size, etc. CFD in CRE IV: Barga, Italy June 19-24,
33 PSD Effect on Fluidization No aggregation and breakage Breakage dominant average size decreases, FB expands CFD in CRE IV: Barga, Italy June 19-24, 2005 Aggregation dominant average size Increases, FB defluidizes 33
34 Volume-Average Mean Diameter Case 1 Case 2 β = 0; a = 0 β = a = m /s; 0.1 s 1 Case 3 N = 2 filled symbols N = 3 empty symbols N = 4 lines β = a = 1s m /s; 1 CFD in CRE IV: Barga, Italy June 19-24,
35 Volume-Average Normalized Moments N = 2 filled symbols N = 3 empty symbols N = 4 lines + N k k mk ( x, t) = n( L;x, t) L dl ω α L α 0 α = 1 CFD in CRE IV: Barga, Italy June 19-24,
36 Extension to Energy/Species Balances ε ρ Thermal energy balance T C + u T = q + H H t sα sα sα psα sα sα sα gα rsα + k ρ L C c k ρ L C T a 3 3 v s α ps T, α v s α ps sα α Changes due to aggregation and breakage Multi-variate DQMOM H = ε R H rsα sα sα, C r R k c P * m sα, C = p[ ] C 2 2 CFD in CRE IV: Barga, Italy June 19-24,
37 1. Introduction Population Balances Coupling with CFD Outline 2. Population Balances in CFD Population Balance Equation Direct Solvers Quadrature Methods 3. Implementation for Gas-Solid Flow Overview of MFIX Polydisperse Solids Model Application of DQMOM 4. Two Open Problems CFD in CRE IV: Barga, Italy June 19-24,
38 Two Open Problems 1. How to extend DQMOM to systems with unknown fluxes at boundaries in phase space? Model problem: pure evaporation How can we estimate it? Estimate flux in DQMOM variables, test with exact solutions: Define vectors: Define cross product : Linear constraint: CFD in CRE IV: Barga, Italy June 19-24,
39 Simple case with monotone flux (N = 2): CFD in CRE IV: Barga, Italy June 19-24,
40 Harder case with multimode flux (N = 2): CFD in CRE IV: Barga, Italy June 19-24,
41 Harder case with N = 3: CFD in CRE IV: Barga, Italy June 19-24,
42 Two Open Problems 2. What is best choice of moments for multivariate DQMOM? Model problem: homogeneous aggregation Aggregation Aggregation + sintering Moments (mass k=1, l=0; area k=0, l=1) CFD in CRE IV: Barga, Italy June 19-24,
43 Choice of moments affects the condition of matrix C B F E D A All choices yield nearly same weights and abscissas Choose moments with lowest condition number? CFD in CRE IV: Barga, Italy June 19-24,
44 Another example: Williams Spray Equation CFD in CRE IV: Barga, Italy June 19-24,
45 Coefficients depend on choice of 5N moments: Condition number of A depends on choice of k, l, m, p In general, A matrix will become singular if 1 < l + m + p Choose l, m, p = (0,1) and vary k to yield 5N distinct moments Number: (k, l, m, p)= 0 Mass: k =1, (l, m, p) = 0 X-Mom: k =1, l = 1 Y-Mom: k = 1, m = 1 Z-Mom: k = 1, p = 1 Is there a general method for choosing moments? CFD in CRE IV: Barga, Italy June 19-24,
Prediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor
Journal of Engineering Science, Vol. 10, 95 105, 2014 Prediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor Vahid Akbari and Mohd.
More informationAggregation and Breakage of Nanoparticle Dispersions in Heterogeneous Turbulent Flows
ggregation and reakage of Nanoparticle Dispersions in Heterogeneous Turbulent Flows M. Soos,, D. Marchisio*, J. Sefcik, M. Morbidelli Swiss Federal Institute of Technology Zürich, Z CH *Politecnico di
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 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 informationPredicting Size Distribution of Dispersed Phase Systems
Predicting Size Distribution of Dispersed Phase Systems Evren Bayraktar Supervisors: Prof. S.Turek (Dortmund), Dr. Ing. Otto Mierka, Dr. Ing. Frank Platte 28.11.2008, Sulzer Page 1 Outline Introduction
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationMODENA. Deliverable 3.2. WP s leader: TU/e. Simulations for foams, dispersion and mixing and developed SW. Principal investigator:
Delivery date: 11-4-2015 MODENA Authors: Patrick Anderson, Martien Hulsen, Christos Mitrias TU/e E-mail : pda@tue.nl Deliverable 3.2 Simulations for foams, dispersion and mixing and developed SW Daniele
More informationSimulation of Particulate Solids Processing Using Discrete Element Method Oleh Baran
Simulation of Particulate Solids Processing Using Discrete Element Method Oleh Baran Outline DEM overview DEM capabilities in STAR-CCM+ Particle types and injectors Contact physics Coupling to fluid flow
More informationprecipitation in a Confined Impinging Jets Reactor by means of
Modeling of nanoparticles precipitation in a Confined Impinging Jets Reactor by means of Computational Fluid Dynamics E. Gavi, D.L. Marchisio, A.A. Barresi Politecnico di Torino, Department of Material
More informationCFD Investigations of Effects of Cohesive Particles Proportion on Fluidization of Binary Particles
Proceedings of the 2 nd World Congress on Momentum, Heat and Mass Transfer (MHMT 17) Barcelona, Spain April 6 8, 2017 Paper No. ICMFHT 122 ISSN: 2371-5316 DOI: 10.11159/icmfht17.122 CFD Investigations
More informationMultiphase Flows. Mohammed Azhar Phil Stopford
Multiphase Flows Mohammed Azhar Phil Stopford 1 Outline VOF Model VOF Coupled Solver Free surface flow applications Eulerian Model DQMOM Boiling Model enhancements Multi-fluid flow applications Coupled
More informationCOMPARISON OF ANALYTICAL SOLUTIONS FOR CMSMPR CRYSTALLIZER WITH QMOM POPULATION BALANCE MODELING IN FLUENT
CHINA PARTICUOLOGY Vol. 3, No. 4, 13-18, 5 COMPARISON OF ANALYTICAL SOLUTIONS FOR CMSMPR CRYSTALLIZER WITH QMOM POPULATION BALANCE MODELING IN FLUENT Bin Wan 1, Terry A. Ring 1, *, Kumar M. Dhanasekharan
More informationSTUDY OF AGGREGATION IN BARIUM SULPHATE PRECIPITATION
15 th International Symposium on Industrial Crystallization Sorrento, Italy, 15-18 September 2002 paper 171 STUDY OF AGGREGATION IN BARIUM SULPHATE PRECIPITATION D. L. Marchisio 1, A. A. Barresi, M. Garbero,
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationCode Verification of Multiphase Flow with MFIX
Code Verification of Multiphase Flow with MFIX A N I R U D D H A C H O U D H A R Y ( A N I R U D D @ V T. E D U ), C H R I S T O P H E R J. ROY ( C J R O Y @ V T. E D U ), ( V I R G I N I A T E C H ) 23
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationOn the Comparison between Population Balance Models for CFD Simulation of Bubble Columns
Chemical and Biological Engineering Publications Chemical and Biological Engineering 2005 On the Comparison between Population Balance Models for CFD Simulation of Bubble Columns Jay Sanyal Fluent Inc.
More informationModeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!
http://www.nd.edu/~gtryggva/cfd-course/! Modeling Complex Flows! Grétar Tryggvason! Spring 2011! Direct Numerical Simulations! In direct numerical simulations the full unsteady Navier-Stokes equations
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 informationPolitecnico di Torino. Porto Institutional Repository
Politecnico di Torino Porto Institutional Repository [Proceeding] Simulation of macromolecule self-assembly in solution: multiscale approach a Original Citation: Lavino, A; Di Pasquale, N; Carbone, P;
More informationC C C C 2 C 2 C 2 C + u + v + (w + w P ) = D t x y z X. (1a) y 2 + D Z. z 2
This chapter provides an introduction to the transport of particles that are either more dense (e.g. mineral sediment) or less dense (e.g. bubbles) than the fluid. A method of estimating the settling velocity
More informationAPPLICATION OF THE CONDITIONAL QUADRATURE METHOD OF MOMENTS FOR THE SIMULATION OF COALESCENCE, BREAKUP AND MASS TRANSFER IN GAS-LIQUID STIRRED TANKS
14 th European Conference on Mixing Warszawa, 10-13 Septemer 2012 APPLICATION OF THE CONDITIONAL QUADRATURE METHOD OF MOMENTS FOR THE SIMULATION OF COALESCENCE, BREAKUP AND MASS TRANSFER IN GAS-LIQUID
More informationComputational Models for Polydisperse Particulate and Multiphase Systems
Computational Models for Polydisperse Particulate and Multiphase Systems Providing a clear description of the theory of polydisperse multiphase flows, with emphasis on the mesoscale modeling approach and
More informationSolution of population balance equations in applications with fine particles: Mathematical modeling and numerical schemes
Chemical and Biological Engineering Publications Chemical and Biological Engineering 11-15-216 Solution of population balance equations in applications with fine particles: Mathematical modeling and numerical
More informationApplications of Computational Fluid Dynamics in the Process Industries. Ahmad Haidari & Peter Spicka Fluent Inc.
Applications of Computational Fluid Dynamics in the Process Industries Ahmad Haidari & Peter Spicka Fluent Inc. 1 Outline Overview of CFD s growth in the process Industry Overview of modeling multiphase
More informationComputational fluid dynamic modeling of fluidized bed polymerization reactors
Graduate Theses and Dissertations Graduate College 2012 Computational fluid dynamic modeling of fluidized bed polymerization reactors Ram Ganesh Rokkam Iowa State University Follow this and additional
More informationMass Transfer in Turbulent Flow
Mass Transfer in Turbulent Flow ChEn 6603 References: S.. Pope. Turbulent Flows. Cambridge University Press, New York, 2000. D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, La Caada CA, 2000.
More informationQuadrature-based moment methods for polydisperse multiphase flow modeling
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 2013 Quadrature-based moment methods for polydisperse multiphase flow modeling Cansheng Yuan Iowa State University
More informationTurbulent Flows. g u
.4 g u.3.2.1 t. 6 4 2 2 4 6 Figure 12.1: Effect of diffusion on PDF shape: solution to Eq. (12.29) for Dt =,.2,.2, 1. The dashed line is the Gaussian with the same mean () and variance (3) as the PDF at
More informationAn Angular Multigrid Acceleration Method for S n Equations with Highly Forward-Peaked Scattering
An Angular Multigrid Acceleration Method for S n Equations with Highly Forward-Peaked Scattering Bruno Turcksin & Jean C. Ragusa & Jim E. Morel Texas A&M University, Dept. of Nuclear Engineering Bruno
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationEulerian models. 2.1 Basic equations
2 Eulerian models In this chapter we give a short overview of the Eulerian techniques for modelling turbulent flows, transport and chemical reactions. We first present the basic Eulerian equations describing
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationMixing and Evaporation of Liquid Droplets Injected into an Air Stream Flowing at all Speeds
Mixing and Evaporation of Liquid Droplets Injected into an Air Stream Flowing at all Speeds F. Moukalled* and M. Darwish American University of Beirut Faculty of Engineering & Architecture Mechanical Engineering
More informationCFD Modelling of Turbulent Mass Transfer in a Mixing Channel
CFD Modelling of Turbulent Mass Transfer in a Mixing Channel Lene K. Hjertager Osenbroch, Bjørn H. Hjertager and Tron Solberg Aalborg University Esbjerg Esbjerg, Denmark Homepage: hugin.aue.auc.dk Prepared
More informationPredictive multiscale computational modeling of nanoparticles in flame reactors
Graduate Theses and Dissertations Graduate College 2012 Predictive multiscale computational modeling of nanoparticles in flame reactors Maulik Mehta Iowa State University Follow this and additional works
More informationMSc. Thesis Project. Simulation of a Rotary Kiln. MSc. Cand.: Miguel A. Romero Advisor: Dr. Domenico Lahaye. Challenge the future
MSc. Thesis Project Simulation of a Rotary Kiln MSc. Cand.: Miguel A. Romero Advisor: Dr. Domenico Lahaye 1 Problem Description What is a Rotary Kiln? A Rotary Kiln is a pyroprocessing device used to raise
More informationOutline. Advances in STAR-CCM+ DEM models for simulating deformation, breakage, and flow of solids
Advances in STAR-CCM+ DEM models for simulating deformation, breakage, and flow of solids Oleh Baran Outline Overview of DEM in STAR-CCM+ Recent DEM capabilities Parallel Bonds in STAR-CCM+ Constant Rate
More informationMicromechanics of Colloidal Suspensions: Dynamics of shear-induced aggregation
: Dynamics of shear-induced aggregation G. Frungieri, J. Debona, M. Vanni Politecnico di Torino Dept. of Applied Science and Technology Lagrangian transport: from complex flows to complex fluids Lecce,
More informationANSYS FLUENT12.0. Population Balance. Module Manual
ANSYS FLUENT12.0 Population Balance Module Manual January 2009 Copyright c 2009 by ANSYS, Inc. All Rights Reserved. No part of this document may be reproduced or otherwise used in any form without express
More informationCFD Simulation of Nanosufur Crystallization Incorporating Population Balance Modeling
Iranica Journal of Energy & Environment 4 (1) Special Issue on Nanotechnology: 36-42, 2013 ISSN 2079-2115 IJEE an Official Peer Reviewed Journal of Babol Noshirvani University of Technology DOI: 10.5829/idosi.ijee.2013.04.01.06
More informationTABLE OF CONTENTS CHAPTER TITLE PAGE
v TABLE OF CONTENTS CHAPTER TITLE PAGE TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES LIST OF SYMBOLS LIST OF APPENDICES v viii ix xii xiv CHAPTER 1 INTRODUCTION 1.1 Introduction 1 1.2 Literature Review
More informationAn open-source quadrature-based population balance solver for OpenFOAM
Chemical and Biological Engineering Publications Chemical and Biological Engineering 2-2-2018 An open-source quadrature-based population balance solver for OpenFOAM Alberto Passalacqua Iowa State University,
More informationMultiscale Modeling of TiO2 Nanoparticle Production in Flame Reactors: Effect of Chemical Mechanism
Chemical and Biological Engineering Publications Chemical and Biological Engineering 2010 Multiscale Modeling of TiO2 Nanoparticle Production in Flame Reactors: Effect of Chemical Mechanism Maulik Mehta
More informationmeters, we can re-arrange this expression to give
Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.
More informationFiltered Two-Fluid Model for Gas-Particle Suspensions. S. Sundaresan and Yesim Igci Princeton University
Filtered Two-Fluid Model for Gas-Particle Suspensions S. Sundaresan and Yesim Igci Princeton University Festschrift for Professor Dimitri Gidaspow's 75th Birthday II Wednesday, November 11, 2009: 3:15
More informationBASIC DESIGN EQUATIONS FOR MULTIPHASE REACTORS
BASIC DESIGN EQUATIONS FOR MULTIPHASE REACTORS Starting Reference 1. P. A. Ramachandran and R. V. Chaudhari, Three-Phase Catalytic Reactors, Gordon and Breach Publishers, New York, (1983). 2. Nigam, K.D.P.
More informationMulti-Scale Modeling of Turbulence and Microphysics in Clouds. Steven K. Krueger University of Utah
Multi-Scale Modeling of Turbulence and Microphysics in Clouds Steven K. Krueger University of Utah 10,000 km Scales of Atmospheric Motion 1000 km 100 km 10 km 1 km 100 m 10 m 1 m 100 mm 10 mm 1 mm Planetary
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationINTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING---
ILASS-Europe 2002 Zaragoza 9 11 September 2002 INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING--- Z. Tapia, A. Chávez e-mail: ztapia@imp.mx Instituto Mexicano del Petróleo Blvd. Adolfo Ruiz Cortines
More informationCalculating equation coefficients
Fluid flow Calculating equation coefficients Construction Conservation Equation Surface Conservation Equation Fluid Conservation Equation needs flow estimation needs radiation and convection estimation
More informationTurbulence Solutions
School of Mechanical, Aerospace & Civil Engineering 3rd Year/MSc Fluids Turbulence Solutions Question 1. Decomposing into mean and fluctuating parts, we write M = M + m and Ũ i = U i + u i a. The transport
More informationOverview of Accelerated Simulation Methods for Plasma Kinetics
Overview of Accelerated Simulation Methods for Plasma Kinetics R.E. Caflisch 1 In collaboration with: J.L. Cambier 2, B.I. Cohen 3, A.M. Dimits 3, L.F. Ricketson 1,4, M.S. Rosin 1,5, B. Yann 1 1 UCLA Math
More informationSolution of population balance equations in applications with fine particles: mathematical modeling and numerical schemes
Solution of population balance equations in applications with fine particles: mathematical modeling and numerical schemes Tan Trung Nguyen, Frédérique Laurent, Rodney Fox, Marc Massot To cite this version:
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationChapter 7 Mixing and Granulation
Chapter 7 Mixing and Granulation 7.1 Mixing and Segregation (Chapter 9) Mixing vs. segregation (1) Types of Mixture * Perfect mixing Random mixing Segregating mixing Figure 9.1 (2) Segregation 1) Causes
More informationNUMERICAL MODELING OF FINE PARTICLE FRACTAL AGGREGATES IN TURBULENT FLOW
THERMAL SCIENCE, Year 2015, Vol. 19, No. 4, pp. 1189-1193 1189 NUMERICAL MODELING OF FINE PARTICLE FRACTAL AGGREGATES IN TURBULENT FLOW by Feifeng CAO a, Zhanhong WAN b,c*, Minmin WANG b, Zhenjiang YOU
More informationFLUENT 6.3. Population Balance Module Manual
FLUENT 6.3 Population Balance Module Manual September 2006 Copyright c 2006 by Fluent Inc. All rights reserved. No part of this document may be reproduced or otherwise used in any form without express
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More informationLesson 6 Review of fundamentals: Fluid flow
Lesson 6 Review of fundamentals: Fluid flow The specific objective of this lesson is to conduct a brief review of the fundamentals of fluid flow and present: A general equation for conservation of mass
More informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationVici Progress Report, March 2013.
Vici Progress Report, March 2013. Nicolás Rivas Project description Two general approaches are used to model granular materials: microscopic, where information is taken from individual particles, and macroscopic,
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
More informationCOMPUTATIONAL STUDY OF PARTICLE/LIQUID FLOWS IN CURVED/COILED MEMBRANE SYSTEMS
COMPUTATIONAL STUDY OF PARTICLE/LIQUID FLOWS IN CURVED/COILED MEMBRANE SYSTEMS Prashant Tiwari 1, Steven P. Antal 1,2, Michael Z. Podowski 1,2 * 1 Department of Mechanical, Aerospace and Nuclear Engineering,
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationTable of Contents. Preface... xiii
Preface... xiii PART I. ELEMENTS IN FLUID MECHANICS... 1 Chapter 1. Local Equations of Fluid Mechanics... 3 1.1. Forces, stress tensor, and pressure... 4 1.2. Navier Stokes equations in Cartesian coordinates...
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationHydrodynamics of Bubble Columns: Turbulence and Population Balance Model
Article Hydrodynamics of Bubble Columns: Turbulence and Population Balance Model Camila Braga Vieira *, and Pierre Proulx ID ID, Giuliana Litrico, Ehsan Askari ID, Gabriel Lemieux ID Chemical Engineering
More informationExplaining and modelling the rheology of polymeric fluids with the kinetic theory
Explaining and modelling the rheology of polymeric fluids with the kinetic theory Dmitry Shogin University of Stavanger The National IOR Centre of Norway IOR Norway 2016 Workshop April 25, 2016 Overview
More informationSteady Flow and its Instability of Gravitational Granular Flow
Steady Flow and its Instability of Gravitational Granular Flow Namiko Mitarai Department of Chemistry and Physics of Condensed Matter, Graduate School of Science, Kyushu University, Japan. A thesis submitted
More informationCHARACTERIZATION OF BRANCHED POLYMERS IN SOLUTION (I)
CHARACTERIZATION OF BRANCHED POLYMERS IN SOLUTION (I) Overview: General Properties of Macromolecules in Solution Molar Mass Dependencies Molar Mass Distributions Generalized Ratios Albena Lederer Leibniz-Institute
More informationMicro-scale modelling of internally
Micro-scale modelling of internally unstable soils Dr Tom Shire School of Engineering, University of Glasgow 1 st September 2017 Outline Internal instability Micro-scale modelling Hydromechanical criteria
More informationLecture 35: The Inertia Tensor
Lecture 35: The Inertia Tensor We found last time that the kinetic energy of a rotating obect was: 1 Trot = ωω i Ii where i, ( I m δ x x x i i, k, i, k So the nine numbers represented by the I i tell us
More informationIntroduction to Heat and Mass Transfer. Week 10
Introduction to Heat and Mass Transfer Week 10 Concentration Boundary Layer No concentration jump condition requires species adjacent to surface to have same concentration as at the surface Owing to concentration
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationChapter 7. Basic Turbulence
Chapter 7 Basic Turbulence The universe is a highly turbulent place, and we must understand turbulence if we want to understand a lot of what s going on. Interstellar turbulence causes the twinkling of
More informationCOMPARISON OF DIFFERENT SUBGRID TURBULENCE MODELS AND BOUNDARY CONDITIONS FOR LARGE-EDDY-SIMULATIONS OF ROOM AIR FLOWS.
7 TH INTRNATINAL CNFRNC N AIR DISTRIBTIN IN RMS, RMVNT 2 pp. 31-36 CMPARISN F DIFFRNT SBGRID TRBLNC MDLS AND BNDARY CNDITINS FR LARG-DDY-SIMLATINS F RM AIR FLWS. D. Müller 1, L. Davidson 2 1 Lehrstuhl
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More informationCALCULATION OF SHOCK STAND-OFF DISTANCE FOR A SPHERE
J. Comput. Fluids Eng. Vol.17, No.4, pp.69-74, 2012. 12 / 69 CALCULATION OF SHOCK STAND-OFF DISTANCE FOR A SPHERE IN NONEQUILIBRIUM HYPERSONIC FLOW M. Ahn Furudate * Dept. of Mechatronics Engineering,
More informationLecture 3: 1. Lecture 3.
Lecture 3: 1 Lecture 3. Lecture 3: 2 Plan for today Summary of the key points of the last lecture. Review of vector and tensor products : the dot product (or inner product ) and the cross product (or vector
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationLattice Boltzmann Methods for Fluid Dynamics
Lattice Boltzmann Methods for Fluid Dynamics Steven Orszag Department of Mathematics Yale University In collaboration with Hudong Chen, Isaac Goldhirsch, and Rick Shock Transient flow around a car LBGK
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationA Simple Turbulence Closure Model
A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: Mean velocity: Turbulent velocity: Gradient operator: Advection operator: V =
More informationEffect of Turbulent Enhancemnt of Collision-coalescence on Warm Rain Formation in Maritime Shallow Convection
Effect of Turbulent Enhancemnt of Collision-coalescence on Warm Rain Formation in Maritime Shallow Convection A. A. WYSZOGRODZKI 1 W. W. GRABOWSKI 1,, L.-P. WANG 2, AND O. AYALA 2 1 NATIONAL CENTER FOR
More informationOn the hydrodynamic diffusion of rigid particles
On the hydrodynamic diffusion of rigid particles O. Gonzalez Introduction Basic problem. Characterize how the diffusion and sedimentation properties of particles depend on their shape. Diffusion: Sedimentation:
More informationDEVELOPMENT OF A NUMERICAL APPROACH FOR SIMULATION OF SAND BLOWING AND CORE FORMATION
TMS (The Minerals, Metals & Materials Society), DEVELOPMENT OF A NUMERICAL APPROACH FOR SIMULATION OF SAND BLOWING AND CORE FORMATION G.F. Yao, C. W. Hirt, and
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 informationSound Pressure Generated by a Bubble
Sound Pressure Generated by a Bubble Adrian Secord Dept. of Computer Science University of British Columbia ajsecord@cs.ubc.ca October 22, 2001 This report summarises the analytical expression for the
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationThe Simulation of Wraparound Fins Aerodynamic Characteristics
The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:
More informationParticle Dynamics: Brownian Diffusion
Particle Dynamics: Brownian Diffusion Prof. Sotiris E. Pratsinis Particle Technology Laboratory Department of Mechanical and Process Engineering, ETH Zürich, Switzerland www.ptl.ethz.ch 1 or or or or nucleation
More informationSTRESS TRANSPORT MODELLING 2
STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1 Introduction In the previous lecture we introduced
More informationComputational fluid dynamics simulation of precipitation processes
Retrospective Theses and Dissertations 2004 Computational fluid dynamics simulation of precipitation processes Liguang Wang Iowa State University Follow this and additional works at: http://lib.dr.iastate.edu/rtd
More informationFOUR-WAY COUPLED SIMULATIONS OF TURBULENT
FOUR-WAY COUPLED SIMULATIONS OF TURBULENT FLOWS WITH NON-SPHERICAL PARTICLES Berend van Wachem Thermofluids Division, Department of Mechanical Engineering Imperial College London Exhibition Road, London,
More informationNumerical Solutions of Population-Balance Models in Particulate Systems
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
More information11.D.2. Collision Operators
11.D.. Collision Operators (11.94) can be written as + p t m h+ r +p h+ p = C + h + (11.96) where C + is the Boltzmann collision operator defined by [see (11.86a) for convention of notations] C + g(p)
More information