Numerical Simulation of the Dispersion of a Plume of Exhaust Gases from Diesel and Petrol Engine Vehicles

Similar documents
EVALUATION OF FORCE COEFFICIENTS FOR A 2-D ANGLE SECTION USING REALIZABLE k-ε TURBULENCE MODEL

Effect of a Vector Wall on the Thermal Field in a SRU Thermal Reactor

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

Transient Numerical of Piston Wind in Subway Station. Haitao Bao

Numerical simulation of a solar chimney power plant in the southern region of Iran

Polymerization Technology Laboratory Course

Diffusion of Heptane in Polyethylene Vinyl Acetate: Modelisation and Experimentation

NUMERICAL SIMULATION AND EXPERIMENTAL INVESTIGATION FOR INDOOR AIR ENVIRONMENT IN AN OFFICE ROOM

Multi-Fuel and Mixed-Mode IC Engine Combustion Simulation with a Detailed Chemistry Based Progress Variable Library Approach

Numerical Studies on Lip Shock Flow Behaviors over Backward Facing Sharp Edge Step with Hybrid RANS-LES

EFFECT OF HEAT FLUX RATIO FROM BOTH SIDE-WALLS ON THERMAL- FLUID FLOW IN CHANNEL

Comparison of LES and RANS in numerical simulation of turbulent non-premixed flame under MILD combustion condition

Lecture 9: Dynamic Properties

Available online at Physics Procedia 32 (2012 )

Solution in semi infinite diffusion couples (error function analysis)

The profile-linear average velocity for the transition in pipes based on the method of LES *

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

NUMERICAL SIMULATION OF A PREMIXED TURBULENT V-SHAPED FLAME

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

COLD GAS ANALYSIS OF A WASTE-GAS INCINERATOR TO ENHANCE MIXING CAPABILITIES USING CFD

Experimental and Numerical Investigation of Temperature Distribution in Room with Displacement Ventilation

A Paper presentation on. Department of Hydrology, Indian Institute of Technology, Roorkee

HEAT FLUX MEASUREMENT OF URBAN BOUNDARY LAYERS IN KYOTO CITY AND ITS PREDICTION BY CFD SIMULATION

A General Lagrangian Approach to Simulate Pollutant Dispersion in Atmosphere for Low-wind Condition

On the turbulence models and turbulent Schmidt number in simulating stratified flows

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2)

INVESTIGATION ON THE PRESSURE MATCHING PERFORMANCE OF THE CONSTANT AREA SUPERSONIC-SUPERSONIC EJECTOR. Jian CHEN, Zhenguo WANG*, Jiping WU, Wanwu XU

Extended Progress Variable and Mixture Fraction Model for Turbulent Premixed Combustion

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

Numerical Simulation on Supersonic Turbulent Flow past Backward Facing Rounded Step Utilizing Hybrid RANS-LES

Computer modelling of technogenic thermal pollution zones in large water bodies

Multiphase CFD Modeling of Trickle-Bed Reactor Hydrodynamics

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

2.1 Constitutive Theory

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

P R = P 0. The system is shown on the next figure:

FI 3103 Quantum Physics

Linear Response Theory: The connection between QFT and experiments

11.2 ESTIMATES OF AIR AND POLLUTANT EXCHANGE FOR STREET CANYONS IN COMBINED WIND- BUOYANCY DRIVEN FLOW USING THE RANS RENORMALIZATION GROUP

Mist/air Cooling in a Two-Pass Rectangular Rotating Channel with 45-deg Angled Rib Turbulators

ELASTIC MODULUS ESTIMATION OF CHOPPED CARBON FIBER TAPE REINFORCED THERMOPLASTICS USING THE MONTE CARLO SIMULATION

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

FTCS Solution to the Heat Equation

PHYS 1443 Section 001 Lecture #4

A Computational Fluid Dynamics Study of Turbulence, Radiation, and Combustion Models for Natural Gas Combustion Burner

VENTILATION AND DEDUSTING OF MELTING SHOPS

Motion of Wavepackets in Non-Hermitian. Quantum Mechanics

On One Analytic Method of. Constructing Program Controls

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

Chapters 2 Kinematics. Position, Distance, Displacement

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

Spurious oscillations and conservation errors in interface-capturing schemes

Mechanics Physics 151

Motion in Two Dimensions

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

Mechanics Physics 151

Thermochemistry study of internal combustion engine

Nonequilibrium models for a multi component reactive distillation column

Transcription: Messenger RNA, mrna, is produced and transported to Ribosomes

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

Anisotropic Behaviors and Its Application on Sheet Metal Stamping Processes

