MODELLING AND COMPUTATION OF IRREGULAR NON- SPHERICAL PARTICLES TRANSPORT IN CONFINED TURBULENT FLOW

Similar documents
Numerical Simulation of Elongated Fibres in Horizontal Channel Flow

Strategy in modelling irregular shaped particle behaviour in confined turbulent flows

Effect of geometry on flow structure and pressure drop in pneumatic conveying of solids along horizontal ducts

STATIC, STAGNATION, AND DYNAMIC PRESSURES

5.1 Heat removal by coolant flow

UNIT V BOUNDARY LAYER INTRODUCTION

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom

Wall treatment in Large Eddy Simulation

1 JAXA Special Pblication JAXA-SP-1-E Small-scale trblence affects flow fields arond a blff body and therefore it governs characteristics of cross-sec

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.

Program Burn Algorithms Based on Detonation Shock Dynamics

Discontinuous Fluctuation Distribution for Time-Dependent Problems

MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA. PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 13 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS

By Dr. Salah Salman. Problem (1)

PARTICLE FLOW MODELLING ON SPIRAL CONCENTRATORS: BENEFITS OF DENSE MEDIA FOR COAL PROCESSING?

IJAET International Journal of Application of Engineering and Technology ISSN: Vol.1 No.1

Chebyshev Affine Arithmetic Based Parametric Yield Prediction Under Limited Descriptions of Uncertainty

FLUID MECHANICS D203 SAE SOLUTIONS TUTORIAL 2 APPLICATIONS OF BERNOULLI SELF ASSESSMENT EXERCISE 3

Computations of the Shock Waves at Hypersonic Velocities Taken into Account the Chemical Reactions that Appear in the Air at High Temperatures

Two identical, flat, square plates are immersed in the flow with velocity U. Compare the drag forces experienced by the SHADED areas.

MEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY

The prediction of turbulence intensities in unsteady flow

A KINEMATIC WAVE MODEL FOR EMERGENCY EVACUATION PLANNING

Modelling with Partial Differential Equations in Metrology

Module 4. Analysis of Statically Indeterminate Structures by the Direct Stiffness Method. Version 2 CE IIT, Kharagpur

FRONT TRACKING FOR A MODEL OF IMMISCIBLE GAS FLOW WITH LARGE DATA

Numerical Simulation of Three Dimensional Flow in Water Tank of Marine Fish Larvae

A random walk model to simulate the atmospheric dispersion of radionuclide

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method

Large Eddy Simulation Of Flow Past A Two-dimensional Hill

DILUTE GAS-LIQUID FLOWS WITH LIQUID FILMS ON WALLS

Lecture contents. Metals: Drude model Conductivity frequency dependence Plasma waves Difficulties of classical free electron model

The Euler-Lagrange Approach for Steady and Unsteady Flows. M. Sommerfeld. www-mvt.iw.uni-halle.de. Title. Zentrum für Ingenieurwissenschaften

A GNF framework for turbulent flow models of drag reducing fluids and proposal for a k ε type closure

Quantitative Characterization of Pressure-related Turbulence Transport Terms using Simultaneous Nonintrusive Pressure and Velocity Measurement*

FEA Solution Procedure

Process Modeling of Wellbore Leakage for GCS Risk Assessment

WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY

FOUR-WAY COUPLED SIMULATIONS OF TURBULENT

Interrogative Simulation and Uncertainty Quantification of Multi-Disciplinary Systems

Chapter 9 Flow over Immersed Bodies

Sources of Non Stationarity in the Semivariogram

Incompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System

Bed-Load Transport Rate Based on the Entrainment Probabilities of Sediment Grains. by Rolling and Lifting. Jun-De Li*

Forces on Particles and Bubbles

Direct Numerical Simulation Of Particle Saltation In Turbulent Channel Flow

The Role of Millionshschkov's Zero-Fourth Cumulant Hypothesis in Homogenous Isotropic Turbulence

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance


CHAPTER 5 INTRODUCTION TO OCEANIC TURBIDITY CURRENTS 5.1 INTRODUCTION

Study of the diffusion operator by the SPH method

Lattice-Boltzmann vs. Navier-Stokes simulation of particulate flows

CFD-Simulation thermoakustischer Resonanzeffekte zur Bestimmung der Flammentransferfunktion

