Rainout Assessment after a Release of Hazardous Liquefied Gas: Effect of Physical and Storage Parameters

Size: px
Start display at page:

Download "Rainout Assessment after a Release of Hazardous Liquefied Gas: Effect of Physical and Storage Parameters"

Transcription

1 Rainou Assessmen afer a Release of Haardous Liquefied Gas: Effec of Physical and Sorage Parameers Valenina Busini 1, Renao Roa 1*, Federica i Vio 2, Sabaino iali, Robero Fiore 3 1 Poliecnico di Milano ip. Chimica, Maeriali e Ingegneria Chimica G. Naa via Mancinelli 7, 2131 Milano, Ialy, renao.roa@polimi.i 2 Tecnimon SpA esign HSE ep. 3 Saipem SpA - B.U. Onshore Snamprogei Cenre of Excellence - Loss Prevenion & Environmen ep. In process indusry i is common pracice o sore and handle dangerous subsances, which are gaseous a amospheric emperaure, as pressuried or refrigeraed liquids. In case of accidenal release, a uid in hese condiions gives rise o a wo-phase discharge. In such condiions, he liquid breaks up in an aerosol cloud, whose behavior in amosphere affecs dispersion disances. The physical phenomena occurring afer he release o he amosphere are mainly hree: expansion o ambien pressure, liquid aomiaion (break-up) and rainou. In paricular, he rainou fracion is a crucial parameer in quaniaive risk analysis and herefore models capable of esimaing wheher i occurs or no are of paramoun imporance. This sudy aims a emphasiing he inuence of some physical parameers and sorage relaed parameers on he rainou fracion, aking ino accoun coninuous and saionary releases of ammonia and LNG (Liquefied Naural Gas). These represen wo differen classes of subsances: gas liquefied by pressure (superheaed) and gas liquefied by refrigeraion (subcooled), which undergo o differen je break-up mechanisms. The resuls obained for hese wo compounds can hen be exended o he whole range of pressurised liquids/liquefied gases. A simple model has been developed o evaluae rainou occurrence, he simple approach here proposed (which has been validaed comparing is resuls o hose of oher simulaion models as well as o some available experimenal daa) allows o invesigae in a simple and effecive way he imporance of he differen physical and sorage parameers involved during rainou. 1. Inroducion angerous maerials ha are gaseous a ambien emperaure are commonly sored as pressuried or refrigeraed liquids. In case of accidenal release, hese liquids give rise o a wo-phase discharge. In such condiions, he liquid breaks up in an aerosol cloud, which behavior in amosphere affecs dispersion disances. The physical phenomena occurring afer he release o he amosphere are mainly hree: expansion o ambien pressure, liquid aomiaion (break-up) and rainou.

2 In paricular, he larges les can lead o he formaion of an evaporaing pool, while he smalles ones will say ino he gas cloud. Undersanding wheher a specific liquid release will lead o a pool formaion or no is of paramoun imporance o define he inle condiions for he cloud dispersion modeling leading o he safey disance esimaion. This is quie cumbersome since he physical phenomena occurring afer he release of a liquefied gas (namely, je expansion up o he ambien pressure, mechanical and/or ashing liquid break-up, and rainou) are no easy o model in deail. The rainou fracion is a parameer of crucial imporance in Quaniaive Risk Analysis (QRA) and herefore lumped parameers mahemaical models have been developed for mechanical and/or ashing liquid break-up, as well as for rainou (Cleary e al., 27; Wilox e al. 27). However, such models, while requiring dedicaed compuer codes (e.g., Wilox and Hol, 2) sill reain many approximaions and heir agreemen wih he available experimenal daa is quesionable. I follows ha for a firs-aemp esimaion of wheher rainou will occur or no (which is he main informaion required for seing up he subsequen simulaion of he amospheric dispersion of he ammable cloud) simpler correlaions would be desirable. Among all he pressuried and refrigeraed liquids in his paper he ammonia and he Liquefied Naural Gas (LNG) has been considered and he available experimenal daa of heir rainou have been used o validae a shor cu mehod for predicing rainou fracion. 2. Shor-cu Rainou Correlaion A simple correlaion o be used for predicing wheher rainou is expeced or no in a given siuaion should be able o accoun for he inuence of he main physical and sorage-relaed parameers inuencing he liquid break up and he subsequen les evaporaion (e.g., laen hea of vaporiaion, vapor pressure, liquid densiy and surface ension, -sie, wind condiion, ambien emperaure, air enrainmen). A simple correlaion has been developed o evaluae rainou occurrence, saring from he observaion ha rainou occurs when he uid released afer a loss of conainmen breaks ou ino les which ake more ime o evaporae (i.e. he evaporaion ime evap ) han o fall on he ground (i.e. igh ime, ). This consideraion leads o he following crieria for rainou exisence: evap evap > 1 < 1 ηrainou ηrainou > = (1) Boh he characerisic evaporaion ime and he characerisic igh ime depend on les diameer. In urn, his characerisic dimension is deermined by he mechanism ha deermines he liquid break up: shear sress (mechanic or aerodynamic breakage) or

