Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis

Save this PDF as:
Size: px
Start display at page:

Download "Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis"

Transcription

1 Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis DAVID KATOSHEVSKI Department of Biotechnology and Environmental Engineering Ben-Gurion niversity of the Negev Beer-Sheva ISRAEL Astract: Basic criteria are estalished for the phenomenon of travel grouping, oth for the possiility of it to occur as well as for the associated characteristic time. These criteria are suggested to e incorporated in the rate of spray evaporation, thus enaling for the estimation of the effect of grouping due to periodic oscillation on the overall vapor production of a multi-size spray. It is shown that the relatively small droplets tend to form groups and thus decrease their evaporation rate. On the other hand, the lager droplets tend not to group, thus maximizing their rate of evaporation. Hence, such travel grouping rings closer the evaporation rates of different droplet size ranges. With respect to the opposite effect, that is, the influence of evaporation on grouping, it is shown that once a grouping or a non-grouping mode of ehavior has een estalished, evaporation does not change this ehavior as time passes or, in other words, droplet s history determines a single mode of time-ehavior. Key-Words: Spray Modeling, Sectional Method, Grouping 1 Introduction The configuration in which the motion of evaporating droplets is influenced y an oscillating flow occurs in many spray applications [1]. One such an application involves fuel-droplet vaporization in an oscillating or pulsating comustion chamer. The spatial spreading of the fuel vapors is a key parameter regarding the characteristics of the comustion, while that spreading of the vapors is primarily influenced y droplet trajectories and y the local concentration of the liquid phase [2]. Thus, if we deal with relatively fast moving multi-size spray droplets we may primarily picture the vapors of droplets to form a narrow cylindrical trail, as a function of their size, marking the droplet trajectory, which later diffuse according to the local vapor concentration gradients. With this picture in mind, for the purpose of the current study we consider a Lagrangian type of approach (y Tamour [3]) to follow the dynamics of multi-size evaporating spray droplets, with special emphasize on the phenomenon we denote here as travel grouping. By this notation we refer to cases in which droplets tend to gather into groups (see for example the work y Bellan and Harstad [4]) as a function of various operating conditions to e later addressed. Figure 1 shows a schematic description of the general configuration. Fig. 1. Schematic description of travel-grouping in oscillating flow field. The effect of grouping on the evaporation has een a suject of various studies (see a partial list in [1]). However, non of them deals with the opposite effect, that is, the effect of evaporation on grouping, an effect to e addressed here. We wish to emphasis that we will not make an attempt here to otain a new formulation for the evaporation of a group of droplets, rather, we will enale criteria to predict when the use of such availale rates should e turned on. Suggestion of modification for the evaporation rate itself will also e addressed, ut only in a qualitative manner.

2 In the route of the current work, the following section deals with the equation of motion of the droplets associated with the effect of the oscillations of the host-gas flow. This asic issue was discussed in a recent work y Katoshevski et al., [5] and thus we will show only the major points with that respect, while a new criteria for grouping will e presented here. Then, we will address the dynamics of the droplets coupled with their evaporation. 2 Droplet motion and grouping/nongrouping criteria For a qualitative investigation of the prolem we consider droplet motion in a one-dimensional periodic fluid flow. As a model of such flow one can apply the following mathematical form [6]: ut, x a sinkx t (1) where u t, x is a fluid velocity at a time t at a location x, a is a the mean flow velocity, is the amplitude of the velocity oscillation, k is the wave numer, and is the angular velocity which is linearly proportional to the oscillation frequency. Without less in generality one can assume that a, 0. The flow has a wave length of L 2 / k along the x-axis and a period T 2 / in time. The wave that is formed propagates along the x-axis in a velocity of w / k. Let us normalize the velocities y w, the distance y k -1 1 and the time y. This leads to the following equation of motion for the droplets (note that from now on we refer to normalized velocities only, and for convenience, omit any indication for designating a normalized character), dud 1 u ud, (2) dt St where St is the Stokes numer which is linear with the square of the droplet diameter. u d denotes the velocity of a droplet. After sustituting the explicit form of the host gas velocity u t, x from Eq. 1 and introducing a new variale x t, the droplet equation of motion ecomes: 2 d 1 d a 1 sin. (3) 2 dt St dt St St Now, replacing time t y ct, where c / St, leads to sin, (4) where the prime denotes a derivative with respect to. In the aove, the parameters and are defined as, 1/ St (5) and a 1 /. (6) In the previous study [4] the analysis of the aove equation has revealed that there are two major modes of droplet grouping ehavior in such a flow configuration that can e descried in the plane. The first mode of ehavior is pronounced y motion of droplets in a fixed travel-group, that is, moving constantly within the original family while the second major mode is the shifting of droplets from one family to the other in addition to reakup of groups. As and are associated with the operating conditions, such as the host gas velocity, the frequency of the oscillations as well as droplet size, we may say that in certain operating conditions we expect to oserve mode I, that is the family type of group motion and in other conditions mode II is expected, that is the non-family type, will e oserved. In addition to these modes of grouping, for certain operating conditions clear non-grouping situations are predicted. It is important to note that Eq. 2 is valid for low droplet Reynolds numers, ased on the velocitydifference etween the host gas velocity and that of the droplet. But, from our calculations we conclude that using a drag coefficient C D as a function of the Reynolds numer for Re>1 leads to grouping tendency which is even greater. In most spray cases the desire is to reakup groups of droplets in order to maximize evaporation. From further analysis of the plane we ring here new criteria for reaching non-grouping, or nearly nongrouping situations. The following relations should e maintained y the operating conditions, for low particle Reynolds numers 1 (7a), OR, 1 AND / 2 (7) and to ensure non-grouping, or nearly non-grouping situations even for a wider range of droplet Reynolds numers we have to maintain the following relation arising from Eqs. 5-7, a 1 NG St 1 (8) where NG stands for Non-Grouping. We note here, that according to our numerical investigation we conclude that grouping is not distinguished already when NG >2 (corresponding to the division y 2 in Eq. 7), and the aove relation is for the sake of generality.

