Basic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos

Similar documents
Turbulent Flow. Turbulent Flow

Diffusion Mass Transfer

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

Turbulent Nonpremixed Flames

Lecture 12. Modeling of Turbulent Combustion

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

A Study on Responces of Premixed Flames of Hydrogen, Methane and Propane having the same Adiabatic Flame Temperature to Flame Stretch

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 6

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson

Appendix II Summary of Important Equations

TURBULENT FLOW A BEGINNER S APPROACH. Tony Saad March

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Equation of State Modeling of Phase Equilibrium in the Low-Density Polyethylene Process

Introduction to Turbulence Modeling

Numerical Study of Propane-Air Mixture Combustion in a Burner Element

1. Why turbulence occur? Hydrodynamic Instability. Hydrodynamic Instability. Centrifugal Instability: Rayleigh-Benard Instability:

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Problem adapted reduced models based on Reaction-Diffusion Manifolds (REDIMs)

Design Equations. ν ij r i V R. ν ij r i. Q n components. = Q f c jf Qc j + Continuous Stirred Tank Reactor (steady-state and constant phase)

CHEMICAL ENGINEERING

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

Journal of Fluid Science and Technology

The Tangential Force Distribution on Inner Cylinder of Power Law Fluid Flowing in Eccentric Annuli with the Inner Cylinder Reciprocating Axially

Effects of Polymer Concentration and Molecular Weight on the Dynamics of Visco-Elasto- Capillary Breakup

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram

Lecture 7 Fluid Systems I. System Analysis Spring

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Airflow and Contaminant Simulation with CONTAM

GENERAL EQUATIONS OF PHYSICO-CHEMICAL

Chemical Equilibrium. Chapter 6 Spontaneity of Reactive Mixtures (gases) Taking into account there are many types of work that a sysem can perform

CFD Simulation of Dense Gas Extraction through Polymeric Membranes

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

Turbulence. Lecture 21. Non-linear Dynamics. 30 s & 40 s Taylor s work on homogeneous turbulence Kolmogorov.

Convection Heat Transfer. Textbook: Convection Heat Transfer. Reference: Convective Heat and Mass Transfer. Convection Heat Transfer

3. Be able to derive the chemical equilibrium constants from statistical mechanics.

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

Mass Transfer Processes

DynaflowC.F.D. conference

Introduction to the lattice Boltzmann method

Mixing in an agitated tubular reactor. J.J. Derksen. School of Engineering, University of Aberdeen, Aberdeen, UK.

Lecture 5.8 Flux Vector Splitting

Modeling of Dynamic Systems

#64. ΔS for Isothermal Mixing of Ideal Gases

Turbulent Flow in Curved Square Duct: Prediction of Fluid flow and Heat transfer Characteristics

Application of the Adjoint Method for Vehicle Aerodynamic Optimization. Dr. Thomas Blacha, Audi AG

Tensor Smooth Length for SPH Modelling of High Speed Impact

The Governing Equations

A STUDY ON THE PHYSICS OF SUPERSONIC MIXING FIELD WITH INJECTION AT DIFFERENT ANGLES

Axial Turbine Analysis

ChE 512: Topic 1 Reactions at a fluid non-porous solid interface. P.A. Ramachandran

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Experimental and Numerical Investigation of Air Suction in Domestic Gas-Burning Heaters to Increase Efficiency

Ph.D. Qualifying Examination in Kinetics and Reactor Design

Energy configuration optimization of submerged propeller in oxidation ditch based on CFD

CHEMICAL REACTIONS AND DIFFUSION

Survey of applications of discrete vortex method in civil engineering

Aerodynamics. Finite Wings Lifting line theory Glauert s method

The Finite Element Method

Energy, Entropy, and Availability Balances Phase Equilibria. Nonideal Thermodynamic Property Models. Selecting an Appropriate Model

and Statistical Mechanics Material Properties

Numerical analysis on the agglomeration behavior of fine particles in plane jet

Mass transfer in multi-component mixtures

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

A NUMERICAL STUDY OF HEAT TRANSFER AND FLUID FLOW IN A BANK OF TUBES WITH INTEGRAL WAKE SPLITTER

DNS of transition in hypersonic boundary-layer flows including high-temperature gas effects

LOW-NOx CATALYTIC COMBUSTION OF HYDROGEN

PART I: MULTIPLE CHOICE (32 questions, each multiple choice question has a 2-point value, 64 points total).

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014

10/9/2003 PHY Lecture 11 1

Macroscopic Momentum Balances

Numerical Analysis of Heat Transfer and Pressure Drop in a Channel Equipped with Triangular Bodies in Side-By-Side Arrangement

Thermodynamics General