3 liquid ash (hermodynamic breakage). The dominaing mechanism depends on he physical characerisics of he liquid as well as on he sorage condiions (e Vaull and King, 1992; Bricard and Friedel, 1998; Wilox e al. 27), which herefore deermine he characerisic le dimension and in urn he wo characerisic imes. evap By comparing hese wo characerisic imes (ha is, when = 1) i is possible o idenify a hreshold le diameer, lim, for rainou: les larger han his hreshold will rainou, and he rainou fracion, η rainou, will be roughly proporional o he difference beween he le diameer and he hreshold diameer, as follow: ( lim ) η rainou = (2) The evaporaion ime for he le can be easily evaluaed inegraing he maerial balance and he energy balance for he evaporaing le: evap 1 ΔH 2 evap ρ liq = 4 ( + sup ) v Sh AB ρ ρ h Nu k ( T T ) (3) is Sauer mean le diameer, ΔH evap where released liquid evaporaion enhalpy, ρ liq released liquid densiy, Sh Sherwood number, AB maerial diffusiviy, ρ air densiy, ρ P v ( T ) MW sup = is he le superficial densiy (MW RT v molecular weigh and R perfec gas consan), h vapor enhalpy, Nu Nussel number, k hermal exchange coefficien, T le emperaure, and T air emperaure. Equaion 3 shows ha he evaporaion ime depends from he le diameer boh direcly and indirecly, for he presence of Sh and Nu. The igh ime of he le can be obained by he inegraion of he momenum equaion for he le along horional and verical direcions, applying he following assumpions: 24 laminar condiions: C =, Re du, consan verical velociy: =, d

4 The igh ime for he le is insead evaluaed as follow: = 2 g( ρ liq ρ ) 18μ + u, (4) Where is he release heigh, g graviy acceleraion, μ he dynamic viscosiy, and verical componen of wind velociy u. Afer having defined he evaporaion ime and he igh ime for he le i is possible o evaluae heir raio, assuming ha: he le velociy ( u ) is consan, he wind direcion is horional ( u, = ), he le emperaure ( T ) is consan, he release direcion is horional. Solving he simple relaion evap =, he limi le diameer is deermined and hen, knowing he iniial le diameer, also he rainou fracion can be esimaed. The iniial le diameer is evaluaed choosing he smalles one beween he wo calculaed considering boh mechanical and aerodynamic uid break-ou. The wo correlaions used o evaluae he diameer of les originaed by mechanical ( mec ) and aerodynamic ( aer ) break-ou are he following: σ Wecri mec = (5) 2 u ρ f u, where σ is he surface ension, he criical Weber number We cri,=12.5, and u f le velociy afer je expansion. = ln (6) 3 3 aer E p where E p [ P v ( T s ) P ] ν s + [ P s P v ( T s )] ν s = ΔH evap (7)

5 is he parial expansion energy, P v is he vapor pressure, T s he sorage emperaure of he released uid, P he pressure a he end of he expansion, P s is he sorage pressure, and ν he specific volume of he released uid a sorage condiions. s 3. Resuls This simple approach has been validaed comparing is resuls o hose of oher simulaion models as well as o some available experimenal daa; in paricular varying he heigh of he release, he sorage emperaure, and he sorage pressure. Moreover, i was applied considering only a mean diameer, bu i would define a limi diameer also in he case of a disribuion of le diameers (i.e., les larger han he lim rainou); of course in his case, he formula should be weighed on he disribuion funcion. For he sake of example, Figure 1 and Figure 2 compare he ypical resuls obained using his simple correlaion for LNG and ammonia respecively, wih he resuls of one of he more deailed models available (Wilox and Hol, 2). We can see ha he simple correlaion is able o represen quie correcly he rainou fracion, and in paricular i is able o foresee he presence or absence of rainou when some physical or operaing parameer (he sorage pressure in his case) changes Rainou fracion rople diameer [mm] (o-lim)/o PHAST rainou fracion rople diameer P L=1, Sorage Pressure [am] Figure 1 Rainou fracion and characerisic les diameer for a 1 m heigh release from a ank a consan emperaure as a funcion of he sorage pressure - LNG

