A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS

Size: px
Start display at page:

Download "A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS"

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.

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 information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow 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 information

Schrödinger s equation.

Schrö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

'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 information

Study of Cooling System with Water Mist Sprayers: Fundamental Examination of Particle Size Distribution and Cooling Effects

Study 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 information

Evaporating droplets tracking by holographic high speed video in turbulent flow

Evaporating 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 information

NUMERICAL SIMULATION OF HEAT AND MASS TRANSFER, MIXTURE FORMATION IN COMBUSTION CHAMBER OF GAS TURBINE ENGINE

NUMERICAL 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 information

Applications of First Order Equations

Applications 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 information

CONTROL CHARTS FOR VARIABLES

CONTROL 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 information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization 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 information

Section 2.7 Derivatives of powers of functions

Section 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 information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE 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 information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 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 information

A Molten Solid Approach for Simulating Urea-Water Solution Droplet Depletion

A 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 information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal 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 information

Linear analysis of a natural circulation driven supercritical water loop

Linear 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 information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 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 information

3.2 Shot peening - modeling 3 PROCEEDINGS

3.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 information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 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 information

SIMULATION OF POROUS MEDIUM COMBUSTION IN ENGINES

SIMULATION 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 information

SPE Copyright 1999, Society of Petroleum Engineers Inc.

SPE 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 information

The Press-Schechter mass function

The 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 information

The total derivative. Chapter Lagrangian and Eulerian approaches

The 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 information

Introduction to the Vlasov-Poisson system

Introduction 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 information

Evaluation of Column Breakpoint and Trajectory for a Plain Liquid Jet Injected into a Crossflow

Evaluation 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 information

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section

The 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 information

Section 2.1 The Derivative and the Tangent Line Problem

Section 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 information

A new identification method of the supply hole discharge coefficient of gas bearings

A 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 information

Aspects of Numerical Modelling of Flash-Boiling Fuel Sprays

Aspects 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 information

Dusty Plasma Void Dynamics in Unmoving and Moving Flows

Dusty 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 information

The effect of nonvertical shear on turbulence in a stably stratified medium

The 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 information

EVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD

EVOLUTION 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 information

ABSTRACT. Analysis of the Initial Spray from Canonical Fire Suppression Nozzles. Professor André Marshall, Department of Fire Protection Engineering

ABSTRACT. 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 information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic 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 information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE 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 information

D DAVID PUBLISHING. Q Energy flux (W) Nu Nusselt number (hd/k) Nomenclature. Greek Symbols

D 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 information

SIMULATION OF DIRECT CONTACT CONDENSATION OF STEAM JETS SUBMERGED IN SUBCOOLED WATER BY MEANS OF A ONE-DIMENSIONAL TWO-FLUID MODEL

SIMULATION 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 information

Practical application of numerical modelling to overbank flows in a compound river channel

Practical 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 information

Exercise 4 - Hydraulic Systems

Exercise 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 information

Model for Dopant and Impurity Segregation During Vapor Phase Growth

Model 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 information

Define each term or concept.

Define 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 information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 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 information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic 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 information

MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING

MULTISCALE 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 information

Study on aero-acoustic structural interactions in fan-ducted system

Study 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 information

Influence the Nozzle Shape on Local Heat Transfer in Impinging Jet

Influence 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 information

Approaches for Predicting Collection Efficiency of Fibrous Filters

Approaches 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 information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d 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 information

water adding dye partial mixing homogenization time

water 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 information

Lagrangian and Hamiltonian Mechanics

Lagrangian 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 information

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies

Chapter-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 information

18 EVEN MORE CALCULUS

18 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 information

The Exact Form and General Integrating Factors

The 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 information

Table of Common Derivatives By David Abraham

Table 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 information

CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić

CHARACTERISTICS 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 information

Homework 7 Due 18 November at 6:00 pm

Homework 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 information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 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 information

An inductance lookup table application for analysis of reluctance stepper motor model

An 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 information

Sparse Reconstruction of Systems of Ordinary Differential Equations

Sparse 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 information

Parametrization of non-convective condensation processes May 1987

Parametrization 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 information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS 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 information

Simulation of jet-flow solid fraction during spray forming

Simulation 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 information

The Principle of Least Action

The 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 information

Basic Differentiation Rules and Rates of Change. The Constant Rule

Basic 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 information

Tutorial Test 5 2D welding robot

Tutorial 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 information

12.11 Laplace s Equation in Cylindrical and

12.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 information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv: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 information

A USER-FRIENDLY TOOL FOR SIMULATING THE TIME-DEPENDENT FIELD DISTRIBUTION IN PE INSULATION ON THE BASIS OF A PHYSICAL APPROACH

A 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 information

Linear First-Order Equations

Linear 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 information

Implicit Differentiation

Implicit 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 information

Quantum Mechanics in Three Dimensions

Quantum 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 information

Modelling of Multiphase Flow in Pressure Swirl Atomizer

Modelling 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 information

An analytical investigation into filmwise condensation on a horizontal tube in a porous medium with suction at the tube surface

An 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 information

Predictive Control of a Laboratory Time Delay Process Experiment

Predictive 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 information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized 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 information

IPMSM Inductances Calculation Using FEA

IPMSM 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 information

Chapter 6: Energy-Momentum Tensors

Chapter 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 information

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

Application 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 information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles 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

θ 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 information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture 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 information

Efficient Macro-Micro Scale Coupled Modeling of Batteries

Efficient 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 information

Chapter 2 Lagrangian Modeling

Chapter 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 information

A 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS

A 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 information

NUMERICAL 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. 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 information

Momentum and Energy. Chapter Conservation Principles

Momentum 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 information

BEYOND 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 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 information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - 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 information

RFSS: Lecture 4 Alpha Decay

RFSS: 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 information

Lie symmetry and Mei conservation law of continuum system

Lie 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 information

Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N.

Modeling 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 information

Stable and compact finite difference schemes

Stable 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 information

Separation of Variables

Separation 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 information

AIR BUBBLE ENTRAINMENT IN HYDRAULIC JUMPS: PHYSICAL MODELING AND SCALE EFFECTS

AIR 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 information

u t v t v t c a u t b a v t u t v t b a

u 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 information

A Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion

A 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 information

5-4 Electrostatic Boundary Value Problems

5-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 information

Essential Considerations for Buckling Analysis

Essential 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 information

COUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS

COUPLING 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 information

V = Flow velocity, ft/sec

V = 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