Lecture 5. Stoichiometry Dimensionless Parameters. Moshe Matalon N X. 00 i. i M i. i=1. i=1

Numerical simulation of the combustion of strong swirling confined reacting flow (natural gas /air) in gas turbine combustor

Natural Gas Conversion in Monolithic Catalysts: Interaction of Chemical Reactions and Transport Phenomena

One-Dimensional Numerical Investigation of a Cylindrical Micro-Combustor Applying Electrohydrodynamics Effect

FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO

Wilbur and Ague 4 WILBUR AND AGUE; APPENDIX DR1. Two-dimensional chemical maps as well as chemical profiles were done at 15 kv using

Introduction to Computational Fluid Dynamics

Assistant professor, Faculty of Mechanical Engineering, Shiraz University, Shiraz, Iran **

a for save as PDF Chemistry 163B Introduction to Multicomponent Systems and Partial Molar Quantities

FEATURES OF TURBULENT TRANSPORT OF MOMENTUM AND HEAT IN STABLY STRATIFIED BOUNDARY LAYERS AND THEIR REPRODUCTION IN ATMOSPHERIC MESOSCALE MODELS

IC Engine Flow Simulation using KIVA code and A Modified Reynolds Stress Turbulence Model

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

2) For a two-dimensional steady turbulent flow in Cartesian coordinates (x,y), with mean velocity components (U,V), write

Numerical Differentiation

LAGRANGIAN MECHANICS

Calculation of Aerodynamic Characteristics of NACA 2415, 23012, Airfoils Using Computational Fluid Dynamics (CFD)

Process Modeling. Improving or understanding chemical process operation is a major objective for developing a dynamic process model

Week 9 Chapter 10 Section 1-5

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

A Computational Viewpoint on Classical Density Functional Theory

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

Angular Momentum and Fixed Axis Rotation. 8.01t Nov 10, 2004

MEEM 3700 Mechanical Vibrations

Simulation of Turbulent Flow Using FEM

Design of a Porous Bluff-Body Disc on Improving the Gas-Mixing Efficiency

Transcription:

Introducton to Toshhsa Ueda School of Scence for Open and Envronmental Systems Keo Unversty, Japan Combuston Mxng and reacton n hgh vscous flud Applcaton of Chaos Keo Unversty 1 Keo Unversty 2 What s reactve flows? The flows wth reactons combuston ar polluton chemcal reactor bo reactor Keo Unversty 3 Keo Unversty 4

Most mportant pont of reactve fluds Bass of mult-component flud flud conssts of reactants and products speces concentraton changes locally and temporally wth heat release or absorpton Chemcal reacton mechansm Flud dynamcs Thermodynamcs Transport phenomena Keo Unversty 5 Keo Unversty 6 Concentraton and fracton Understandng reactve flows Concentraton: amount / unt volume mass concentraton m [kg/m 3 ] molar concentraton c [mol/m 3 ] Fracton: rato of concentraton mass fracton Y = m m j molar fracton X = c c j flow moton velocty and pressure u, v, w, p energy transport temperature T speces transport mass fracton of chemcal speces Y Keo Unversty 7 Keo Unversty 8

Governng equatons Conservaton equatons of physcal moton Conservaton equatons of physcal moton mass, momentum, energy, chemcal speces Equatons of chemcal varatons Chemcal reacton formula Rate equatons for chemcal reacton concentraton dependence temperature dependence Keo Unversty 9 Mass: Momentum: energy: chemcal speces: ρ ρu ρv ρw + + + = 0 t x y z Du p τxx τyx τzx ρ = + + + Dt x x y z Dh Dp ρ = q + ρ U h + Dt Dt dh = C dt DY ρ = ( ρ U ) + w Dt + Φ ( : for prefect gases) p Keo Unversty 10 Reacton model Reacton model CH 2 H O + CO 4 + O2 2 2 Reacton rate equaton 3 dcch4 a b ω [ m s] = = kc c Concentraton : Temperature : CH 2 4 kmol/ CH 4 O2 dt E = a = n E a k A'exp AT exp R0T R0T Keo Unversty 11 Passve scalar or actve scalar? Passve scalar : scalar varables passve for flow moton Actve scalar : scalar varables actve for flow moton constant propertes Passve scalar varable propertes Actve scalar Keo Unversty 12