6 Rainou fracion rople diameer [mm] (o-lim)/o PHAST rainou fracion rople diameer P L =6, Sorage Pressure [am] Figure 2 Rainou fracion and characerisic les diameer for a 1 m heigh release from a ank a consan emperaure as a funcion of he sorage pressure - ammonia References Bricard, P., Friedel, L., 1998, Two-phase je dispersion. J. Haardous Maerials, 59, 287. Cleary, V.M., Bowen, P.J., Wilox, H.W.M., 27, Flashing liquid jes and wo-phase le dispersion. I. Experimens for derivaion of le aomisaion correlaions. J. Haardous Maerials, 142, 786. e Vaull, G. E., King, J. A., 1992, Similariy scaling of le evaporaion and liquid rain-ou following he release of a superheaed ashing liquid o he environmen, For Presenaion a he 85h Annual Meeing, Air & Wase Managemen Assoc., Kansas Ciy, MO, June Wilox, H.W.M., Harper, M., Bowenb, P., Cleary, V.M., 27, Flashing liquid jes and wo-phase le dispersion. II. Comparison and validaion of le sie and rainou formulaions. J. of Haardous Maerials, 142, 797. Wilox, H.W.M., Hol, A., 2, Unified ispersion Model, in: Technical Reference Manual 6., NV, London.

Turbulence in Fluids. Plumes and Thermals. Benoit Cushman-Roisin Thayer School of Engineering Dartmouth College

Turbulence in Fluids. Plumes and Thermals. Benoit Cushman-Roisin Thayer School of Engineering Dartmouth College Turbulence in Fluids Plumes and Thermals enoi Cushman-Roisin Thayer School of Engineering Darmouh College Why do hese srucures behave he way hey do? How much mixing do hey accomplish? 1 Plumes Plumes are

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

ENVIRONMENTAL FLUID MECHANICS

ENVIRONMENTAL FLUID MECHANICS ENVIONMENTAL FLUID MECHANICS Plumes & Thermals Why do hese srucures behave he way hey do? How much mixing do hey accomplish? enoi Cushman-oisin Thayer School of Engineering Darmouh College hp://hayer.darmouh.edu/~cushman/books/efm/chap1.pdf

More information

5.2. The Natural Logarithm. Solution

5.2. The Natural Logarithm. Solution 5.2 The Naural Logarihm The number e is an irraional number, similar in naure o π. Is non-erminaing, non-repeaing value is e 2.718 281 828 59. Like π, e also occurs frequenly in naural phenomena. In fac,

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

Sub Module 2.6. Measurement of transient temperature

Sub Module 2.6. Measurement of transient temperature Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Scientific Herald of the Voronezh State University of Architecture and Civil Engineering. Construction and Architecture

Scientific Herald of the Voronezh State University of Architecture and Civil Engineering. Construction and Architecture Scienific Herald of he Voronezh Sae Universiy of Archiecure and Civil Engineering. Consrucion and Archiecure UDC 625.863.6:551.328 Voronezh Sae Universiy of Archiecure and Civil Engineering Ph. D. applican

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

Chapter 12: Velocity, acceleration, and forces

Chapter 12: Velocity, acceleration, and forces To Feel a Force Chaper Spring, Chaper : A. Saes of moion For moion on or near he surface of he earh, i is naural o measure moion wih respec o objecs fixed o he earh. The 4 hr. roaion of he earh has a measurable

More information

5 The fitting methods used in the normalization of DSD

5 The fitting methods used in the normalization of DSD The fiing mehods used in he normalizaion of DSD.1 Inroducion Sempere-Torres e al. 1994 presened a general formulaion for he DSD ha was able o reproduce and inerpre all previous sudies of DSD. The mehodology

More information

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

Comparing Means: t-tests for One Sample & Two Related Samples

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

Ch.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