T q (heat generation) Figure 0-1: Slab heated with constant source 2 = q k

Robustness Experiments with Two Variance Components

INFLUENCE OF MODELING OF TURBULENCE IN THE FLOW PARAMETERS WITHIN A FOOD OVEN USING THE OPENFOAM

Thermal fatigue estimation due to thermal stratification in the RCS branch line using one-way FSI scheme

Molecular Dynamics Simulation Study forgtransport Properties of Diatomic Liquids

SHOALING OF NONLINEAR INTERNAL WAVES ON A UNIFORMLY SLOPING BEACH

Genetic Algorithm in Parameter Estimation of Nonlinear Dynamic Systems

Measurement of liquid holdup and axial dispersion in trickle bed reactors using radiotracer technique

APPLICATION OF THE SCHWARZ-CHRISTOFFEL TRANSFORMATION IN SOLVING TWO-DIMENSIONAL TURBULENT FLOWS IN COMPLEX GEOMETRIES

Single-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method

Comprehensive Integrated Simulation and Optimization of LPP for EUV Lithography Devices

ARWtr 2004 Modern Transformers October. Vigo Spain Transformers

CST036. -SST Turbulence Model for Separated Particle-Laden Flows

WiH Wei He

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

WAVE BREAKING USING A ROLLER APPROACH IN A HYBRID FINITE-VOLUME FINITE-DIFFERENCE BOUSSINESQ-TYPE MODEL

Effect of Terrain on Mass and Heat Transport of Smoke Plume Based on Geometric Properties of Digital Elevation Model and Satellite Imagery

( ) () we define the interaction representation by the unitary transformation () = ()

TSS = SST + SSE An orthogonal partition of the total SS

Variants of Pegasos. December 11, 2009

A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION

Implementation of Quantized State Systems in MATLAB/Simulink

Numerical analysis on the distribution features of velocity circulation of axial-flow pump

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model

First-order piecewise-linear dynamic circuits

Density Matrix Description of NMR BCMB/CHEM 8190

Notes on the stability of dynamic systems and the use of Eigen Values.

National Exams December 2015 NOTES: 04-BS-13, Biology. 3 hours duration

Chapter 2 Linear dynamic analysis of a structural system

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :

CFD modeling of precipitation of nanoparticles in Confined Impinging Jet Reactors

Graduate Macroeconomics 2 Problem set 5. - Solutions

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

Tools for Analysis of Accelerated Life and Degradation Test Data

Robustness of DEWMA versus EWMA Control Charts to Non-Normal Processes

Transcription:

World Academy of Scence, Engneerng and Technology 67 01 Numercal Smulaon of he Dsperson of a Plume of Exhaus Gases from Desel and Perol Engne Vehcles H. ZAHLOUL, and M. MERIEM-BENZIANE Absrac The obecve of hs wor s o sudy numercally he dsperson characerscs of desel engne vehcle polluans a an dle condon n an dealzed and smplfed envronmen. The Fluen code wh he - eddy dsspaon urbulen model was used for he numercal smulaon of he parculae concenraons of exhaus plume NOx, CO and HC. The smulaed resuls mach very well wh he expermenal ones close o he source of emsson bu decay o he amben concenraons much slower. Despe he dfference n he nal concenraons, he polluans decayed exponenally n all drecons. The effecs of he nal emsson concenraon and ex velocy have been nvesgaed. Keywords Concenraon, dsperson, exhaus plume, polluan I. INTRODUCTION MISSIONS from moor vehcles are he man source of Ear polluon n mos urban ces. The dsperson of polluans s hghly dependen on local weaher and drvng condons. A he cy level, he ncrease n urbulence and he dmnuon of he sably of he boundary layer over urban areas have resuled n changng he dsperson of polluans from auomobles and road raffc. One mporan feaure of he urbulen flow s s ably o ranspor and mx he polluan much more effcenly han a lamnar flow. Idlng emssons are consdered more harmful because he polluans canno be dspersed by he wae creaed by a movng vehcle. A specfc dle mode s one n whch he engne vehcle s pared along he sree wh he moor swched on. They have been varous numercal and expermenal sudes on he dsperson of polluans emed by moor vehcles. Mos of he expermenal ones were carred ou eher n wnd unnels or on roadsde [1], []. For he numercal nvesgaons, hey were concenraed on he developmen of a numercal model for smulang he dsperson process In hs wor, he FLUENT code wh he - urbulen model was used o sudy he dsperson of polluan emed by a moor vehcle whch s operang a dle under dealzed and smplfed condons. The obaned resuls were analyzed and compared o he expermenal ones [3]. H. Zahloul s wh he Deparmen of Mechancal Engneerng, Unversy Hassba Benboual of Chlef, Algera (phone : 13-771797990; fax : 13-771794; e-mal : zahloulh@yahoo.fr) M. Merem Benzane s wh he Deparmen of Mechancal Engneerng, Unversy Hassba Benboual of Chlef, Algera, (e-mal : mbmadd001@yahoo.fr). II. GOVERNING EQUATIONS In order o descrbe he process of he dsperson of exhaus polluans n amosphere accuraely, he effecs of urbulen mxng, convecon, dffuson, emperaure varaons and speces ranspor mus be consdered. The governng equaons for he problem suded can be wren as [4]. A. Connuy equaon : ρ r + dv( u ) = 0, (1) Where ρ s he flud average densy; U s he mean velocy componen n he drecon of coordnae. Momenum equaon : U U ρu U + pδ + U 3 δ + ρg + () ( ρu U ) = 0 Where p s he flud average pressure, s he flud lamnar vscosy, g s he gravaonal acceleraon n he drecon of coordnae. In hs sudy, he Reynolds sresses n Eq. () can be wren as U U U U U (3) ( ρ ) = + δ x x 3 Where s he flud urbulen vscosy and can be expressed as = c (4) / Where c s a consan of he model. and can be solved by he urbulen energy ( ) equaon and he urbulen dsspaon ( ) equaon : B. Turbulen energy( ) equaon : ( ) U U ρ Γ + U U + ρ = 0 (5) 885

