Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow Theoretical Analysis


 Clemence Cobb
 9 months ago
 Views:
Transcription
1 Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow Theoretical Analysis DAVID KATOSHEVSKI Department of Biotechnology and Environmental Engineering BenGurion niversity of the Negev BeerSheva 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 multisize 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 nongrouping 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 timeehavior. KeyWords: 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 fueldroplet 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 multisize 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 multisize 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 travelgrouping 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 hostgas 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 onedimensional 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 xaxis and a period T 2 / in time. The wave that is formed propagates along the xaxis 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 travelgroup, 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 nonfamily type, will e oserved. In addition to these modes of grouping, for certain operating conditions clear nongrouping 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 nongrouping, 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 nongrouping, or nearly nongrouping situations even for a wider range of droplet Reynolds numers we have to maintain the following relation arising from Eqs. 57, a 1 NG St 1 (8) where NG stands for NonGrouping. 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 nongrouping 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/nongrouping 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 wavelength (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 nongrouping 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 xt) 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 cleargrouping 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 nonevaporating droplets of 30 microns in diameter, and 1000 Hz.
4 Fig. 3. Nongrouping mode, for nonevaporating droplets of 30 microns in diameter, and 1800 Hz. Fig. 4. Nongrouping mode, for nonevaporating 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 monodisperse 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 multisized 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 monosize spray, which can e used also for the case where similarsized groups of droplets of a multisize 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 RanzMarshal 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 coworkers (Sirignano, 1999, and references therein) for evaporation and comustion of nonisolate 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 multisized 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 monosized 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 Dsquare 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 nongrouping droplets. Thus, for large droplets, which initially do not form a travelgroup, evaporation and coincidently reduction of size, will not lead to grouping ehavior. Now, we add to the aove the effect of the slowdown 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, slowdown 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 Dsquare 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 suindex d denotes dropletsize) 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 sizesgroups together as a function of time m(t) can e evaluated y Eqs. (14) and (18). As a firstorder 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 multisize spray droplets are sujected to an oscillatory flow field. The current work also points out that once grouping or nongrouping 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 BenGurion 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 BurkeSchumann 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 Itamar Hite, Michal Ruzal Mendelevich,, Tal Shakked, Prof. Eran Sher and Prof. David Katoshevski BenGurion niversity of the Negev BeerSheva
More informationModeling Droplet Breakup Processes in Biofuel Diesel Engines under Microexplosion Conditions
ILASS Americas, th Annual Conference on Liquid Atomization and Spray Systems, Chicago, IL, May 7 Modeling Droplet Breakup Processes in Biofuel Diesel Engines under Microexplosion Conditions K. T. Wang
More informationOn 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: 19966 Vol. 6, Issue 1 (June 11) pp. 7 (Previously, Vol. 6, Issue 11, pp. 1767 178) Applications and Applied Mathematics: An International Journal
More informationSection 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 MassSpring Systems In this section we will investigate two applications of our work in Section 8.4. First, we will consider
More information1D 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 informationAvailable online at Energy Procedia 100 (2009) (2008) GHGT9
Availale online at www.sciencedirect.com Energy Procedia (29) (28) 655 66 Energy Procedia www.elsevier.com/locate/procedia www.elsevier.com/locate/xxx GHGT9 Precomustion CO 2 capture for IGCC plants
More informationLecture 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 informationEffect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium
IOSR Journal of Mathematics (IOSRJM) eissn: 78578, pissn: 39765X. Volume, Issue Ver. III (Jan.  Fe. 06), 565 www.iosrjournals.org Effect of Uniform Horizontal Magnetic Field on Thermal Instaility
More informationTopics in Other Lectures Droplet Groups and Array Instability of Injected Liquid Liquid FuelFilms
Lecture Topics Transient Droplet Vaporization Convective Vaporization Liquid Circulation Transcritical Thermodynamics Droplet Drag and Motion Spray Computations Turbulence Effects Topics in Other Lectures
More informationLuis 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 [JIA01] 1. Introduction The evolution of VLSI chips toward larger die sizes and faster clock speeds makes
More informationRobot 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 informationA 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 informationCAV2001:SESSIONA
CAV001:SESSIONA1.004 1 Effects of model size and free stream nuclei on tip vortex cavitation inception scaling Young Shen*, Georges Chahine**, ChaoTsung Hsiao**, and Stuart Jessup* *NSWCCD, West Bethesda,
More informationMath 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 informationNonsteady interior ballistics of cylindricalgrain solid rocket motors
Computational Ballistics II 281 Nonsteady interior allistics of cylindricalgrain solid rocket motors D. R. Greatrix Department of Aerospace Engineering, Ryerson University, Canada Astract A numerical
More informationBuckling 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 informationEFFECTS 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 informationCritical 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 informationLecture 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 (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 informationPARAMETER 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@RoseHulman.edu ABSTRACT In this paper a freely rotating pendulum suspended from
More informationarxiv: v1 [hepph] 20 Dec 2012
Decemer 21, 2012 UMDPP012027 Using Energy Peaks to Count Dark Matter Particles in Decays arxiv:1212.5230v1 [hepph] 20 Dec 2012 Kaustuh Agashe a, Roerto Franceschini a, Doojin Kim a, and Kyle Wardlow
More informationFieldinduced magnetization of a freeelectron gas in thin films
Materials SciencePoland, Vol. 5, No., 007 ieldinduced magnetization of a freeelectron gas in thin films T. BALCERZAK 1* 1 Department of Solid State Physics, University of Łódź, Pomorska 149/153, 9036
More informationcen29305_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 informationTheoretical 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 informationConcentration 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 informationPair diffusion in quasione and quasitwodimensional binary colloid suspensions
THE JOURNAL OF CHEMICAL PHYSICS 126, 134908 2007 Pair diffusion in quasione and quasitwodimensional binary colloid suspensions David T. Valley, Stuart A. Rice, Bianxiao Cui, a and Hau My Ho Department
More informationTOPICAL PROBLEMS OF FLUID MECHANICS 17 ONEDIMENSIONAL TEMPERATURE DISTRIBUTION OF CONDENSING ANNULAR FINS OF DIFFERENT PROFILES
TOPICAL PROBLEMS OF FLUID MECHANICS 17 ONEDIMENSIONAL 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 informationSEG/New Orleans 2006 Annual Meeting. Nonorthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University
Nonorthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University SUMMARY Wavefield extrapolation is implemented in nonorthogonal Riemannian spaces. The key component is the development
More informationMixing 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. SapirKatiraie The BenGurion University of the Negev, BeerSheva, Israel G. Hazak
More informationOptimal 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 AlcatelLucent Bell Las, 916 Nozay, France Department of Telecommunications,
More informationMotion 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 informationHot 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 informationThreering tautomerism of the 2 isoxazoline  2 pyrazoline  1,3,4thiadiazine system
Threering tautomerism of the 2 isoxazoline  2 pyrazoline  1,3,4thiadiazine system Andrei Yu. Ershov, Maksim V. Mokeev, Elena V. Beloorodova, and Alexander V. Grianov Institute of Macromolecular Compounds
More informationDynamical Systems Solutions to Exercises
Dynamical Systems Part 56 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 informationOnedimensional Spray Combustion Optimization with a Sequential Linear Quadratic Algorithm
Onedimensional 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 informationBOUSSINESQTYPE MOMENTUM EQUATIONS SOLUTIONS FOR STEADY RAPIDLY VARIED FLOWS. Yebegaeshet T. Zerihun 1 and John D. Fenton 2
ADVANCES IN YDROSCIENCE AND ENGINEERING, VOLUME VI BOUSSINESQTYPE MOMENTUM EQUATIONS SOLUTIONS FOR STEADY RAPIDLY VARIED FLOWS Yeegaeshet T. Zerihun and John D. Fenton ABSTRACT The depth averaged SaintVenant
More informationChaos 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 informationNumerical Studies of Droplet Deformation and Breakup
ILASS Americas 14th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2001 Numerical Studies of Droplet Deformation and Breakup B. T. Helenbrook Department of Mechanical and
More informationarxiv:physics/ v1 [physics.plasmph] 7 Apr 2006
Europhysics Letters PREPRINT arxiv:physics/6455v [physics.plasmph] 7 Apr 26 Olique electromagnetic instailities for an ultra relativistic electron eam passing through a plasma A. Bret ETSI Industriales,
More informationMathematical 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 informationAMS 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 informationFLOW AND HEATTRANSFER MODELLING OF THREEDIMENSIONAL JET IMPINGEMENT ON A CONCAVE SURFACE
FLOW AND HEATTRANSFER MODELLING OF THREEDIMENSIONAL JET IMPINGEMENT ON A CONCAVE SURFACE Tim J. Craft Hector Iacovides Nor A. Mostafa School of Mechanical Aerospace & Civil Engineering The University
More informationCharacteristics of flow over the free overfall of triangular channel
MATEC We of Conferences 16, 03006 (018) BCEE3017 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 informationP = ρ{ 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 informationA Stable and Convergent Finite Difference Scheme for 2D Incompressible Nonlinear Viscoelastic Fluid Dynamics Problem
Applied and Computational Mathematics 2018; (1): 1118 http://www.sciencepulishinggroup.com/j/acm doi: 10.11648/j.acm.2018001.12 ISSN: 23285605 (Print); ISSN: 23285613 (Online) A Stale and Convergent
More informationModule 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 informationCombustion 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 informationTP 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 informationCULAR 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 informationEstimation 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 informationragsdale (zdr82) HW7 ditmire (58335) 1 The magnetic force is
ragsdale (zdr8) HW7 ditmire (585) This printout should have 8 questions. Multiplechoice questions ma continue on the net column or page find all choices efore answering. 00 0.0 points A wire carring
More informationOPTIMAL WAVELENGTH SELECTION ALGORITHM OF NONSPHERICAL PARTICLE SIZE DISTRIBUTION BASED ON THE LIGHT EXTINCTION DATA
THERMAL SCIENCE, Year 202, Vol. 6, No. 5, pp. 353357 353 OPTIMAL WAVELENGTH SELECTION ALGORITHM OF NONSPHERICAL PARTICLE SIZE ISTRIBUTION BASE ON THE LIGHT EXTINCTION ATA y Hong TANG * College of Metrology
More informationSuperluminal 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 informationNotes 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 informationVelocity crosscorrelations and atomic momentum transfer in simple liquids with different potential cores
PHYSICAL REVIEW E VOLUME 62, NUMBER 1 JULY 2000 Velocity crosscorrelations and atomic momentum transfer in simple liquids with different potential cores A. Verdaguer and J. A. Padró Departament de Física
More informationSTUDY ON BUBBLE BEHAVIOR OF INERT GASES AT ENTRANCE NOZZLE IN SODIUMCOOLED FAST REACTOR
NTHAS8: The Eighth JapanKorea Symposium on Nuclear Thermal Hydraulics and Safety Beppu, Japan, Decemer 912, 2012 Paper Numer N8P1081 STUDY ON BUBBLE BEHAVIOR OF INERT GASES AT ENTRANCE NOZZLE IN SODIUMCOOLED
More informationA 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. 137145 A COMPUTATIONAL APPROACH TO DESIGN CODES FOR TSUNAMI RESISTING COASTAL STRUCTURES Christopher Koutitas* and
More informationPHY451, 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 informationSupporting 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 informationModeling 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 informationMATH 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 informationMinimizing 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 SouthWestern University 2700 Blagoevgrad
More informationASEISMIC 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 informationDeterminants 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 informationSlope 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 informationSpray 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 liquidphase
More informationThe 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 informationRATIONAL EXPECTATIONS AND THE COURNOTTHEOCHARIS PROBLEM
RATIONAL EXPECTATIONS AND THE COURNOTTHEOCHARIS 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 informationON 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 informationGravitational Model of the Three Elements Theory: Mathematical Explanations
Journal of odern Physics, 3, 4, 735 http://dxdoiorg/436/jmp34738 Pulished Online July 3 (http://wwwscirporg/journal/jmp) Gravitational odel of the Three Elements Theory: athematical Explanations Frederic
More informationAppendix (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 informationNumerical Study of Heat Propagation in Living Tissue Subjected to Instantaneous Heating
Indian Journal of Biomechanics: Special Issue (NCBM 78 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 informationwe make slices perpendicular to the xaxis. 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 informationExact 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 informationEstimation of fracture parameters from reflection seismic data Part II: Fractured models with orthorhombic symmetry
GEOPHYSICS, VOL. 65, NO. 6 (NOVEMBERDECEMBER 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 informationTIGHT BOUNDS FOR THE FIRST ORDER MARCUM QFUNCTION
TIGHT BOUNDS FOR THE FIRST ORDER MARCUM QFUNCTION Jiangping Wang and Dapeng Wu Department of Electrical and Computer Engineering University of Florida, Gainesville, FL 3611 Correspondence author: Prof.
More informationModeling of Water Flows around a Circular Cylinder with the SPH Method
Archives of HydroEngineering and Environmental Mechanics Vol. 61 (2014), No. 1 2, pp. 39 60 DOI: 10.1515/heem20150003 IBW PAN, ISSN 1231 3726 Modeling of Water Flows around a Circular Cylinder with
More informationAir 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 informationIntroduction 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 informationEstimating a Finite Population Mean under Random NonResponse in Two Stage Cluster Sampling with Replacement
Open Journal of Statistics, 07, 7, 834848 http://www.scirp.org/journal/ojs ISS Online: 6798 ISS Print: 678X Estimating a Finite Population ean under Random onresponse in Two Stage Cluster Sampling
More informationInternational 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 informationQUADRATIC 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 informationMicrowave Absorption by Lightinduced Free Carriers in Silicon
Microwave Asorption y Lightinduced Free Carriers in Silicon T. Sameshima and T. Haa Tokyo University of Agriculture and Technology, Koganei, Tokyo 1848588, Japan Email address: tsamesim@cc.tuat.ac.jp
More informationCondition 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 informationSpiking 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 endpoint, where
More informationExact Shape Functions for Timoshenko Beam Element
IOSR Journal of Computer Engineering (IOSRJCE) eiss: 7866,pISS: 78877, Volume 9, Issue, Ver. IV (May  June 7), PP  www.iosrjournals.org Exact Shape Functions for Timoshenko Beam Element Sri Tudjono,
More informationMultiattribute seismic analysis on AVOderived parameters A case study
Multiattriute seismic analysis on AVOderived parameters A case study SATINDER CHOPRA, Core La Reservoir Technologies, Calgary, Canada DOUG PRUDEN, GEDCO, Calgary, Canada Prospecting for reservoir zones
More informationModal analysis of waveguide using method of moment
HAIT Journal of Science and Engineering B, Volume x, Issue x, pp. xxxxxx Copyright C 27 Holon Institute of Technology Modal analysis of waveguide using method of moment Arti Vaish and Harish Parthasarathy
More informationS. 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 nonuniform velocity fields: Assessment of the significance of intrinsic instabilities in turbulent
More informationUNSTEADY 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 informationA 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 informationDifferential Flatnessbased Kinematic and Dynamic Control of a Differentially Driven Wheeled Mobile Robot
Differential Flatnessased 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 informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of OneDimensional Gas Flow The dissipative processes  viscosity (internal friction) and heat conduction  are connected with existence of molecular
More informationPairwise Interaction Extended PointParticle (PIEP) Model for dropletladen flows: Towards application to the midfield of a spray
Pairwise Interaction Extended PointParticle (PIEP) Model for dropletladen flows: Towards application to the midfield of a spray Georges Akiki, Kai Liu and S. Balachandar * Department of Mechanical &
More informationIMPLEMENTATION 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 informationOptimal 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 peertopeer networks. Chord is an undirected
More informationMotivation: Can the equations of physics be derived from informationtheoretic principles?
10. Physics from Fisher Information. Motivation: Can the equations of physics e derived from informationtheoretic principles? I. Fisher Information. Task: To otain a measure of the accuracy of estimated
More informationComments on A Time Delay Controller for Systems with Uncertain Dynamics
Comments on A Time Delay Controller for Systems with Uncertain Dynamics QingChang Zhong Dept. of Electrical & Electronic Engineering Imperial College of Science, Technology, and Medicine Exhiition Rd.,
More informationChapter 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