Ch.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC Ch.. Group Work Unis Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Uni 2 Jusify wheher he following saemens are rue or false: a) Two sreamlines, corresponding o a same insan of ime, can never inersec

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES

THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES Kragujevac J. Sci. 3 () 7-4. UDC 53.5:536. 4 THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES Hazem A. Aia Dep. of Mahemaics, College of Science,King Saud Universiy

More information

Cosmic String Loop Distribution with a Gravitational Wave Cutoff

Cosmic String Loop Distribution with a Gravitational Wave Cutoff Ouline Inroducion Towards he Cosmological Aracor Cosmic Sring Loop Disribuion wih a Graviaional Wave Cuoff Larissa Lorenz Insiue of Mahemaics and Physics Cenre for Cosmology, Paricle Physics and Phenomenology

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Q.1 Define work and its unit?

Q.1 Define work and its unit? CHP # 6 ORK AND ENERGY Q.1 Define work and is uni? A. ORK I can be define as when we applied a force on a body and he body covers a disance in he direcion of force, hen we say ha work is done. I is a scalar

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

MECHANICAL PROPERTIES OF FLUIDS NCERT

MECHANICAL PROPERTIES OF FLUIDS NCERT Chaper Ten MECHANICAL PROPERTIES OF FLUIDS MCQ I 10.1 A all cylinder is filled wih iscous oil. A round pebble is dropped from he op wih zero iniial elociy. From he plo shown in Fig. 10.1, indicae he one

More information

Development of a new metrological model for measuring of the water surface evaporation Tovmach L. Tovmach Yr. Abstract Introduction

Development of a new metrological model for measuring of the water surface evaporation Tovmach L. Tovmach Yr. Abstract Introduction Developmen of a new merological model for measuring of he waer surface evaporaion Tovmach L. Tovmach Yr. Sae Hydrological Insiue 23 Second Line, 199053 S. Peersburg, Russian Federaion Telephone (812) 323

More information

A NUMERICAL STUDY OF THE EFFECT OF CONTACT ANGLE ON THE DYNAMICS OF A SINGLE BUBBLE DURING POOL BOILING

A NUMERICAL STUDY OF THE EFFECT OF CONTACT ANGLE ON THE DYNAMICS OF A SINGLE BUBBLE DURING POOL BOILING Proceedings of IMECE ASME Inernaional Mechanical Engineering Congress & Exposiion November 17-,, New Orleans, Louisiana IMECE-33876 A NUMERICAL STUDY OF THE EFFECT OF CONTACT ANGLE ON THE DYNAMICS OF A

More information

not to be republished NCERT MATHEMATICAL MODELLING Appendix 2 A.2.1 Introduction A.2.2 Why Mathematical Modelling?

not to be republished NCERT MATHEMATICAL MODELLING Appendix 2 A.2.1 Introduction A.2.2 Why Mathematical Modelling? 256 MATHEMATICS A.2.1 Inroducion In class XI, we have learn abou mahemaical modelling as an aemp o sudy some par (or form) of some real-life problems in mahemaical erms, i.e., he conversion of a physical

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion

More information

ASSESSMENT OF BUOYANCY-CORRECTED TURBULENCE MODELS FOR THERMAL PLUMES

ASSESSMENT OF BUOYANCY-CORRECTED TURBULENCE MODELS FOR THERMAL PLUMES Engineering Applicaions of Compuaional Fluid Mechanics Vol. 7, No., pp. 9 49 (1) ASSESSMENT OF BUOYANCY-CORRECTED TURBULENCE MODELS FOR THERMAL PLUMES Raesh Kumar and Anupam Dewan * * Deparmen of Applied

More information

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he

More information

1 Model equations and parameters

1 Model equations and parameters Supporing Informaion: Pre-combusion capure by PSA: Comparison of laboraory PSA Experimens and Simulaions (Indusrial Engineering Chemisry Research) J. Schell, N. Casas, D. Marx and M. Mazzoi ETH Zurich,

More information

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

More information

Non-uniform circular motion *

Non-uniform circular motion * OpenSax-CNX module: m14020 1 Non-uniform circular moion * Sunil Kumar Singh This work is produced by OpenSax-CNX and licensed under he Creaive Commons Aribuion License 2.0 Wha do we mean by non-uniform

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

5.1 - Logarithms and Their Properties

5.1 - Logarithms and Their Properties Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We

More information

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

More information

Verification of a CFD benchmark solution of transient low Mach number flows with Richardson extrapolation procedure 1