World Academy of Scence, Engneerng and Technology 67 01 C. Turbulen dsspaon equaon : ( ) U U Γ + U ρu C + C ρ 0 1 = (6) Where C 1 and C are consans of he model, Γ and be deermned va Bousssnesq approxmaons. D. Energy conservaons equaon Γ can T Dp u ( ρh) + ( ρuh) = ( + ) Σ( h J, ) + + τ = 0 D u Where h s he sac enhalpy, s he molecular conducvy, s he effecve molecular conducvy and can be expressed as = / Pr, s ' he dffuson flux J, of he h chemcal speces due o he concenraons gradens of he h speces. E. Speces ranspor equaon: J ', ( C ) + ( uc ) = + R, (8) ρ Where C s he mean concenraon of he h speces and R s he source erm due o he chemcal reacons. The ranspor equaons descrbed n hs secon were dscrezed by he fne volume mehod and he compuaons were performed usng he FLUENT code. III. NUMERICAL SIMULATION In hs par, we focus on he feasbly of usng a numercal code for flud dynamcs o smulae he concenraons of gaseous polluans (NOx, CO, HC). Then, we compare he resuls of our calculaons wh measuremens aen from he leraure [3] A. Expermenal daa [3] (7) Fg. 1 Schemac of caresan coordnae sysem and vehcular exhaus plume geomery Expermens were carred ou n a sem-enclosed space ha has suffcen space for he plume o dsperse freely n all drecons. The expermenal se permed repeaable measuremens n he vehcular exhaus plume under realsc and reproducble condons.the expermens were also performed wh he followng parameers: - Sable weaher condons (wea wnd) and n he evenng o avod convecon curren due o he hea of he sun. - Temperaure: 0- C -Relave humdy: 50-60% -No movng vehcle bu wh runnng moor Inal condons: a) Enry of flud: All he properes of he flud a he nle are nown [1]. b) Speed: - The ex velocy of polluan = 4.8 m / s. - The ar velocy n sable weaher condon = 0.3 m / s. c) Temperaure: - The emperaure a he oule of he exhaus ppe s 380 K. - The emperaure of he compuaonal doman s 93 K. d) Concenraon: -producs of he exhaus gases from he combuson of desel and perol engnes are gven by he followng able: The expermens were carred ou o gan nsgh no he evoluon of he vehcular polluans nsgh he feld regon of he exhaus plume of wo moor vehcles (Desel and Perol) under dle condons and oban daa for esng he performance of he numercal mode. In he expermens, he es vehcles were operaed under dle condon on he es se and experenced no prevalng wnd. The caresann coordnae sysem (x, y, z) s fxed on he ground, as shown n Fg. 1. The orgn of he coordnae sysem s on he cenerlne of he alppe ex of he vehcle. The hegh of he alppe above he ground s H. 886