3 Note that the aove relation is presented here for the first time and was not otained in the previous study y Katoshevski et al., [5]. Keeping in mind that the velocities are normalized y w / k and that the Stokes numer is linear with the frequency, then, if one increases the frequency NG will also increase and the tendency will e toward non-grouping situations. The same ehavior of tendency to group is true for an increase in droplet size, or an increase in droplet liquid density, oth causing an increase in St. The aove criteria can e used for correcting the evaporation rates according to grouping/non-grouping situations. This will e addressed in the next section, while here we first would like to demonstrate some related features of the grouping ehavior as shown in Figs Figure 2 presents a clear grouping case, in which the droplets of 30 microns in diameter that are initially dispersed along a distance of one wave-length (L in length) along the x axis form two travel groups. For that case the value of NG is Then, if we increase the frequency, so that (which is 2 times the frequency) increases from 1000 Hz to 1800 Hz, NG ecomes 1.04, and one of the two groups splits up into two, that is, a secondary reakup occurs of the overall group, as shown in Fig. 3. Our calculations show that, a further increase of to 2000 Hz leads to an increase of NG to 1.11 and causes another stage of group reakup. Now, we jump a few steps ahead and increase the droplet size to 70 microns and to 6000 Hz in order to show a clear non-grouping case, as presented in Fig. 4. In that case NG is well over unity (equals 4.45) as in the criterion of Eq.8. In order to account for the grouping effect on evaporation, which is the topic of the next section, it will e useful to achieve an approximate analytical solution of the equation of droplet motion under the grouping criteria. For such an approximation, ased on Eqs. 7a-, we may consider the grouping criterion as, and for an asymptotic solution account for and define Equation 4 then ecomes sin (9) and the solution can e represented asymptotically as (0) (1)... (10) (0) The first approximation will e const., which implies from the definition of x-t) that the group velocity is aout the same as that of the wave propagation. Thus, we may conclude that for the sake of descriing the clear-grouping situation we can account for the dimensional velocity of the groups as equal to G / k, where the index G denotes group. This value of G is also deduced from the trajectories presented in Fig. 2. In order to satisfy the initial condition of the particle initial position we set (0) (0). For the correction term which is in the order of we allow ourselves to linearize the sinus term and otain the following approximation, ( 0) t sin t (11) St St From the aove, one can estimate the dimensional time which passes until grouping is clearly estalished, which is, St tg tˆ (12) where t G is the dimensional characteristic time for grouping to occur, where tˆ is a constant that equals 2. For example, with respect to the parameters of Fig. 2, this corresponds to the dimensional time of sec, which agrees well with the trajectories shown in that figure. Fig. 2. Trajectories of droplets in the grouping mode, for non-evaporating droplets of 30 microns in diameter, and 1000 Hz.

4 Fig. 3. Non-grouping mode, for non-evaporating droplets of 30 microns in diameter, and 1800 Hz. Fig. 4. Non-grouping mode, for non-evaporating droplets of 70 microns in diameter, 6000 Hz. 3 Travel grouping and evaporation coupling We now turn to suggest how to incorporate the effect of grouping on evaporation, and visa- versa, that is, the effect of evaporation on grouping. For the frame of our qualitative description we choose to employ the Lagrangian method [3, 7], for a mono-disperse spray, dq C Q (13) dt where Q is the total mass of droplets in a control volume. C is an evaporation coefficient. In case a multi-sized spray is addressed, the overall size range is divided into size sections, and the aove equation will then incorporate an additional term accounting for droplets that move from one size section to the other as their size is reduced due to evaporation [3]. For the purpose of the current study we treat a mono-size spray, which can e used also for the case where similar-sized groups of droplets of a multi-size spray move closely together. The evaporation coefficient is general a function of the droplet size, the local temperature, the liquid specifications, and can also incorporate the effect of the velocity difference etween that of the droplet and that of the host gas y a Ranz-Marshal type of correlation [8]. As mentioned aove, in the frame of our approximations we use a Lagrangian type of an approach [3], and for the estimation of overall production of the vapors we may write, dm C Q (14) dt When dealing with the production of vapors y specific range of droplet sizes that may travel within a group we will sum up the contriutions of the various groups to the overall vapor production. In general, we can estimate the overall vapor production y dividing the time frame into two periods which are (I) up to t G and (II) after t G. In the second period, the fact that grouping is accounted for will hinder evaporation comparing to the evaporation rate from isolated droplets. This effect on evaporation will e accounted for y introducing a theoretical correction to the evaporation coefficient. The correction for the evaporation coefficient due to the grouping is pronounced y multiplying these coefficients y a function E, where E( G) 1, if NG 1 and t tg E (15) 1 otherwise and G is the known correction term suggested y Chiu and co-workers (Sirignano, 1999, and references therein) for evaporation and comustion of non-isolate droplets. Based on studies on group comustion (see in Ref. [1]) G is primarily a function of the Lewis numer, Le, the total numer of droplets in the group N, the averaged droplet radius r and the spacing etween the centers of the droplets s, 2 / 3 G ~ Le N ( r / s) (16) Thus, in addition of eing a function of the sectional size range, the liquid properties, and the temperature [3], the sectional evaporation coefficients may also e

5 adjusted y the aove E function. Adding the Ranz- Marshal correlation [8] denoted here as R M, we may write C ~ RM (Re, Pr) E( G, NG, tg ) (17) where Re and Pr are the Reynolds and Prandtl numers, respectively. Hence, in a formal way of representation, the overall vapor production from a multi-sized spray undergoing grouping and is eing followed numerically, may e estimated y tg t m( t) C Q dt E C Q dt (18) o t G where the evaporation coefficient C is for nongrouping situation ut includes the effect of the Ranz- Marshal correlation [8]. When a spray flame is to e considered (as for example, y Greenerg [7]), which is not the case in the present work, it will e important to add the effect of diffusion in order to descrie the time evolution of the spatial vapor distriution in the comustion chamer, along with the comustion instaility which is a result of flow oscillations. With respect to droplet coalescence, the same treatment can e applied, that is, as grouping will enhance the sectional rates of coalescence, these rates [3] should e adjusted during the time greater than t G. 4 Demonstration and conclusions For demonstration, let us consider the case of two mono-sized sprays of 10 m and 70 m in diameter. We have calculated the trajectories of these droplets, without evaporation, and then added the effect of droplet size reduction, employing the commonly used D-square law of evaporation (and we modify Eq. (1) so that the time derivative will account for droplet size reduction). For a qualitative demonstration we have the square of the droplet diameter reduced dramatically during the time of 0.08 seconds, to e consistent with the time scale of the former figures. Thus, in Fig. 5a the droplet size reduces from 10 microns to 1 micron. Comparing to our nonevaporating calculations we find that the reduction in size does not affect the grouping ehavior of these small droplets. This coincides with the aove mentioned conclusion, that smaller droplets have higher tendency to form groups. Thus, after grouping has occurred, a further reduction in size will not reak that group. With respect to the larger droplets (Fig. 5), their diameter reduces from 70 microns to 6 microns during the same period of time of 0.08 seconds. Although the size of these relatively large droplets reduces consideraly with time as they evaporate, our calculations show that the tendency of these droplets to group is negligile, and this is due to their history as non-grouping droplets. Thus, for large droplets, which initially do not form a travel-group, evaporation and coincidently reduction of size, will not lead to grouping ehavior. Now, we add to the aove the effect of the slow-down in evaporation as grouping occurs. This is done y employing the function E(G) of Eq. 15. Grouping is relevant to the small droplets as already mentioned. Thus, we do not apply the function E(G) on the large droplets of Fig. 5. As already concluded, from Fig. 5, the rate of evaporation does not change the grouping ehavior of the small droplets once grouping has een estalished. Hence, slow-down in evaporation rate will not change the trajectories of Fig. 5a. What is then left is to estimate the change in the vapor production, or the change in evaporation rate of the group of small droplets compared with the same numer of isolated ones. Laowsky [9] has shown that the D-square law essentially descries the evaporation rate of nonisolated droplets, ut with a decrease factor, represented here y the function E(G). That study, as well as others [1], have shown that depending on spacing etween droplets, and operating conditions, the reduction in the evaporation rate can e sustantially larger than 50%. As our simple mathematical modeling can represent only a qualitative picture, for that frame of work we can use such an estimation for the purpose of a qualitative demonstration. In addition, from the study y Tamour [3] regarding the aove mentioned sectional approach we may estimate that the averaged sectional evaporation coefficient C d (where the su-index d denotes droplet-size) of droplet initially with size of 10 microns is 3 times larger than the corresponding rate of a 70 micron droplet, when they are isolated. Hence grouping can ring the evaporation rates of these two classes of sizes to very close values. The amounts of vapors produced from two sizes-groups together as a function of time m(t) can e evaluated y Eqs. (14) and (18). As a first-order approximation, accounting for t>t G we otain, m( t) Q Q (0) 1 ~ (0)! 1 Exp( C70 t) " ~ ~ Exp( C t )! 2 Exp( EC ( t )" 10 G 10 t G (19)