Verification of a CFD benchmark solution of transient low Mach number flows with Richardson extrapolation procedure 1 Verificaion of a CFD benchmark soluion of ransien low Mach number flows wih Richardson exrapolaion procedure S. Beneboula, S. Gounand, A. Beccanini and E. Suder DEN/DANS/DMS/STMF Commissaria à l Energie

More information

V AK (t) I T (t) I TRM. V AK( full area) (t) t t 1 Axial turn-on. Switching losses for Phase Control and Bi- Directionally Controlled Thyristors

V AK (t) I T (t) I TRM. V AK( full area) (t) t t 1 Axial turn-on. Switching losses for Phase Control and Bi- Directionally Controlled Thyristors Applicaion Noe Swiching losses for Phase Conrol and Bi- Direcionally Conrolled Thyrisors V AK () I T () Causing W on I TRM V AK( full area) () 1 Axial urn-on Plasma spread 2 Swiching losses for Phase Conrol

More information

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

More information

SPH3U: Projectiles. Recorder: Manager: Speaker:

SPH3U: Projectiles. Recorder: Manager: Speaker: SPH3U: Projeciles Now i s ime o use our new skills o analyze he moion of a golf ball ha was ossed hrough he air. Le s find ou wha is special abou he moion of a projecile. Recorder: Manager: Speaker: 0

More information

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an

More information

Displacement ( x) x x x

Displacement ( x) x x x Kinemaics Kinemaics is he branch of mechanics ha describes he moion of objecs wihou necessarily discussing wha causes he moion. 1-Dimensional Kinemaics (or 1- Dimensional moion) refers o moion in a sraigh

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

MATHEMATICAL MODELING OF THE TRACTOR-GRADER AGRICULTURAL SYSTEM CINEMATIC DURING LAND IMPROVING WORKS

MATHEMATICAL MODELING OF THE TRACTOR-GRADER AGRICULTURAL SYSTEM CINEMATIC DURING LAND IMPROVING WORKS Bullein of he Transilvania Universiy of Braşov Series II: Foresry Wood Indusry Agriculural Food Engineering Vol. 5 (54) No. 1-2012 MATHEMATICA MODEING OF THE TRACTOR-GRADER AGRICUTURA SYSTEM CINEMATIC

More information

6.003 Homework 1. Problems. Due at the beginning of recitation on Wednesday, February 10, 2010.

6.003 Homework 1. Problems. Due at the beginning of recitation on Wednesday, February 10, 2010. 6.003 Homework Due a he beginning of reciaion on Wednesday, February 0, 200. Problems. Independen and Dependen Variables Assume ha he heigh of a waer wave is given by g(x v) where x is disance, v is velociy,

More information

Heat Transfer. Revision Examples

Heat Transfer. Revision Examples Hea Transfer Revision Examples Hea ransfer: energy ranspor because of a emperaure difference. Thermal energy is ransferred from one region o anoher. Hea ranspor is he same phenomena lie mass ransfer, momenum

More information

USP. Surplus-Production Models

USP. Surplus-Production Models USP Surplus-Producion Models 2 Overview Purpose of slides: Inroducion o he producion model Overview of differen mehods of fiing Go over some criique of he mehod Source: Haddon 2001, Chaper 10 Hilborn and

More information

Sterilization D Values

Sterilization D Values Seriliaion D Values Seriliaion by seam consis of he simple observaion ha baceria die over ime during exposure o hea. They do no all live for a finie period of hea exposure and hen suddenly die a once,

More information

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070

More information

7 The Itô/Stratonovich dilemma

7 The Itô/Stratonovich dilemma 7 The Iô/Sraonovich dilemma The dilemma: wha does he idealizaion of dela-funcion-correlaed noise mean? ẋ = f(x) + g(x)η() η()η( ) = κδ( ). (1) Previously, we argued by a limiing procedure: aking noise

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

5 Dilute systems. 5.1 Weight, drag and Particle Reynolds number

5 Dilute systems. 5.1 Weight, drag and Particle Reynolds number 5 Dilue sysems This chaper considers he behaviour of a single paricle suspended in a fluid. In pracice, he equaions and principles described are used o undersand how a number of paricles behave, provided

More information

ψ(t) = V x (0)V x (t)

ψ(t) = V x (0)V x (t) .93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in

More information

Testing for a Single Factor Model in the Multivariate State Space Framework

