A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS
|
|
- Camron Allen
- 5 years ago
- Views:
Transcription
1 International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical an Prouction Engineering Nanyang Technological niversity, 5 Nanyang Avenue, Singapore (Receive February ; Accepte September ) ABSTRACT Analytical moel is elope to preict smoe layer temperatures after activation of sprinler or water mist system. The governing equations are solve to etermine water roplet trajectory, temperature history an aporation rate. As a result, a rate of heat absorption by water spray an temperature history of the hot smoe layer can be preicte. A goo agreement is emonstrate between the analytical moel an the results of full CF simulations.. INTROCTION The problem of preicting smoe layer temperature is typically arising in many fire safety applications. A part of this general problem is a preiction of smoe layer behaviour upon activation of fire suppression systems. Typical examples of such systems are water sprinlers an water mist nozzles. Heat absorption by water spray generally leas to the rop of the overall temperature insie the compartment. Quantitative preiction of this effect is quite important for estimation of attainable conitions insie builings uring fire ents. Primary ifferences between water sprinler an water mist systems are in water ischarge rates an roplet size istribution. Water spray for conventional sprinler is relatively coarse, with significant amount of roplets larger than mm. In contrast, mist spray is efine by NFPA as the spray containing 99% of its volume in the roplets smaller than mm. Recent interest in water mist application is mostly riven by phasing out of halons, an the nee to fin a replacement which can act as volumetric suppression agent. ue to very fine ispersion, mist is capable of flooing the room, although it oes not have as goo flooing properties as gaseous agents. ue to high surface area of the roplet phase, water mist aporates very effectively, an therefore water ischarge quantities may be ept much smaller, compare to conventional sprinlers. Water mist spray oes not generally penetrate to the surface of burning material (ue to low momentum of iniviual roplets), an suppress flame in gaseous phase. In contrast, sprinler sprays are most effective for unshiele fuel surfaces, as roplets easily penetrate to burning material, cool the surface, an suppress pyrolysis process. Computational Flui ynamics (CF) moelling of water sprinler sprays is well-establishe [-3]. Hower, this technique is relatively expensive. Typically, tens of thousans of trajectories are trace to represent water spray statistically. For real fire engineering applications, this may not be necessary in all situations. Basic unerstaning of water spray behaviour can be achie using simple estimations of roplet ynamics, similar to consierations applie in the pioneering wor of Chow an Tang [4]. The moel escribe in the present paper provies a more accurate estimation of roplet interaction with the hot layer of combustion proucts. In contrast to the paper [4], heat transfer problem for the moving roplet is solve in aition to the ynamic problem. The major avancement of the moel is the analytical solution for the roplet motion, obtaine for the case of real rag coefficient correlations, applicable for spherical particles moving in a turbulent gas stream. The equations an the metho to obtain analytical solutions are iscusse first. From the solutions, the roplet trajectories an heat absorption rates for iniviual roplets can be preicte. Such preictions are illustrate for various roplet sizes an ifferent temperatures of the hot layer. Finally, the moel is utilize to mae preictions of temperature histories in compartments after activation of water sprinlers. The results from the present analytical moel are compare to the CF simulations of the similar problem. 95
2 International Journal on Engineering Performance-Base Fire Coes. RESLTS AN ISCSSION. Sprinler Spray ynamics within the Smoe Layer Consier an interaction of a sprinler spray with the hot layer of certain thicness an temperature T (Fig. ). It is assume that absorption of heat by the spray within the layer is uniform. This assumption is reasonably accurate if the with of the spray is comparable with the length scale of the fire ceiling jet, as has been emonstrate in [7]. The roplets injecte into the hot smoe layer unergo two major stages: () heat-up stage an () aporation stage. The overall heat absorption by a single roplet epens on its resience time within the layer. If the resience time is large, then the roplet may be completely aporate. roplet ynamics an heat transfer equations must be solve simultaneously to obtain roplet trajectory an temperature history. In orer to elop a simplifie analytical solution, let us consier the stages () an () separately... Heat-up stage roplet ynamics: uring the heat-up stage, the roplet is consiere as having constant mass, i.e. there is no aporation. This assumption can be shown to be quite accurate [5]. The general equation of the roplet motion through a gas layer (neglecting forces other than rag an gravity) may be written as the following set of two ifferential equations for the velocity components [6]: m t V m t π ρc 8 π ρc 8 + V + V V + mg () For a typical conventional sprinler the Volumetric Mean iameter is of the orer of to mm. Taing the characteristic roplet size as mm, Reynols particle number for the roplets exiting orifice at velocities of to ms - may be estimate to be Re ~ 7 to 5. In this range the rag coefficient for spherical particles is a wea function of the particle Reynols number [6] an varies between.47 an.6. Therefore, it may be represente by a constant average value with a reasonable accuracy. It is obvious that en with the constant rag coefficient, the equations () are still couple. Hower, since roplets hit the sprinler eflector, their vertical velocity is usually small compare to the tangential velocity. Since the rag coefficient is assume to be constant, then in the case V <<, equations () become ecouple an may be solve separately. The system is simplifie as follows: m t V m t π ρc 8 π ρc 8 V + mg () Sprinler hea roplet trajectories Smoe layer Opening Fig. : Schematic view of smoe layer an roplet trajectories in a compartment 96
3 International Journal on Engineering Performance-Base Fire Coes The first equation contains only tangential component, an can be easily solve by separating variables. The solution has the following form: π = C t + ρ 8m (3) Substitution of this expression to the secon (V component) equation yiels: V π π m ρc C t V + mg t 8 + ρ 8m (4) To fin solution of (4) one nees to solve homogeneous equation first (which is one by separating variables). Then a metho of variation of parameters may be applie to fin solution of the non-homogeneous equation. The final solution of equation (4) can be written as: g κ V κ g V = + ( t + κ) (5) t + κ where 4 ρ p κ = 3C ρ (6) If the roplet leaves the smoe layer uring its heat-up stage, then the following equation may be easily erive to estimate the roplet exit time t * from the layer: * ( t + κ) * g t g * κ V κ ln + + t = (7) κ 4 Heat absorption: The heat transfer equation for the roplet heat-up stage may be written as: 4 3 T Nu π ( / ) ρpcp = π (Tsm T ) (8) 3 t The correlation for the Nusselt number in the case of spherical roplets moving through the air is taen in the form [6]: / / 3 Nu =.+.6 Re Pr (9) It is seen from this expression that the Nusselt number is a rather wea function of the particle Reynols number (Nu ~ Re / ) an may be represente with sufficient accuracy by its constant average value. roplet temperature at any moment uring the heat-up stage is obtaine by straightforwar integration of equation (8): 6Nu ρ C () t T = exp t (T T ) T sm p p () Expression () allows an estimation of the egree of a single roplet heat-up to be mae. The total heat absorption by the spray can be calculate provie the roplet size istribution in the spray is nown. Taing into account the temperature change along the roplet trajectory (), the heat loss rate to the spray from the layer, which has the temperature T sm is given by: 6Nu Q& * spr = ψ& c f (s) p exp t pc p ρ sm () () s s ( T T ) where ψ& is the water ischarge rate, f(s) is the probability ensity function for the roplet size istribution. The time t * (s) correspons to the en of the heat-up stage for a roplet of a given size. At this time, either roplet aporation begins, or the roplet exits from the smoe layer... Evaporation stage roplet ynamics: Generally, water roplets aporate as they move through the hot gas. Once the roplet is injecte into the hot layer, it will experience a perio of unsteay aporation, followe by a steay-state aporation regime, characterise by approximately constant surface roplet temperature ( wet bulb temperature ) an constant surface regression rate [5]. At high environment temperatures, wet bulb temperature approaches boiling temperature, so it may be reasonably assume that the steay-state aporation occurs at the surface boiling temperature. Equations of motion () apply for both heat-up an aporation stages. They can be re-written in the equivalent form: 3ρC t 4ρ p V 3ρC V + g t 4ρ p () 97
4 International Journal on Engineering Performance-Base Fire Coes uring the heat-up phase, the roplet iameter oes not change, an the solution is the same as for the set of equations (), that is given by the formulae (3),(5) an (6). uring the steay-state aporation, hower, the equations of motion () must be solve with the roplet iameter varying as a function of time. This function is taen as the theoretically an experimentally observe - law, which implies (t) = (t) t (3) The parameter is the so-calle aporation constant, an may be taen in the form [5]: ( T T ) 8λ g C g sm w = + ln (4) ρ pcg Hfg With the iameter regression rate (4), the equations () can still be integrate analytically. The solution can be written in the following form: 3ρC = ρp V = t + g α t 3 α t t + K where the parameters K an α are given by: K = α g + 3 ρ p α = + 3ρC ( α ) V (5) (6) ( an V are the roplet velocity components at the en of the heat-up stage, an is the initial roplet iameter). Heat absorption: For calculation of the total heat loss to the spray we note that the heat loss ue to aporation of a single roplet is: Q p ( T T ) + H fg m w = m C (7) where the roplet mass loss at a given time is calculate as: ( ) t 3/ 3 m = πρp (8) 6 The total heat loss rate to the spray may then be obtaine as: Q& spr = ψ& C p ( ) ** T T + H t (s) f (s) s w..3 Test calculations fg s 3 / (9) For the purpose of calculations, the roplet size istribution function of the spray is iscretize, an the roplets, therefore, are ivie between finite number of groups. Each group contains roplets of the same iameter. Calculation proceure is esigne in such a way that the equations for both heat-up an aporation stages are solve consecutively along roplet trajectory for one representative roplet from each of the groups. In general, each roplet unergoes both (heat-up an aporation) stages, an the total heat absorption by the roplet epens on its resience time in the smoe layer. The resience time for any particular roplet is etermine using the following equation: * t V ** t () s s V () s s = + () where V (s) an V (s) are taen from the solutions (5) an (5), respectively, is the smoe layer thicness. The resience time t res is equal to: t res = t * + t ** () where t * an t ** are the times for the heat-up an aporation stages, respectively. Equation () oes not have solution if roplet aporates completely in the layer. In that case, the resience time is taen as t * + t where t is the time for complete aporation of the roplet: t = () The total heat absorption rate is obtaine by multiplying the heat absorption rate for a single roplet by a number of roplets in the respective group, an summing up over all groups in the iscretize spray istribution. 98
5 International Journal on Engineering Performance-Base Fire Coes The typical calculate roplet trajectories are presente in Figs. to 4. The ifference in the roplet behavior can be observe for the ifferent temperatures of the smoe layer, ranging from 3 to 8 K. The trajectories are presente for ifferent representative roplet iameters an ifferent initial velocities. The smoe layer thicness is taen as m in Figs. to 4. Generally, the roplets below 5 µm are affecte by the increase temperature of the smoe. Compare to the ambient conitions (3 K, Fig. 4), fine roplets of to µm in iameter aporate quite quicly in the smoe layer with the temperatures above 5 K. This may be seen from Figs. an 3. A more important parameter for fire engineering applications is the rate of heat absorption by the whole spray. In orer to valiate the performance of the analytical moel, the results are compare to the CF calculations, which were performe solving exact equations for the two-phase (ropletgas) flow [7]. The results are presente for the AM4 spray nozzle. The roplet size istribution for such a nozzle [8] is shown in Fig. 5. In this case, the meium roplet iameter in the spray is of the orer of 3 µm. Heat absorption rate is plotte in Fig. 6 as a function of the hot layer temperature. It appears that the agreement between the CF an the simplifie moel is quite goo. Note that the heat absorption rate approaches the plateau at high smoe layer temperatures. This plateau correspons to the complete aporation of the spray. Also shown in Fig. 6 is a convective part of the total absorption rate. As expecte for water mist spray, major contribution to heat absorption is from vaporization process. X / m Y / m roplet iameters an initial conitions: : = 3 µm; = ms - ; V = ms - : = 5 µm; = ms - ; V = ms - : = µm; = ms - ; V = ms - : = µm; = ms - ; V = 4 ms - Fig. : Selecte roplet trajectories in a 8 K smoe layer (layer epth is m) Y / m X / m roplet iameters an initial conitions: : = 3 µm; = ms - ; V = ms - : = 5 µm; = ms - ; V = ms - : = µm; = ms - ; V = ms - : = µm; = ms - ; V = 4 ms - Fig. 3: Selecte roplet trajectories in a 5 K smoe layer (layer epth is m) 99
6 International Journal on Engineering Performance-Base Fire Coes X / m Y / m roplet iameters an initial conitions: : = 3 µm; = ms - ; V = ms - : = 5 µm; = ms - ; V = ms - : = µm; = ms - ; V = ms - : = µm; = ms - ; V = 4 ms - Fig. 4: Selecte roplet trajectories in a 3 K smoe layer (layer epth is m). Number fraction of roplets roplet iameter /mm Fig. 5: roplet size istribution for the AM4 nozzle [8] Total absorption rate / W Smoe layer temperature /K present analytical moel _ CF preictions [7]. convective absorption rate, as preicte by analytical moel Fig. 6: Heat absorption rate preictions for the AM4 nozzle (layer epth is m; water ischarge rate is equal to ψ& = 8 liters/min)
7 International Journal on Engineering Performance-Base Fire Coes The presente moel can also be use to preict temperature history of the smoe layer in compartments. This is accomplishe solving the simple energy conservation equation for the smoe layer: M sm c pg T t sm spr ( T ) Q& (3) sm where M sm is the mass of the smoe layer. For the solution of equation (3), the rate of heat absorption Q& spr ( T sm ) must be upate as the solution progresses using the escribe calculation proceure (equations () to ()). Typical results are presente in Figs. 7 an 8. The first of these figures plots smoe layer temperature versus time for ifferent water ischarge rates. The water ischarge rate of ψ& = 8 liters/min correspons to the nominal operating conitions for the AM4 nozzle. The other curves in Fig. 7 correspon to the reuce water ischarge rates. At the initial stage after sprinler activation, the temperature ecrease rate is roughly linearly proportional to the water ischarge rate. All curves in Fig. 7 approach the ambient temperature lel (3 K) at infinitely long times. The effect of compartment size on temperature ecrease rate can be estimate from Fig. 8. Sample results are presente for the three ifferent floor areas of the room. The water nozzle operates at its nominal lel of ψ& = 8 liters/min. The rooms with smaller floor area emonstrate higher temperature rop rate, until the temperature lels off to the ambient lel. Typically, the thermal time constant of the smoe layer is of the orer of min. Finally, typical roplet heat-up times (time require for the roplet to heat to boiling temperature) are presente in Fig. 9. The heat-up times are presente for ifferent smoe layer temperatures, an for seral roplet sizes, which contain most of the AM4 spray volume. 3. CONCLSIONS A simple analytical solution has been presente for water roplet motion an aporation rates in a hot smoe layer. The moel allows calculation of the heat absorption rates for water sprays an smoe layer temperature histories to be performe. As a calculation test, the results are presente for iniviual representative roplet trajectories in hot layers of ifferent temperatures. The calculations of the total heat absorption rate by the spray is also presente, an a goo agreement is foun between the simplifie moel an the results from comprehensive CF simulations. The agreement between the analytical an CF moeling emonstrates that the present moel can be applie for fire engineering calculations as a quic an economical alternative to full CF simulations. Smoe layer temperature /K : ψ& = 8 liters/min : ψ& = liters/min 3 : ψ& = 4 liters/min 4 : ψ& = 8 liters/min Time / s Fig. 7: Smoe layer temperature histories for ifferent water ischarge rates (roplet size istribution correspons to the AM4 nozzle; layer epth is m)
8 International Journal on Engineering Performance-Base Fire Coes Smoe layer temperature /K Time / s : 5 m x 5 m compartment : 4 m x 4 m compartment 3 : 3 m x 3 m compartment Fig. 8: Effect of the compartment floor area on the smoe temperature rop rate (roplet size istribution correspons to the AM4 nozzle; layer epth is m; water ischarge rate is equal to ψ& = 8 liters/min ) Heat-up time / s.4. roplet iameters:. : = 3 µm.8 : = 5 µm.6 : = 36 µm Smoe layer temperature / K Fig. 9: Heat-up time for various roplets NOMENCLATRE c p specific heat C rag coefficient roplet iameter f (s) pf for roplet size istribution H fg latent heat of water vaporisation m roplet mass Nu Nusselt number Q & rate of heat absorption Re Reynols number t time T temperature tangential velocity component of the roplet V vertical velocity component of the roplet y coorinate normal to the ceiling Gree symbols ϑ ρ ψ& smoe layer thicness imensionless temperature ensity water ischarge rate Subscripts initial g gas w wet bulb temperature p particle spr spray sm smoe
9 International Journal on Engineering Performance-Base Fire Coes REFERENCES. W.K. Chow an Y.L. Cheung, Simulation of sprinler - Hot layer interaction using a fiel moel, Fire an Materials, Vol. 8, pp (994).. V. Novozhilov,.J.E. Harvie, A.R. Green an J.H. Kent, A computational flui ynamic moel of fire burning rate an extinction by water sprinler, Combustion Science an Technology, Vol. 3, pp (997). 3. V. Novozhilov,.J.E. Harvie, J.H. Kent, V.B. Apte an. Pearson, A computational flui ynamics stuy of woo fire extinguishment by water sprinler, Fire Safety Journal, Vol. 9, pp (997). 4. W.K. Chow an A.C. Tang, Experimental stuies on sprinler water spray - Smoe layer interaction, Journal of Applie Fire Science, Vol. 4, No. 3, pp (994-95). 5. A.H. Lefebvre, Atomization an sprays, Hemisphere, New Yor (989). 6. G.M. Faeth, Evaporation an combustion of sprays, Progress in Energy an Combustion Science, Vol. 9, pp. -76 (983) 7. K. Nyanina an V. Novozhilov, Preiction of fire environment uner water mist suppression conitions by two-phase zone moel, Journal of Applie Fire Science, Vol., No. 3, pp. 93-4, (-). 8. R. Thomas, Low pressure water mist performance in open area protection, Fire Australia, pp. - (997). 9. Grinnell Fire Protection Proucts & Equipment, Catalog GPC 4.97, SA (997)
To understand how scrubbers work, we must first define some terms.
SRUBBERS FOR PARTIE OETION Backgroun To unerstan how scrubbers work, we must first efine some terms. Single roplet efficiency, η, is similar to single fiber efficiency. It is the fraction of particles
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationStudy of Cooling System with Water Mist Sprayers: Fundamental Examination of Particle Size Distribution and Cooling Effects
Buil Simul (28) 1: 95 11 DOI 1.17/s12273-8-8115-y RESEARCH ARTICLE Stuy of Cooling System with Water Mist Sprayers: Funamental Examination of Particle Size Distribution an Cooling Effects Hieki Yamaa 1,
More informationEvaporating droplets tracking by holographic high speed video in turbulent flow
Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es
More informationNUMERICAL SIMULATION OF HEAT AND MASS TRANSFER, MIXTURE FORMATION IN COMBUSTION CHAMBER OF GAS TURBINE ENGINE
MATHEMATICA MONTISNIGRI Vol XL (2017) MATHEMATICAL MODELING NUMERICAL SIMULATION OF HEAT AND MASS TRANSFER, MIXTURE FORMATION IN COMBUSTION CHAMBER OF GAS TURBINE ENGINE A.A. SVIRIDENKOV, P.D. TOKTALIEV
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More informationCONTROL CHARTS FOR VARIABLES
UNIT CONTOL CHATS FO VAIABLES Structure.1 Introuction Objectives. Control Chart Technique.3 Control Charts for Variables.4 Control Chart for Mean(-Chart).5 ange Chart (-Chart).6 Stanar Deviation Chart
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationSection 2.7 Derivatives of powers of functions
Section 2.7 Derivatives of powers of functions (3/19/08) Overview: In this section we iscuss the Chain Rule formula for the erivatives of composite functions that are forme by taking powers of other functions.
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationA Molten Solid Approach for Simulating Urea-Water Solution Droplet Depletion
ILASS Americas 7th Annual Conference on Liqui Atomization an Spray Systems, Raleigh, NC, May 015 A Molten Soli Approach for Simulating Urea-Water Solution Droplet Depletion Shaoping Quan* 1, Mingjie Wang
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationLinear analysis of a natural circulation driven supercritical water loop
TU Delft Bachelor Thesis Linear analysis of a natural circulation riven supercritical water loop D J van er Ham 4285816 supervise by Dr. Ir. M. Rohe July 3, 216 Nomenclature Symbol Units Description A
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More information3.2 Shot peening - modeling 3 PROCEEDINGS
3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationSIMULATION OF POROUS MEDIUM COMBUSTION IN ENGINES
SIMULATION OF POROUS MEDIUM COMBUSTION IN ENGINES Jan Macek, Miloš Polášek Czech Technical University in Prague, Josef Božek Research Center Introuction Improvement of emissions from reciprocating internal
More informationSPE Copyright 1999, Society of Petroleum Engineers Inc.
SPE 664 Effect of Flow Through a Choke Valve on Emulsion Stability M.J. van er Zane, SPE, K.R. van Heuven, J.H. Muntinga, SPE, an W.M.G.T. van en Broek, SPE, Delft University of Technology Copyright 1999,
More informationThe Press-Schechter mass function
The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationEvaluation of Column Breakpoint and Trajectory for a Plain Liquid Jet Injected into a Crossflow
ILASS Americas, 1 st Annual Conference on Liqui Atomization an Spray Systems, Orlano, Floria, May 008 Evaluation of Column Breakpoint an Trajectory for a Plain Liqui Jet Injecte into a Crossflow S.M. Thawley,
More informationThe influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section
IOP Conference Series: Materials Science an Engineering PAPER OPEN ACCESS The influence of the equivalent hyraulic iameter on the pressure rop preiction of annular test section To cite this article: A
More informationSection 2.1 The Derivative and the Tangent Line Problem
Chapter 2 Differentiation Course Number Section 2.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan
More informationA new identification method of the supply hole discharge coefficient of gas bearings
Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico
More informationAspects of Numerical Modelling of Flash-Boiling Fuel Sprays
Copyright 15 SAE International 15-4-463 Aspects of Numerical Moelling of Flash-Boiling Fuel Sprays Christopher Price, Arash Hamzehloo an Pavlos Aleiferis University College Lonon, UK Davi Richarson Jaguar
More informationDusty Plasma Void Dynamics in Unmoving and Moving Flows
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationEVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD
137 THERMAL SCIENCE, Year 1, Vol. 16, No. 5, pp. 137-1376 EVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD by Fu-Jun GAN a an Jian-Zhong LIN a,b * a Department
More informationABSTRACT. Analysis of the Initial Spray from Canonical Fire Suppression Nozzles. Professor André Marshall, Department of Fire Protection Engineering
ABSTRACT Title of Document: Analysis of the Initial Spray from Canonical Fire Suppression Nozzles Ning Ren, Master of Science, 007 Directe By: Professor Anré Marshall, Department of Fire Protection Engineering
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationD DAVID PUBLISHING. Q Energy flux (W) Nu Nusselt number (hd/k) Nomenclature. Greek Symbols
Journal of Energy an Power Engineering 9 (2015) 222-231 oi: 1017265/1934-8975/201502012 D DAVID PUBLISHING Effect of Injection Parameters on Evaporation an Thermolysis Characteristics of UWS (Urea-Water-Solution)
More informationSIMULATION OF DIRECT CONTACT CONDENSATION OF STEAM JETS SUBMERGED IN SUBCOOLED WATER BY MEANS OF A ONE-DIMENSIONAL TWO-FLUID MODEL
HEFAT014 10 th International Conference on Heat Transfer, Flui Mechanics an Thermoynamics 14 16 July 014 Orlano, Floria SIMULATION OF DIRECT CONTACT CONDENSATION OF STEAM JETS SUBMERGED IN SUBCOOLED WATER
More informationPractical application of numerical modelling to overbank flows in a compound river channel
Practical application of numerical moelling to overbank flows in a compoun river channel Moreta, P.M. Lecturer, School of Computing an Engineering, University of West Lonon, Ealing, W5 5RF Lonon, UK Rotimi
More informationExercise 4 - Hydraulic Systems
Exercise 4 - Hyraulic Systems 4.1 Hyraulic Systems Hyraulic systems are, in general, escribe by the Navier-Stokes equations as you might have learne in flui ynamics courses. In orer to simplify the moeling
More informationModel for Dopant and Impurity Segregation During Vapor Phase Growth
Mat. Res. Soc. Symp. Proc. Vol. 648, P3.11.1-7 2001 Materials Research Society Moel for Dopant an Impurity Segregation During Vapor Phase Growth Craig B. Arnol an Michael J. Aziz Division of Engineering
More informationDefine each term or concept.
Chapter Differentiation Course Number Section.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationMULTISCALE FRICTION MODELING FOR SHEET METAL FORMING
MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of
More informationStudy on aero-acoustic structural interactions in fan-ducted system
Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,
More informationInfluence the Nozzle Shape on Local Heat Transfer in Impinging Jet
Issue 6, Volume 6, 12 Influence the Nozzle Shape on Local Heat Transfer in Impinging Jet M. Attalla an M. S. Ahme 1 Abstract The local Nusselt number istributions of circular nozzle on a heate flat plate
More informationApproaches for Predicting Collection Efficiency of Fibrous Filters
Volume 5, Issue, Summer006 Approaches for Preicting Collection Efficiency of Fibrous Filters Q. Wang, B. Maze, H. Vahei Tafreshi, an B. Poureyhimi Nonwovens Cooperative esearch Center, North Carolina State
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationwater adding dye partial mixing homogenization time
iffusion iffusion is a process of mass transport that involves the movement of one atomic species into another. It occurs by ranom atomic jumps from one position to another an takes place in the gaseous,
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationChapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies
hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation
More information18 EVEN MORE CALCULUS
8 EVEN MORE CALCULUS Chapter 8 Even More Calculus Objectives After stuing this chapter you shoul be able to ifferentiate an integrate basic trigonometric functions; unerstan how to calculate rates of change;
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationCHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić
Sensors an Actuators A: Physical 168 (211) 149-154 oi: 1.116/j.sna.211..7 211 Elsevier B.V. CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS Jože Kutin *, Ivan Bajsić Laboratory of
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationAn inductance lookup table application for analysis of reluctance stepper motor model
ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR
More informationSparse Reconstruction of Systems of Ordinary Differential Equations
Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA
More informationParametrization of non-convective condensation processes May 1987
Parametrization of non-convective conensation processes May 987 By M. Tietke European Centre for Meium-Range Weather Forecasts Table of contents. Thermoynamics of moist air 2. Clou physical processes 2.
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationSimulation of jet-flow solid fraction during spray forming
International Journal of Minerals, Metallurgy an Materials Volume 24, Number 6, June 2017, Page 657 DOI: 10.1007/s12613-017-1448-8 Simulation of jet-flow soli fraction uring spray forming Zi-qiang Pi 1),
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationBasic Differentiation Rules and Rates of Change. The Constant Rule
460_00.q //04 4:04 PM Page 07 SECTION. Basic Differentiation Rules an Rates of Change 07 Section. The slope of a horizontal line is 0. Basic Differentiation Rules an Rates of Change Fin the erivative of
More informationTutorial Test 5 2D welding robot
Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationarxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003
Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep
More informationA USER-FRIENDLY TOOL FOR SIMULATING THE TIME-DEPENDENT FIELD DISTRIBUTION IN PE INSULATION ON THE BASIS OF A PHYSICAL APPROACH
A USER-FRIENLY TOOL FOR SIMULATING THE TIME-EPENENT FIEL ISTRIBUTION IN PE INSULATION ON THE BASIS OF A PHYSICAL APPROACH Séverine LE ROY, Laboratoire Plasma et Conversion Energie, (France), severine.leroy@laplace.univ-tlse.fr
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationModelling of Multiphase Flow in Pressure Swirl Atomizer
International Journal of Engineering an Technology olume 6 No.5, May, 216 Moelling of Multiphase Flow in Pressure wirl Atomizer.K. Ameorme chool of Mechanical Engineering, Institute of Thermofluis an Combustion,
More informationAn analytical investigation into filmwise condensation on a horizontal tube in a porous medium with suction at the tube surface
Heat Mass Transfer (29) 45:355 361 DOI 1.17/s231-8-436-y ORIGINAL An analytical investigation into filmwise conensation on a horizontal tube in a porous meium with suction at the tube surface Tong Bou
More informationPredictive Control of a Laboratory Time Delay Process Experiment
Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationIPMSM Inductances Calculation Using FEA
X International Symposium on Inustrial Electronics INDEL 24, Banja Luka, November 68, 24 IPMSM Inuctances Calculation Using FEA Dejan G. Jerkan, Marko A. Gecić an Darko P. Marčetić Department for Power,
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationA 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS
A 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS Peifeng Ma 1 an Ole Secher Masen Accurate preiction of current velocity an bottom shear stress, which both can be significantly influence by
More informationNUMERICAL STUDY OF THERMAL RADIATIONS AND THERMAL STRATIFICATION MECHANISMS IN MHD CASSON FLUID FLOW. and Sardar Muhammad BILAL c
NUMERICAL STUDY OF THERMAL RADIATIONS AND THERMAL STRATIFICATION MECHANISMS IN MHD CASSON FLUID FLOW Khalil Ur REHMAN b c * Noor Ul SABA b Iffat ZEHRA c Muhamma Yousaf MALIK ab an Sarar Muhamma BILAL c
More informationMomentum and Energy. Chapter Conservation Principles
Chapter 2 Momentum an Energy In this chapter we present some funamental results of continuum mechanics. The formulation is base on the principles of conservation of mass, momentum, angular momentum, an
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationRFSS: Lecture 4 Alpha Decay
RFSS: Lecture 4 Alpha Decay Reaings Nuclear an Raiochemistry: Chapter 3 Moern Nuclear Chemistry: Chapter 7 Energetics of Alpha Decay Geiger Nuttall base theory Theory of Alpha Decay Hinrance Factors Different
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationModeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N.
Submitte to the Journal of Rheology Moeling the effects of polyispersity on the viscosity of noncolloial har sphere suspensions Paul M. Mwasame, Norman J. Wagner, Antony N. Beris a) epartment of Chemical
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationAIR BUBBLE ENTRAINMENT IN HYDRAULIC JUMPS: PHYSICAL MODELING AND SCALE EFFECTS
AIR BUBBLE ENTRAINMENT IN HYDRAULIC JUMPS: PHYSICAL MODELING AND SCALE EFFECTS Hubert CHANSON Professor in Civil Engineering The University of Queenslan Brisbane QLD 4072 Australia Ph.: (6 7) 3365 463
More informationu t v t v t c a u t b a v t u t v t b a
Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying
More informationA Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion
Open Journal o Flui Dynamics, 5, 5, 5-57 Publishe Online March 5 in SciRes. http://www.scirp.org/journal/oj http://x.oi.org/.46/oj.5.57 A Short Note on Sel-Similar Solution to Unconine Flow in an Aquier
More information5-4 Electrostatic Boundary Value Problems
11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions
More informationEssential Considerations for Buckling Analysis
Worlwie Aerospace Conference an Technology Showcase, Toulouse, France, Sept. 24-26, 2001 Essential Consierations for Buckling Analysis 2001-120 Sang H. Lee MSC.Software Corporation, 2 MacArthur Place,
More informationCOUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND
More informationV = Flow velocity, ft/sec
1 Drag Coefficient Preiction Chapter 1 The ieal force acting on a surface positione perpenicular to the airflow is equal to a ynamic pressure, enote by q, times the area of that surface. Dynamic pressure
More information