6 where the tilde represents an averaged value of vaporization rate. The caliration of the coefficients in Eq. (19) can e done when experiments are involved. In conclusion, the smaller size droplets tend to form groups and thus decrease their evaporation rate, while the lager droplets do not tend to group, thus maximizing their rate of evaporation. This may ring closer the evaporation rates of different size ranges. The criterion for travel grouping introduced here, along with the characteristic time until grouping may occur can e of use in cases where multi-size spray droplets are sujected to an oscillatory flow field. The current work also points out that once grouping or non-grouping mode of ehavior has een estalished, nor evaporation nor time will cause a switch of that mode. a Acknowledgment This research was supported y the fund for promotion of research at the Ben-Gurion niversity of the Negev. References: [1] W. A. Sirignano, Fluid Dynamics and Transport of Droplets and Sprays, Camridge niversity Press, [2] N. A. Chigier, Progress in. Energy and Comustion Science, Vol. 2, 1976, pp [3] Y. Tamour, A Lagrangian Sectional Approach for Simulating Droplet Size Distriution of Vaporizing Fuel Sprays in a Turulent Jet, Comustion and Flame, Vol. 60, 1985, pp [4] J. Bellan, and K. Harstad, The Dynamics of Dense and Dilute Clusters of Drops Evaporating in Large Coherent Vortices, Proceedings of the Comustion Institute, Vol. 23, 1990, pp [5] D. Katoshevski, Z. Dodin, and G. Ziskind, Aerosol Particle Clustering in Oscillating Flows- Mathematical Analysis, Atomization and Sprays, in press, [6] S. C. Crow, F. H. Champagne, Orderly Structure in Jet Turulence, Journal of Fluid Mechanics, Vol. 48, 1971, pp [7] J. B. Greenerg, The Burke-Schumann Diffusion Flame Revisited- with Fuel Spray Injection, Comustion and Flame, Vol. 77, 1989, pp [8] V. E. Ranz, and W. R. Marshall, Evaporation from Drops, Chemical Engineering Progress, Vol. 48, 1952, pp [9] M. Laowsky, Calculation of Burning Rates of Interacting Fuel Droplets, Comustion Science and Technology, Vol. 22, 1980, pp Fig. 5. Trajectories of evaporating droplets: (a) 10 microns, and () 70 microns Hz.

Induced Particle Grouping and Coagulation in Diesel Exhaust

Induced Particle Grouping and Coagulation in Diesel Exhaust Induced Particle Grouping and Coagulation in Diesel Exhaust Itamar Hite, Michal Ruzal Mendelevich,, Tal Shakked, Prof. Eran Sher and Prof. David Katoshevski Ben-Gurion niversity of the Negev Beer-Sheva

More information

Modeling Droplet Breakup Processes in Bio-fuel Diesel Engines under Micro-explosion Conditions

Modeling Droplet Breakup Processes in Bio-fuel Diesel Engines under Micro-explosion Conditions ILASS Americas, th Annual Conference on Liquid Atomization and Spray Systems, Chicago, IL, May 7 Modeling Droplet Breakup Processes in Bio-fuel Diesel Engines under Micro-explosion Conditions K. T. Wang

More information

On Temporal Instability of Electrically Forced Axisymmetric Jets with Variable Applied Field and Nonzero Basic State Velocity

On Temporal Instability of Electrically Forced Axisymmetric Jets with Variable Applied Field and Nonzero Basic State Velocity Availale at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 19-966 Vol. 6, Issue 1 (June 11) pp. 7 (Previously, Vol. 6, Issue 11, pp. 1767 178) Applications and Applied Mathematics: An International Journal

More information

Section 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)).

Section 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)). Difference Equations to Differential Equations Section 8.5 Applications: Pendulums Mass-Spring Systems In this section we will investigate two applications of our work in Section 8.4. First, we will consider

More information

1D spirals: is multi stability essential?

1D spirals: is multi stability essential? 1D spirals: is multi staility essential? A. Bhattacharyay Dipartimento di Fisika G. Galilei Universitá di Padova Via Marzolo 8, 35131 Padova Italy arxiv:nlin/0502024v2 [nlin.ps] 23 Sep 2005 Feruary 8,

More information

Available online at Energy Procedia 100 (2009) (2008) GHGT-9

Available online at   Energy Procedia 100 (2009) (2008) GHGT-9 Availale online at www.sciencedirect.com Energy Procedia (29) (28) 655 66 Energy Procedia www.elsevier.com/locate/procedia www.elsevier.com/locate/xxx GHGT-9 Pre-comustion CO 2 capture for IGCC plants

More information

Lecture 9 Laminar Diffusion Flame Configurations

Lecture 9 Laminar Diffusion Flame Configurations Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.

More information

Effect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium

Effect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue Ver. III (Jan. - Fe. 06), 5-65 www.iosrjournals.org Effect of Uniform Horizontal Magnetic Field on Thermal Instaility

More information

Topics in Other Lectures Droplet Groups and Array Instability of Injected Liquid Liquid Fuel-Films

Topics in Other Lectures Droplet Groups and Array Instability of Injected Liquid Liquid Fuel-Films Lecture Topics Transient Droplet Vaporization Convective Vaporization Liquid Circulation Transcritical Thermodynamics Droplet Drag and Motion Spray Computations Turbulence Effects Topics in Other Lectures

More information

Luis Manuel Santana Gallego 100 Investigation and simulation of the clock skew in modern integrated circuits. Clock Skew Model

Luis Manuel Santana Gallego 100 Investigation and simulation of the clock skew in modern integrated circuits. Clock Skew Model Luis Manuel Santana Gallego 100 Appendix 3 Clock Skew Model Xiaohong Jiang and Susumu Horiguchi [JIA-01] 1. Introduction The evolution of VLSI chips toward larger die sizes and faster clock speeds makes

More information

Robot Position from Wheel Odometry

Robot Position from Wheel Odometry Root Position from Wheel Odometry Christopher Marshall 26 Fe 2008 Astract This document develops equations of motion for root position as a function of the distance traveled y each wheel as a function

More information

A theory for one dimensional asynchronous waves observed in nonlinear dynamical systems