Testing for a Single Factor Model in the Multivariate State Space Framework esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics

More information

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran Curren Trends in Technology and Science ISSN : 79-055 8hSASTech 04 Symposium on Advances in Science & Technology-Commission-IV Mashhad, Iran A New for Sofware Reliabiliy Evaluaion Based on NHPP wih Imperfec

More information

APPM 2360 Homework Solutions, Due June 10

APPM 2360 Homework Solutions, Due June 10 2.2.2: Find general soluions for he equaion APPM 2360 Homework Soluions, Due June 10 Soluion: Finding he inegraing facor, dy + 2y = 3e µ) = e 2) = e 2 Muliplying he differenial equaion by he inegraing

More information

2) Of the following questions, which ones are thermodynamic, rather than kinetic concepts?

2) Of the following questions, which ones are thermodynamic, rather than kinetic concepts? AP Chemisry Tes (Chaper 12) Muliple Choice (40%) 1) Which of he following is a kineic quaniy? A) Enhalpy B) Inernal Energy C) Gibb s free energy D) Enropy E) Rae of reacion 2) Of he following quesions,

More information

Keywords: thermal stress; thermal fatigue; inverse analysis; heat conduction; regularization

Keywords: thermal stress; thermal fatigue; inverse analysis; heat conduction; regularization Proceedings Inverse Analysis for Esimaing Temperaure and Residual Sress Disribuions in a Pipe from Ouer Surface Temperaure Measuremen and Is Regularizaion Shiro Kubo * and Shoki Taguwa Deparmen of Mechanical

More information

Unsteady Flow Problems

Unsteady Flow Problems School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.

More information