Passve: constant property Actve: varable property Mass: Momentum: ρ ρu ρv ρw + + + = 0 t x y z Du p τxx τyx τzx ρ = + + + Dt x x y z energy: u, v, w, p Dh Dp ρ = q + ρ U h + + Φ Dt Dt DY chemcal ρ = ( ρ U ) + w Dt speces: Keo Unversty 13 Mass: Momentum: u, v, w, p ρ ρu ρv ρw + + + = 0 t x y z Du p τxx τyx τzx ρ = + + + Dt x x y z energy: Dh Dp ρ = q + ρ U h + + Φ Dt Dt DY chemcal ρ = ( ρ U ) + w Dt speces: (Y, T ), (Y, T ) Keo Unversty 14 Non-dmensonal parameter for reactve flows τ f Damköhler number : Da = τ f :characterstc tme of flud moton (mxng) c :characterstc tme of chemcal reacton c Damköhler number τ >> Da f τ c τ f Da = τ The chemcal reacton s much faster than the flud moton, e.g. engnes, hgh temperature chemcal reactors Flud moton (mxng) determnes the progress n chemcal reacton the chemcal reacton nstantly occurs n reacton surface, as soon as reactants meet c Keo Unversty 15 Reactant A Reactant B Product Reacton surface Keo Unversty 16

Damköhler number τ << Da 0 f τ c τ f Da = τ The flud moton s much faster than the chemcal reacton, e.g. low temperature chemcal reactors, boreactors Rate of chemcal reacton determnes the progress n chemcal reacton Reactant A Reactant B Chemcal reacton gradually occurs after Reactants mx well c Product Mxng Mxng due to Molecular dffuson Molecular level mxng necessary for chemcal reacton τ = du µ dy du d ( ρu Newton s s law for momentum transfer τ = µ = ν ) dy dy dt d ρc p Fourer's law for heat conducton q = k = α Fck s law for speces dffuson j = dy dm D dy ( T ) dy Keo Unversty 17 Keo Unversty 18 Mxng Mxng due to advecton Turbulent advecton (Turbulent mxng) mxng by mult-scale vortex moton Chaotc advecton (Chaotc mxng) mxng by stretchng and foldng Mxng When the vscosty s low Mxng due to turbulence s effectve When the vscosty s hgh It s not easy to make the flow turbulent It s necessary to enhance the mxng even n lamnar flow condton Chaotc mxng Keo Unversty 19 Keo Unversty 20

Non-steady moton opens new world Perodc moton precse control s possble Turbulent moton transport, such as mxng, s enhanced most reactve systems use turbulent moton engnes furnaces reactors Keo Unversty 21 Non-steady moton opens new world Chaotc advecton enhance lamnar mxng n low Reynolds number area low velocty heat sland ssue low shear bo reactor small area mcro reactor hgh vscous fluds polymer technology non-newtonan fluds artfcal organs Keo Unversty 22 Combuston Combuston plays a key role n many energy devces Combuston Engnes Furnaces Dynamcs of Combuston s a future technology for mprovng combuston Keo Unversty 23 Bunsen flame wth burner rotaton Keo Unversty 24

Combuston Combuston over a Clathrate Hydrate Combuston Combuston over a hydrate can be predcted numercally The flame front propagetes upstream drecton Ar Hydrate Keo Unversty 25 Keo Unversty 26 Mxng and reacton n hgh vscous fluds Test secton 190 50 Mxng and reacton n hgh vscous flud Outer cylnder Inner cylnder Keo Unversty 27 o Reactants Workng flud (glycerne) Sde vew Keo Unversty 28

Mxng and reacton n hgh vscous fluds Parameters Workng flud Glycern 1300ml (98wt%) Rato of rad Eccentrcty R R o ε = 0.26 d R o R + [ Fe( CN) ] KFe[ Fe( CN) ] + 3 + Chemcal reacton Fe 3 + K 4 6 6 K Yellow Red Dark blue Mxng and reacton n hgh vscous fluds Flow Feld stretchng chaotc regon Chaotc regon Regular regon foldng regular regon The mxng occurs well The mxng occurs scarcely Keo Unversty 29 Keo Unversty 30 Mxng and reacton n hgh vscous fluds Effect of ntal locatons Both reactants are ntally located n the regular regon Mxng and reacton n hgh vscous fluds Effect of ntal locatons regular-regular regular-chaos chaos-chaos One reactant s n regular, the other n chaotc regon Both reactants are ntally located n the chaotc regon Product Keo Unversty 31 Keo Unversty 32

Applcaton of Chaos Applcaton of Chaos Non-element mxer and reactor Non-element mxer and reactor s realzed by usng chaotc advecton (chaotc mxng) Hgh vscosty mxng and reacton n lamnar flow condton Mcro scale and low speed Mcro, mllmeter scale machnes Low shear polymer flud, bo flud Keo Unversty 33 Keo Unversty 34 Applcaton of Chaos Summary ø10 50 50 50 50 50 50 Stretchng foldng Combuston Mxng and reacton n lamnar flow condton Mcro, mllmeter scale machnes Applcaton of chaos Dffuson normal to the flow drecton s an mportant, nterestng and challengng feld Keo Unversty 35 Keo Unversty 36