A theory for one dimensional asynchronous waves observed in nonlinear dynamical systems A theory for one dimensional asynchronous waves oserved in nonlinear dynamical systems arxiv:nlin/0510071v1 [nlin.ps] 28 Oct 2005 A. Bhattacharyay Dipartimento di Fisika G. Galilei Universitá di Padova

More information

CAV2001:SESSIONA

CAV2001:SESSIONA CAV001:SESSIONA1.004 1 Effects of model size and free stream nuclei on tip vortex cavitation inception scaling Young Shen*, Georges Chahine**, Chao-Tsung Hsiao**, and Stuart Jessup* *NSWCCD, West Bethesda,

More information

Math 216 Second Midterm 28 March, 2013

Math 216 Second Midterm 28 March, 2013 Math 26 Second Midterm 28 March, 23 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that

More information

Nonsteady interior ballistics of cylindrical-grain solid rocket motors

Nonsteady interior ballistics of cylindrical-grain solid rocket motors Computational Ballistics II 281 Nonsteady interior allistics of cylindrical-grain solid rocket motors D. R. Greatrix Department of Aerospace Engineering, Ryerson University, Canada Astract A numerical

More information

Buckling Behavior of Long Symmetrically Laminated Plates Subjected to Shear and Linearly Varying Axial Edge Loads

Buckling Behavior of Long Symmetrically Laminated Plates Subjected to Shear and Linearly Varying Axial Edge Loads NASA Technical Paper 3659 Buckling Behavior of Long Symmetrically Laminated Plates Sujected to Shear and Linearly Varying Axial Edge Loads Michael P. Nemeth Langley Research Center Hampton, Virginia National

More information

EFFECTS OF STRONG TEMPERATURE GRADIENT ON A COMPRESSIBLE TURBULENT CHANNEL FLOW

EFFECTS OF STRONG TEMPERATURE GRADIENT ON A COMPRESSIBLE TURBULENT CHANNEL FLOW th International Symposium on Turulence and Shear Flo Phenomena (TSFP, Chicago, USA, July, 7 EFFECTS OF STRONG TEMPERATURE GRADIENT ON A COMPRESSIBLE TURBULENT CHANNEL FLOW Mitsuhiro Nagata Mechanical

More information

Critical value of the total debt in view of the debts. durations

Critical value of the total debt in view of the debts. durations Critical value of the total det in view of the dets durations I.A. Molotov, N.A. Ryaova N.V.Pushov Institute of Terrestrial Magnetism, the Ionosphere and Radio Wave Propagation, Russian Academy of Sciences,

More information

Lecture 12. Droplet Combustion Spray Modeling. Moshe Matalon

Lecture 12. Droplet Combustion Spray Modeling. Moshe Matalon Lecture 12 Droplet Combustion Spray Modeling Spray combustion: Many practical applications liquid fuel is injected into the combustion chamber resulting in fuel spray. Spray combustion involves many physical

More information

#A50 INTEGERS 14 (2014) ON RATS SEQUENCES IN GENERAL BASES

#A50 INTEGERS 14 (2014) ON RATS SEQUENCES IN GENERAL BASES #A50 INTEGERS 14 (014) ON RATS SEQUENCES IN GENERAL BASES Johann Thiel Dept. of Mathematics, New York City College of Technology, Brooklyn, New York jthiel@citytech.cuny.edu Received: 6/11/13, Revised:

More information

PARAMETER IDENTIFICATION, MODELING, AND SIMULATION OF A CART AND PENDULUM

PARAMETER IDENTIFICATION, MODELING, AND SIMULATION OF A CART AND PENDULUM PARAMETER IDENTIFICATION, MODELING, AND SIMULATION OF A CART AND PENDULUM Erin Bender Mechanical Engineering Erin.N.Bender@Rose-Hulman.edu ABSTRACT In this paper a freely rotating pendulum suspended from

More information

arxiv: v1 [hep-ph] 20 Dec 2012

arxiv: v1 [hep-ph] 20 Dec 2012 Decemer 21, 2012 UMD-PP-012-027 Using Energy Peaks to Count Dark Matter Particles in Decays arxiv:1212.5230v1 [hep-ph] 20 Dec 2012 Kaustuh Agashe a, Roerto Franceschini a, Doojin Kim a, and Kyle Wardlow

More information

Field-induced magnetization of a free-electron gas in thin films

Field-induced magnetization of a free-electron gas in thin films Materials Science-Poland, Vol. 5, No., 007 ield-induced magnetization of a free-electron gas in thin films T. BALCERZAK 1* 1 Department of Solid State Physics, University of Łódź, Pomorska 149/153, 90-36

More information

cen29305_ch08.qxd 11/30/05 3:05 PM Page 451 INTERNAL FORCED CONVECTION CHAPTER 8 Liquid or gas flow through pipes or ducts is commonly used in heating

cen29305_ch08.qxd 11/30/05 3:05 PM Page 451 INTERNAL FORCED CONVECTION CHAPTER 8 Liquid or gas flow through pipes or ducts is commonly used in heating cen29305_ch08.qxd 11/30/05 3:05 PM Page 451 INTERNAL FORCED CONVECTION CHAPTER 8 Liquid or gas flow through pipes or ducts is commonly used in heating and cooling applications. The fluid in such applications

More information

Theoretical Developments in Group Combustion of Droplets and Sprays

Theoretical Developments in Group Combustion of Droplets and Sprays Theoretical Developments in Group Combustion of Droplets and Sprays William A. Sirignano University of California, Irvine Collaborations: Guang Wu, current student; Randall Imaoka, former student, US Navy;

More information

Concentration of magnetic transitions in dilute magnetic materials

Concentration of magnetic transitions in dilute magnetic materials Journal of Physics: Conference Series OPEN ACCESS Concentration of magnetic transitions in dilute magnetic materials To cite this article: V I Beloon et al 04 J. Phys.: Conf. Ser. 490 065 Recent citations

More information

Pair diffusion in quasi-one- and quasi-two-dimensional binary colloid suspensions

Pair diffusion in quasi-one- and quasi-two-dimensional binary colloid suspensions THE JOURNAL OF CHEMICAL PHYSICS 126, 134908 2007 Pair diffusion in quasi-one- and quasi-two-dimensional binary colloid suspensions David T. Valley, Stuart A. Rice, Bianxiao Cui, a and Hau My Ho Department

More information

TOPICAL PROBLEMS OF FLUID MECHANICS 17 ONE-DIMENSIONAL TEMPERATURE DISTRIBUTION OF CONDENSING ANNULAR FINS OF DIFFERENT PROFILES

TOPICAL PROBLEMS OF FLUID MECHANICS 17 ONE-DIMENSIONAL TEMPERATURE DISTRIBUTION OF CONDENSING ANNULAR FINS OF DIFFERENT PROFILES TOPICAL PROBLEMS OF FLUID MECHANICS 17 ONE-DIMENSIONAL TEMPERATURE DISTRIBUTION OF CONDENSING ANNULAR FINS OF DIFFERENT PROFILES A. Bouraaa 1, 2, M. Saighi 2, K. Salhi 1, A. Hamidat 1 and M. M. Moundi

More information

SEG/New Orleans 2006 Annual Meeting. Non-orthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University

SEG/New Orleans 2006 Annual Meeting. Non-orthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University Non-orthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University SUMMARY Wavefield extrapolation is implemented in non-orthogonal Riemannian spaces. The key component is the development

More information

Mixing at the External Boundary of a Submerged Turbulent Jet

Mixing at the External Boundary of a Submerged Turbulent Jet Mixing at the External Boundary of a Submerged Turbulent Jet A. Eidelman, T. Elperin, N. Kleeorin, I. Rogachevskii, I. Sapir-Katiraie The Ben-Gurion University of the Negev, Beer-Sheva, Israel G. Hazak

More information

Optimal Area Power Efficiency in Cellular Networks

Optimal Area Power Efficiency in Cellular Networks Optimal Area Power Efficiency in Cellular Networks Bhanukiran Peraathini 1, Marios Kountouris, Mérouane Deah and Alerto Conte 1 1 Alcatel-Lucent Bell Las, 916 Nozay, France Department of Telecommunications,

More information

Motion Analysis of Euler s Disk

Motion Analysis of Euler s Disk Motion Analysis of Euler s Disk Katsuhiko Yaada Osaka University) Euler s Disk is a nae of a scientific toy and its otion is the sae as a spinning coin. In this study, a siple atheatical odel is proposed