World Academy of Scence, Engneerng and Technology 67 01 TABLE II Gaseous polluans concenraons for each vehcle [3] Engne NO X (ppm) CO (ppm) HC (ppm) VEHICLE A DIESEL 140 60 4 VEHICLE B PETROL 6 80 40 In he process of smulaon, he dmensons of he negraed doman were aen o be 15m x 6m x 3m. The dameer of he alppe was 0.03m and he hegh of he alppe above he ground was 0.3m. The exernal boundary of he negraed doman was defned as free sream boundary. Fg. 3-a Comparson of NOx concenraon beween he expermenal daa and smulaed resuls IV. RESULTS AND DISCUSSION A. Concenraons of CO, NOx and HC of vehcle A (Desel) Fg 3-b Dsrbuon of NOx concenraon [ppm] Fg. -a Comparson of CO concenraon beween he expermenal daa and smulaed resuls Fg. 4 Comparson of HC concenraon beween he daa and smulaed resuls expermenal Fg -b Dsrbuon of CO concenraon [ppm] We noce ha he smulaed resuls show a good agreemen wh he expermenal daa Fg. a and fg. 3a show ha he concenraons of CO and NOx decreased rapdly along he downsream dsance of he vehcular exhaus plume. The concenraons of CO and NOx a he alppe ex were 55 and 140 by volume respecvely and dropped exponenally o less han 1 ppm whn a dsance of 3.0m. Smlarly, n fg. 4, he HC concenraon decreased rapdly along he downsream dsance. Ths concenraon was 4.5 and dropped exponenally o less han 1 ppm from a dsance of m. 887

World Academy of Scence, Engneerng and Technology 67 01 B. Concenraons of CO and NOx of vehcle B (perol) Fg. 5 Comparson of CO concenraon beween he expermenal daa and smulaed resuls Fg. 8 Comparson of HC concenraon beween vehcle A and B Fg. 7 shows ha he desel vehcle (vehcle A)has he hghes level of NOx concenraon, whle he perol engne vehcle B) has he lowes of hese polluans. Furhermore, Fg. 8 shows ha he desel engne has he lowes concenraon of HC emsson whle he perol engne has he hghes of hese polluans. I can be concluded ha here are sgncan dfferences of polluan concenraons n he exhaus plume for dfferen vehcles under dle condon and he dfference depends on he ype of he vehcle, he vehcle echnologe level and he condons of he vehcle. Fg. 6 Comparson of NOx concenraon beween he expermenal daa and smulaed resul D. Smulaon resul of flow velocy of he exhaus plume We noce ha he smulaed resuls show a good agreemen wh he expermenal daa Fg. 5 and fg. 6 show ha he concenraons of CO and NOx decreased rapdly along he downsream dsance of he vehcular exhaus plume. The concenraons of of CO and NOx a he alppe ex were 57 and 5.5 by volume respecvely and dropped exponenally o less han 1 ppm whn a dsance of 3.5m C. Comparson of NOx and HC concenraons beween vehcle A and B (numercal and expermenal) Fg. 9 Flow velocy profle along X axs Fg 10 Dsrbuon of flow velocy of he exhaus plume Fg. 7 Comparson of NOx concenraon beween Vehcle A and B Fg. 9 and fg. 10 show ha he evoluon of polluan plume velocy n he longudnal drecon X. One noces ha here s a sudden decrease of he flow speed of exhaus gases. I decreases exponenally wh a value of 4.8m/s a he ex of he alppe o reach he ar amben speed for a dsance of m. 888

World Academy of Scence, Engneerng and Technology 67 01 V. CONCLUSION Idlng emssons of a vehcle s consdered more harmful on human healh because he polluans canno be dspersed by he wae creaed by a movng vehcle. Ths paper nvesgaes numercally he dsperson characerscs of he exhaus plume when a vehcle s operang under dle condon by consderng he dsperson of wo vehcles (desel and perol) from s pon of emsson o some dsance downsream. Also, he CFD code FLUENT was used o analyze he dsperson of gaseous polluan usng he - eddy dsspaon model. The model was valdaed wh NOx, CO and HC concenraons measured expermenally. The model was able o smulae closely he expermenal resuls close o he source bu he smulaed resuls oo a longer dsance o decay o amben values. REFERENCES [1] K. T., Knapp, S. B., Teada, Cenral Calorna vehcle parculae emsson sudy, Fnal repor US, Envronnemenal Proecon Agency, 000, MD-46 [] D. B. Kelson, W. F. Was, Fne parcle emsson on Mnnoesa hghways, Fnal repor, Mnnesoa, 001, MN/RC-001-1 [3] Z. Nng, C. S., Cheung, Y. Lu, Expermenal and numercal sudy of he dperson of moor vehcle polluans under dle condon, Amospherc envronmen, vol.39, 005, pp 7880-7893 [4] Km e all, On he predcon of concenraon varaons n a dspersng heavy-duy ruc exhaus plume usng - urbulen closure, Amospherc envronmen, vol.35, 001, pp 567-575 889