Suggested Practice Problems (set #2) for the Physics Placement Test

Suggested Practice Problems (set #2) for the Physics Placement Test Deparmen of Physics College of Ars and Sciences American Universiy of Sharjah (AUS) Fall 014 Suggesed Pracice Problems (se #) for he Physics Placemen Tes This documen conains 5 suggesed problems ha are

More information

EE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:

EE650R: Reliability Physics of Nanoelectronic Devices Lecture 9: EE65R: Reliabiliy Physics of anoelecronic Devices Lecure 9: Feaures of Time-Dependen BTI Degradaion Dae: Sep. 9, 6 Classnoe Lufe Siddique Review Animesh Daa 9. Background/Review: BTI is observed when he

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8. Kinemaics Vocabulary Kinemaics and One Dimensional Moion 8.1 WD1 Kinema means movemen Mahemaical descripion of moion Posiion Time Inerval Displacemen Velociy; absolue value: speed Acceleraion Averages

More information

Lab #2: Kinematics in 1-Dimension

Lab #2: Kinematics in 1-Dimension Reading Assignmen: Chaper 2, Secions 2-1 hrough 2-8 Lab #2: Kinemaics in 1-Dimension Inroducion: The sudy of moion is broken ino wo main areas of sudy kinemaics and dynamics. Kinemaics is he descripion

More information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

Motion along a Straight Line

Motion along a Straight Line chaper 2 Moion along a Sraigh Line verage speed and average velociy (Secion 2.2) 1. Velociy versus speed Cone in he ebook: fer Eample 2. Insananeous velociy and insananeous acceleraion (Secions 2.3, 2.4)

More information

Inventory Control of Perishable Items in a Two-Echelon Supply Chain

Inventory Control of Perishable Items in a Two-Echelon Supply Chain Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan

More information

Failure of the work-hamiltonian connection for free energy calculations. Abstract

Failure of the work-hamiltonian connection for free energy calculations. Abstract Failure of he work-hamilonian connecion for free energy calculaions Jose M. G. Vilar 1 and J. Miguel Rubi 1 Compuaional Biology Program, Memorial Sloan-Keering Cancer Cener, 175 York Avenue, New York,

More information

Viscoelastic Catenary

Viscoelastic Catenary Viscoelasic Caenary Anshuman Roy 1 Inroducion This paper seeks o deermine he shape of a hin viscoelasic fluid filamen as i sags under is own weigh. The problem is an exension of he viscous caenary [1]

More information

NUMERICAL SIMULATION OF A LIQUID SODIUM TURBULENT FLOW OVER A BACKWARD FACING STEP WITH A FOUR PARAMETER LOGARITHMIC TURBULENCE MODEL

NUMERICAL SIMULATION OF A LIQUID SODIUM TURBULENT FLOW OVER A BACKWARD FACING STEP WITH A FOUR PARAMETER LOGARITHMIC TURBULENCE MODEL h European Conference on Compuaional Mechanics (ECCM ) 7h European Conference on Compuaional Fluid Dynamics (ECFD 7) 11 June, Glasgow, UK NUMERICAL SIMULATION OF A LIQUID SODIUM TURBULENT FLOW OVER A BACKWARD

More information

Summary of shear rate kinematics (part 1)

Summary of shear rate kinematics (part 1) InroToMaFuncions.pdf 4 CM465 To proceed o beer-designed consiuive equaions, we need o know more abou maerial behavior, i.e. we need more maerial funcions o predic, and we need measuremens of hese maerial

More information

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The

More information

Math 10B: Mock Mid II. April 13, 2016

Math 10B: Mock Mid II. April 13, 2016 Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.

More information

t 2 B F x,t n dsdt t u x,t dxdt

t 2 B F x,t n dsdt t u x,t dxdt Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.

More information

k 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series

k 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series Final Review A Puzzle... Consider wo massless springs wih spring consans k 1 and k and he same equilibrium lengh. 1. If hese springs ac on a mass m in parallel, hey would be equivalen o a single spring

More information

Physics 101 Fall 2006: Exam #1- PROBLEM #1

Physics 101 Fall 2006: Exam #1- PROBLEM #1 Physics 101 Fall 2006: Exam #1- PROBLEM #1 1. Problem 1. (+20 ps) (a) (+10 ps) i. +5 ps graph for x of he rain vs. ime. The graph needs o be parabolic and concave upward. ii. +3 ps graph for x of he person

More information

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H. ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models

More information

Interpretation of special relativity as applied to earth-centered locally inertial

Interpretation of special relativity as applied to earth-centered locally inertial Inerpreaion of special relaiviy as applied o earh-cenered locally inerial coordinae sysems in lobal osiioning Sysem saellie experimens Masanori Sao Honda Elecronics Co., Ld., Oyamazuka, Oiwa-cho, Toyohashi,

More information

MOMENTUM CONSERVATION LAW

MOMENTUM CONSERVATION LAW 1 AAST/AEDT AP PHYSICS B: Impulse and Momenum Le us run an experimen: The ball is moving wih a velociy of V o and a force of F is applied on i for he ime inerval of. As he resul he ball s velociy changes

More information

Second Law. first draft 9/23/04, second Sept Oct 2005 minor changes 2006, used spell check, expanded example

Second Law. first draft 9/23/04, second Sept Oct 2005 minor changes 2006, used spell check, expanded example Second Law firs draf 9/3/4, second Sep Oc 5 minor changes 6, used spell check, expanded example Kelvin-Planck: I is impossible o consruc a device ha will operae in a cycle and produce no effec oher han

More information

NUMERICAL INVESTIGATION OF STROUHAL FREQUENCIES OF TWO STAGGERED BLUFF BODIES

NUMERICAL INVESTIGATION OF STROUHAL FREQUENCIES OF TWO STAGGERED BLUFF BODIES NUMERICAL INVESTIGATION OF STROUHAL FREQUENCIES OF TWO STAGGERED BLUFF BODIES Eswaran M 1, P. Goyal, Anu Dua, G.R. Reddy, R. K. Singh and K.K. Vaze Bhabha Aomic Research Cenre, Mumbai, India. 1 Corresponding

More information

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18 A Firs ourse on Kineics and Reacion Engineering lass 19 on Uni 18 Par I - hemical Reacions Par II - hemical Reacion Kineics Where We re Going Par III - hemical Reacion Engineering A. Ideal Reacors B. Perfecly

More information

Polymerization Technology Laboratory

Polymerization Technology Laboratory Versuch eacion Calorimery Polymerizaion Technology Laboraory eacion Calorimery 1. Subjec Isohermal and adiabaic emulsion polymerizaion of mehyl mehacrylae in a bach reacor. 2. Theory 2.1 Isohermal and

More information

AP Chemistry--Chapter 12: Chemical Kinetics

AP Chemistry--Chapter 12: Chemical Kinetics AP Chemisry--Chaper 12: Chemical Kineics I. Reacion Raes A. The area of chemisry ha deals wih reacion raes, or how fas a reacion occurs, is called chemical kineics. B. The rae of reacion depends on he

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information