More information

Hot Electron Production in Plasmas Illuminated by Intense Lasers

Hot Electron Production in Plasmas Illuminated by Intense Lasers JETP Letters, Vol. 81, No. 1, 005, pp. 1 5. From Pis ma v Zhurnal Éksperimental noœ i Teoreticheskoœ Fiziki, Vol. 81, No. 1, 005, pp. 7. Original English Text Copyright 005 y Balakin, Fraiman, Fisch. Hot

More information

Three-ring tautomerism of the 2 -isoxazoline - 2 -pyrazoline - 1,3,4-thiadiazine system

Three-ring tautomerism of the 2 -isoxazoline - 2 -pyrazoline - 1,3,4-thiadiazine system Three-ring tautomerism of the 2 -isoxazoline - 2 -pyrazoline - 1,3,4-thiadiazine system Andrei Yu. Ershov, Maksim V. Mokeev, Elena V. Beloorodova, and Alexander V. Grianov Institute of Macromolecular Compounds

More information

Dynamical Systems Solutions to Exercises

Dynamical Systems Solutions to Exercises Dynamical Systems Part 5-6 Dr G Bowtell Dynamical Systems Solutions to Exercises. Figure : Phase diagrams for i, ii and iii respectively. Only fixed point is at the origin since the equations are linear

More information

One-dimensional Spray Combustion Optimization with a Sequential Linear Quadratic Algorithm

One-dimensional Spray Combustion Optimization with a Sequential Linear Quadratic Algorithm One-dimensional Spray Combustion Optimization with a Sequential Linear Quadratic Algorithm Justin A. Sirignano, Luis Rodriguez, Athanasios Sideris, and William A. Sirignano Department of Mechanical and

More information

BOUSSINESQ-TYPE MOMENTUM EQUATIONS SOLUTIONS FOR STEADY RAPIDLY VARIED FLOWS. Yebegaeshet T. Zerihun 1 and John D. Fenton 2

BOUSSINESQ-TYPE MOMENTUM EQUATIONS SOLUTIONS FOR STEADY RAPIDLY VARIED FLOWS. Yebegaeshet T. Zerihun 1 and John D. Fenton 2 ADVANCES IN YDRO-SCIENCE AND ENGINEERING, VOLUME VI BOUSSINESQ-TYPE MOMENTUM EQUATIONS SOLUTIONS FOR STEADY RAPIDLY VARIED FLOWS Yeegaeshet T. Zerihun and John D. Fenton ABSTRACT The depth averaged Saint-Venant

More information

Chaos and Dynamical Systems

Chaos and Dynamical Systems Chaos and Dynamical Systems y Megan Richards Astract: In this paper, we will discuss the notion of chaos. We will start y introducing certain mathematical concepts needed in the understanding of chaos,

More information

Numerical Studies of Droplet Deformation and Break-up

Numerical Studies of Droplet Deformation and Break-up ILASS Americas 14th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2001 Numerical Studies of Droplet Deformation and Break-up B. T. Helenbrook Department of Mechanical and

More information

arxiv:physics/ v1 [physics.plasm-ph] 7 Apr 2006

arxiv:physics/ v1 [physics.plasm-ph] 7 Apr 2006 Europhysics Letters PREPRINT arxiv:physics/6455v [physics.plasm-ph] 7 Apr 26 Olique electromagnetic instailities for an ultra relativistic electron eam passing through a plasma A. Bret ETSI Industriales,

More information

Mathematical Ideas Modelling data, power variation, straightening data with logarithms, residual plots

Mathematical Ideas Modelling data, power variation, straightening data with logarithms, residual plots Kepler s Law Level Upper secondary Mathematical Ideas Modelling data, power variation, straightening data with logarithms, residual plots Description and Rationale Many traditional mathematics prolems

More information

AMS 147 Computational Methods and Applications Lecture 13 Copyright by Hongyun Wang, UCSC

AMS 147 Computational Methods and Applications Lecture 13 Copyright by Hongyun Wang, UCSC Lecture 13 Copyright y Hongyun Wang, UCSC Recap: Fitting to exact data *) Data: ( x j, y j ), j = 1,,, N y j = f x j *) Polynomial fitting Gis phenomenon *) Cuic spline Convergence of cuic spline *) Application

More information

FLOW AND HEAT-TRANSFER MODELLING OF THREE-DIMENSIONAL JET IMPINGEMENT ON A CONCAVE SURFACE

FLOW AND HEAT-TRANSFER MODELLING OF THREE-DIMENSIONAL JET IMPINGEMENT ON A CONCAVE SURFACE FLOW AND HEAT-TRANSFER MODELLING OF THREE-DIMENSIONAL JET IMPINGEMENT ON A CONCAVE SURFACE Tim J. Craft Hector Iacovides Nor A. Mostafa School of Mechanical Aerospace & Civil Engineering The University

More information

Characteristics of flow over the free overfall of triangular channel

Characteristics of flow over the free overfall of triangular channel MATEC We of Conferences 16, 03006 (018) BCEE3-017 https://doi.org/10.101/matecconf/0181603006 Characteristics of flow over the free overfall of triangular channel Raad Irzooki 1,*, Safa Hasan 1 Environmental

More information

P = ρ{ g a } + µ 2 V II. FLUID STATICS

P = ρ{ g a } + µ 2 V II. FLUID STATICS II. FLUID STATICS From a force analysis on a triangular fluid element at rest, the following three concepts are easily developed: For a continuous, hydrostatic, shear free fluid: 1. Pressure is constant

More information

A Stable and Convergent Finite Difference Scheme for 2D Incompressible Nonlinear Viscoelastic Fluid Dynamics Problem