Proceedings of the 51st Anniversary Conference of KSME PHYSICAL MODELING OF ATMOSPHERIC FLOW AND ENVIRONMENTAL APPLICATIONS B. R.

Gravity Waves in Shear and Implications for Organized Convection

Mathematical and Numerical Modeling of Tsunamis in Nearshore Environment: Present and Future

Steady State and Transient Thermal Analysis of Switched Reluctance Machine

PhysicsAndMathsTutor.com

A priori error analysis for transient problems using Enhanced Velocity approach in the discrete-time setting.

Fluid Dynamics. Type of Flows Continuity Equation Bernoulli Equation Steady Flow Energy Equation Applications of Bernoulli Equation

Computational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model

Experimental Study of an Impinging Round Jet

Shock wave structure for Generalized Burnett Equations

Appendix A: The Fully Developed Velocity Profile for Turbulent Duct Flows

Turbulence and boundary layers

Prandl established a universal velocity profile for flow parallel to the bed given by

Turbulence Deposition

QGSJET-II: physics, recent improvements, and results for air showers

EXPERIMENTAL STUDY OF PARAMETERS INFLUENCING THE PICKUP VELOCITY IN PNEUMATIC CONVEYING SYSTEMS

The analysis of scoring matrix in voting system with multiple criteria and multiple places

Euler-Lagrange modeling of a gas-liquid stirred reactor with consideration. of bubble breakage and coalescence

Consistent Numerical Model for Wind Buffeting Analysis of Long-Span Bridges

Lewis number and curvature effects on sound generation by premixed flame annihilation

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers

Sound speed profile structure and variability measured over flat terrain

Proceedings of the ASME nd International Conference on Ocean, Offshore and Arctic Engineering OMAE2013 June 9-14, 2013, Nantes, France

arxiv: v1 [physics.flu-dyn] 4 Sep 2013

Viscous Dissipation and Heat Absorption effect on Natural Convection Flow with Uniform Surface Temperature along a Vertical Wavy Surface

Applicability Limits of Operational Modal Analysis to Operational Wind Turbines

Kragujevac J. Sci. 34 (2012) UDC 532.5: :537.63

FEA Solution Procedure

High speed analysis of high pressure combustion in a constant volume cell

Numerical Modeling of Powder Flow during Coaxial Laser Direct Metal Deposition Comparison between Ti-6Al-4V Alloy and Stainless Steel 316L

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS

LOCAL DEPOSITION SITES OF DRUG PARTICLES IN A HUMAN NASAL CAVITY

Outline. Model Predictive Control: Current Status and Future Challenges. Separation of the control problem. Separation of the control problem

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently

On the Optimization of Numerical Dispersion and Dissipation of Finite Difference Scheme for Linear Advection Equation

Dynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining

Pressure coefficient evaluation on the surface of the SONDA III model tested in the TTP Pilot Transonic Wind Tunnel

ON THE SIMILARITIES OF THE ENGINEERING AND ATMOSPHERIC BOUNDARY LAYERS

Numerical Simulation of Particle Concentration in a Gas Cyclone Separator *

An effect of the averaging time on maximum mean wind speeds during tropical cyclone

FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL RESPONSE OF A CANTILEVER MODEL

Regression Analysis of Octal Rings as Mechanical Force Transducers

FEA Solution Procedure

Uncertainties of measurement

Modeling Effort on Chamber Clearing for IFE Liquid Chambers at UCLA

Transcription:

Martin-Lther-Universität Halle-Wittenberg 13 th Int. onf. Mltihase low in Indstrial Plants, MIP014 Setember 17-19, 014, Sestri-Levante, Italy MODELLING AND OMPUTATION O IRREGULAR NON- SPHERIAL PARTILES TRANSPORT IN ONINED TURBULENT LOW Santiago Laín b, Martin Sommerfeld a,b, Brian Qintero b and Zeeshan Qadir a a Martin-Lther-Universität Halle-Wittenberg, Zentrm für Ingenierwissenschaften,, 06099 Halle (Saale), Germany. www-mvt.iw.ni-halle.de b Universidad Atonoma de Occidente (UAO), lid Mechanics Research Gro, ali, olombia.