A Stable and Convergent Finite Difference Scheme for 2D Incompressible Nonlinear Viscoelastic Fluid Dynamics Problem Applied and Computational Mathematics 2018; (1): 11-18 http://www.sciencepulishinggroup.com/j/acm doi: 10.11648/j.acm.2018001.12 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Stale and Convergent

More information

Module 9: Further Numbers and Equations. Numbers and Indices. The aim of this lesson is to enable you to: work with rational and irrational numbers

Module 9: Further Numbers and Equations. Numbers and Indices. The aim of this lesson is to enable you to: work with rational and irrational numbers Module 9: Further Numers and Equations Lesson Aims The aim of this lesson is to enale you to: wor with rational and irrational numers wor with surds to rationalise the denominator when calculating interest,

More information

Combustion MATHEMATICAL MODEL FOR TRANSIENT. S. M. Frolov Λ,F.S.Frolov Λ, and B. Basara y

Combustion MATHEMATICAL MODEL FOR TRANSIENT. S. M. Frolov Λ,F.S.Frolov Λ, and B. Basara y Combustion MATHEMATICAL MODEL FOR TRANSIENT DROPLET VAPORIZATION S. M. Frolov Λ,F.S.Frolov Λ, and B. Basara y Λ N. N. Semenov Institute of Chemical Physics Russian Academy of Sciences Moscow, Russia y

More information

TP A.18 Distance required for stun and natural roll to develop for different tip offsets

TP A.18 Distance required for stun and natural roll to develop for different tip offsets technical proof technical proof TP A.18 Distance required for stun and natural roll to develop for different tip offsets supporting: The Illustrated Principles of Pool and Billiards http://illiards.colostate.edu

More information

CULAR CATTERING. During the first quarter of this century it became. h~,, ~ " " " " " L. Monchick