1. Motivation Martin-Lther-Universität Halle-Wittenberg Otline. Aerodynamic coefficients correlations 3. Wall collision arameters correlations 4. Nmerical aroach Motion eqations Stochastic modelling aroach Elementary rocess Drolast articles volme eqivalent mean diameter 40μm 5. omarison with exeriments 6. Smmary and onclsions Qartz sand volme eqivalent mean diameter 185µm

Martin-Lther-Universität Halle-Wittenberg Nmerical redictions of disersed article-laden flows sally assme sherical articles, bt in ractice articles are often non-sherical. Modelling the behavior of non-sherical articles is essential for a reliable nmerical calclation of (confined) article-laden rocesses. Lagrangian treatment of non-sherical articles Reglarly shaed non-sherical articles: Additionally solve for the article orientation and torqe Reqires information on the resistance coefficients in deendence of article orientation (e.g. Hölzer and Sommerfeld 009) Irreglar non-sherical articles: Statistical generation of the forces acting on the articles Reqires information on the PD`s of the resistance coefficients for random orientation of the articles in the flow (LBM) Re = 90 c D 1.5 Re Pa =90 1.4 1.3 1. 1.1 1.0 0 10 0 30 40 50 60 70 80 90 [degree] Additionally, the wall collision rocess is treated in a stochastic manner, based on statistical measrements of the restittion and friction coefficients.

Martin-Lther-Universität Halle-Wittenberg Resistance oefficients Determination by Lattice Boltzmann Method O U T L E T Symmetry Symmetry I N L E T P D D A U 0 P L L A U 0 0 P P T D A U T Drag Lift Pitching Torqe d l T (D VES,hll = 100 m, = 0.87) Simlations for 4 different articles = 0.87; for each Re nmber the article was randomly rotated to 71 ositions; derivation of PD`s of the coefficients omtational Domain

mean resistance coefficients [ - ] rms resistance coefficients [ - ] Martin-Lther-Universität Halle-Wittenberg orrelations Resistance oefficients D D 10 L M 1 L M 1 0.1 0.1 1 10 100 Re [ - ] x 4 Re B 1 ARe 1 Re 00 1 D / Re 1 10 100 Re [ - ] Parameter A B D D 0.56050 0.75899 0.64044 7.45770 s D -0.99615-0.005536 0.1343 0.111044 L -0.96964-0.031363 0.374803 51.7880 s L -0.966977-0.00794 0.114544 3.1660 T -0.973699-0.011175 0.116584 3.05970 s T -0.98897-0.004679 0.058604 3.158780

Martin-Lther-Universität Halle-Wittenberg Wall collision irreglar non-sherical articles (imaging exeriments) Side view To view

Normal Restittion oefficient [ - ] riction oefficient [ - ] 1.0 Martin-Lther-Universität Halle-Wittenberg Wall collision coefficients (exeriments) 1.4 0.8 0.6 1. 1.0 0.8 measrements S&H 1999 model assmtion 0.4 0. measrements mean model assmtion 0.0 0 10 0 30 40 50 60 70 80 90 Imact Angle [Degree] 0.6 0.4 0. 0.0 0 10 0 30 40 50 Imact Angle [Degree] 1 0.0 5 A e ( ) ( ) 1 B n 0.5 5 A = 0.5048, B = 3.8069, = 0.3916 s ( ) 0.07

Martin-Lther-Universität Halle-Wittenberg Nmerical Method: Eler Lagrange aroach lid Phase: Reynolds-averaged Navier-Stokes eqations in connection with the Reynolds Stress trblence model. General form of the 3D conservation eqations:

Martin-Lther-Universität Halle-Wittenberg Nmerical Method: Eler Lagrange aroach Particle Phase: simltaneos, time-deendent tracking of a large nmber of articles i i dt dx Li Di i i g m dt d m 1 i i T dt d I D D c m D ) ( 4 3 n D L L ˆ 4 1 ˆ' 8 1 3 n D T T P t ˆ 0 ˆ ˆ n t res res t ˆ' 0 ˆ' ' ˆ n t L D res

onts [ - ] Martin-Lther-Universität Halle-Wittenberg Stochastic oefficients Determination irst aroach: Gassian random rocess 80 Instantaneos flow coefficients, D, L, T 70 60 50 Re = 100 x (Re, t) x (Re ) s x (Re ) N(0,1) 40 30 0 10 Instantaneos restittion coefficient, e n : 0 0.8 1.0 1. 1.4 1.6 1.8 D [ - ] e n (, t) e ( ) s ( ) N(0,1) Instantaneos friction coefficient, : N (0, 1): standard gassian random rocess n (, t) ( ) s ( ) N(0,1) e n

Martin-Lther-Universität Halle-Wittenberg Particle-rogh wall collision momentm eqations + olomb s law of friction with modified coefficients = 1.0 Elementary rocesses 1 1 Inter-article collisions (Sommerfeld, 001) Based on the fictitios collision artner concet soltion of the imlse eqations olomb`s law of friction obliqe inelastic collision (Hard Shere Model) 1 L rel collision cylinder 1 Aroximation: inter-article collisions of isometric articles with high rotation rates are treated as sherical articles

Martin-Lther-Universität Halle-Wittenberg hannel flow laden with non-sherical articles (Kssin, 004) hannel length: 6000 mm hannel heigth: 35 mm Mean conveying velocity: 0 m/s Measring location: 5800 mm Qartz articles Eqivalent mean D = 185 m = 650 kg/m 3 Measred shericity = 0.8 Mass loading: 0.1

Normalised height [ - ] Martin-Lther-Universität Halle-Wittenberg hannel low Reslts (1) 1.0 0.8 0.6 0.4 0. Ex. articles Nm. L&S 007 Nm. non-sherical no coll. Nm. non-sherical coll. Nm. sherical coll. 0.0 0. 0.4 0.6 0.8 1.0 1. 1.4 1.6 Normalised Particle Nmber lx [ - ] There is not a large difference in the nmber flx rofiles for the varios cases.

Normalised height [ - ] Martin-Lther-Universität Halle-Wittenberg hannel low Reslts () 1.0 0.8 0.6 0.4 0. Ex. lid Ex. Particles Nm. lid Nm. L&S 007 Nm. non-sherical no coll. Nm. non-sherical coll. Nm. sherical coll. 0.0 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 Normalised Mean Velocity [ - ] or non-sherical articles comared to sherical ones the mean velocity near the bottom is lower.

Normalised Height [ - ] Martin-Lther-Universität Halle-Wittenberg hannel low Reslts (3) 1.0 0.8 0.6 0.4 0. Ex. lid Ex. Particles Nm. lid Nm. L&S 007 Nm. non-sherical no coll. Nm. non-sherical coll. Nm. sherical coll. 0.0 0.04 0.06 0.08 0.10 0.1 0.14 Normalised Streamwise rms [ - ] The stream-wise article velocity flctation is remarkably higher comared to sherical articles.

Normalised Height [ - ] Martin-Lther-Universität Halle-Wittenberg 1.0 hannel low Reslts (4) 0.8 0.6 0.4 0. Ex. lid Ex. Particles Nm. lid Nm. L&S 007 Nm. non-sherical no coll. Nm. non-sherical coll. Nm. sherical coll. 0.0 0.01 0.0 0.03 0.04 0.05 0.06 Normalised Vertical rms [ - ] The vertical article velocity flctation is remarkably lower than for sherical articles.

Martin-Lther-Universität Halle-Wittenberg Smmary and conclsions Stochastic models for flid dynamic resistance coefficients of nonsherical isometric articles have been constrcted from PD`s obtained with DNS-LBM simlations (Qadir & Sommerfeld, 013) Stochastic models for article-wall collision arameters (PD`s of normal restittion and friction coefficients) for non-sherical isometric articles have been constrcted from imaging exeriments (Saeed, 013) Pnematic conveying of qartz non-sherical isometric articles has been considered for validating the stochastic models. omarison with measrements of Kssin (004) Reasonable agreement fond for article mean and flctating velocities as well as for article nmber flx imroving revios reslts of Laín & Sommerfeld (007) Rotational and shear lift forces coefficients derived from DNS- LBM are exected to imrove the resent reslts

Martin-Lther-Universität Halle-Wittenberg Acknowledgement The financial sort of the resent stdies by the Detsche orschngs Gemeinschaft nder contract SO04/39-1 and Universidad Atónoma de Occidente throgh the joint roject Transort von nicht-shärischen Partikeln in trblenten Innenströmngen is grateflly acknowledged. Also, the sort of DAAD & olciencias throgh the German- olombian exchange mobility rogram is grateflly acknowledged.