CULAR CATTERING. During the first quarter of this century it became. h~,, ~      L. Monchick The characteristic interference phenomena of wave mechanics usually can e inferred from experiment only in an indirect manner. Molecular and atomic scattering exhiit this effect directly, ecause thf( triution

More information

Estimation of Hottest Spot Temperature in Power Transformer Windings with NDOF and DOF Cooling

Estimation of Hottest Spot Temperature in Power Transformer Windings with NDOF and DOF Cooling Transactions D: Computer Science & Engineering and Electrical Engineering Vol. 16, No. 2, pp. 163{170 c Sharif University of Technology, Decemer 2009 Research Note Estimation of Hottest Spot Temperature

More information

ragsdale (zdr82) HW7 ditmire (58335) 1 The magnetic force is

ragsdale (zdr82) HW7 ditmire (58335) 1 The magnetic force is ragsdale (zdr8) HW7 ditmire (585) This print-out should have 8 questions. Multiple-choice questions ma continue on the net column or page find all choices efore answering. 00 0.0 points A wire carring

More information

OPTIMAL WAVELENGTH SELECTION ALGORITHM OF NON-SPHERICAL PARTICLE SIZE DISTRIBUTION BASED ON THE LIGHT EXTINCTION DATA

OPTIMAL WAVELENGTH SELECTION ALGORITHM OF NON-SPHERICAL PARTICLE SIZE DISTRIBUTION BASED ON THE LIGHT EXTINCTION DATA THERMAL SCIENCE, Year 202, Vol. 6, No. 5, pp. 353-357 353 OPTIMAL WAVELENGTH SELECTION ALGORITHM OF NON-SPHERICAL PARTICLE SIZE ISTRIBUTION BASE ON THE LIGHT EXTINCTION ATA y Hong TANG * College of Metrology

More information

Superluminal Hidden Communication as the Underlying Mechanism for Quantum Correlations: Constraining Models

Superluminal Hidden Communication as the Underlying Mechanism for Quantum Correlations: Constraining Models 38 Brazilian Journal of Physics, vol. 35, no. A, June, 005 Superluminal Hidden Communication as the Underlying Mechanism for Quantum Correlations: Constraining Models Valerio Scarani and Nicolas Gisin

More information

Notes to accompany Continuatio argumenti de mensura sortis ad fortuitam successionem rerum naturaliter contingentium applicata

Notes to accompany Continuatio argumenti de mensura sortis ad fortuitam successionem rerum naturaliter contingentium applicata otes to accompany Continuatio argumenti de mensura sortis ad fortuitam successionem rerum naturaliter contingentium applicata Richard J. Pulskamp Department of Mathematics and Computer Science Xavier University,

More information

Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores

Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores PHYSICAL REVIEW E VOLUME 62, NUMBER 1 JULY 2000 Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores A. Verdaguer and J. A. Padró Departament de Física

More information

STUDY ON BUBBLE BEHAVIOR OF INERT GASES AT ENTRANCE NOZZLE IN SODIUM-COOLED FAST REACTOR

STUDY ON BUBBLE BEHAVIOR OF INERT GASES AT ENTRANCE NOZZLE IN SODIUM-COOLED FAST REACTOR NTHAS8: The Eighth Japan-Korea Symposium on Nuclear Thermal Hydraulics and Safety Beppu, Japan, Decemer 9-12, 2012 Paper Numer N8P1081 STUDY ON BUBBLE BEHAVIOR OF INERT GASES AT ENTRANCE NOZZLE IN SODIUM-COOLED

More information

A COMPUTATIONAL APPROACH TO DESIGN CODES FOR TSUNAMI- RESISTING COASTAL STRUCTURES

A COMPUTATIONAL APPROACH TO DESIGN CODES FOR TSUNAMI- RESISTING COASTAL STRUCTURES ISET Journal of Earthquake Technology, Paper No. 461, Vol. 42, No. 4, Decemer 25, pp. 137-145 A COMPUTATIONAL APPROACH TO DESIGN CODES FOR TSUNAMI- RESISTING COASTAL STRUCTURES Christopher Koutitas* and

More information

PHY451, Spring /5

PHY451, Spring /5 PHY451, Spring 2011 Notes on Optical Pumping Procedure & Theory Procedure 1. Turn on the electronics and wait for the cell to warm up: ~ ½ hour. The oven should already e set to 50 C don t change this

More information

Supporting Information. Dynamics of the Dissociating Uracil Anion. Following Resonant Electron Attachment

Supporting Information. Dynamics of the Dissociating Uracil Anion. Following Resonant Electron Attachment Supporting Information Dynamics of the Dissociating Uracil Anion Following Resonant Electron Attachment Y. Kawarai,, Th. Weer, Y. Azuma, C. Winstead, V. McKoy, A. Belkacem, and D.S. Slaughter, Department

More information

Modeling of Humidification in Comsol Multiphysics 4.4

Modeling of Humidification in Comsol Multiphysics 4.4 Modeling of Humidification in Comsol Multiphysics 4.4 Indrajit Wadgaonkar *1 and Suresh Arikapudi 1 1 Tata Motors Ltd. Pimpri, Pune, India, 411018. *Corresponding author: Indrajit Wadgaonkar, Tata Motors

More information

MATH 723 Mathematical Neuroscience Spring 2008 Instructor: Georgi Medvedev

MATH 723 Mathematical Neuroscience Spring 2008 Instructor: Georgi Medvedev MATH 723 Mathematical Neuroscience Spring 28 Instructor: Georgi Medvedev 2 Lecture 2. Approximate systems. 2. Reduction of the HH model to a 2D system. The original HH system consists of 4 differential

More information

Minimizing a convex separable exponential function subject to linear equality constraint and bounded variables

Minimizing a convex separable exponential function subject to linear equality constraint and bounded variables Minimizing a convex separale exponential function suect to linear equality constraint and ounded variales Stefan M. Stefanov Department of Mathematics Neofit Rilski South-Western University 2700 Blagoevgrad

More information

ASEISMIC DESIGN OF TALL STRUCTURES USING VARIABLE FREQUENCY PENDULUM OSCILLATOR

ASEISMIC DESIGN OF TALL STRUCTURES USING VARIABLE FREQUENCY PENDULUM OSCILLATOR ASEISMIC DESIGN OF TALL STRUCTURES USING VARIABLE FREQUENCY PENDULUM OSCILLATOR M PRANESH And Ravi SINHA SUMMARY Tuned Mass Dampers (TMD) provide an effective technique for viration control of flexile

More information

Determinants of generalized binary band matrices

Determinants of generalized binary band matrices Determinants of generalized inary and matrices Dmitry Efimov arxiv:17005655v1 [mathra] 18 Fe 017 Department of Mathematics, Komi Science Centre UrD RAS, Syktyvkar, Russia Astract Under inary matrices we

More information

Slope and wind effects on fire propagation

Slope and wind effects on fire propagation CSIRO PUBLISHING www.pulish.csiro.au/journals/ijwf International Journal of Wildland Fire, 24, 13, 143 156 Slope and wind effects on fire propagation Domingos X. Viegas Centro de Estudos sore Incêndios

More information

Spray evaporation model sensitivities

Spray evaporation model sensitivities Center for Turbulence Research Annual Research Briefs 20 23 Spray evaporation model sensitivities By Shashank, E. Knudsen AND H. Pitsch. Motivation and objective The energy density of solid- and liquid-phase

More information

The Role of the Graphing Calculator in the Qualitative Analysis of Function

The Role of the Graphing Calculator in the Qualitative Analysis of Function The Role of the Graphing Calculator in the Qualitative Analysis of Function Li Hong lihong2209@sina.com Beijing No.22 High School China Astract: This paper focuses on several cases of student s qualitative

More information

RATIONAL EXPECTATIONS AND THE COURNOT-THEOCHARIS PROBLEM

RATIONAL EXPECTATIONS AND THE COURNOT-THEOCHARIS PROBLEM RATIONAL EXPECTATIONS AND THE COURNOT-THEOCHARIS PROBLEM TÖNU PUU Received 18 April 006; Accepted 1 May 006 In dynamic models in economics, often rational expectations are assumed. These are meant to show

More information

ON THE COMPARISON OF BOUNDARY AND INTERIOR SUPPORT POINTS OF A RESPONSE SURFACE UNDER OPTIMALITY CRITERIA. Cross River State, Nigeria

ON THE COMPARISON OF BOUNDARY AND INTERIOR SUPPORT POINTS OF A RESPONSE SURFACE UNDER OPTIMALITY CRITERIA. Cross River State, Nigeria ON THE COMPARISON OF BOUNDARY AND INTERIOR SUPPORT POINTS OF A RESPONSE SURFACE UNDER OPTIMALITY CRITERIA Thomas Adidaume Uge and Stephen Seastian Akpan, Department Of Mathematics/Statistics And Computer

More information

Gravitational Model of the Three Elements Theory: Mathematical Explanations

Gravitational Model of the Three Elements Theory: Mathematical Explanations Journal of odern Physics, 3, 4, 7-35 http://dxdoiorg/436/jmp34738 Pulished Online July 3 (http://wwwscirporg/journal/jmp) Gravitational odel of the Three Elements Theory: athematical Explanations Frederic

More information

Appendix (Supplementary Material)

Appendix (Supplementary Material) Appendix (Supplementary Material) Appendix A. Real and Imaginary Parts of the complex modulus E*(f), and E*(f) for the Pyramidal Proe Tip The magnitude of the complex modulus, E*(f), and the phase angle,

More information

Numerical Study of Heat Propagation in Living Tissue Subjected to Instantaneous Heating

Numerical Study of Heat Propagation in Living Tissue Subjected to Instantaneous Heating Indian Journal of Biomechanics: Special Issue (NCBM 7-8 March 9) Numerical Study of Heat Propagation in Living Tissue Sujected to Instantaneous Heating P. R. Sharma 1, Sazid Ali, V. K. Katiyar 1 Department

More information

we make slices perpendicular to the x-axis. If the slices are thin enough, they resemble x cylinders or discs. The formula for the x

we make slices perpendicular to the x-axis. If the slices are thin enough, they resemble x cylinders or discs. The formula for the x Math Learning Centre Solids of Revolution When we rotate a curve around a defined ais, the -D shape created is called a solid of revolution. In the same wa that we can find the area under a curve calculating

More information

Exact Free Vibration of Webs Moving Axially at High Speed

Exact Free Vibration of Webs Moving Axially at High Speed Eact Free Viration of Wes Moving Aially at High Speed S. HATAMI *, M. AZHARI, MM. SAADATPOUR, P. MEMARZADEH *Department of Engineering, Yasouj University, Yasouj Department of Civil Engineering, Isfahan

More information

Estimation of fracture parameters from reflection seismic data Part II: Fractured models with orthorhombic symmetry

Estimation of fracture parameters from reflection seismic data Part II: Fractured models with orthorhombic symmetry GEOPHYSICS, VOL. 65, NO. 6 (NOVEMBER-DECEMBER 2000); P. 1803 1817, 5 FIGS., 1 TABLE. Estimation of fracture parameters from reflection seismic data Part II: Fractured models with orthorhomic symmetry Andrey

More information

TIGHT BOUNDS FOR THE FIRST ORDER MARCUM Q-FUNCTION

TIGHT BOUNDS FOR THE FIRST ORDER MARCUM Q-FUNCTION TIGHT BOUNDS FOR THE FIRST ORDER MARCUM Q-FUNCTION Jiangping Wang and Dapeng Wu Department of Electrical and Computer Engineering University of Florida, Gainesville, FL 3611 Correspondence author: Prof.

More information

Modeling of Water Flows around a Circular Cylinder with the SPH Method

Modeling of Water Flows around a Circular Cylinder with the SPH Method Archives of Hydro-Engineering and Environmental Mechanics Vol. 61 (2014), No. 1 2, pp. 39 60 DOI: 10.1515/heem-2015-0003 IBW PAN, ISSN 1231 3726 Modeling of Water Flows around a Circular Cylinder with

More information

Air and Heat Flow through Large Vertical Openings

Air and Heat Flow through Large Vertical Openings Air and Heat Flow through Large Vertical Openings J L M Hensen *, J van der Maas #, A Roos * * Eindhoven University of Technology # Ecole Polytechnique Federale de Lausanne After a short description of

More information

Introduction to Acoustics. Phil Joseph

Introduction to Acoustics. Phil Joseph Introduction to Acoustics Phil Joseph INTRODUCTION TO ACOUSTICS Sound and Noise Sound waves Frequency, wavelength and wavespeed Point sources Sound power and intensity Wave reflection Standing waves Measures

More information

Estimating a Finite Population Mean under Random Non-Response in Two Stage Cluster Sampling with Replacement

Estimating a Finite Population Mean under Random Non-Response in Two Stage Cluster Sampling with Replacement Open Journal of Statistics, 07, 7, 834-848 http://www.scirp.org/journal/ojs ISS Online: 6-798 ISS Print: 6-78X Estimating a Finite Population ean under Random on-response in Two Stage Cluster Sampling

More information

International Journal of Thermal Sciences

International Journal of Thermal Sciences International Journal of Thermal Sciences 49 (21) 643 652 Contents lists availale at ScienceDirect International Journal of Thermal Sciences journal homepage: www.elsevier.com/locate/ijts Investigation

More information

QUADRATIC EQUATIONS EXPECTED BACKGROUND KNOWLEDGE

QUADRATIC EQUATIONS EXPECTED BACKGROUND KNOWLEDGE 6 QUADRATIC EQUATIONS In this lesson, you will study aout quadratic equations. You will learn to identify quadratic equations from a collection of given equations and write them in standard form. You will

More information

Microwave Absorption by Light-induced Free Carriers in Silicon

Microwave Absorption by Light-induced Free Carriers in Silicon Microwave Asorption y Light-induced Free Carriers in Silicon T. Sameshima and T. Haa Tokyo University of Agriculture and Technology, Koganei, Tokyo 184-8588, Japan E-mail address: tsamesim@cc.tuat.ac.jp

More information

Condition number of the BEM matrix arising from the Stokes equations in 2D

Condition number of the BEM matrix arising from the Stokes equations in 2D Condition numer of the BEM matrix arising from the Stokes equations in D W. Dijkstra a, R.M.M. Mattheij Department of Mathematics and Computing Science, Eindhoven University of Technology, P.O. Box 53,

More information

Spiking problem in monotone regression : penalized residual sum of squares

Spiking problem in monotone regression : penalized residual sum of squares Spiking prolem in monotone regression : penalized residual sum of squares Jayanta Kumar Pal 12 SAMSI, NC 27606, U.S.A. Astract We consider the estimation of a monotone regression at its end-point, where

More information

Exact Shape Functions for Timoshenko Beam Element

Exact Shape Functions for Timoshenko Beam Element IOSR Journal of Computer Engineering (IOSR-JCE) e-iss: 78-66,p-ISS: 78-877, Volume 9, Issue, Ver. IV (May - June 7), PP - www.iosrjournals.org Exact Shape Functions for Timoshenko Beam Element Sri Tudjono,

More information

Multiattribute seismic analysis on AVO-derived parameters A case study

Multiattribute seismic analysis on AVO-derived parameters A case study Multiattriute seismic analysis on AVO-derived parameters A case study SATINDER CHOPRA, Core La Reservoir Technologies, Calgary, Canada DOUG PRUDEN, GEDCO, Calgary, Canada Prospecting for reservoir zones

More information

Modal analysis of waveguide using method of moment

Modal analysis of waveguide using method of moment HAIT Journal of Science and Engineering B, Volume x, Issue x, pp. xxx-xxx Copyright C 27 Holon Institute of Technology Modal analysis of waveguide using method of moment Arti Vaish and Harish Parthasarathy

More information

S. Kadowaki, S.H. Kim AND H. Pitsch. 1. Motivation and objectives

S. Kadowaki, S.H. Kim AND H. Pitsch. 1. Motivation and objectives Center for Turbulence Research Annual Research Briefs 2005 325 The dynamics of premixed flames propagating in non-uniform velocity fields: Assessment of the significance of intrinsic instabilities in turbulent

More information

UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION

UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume, Numer 4/0, pp. 9 95 UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION

More information

A Sufficient Condition for Optimality of Digital versus Analog Relaying in a Sensor Network

A Sufficient Condition for Optimality of Digital versus Analog Relaying in a Sensor Network A Sufficient Condition for Optimality of Digital versus Analog Relaying in a Sensor Network Chandrashekhar Thejaswi PS Douglas Cochran and Junshan Zhang Department of Electrical Engineering Arizona State

More information

Differential Flatness-based Kinematic and Dynamic Control of a Differentially Driven Wheeled Mobile Robot

Differential Flatness-based Kinematic and Dynamic Control of a Differentially Driven Wheeled Mobile Robot Differential Flatness-ased Kinematic and Dynamic Control of a Differentially Driven Wheeled Moile Root Chin Pei Tang, IEEE Memer Astract In this paper, we present an integrated motion planning and control

More information

n v molecules will pass per unit time through the area from left to

n v molecules will pass per unit time through the area from left to 3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular

More information

Pairwise Interaction Extended Point-Particle (PIEP) Model for droplet-laden flows: Towards application to the mid-field of a spray

Pairwise Interaction Extended Point-Particle (PIEP) Model for droplet-laden flows: Towards application to the mid-field of a spray Pairwise Interaction Extended Point-Particle (PIEP) Model for droplet-laden flows: Towards application to the mid-field of a spray Georges Akiki, Kai Liu and S. Balachandar * Department of Mechanical &

More information

IMPLEMENTATION OF FULLY COUPLED HEAT AND MASS TRANSPORT MODEL TO DETERMINE TEMPERTATURE AND MOISTURE STATE AT ELEVATED TEMPERATURES.

IMPLEMENTATION OF FULLY COUPLED HEAT AND MASS TRANSPORT MODEL TO DETERMINE TEMPERTATURE AND MOISTURE STATE AT ELEVATED TEMPERATURES. 11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM ) 6th European Conference on Computational Fluid Dynamics (ECFD I) IMPLEMENTATION OF FULLY

More information

Optimal Routing in Chord

Optimal Routing in Chord Optimal Routing in Chord Prasanna Ganesan Gurmeet Singh Manku Astract We propose optimal routing algorithms for Chord [1], a popular topology for routing in peer-to-peer networks. Chord is an undirected

More information

Motivation: Can the equations of physics be derived from information-theoretic principles?

Motivation: Can the equations of physics be derived from information-theoretic principles? 10. Physics from Fisher Information. Motivation: Can the equations of physics e derived from information-theoretic principles? I. Fisher Information. Task: To otain a measure of the accuracy of estimated

More information

Comments on A Time Delay Controller for Systems with Uncertain Dynamics

Comments on A Time Delay Controller for Systems with Uncertain Dynamics Comments on A Time Delay Controller for Systems with Uncertain Dynamics Qing-Chang Zhong Dept. of Electrical & Electronic Engineering Imperial College of Science, Technology, and Medicine Exhiition Rd.,

More information

Chapter 2 Canonical Correlation Analysis

Chapter 2 Canonical Correlation Analysis Chapter 2 Canonical Correlation Analysis Canonical correlation analysis CCA, which is a multivariate analysis method, tries to quantify the amount of linear relationships etween two sets of random variales,

More information